Showing 1–50 of 75 results for author: Fernandez-Lopez, M

Magnetic Fields in Massive Star-forming Regions (MagMaR) IV: Tracing the Magnetic Fields in the O-type protostellar system IRAS 16547$-$4247

Authors: Luis A. Zapata, Manuel Fernández-López, Patricio Sanhueza, Josep M. Girart, Luis F. Rodríguez, Paulo Cortes, Koch Patrick, María T. Beltrán, Kate Pattle, Henrik Beuther, Piyali Saha, Wenyu Jiao, Fengwei Xu, Xing Walker Lu, Fernando Olguin, Shanghuo Li, Ian W. Stephens, Ji-hyun Kang, Yu Cheng, Spandan Choudhury, Kaho Morii, Eun Jung Chung, Jia-Wei Wang, Jihye Hwang, A-Ran Lyo , et al. (2 additional authors not shown)

Abstract: The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4… ▽ More The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) 1.2 mm full polarized continuum, and H$^{13}$CO$^+$(3$-$2), CS(5$-$4), and HN$^{13}$C(3$-$2) line observations with a high angular resolution ($\sim$0.4$''$ or 1100 au). In the 1.2 mm continuum emission, we reveal a dusty envelope surrounding the massive protostars, IRAS16547-E and IRAS16547-W, with dimensions of $\sim$10,000 au. This envelope has a bi-conical structure likely carved by the powerful thermal radio jet present in region. The magnetic fields vectors follow very-well the bi-conical envelope. The polarization fraction is $\sim$2.0\% in this region. Some of these vectors seem to converge to IRAS 16547-E, and IRAS 16547-W, the most massive protostars. Moreover, the velocity fields revealed from the spectral lines H$^{13}$CO$^+$(3$-$2), and HN$^{13}$C(3$-$2) show velocity gradients with a good correspondence with the magnetic fields, that maybe are tracing the cavities of molecular outflows or maybe in some parts infall. We derived a magnetic field strength in some filamentary regions that goes from 2 to 6.1\,mG. We also find that the CS(5$-$4) molecular line emission reveals multiple outflow cavities or bow-shocks with different orientations, some of which seem to follow the NW-SE radio thermal jet. △ Less

Submitted 19 August, 2024; originally announced August 2024.

Comments: Accepted by the Astrophysical Journal, 13 pages

The Class 0 protostars in Orion: Characterizing the properties of their magnetized envelopes

Authors: B. Huang, J. M. Girart, I. W. Stephens, M. Fernandez-Lopez, J. J. Tobin, P. Cortes, N. M. Murillo, P. C. Myers, S. Sadavoy, Q. Zhang, H. G. Arce, J. M. Carpenter, W. Kwon, V. J. M. Le Gouellec, Z. -Y. Li, L. W. Looney, T. Megeath, E. G. Cox, N. Karnath, D. Segura-Cox

Abstract: We present a study connecting the physical properties of protostellar envelopes to the morphology of the envelope-scale magnetic field. We use the ALMA polarization observations of 55 young prtostars at 0.87 mm on $\sim400-3000$ au scales from the {\em B}-field Orion Protostellar Survey (BOPS) to infer the envelope-scale magnetic field and both dust and gas emission on comparable scales to measure… ▽ More We present a study connecting the physical properties of protostellar envelopes to the morphology of the envelope-scale magnetic field. We use the ALMA polarization observations of 55 young prtostars at 0.87 mm on $\sim400-3000$ au scales from the {\em B}-field Orion Protostellar Survey (BOPS) to infer the envelope-scale magnetic field and both dust and gas emission on comparable scales to measure the envelope properties. We find that the protostellar envelopes with compact polarized dust emission tend to have lower envelope masses, than the sources with more extended envelopes. We also find that protostars showing hourglass-field morphologies tend to have lower velocity dispersions in their envelopes, whereas systems with spiral-field morphologies have higher velocity dispersion. Combining with the disk properties taken from the Orion VLA/ALMA Nascent Disk and Multiplicity (VANDAM) survey, we connect envelope properties to fragmentation. Our results suggest that envelope mass may not correlate with fragmentation, whereas turbulence appears to promote fragmentation. On the other hand, we find that fragmentation is suppressed in systems with pinched magnetic fields, suggesting that the magnetic field play a role on providing additional support against gravitational collapse, and the formation of an hourglass-like field may coincide with enhanced magnetic braking that removes angular momentum and hinders the formation of embedded disks. Nevertheless, significant misalignment between magnetic field and outflow axes tends to reduce magnetic braking, leading to the formation of larger disks. △ Less

Submitted 28 July, 2024; originally announced July 2024.

Comments: 18 pages, 6 figures, submitted to ApJ

ALMA-IMF XV: The core mass function in the high-mass star-formation regime

Authors: F. Louvet, P. Sanhueza, A. Stutz, A. Men'shchikov, F. Motte, R. Galván-Madrid, S. Bontemps, Y. Pouteau, A. Ginsburg, T. Csengeri, J. Di Francesco, P. Dell'Ova, M. González, P. Didelon, J. Braine, N. Cunningham, B. Thomasson, P. Lesaffre, P. Hennebelle, M. Bonfand, A. Gusdorf, R. H. Álverez-Gutiérrez, T. Nony, G. Busquet, F. Olguin , et al. (16 additional authors not shown)

Abstract: The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Halpha emission to estimate the star formation rate of a system and interpret unresolved star clusters across the universe. So far, there is little firm evidence of large-scale variations of the I… ▽ More The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Halpha emission to estimate the star formation rate of a system and interpret unresolved star clusters across the universe. So far, there is little firm evidence of large-scale variations of the IMF, which is thus generally considered universal. Stars form from cores and it is now possible to estimate core masses and compare the core mass function (CMF) with the IMF, which it presumably produces. The goal of the ALMA-IMF large program is to measure the core mass function at high linear resolution (2700 au) in 15 typical Milky Way protoclusters spanning a mass range of 2500 to 32700 Msun. In this work, we used two different core extraction algorithms to extract about 680 gravitationally bound cores from these 15 protoclusters. We adopt per core temperature using the temperature estimate from the PPMAP Bayesian method. A power-law fit to the CMF of the sub-sample of cores above the 1.64 Msun completeness limit, 330 cores, through the maximum likelihood estimate technique yields a slope of 1.97 +/- 0.06, significantly flatter than the 2.35 Salpeter slope. Assuming a self-similar mapping between the CMF and the IMF, this result implies that these 15 high-mass protoclusters will generate atypical IMFs. This sample is the largest to date produced and analysed self-consistently, derived at matched physical resolution, with per-core temperature estimates and cores as massive as 150 Msun. We provide the raw source extraction catalogues and the source derived size, temperature, mass, and spectral indices in the 15 protoclusters. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Magnetic Fields in Massive Star-forming Regions (MagMaR): Unveiling an Hourglass Magnetic Field in G333.46-0.16 using ALMA

Authors: Piyali Saha, Patricio Sanhueza, Marco Padovani, Josep M. Girart, Paulo Cortes, Kaho Morii, Junhao Liu, A. Sanchez-Monge, Daniele Galli, Shantanu Basu, Patrick M. Koch, Maria T. Beltran, Shanghuo Li, Henrik Beuther, Ian W. Stephens, Fumitaka Nakamura, Qizhou Zhang, Wenyu Jiao, M. Fernandez-Lopez, Jihye Hwang, Eun Jung Chung, Kate Pattle, Luis A. Zapata, Fengwei Xu, Fernando A. Olguin , et al. (11 additional authors not shown)

Abstract: The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as par… ▽ More The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{POS}$, obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis-Chandrasekhar-Fermi method, the magnetic field strength towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Badminton Birdie-Like Aerodynamic Alignment of Drifting Dust Grains by Subsonic Gaseous Flows in Protoplanetary Disks

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Haifeng Yang, Leslie W. Looney, Ian W. Stephens, Manuel Fernández-López, Rachel E. Harrison

