-
Rotation curves in protoplanetary disks with thermal stratification
Authors:
Paola Martire,
Cristiano Longarini,
Giuseppe Lodato,
Giovanni P. Rosotti,
Andrew Winter,
Stefano Facchini,
Caitlyn Hardiman,
Myriam Benisty,
Jochen Stadler,
Andrés F. Izquierdo,
Leonardo Testi
Abstract:
In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal strat…
▽ More
In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal stratification on the azimuthal velocity of protoplanetary disks. High resolution gas observations are showing velocity differences between CO isotopologues, which cannot be adequately explained with vertically isothermal models. The aim of this work is to determine whether a stratified model can explain this discrepancy. We analytically solve the hydrostatic equilibrium for a stratified disk and we derive the azimuthal velocity. We test the model with SPH numerical simulations and then we use it to fit for star mass, disk mass and scale radius of the sources in the MAPS sample. In particular, we use 12CO and 13CO datacubes.
△ Less
Submitted 11 March, 2024; v1 submitted 19 February, 2024;
originally announced February 2024.
-
Massive clumps in W43-main: Structure formation in an extensively shocked molecular cloud
Authors:
Yuxin Lin,
Friedrich Wyrowski,
Hauyu Baobab Liu,
Yan Gong,
Olli Sipilä,
Andrés F. Izquierdo,
Timea Csengeri,
Adam Ginsburg,
Guang-Xing Li,
Silvia Spezzano,
Jaime E. Pineda,
Silvia Leurini,
Paola Caselli,
Karl M. Menten
Abstract:
W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm…
▽ More
W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm with an angular resolution of $\sim$0.1 pc towards five massive clumps in W43 main. We use CH$_{3}$CCH and H$_{2}$CS lines to trace the extended gas temperature and CH$_{3}$OH lines to probe the volume density of the dense gas ($\gtrsim$10$^{5}$ cm$^{-3}$). The emission of SiO (2-1) is extensive across the region ($\sim$4 pc) and is mostly contained within a low-velocity regime, hinting at a large-scale origin of the shocks. The position-velocity maps of multiple tracers show systematic spatio-kinematic offsets supporting the cloud-cloud collision/merging scenario. We identify an additional extended velocity component in CCH emission, which coincides with one of the velocity components of the larger scale $^{13}$CO (2-1) emission, likely representing an outer, less dense gas layer in the cloud merging process. We find that the V-shaped, asymmetric SiO wings are tightly correlated with localised gas density enhancements, which is direct evidence of dense gas formation and accumulation in shocks. We resolve two categories of NH$_{2}$D cores: ones exhibiting only subsonic to transonic velocity dispersion, and the others with an additional supersonic velocity dispersion. The centroid velocities of the latter cores are correlated with the shock front seen by SiO. The kinematics of the $\sim$0.1 pc NH$_{2}$D cores are heavily imprinted by shock activities, and may represent a population of early-stage cores forming around the shock interface.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
High turbulence in the IM Lup protoplanetary disk: Direct observational constraints from CN and C$_2$H emission
Authors:
T. Paneque-Carreño,
A. F. Izquierdo,
R. Teague,
A. Miotello,
E. A. Bergin,
R. Loomis,
E. F. van Dishoeck
Abstract:
Constraining turbulence in disks is key towards understanding their evolution through the transport of angular momentum. Until now measurements of high turbulence have remained elusive and methods for estimating turbulence relay mostly on complex radiative transfer models of the data. Using the disk emission from IM Lup, a source proposed to be undergoing Magneto-Rotational Instabilities (MRI) and…
▽ More
Constraining turbulence in disks is key towards understanding their evolution through the transport of angular momentum. Until now measurements of high turbulence have remained elusive and methods for estimating turbulence relay mostly on complex radiative transfer models of the data. Using the disk emission from IM Lup, a source proposed to be undergoing Magneto-Rotational Instabilities (MRI) and possibly have high turbulence values in the upper disk layers, we present a new way of directly measuring turbulence without need of radiative transfer or thermochemical models. Through the characterization of the CN and C$_2$H emission in IM Lup, we aim to connect the information on the vertical and thermal structure of a particular disk region to derive turbulence at that location. By using an optically thin tracer it is possible to directly measure turbulence from the non-thermal broadening of the line. The vertical layers of the CN and C$_2$H emission are traced directly from the channel maps using ALFAHOR. By comparing their position to that of optically thick CO observations we are able to characterize the kinetic temperature of the emitting region. Using a simple parametric model of the line intensity with DISCMINER we accurately measure the emission linewidth and separate the thermal and non-thermal components. Assuming that the non-thermal component is fully turbulent, we are able to directly estimate the turbulent motions at the studied radial and vertical location of CN emission. IM Lup shows high turbulence of Mach 0.4-0.6 at $z/r \sim$ 0.25. Considering previous estimates of low turbulence near the midplane, this may indicate a vertical gradient in the disk turbulence, which is a key prediction of MRI studies. CN and C$_2$H are both emitting from a localized upper disk region at $z/r =$0.2-0.3, in agreement with thermochemical models.
△ Less
Submitted 6 February, 2024; v1 submitted 7 December, 2023;
originally announced December 2023.
-
The Disc Miner II: Revealing Gas substructures and Kinematic signatures from Planet-disc interaction through line profile analysis
Authors:
Andres F. Izquierdo,
Leonardo Testi,
Stefano Facchini,
Giovanni P. Rosotti,
Ewine van Dishoeck,
Lisa Wölfer,
Teresa Paneque-Carreño
Abstract:
[Abridged] The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the DISCMINER modelling fr…
▽ More
[Abridged] The aim of this work is to identify potential signatures from planet-disc interaction in the circumstellar discs around MWC 480, HD 163296, AS 209, IM Lup, and GM Aur, through the study of molecular lines observed as part of the ALMA large program MAPS. Extended and localised perturbations in velocity, line width, and intensity have been analysed jointly using the DISCMINER modelling framework, in three bright CO lines, 12CO, 13CO, and C18O $J=2-1$, to provide a comprehensive summary of the kinematic and column density substructures that planets might be actively sculpting in these discs. We find convincing evidence for the presence of four giant planets located at wide orbits in three of the discs in the sample: two around HD 163296, one in MWC 480, and one in AS 209. One of the planet candidates in HD 163296, P94, previously associated with velocity signatures detected in lower velocity resolution 12CO data, is confirmed and linked to localised velocity and line width perturbations in 13CO and C18O too. We highlight that line widths are also powerful tracers of planet-forming sites as they are sensitive to turbulent motions triggered by planet-disc interactions. In MWC 480, we identified non-axisymmetric line width enhancements around the radial separation of candidate planet-driven buoyancy spirals, which we used to narrow the location of the possible planet to an orbital radius of $R=245$ au and $\rm{PA}=193^\circ$. In the disc of AS 209, we found excess 12CO line widths centred at $R=210$ au, $\rm{PA}=151^\circ$, spanning around the immediate vicinity of a circumplanetary disc candidate proposed previously, which further supports its presence. Our simultaneous analysis of multiple tracers and observables aims to lay the groundwork for robust studies of molecular line properties focused on the search for young planets in discs.
△ Less
Submitted 15 May, 2023; v1 submitted 7 April, 2023;
originally announced April 2023.
-
A kinematically detected planet candidate in a transition disk
Authors:
Jochen Stadler,
Myriam Benisty,
Andrés F. Izquierdo,
Stefano Facchini,
Richard Teague,
Nicolas Kurtovic,
Paola Pinilla,
Jaehan Bae,
Megan Ansdell,
Ryan Loomis,
Satoshi Mayama,
Laura M. Pérez,
Leonardo Testi
Abstract:
Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps an…
▽ More
Transition disks are protoplanetary disks with inner cavities possibly cleared by massive companions, which makes them prime targets to observe at high resolution to map their velocity structure. We present ALMA Band 6 dust and gas observations of the transition disk around RXJ1604.3-2130 A, known to feature nearly symmetric shadows in scattered light. We studied the $^{12}$CO line channel maps and moment maps of the line of sight velocity and peak intensity. We fitted a Keplerian model of the channel-by-channel emission to study line profile differences and produced deprojected radial profiles for all velocity components. The $^{12}$CO emission shows a cavity inwards of $\sim$56 au and within the dust continuum ring at 81 au. The azimuthal brightness variations in the $^{12}$CO line and dust continuum are broadly aligned with the shadows detected in scattered-light observations. We find a strong localized non-Keplerian feature toward the west within the continuum ring (at $R=41\pm 10$ au and $PA=280\pm 2 ^\circ$). A tightly wound spiral is also detected which extends over 300$^\circ$ in azimuth, possibly connected to the localized non-Keplerian feature. Finally, a bending of the iso-velocity contours within the gas cavity indicates a highly perturbed inner region, possibly related to the presence of a misaligned inner disk. While broadly aligned with the scattered-light shadows, the localized non-Keplerian feature cannot be solely due to changes in temperature. Instead, we interpret the kinematical feature as tracing a massive companion located at the edge of the dust continuum ring. We speculate that the spiral is caused by buoyancy resonances driven by planet-disk interactions. However, this potential planet at $\sim$41 au cannot explain the gas-depleted cavity, the low accretion rate, and the misaligned inner disk, which suggests the presence of another companion closer in.
△ Less
Submitted 11 January, 2023; v1 submitted 4 January, 2023;
originally announced January 2023.
-
Directly tracing the vertical stratification of molecules in protoplanetary disks
Authors:
T. Paneque-Carreño,
A. Miotello,
E. F. van Dishoeck,
B. Tabone,
A. F. Izquierdo,
S. Facchini
Abstract:
We aim to directly trace the vertical location of the emitting surface of multiple molecular tracers in protoplanetary disks. Our sample of disks includes Elias 2-27, WaOph 6 and the sources targeted by the MAPS ALMA Large Program. The set of molecules studied include CO isotopologues in various transitions, HCN, CN, H2CO, HCO+, C2H and c-C3H2. The vertical emitting region is determined directly f…
▽ More
We aim to directly trace the vertical location of the emitting surface of multiple molecular tracers in protoplanetary disks. Our sample of disks includes Elias 2-27, WaOph 6 and the sources targeted by the MAPS ALMA Large Program. The set of molecules studied include CO isotopologues in various transitions, HCN, CN, H2CO, HCO+, C2H and c-C3H2. The vertical emitting region is determined directly from the channel maps, implementing accurate masking of the channel emission to recover the vertical location of the emission surface even at large radial distances from the star and for low-SNR lines. The vertical location of the emitting layer is obtained for 4-10 lines in each disk. IM Lup, HD163296 and MWC 480 12CO and 13CO show vertical modulations, which are coincident with dust gaps and kinematical perturbations. We also present estimates of the gas pressure scale height in the disks from the MAPS sample. Compared to physical-chemical models we find good agreement with the vertical location of CO isotopologues. In HD 163296 CN and HCN trace a similar intermediate layer, for the other disks, the UV flux tracers and the vertical profiles of HCN and C2H are lower than predicted in theoretical models. HCN and H2CO show a highly structured vertical profile, possibly indicative of different formation pathways. It is possible to trace the vertical locations of multiple molecular species that trace a wide variety of physical and chemical disk properties. The distribution of CO isotopologues are found at a wide range of vertical heights $z/r = $ 0.5-0.05. Other molecular lines are mostly found at $z/r \leq $0.15. The vertical layering of molecules is in agreement with theory in some systems, but not in all, therefore dedicated chemical-physical models are needed to further study and understand the emission surfaces.
△ Less
Submitted 3 October, 2022;
originally announced October 2022.
-
Kinematics and brightness temperatures of transition discs -- A survey of gas substructures as seen with ALMA
Authors:
Lisa Wölfer,
Stefano Facchini,
Nienke van der Marel,
Ewine F. van Dishoeck,
Myriam Benisty,
Alexander J. Bohn,
Logan Francis,
Andrés F. Izquierdo,
Richard D. Teague
Abstract:
In recent years, high-angular-resolution observations of the dust and gas in circumstellar discs have revealed a variety of morphologies, naturally triggering the question of whether these substructures are driven by forming planets. While it remains difficult to directly image embedded planets, a promising method to distinguish disc-shaping mechanisms is to study the gas kinematics as characteris…
▽ More
In recent years, high-angular-resolution observations of the dust and gas in circumstellar discs have revealed a variety of morphologies, naturally triggering the question of whether these substructures are driven by forming planets. While it remains difficult to directly image embedded planets, a promising method to distinguish disc-shaping mechanisms is to study the gas kinematics as characterising deviations from Keplerian rotation can be used to probe underlying perturbations such as planets. Creating spiral structures, the latter can also be traced in the brightness temperature. Here we analyse the brightness temperatures and kinematics of a sample of 36 transition discs observed with ALMA to search for substructures possibly tracing companions. We use archival Band 6 and 7 observations of different CO isotopologues and fit Keplerian disc models to the velocity fields. After subtraction of an azimuthally averaged brightness temperature and Keplerian rotation model from the data, we find significant substructures in both residuals of eight sources. Other sources show tentative features, while half of our sample does not show any substructures that may be indicative of planet-disc interactions. For the first time, we compare the substructures from our analysis with various other indicators of planets. About 20% of discs show strong features such as arcs or spirals, possibly associated with the presence of planets, while the majority do not present as clear planet-driven signatures. Almost all discs that exhibit spirals in near-infrared scattered light show at least tentative features in the CO data. The present data are able to reveal only very massive bodies and a lack of features may suggest that, if there are planets, they are of lower mass (<1-3Mj) or located closer to the star within deep cavities. Dedicated observations and modelling efforts are needed to confirm such scenarios.
△ Less
Submitted 19 August, 2022;
originally announced August 2022.
-
A giant planet shaping the disk around the very low-mass star CIDA 1
Authors:
P. Curone,
A. F. Izquierdo,
L. Testi,
G. Lodato,
S. Facchini,
A. Natta,
P. Pinilla,
N. T. Kurtovic,
C. Toci,
M. Benisty,
M. Tazzari,
F. Borsa,
M. Lombardi,
C. F. Manara,
E. Sanchis,
L. Ricci
Abstract:
(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, s…
▽ More
(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, show substructures hinting at the presence of a massive planet. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its $^{12}$CO (J=3-2) and $^{13}$CO (J=3-2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet. We model the protoplanetary disk with a set of hydrodynamical simulations, varying the mass and locations of the embedded planet. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations treating the images as the actual ALMA observations. Our models indicate that a planet with a minimum mass of $\sim1.4\,\text{M}_\text{Jup}$ orbiting at a distance of $\sim 9-10$ au can explain the morphology and location of the observed dust ring at Band 7 and Band 4. We can reproduce the low spectral index ($\sim 2$) observed where the dust ring is detected. Our synthetic images reproduce the morphology of the $^{12}$CO and $^{13}$CO observed channel maps where the cloud absorption allowed a detection. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of $\sim4 - 8\,\text{M}_\text{Jup}$. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars.
△ Less
Submitted 17 August, 2022; v1 submitted 20 May, 2022;
originally announced May 2022.
-
A new planet candidate detected in a dust gap of the disc around HD 163296 through localised kinematic signatures. An observational validation of the Discminer
Authors:
Andres F. Izquierdo,
Stefano Facchini,
Giovanni P. Rosotti,
Ewine F. van Dishoeck,
Leonardo Testi
Abstract:
We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the c…
▽ More
We report the robust detection of coherent, localised deviations from Keplerian rotation possibly associated with the presence of two giant planets embedded in the disc around HD 163296. The analysis is performed using the DISCMINER channel-map modelling framework on $^{12}$CO $J=2-1$ DSHARP data. Not only orbital radius, but also azimuth of the planets are retrieved by our technique. One of the candidate planets, detected at $R=94\pm6$ au, $φ=50\pm3^\circ$ (P94), is near the centre of one of the gaps in dust continuum emission, and is consistent with a planet mass of 1 $\rm M_{\rm Jup}$. The other planet, located at $R=261\pm4$ au, $φ=57\pm1^\circ$ (P261), is in the region where a velocity kink was previously observed in $^{12}$CO channel maps. Also, we provide a simultaneous description of the height and temperature of the upper and lower emitting surfaces of the disc, and propose the line width as a solid observable to track gas substructure. Using azimuthally averaged line width profiles we detect gas gaps at $R=38$ au, $R=88$ au, and $R=136$ au, closely matching the location of their dust and kinematical counterparts. Furthermore, we observe strong azimuthal asymmetries in line widths around the gas gap at $R=88$ au, possibly linked to turbulent motions driven by the P94 planet. Our results confirm that the DISCMINER is capable of finding localised, otherwise unseen velocity perturbations thanks to its robust statistical framework, but also that it is well suited for studies of the gas properties and vertical structure of protoplanetary discs.
△ Less
Submitted 7 January, 2022; v1 submitted 11 November, 2021;
originally announced November 2021.
-
Zooming into the Collimation Zone in a Massive Protostellar Jet
Authors:
Carlos Carrasco-González,
Alberto Sanna,
Adriana Rodríguez-Kamenetzky,
Luca Moscadelli,
Melvin Hoare,
José M. Torrelles,
Roberto Galván-Madrid,
Andrés F. Izquierdo
Abstract:
Protostellar jets have a fundamental role at the earliest evolution of protostars of all masses. In the case of low-mass (<8 Msun) protostars, strong observational evidence exists that the launching and collimation is due to the X- and/or disk-wind mechanisms. In these models, it is the protostar/disk system that creates all the necessary conditions to launch and collimate the jets near the protos…
▽ More
Protostellar jets have a fundamental role at the earliest evolution of protostars of all masses. In the case of low-mass (<8 Msun) protostars, strong observational evidence exists that the launching and collimation is due to the X- and/or disk-wind mechanisms. In these models, it is the protostar/disk system that creates all the necessary conditions to launch and collimate the jets near the protostar via strong magnetic fields. The origin of jets from more massive protostars has been investigated much less, in part because of the difficulty of resolving the collimation zone in these more distant objects. Here we present the highest angular resolution observations of a jet powered by a massive protostar, the Cep A HW2 radio jet. We imaged the radio emission at projected distances of only ~20 au from the protostar, resolving the innermost 100 au of a massive protostellar jet for the first time. The morphology of the radio jet emission in this massive object is very different than what is usually observed in jets from low-mass protostars. We found that the outflowing material in HW2 has two components: a wide-angle wind launched from the protostar/disk system, and a highly collimated jet starting at 20-30 au from the protostar. We discuss two possible scenarios: an extension of the classical disk-wind to a massive protostar, or external collimation of a wide-angle wind. These results have important consequences for our understanding of how stars of different masses are formed.
△ Less
Submitted 2 June, 2021;
originally announced June 2021.
-
The Disc Miner I: A statistical framework to detect and quantify kinematical perturbations driven by young planets in discs
Authors:
Andres F. Izquierdo,
Leonardo Testi,
Stefano Facchini,
Giovanni P. Rosotti,
Ewine F. van Dishoeck
Abstract:
[Abridged] The study of disc kinematics has recently opened up as a promising method to detect unseen planets. However, a systematic, statistically meaningful analysis of such an approach remains missing. The aim of this work is to devise an automated, statistically robust technique to identify kinematical perturbations induced by the presence of planets in a gas disc, and to accurately infer thei…
▽ More
[Abridged] The study of disc kinematics has recently opened up as a promising method to detect unseen planets. However, a systematic, statistically meaningful analysis of such an approach remains missing. The aim of this work is to devise an automated, statistically robust technique to identify kinematical perturbations induced by the presence of planets in a gas disc, and to accurately infer their location. For this purpose, we produce hydro simulations of planet-disc interactions with different planet masses, 0.3, 1.0 and 3.0 $M_{Jup}$, at a radius of $R=100$ au in the disc, and perform radiative transfer calculations of CO to simulate observables for 13 planet azimuths. Using the DISCMINER package, we fit the synthetic data cubes with a Keplerian model of the channel-by-channel emission to study line profile differences, including deviations from Keplerian rotation. The detection technique, based on line centroid differences, captures localised planet-driven perturbations, and can distinguish them from axisymmetric velocity perturbations. The method can detect all three simulated planets, at all azimuths, with an average accuracy of $\pm3^\circ$ in azimuth and $\pm8$ au in radius. Owing to disc structure and line-of-sight projection effects, planets at azimuths close to $\pm45^\circ$ yield the highest velocity fluctuations, whereas those at limiting cases, $0^\circ$ and $\pm90^\circ$, drive the lowest. The observed peak velocities range within 40$-$70, 70$-$170 and 130$-$450 m s$^{-1}$ for 0.3, 1.0 and 3.0 $M_{Jup}$ planets. Our analysis indicates that the variance of peak velocities is boosted near planets due to organised gas motions prompted by their localised gravitational well. We propose an approach that exploits this velocity coherence to provide, for the first time, statistically significant detections of localised planet-driven perturbations in the gas disc kinematics.
△ Less
Submitted 26 June, 2021; v1 submitted 19 April, 2021;
originally announced April 2021.
-
The Cloud Factory II: Gravoturbulent Kinematics of Resolved Molecular Clouds in a Galactic Potential
Authors:
Andres F. Izquierdo,
Rowan J. Smith,
Simon C. O. Glover,
Ralf S. Klessen,
Robin G. Tress,
Mattia C. Sormani,
Paul C. Clark,
Ana Duarte-Cabral,
Catherine Zucker
Abstract:
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are…
▽ More
We present a statistical analysis of the gravoturbulent velocity fluctuations in molecular cloud complexes extracted from our "Cloud Factory" galactic-scale ISM simulation suite. For this purpose, we produce non-LTE $^{12}$CO J=1-0 synthetic observations and apply the Principal Component Analysis (PCA) reduction technique on a representative sample of cloud complexes. The velocity fluctuations are self-consistently generated by different physical mechanisms at play in our simulations, which include galactic-scale forces, gas self-gravity, and supernova feedback. The statistical analysis suggests that, even though purely gravitational effects are necessary to reproduce standard observational laws, they are not sufficient in most cases. We show that the extra injection of energy from supernova explosions plays a key role in establishing the global turbulent field and the local dynamics and morphology of molecular clouds. Additionally, we characterise structure function scaling parameters as a result of cloud environmental conditions: some of the complexes are immersed in diffuse (inter-arm) or dense (spiral-arm) environments, and others are influenced by embedded or external supernovae. In quiescent regions, we obtain time-evolving trajectories of scaling parameters driven by gravitational collapse and supersonic turbulent flows. Our findings suggests that a PCA-based statistical study is a robust method to diagnose the physical mechanisms that drive the gravoturbulent properties of molecular clouds. Also, we present a new open source module, the PCAFACTORY, which smartly performs PCA to extract velocity structure functions from simulated or real data of the ISM in a user-friendly way. Software DOI: 10.5281/zenodo.3822718
△ Less
Submitted 4 November, 2020;
originally announced November 2020.
-
The history of dynamics and stellar feedback revealed by the HI filamentary structure in the disk of the Milky Way
Authors:
J. D. Soler,
H. Beuther,
J. Syed,
Y. Wang,
L. D. Anderson,
S. C. O. Glover,
P. Hennebelle,
M. Heyer,
Th. Henning,
A. F. Izquierdo,
R. S. Klessen,
H. Linz,
N. M. McClure-Griffiths,
J. Ott,
S. E. Ragan,
M. Rugel,
N. Schneider,
R. J. Smith,
M. C. Sormani,
J. M. Stil,
R. Treß,
J. S. Urquhart
Abstract:
We present a study of the filamentary structure in the emission from the neutral atomic hydrogen (HI) at 21 cm across velocity channels in the 40'' and 1.5-km/s resolution position-position-velocity cube resulting from the combination of the single-dish and interferometric observations in The HI/OH/Recombination (THOR) line survey. Using the Hessian matrix method in combination with tools from cir…
▽ More
We present a study of the filamentary structure in the emission from the neutral atomic hydrogen (HI) at 21 cm across velocity channels in the 40'' and 1.5-km/s resolution position-position-velocity cube resulting from the combination of the single-dish and interferometric observations in The HI/OH/Recombination (THOR) line survey. Using the Hessian matrix method in combination with tools from circular statistics, we find that the majority of the filamentary structures in the HI emission are aligned with the Galactic plane. Part of this trend can be assigned to long filamentary structures that are coherent across several velocity channels. However, we also find ranges of Galactic longitude and radial velocity where the HI filamentary structures are preferentially oriented perpendicular to the Galactic plane. These are located (i) around the tangent point of the Scutum spiral arm, $l \approx 28^{\circ}$ and $v_{\rm LSR}\approx 100$ km/s, (ii) toward $l \approx 45^{\circ}$ and $v_{\rm LSR}\approx 50$ km/s, (iii) around the Riegel-Crutcher cloud, and (iv) toward the terminal velocities. Comparison with numerical simulations indicates that the prevalence of horizontal filamentary structures is most likely the result of the large-scale dynamics and that vertical structures identified in (i) and (ii) may arise from the combined effect of supernova (SN) feedback and strong magnetic fields. The vertical filamentary structures in (iv) can be related to the presence of clouds from extra-planar HI gas falling back into the Galactic plane after being expelled by SNe. Our results indicate that a systematic characterization of the emission morphology toward the Galactic plane provides an unexplored link between the observations and the dynamical behaviour of the interstellar medium, from the effect of large-scale Galactic dynamics to the Galactic fountains driven by SNe.
△ Less
Submitted 10 September, 2020; v1 submitted 14 July, 2020;
originally announced July 2020.
-
The Cloud Factory I: Generating resolved filamentary molecular clouds from galactic-scale forces
Authors:
Rowan J. Smith,
Robin G. Treß,
Mattia C. Sormani,
Simon C. O. Glover,
Ralf S. Klessen,
Paul C. Clark,
Andrés F. Izquierdo,
Ana Duarte Cabral,
Catherine Zucker
Abstract:
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitati…
▽ More
We introduce a new suite of simulations, "The Cloud Factory", which self-consistently forms molecular cloud complexes at high enough resolution to resolve internal substructure (up to 0.25 Msol in mass) all while including galactic-scale forces. We use a version of the Arepo code modified to include a detailed treatment of the physics of the cold molecular ISM, and an analytical galactic gravitational potential for computational efficiency. The simulations have nested levels of resolution, with the lowest layer tied to tracer particles injected into individual cloud complexes. These tracer refinement regions are embedded in the larger simulation so continue to experience forces from outside the cloud. This allows the simulations to act as a laboratory for testing the effect of galactic environment on star formation. Here we introduce our method and investigate the effect of galactic environment on filamentary clouds. We find that cloud complexes formed after a clustered burst of feedback, have shorter lengths and are less likely to fragment compared to quiescent clouds (e.g. the Musca filament) or those dominated by the galactic potential (e.g. Nessie). Spiral arms and differential rotation preferentially align filaments, but strong feedback randomises them. Long filaments formed within the cloud complexes are necessarily coherent with low internal velocity gradients, which has implications for the formation of filamentary star-clusters. Cloud complexes formed in regions dominated by supernova feedback have fewer star-forming cores, and these are more widely distributed. These differences show galactic-scale forces can have a significant impact on star formation within molecular clouds.
△ Less
Submitted 18 November, 2019; v1 submitted 13 November, 2019;
originally announced November 2019.
-
On the Effects of Self-Obscuration in the (Sub-)Millimeter Spectral Indices and Appearance of Protostellar Disks
Authors:
Roberto Galván-Madrid,
Hauyu Baobab Liu,
Andrés F. Izquierdo,
Anna Miotello,
Bo Zhao,
Carlos Carrasco-González,
Susana Lizano,
Luis F. Rodríguez
Abstract:
In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $α_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclina…
▽ More
In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $α_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $ν$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $β$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $β=1.7$ gives better results than $β=1$. A wide distribution of disk masses is also needed to produce some cases with $α<2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account.
△ Less
Submitted 1 October, 2018;
originally announced October 2018.
-
Radiative transfer modelling of W33A MM1: 3-D structure and dynamics of a complex massive star forming region
Authors:
Andrés F. Izquierdo,
Roberto Galván-Madrid,
Luke T. Maud,
Melvin G. Hoare,
Katharine G. Johnston,
Eric R. Keto,
Qizhou Zhang,
Willem-Jan de Wit
Abstract:
We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represen…
▽ More
We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.