Abstract: Recent (sub)millimeter polarization observations of protoplanetary disks reveal toroidally aligned, effectively prolate dust grains large enough (at least ~100 $μ$m) to efficiently scatter millimeter light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas-dust relative motion when the grain's center of mass is offset from… ▽ More Recent (sub)millimeter polarization observations of protoplanetary disks reveal toroidally aligned, effectively prolate dust grains large enough (at least ~100 $μ$m) to efficiently scatter millimeter light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas-dust relative motion when the grain's center of mass is offset from its geometric center, analogous to a badminton birdie's alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric center in the flow direction relative to its center of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the center of mass can be shifted away from the spheroid's geometric center. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the center of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disk model shows that the alignment timescales are shorter than or comparable to the orbital time. The grain alignment direction in a disk depends on the disk (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas disks and for small St grains near the gas pressure extrema, such as rings and gaps. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 21 pages, 12 figures, submitted to MNRAS

ALMA-IMF XII: Point-process mapping of 15 massive protoclusters

Authors: P. Dell'Ova, F. Motte, A. Gusdorf, Y. Pouteau, A. Men'shchikov, D. Diaz-Gonzalez, R. Galván-Madrid, P. Lesaffre, P. Didelon, A. M. Stutz, A. P. M. Towner, K. Marsh, A. Whitworth, M. Armante, M. Bonfand, T. Nony, M. Valeille-Manet, S. Bontemps, T. Csengeri, N. Cunningham, A. Ginsburg, F. Louvet, R. H. Alvarez-Gutierrez, N. Brouillet, J. Salinas , et al. (7 additional authors not shown)

Abstract: A crucial aspect in addressing the challenge of measuring the core mass function, that is pivotal for comprehending the origin of the initial mass function, lies in constraining the temperatures of the cores. We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. High angular resolution mapping is required to captur… ▽ More A crucial aspect in addressing the challenge of measuring the core mass function, that is pivotal for comprehending the origin of the initial mass function, lies in constraining the temperatures of the cores. We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. High angular resolution mapping is required to capture the properties of protostellar and pre-stellar cores and to effectively separate them from larger features, such as dusty filaments. We employed the point process mapping (PPMAP) technique, enabling us to perform spectral energy distribution fitting of far-infrared and submillimeter observations across the 15 ALMA-IMF fields, at an unmatched 2.5" angular resolution. By combining the modified blackbody model with near-infrared data, we derived bolometric luminosity maps. We estimated the errors impacting values of each pixel in the temperature, column density, and luminosity maps. Subsequently, we employed the extraction algorithm getsf on the luminosity maps in order to detect luminosity peaks and measure their associated masses. We obtained high-resolution constraints on the luminosity, dust temperature, and mass of protoclusters, that are in agreement with previously reported measurements made at a coarser angular resolution. We find that the luminosity-to-mass ratio correlates with the evolutionary stage of the studied regions, albeit with intra-region variability. We compiled a PPMAP source catalog of 313 luminosity peaks using getsf on the derived bolometric luminosity maps. The PPMAP source catalog provides constraints on the mass and luminosity of protostars and cores, although one source may encompass several objects. Finally, we compare the estimated luminosity-to-mass ratio of PPMAP sources with evolutionary tracks and discuss the limitations imposed by the 2.5" beam. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 35 pages, 19 figures, 5 tables Accepted by A&A

ALMA-IMF XIV: Free-Free Templates Derived from H$41α$ and Ionized Gas Content in Fifteen Massive Protoclusters

Authors: Roberto Galván-Madrid, Daniel J. Díaz-González, Frédérique Motte, Adam Ginsburg, Nichol Cunningham, Karl M. Menten, Mélanie Armante, Mélisse Bonfand, Jonathan Braine, Timea Csengeri, Pierre Dell'Ova, Fabien Louvet, Thomas Nony, Rudy Rivera-Soto, Patricio Sanhueza, Amelia M. Stutz, Friedrich Wyrowski, Rodrigo H. Álvarez-Gutiérrez, Tapas Baug, Sylvain Bontemps, Leonardo Bronfman, Manuel Fernández-López, Antoine Gusdorf, Atanu Koley, Hong-Li Liu , et al. (3 additional authors not shown)

Abstract: We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen i… ▽ More We use the H$41α$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41α$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41α$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $ΔV_\mathrm{dyn}$ with $D$ and $Q_0$. $ΔV_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

MagMar III -- Resisting the Pressure, Is the Magnetic Field Overwhelmed in NGC6334I?

Authors: Paulo C. Cortes, Josep M. Girart, Patricio Sanhueza, Junhao Liu, Sergio Martin, Ian W. Stephens, Henrik Beuther, Patrick M. Koch, M. Fernandez-Lopez, Alvaro Sanchez-Monge, Jia-Wei Wang, Kaho Morii, Shanghuo Li, Piyali Saha, Qizhou Zhang, David Rebolledo, Luis A. Zapata, Ji-hyun Kang, Wenyu Jiao, Jongsoo Kim, Yu Cheng, Jihye Hwang, Eun Jung Chung, Spandan Choudhury, A-Ran Lyo , et al. (1 additional authors not shown)

Abstract: We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a pr… ▽ More We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Accepted for Publication at the Astrophysical Journal

Protoplanetary Disk Polarization at Multiple Wavelengths: Are Dust Populations Diverse?

Authors: Rachel E. Harrison, Zhe-Yu Daniel Lin, Leslie W. Looney, Zhi-Yun Li, Haifeng Yang, Ian Stephens, Manuel Fernández-López

Abstract: Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $μ$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can… ▽ More Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $μ$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $μ$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $μ$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $μ$m ) near the disk midplane and small grains ($a_{max}$ = 140 $μ$m) above and below the midplane. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 22 pages, 12 figures

ALMA-IMF XIII: N$_2$H$^+$ kinematic analysis on the intermediate protocluster G353.41

Authors: R. H. Álvarez-Gutiérrez, A. M. Stutz, N. Sandoval-Garrido, F. Louvet, F. Motte, R. Galván-Madrid, N. Cunningham, P. Sanhueza, M. Bonfand, S. Bontemps, A. Gusdorf, A. Ginsburg, T. Csengeri, S. D. Reyes, J. Salinas, T. Baug, L. Bronfman, G. Busquet, D. J. Díaz-González, M. Fernandez-Lopez, A. Guzmán, A. Koley, H. -L. Liu, F. A. Olguin, M. Valeille-Manet , et al. (1 additional authors not shown)

Abstract: The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with an spatial resolution $\sim$0.02 pc. G353.41, at a distance of $\sim$2 kpc, has a mass of $\sim$2500 M$_{\odot}$ within $1.3\times1.3$ pc$^2$. We extract… ▽ More The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with an spatial resolution $\sim$0.02 pc. G353.41, at a distance of $\sim$2 kpc, has a mass of $\sim$2500 M$_{\odot}$ within $1.3\times1.3$ pc$^2$. We extract the N$_2$H$^+$ isolated line component and we decompose it by fitting up to 3 Gaussian velocity components. This allows us to identify velocity structures that are impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients on large ($\sim$1 pc) and small scales ($\sim$0.2 pc). We find good agreement between the N$_2$H$^+$ velocities and the previously reported DCN core velocities, suggesting that cores are kinematically coupled to the dense gas in which they form. We measure 9 converging ``V-shaped'' velocity gradients ($\sim20$ km/s/pc), located in filaments, which are sometimes associated with cores near their point of convergence. The average inflow timescale is $\sim67$ kyr, or about twice the free-fall time of cores in the same area ($\sim33$ kyr) but substantially shorter than protostar lifetime estimates ($\sim$0.5 Myr). We derive mass accretion rates in the range of $(0.35-8.77)\,\times\,10^{-4}$ M$_{\odot}$/yr. This feeding might lead to further filament collapse and formation of new cores. We suggest that the protocluster is collapsing on large scales, but the velocity signature of collapse is slow compared to pure free-fall. These data are consistent with a comparatively slow global protocluster contraction under gravity, and faster core formation within, suggesting the formation of multiple generations of stars over the protocluster lifetime. △ Less

Submitted 13 June, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Accepted in A&A, 22 pages, 21 figures, 1 interactive figure, 4 tables

Journal ref: A&A 689, A74 (2024)

ALMA-IMF XI: The sample of hot core candidates A rich population of young high-mass proto-stars unveiled by the emission of methyl formate

Authors: M. Bonfand, T. Csengeri, S. Bontemps, N. Brouillet, F. Motte, F. Louvet, A. Ginsburg, N. Cunningham, R. Galván-Madrid, F. Herpin, F. Wyrowski, M. Valeille-Manet, A. M. Stutz, J. Di Francesco, A. Gusdorf, M. Fernández-López, B. Lefloch, H-L. Liu, P. Sanhueza, R. H. Álvarez-Gutiérrez, F. Olguin, T. Nony, A. Lopez-Sepulcre, P. Dell'Ova, Y. Pouteau , et al. (6 additional authors not shown)

Abstract: Sites associated with high-mass star and cluster formation exhibit a so-called hot core phase, characterized by high temperatures and column densities of complex organic molecules. We built a comprehensive census of hot core candidates towards the ALMA-IMF protoclusters based on the detection of two CH3OCHO emission lines at 216.1 GHz. We used the source extraction algorithm GExt2D to identify pea… ▽ More Sites associated with high-mass star and cluster formation exhibit a so-called hot core phase, characterized by high temperatures and column densities of complex organic molecules. We built a comprehensive census of hot core candidates towards the ALMA-IMF protoclusters based on the detection of two CH3OCHO emission lines at 216.1 GHz. We used the source extraction algorithm GExt2D to identify peaks of methyl formate (CH3OCHO) emission that is a complex species commonly observed towards sites of star formation. We built up a catalog of 76 hot core candidates with masses ranging from about 0.2 to 80 Msun , of which 56 are new detections. A large majority of these objects are compact, rather circular, with deconvolved FWHM sizes of about 2300 au on average. About 30% of our sample of methyl formate sources have core masses above 8 Msun within sizes ranging from about 1000 au to 13400 au, which well correspond to archetypical hot cores. The origin of the CH3OCHO emission toward the lower-mass cores can be explained by a mixture of contribution from shocks, or may correspond to objects in a more evolved state, i.e. beyond the hot core stage. We find that the fraction of hot core candidates increases with the core mass. The large fraction of hot core candidates towards the most massive cores suggests that they rapidly enter the hot core phase and feedback effects from the forming protostar(s) impact their environment on short time-scales. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Journal ref: A&A 687, A163 (2024)

Turbulent vortex with moderate dust settling probed by scattering-induced polarization in the IRS 48 system

Authors: Haifeng Yang, Manuel Fernández-López, Zhi-Yun Li, Ian W. Stephens, Leslie W. Looney, Zhe-Yu Daniel Lin, Rachel Harrison

Abstract: We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $μ$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale h… ▽ More We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $μ$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale heights in this system. We model the polarization observations, focusing on the effects of dust scale heights. We find that the interplay between the inclination-induced polarization and the polarization arising from radiation anisotropy in the crescent determines the observed polarization; the anisotropy is controlled by the dust optical depth along the midplane, which is, in turn, determined by the dust scale height in the vertical direction. We find that the dust grains can neither be completely settled nor well mixed with the gas. The completely settled case produces little radial displacement between the total and polarized flux, while the well-mixed case produces an azimuthal pattern in the outer (radial) edge of the crescent that is not observed. Our best model has a gas-to-dust scale height ratio of 2, and can reproduce both the radial displacement and the azimuthal displacement between the total and polarized flux. We infer an effective turbulence $α$ parameter of approximately $0.0001-0.005$. The scattering-induced polarization provides insight into a turbulent vortex with a moderate level of dust settling in the IRS 48 system, which is hard to achieve otherwise. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 9 pages in the main text, 5 figures, accepted by ApJ

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

ALMA-IMF X -- The core population in the evolved W33-Main (G012.80) protocluster

Authors: M. Armante, A. Gusdorf, F. Louvet, F. Motte, Y. Pouteau, P. Lesaffre, R. Galván-Madrid, P. Dell'Ova, M. Bonfand, T. Nony, N. Brouillet, N. Cunningham, A. Ginsburg, A. Men'shchikov, S. Bontemps, D. Díaz González, T. Csengeri, M. Fernández-López, M. González, F. Herpin, H. -L. Liu, P. Sanhueza, A. M. Stutz, M. Valeille-Manet

Abstract: We aimed to measure the CMF in the evolved W33-Main star-forming protocluster to compare it with CMF recently obtained in other Galactic star-forming regions, including the ones included in the ALMA-IMF program. We used observations from the ALMA-IMF large program: 2'x2' maps of emission from the continuum and selected lines at 1.3mm and 3mm observed by the ALMA 12m only antennas. Our angular reso… ▽ More We aimed to measure the CMF in the evolved W33-Main star-forming protocluster to compare it with CMF recently obtained in other Galactic star-forming regions, including the ones included in the ALMA-IMF program. We used observations from the ALMA-IMF large program: 2'x2' maps of emission from the continuum and selected lines at 1.3mm and 3mm observed by the ALMA 12m only antennas. Our angular resolution was typically 1'', that is 2400au at a distance of 2.4kpc. The lines we analysed are CO(2-1), SiO(5-4), N2H+(1-0), H41alpha as well as He41alpha blended with C41alpha. We built a census of dense cores in the region, and we measured the associated CMF based on a core-dependent temperature value. We confirmed the 'evolved' status of W33-Main by identifiying three HII regions within the field, and to a lesser extent based on the number and extension of N2H+ filaments. We produced a filtered core catalog of 94 candidates, that we refined to take into account the contamination of the continuum by free-free and line emission, obtaining 80 cores with masses that range from 0.03 to 13.2Msun. We fitted the resulting high-mass end of the CMF with a single power law of the form N(log(M)) ~ M^alpha, obtaining alpha = -1.44(+0.16)(-0.22), slightly steeper but consistent with the Salpeter index. We categorized our cores in pre- and protostellar, mostly based on outlow activity and hot core nature. We found the prestellar CMF to be steeper than a Salpeter-like distribution, and the protostellar CMF to be slightly top heavy. We found a higher proportion of cores within the HII regions and their surroundings than in the rest of the field. We also found that the cores' masses were rather low (maximum mass of 13Msun). △ Less

Submitted 17 January, 2024; originally announced January 2024.

Journal ref: A&A 686, A122 (2024)

LLAMA Millimeter and Submillimeter Observatory. Update on its Science Opportunities

Authors: M. Fernandez-Lopez, P. Benaglia, S. Cichowolski, F. S. Correra, G. Cristiani, T. P. Dominici, N. Duronea, G. Gimenez de Castro, J. R. D. Lepine, I. F. Mirabel, J. P. Raulin, H. Saldano, L. Suad, C. Valotto

Abstract: The Large Latin American Millimeter Array (LLAMA for short) is a joint scientific and technological undertaking of Argentina and Brazil whose goal is to install and to operate an observing facility capable of performing observations of the Universe at millimeter and sub-millimeter wavelengths. It will consist of a 12m ALMA-like antenna with the addition of two Nasmyth cabins. LLAMA is located at 4… ▽ More The Large Latin American Millimeter Array (LLAMA for short) is a joint scientific and technological undertaking of Argentina and Brazil whose goal is to install and to operate an observing facility capable of performing observations of the Universe at millimeter and sub-millimeter wavelengths. It will consist of a 12m ALMA-like antenna with the addition of two Nasmyth cabins. LLAMA is located at 4850m above sea level in the Puna Saltenia, in the northwest region of Argentina. When completed, LLAMA will be equipped with six ALMA receivers covering Bands 1, 2+3, 5, 6, 7, and 9, which will populate the two Nasmyth cabins. We summarize here the main ideas related with the Science that LLAMA could accomplish on different astronomical topics, gathered from the experience of a group of international experts on each field. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 11 pages, contributed paper to the workshop "Prospects for low-frequency radio astronomy in South America", held in Buenos Aires, Nov 2022; to appear in the RevMexAA-SC

Multiple shells driven by disk winds: ALMA observations in the HH 30 outflow

Authors: J. A. López-Vázquez, Chin-Fei Lee, M. Fernández-López, Fabien Louvet, O. Guerra-Alvarado, Luis A. Zapata

Abstract: We present archive Atacama Large Millimeter/Submillimeter Array (ALMA) Band 6 observations of the $^{13}$CO (J=2-1) and $^{12}$CO (J=2-1) molecular line emission of the protostellar system associated with HH 30. The $^{13}$CO molecular line shows the accretion disk while the molecular outflow is traced by the emission of the $^{12}$CO molecular line. We estimated a dynamical mass for the central o… ▽ More We present archive Atacama Large Millimeter/Submillimeter Array (ALMA) Band 6 observations of the $^{13}$CO (J=2-1) and $^{12}$CO (J=2-1) molecular line emission of the protostellar system associated with HH 30. The $^{13}$CO molecular line shows the accretion disk while the molecular outflow is traced by the emission of the $^{12}$CO molecular line. We estimated a dynamical mass for the central object of $0.45\pm0.14$ M$_\odot$, and a mass for the molecular outflow of $1.83\pm0.19\times10^{-4}$ M$_\odot$. The molecular outflow presents an internal cavity as well as multiple outflowing shell structures. We distinguish three different shells with constant expansion ($\sim4-6$ km s$^{-1}$) and possible rotation signatures ($\leq0.5$ km s$^{-1}$). We find that the shells can be explained by magnetocentrifugal disk winds with launching radii $R_\mathrm{launch}\lesssim4$ au and a small magnetic lever arm $λ\sim1.6-1.9$. The multiple shell structure may be the result of episodic ejections of the material from the accretion disk associated with three different epochs with dynamical ages of $497\pm15$ yr, $310\pm9$ yr, and $262\pm11$ yr for the first, second, and third shells, respectively. The outermost shell was ejected $187\pm17$ yr before the medium shell, while the medium shell was launched $48\pm14$ yr before the innermost shell. Our estimations of the linear and angular momentum rates of the outflow as well as the accretion luminosity are consistent with the expected values if the outflow of HH 30 is produced by a wide-angle disk wind. △ Less

Submitted 8 December, 2023; v1 submitted 5 December, 2023; originally announced December 2023.

Aligned Grains and Scattered Light Found in Gaps of Planet-Forming Disk

Authors: Ian W. Stephens, Zhe-Yu Daniel Lin, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Haifeng Yang, Rachel Harrison, Akimasa Kataoka, Carlos Carrasco-Gonzalez, Satoshi Okuzumi, Ryo Tazaki

Abstract: Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields a… ▽ More Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $μ$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Published in Nature

ALMA-IMF VIII -- Combination of Interferometric Continuum Images with Single-Dish Surveys and Structural Analysis of Six Protoclusters

Authors: Daniel J. Díaz-González, Roberto Galván-Madrid, Adam Ginsburg, Frédérique Motte, Pierre Dell'Ova, Stan Kurtz, Nichol Cunningham, Amelia M. Stutz, Fabien Louvet, Timea Csengeri, Manuel Fernández-López, Patricio Sanhueza, Thomas Nony, Rudy Rivera-Soto, Rodrigo H. Álvarez-Gutiérrez, Melanie Armante, Melisse Bonfand, Sylvain Bontemps, Antoine Gusdorf, Hong-Li Liu

Abstract: We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the… ▽ More We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the combined images, and performed further analysis in the six fields with combination in both bands (G012.80, W43-MM1, W43-MM2, W43-MM3, W51-E, W51-IRS2). In these fields, we calculated spectral index maps and used them to separate regions dominated by dust or free-free emission, and then performed further structural analysis. We report the basic physical parameters of the dust-dominated (column densities, masses) and ionized (emission measures, hydrogen ionization photon rates) structures. We also searched for multi-scale relations in the dust-dominated structures across the analyzed fields, finding that the fraction of mass in dendrogram leaves (which we label as "Leaf Mass Eficiency", LME) as a function of molecular gas column density follows a similar trend: a rapid, exponential-like growth, with maximum values approaching 100% in most cases. The observed behaviour of the LME with gas column is tentatively interpreted as an indicator of large star formation activity within the ALMA-IMF protoclusters. W51-E and G012.80 stand out as cases with comparatively large and reduced potential for further star formation, respectively. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: Accepted to The Astrophysical Journal Supplement

ALMA-IMF IX: Catalog and Physical Properties of 315 SiO Outflow Candidates in 15 Massive Protoclusters

Authors: A. P. M. Towner, A. Ginsburg, P. Dell'Ova, A. Gusdorf, S. Bontemps, T. Csengeri, R. Galván-Madrid, F. K. Louvet, F. Motte, P. Sanhueza, A. M. Stutz, J. Bally, T. Baug, H. R. V. Chen, N. Cunningham, M. Fernández-López, H. -L. Liu, X. Lu, T. Nony, M. Valeille-Manet, B. Wu, R. H. Álvarez-Gutiérrez, M. Bonfand, J. Di Francesco, Q. Nguyen-Luong , et al. (2 additional authors not shown)

Abstract: We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,00… ▽ More We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,000 years, yielding median mass, momentum, and energy rates of $\dot{M}$ = 10$^{-4.4}$ M$_{\odot}$ yr$^{-1}$, $\dot{P}$ = 10$^{-3.2}$ M$_{\odot}$ km/s yr$^{-1}$, and $\dot{E}$ = 1 L$_{\odot}$. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (~3$-$5$σ$), and no correlations above 3$σ$ with clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass - which has been previously described in the literature - may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15-30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force-bolometric luminosity relationship in the literature, and our data extend this relationship up to L $\geq$ 10$^6$ L$_{\odot}$ and $\dot{P}$ $\geq$ 1 M$_{\odot}$ km/s yr$^{-1}$. Our lack of correlation with clump L/M is inconsistent with models of protocluster formation in which all protostars start forming at the same time. △ Less

Submitted 27 October, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 46 pages, 14 figures, 10 tables. This publication has an associated Zenodo entry, which can be found here: https://zenodo.org/records/8350595

One, Two, Three ... An Explosive Outflow in IRAS 12326$-$6245 revealed by ALMA

Authors: Luis A. Zapata, Manuel Fernández-López, Silvia Leurini, Estrella Guzmán Ccolque, Luis F. Rodriguez, Aina Palau, Karl M. Menten, Friedrich Wyrowski

Abstract: In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using… ▽ More In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using Band 6 sensitive and high angular resolution ($\sim$0.2$"$) Atacama Large Millimeter/Submillimeter Array (ALMA) observations. We find over 10 molecular and collimated well-defined streamers, with Hubble-Lemaitre like expansion motions, and pointing right to the center of a dusty and molecular shell (reported for the first time here) localized in the northern part of the UCHII region known as G301.1A. The estimated kinematic age, and energy for the explosion are $\sim$700 yrs, and 10$^{48}$ erg, respectively. Taking into account the recently reported explosive outflows together with IRAS 12326$-$6245, we estimate an event rate of once every 90 yr in our Galaxy, similar to the formation rate of massive stars. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Accepted to The Astrophysical Journal Letters

Panchromatic (Sub)millimeter Polarization Observations of HL Tau Unveil Aligned Scattering Grains

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Ian W. Stephens, Manuel Fernández-López, Carlos Carrasco-González, Claire J. Chandler, Alice Pasetto, Leslie W. Looney, Haifeng Yang, Rachel E. Harrison, Sarah I. Sadavoy, Thomas Henning, A. Meredith Hughes, Akimasa Kataoka, Woojin Kwon, Takayuki Muto, Dominique Segura-Cox

Abstract: Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $λ$. The polarization m… ▽ More Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $λ$. The polarization morphology at Band 7 is predominantly parallel to the disk minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ~0.2" (29 au) polarization observations at Q-Band (7.0 mm) using the Karl G. Jansky Very Large Array (VLA) and at Bands 4 (2.1 mm), 5 (1.5 mm), and 7 using ALMA, consolidating HL Tau's position as the protoplanetary disk with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 continue to follow the morphological transition with wavelength previously identified in Bands 3, 6, and 7. Based on the azimuthal variation, we decompose the polarization into contributions from scattering ($s$) and thermal emission ($t$). We find that $s$ decreases slowly with increasing $λ$, and $t$ increases more rapidly with $λ$ which are expected from optical depth effects of toroidally aligned, scattering prolate grains. The relatively weak $λ$ dependence of $s$ is consistent with large, porous grains. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 22 pages, 13 figures, submitted to MNRAS

Disk and Envelope Streamers of the GGD27-MM1 Massive Protostar

Authors: M. Fernández-López, J. M. Girart, J. A. López-Vázquez, R. Estalella, G. Busquet, S. Curiel, N. Añez-López

Abstract: We present new Atacama Large (sub)Millimeter Array 0.98 mm observations of the continuum emission and several molecular lines toward the high-mass protostellar system GGD27-MM1, driving the HH 80-81 radio-jet. The detailed analysis of the continuum and the CH$_3$CN molecular emission allows us to separate the contributions from the dust content of the disk (extending up to 190 au), the molecular c… ▽ More We present new Atacama Large (sub)Millimeter Array 0.98 mm observations of the continuum emission and several molecular lines toward the high-mass protostellar system GGD27-MM1, driving the HH 80-81 radio-jet. The detailed analysis of the continuum and the CH$_3$CN molecular emission allows us to separate the contributions from the dust content of the disk (extending up to 190 au), the molecular content of the disk (extending from 140 to 360 au), and the content of the envelope, revealing the presence of several possible accretion streamers (also seen in other molecular tracers, such as CH$_3$OH). We analyze the physical properties of the system, producing temperature and column density maps, and radial profiles for the disk and the envelope. We qualitatively reproduce the trajectories and line-of-sight velocities of the possible streamers using a theoretical model approach. An ad-hoc model of a flared disk comprising a hot dust disk embedded in cold gas fits the H$_2$S emission, which revealed the molecular disk as crescent-shape with a prominent central absorption. Another fit to the central absorption spectrum suggests that the absorption is probably caused by different external cold layers from the envelope or the accretion streamers. Finally, the analysis of the rotation pattern of the different molecular transitions in the molecular disk, suggests that there is an inner zone devoid of molecular content. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: 36 pages, 28 figures, 6 tables. Accepted by the ApJ

ALMA-IMF. VII. First release of the full spectral line cubes: Core kinematics traced by DCN J=(3-2)

Authors: N. Cunningham, A. Ginsburg, R. Galván-Madrid, F. Motte, T. Csengeri, A. M. Stutz, M. Fernández-López, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, S. Bontemps, J. Braine, N. Brouillet, G. Busquet, D. J. Díaz-González, J. Di Francesco, A. Gusdorf, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, X. Lu, L. Maud, T. Nony , et al. (8 additional authors not shown)

Abstract: ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), prod… ▽ More ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), produced from the combination of two ALMA 12m-array configurations. The data include 12 spectral windows, with eight at 1.3mm and four at 3mm. The broad spectral coverage of ALMA-IMF (~6.7 GHz bandwidth coverage per field) hosts a wealth of simple atomic, molecular, ionised, and complex organic molecular lines. We describe the line cube calibration done by ALMA and the subsequent calibration and imaging we performed. We discuss our choice of calibration parameters and optimisation of the cleaning parameters, and we demonstrate the utility and necessity of additional processing compared to the ALMA archive pipeline. As a demonstration of the scientific potential of these data, we present a first analysis of the DCN (3-2) line. We find that DCN traces a diversity of morphologies and complex velocity structures, which tend to be more filamentary and widespread in evolved regions and are more compact in the young and intermediate-stage protoclusters. Furthermore, we used the DCN (3-2) emission as a tracer of the gas associated with 595 continuum cores across the 15 protoclusters, providing the first estimates of the core systemic velocities and linewidths within the sample. We find that DCN (3-2) is detected towards a higher percentage of cores in evolved regions than the young and intermediate-stage protoclusters and is likely a more complete tracer of the core population in more evolved protoclusters. The full ALMA 12m-array cubes for the ALMA-IMF Large Program are provided with this DR1 release. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 75 pages (21 main body; 54 appendix), 37 figures. The ALMA-IMF DR1 line release is hosted at https://dataverse.harvard.edu/dataverse/alma-imf

Journal ref: A&A 678, A194 (2023)

Eccentric Dust Ring in the IRS 48 Transition Disk

Authors: Haifeng Yang, Manuel Fernandez-Lopez, Zhi-Yun Li, Ian W. Stephens, Leslie W. Looney, Zhe-Yu Daniel Lin, Rachel Harrison

Abstract: Crescent-shaped structures in transition disks hold the key to studying the putative companions to the central stars. The dust dynamics, especially that of different grain sizes, is important to understanding the role of pressure bumps in planet formation. In this work, we present deep dust continuum observation with high resolution towards the Oph IRS 48 system. For the first time, we are able to… ▽ More Crescent-shaped structures in transition disks hold the key to studying the putative companions to the central stars. The dust dynamics, especially that of different grain sizes, is important to understanding the role of pressure bumps in planet formation. In this work, we present deep dust continuum observation with high resolution towards the Oph IRS 48 system. For the first time, we are able to significantly trace and detect emission along $95\%$ of the ring crossing the crescent-shaped structure. The ring is highly eccentric with an eccentricity of $0.27$. The flux density contrast between the peak of the flux and its counter part along the ring is $\sim 270$. In addition, we detect a compact emission toward the central star. If the emission is an inner circumstellar disk inside the cavity, it has a radius of at most a couple of astronomical units with a dust mass of $1.5\times 10^{-8}\rm\, M_\odot$, or $0.005\rm\, M_\oplus$. We also discuss the implications of the potential eccentric orbit on the proper motion of the crescent, the putative secondary companion, and the asymmetry in velocity maps. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 9 pages, 5 figures, accepted for publication on AAS Journals

ALMA-IMF. V. Prestellar and protostellar core populations in the W43 cloud complex

Authors: T. Nony, R. Galvan-Madrid, F. Motte, Y. Pouteau, N. Cunningham, F. Louvet, A. M. Stutz, B. Lefloch, S. Bontemps, N. Brouillet, A. Ginsburg, I. Joncour, F. Herpin, P. Sanhueza, T. Csengeri, A. P. M. Towner, M. Bonfand, M. Fernández-López, T. Baug, L. Bronfman, G. Busquet, J. Di Francesco, A. Gusdorf, X. Lu, F. Olguin , et al. (2 additional authors not shown)

Abstract: The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and… ▽ More The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $α=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $α=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores. △ Less

Submitted 15 March, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

Journal ref: A&A 674, A75 (2023)

ALMA-IMF VI -- Investigating the origin of stellar masses: Core mass function evolution in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, M. Gonzalez, I. Joncour, J. -F. Robitaille, G. Busquet, R. Galvan-Madrid, A. Gusdorf, P. Hennebelle, A. Ginsburg, T. Csengeri, P. Sanhueza, P. Dell'Ova, A. M. Stutz, A. P. M. Towner, N. Cunningham, F. Louvet, A. Men'shchikov, M. Fernandez-Lopez, N. Schneider, M. Armante, J. Bally, T. Baug, M. Bonfand , et al. (13 additional authors not shown)

Abstract: Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present… ▽ More Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present study focuses on the W43-MM2&MM3 mini-starburst, whose CMF has recently been found to be top-heavy with respect to the Salpeter slope. W43-MM2&MM3 harbors a rich cluster that contains a statistically significant number of cores, which was previously characterized in Paper III. We applied a multi-scale decomposition technique to the ALMA 1.3 mm and 3 mm continuum images to define six subregions. For each subregion we characterized the high column density probability distribution function, n-PDF, and the shape of the cloud gas using the 1.3 mm image. Using the core catalog, we investigate correlations between the CMF and cloud and core properties. We classify the subregions into different stages of evolution, from quiescent to burst to post-burst, based on the surface number density of cores, number of outflows, and UCHii presence. The high-mass end of the subregion CMFs varies from being close to the Salpeter slope (quiescent) to top-heavy (burst and post-burst). Moreover, the second tail of the n-PDF varies from steep, to flat like observed for the high mass star-forming clouds. We found that subregions with flat second n-PDF tails display top-heavy CMFs. The CMF may evolve from Salpeter to top-heavy throughout the star formation process from the quiescent to the burst phase. This scenario raises the question of if the CMF might revert again to Salpeter as the cloud approaches the end of its star formation stage, a hypothesis that remains to be tested. △ Less

Submitted 20 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A (december, 5th 2022) after language editing; 25 pages, 14 figures

Journal ref: A&A 674, A76 (2023)

Possible Explosive Dispersal Outflow in IRAS 16076-5134 revealed with ALMA

Authors: Estrella Guzmán Ccolque, Manuel Fernández-López, Luis A. Zapata, Tapas Baug

Abstract: We present 0.9 mm continuum and CO (3-2) line emission observations retrieved from the Atacama Large Millimeter/submillimeter Array (ALMA) archive toward the high-mass star formation region IRAS 16076-5134. We identify fourteen dense cores with masses between 0.3 to 22 M$_{\odot}$. We find an ensemble of filament-like CO (3-2) ejections from -62 to +83 km s$^{-1}$ that appear to arise radially fro… ▽ More We present 0.9 mm continuum and CO (3-2) line emission observations retrieved from the Atacama Large Millimeter/submillimeter Array (ALMA) archive toward the high-mass star formation region IRAS 16076-5134. We identify fourteen dense cores with masses between 0.3 to 22 M$_{\odot}$. We find an ensemble of filament-like CO (3-2) ejections from -62 to +83 km s$^{-1}$ that appear to arise radially from a common central position, close to the dense core MM8. The ensemble of filaments, has a quasi-isotropic distribution in the plane of the sky. The radial velocity of several filaments follow a linear velocity gradient, incresing from a common origin. Considering the whole ensemble of filaments, we estimate its total mass to be 138 and 216 M$_{\odot}$ from its CO emission, for 70 K and 140 K respectively. Also, assuming a constant velocity expansion of the filaments (of 83 km s$^{-1}$) we estimate the dynamical age of the outflowing material (3500 years), its momentum (~10$^{4}$ M$_{\odot}$ km s$^{-1}$) and its kinetic energy (~10$^{48-49}$ erg). The morphology and kinematics presented by the filaments suggest the presence of a dispersal outflow with explosive characteristics in IRAS 16076-5134. In addition, we make a raw estimate of the lower limit of the frequency rate of the explosive dispersal outflows in the Galaxy (one every 110 years) considering constant star formation rate and efficiency with respect to the galactocentric radius of the Galaxy. This may imply a comparable rate of dispersal outflows and supernovae (approximately one every 50 years), which may be important for the energy budget of the Interstellar Medium and the link between dispersal outflows and high-mass star formation. △ Less

Submitted 25 August, 2022; originally announced August 2022.

Comments: 13 pages, 8 figures. ApJ accepted

ALMA-IMF IV -- A comparative study of the main hot cores in W43-MM1: detection, temperature and molecular composition

Authors: N. Brouillet, D. Despois, J. Molet, T. Nony, F. Motte, A. Gusdorf, F. Louvet, S. Bontemps, F. Herpin, M. Bonfand, T. Csengeri, A. Ginsburg, N. Cunningham, R. Galvan-Madrid, L. Maud, G. Busquet, L. Bronfman, M. Fernandez-Lopez, D. L. Jeff, B. Lefloch, Y. Pouteau, P. Sanhueza, A. M. Stutz, M. Valeille-Manet

Abstract: W43-MM1 is a young region, very rich in terms of high-mass star formation. We aim to systematically identify the massive cores which contain a hot core and compare their molecular composition. We used ALMA high-spatial resolution (2500 au) data of W43-MM1 to identify line-rich protostellar cores and make a comparative study of their temperature and molecular composition. The identification of hot… ▽ More W43-MM1 is a young region, very rich in terms of high-mass star formation. We aim to systematically identify the massive cores which contain a hot core and compare their molecular composition. We used ALMA high-spatial resolution (2500 au) data of W43-MM1 to identify line-rich protostellar cores and make a comparative study of their temperature and molecular composition. The identification of hot cores is based on both the spatial distribution of the complex organic molecules and the contribution of molecular lines relative to the continuum intensity. We rely on the analysis of CH3CN and CH3CCH to estimate the temperatures of the selected cores. Finally, we rescale the spectra of the different hot cores based on their CH3OCHO line intensities to directly compare the detections and line intensities of the other species. W43-MM1 turns out to be a region rich in massive hot cores with at least 1 less massive and 7 massive hot cores. The excitation temperature of CH3CN is of the same order for all of them (120-160 K). There is a factor of up to 30 difference in the intensity of the complex organic molecules (COMs) lines. However the molecular emission of the hot cores appears to be the same within a factor 2-3. This points towards both a similar chemical composition and excitation of most of the COMs over these massive cores, which span about an order of magnitude in core mass. In contrast, CH3CCH emission is found to preferentially trace more the envelope, with a temperature ranging from 50 K to 90 K. Core 1, the most massive hot core of W43-MM1, shows a richer line spectrum than the other cores. In core 2, the emission of O-bearing molecules does not peak at the dust continuum core center; the blue and red shifted emission correspond to the outflow lobes, suggesting a formation via the sublimation of the ice mantles through shocks or UV irradiation on the walls of the cavity. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 19 pages, 17 figures. Accepted for publication in A&A

Journal ref: A&A 665, A140 (2022)

(Sub)millimeter Dust Polarization of Protoplanetary Disks from Scattering by Large Millimeter-Sized Irregular Grains

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Haifeng Yang, Olga Muñoz, Leslie Looney, Ian Stephens, Charles L. H. Hull, Manuel Fernández-López, Rachel Harrison

Abstract: The size of dust grains, $a$, is key to the physical and chemical processes in circumstellar disks, but observational constraints of grain size remain challenging. (Sub)millimeter continuum observations often show a percent-level polarization parallel to the disk minor axis, which is generally attributed to scattering by $\sim 100μ$m-sized spherical grains (with a size parameter… ▽ More The size of dust grains, $a$, is key to the physical and chemical processes in circumstellar disks, but observational constraints of grain size remain challenging. (Sub)millimeter continuum observations often show a percent-level polarization parallel to the disk minor axis, which is generally attributed to scattering by $\sim 100μ$m-sized spherical grains (with a size parameter $x \equiv 2πa / λ< 1$, where $λ$ is the wavelength). Larger spherical grains (with $x$ greater than unity) would produce opposite polarization direction. However, the inferred size is in tension with the opacity index $β$ that points to larger mm/cm-sized grains. We investigate the scattering-produced polarization by large irregular grains with a range of $x$ greater than unity with optical properties obtained from laboratory experiments. Using the radiation transfer code, RADMC-3D, we find that large irregular grains still produce polarization parallel to the disk minor axis. If the original forsterite refractive index in the optical is adopted, then all samples can produce the typically observed level of polarization. Accounting for the more commonly adopted refractive index using the DSHARP dust model, only grains with $x$ of several (corresponding to $\sim$mm-sized grains) can reach the same polarization level. Our results suggest that grains in disks can have sizes in the millimeter regime, which may alleviate the tension between the grain sizes inferred from scattering and other means. Additionally, if large irregular grains are not settled to the midplane, their strong forward scattering can produce asymmetries between the near and far side of an inclined disk, which can be used to infer their presence. △ Less

Submitted 6 March, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 14 pages, 9 figures, published by MNRAS

Catching the Butterfly and the Homunculus of $η$ Carinae with ALMA

Authors: Luis A. Zapata, Laurent Loinard, Manuel Fernández-López, Jesús A. Toalá, Ricardo F. González, Luis F. Rodrí guez, Theodore R. Gull, Patrick W. Morris, Karl M. Menten, Tomasz Kamiński

Abstract: The nature and origin of the molecular gas component located in the circumstellar vicinity of $η$ Carinae are still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) CO(3$-$2) observations with a high angular resolution ($\sim$0.15$''$), and a great sensitivity that are employed to reveal the origin of this component in $η$ Carinae. These ob… ▽ More The nature and origin of the molecular gas component located in the circumstellar vicinity of $η$ Carinae are still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) CO(3$-$2) observations with a high angular resolution ($\sim$0.15$''$), and a great sensitivity that are employed to reveal the origin of this component in $η$ Carinae. These observations reveal much higher velocity ($-$300 to $+$270 km s$^{-1}$) blue and redshifted molecular thermal emission than previously reported, which we associate with the lobes of the Homunculus Nebula, and that delineates very well the innermost contours of the red- and blue-shifted lobes likely due by limb brightening. The inner contour of the redshifted emission was proposed to be a {\it disrupted torus}, but here we revealed that it is at least part of the molecular emission originated from the lobes and/or the expanding equatorial skirt. On the other hand, closer to systemic velocities ($\pm$100 km s$^{-1}$), the CO molecular gas traces an inner butterfly-shaped structure that is also revealed at NIR and MIR wavelengths as the region in which the shielded dust resides. The location and kinematics of the molecular component indicate that this material has formed after the different eruptions of $η$ Carinae. △ Less

Submitted 26 May, 2022; originally announced May 2022.

Comments: Accepted to Astrophysical Journal. 3D animation -> https://drive.google.com/file/d/1MRWaTkij-Lf_uzTfamJfHeGvgzjDe6Oq/view?usp=sharing

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software p… ▽ More The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

Journal ref: A&A 664, A26 (2022)

Thermal Emission and Scattering by Aligned Grains: Plane-Parallel Model and Application to Multiwavelength Polarization of the HL Tau Disk

Authors: Zhe-Yu Daniel Lin, Zhi-Yun Li, Haifeng Yang, Ian Stephens, Leslie Looney, Rachel Harrison, Manuel Fernández-López

Abstract: Telescopes are now able to resolve dust polarization across circumstellar disks at multiple wavelengths, allowing the study of the polarization spectrum. Most disks show clear evidence of dust scattering through their unidirectional polarization pattern typically at the shorter wavelength of $\sim 870 μ$m. However, certain disks show an elliptical pattern at $\sim 3$mm, which is likely due to alig… ▽ More Telescopes are now able to resolve dust polarization across circumstellar disks at multiple wavelengths, allowing the study of the polarization spectrum. Most disks show clear evidence of dust scattering through their unidirectional polarization pattern typically at the shorter wavelength of $\sim 870 μ$m. However, certain disks show an elliptical pattern at $\sim 3$mm, which is likely due to aligned grains. With HL Tau, its polarization pattern at $\sim 1.3$mm shows a transition between the two patterns making it the first example to reveal such transition. We use the T-matrix method to model elongated dust grains and properly treat scattering of aligned non-spherical grains with a plane-parallel slab model. We demonstrate that a change in optical depth can naturally explain the polarization transition of HL Tau. At low optical depths, the thermal polarization dominates, while at high optical depths, dichroic extinction effectively takes out the thermal polarization and scattering polarization dominates. Motivated by results from the plane-parallel slab, we develop a simple technique to disentangle thermal polarization of the aligned grains $T_{0}$ and polarization due to scattering $S$ using the azimuthal variation of the polarization fraction. We find that, with increasing wavelength, the fractional polarization spectrum of the scattering component $S$ decreases, while the thermal component $T_{0}$ increases, which is expected since the optical depth decreases. We find several other sources similar to HL Tau that can be explained by azimuthally aligned scattering prolate grains when including optical depth effects. In addition, we explore how spirally aligned grains with scattering can appear in polarization images. △ Less

Submitted 25 April, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 27 pages, 21 figures, accepted by MNRAS

ALMA-IMF II -- investigating the origin of stellar masses: Continuum Images and Data Processing

Authors: A. Ginsburg, T. Csengeri, R. Galván-Madrid, N. Cunningham, R. H. Álvarez-Gutiérrez, T. Baug, M. Bonfand, S. Bontemps, G. Busquet, D. J. Díaz-González, M. Fernández-López, A. Guzmán, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, L. Maud, F. Motte, F. Nakamura, T. Nony, F. A. Olguin, Y. Pouteau, P. Sanhueza, A. M. Stutz, A. P. M. Towner , et al. (27 additional authors not shown)

Abstract: We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Ar… ▽ More We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Array (ALMA) observatory and the subsequent calibration and imaging we performed. The image products are combinations of multiple 12m array configurations created from a selection of the observed bandwidth using multi-term, multi-frequency synthesis imaging and deconvolution. The data products are self-calibrated and exhibit substantial noise improvements over the images produced from the delivered data. We compare different choices of continuum selection, calibration parameters, and image weighting parameters, demonstrating the utility and necessity of our additional processing work. Two variants of continuum selection are used and will be distributed: the "best-sensitivity" data, which include the full bandwidth, including bright emission lines that contaminate the continuum, and "cleanest", which select portions of the spectrum that are unaffected by line emission. We present a preliminary analysis of the spectral indices of the continuum data, showing that the ALMA products are able to clearly distinguish free-free emission from dust emission, and that in some cases we are able to identify optically thick emission sources. The data products are made public with this release. △ Less

Submitted 13 May, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: Data released on Zenodo at https://zenodo.org/record/5702966. Accepted for publication in A&A

Journal ref: A&A 662, A9 (2022)

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

Authors: F. Motte, S. Bontemps, T. Csengeri, Y. Pouteau, F. Louvet, A. M. Stutz, N. Cunningham, A. López-Sepulcre, N. Brouillet, R. Galván-Madrid, A. Ginsburg, L. Maud, A. Men'shchikov, F. Nakamura, T. Nony, P. Sanhueza, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, G. Busquet, E. Chapillon, D. Díaz-González, M. Fernández-López, A. E. Guzmán , et al. (39 additional authors not shown)

Abstract: The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2… ▽ More The ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 27 pages, 9 figures, accepted for publication by A&A

Journal ref: A&A 662, A8 (2022)

Magnetic Fields in Massive Star-Forming Regions (MagMaR) II. Tomography Through Dust and Molecular Line Polarization in NGC 6334I(N)

Authors: Paulo C. Cortes, Patricio Sanhueza, Martin Houde, Sergio Martin, Charles L. H. Hull, Josep M. Girart, Qizhou Zhang, Manuel Fernandez-Lopez, Luis A. Zapata, Ian W. Stephens, Hua-bai Li, Benjamin Wu, Fernando Olguin, Xing Lu, Andres E. Guzman, Fumitaka Nakamura

Abstract: Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to th… ▽ More Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to the field. By considering previous findings, the field retains a pinched shape which can be traced to clump length-scales from the envelope scales traced by ALMA, suggesting that the field is dynamically important across multiple length-scales in this region. The CS total intensity emission is found to be optically thick ($τ_{\mathrm{CS}} = 32 \pm 12$) while the C$^{33}$S emission appears to be optically thin ($τ_{\mathrm{C^{33}S}} = 0.1 \pm 0.01$). This suggests that sources of anisotropy other than large velocity gradients, i.e. anisotropies in the radiation field are required to explain the polarized emission from CS seen by ALMA. By using four variants of the Davis-Chandrasekhar-Fermi technique and the angle dispersion function methods (ADF), we obtain an average of estimates for the magnetic field strength onto the plane of the sky of $\left< \mathrm{B}_{\mathrm{pos}} \right> = 16$ mG from the dust and $\left< \mathrm{B}_{\mathrm{pos}} \right> \sim 2$ mG from the CS emission, where each emission traces different molecular hydrogen number densities. This effectively enables a tomographic view of the magnetic field within a single ALMA observation. △ Less

Submitted 19 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ

Gravity Driven Magnetic Field at ~1000 au Scales in High-mass Star Formation

Authors: Patricio Sanhueza, Josep Miquel Girart, Marco Padovani, Daniele Galli, Charles L. H. Hull, Qizhou Zhang, Paulo Cortes, Ian W. Stephens, Manuel Fernandez-Lopez, James M. Jackson, Pau Frau, Patrick M. Kock, Benjamin Wu, Luis A. Zapata, Fernando Olguin, Xing Lu, Andrea Silva, Ya-Wen Tang, Takeshi Sakai, Andres E. Guzman, Ken'ichi Tatematsu, Fumitaka Nakamura, Huei-Ru Vivien Chen

Abstract: A full understanding of high-mass star formation requires the study of one of the most elusive components of the energy balance in the interstellar medium: magnetic fields. We report ALMA 1.2 mm, high-resolution (700 au) dust polarization and molecular line observations of the rotating hot molecular core embedded in the high-mass star-forming region IRAS 18089-1732. The dust continuum emission and… ▽ More A full understanding of high-mass star formation requires the study of one of the most elusive components of the energy balance in the interstellar medium: magnetic fields. We report ALMA 1.2 mm, high-resolution (700 au) dust polarization and molecular line observations of the rotating hot molecular core embedded in the high-mass star-forming region IRAS 18089-1732. The dust continuum emission and magnetic field morphology present spiral-like features resembling a whirlpool. The velocity field traced by the H13CO+ (J=3-2) transition line reveals a complex structure with spiral filaments that are likely infalling and rotating, dragging the field with them. We have modeled the magnetic field and find that the best model corresponds to a weakly magnetized core with a mass-to-magnetic-flux ratio (lambda) of 8.38. The modeled magnetic field is dominated by a poloidal component, but with an important contribution from the toroidal component that has a magnitude of 30% of the poloidal component. Using the Davis-Chandrasekhar-Fermi method, we estimate a magnetic field strength of 3.5 mG. At the spatial scales accessible to ALMA, an analysis of the energy balance of the system indicates that gravity overwhelms turbulence, rotation, and the magnetic field. We show that high-mass star formation can occur in weakly magnetized environments, with gravity taking the dominant role. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Accepted for publications in ApJL (12 pages, 5 figures, Appendix)

Magnetic Fields in Massive Star-Forming Regions (MagMaR) I. Linear Polarized Imaging of the UCHII Region G5.89-0.39

Authors: M. Fernández-López, P. Sanhueza, L. A. Zapata, I. Stephens, C. Hull, Q. Zhang, J. M. Girart, P. M. Koch, P. Cortés, A. Silva, K. Tatematsu, F. Nakamura, A. E. Guzmán, Q. Nguyen Luong, E. Guzmán Ccolque, Y. -W. Tang, V. Chen

Abstract: We report 1.2 mm polarized continuum emission observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the high-mass star formation region G5.89-0.39. The observations show a prominent 0.2 pc north-south filamentary structure. The UCHII in G5.89-0.39 breaks the filament in two pieces. Its millimeter emission shows a dusty belt with a mass of 55-115 M$_{\odot}$ a… ▽ More We report 1.2 mm polarized continuum emission observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the high-mass star formation region G5.89-0.39. The observations show a prominent 0.2 pc north-south filamentary structure. The UCHII in G5.89-0.39 breaks the filament in two pieces. Its millimeter emission shows a dusty belt with a mass of 55-115 M$_{\odot}$ and 4,500 au in radius, surrounding an inner part comprising mostly ionized gas with a dust emission only accounting about 30% of the total millimeter emission. We also found a lattice of convex arches which may be produced by dragged dust and gas from the explosive dispersal event involving the O5 Feldt's star. The north-south filament has a mass between 300-600 M$_{\odot}$ and harbours a cluster of about 20 millimeter envelopes with a median size and mass of 1700 au and 1.5 M$_{\odot}$, respectively, some of which are already forming protostars. We interpret the polarized emission in the filament as mainly coming from magnetically aligned dust grains. The polarization fraction is ~4.4% in the filaments and 2.1% at the shell. The magnetic fields are along the North Filament and perpendicular to the South Filament. In the Central Shell, the magnetic fields are roughly radial in a ring surrounding the dusty belt between 4,500 and 7,500 au, similar to the pattern recently found in the surroundings of Orion BN/KL. This may be an independent observational signpost of explosive dispersal outflows and should be further investigated in other regions. △ Less

Submitted 27 April, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted in ApJ. 26 pages, 17 figures

Outflows, envelopes, and disks as evolutionary indicators in Lupus YSOs

Authors: M. M. Vazzano, M. Fernández-López, A. Plunkett, I. de Gregorio-Monsalvo, A. Santamaría-Miranda, S. Takahashi, C. Lopez

Abstract: By studying 7 objects in the Lupus clouds we aim to test if a coherence exists between commonly used evolutionary tracers. We present ALMA observations of the continuum and molecular line emission that probe the dense gas and dust of cores and their associated molecular outflows. Our source selection in a common environment allows for a consistent comparison across different evolutionary stages.… ▽ More By studying 7 objects in the Lupus clouds we aim to test if a coherence exists between commonly used evolutionary tracers. We present ALMA observations of the continuum and molecular line emission that probe the dense gas and dust of cores and their associated molecular outflows. Our source selection in a common environment allows for a consistent comparison across different evolutionary stages. The quality of the ALMA molecular data allows us to reveal the nature of the molecular outflows by studying their morphology and kinematics. The images in IRAS15398-3359 appear to show that it drives a precessing episodic jet-driven outflow with at least 4 ejections separated by periods of time between 50 and 80 years, while data in IRAS16059-3857 show similarities with a wide-angle wind model also showing signs of being episodic. The outflow of J160115-41523 could be better explain with the wide-angle wind model as well, but new observations are needed to explore its nature. We find that the most common evolutionary tracers are useful for broad evolutionary classifications, but are not consistent with each other to provide enough granularity to disentangle different evolutionary stage of sources that belong to the same Class. Outflow properties used as protostellar age tracers (mass, momentum, energy, opening angle) may suffer from differences in the nature of each outflow, thus detailed observations are needed to refine evolutionary classifications. We found both AzTEC-lup1-2 and AzTEC-lup3-5 to be in the pre-stellar stage, although the latter could be more evolved. IRAS15398-3359, IRAS16059-3857 and J160115-41523, which have clearly detected outflows, are Class 0 sources, although we are not able to determine which is younger and which is older. Sz102 and Merin28 are the most evolved sources and show signs of having associated flows, not as well traced by CO as for the younger sources. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Journal ref: A&A 648, A41 (2021)

ALMA CN Zeeman Observations of AS 209: Limits on Magnetic Field Strength and Magnetically Driven Accretion Rate

Authors: Rachel E. Harrison, Leslie W. Looney, Ian W. Stephens, Zhi-Yun Li, Richard Teague, Richard M. Crutcher, Haifeng Yang, Erin Cox, Manuel Fernández-López, Hiroko Shinnaga

Abstract: While magnetic fields likely play an important role in driving the evolution of protoplanetary disks through angular momentum transport, observational evidence of magnetic fields has only been found in a small number of disks. Although dust continuum linear polarization has been detected in an increasing number of disks, its pattern is more consistent with that from dust scattering than from magne… ▽ More While magnetic fields likely play an important role in driving the evolution of protoplanetary disks through angular momentum transport, observational evidence of magnetic fields has only been found in a small number of disks. Although dust continuum linear polarization has been detected in an increasing number of disks, its pattern is more consistent with that from dust scattering than from magnetically aligned grains in the vast majority of cases. Continuum linear polarization from dust grains aligned to a magnetic field can reveal information about the magnetic field's direction, but not its strength. On the other hand, observations of circular polarization in molecular lines produced by Zeeman splitting offer a direct measure of the line-of-sight magnetic field strength in disks. We present upper limits on the net toroidal and vertical magnetic field strengths in the protoplanetary disk AS 209 derived from Zeeman splitting observations of the CN 2-1 line. The 3$σ$ upper limit on the net line-of-sight magnetic field strength in AS 209 is 5.0 mG on the redshifted side of the disk and 4.2 mG on the blueshifted side of the disk. Given the disk's inclination angle, we set a 3$σ$ upper limit on the net toroidal magnetic field strength of 8.7 and 7.3 mG for the red and blue sides of the disk, respectively, and 6.2 and 5.2 mG on the net vertical magnetic field on the red and blue sides of the disk. If magnetic disk winds are a significant mechanism of angular momentum transport in the disk, magnetic fields of a strength close to the upper limits would be sufficient to drive accretion at the rate previously inferred for regions near the protostar. △ Less

Submitted 22 June, 2021; v1 submitted 5 January, 2021; originally announced January 2021.

Comments: June 22, 2021: Corrected disk size scale bar in Figures 1, 2 and 5. Corrected percent polarization scale bar label in Figure 5

Confirming the Explosive Outflow in G5.89 with ALMA

Authors: Luis A. Zapata, Paul T. P. Ho, Manuel Fernández-López, Estrella Guzmán Ccolque, Luis F. Rodriguez, José Reyes-Valdés, John Bally, Aina Palau, Masao Saito, Patricio Sanhueza, P. R. Rivera-Ortiz, A. Rodriguez-González

Abstract: The explosive molecular outflow detected decades ago in the Orion BN/KL region of massive star formation was considered to be a bizarre event. This belief was strengthened by the non detection of similar cases over the years with the only exception of the marginal case of DR21. Here, we confim a similar explosive outflow associated with the UCH$_{\rm II}$ region G5.89$-$0.39 that indicates that th… ▽ More The explosive molecular outflow detected decades ago in the Orion BN/KL region of massive star formation was considered to be a bizarre event. This belief was strengthened by the non detection of similar cases over the years with the only exception of the marginal case of DR21. Here, we confim a similar explosive outflow associated with the UCH$_{\rm II}$ region G5.89$-$0.39 that indicates that this phenomenon is not unique to Orion or DR21. Sensitive and high angular resolution ($\sim$ 0.1$''$) ALMA CO(2$-$1) and SiO(5$-$4) observations show that the molecular outflow in the massive star forming region G5.89$-$0.39 is indeed an explosive outflow with an age of about 1000 yrs and a liberated kinetic energy of 10$^{46-49}$ erg. Our new CO(2$-$1) ALMA observations revealed over 30 molecular filaments, with Hubble-like expansion motions, pointing to the center of UCH$_{\rm II}$ region. In addition, the SiO(5$-$4) observations reveal warmer and strong shocks very close to the origin of the explosion, confirming the true nature of the flow. A simple estimation for the occurrence of these explosive events during the formation of the massive stars indicates an event rate of once every $\sim$100 yrs, which is close to the supernovae rate. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: Accepted by ApJ Letters. The animation is located here: https://drive.google.com/file/d/1a1c9JMwm0Mf8CxlfuEzyh6XmCOsgOyMA/view?usp=sharing

Four annular structures in a protostellar disk less than 500,000 years old

Authors: Dominique M. Segura-Cox, Anika Schmiedeke, Jaime E. Pineda, Ian W. Stephens, Manuel Fernández-López, Leslie W. Looney, Paola Caselli, Zhi-Yun Li, Lee G. Mundy, Woojin Kwon, Robert J. Harris

Abstract: Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when… ▽ More Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when the protostar is still embedded in a larger scale dense envelope of gas and dust. While younger disks likely play an important role in the onset of planet formation, only within the past decade have detailed properties of disks in the youngest star-forming phases begun to be observed. Here we present 1.3 mm dust emission observations with 5 au resolution that show four annular substructures in the disk of the young (<500,000 years) protostar IRS 63. IRS 63, a single Class I source located in the nearby Ophiuchus molecular cloud (at a distance of 144 pc), is one of the brightest Class I protostars at (sub)millimeter wavelengths that also has a relatively large disk (>50 au). Multiple annular substructures observed towards disks at young times can act as an early foothold for dust grain growth, a prerequisite of planet formation. Whether planets already exist or not in the disk of IRS 63, it is clear that the planet formation process begins in the young protostellar phases, earlier than predicted by current planet formation theories. △ Less

Submitted 7 October, 2020; originally announced October 2020.

Comments: Accepted for publication in Nature, appearing in the 8 Oct 2020 issue, DOI: 10.1038/s41586-020-2779-6

Low Level Carbon Monoxide Line Polarization in two Protoplanetary Disks: HD 142527 and IM Lup

Authors: Ian W. Stephens, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Richard Teague

Abstract: Magnetic fields are expected to play an important role in accretion processes for circumstellar disks. Measuring the magnetic field morphology is difficult, especially since polarimetric (sub)millimeter continuum observations may not trace fields in most disks. The Goldreich-Kylafis (GK) effect suggests that line polarization is perpendicular or parallel to the magnetic field direction. We attempt… ▽ More Magnetic fields are expected to play an important role in accretion processes for circumstellar disks. Measuring the magnetic field morphology is difficult, especially since polarimetric (sub)millimeter continuum observations may not trace fields in most disks. The Goldreich-Kylafis (GK) effect suggests that line polarization is perpendicular or parallel to the magnetic field direction. We attempt to observe CO(2-1), $^{13}$CO(2-1), and C$^{18}$O(2-1) line polarization toward HD 142527 and IM Lup, which are large, bright protoplanetary disks. We use spatial averaging and spectral integration to search for signals in both disks, and detect a potential CO(2-1) Stokes $Q$ signal toward both disks. The total CO(2-1) polarization fractions are 1.57 $\pm$ 0.18% and 1.01 $\pm$ 0.10% for HD 142527 and IM Lup, respectively. Our Monte Carlo simulations indicate that these signals are marginal. We also stack Stokes parameters based on the Keplerian rotation, but no signal was found. Across the disk traced by dust of HD 142527, the 3$σ$ upper limits for $P_{\text{frac}}$ at 0.5$^{\prime\prime}$ ($\sim$80 au) resolution are typically less than 3% for CO(2-1) and $^{13}$CO(2-1) and 4% for C$^{18}$O(2-1). For IM Lup, the 3$σ$ upper limits for these three lines are typically less than 3%, 4%, and 12%, respectively. Upper limits based on our stacking technique are up to a factor of $\sim$10 lower, though stacking areas can potentially average out small-scale polarization structure. We also compare our continuum polarization at 1.3 mm to observations at 870 $μ$m from previous studies. The polarization in the northern dust trap of HD 142527 shows a significant change in morphology and an increase in $P_{\text{frac}}$ as compared to 870 $μ$m. For IM Lup, the 1.3 mm polarization may be more azimuthal and has a higher $P_{\text{frac}}$ than at 870 $μ$m. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: Accepted to ApJ

Tidal Interaction between the UX Tauri Disk A/C System Revealed by ALMA

Authors: Luis A. Zapata, Luis F. Rodrí guez, Manuel Fernández-López, Aina Palau, Robert Estalella, Mayra Osorio, Guillem Anglada, Nuria Huelamo

Abstract: We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$4… ▽ More We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$σ$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: Accepted by The Astrophysical Journal

Modeling the accretion disk around the high-mass protostar GGD 27-MM1

Authors: N. Añez-Lopez, M. Osorio, G. Busquet, J. M. Girart, E. Macıas, C. Carrasco-Gonzalez, S. Curiel, R. Estalella, M. Fernandez-Lopez, R. Galvan-Madrid, J. Kwon, J. M. Torrelles

Abstract: Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present r… ▽ More Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (Rdisk~170 au) and massive (5Msun), about 20% of the stellar mass of 20 Msun. We estimate the total dynamical mass of the star-disk system from the molecular line emission finding a range between 21 and 30 Msun, which is consistent with our model. We fit the density and temperature structures found by our model with power law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27-MM1 disk with that found in low- and intermediate-mass stars and discuss possible implications on the water snow line. We have also carried about a study of the distance based on Gaia DR2 data and the population of young stellar objects (YSOs) in this region, and from the extinction maps. We conclude that the source distance is in within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc). △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: 20 pages, 15 figures

Radio observations of evaporating objects in the Cygnus OB2 region

Authors: N. L. Isequilla, M. Fernandez-Lopez, P. Benaglia, C. H. Ishwara-Chandra, S. del Palacio

Abstract: We present observations of the Cygnus OB2 region obtained with the Giant Metrewave Radio Telescope (GMRT) at the frequencies of 325 MHz and 610 MHz. In this contribution we focus on the study of proplyd-like objects (also known as free-floating Evaporating Gas Globules or frEGGs) that typically show an extended cometary morphology. We identify eight objects previously studied at other wavelengths… ▽ More We present observations of the Cygnus OB2 region obtained with the Giant Metrewave Radio Telescope (GMRT) at the frequencies of 325 MHz and 610 MHz. In this contribution we focus on the study of proplyd-like objects (also known as free-floating Evaporating Gas Globules or frEGGs) that typically show an extended cometary morphology. We identify eight objects previously studied at other wavelengths and derive their physical properties by obtaining their optical depth at radio-wavelengths. Using their geometry and the photoionization rate needed to produce their radio-continuum emission, we find that these sources are possibly ionized by a contribution of the stars Cyg OB2 #9 and Cyg OB2 #22. Spectral index maps of the eight frEGGs were constructed, showing a flat spectrum in radio frequencies in general. We interpret these as produced by optically-thin ionized gas, although it is possible that a combination of thermal emission, not necessarily optically thin, produced by a diffuse gas component and the instrument response (which detects more diffuse emission at low frequencies) can artificially generate negative spectral indices. In particular, for the case of the Tadpole we suggest that the observed emission is not of non-thermal origin despite the presence of regions with negative spectral indices in our maps. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: Acepted Astronomy & Astrophysics, 9 pages, 4 figures

G5.89: An Explosive Outflow Powered by a Proto-Stellar Merger?

Authors: Luis A. Zapata, Paul T. P. Ho, Estrella Guzman-Ccolque, Manuel Fernandez-Lopez, Luis F. Rodriguez, John Bally, Patricio Sanhueza, Aina Palau, Masao Saito

Abstract: The explosive outflows are a newly-discovered family of molecular outflows associated with high-mass star forming regions. Such energetic events are possibly powered by the release of gravitational energy related with the formation of a (proto)stellar merger or a close stellar binary. Here, we present sensitive and high angular resolution observations (0.85$''$) archival CO(J=3-2) observations car… ▽ More The explosive outflows are a newly-discovered family of molecular outflows associated with high-mass star forming regions. Such energetic events are possibly powered by the release of gravitational energy related with the formation of a (proto)stellar merger or a close stellar binary. Here, we present sensitive and high angular resolution observations (0.85$''$) archival CO(J=3-2) observations carried out with the Submillimeter Array (SMA) of the high-mass star forming region G5.89$-$0.39 that reveal the possible presence of an explosive outflow. We find six well-defined and narrow straight filament-like ejections pointing back approximately to the center of an expanding molecular and ionized shell located at the center of this region. These high velocity ($-$120 to $+$100 km s$^{-1}$) filaments follow a Hubble-like velocity law with the radial velocities increasing with the projected distance. The estimated kinematical age of the filaments is about of 1000 yrs, a value similar to the dynamical age found for the expanding ionized shell. G5.89 is the thus the third explosive outflow reported in the galaxy (together with Orion BN-KL and DR21) and argues in favor of the idea that this is a frequent phenomenon. In particular, explosive outflows, in conjunction with runaway stars, demonstrate that dynamical interactions in such groups are a very important ingredient in star formation. △ Less

Submitted 8 April, 2019; originally announced April 2019.

Comments: Accepted by MNRAS Letters

Polarization in Disks

Authors: Ian Stephens, Zhi-Yun Li, Haifeng Yang, Akimasa Kataoka, Leslie W. Looney, Charles L. H. Hull, Manuel Fernández-López, Sarah I. Sadavoy, Woojin Kwon, Satoshi Ohashi, Ryo Tazaki, Dan Li, Thiem Hoang, Gesa H. -M. Bertrang, Carlos Carrasco-González, William R. F. Dent, Satoko Takahashi, Francesca Bacciotti, Felipe O. Alves, Josep M. Girart, Qizhou Zhang, Ramprasad Rao, Adriana Pohl, Marco Padovani, Daniele Galli , et al. (2 additional authors not shown)

Abstract: Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometrie… ▽ More Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Unveiling a cluster of protostellar disks around the massive protostar GGD27 MM1

Authors: G. Busquet, J. M. Girart, R. Estalella, M. Fernández-López, R. Galván-Madrid, G. Anglada, C. Carrasco-González, N. Añez-López, S. Curiel, M. Osorio, L. F. Rodríguez, J. M. Torrelles

Abstract: We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicat… ▽ More We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicates that GGD27 displays moderate subclustering. This result combined with the dynamical timescale of the radio jet (10000 years) suggests the youthfulness of the cluster. The lack of disk mass segregation signatures may support this too. We found a clear paucity of disks with Rdisk >100 au. The median value of the radius is 34 au, smaller than the median of 92 au for Taurus but comparable to the value found in Ophiuchus and in the Orion Nebula Cluster. In GGD27 there is no evidence of a distance-dependent disk mass distribution (i. e., disk mass depletion due to external photoevaporation), most likely due to the cluster youth. There is a clear deficit of disks for distances <0.02 pc. Only for distances >0.04 pc stars can form larger and more massive disks, suggesting that dynamical interactions far from the cluster center are weaker, although the small disks found could be the result of disk truncation. This work demonstrates the potential to characterize disks from low-mass YSOs in distant and massive (still deeply embedded) clustered environments. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

Journal ref: A&A 623, L8 (2019)

An Asymmetric Keplerian Disk Surrounding the O-type Protostar IRAS16547$-$4247

Authors: Luis A. Zapata, Guido Garay, Aina Palau, Luis F. Rodriguez, Manuel Fernandez-Lopez, Robert Estalella, Andres Guzman

Abstract: During the last decades, a great interest has emerged to know if even the most massive stars in our galaxy (namely the spectral O-type stars) are formed in a similar manner as the low- and intermediate-mass stars, that is, through the presence of accreting disks and powerful outflows. Here, using sensitive observations of the Atacama Large Millimeter/Submillimeter Array (ALMA), we report a resolve… ▽ More During the last decades, a great interest has emerged to know if even the most massive stars in our galaxy (namely the spectral O-type stars) are formed in a similar manner as the low- and intermediate-mass stars, that is, through the presence of accreting disks and powerful outflows. Here, using sensitive observations of the Atacama Large Millimeter/Submillimeter Array (ALMA), we report a resolved Keplerian disk (with fifteen synthesized beams across its major axis) surrounding the deeply embedded O-type protostar IRAS16547$-$4247. The disk shows some asymmetries that could arise because of the disk is unstable and fragmenting or maybe because of different excitation conditions within the disk. The enclosed mass estimated from the disk Keplerian radial velocities is 25$\pm$3 M$_\odot$. The molecular disk is at the base of an ionized thermal radio jet and is approximately perpendicular to the jet axis orientation. We additionally find the existence of a binary system of compact dusty objects at the center of the accreting disk, which indicates the possible formation of an O-type star and a companion of lower mass. This is not surprising due to the high binary fraction reported in massive stars. Subtracting the contribution of the dusty disk plus the envelope and the companion, we estimated a mass of 20 M$_\odot$ for the central star. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: Accepted by The Astrophysical Journal

ALMA Observations of Polarized 872 $μ$m Dust Emission from the Protostellar Systems VLA 1623 and L1527

Authors: Robert J. Harris, Erin G. Cox, Leslie W. Looney, Zhi-Yun Li, Haifeng Yang, Manuel Fernández-López, Woojin Kwon, Sarah Sadavoy, Dominique Segura-Cox, Ian Stephens, John Tobin

Abstract: We present high-sensitivity ($σ_I\sim0.2-0.5$ mJy, $σ_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $μ$m dust emission from the young multiple system VLA 1623 in $ρ$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W… ▽ More We present high-sensitivity ($σ_I\sim0.2-0.5$ mJy, $σ_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $μ$m dust emission from the young multiple system VLA 1623 in $ρ$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W, and L1527 strongly in the polarized emission, at the $\sim1-3\%$ level. We spatially resolve VLA 1623A into two sources, VLA 1623Aa and VLA 1623Ab, separated by $\sim30$ au and located within a cavity of radius $\sim50$ au within the circumbinary Keplerian disk, as well as the edge-on disk of VLA 1623W. The polarization angle of the emission is uniform across each protostellar source and nearly coincides with each disk's minor axis. The offsets between the minor axis and the polarization angle are not uniformly distributed at the $P\lesssim2\times10^{-4}$ level. The circumbinary disk surrounding VLA 1623Aab is azimuthally symmetrically polarized. Each compact source's emission is partially optically thick ($τ\gtrsim1$) at 872 $μ$m, complicating interpretations of polarization involving aligned grains. We find evidence against alignment by radiative flux in each source, particularly in the edge-on VLA 1623W and L1527. We detect astrometric offsets between the polarized emission and the total intensity in VLA 1623Aa, VLA 1623Ab, and VLA 1623B, as predicted if self-scattering in the optically thick limit operates. We conclude that self-scattering is likely responsible for disk-scale polarization at 872 $μ$m in these systems. △ Less

Submitted 22 May, 2018; originally announced May 2018.

Comments: 22 pages, 7 figures, accepted to Astrophysical Journal