-
Filamentary Hierarchies and Superbubbles: Galactic Multiscale MHD Simulations of GMC to Star Cluster Formation
Authors:
Bo Zhao,
Ralph E. Pudritz,
Rachel Pillsworth,
Hector Robinson,
James Wadsley
Abstract:
There is now abundant observational evidence that star formation is a highly dynamical process that connects filament hierarchies and supernova feedback from galaxy scale kpc filaments and superbubbles, to giant molecular clouds (GMCs) on 100 pc scales and star clusters (1 pc). Here we present galactic multi-scale MHD simulations that track the formation of structure from galactic down to sub pc s…
▽ More
There is now abundant observational evidence that star formation is a highly dynamical process that connects filament hierarchies and supernova feedback from galaxy scale kpc filaments and superbubbles, to giant molecular clouds (GMCs) on 100 pc scales and star clusters (1 pc). Here we present galactic multi-scale MHD simulations that track the formation of structure from galactic down to sub pc scales in a magnetized, Milky Way like galaxy undergoing supernova driven feedback processes. We do this by adopting a novel zoom-in technique that follows the evolution of typical 3-kpc sub regions without cutting out the surrounding galactic environment, allowing us to reach 0.28 pc resolution in the individual zoom-in regions. We find a wide range of morphologies and hierarchical structure including superbubbles, turbulence, kpc atomic gas filaments hosting multiple GMC condensations that are often associated with superbubble compression; down to smaller scale filamentary GMCs and star cluster regions within them. Gas accretion and compression ultimately drive filaments over a critical, scale - dependent, line mass leading to gravitational instabilities that produce GMCs and clusters. In quieter regions, galactic shear can produce filamentary GMCs within flattened, rotating disk-like structures on 100 pc scales. Strikingly, our simulations demonstrate the formation of helical magnetic fields associated with the formation of these disk like structures.
△ Less
Submitted 28 May, 2024;
originally announced May 2024.
-
Observations of high-order multiplicity in a high-mass stellar protocluster
Authors:
Shanghuo Li,
Patricio Sanhueza,
Henrik Beuther,
Huei-Ru Vivien Chen,
Rolf Kuiper,
Fernando A. Olguin,
Ralph E. Pudritz,
Ian W. Stephens,
Qizhou Zhang,
Fumitaka Nakamura,
Xing Lu,
Rajika L. Kuruwita,
Takeshi Sakai,
Thomas Henning,
Kotomi Taniguchi,
Fei Li
Abstract:
The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mas…
▽ More
The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23--0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation.
△ Less
Submitted 12 January, 2024;
originally announced January 2024.
-
Necessary conditions for the formation of filaments and star clusters in the cold neutral medium
Authors:
Rachel Pillsworth,
Ralph E. Pudritz
Abstract:
Star formation takes place in filamentary molecular clouds which arise by physical processes that take place in the cold, neutral medium (CNM). We address the necessary conditions for this diffuse ($n \approx 30$ cm$^{-3}$), cold (T $\approx$ 60 K), magnetized gas undergoing shock waves and supersonic turbulence, to produce filamentary structures capable of fragmenting into cluster forming regions…
▽ More
Star formation takes place in filamentary molecular clouds which arise by physical processes that take place in the cold, neutral medium (CNM). We address the necessary conditions for this diffuse ($n \approx 30$ cm$^{-3}$), cold (T $\approx$ 60 K), magnetized gas undergoing shock waves and supersonic turbulence, to produce filamentary structures capable of fragmenting into cluster forming regions. Using RAMSES and a magnetized CNM environment as our initial conditions, we simulate a 0.5 kpc turbulent box to model a uniform gas with magnetic field strength of 7 $μG$, varying the 3D velocity dispersion via decaying turbulence. We use a surface density of $320 M_{\odot} pc^{-2}$, representative of the inner 4.0 kpc CMZ of the Milky Way and typical luminous galaxies. Filamentary molecular clouds are formed dynamically via shocks within a narrow range of velocity dispersions in the CNM of 5 - 10 km/s with a preferred value at 8 km/s. Cluster sink particles appear in filaments which exceed their critical line mass, occurring optimally for velocity dispersions of 8 km/s. Tracking the evolution of magnetic fields, we find that they lead to double the dense star forming gas than in purely hydro runs. Perpendicular orientations between magnetic field and filaments can increase the accretion rates onto filaments and hence their line masses. Because magnetic fields help support gas, MHD runs result in average temperatures an order of magnitude higher than unmagnetized counterparts. Finally, we find magnetic fields delay the onset of cluster formation by $\propto 0.4$ Myr.
△ Less
Submitted 3 January, 2024;
originally announced January 2024.
-
Density distributions, magnetic field structures and fragmentation in high-mass star formation
Authors:
H. Beuther,
C. Gieser,
J. D. Soler,
Q. Zhang,
R. Rao,
D. Semenov,
Th. Henning,
R. Pudritz,
T. Peters,
P. Klaassen,
M. T. Beltran,
A. Palau,
T. Moeller,
K. G. Johnston,
H. Zinnecker,
J. Urquhart,
R. Kuiper,
A. Ahmadi,
A. Sanchez-Monge,
S. Feng,
S. Leurini,
S. E. Ragan
Abstract:
Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures.
Results…
▽ More
Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures.
Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE
Authors:
Aida Ahmadi,
H. Beuther,
F. Bosco,
C. Gieser,
S. Suri,
J. C. Mottram,
R. Kuiper,
Th. Henning,
Á. Sánchez-Monge,
H. Linz,
R. E. Pudritz,
D. Semenov,
J. M. Winters,
T. Möller,
M. T. Beltrán,
T. Csengeri,
R. Galván-Madrid,
K. G. Johnston,
E. Keto,
P. D. Klaassen,
S. Leurini,
S. N. Longmore,
S. L. Lumsden,
L. T. Maud,
L. Moscadelli
, et al. (6 additional authors not shown)
Abstract:
The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in…
▽ More
The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars.
△ Less
Submitted 3 May, 2023; v1 submitted 28 April, 2023;
originally announced May 2023.
-
Combined Effects of Disk Winds and Turbulence-Driven Accretion on Planet Populations
Authors:
Matthew Alessi,
Ralph E. Pudritz
Abstract:
Recent surveys show that protoplanetary disks have lower levels of turbulence than expected based on their observed accretion rates. A viable solution to this is that magnetized disk winds dominate angular momentum transport. This has several important implications for planet formation processes. We compute the physical and chemical evolution of disks and the formation and migration of planets und…
▽ More
Recent surveys show that protoplanetary disks have lower levels of turbulence than expected based on their observed accretion rates. A viable solution to this is that magnetized disk winds dominate angular momentum transport. This has several important implications for planet formation processes. We compute the physical and chemical evolution of disks and the formation and migration of planets under the combined effects of angular momentum transport by turbulent viscosity and disk winds. We take into account the critical role of planet traps to limit Type I migration in all of these models and compute thousands of planet evolution tracks for single planets drawn from a distribution of initial disk properties and turbulence strengths. We do not consider multi-planet models nor include N-body planet-planet interactions. Within this physical framework we find that populations with a constant value disk turbulence and winds strength produce mass-semimajor axis distributions in the M-a diagram with insufficient scatter to compare reasonably with observations. However, populations produced as a consequence of sampling disks with a distribution of the relative strengths of disk turbulence and winds fit much better. Such models give rise to a substantial super Earth population at orbital radii 0.03-2 AU, as well as a clear separation between the produced hot Jupiter and warm Jupiter populations. Additionally, this model results in a good comparison with the exoplanetary mass-radius distribution in the M-R diagram after post-disk atmospheric photoevaporation is accounted for.
△ Less
Submitted 4 July, 2022;
originally announced July 2022.
-
Possible Ribose Synthesis in Carbonaceous Planetesimals
Authors:
Klaus Paschek,
Kai Kohler,
Ben K. D. Pearce,
Kevin Lange,
Thomas K. Henning,
Oliver Trapp,
Ralph E. Pudritz,
Dmitry A. Semenov
Abstract:
The origin of life might be sparked by the polymerization of the first RNA molecules in Darwinian ponds during wet-dry cycles. The key life-building block ribose was found in carbonaceous chondrites. Its exogenous delivery onto the Hadean Earth could be a crucial step toward the emergence of the RNA world. Here, we investigate the formation of ribose through a simplified version of the formose rea…
▽ More
The origin of life might be sparked by the polymerization of the first RNA molecules in Darwinian ponds during wet-dry cycles. The key life-building block ribose was found in carbonaceous chondrites. Its exogenous delivery onto the Hadean Earth could be a crucial step toward the emergence of the RNA world. Here, we investigate the formation of ribose through a simplified version of the formose reaction inside carbonaceous chondrite parent bodies. Following up on our previous studies regarding nucleobases with the same coupled physico-chemical model, we calculate the abundance of ribose within planetesimals of different sizes and heating histories. We perform laboratory experiments using catalysts present in carbonaceous chondrites to infer the yield of ribose among all pentoses (5Cs) forming during the formose reaction. These laboratory yields are used to tune our theoretical model that can only predict the total abundance of 5Cs. We found that the calculated abundances of ribose were similar to the ones measured in carbonaceous chondrites. We discuss the possibilities of chemical decomposition and preservation of ribose and derived constraints on time and location in planetesimals. In conclusion, the aqueous formose reaction might produce most of the ribose in carbonaceous chondrites. Together with our previous studies on nucleobases, we found that life-building blocks of the RNA world could be synthesized inside parent bodies and later delivered onto the early Earth.
△ Less
Submitted 4 April, 2022;
originally announced April 2022.
-
Towards RNA life on Early Earth: From atmospheric HCN to biomolecule production in warm little ponds
Authors:
Ben K. D. Pearce,
Karan Molaverdikhani,
Ralph E. Pudritz,
Thomas Henning,
Kaitlin E. Cerrillo
Abstract:
The origin of life on Earth involves the early appearance of an information-containing molecule such as RNA. The basic building blocks of RNA could have been delivered by carbon-rich meteorites, or produced in situ by processes beginning with the synthesis of hydrogen cyanide (HCN) in the early Earth's atmosphere. Here, we construct a robust physical and non-equilibrium chemical model of the early…
▽ More
The origin of life on Earth involves the early appearance of an information-containing molecule such as RNA. The basic building blocks of RNA could have been delivered by carbon-rich meteorites, or produced in situ by processes beginning with the synthesis of hydrogen cyanide (HCN) in the early Earth's atmosphere. Here, we construct a robust physical and non-equilibrium chemical model of the early Earth atmosphere. The atmosphere is supplied with hydrogen from impact degassing of meteorites, sourced with water evaporated from the oceans, carbon dioxide from volcanoes, and methane from undersea hydrothermal vents, and in which lightning and external UV-driven chemistry produce HCN. This allows us to calculate the rain-out of HCN into warm little ponds (WLPs). We then use a comprehensive sources and sinks numerical model to compute the resulting abundances of nucleobases, ribose, and nucleotide precursors such as 2-aminooxazole resulting from aqueous and UV-driven chemistry within them. We find that at 4.4 bya (billion years ago) the limits of adenine concentrations in ponds for habitable surfaces is 0.05$μ$M in the absence of seepage. These concentrations can be maintained for over 100 Myr. Meteorite delivery of adenine to WLPs can provide boosts in concentration by 2-3 orders of magnitude, but these boosts deplete within months by UV photodissociation, seepage, and hydrolysis. The early evolution of the atmosphere is dominated by the decrease of hydrogen due to falling impact rates and atmospheric escape, and the rise of oxygenated species such as OH from H2O photolysis. Our work points to an early origin of RNA on Earth within ~200 Myr of the Moon-forming impact.
△ Less
Submitted 28 January, 2022; v1 submitted 3 January, 2022;
originally announced January 2022.
-
Meteorites and the RNA World: Synthesis of Nucleobases in Carbonaceous Planetesimals and the Role of Initial Volatile Content
Authors:
Klaus Paschek,
Dmitry A. Semenov,
Ben K. D. Pearce,
Kevin Lange,
Thomas K. Henning,
Ralph E. Pudritz
Abstract:
Prebiotic molecules, fundamental building blocks for the origin of life, have been found in carbonaceous chondrites. The exogenous delivery of these organic molecules onto the Hadean Earth could have sparked the polymerization of the first RNA molecules in Darwinian ponds during wet-dry cycles. Here, we investigate the formation of the RNA and DNA nucleobases adenine, uracil, cytosine, guanine, an…
▽ More
Prebiotic molecules, fundamental building blocks for the origin of life, have been found in carbonaceous chondrites. The exogenous delivery of these organic molecules onto the Hadean Earth could have sparked the polymerization of the first RNA molecules in Darwinian ponds during wet-dry cycles. Here, we investigate the formation of the RNA and DNA nucleobases adenine, uracil, cytosine, guanine, and thymine inside parent body planetesimals of carbonaceous chondrites. An up-to-date thermochemical equilibrium model coupled with a 1D thermodynamic planetesimal model is used to calculate the nucleobase concentrations. Different from the previous study (Pearce & Pudritz 2016), we assume initial volatile concentrations more appropriate for the formation zone of carbonaceous chondrite parent bodies. This represents more accurately cosmochemical findings that these bodies have formed inside the inner, $\sim 2\mathrm{-}5\,\mathrm{au}$, warm region of the solar system. Due to these improvements, our model represents the concentrations of adenine and guanine measured in carbonaceous chondrites. Our model did not reproduce per se the measurements of uracil, cytosine, and thymine in these meteorites. This can be explained by transformation reactions between nucleobases and potential decomposition of thymine. The synthesis of prebiotic organic matter in carbonaceous asteroids could be well explained by a combination of i) radiogenic heating, ii) aqueous chemistry involving a few key processes at a specific range of radii inside planetesimals where water can exist in the liquid phase, and iii) a reduced initial volatile content (H$_2$, CO, HCN, CH$_2$O) of the protoplanetary disk material in the parent body region compared to the outer region of comets.
△ Less
Submitted 11 November, 2022; v1 submitted 16 December, 2021;
originally announced December 2021.
-
Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories
Authors:
Jessica Speedie,
Ralph E. Pudritz,
Alex J. Cridland,
Farzana Meru,
Richard A. Booth
Abstract:
ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed thro…
▽ More
ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed throughout the disk, assuming either a high or low turbulent $α$ viscosity. We have discovered that there is a bifurcation in the migration history of forming planets as a consequence of varying the disk viscosity. In our high viscosity disks, inward migration prevails and yields compact planetary systems, tempered only by planet trapping at the water iceline around 5 au. In our lower viscosity models however, low mass planets can migrate outward to twice their initial orbital radii, driven by a radially extended region of strong outward-directed corotation torques located near the heat transition (where radiative heating of the disk by the star is comparable to viscous heating) - before eventually migrating inwards. We derive analytic expressions for the planet mass at which the corotation torque dominates, and find that this "corotation mass" scales as $M_{\rm p, corot} \sim α^{2/3}$. If disk winds dominate the corotation torque, the corotation mass scales linearly with wind strength. We propose that the observed bifurcation in disk demographics into a majority of compact dust disks and a minority of extended ring/gap systems is a consequence of a distribution of viscosity across the disk population.
△ Less
Submitted 4 February, 2022; v1 submitted 2 November, 2021;
originally announced November 2021.
-
Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953
Authors:
C. Gieser,
H. Beuther,
D. Semenov,
S. Suri,
J. D. Soler,
H. Linz,
J. Syed,
Th. Henning,
S. Feng,
T. Möller,
A. Palau,
J. M. Winters,
M. T. Beltrán,
R. Kuiper,
L. Moscadelli,
P. Klaassen,
J. S. Urquhart,
T. Peters,
S. N. Longmore,
Á. Sánchez-Monge,
R. Galván-Madrid,
R. E. Pudritz,
K. G. Johnston
Abstract:
We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from t…
▽ More
We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups.
△ Less
Submitted 14 October, 2021; v1 submitted 5 October, 2021;
originally announced October 2021.
-
Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3
Authors:
S. Suri,
H. Beuther,
C. Gieser,
A. Ahmadi,
Á. Sánchez-Monge,
J. M. Winters,
H. Linz,
Th. Henning,
M. T. Beltrán,
F. Bosco,
R. Cesaroni,
T. Csengeri,
S. Feng,
M. G. Hoare,
K. G. Johnston,
P. Klaasen,
R. Kuiper,
S. Leurini,
S. Longmore,
S. Lumsden,
L. Maud,
L. Moscadelli,
T. Möller,
A. Palau,
T. Peters
, et al. (7 additional authors not shown)
Abstract:
Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk ar…
▽ More
Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources.
△ Less
Submitted 10 September, 2021;
originally announced September 2021.
-
Fragmentation and kinematics in high-mass star formation: CORE-extension targeting two very young high-mass star-forming regions
Authors:
H. Beuther,
C. Gieser,
S. Suri,
H. Linz,
P. Klaassen,
D. Semenov,
J. M. Winters,
Th. Henning,
J. D. Soler,
J. S. Urquhart,
J. Syed,
S . Feng,
T. Moeller,
M. T. Beltran,
A. Sanchez-Monge,
S. N. Longmore,
T. Peters,
J. Ballesteros-Paredes,
P. Schilke,
L. Moscadelli,
A. Palau,
R. Cesaroni,
S. Lumsden,
R. Pudritz,
F. Wyrowski
, et al. (2 additional authors not shown)
Abstract:
Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we…
▽ More
Context: The formation of high-mass star-forming regions from their parental gas cloud and the subsequent fragmentation processes lie at the heart of star formation research. Aims: We aim to study the dynamical and fragmentation properties at very early evolutionary stages of high-mass star formation. Methods: Employing the NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30m telescope, we observed two young high-mass star-forming regions, ISOSS22478 and ISOSS23053, in the 1.3mm continuum and spectral line emission at a high angular resolution (~0.8''). Results: We resolved 29 cores that are mostly located along filament-like structures. Depending on the temperature assumption, these cores follow a mass-size relation of approximately M~r^2.0, corresponding to constant mean column densities. However, with different temperature assumptions, a steeper mass-size relation up to M~r^3.0, which would be more likely to correspond to constant mean volume densities, cannot be ruled out. The correlation of the core masses with their nearest neighbor separations is consistent with thermal Jeans fragmentation. We found hardly any core separations at the spatial resolution limit, indicating that the data resolve the large-scale fragmentation well. Although the kinematics of the two regions appear very different at first sight - multiple velocity components along filaments in ISOSS22478 versus a steep velocity gradient of more than 50km/s/pc in ISOSS23053 - the findings can be explained within the framework of a dynamical cloud collapse scenario. Conclusions: While our data are consistent with a dynamical cloud collapse scenario and subsequent thermal Jeans fragmentation, the importance of additional environmental properties, such as the magnetization of the gas or external shocks triggering converging gas flows, is nonetheless not as well constrained and would require future investigation.
△ Less
Submitted 6 April, 2021;
originally announced April 2021.
-
The physical and chemical structure of high-mass star-forming regions. Unraveling chemical complexity with the NOEMA large program "CORE"
Authors:
C. Gieser,
H. Beuther,
D. Semenov,
A. Ahmadi,
S. Suri,
T. Möller,
M. T. Beltran,
P. Klaassen,
Q. Zhang,
J. S. Urquhart,
Th. Henning,
S. Feng,
R. Galván-Madrid,
V. de Souza Magalhães,
L. Moscadelli,
S. Longmore,
S. Leurini,
R. Kuiper,
T. Peters,
K. M. Menten,
T. Csengeri,
G. Fuller,
F. Wyrowski,
S. Lumsden,
Á. Sánchez-Monge
, et al. (8 additional authors not shown)
Abstract:
We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profil…
▽ More
We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again.
△ Less
Submitted 23 February, 2021;
originally announced February 2021.
-
Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind
Authors:
L. Moscadelli,
H. Beuther,
A. Ahmadi,
C. Gieser,
F. Massi,
R. Cesaroni,
Á. Sánchez-Monge,
F. Bacciotti,
M. T. Beltrán,
T. Csengeri,
R. Galván-Madrid,
Th. Henning,
P. D. Klaassen,
R. Kuiper,
S. Leurini,
S. N. Longmore,
L. T. Maud,
T. Möller,
A. Palau,
T. Peters,
R. E. Pudritz,
A. Sanna,
D. Semenov,
J. S. Urquhart,
J. M. Winters
, et al. (1 additional authors not shown)
Abstract:
In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we deri…
▽ More
In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively.
△ Less
Submitted 9 February, 2021;
originally announced February 2021.
-
CRAHCN-O: A Consistent Reduced Atmospheric Hybrid Chemical Network Oxygen Extension for Hydrogen Cyanide and Formaldehyde Chemistry in CO2-, N2-, H2O-, CH4-, and H2-Dominated Atmospheres
Authors:
Ben K. D. Pearce,
Paul W. Ayers,
Ralph E. Pudritz
Abstract:
Hydrogen cyanide (HCN) and formaldehyde (H2CO) are key precursors to biomolecules such as nucleobases and amino acids in planetary atmospheres; However, many reactions which produce and destroy these species in atmospheres containing CO2 and H2O are still missing from the literature. We use a quantum chemistry approach to find these missing reactions and calculate their rate coefficients using can…
▽ More
Hydrogen cyanide (HCN) and formaldehyde (H2CO) are key precursors to biomolecules such as nucleobases and amino acids in planetary atmospheres; However, many reactions which produce and destroy these species in atmospheres containing CO2 and H2O are still missing from the literature. We use a quantum chemistry approach to find these missing reactions and calculate their rate coefficients using canonical variational transition state theory and Rice-Ramsperger-Kassel-Marcus/master equation theory at the BHandHLYP/aug-cc-pVDZ level of theory. We calculate the rate coefficients for 126 total reactions, and validate our calculations by comparing with experimental data in the 39% of available cases. Our calculated rate coefficients are most frequently within an factor of 2 of experimental values, and generally always within an order of magnitude of these values. We discover 45 previously unknown reactions, and identify 6 from this list that are most likely to dominate H2CO and HCN production and destruction in planetary atmospheres. We highlight $^1$O + CH3 $\rightarrow$ H2CO + H as a new key source, and H2CO + $^1$O $\rightarrow$ HCO + OH as a new key sink, for H2CO in upper planetary atmospheres. In this effort, we develop an oxygen extension to our consistent reduced atmospheric hybrid chemical network (CRAHCN-O), building off our previously developed network for HCN production in N2-, CH4-, and H2-dominated atmospheres (CRAHCN). This extension can be used to simulate both HCN and H2CO production in atmospheres dominated by any of CO2, N2, H2O, CH4, and H2.
△ Less
Submitted 25 September, 2020;
originally announced September 2020.
-
Connecting planet formation and astrochemistry: C/O and N/O of warm giant planets and Jupiter-analogs
Authors:
Alex J. Cridland,
Ewine F. van Dishoeck,
Matthew Alessi,
Ralph E. Pudritz
Abstract:
(Abridged) The chemical composition of planetary atmospheres has long been thought to store information regarding where and when a planet accretes its material. Predicting this chemical composition theoretically is a crucial step in linking observational studies to the underlying physics that govern planet formation. As a follow-up to a study of hot Jupiters in our previous work, we present a popu…
▽ More
(Abridged) The chemical composition of planetary atmospheres has long been thought to store information regarding where and when a planet accretes its material. Predicting this chemical composition theoretically is a crucial step in linking observational studies to the underlying physics that govern planet formation. As a follow-up to a study of hot Jupiters in our previous work, we present a population of warm Jupiters (semi-major axis between 0.5-4 AU) extracted from the same planetesimal formation population synthesis model as used in our previous work. We compute the astrochemical evolution of the protoplanetary disks included in this population to predict the carbon-to-oxygen (C/O) and nitrogen-to-oxygen (N/O) ratio evolution of the disk gas, ice, and refractory sources, the accretion of which greatly impacts the resulting C/O and N/O in the atmosphere of giant planets. We confirm that the main sequence (between accreted solid mass and atmospheric C/O) we found previously is largely reproduced by the presented population of synthetic warm Jupiters. And as a result, the majority of the population fall along the empirically derived mass-metallicity relation when the natal disk has solar or lower metallicity. Planets forming from disks with high metallicity ([Fe/H] $>$ 0.1) result in more scatter in chemical properties which could explain some of the scatter found in the mass-metallicity relation. Combining predicted C/O and N/O ratios shows that Jupiter does not fall among our population of synthetic planets, suggesting that it likely did not form in the inner 5 AU of the solar system before proceeding into a Grand Tack. This result is consistent with recent analysis of the chemical composition of Jupiter's atmosphere which suggests that it accreted most of its heavy element abundance farther than tens of AU away from the Sun.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
HCN production in Titan's Atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling
Authors:
Ben K. D. Pearce,
Karan Molaverdikhani,
Ralph E. Pudritz,
Thomas Henning,
Eric Hébrard
Abstract:
Hydrogen cyanide (HCN) is a critical reactive source of nitrogen for building key biomolecules relevant for the origin of life. Still, many HCN reactions remain uncharacterized by experiments and theory, and the complete picture of HCN production in planetary atmospheres is not fully understood. To improve this situation, we develop a novel technique making use of computational quantum chemistry,…
▽ More
Hydrogen cyanide (HCN) is a critical reactive source of nitrogen for building key biomolecules relevant for the origin of life. Still, many HCN reactions remain uncharacterized by experiments and theory, and the complete picture of HCN production in planetary atmospheres is not fully understood. To improve this situation, we develop a novel technique making use of computational quantum chemistry, experimental data, and atmospheric numerical simulations. First, we use quantum chemistry simulations to explore the entire field of possible reactions for a list of primary species in N2-, CH4-, and H2-dominated atmospheres. In this process, we discover 33 new reactions with no previously known rate coefficients. From here, we develop a consistent reduced atmospheric hybrid chemical network (CRAHCN) containing experimental values when available, and our calculated rate coefficients otherwise. Next, we couple CRAHCN to a 1D chemical kinetic model (ChemKM) to compute the HCN abundance as a function of atmospheric depth on Titan. Our simulated atmospheric HCN profile agrees very well with the Cassini observations. CRAHCN contains 104 reactions however nearly all of the simulated atmospheric HCN profile can be obtained using a scaled down network of only 19 dominant reactions. From here, we form a complete picture of HCN chemistry in Titan's atmosphere, from the dissociation of the main atmospheric species, down to the direct production of HCN along 4 major channels. One of these channels was first discovered and characterized in Pearce et al. (2019) and this work.
△ Less
Submitted 10 August, 2020;
originally announced August 2020.
-
Formation of Planetary Populations III: Core Composition & Atmospheric Evaporation
Authors:
Matthew Alessi,
Julie Inglis,
Ralph E. Pudritz
Abstract:
The exoplanet mass radius diagram reveals that super Earths display a wide range of radii, and therefore mean densities, at a given mass. Using planet population synthesis models, we explore the key physical factors that shape this distribution: planets' solid core compositions, and their atmospheric structure. For the former, we use equilibrium disk chemistry models to track accreted minerals ont…
▽ More
The exoplanet mass radius diagram reveals that super Earths display a wide range of radii, and therefore mean densities, at a given mass. Using planet population synthesis models, we explore the key physical factors that shape this distribution: planets' solid core compositions, and their atmospheric structure. For the former, we use equilibrium disk chemistry models to track accreted minerals onto planetary cores throughout formation. For the latter, we track gas accretion during formation, and consider photoevaporation-driven atmospheric mass loss to determine what portion of accreted gas escapes after the disk phase. We find that atmospheric stripping of Neptunes and sub-Saturns at small orbital radii ($\lesssim$0.1AU) plays a key role in the formation of short-period super Earths. Core compositions are strongly influenced by the trap in which they formed. We also find a separation between Earth-like planet compositions at small orbital radii $\lesssim$0.5AU and ice-rich planets (up to 50\% by mass) at larger orbits $\sim$1AU. This corresponds well with the Earth-like mean densities inferred from the observed position of the low-mass planet radius valley at small orbital periods. Our model produces planet radii comparable to observations at masses $\sim$1-3M$_\oplus$. At larger masses, planets' accreted gas significantly increases their radii to be larger than most of the observed data. While photoevaporation, affecting planets at small orbital radii $\lesssim$0.1AU, reduces a subset of these planets' radii and improves our comparison, most planets in our computed populations are unaffected due to low FUV fluxes as they form at larger separations.
△ Less
Submitted 13 July, 2020;
originally announced July 2020.
-
Formation of Planetary Populations II: Effects of Initial Disk Size & Radial Dust Drift
Authors:
Matthew Alessi,
Ralph E. Pudritz,
Alex J. Cridland
Abstract:
Recent ALMA observations indicate that while a range of disk sizes exist, typical disk radii are small, and that radial dust drift affects the distribution of solids in disks. Here we explore the consequences of these features in planet population synthesis models. A key feature of our model is planet traps - barriers to otherwise rapid type-I migration of forming planets - for which we include th…
▽ More
Recent ALMA observations indicate that while a range of disk sizes exist, typical disk radii are small, and that radial dust drift affects the distribution of solids in disks. Here we explore the consequences of these features in planet population synthesis models. A key feature of our model is planet traps - barriers to otherwise rapid type-I migration of forming planets - for which we include the ice line, heat transition, and outer edge of the dead zone. We find that the ice line plays a fundamental role in the formation of warm Jupiters. In particular, the ratio of super Earths to warm Jupiters formed at the ice line depend sensitively on the initial disk radius. Initial gas disk radii of $\sim$50 AU results in the largest super Earth populations, while both larger and smaller disk sizes result in the ice line producing more gas giants near 1 AU. This transition between typical planet class formed at the ice line at various disk radii confirms that planet formation is fundamentally linked to disk properties (in this case, disk size), and is a result that is only seen when dust evolution effects are included in our models. Additionally, we find that including radial dust drift results in the formation of more super Earths between 0.1 - 1 AU, having shorter orbital radii than those produced in models where dust evolution effects are not included.
△ Less
Submitted 28 January, 2020;
originally announced January 2020.
-
The Role of Magnetic Fields in Protostellar Outflows and Star Formation
Authors:
Ralph E. Pudritz,
Tom P. Ray
Abstract:
The role of outflows in the formation of stars and the protostellar disks that generate them is a central question in astrophysics. Outflows are associated with star formation across the entire stellar mass spectrum. In this review, we describe the observational, theoretical, and computational advances on magnetized outflows, and their role in the formation of disks and stars of all masses in turb…
▽ More
The role of outflows in the formation of stars and the protostellar disks that generate them is a central question in astrophysics. Outflows are associated with star formation across the entire stellar mass spectrum. In this review, we describe the observational, theoretical, and computational advances on magnetized outflows, and their role in the formation of disks and stars of all masses in turbulent, magnetized clouds. The ability of torques exerted on disks by magnetized winds to efficiently extract and transport disk angular momentum was developed in early theoretical models and confirmed by a variety of numerical simulations. The recent high resolution ALMA observations of disks and outflows now confirm several key aspects of these ideas, e.g. that jets rotate and originate from large regions of their underlying disks. New insights on accretion disk physics show that magneto-rotational instability (MRI) turbulence is strongly damped, leaving magnetized disk winds as the dominant mechanism for transporting disk angular momentum. This has major consequences for star formation, as well as planet formation. Outflows also play an important role in feedback processes particularly in the birth of low mass stars and cluster formation. Despite being almost certainly fundamental to their production and focusing, magnetic fields in outflows in protostellar systems, and even in the disks, are notoriously difficult to measure. Most methods are indirect and lack precision, as for example, when using optical/near-infrared line ratios. Moreover, in those rare cases where direct measurements are possible - where synchrotron radiation is observed, one has to be very careful in interpreting derived values. Here we also explore what is known about magnetic fields from observations, and take a forward look to the time when facilities such as SPIRou and the SKA are in routine operation.
△ Less
Submitted 11 December, 2019;
originally announced December 2019.
-
The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks
Authors:
James Di Francesco,
Helen Kirk,
Doug Johnstone,
Ralph Pudritz,
Shantanu Basu,
Sarah Sadavoy,
Laura Fissel,
Lewis Knee,
Mehrnoosh Tahani,
Rachel Friesen,
Simon Coudé,
Erik Rosolowsky,
Nienke van der Marel,
Michel Fich,
Christine Wilson,
Chris Matzner,
Ruobing Dong,
Brenda Matthews,
Gerald Schieven
Abstract:
Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length sc…
▽ More
Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.
△ Less
Submitted 9 November, 2019; v1 submitted 5 November, 2019;
originally announced November 2019.
-
Connecting planet formation and astrochemistry: A main sequence for C/O in hot-exoplanetary atmospheres
Authors:
Alex J. Cridland,
Ewine F. van Dishoeck,
Matthew Alessi,
Ralph E. Pudritz
Abstract:
To understand the role that planet formation history has on the observable atmospheric carbon-to-oxygen ratio (C/O) we have produced a population of astrochemically evolving protoplanetary disks. Based on the parameters used in a pre-computed population of growing planets their combination allows us to trace the molecular abundances of the gas that is being collected into planetary atmospheres. We…
▽ More
To understand the role that planet formation history has on the observable atmospheric carbon-to-oxygen ratio (C/O) we have produced a population of astrochemically evolving protoplanetary disks. Based on the parameters used in a pre-computed population of growing planets their combination allows us to trace the molecular abundances of the gas that is being collected into planetary atmospheres. We include atmospheric pollution of incoming (icy) planetesimals as well as the effect of refractory carbon erosion noted to exist in our own solar system. We find that the carbon and oxygen content of Neptune-mass planets are determined primarily through solid accretion and result in more oxygen-rich (by roughly two orders of magnitude) atmospheres than Hot-Jupiters, whose C/O are primarily determined by gas accretion. Generally we find a `main-sequence' between the fraction of planetary mass accreted through solid accretion and the resulting atmospheric C/O - with planets of higher solid accretion fraction having lower C/O. Hot-Jupiters whose atmospheres have been chemically characterized agree well with our population of planets, and our results suggest that Hot Jupiter formation typically begins near the water ice line. Lower mass hot-Neptunes are observed to be much more carbon-rich (with 0.33 $\lesssim$ C/O $\lesssim$ 1) than is found in our models (C/O $\sim 10^{-2}$), and suggest that some form of chemical processing may affect their observed C/O over the few Gyrs between formation and observation. Our population reproduces the general mass-metallicity trend of the solar system and qualitatively reproduces the C/O-metallicity anti-correlation that has been inferred for the population of characterized exoplanetary atmospheres.
△ Less
Submitted 29 October, 2019;
originally announced October 2019.
-
Chemical complexity in high-mass star formation: An observational and modeling case study of the AFGL 2591 VLA 3 hot core
Authors:
C. Gieser,
D. Semenov,
H. Beuther,
A. Ahmadi,
J. C. Mottram,
Th. Henning,
M. Beltran,
L. T. Maud,
F. Bosco,
S. Leurini,
T. Peters,
P. Klaassen,
R. Kuiper,
S. Feng,
J. S. Urquhart,
L. Moscadelli,
T. Csengeri,
S. Lumsden,
J. M. Winters,
S. Suri,
Q. Zhang,
R. Pudritz,
A. Palau,
K. M. Menten,
R. Galvan-Madrid
, et al. (8 additional authors not shown)
Abstract:
We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observ…
▽ More
We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.
△ Less
Submitted 11 October, 2019;
originally announced October 2019.
-
LRP2020: Signposts of planet formation in protoplanetary disks
Authors:
Nienke van der Marel,
Ruobing Dong,
Ralph Pudritz,
James Wadsley,
Aaron Boley,
Eve Lee,
Mohamad Ali-Dib,
Brenda Matthews,
Christian Marois,
Henry Ngo
Abstract:
Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary…
▽ More
Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.
△ Less
Submitted 30 September, 2019;
originally announced October 2019.
-
Global dynamics of the interstellar medium in magnetised disc galaxies
Authors:
Bastian Körtgen,
Robi Banerjee,
Ralph E. Pudritz,
Wolfram Schmidt
Abstract:
Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stron…
▽ More
Magnetic fields are an elemental part of the interstellar medium in galaxies. However, their impact on gas dynamics and star formation in galaxies remains controversial. We use a suite of global magnetohydrodynamical simulations of isolated disc galaxies to study the influence of magnetic fields on the diffuse and dense gas in the discs. We find that the magnetic field acts in multiple ways. Stronger magnetised discs fragment earlier due to the shorter growth time of the Parker instability. Due to the Parker instability in the magnetised discs we also find cold ($T<50\,\mathrm{K}$) and dense ($n\sim10^3-10^4\,\mathrm{cm}^{-3}$) gas several hundred pc above/below the midplane without any form of stellar feedback. In addition, magnetic fields change the fragmentation pattern. While in the hydrodynamical case, the disc breaks up into ring-like structures, magnetised discs show the formation of filamentary entities that extent both in the azimuthal and radial direction. These kpc scale filaments become magnetically (super-)critical very quickly and allow for the rapid formation of massive giant molecular clouds. Our simulations suggest that major differences in the behaviour of star formation - due to a varying magnetisation - in galaxies could arise.
△ Less
Submitted 4 September, 2019;
originally announced September 2019.
-
Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE
Authors:
F. Bosco,
H. Beuther,
A. Ahmadi,
J. C. Mottram,
R. Kuiper,
H. Linz,
L. Maud,
J. M. Winters,
T. Henning,
S. Feng,
T. Peters,
D. Semenov,
P. D. Klaassen,
P. Schilke,
J. S. Urquhart,
M. T. Beltrán,
S. L. Lumsden,
S. Leurini,
L. Moscadelli,
R. Cesaroni,
Á. Sánchez-Monge,
A. Palau,
R. Pudritz,
F. Wyrowski,
S. Longmore
, et al. (1 additional authors not shown)
Abstract:
The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. U…
▽ More
The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.
△ Less
Submitted 9 July, 2019;
originally announced July 2019.
-
IRAS23385+6053: An embedded massive cluster in the making
Authors:
R. Cesaroni,
H. Beuther,
A. Ahmadi,
M. T. Beltran,
T. Csengeri,
R. Galvan-Madrid,
C. Gieser,
T. Henning,
K. G. Johnston,
P. D. Klaassen,
R. Kuiper,
S. Leurini,
H. Linz,
S. Longmore,
S. L. Lumsden,
L. T. Maud,
L. Moscadelli,
J. C. Mottram,
A. Palau,
T. Peters,
R. E. Pudritz,
A. Sanchez-Monge,
P. Schilke,
D. Semenov,
S. Suri
, et al. (4 additional authors not shown)
Abstract:
This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additi…
▽ More
This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.
△ Less
Submitted 27 May, 2019;
originally announced May 2019.
-
A Consistent Reduced Network for HCN Chemistry in Early Earth and Titan Atmospheres: Quantum Calculations of Reaction Rate Coefficients
Authors:
Ben K. D. Pearce,
Paul W. Ayers,
Ralph E. Pudritz
Abstract:
HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate…
▽ More
HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate rate coefficients at 298 K using canonical variational transition state theory (CVT) paired with computational quantum chemistry simulations at the BHandHLYP/augcc-pVDZ level of theory. We also calculate the temperature dependence of the rate coefficients for the reactions that have barriers from 50 to 400 K. We present 15 new reaction rate coefficients with no previously known value; 93% of our calculated coefficients are within an order of magnitude of the nearest experimental or recommended values. Above 320 K, the rate coefficient for the new reaction H2CN -> HCN + H dominates. Contrary to experiments, we find the HCN reaction pathway, N + CH3 -> HCN + H2, to be inefficient and suggest that the experimental rate coefficient actually corresponds to an indirect pathway, through the H2CN intermediate. We present CVT using energies computed with density functional theory as a feasible and accurate method for calculating a large network of rate coefficients of small-molecule reactions.
△ Less
Submitted 14 February, 2019;
originally announced February 2019.
-
Physics of Planet Trapping with Applications to HL Tau
Authors:
Alex J. Cridland,
Ralph E. Pudritz,
Matthew Alessi
Abstract:
We explore planet formation in the HL Tau disk and possible origins of the prominent gaps and rings observed by ALMA. We investigate whether dust gaps are caused by dynamically trapped planetary embryos at the ice lines of abundant volatiles. The global properties of the HL Tau disk (total mass, size) at its current age are used to constrain an evolving analytic disk model describing its temperatu…
▽ More
We explore planet formation in the HL Tau disk and possible origins of the prominent gaps and rings observed by ALMA. We investigate whether dust gaps are caused by dynamically trapped planetary embryos at the ice lines of abundant volatiles. The global properties of the HL Tau disk (total mass, size) at its current age are used to constrain an evolving analytic disk model describing its temperature and density profiles. By performing a detailed analysis of the planet-disk interaction for a planet near the water ice line including a rigorous treatment of the dust opacity, we confirm that water is sufficiently abundant (1.5x10^-4 molecules per H) to trap planets at its ice line due to an opacity transition. When the abundance of water is reduced by 50% planet trapping disappears. We extend our analysis to other planet traps: the heat transition, dead zone edge, and the CO_2 ice line and find similar trapping. The formation of planets via planetesimal accretion is computed for dynamically trapped embryos at the water ice line, dead zone, and heat transition. The end products orbit in the inner disk (R < 3 AU), unresolved by ALMA, with masses that range between sub-Earth to 5 Jupiter masses. While we find that the dust gaps correspond well with the radial positions of the CO_2 , CH_4 , and CO ice lines, the planetesimal accretion rates at these radii are too small to build large embryos within 1 Myr.
△ Less
Submitted 3 January, 2019;
originally announced January 2019.
-
Meteoritic Abundances of Fatty Acids and Potential Reaction Pathways in Planetesimals
Authors:
James C. -Y. Lai,
Ben K. D. Pearce,
Ralph E. Pudritz,
Drake Lee
Abstract:
The origin of fatty acids on the prebiotic Earth is important as they likely formed the encapsulating membranes of the first protocells. Carbon-rich meteorites (i.e., carbonaceous chondrites) such as Murchison and Tagish Lake are well known to contain these molecules, and their delivery to the early planet by intense early meteorite bombardments constitutes a key prebiotic source. We collect the f…
▽ More
The origin of fatty acids on the prebiotic Earth is important as they likely formed the encapsulating membranes of the first protocells. Carbon-rich meteorites (i.e., carbonaceous chondrites) such as Murchison and Tagish Lake are well known to contain these molecules, and their delivery to the early planet by intense early meteorite bombardments constitutes a key prebiotic source. We collect the fatty acid abundances measured in various carbonaceous chondrites from the literature and analyze them for patterns and correlations. Fatty acids in meteorites include straight-chain and branched-chain monocarboxylic and dicarboxylic acids up to 12 carbons in length---fatty acids with at least 8 carbons are required to form vesicles, and modern cell membranes employ lipids with ~12--20 carbons. To understand the origin of meteoritic fatty acids, we search the literature for abiotic fatty acid reaction pathways and create a candidate list of 11 reactions that could potentially produce these fatty acids in meteorite parent bodies. Straight-chain monocarboxylic acids (SCMA) are the dominant fatty acids in meteorites, followed by branched-chain monocarboxylic acids (BCMA). SCMA are most abundant in CM2 and Tagish Lake (ungrouped) meteorites, ranging on average from 10$^2$ ppb to 4x10$^5$ ppb, and 10$^4$ ppb to 5x10$^6$ ppb, respectively. In CM, CV, and Tagish Lake meteorites, SCMA abundances generally decrease with increasing carbon chain length. Conversely, SCMA abundances in CR meteorites peak at 5 and 6 carbons in length, and decrease on either side of this peak. This unique CR fatty acid distribution may hint at terrestrial contamination, or that certain fatty acid reactions mechanisms are active in different meteorite parent bodies (planetesimals). We identify Fischer-Tropsch-type synthesis as the most promising pathway for further analysis in the production of fatty acids in planetesimals.
△ Less
Submitted 25 September, 2018;
originally announced September 2018.
-
Constraining the Time Interval for the Origin of Life on Earth
Authors:
Ben K. D. Pearce,
Andrew S. Tupper,
Ralph E. Pudritz,
Paul G. Higgs
Abstract:
Estimates of the time at which life arose on Earth make use of two types of evidence. First, astrophysical and geophysical studies provide a timescale for the formation of Earth and the Moon, for large impact events on early Earth, and for the cooling of the early magma ocean. From this evidence, we can deduce a habitability boundary, which is the earliest point at which Earth became habitable. Se…
▽ More
Estimates of the time at which life arose on Earth make use of two types of evidence. First, astrophysical and geophysical studies provide a timescale for the formation of Earth and the Moon, for large impact events on early Earth, and for the cooling of the early magma ocean. From this evidence, we can deduce a habitability boundary, which is the earliest point at which Earth became habitable. Second, biosignatures in geological samples, including microfossils, stromatolites, and chemical isotope ratios, provide evidence for when life was actually present. From these observations we can deduce a biosignature boundary, which is the earliest point at which there is clear evidence that life existed. Studies with molecular phylogenetics and records of the changing level of oxygen in the atmosphere give additional information that helps to determine the biosignature boundary. Here, we review the data from a wide range of disciplines to summarize current information on the timings of these two boundaries. The habitability boundary could be as early as 4.5 Ga, the earliest possible estimate of the time at which Earth had a stable crust and hydrosphere, or as late as 3.9 Ga, the end of the period of heavy meteorite bombardment. The lack of consensus on whether there was a late heavy meteorite bombardment that was significant enough to prevent life is the largest uncertainty in estimating the time of the habitability boundary. The biosignature boundary is more closely constrained. Evidence from carbon isotope ratios and stromatolite fossils both point to a time close to 3.7 Ga. Life must have emerged in the interval between these two boundaries. The time taken for life to appear could, therefore, be within 200 Myr or as long as 800 Myr.
△ Less
Submitted 28 August, 2018;
originally announced August 2018.
-
On the Origin of Multiple Populations During Massive Star Cluster Formation
Authors:
Corey S. Howard,
Ralph E. Pudritz,
Alison Sills,
William E. Harris
Abstract:
We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accret…
▽ More
We investigate the possibility that multiple populations in globular clusters arise as a natural by-product of massive star-cluster formation. We use 3D radiative hydrodynamics simulations for the formation of young massive clusters to track their chemical self-enrichment during their first 5 Myr. These clusters form embedded within filamentary Giant Molecular Clouds by a combination of gas accretion and rapid merging of protoclusters. Chemical enrichment is a dynamic process happening as the young cluster assembles, so that the original (1P) and enriched (2P) subpopulations of stars form almost simultaneously. Here we test two simple and opposite extremes for the injection of enriched material into the intracluster gas: we assume either continuous injection in a way that tracks the star formation rate; or sudden injection by a single instantaneous event. Using helium abundance as a proxy for the enrichment, we find that realistic multiple population features can be reproduced by injecting a total helium mass amounting to a few percent of the cluster's total mass. The differences in individual growth histories can lead to widely differing 1P/2P outcomes. These models suggest that dual or multiple populations can emerge rapidly in massive star clusters undergoing the typical mode of star cluster formation.
△ Less
Submitted 5 April, 2019; v1 submitted 21 August, 2018;
originally announced August 2018.
-
A universal route for the formation of massive star clusters in giant molecular clouds
Authors:
Corey S. Howard,
Ralph E. Pudritz,
William E. Harris
Abstract:
Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show tha…
▽ More
Young massive star clusters (YMCs, with M $\geq$10$^4$ M$_{\odot}$) are proposed modern-day analogues of the globular clusters (GCs) that were products of extreme star formation in the early universe. The exact conditions and mechanisms under which YMCs form remain unknown -- a fact further complicated by the extreme radiation fields produced by their numerous massive young stars. Here we show that GC-sized clusters are naturally produced in radiation-hydrodynamic simulations of isolated 10$^7$ M$_{\odot}$ Giant Molecular Clouds (GMCs) with properties typical of the local universe, even under the influence of radiative feedback. In all cases, these massive clusters grow to GC-level masses within 5 Myr via a roughly equal combination of filamentary gas accretion and mergers with several less massive clusters. Lowering the heavy-element abundance of the GMC by a factor of 10 reduces the opacity of the gas to radiation and better represents the high-redshift formation conditions of GCs. This results in higher gas accretion leading to a mass increase of the largest cluster by a factor of ~4. When combined with simulations of less massive GMCs (10$^{4-6}$ M$_{\odot}$), a clear relation emerges between the maximum YMC mass and the mass of the host GMC. Our results demonstrate that YMCs, and potentially GCs, are a simple extension of local cluster formation to more massive clouds and do not require suggested exotic formation scenarios.
△ Less
Submitted 21 August, 2018;
originally announced August 2018.
-
Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE
Authors:
A. Ahmadi,
H. Beuther,
J. C. Mottram,
F. Bosco,
H. Linz,
Th. Henning,
J. M. Winters,
R. Kuiper,
R. Pudritz,
Á. Sánchez-Monge,
E. Keto,
M. Beltran,
S. Bontemps,
R. Cesaroni,
T. Csengeri,
S. Feng,
R. Galvan-Madrid,
K. G. Johnston,
P. Klaassen,
S. Leurini,
S. N. Longmore,
S. Lumsden,
L. T. Maud,
K. M. Menten,
L. Moscadelli
, et al. (8 additional authors not shown)
Abstract:
The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.3…
▽ More
The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.
△ Less
Submitted 1 August, 2018;
originally announced August 2018.
-
The Origin of Filamentary Star Forming Clouds in Magnetised Galaxies
Authors:
Bastian Körtgen,
Robi Banerjee,
Ralph E. Pudritz,
Wolfram Schmidt
Abstract:
Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from…
▽ More
Observations show that galaxies and their interstellar media are pervaded by strong magnetic fields with energies in the diffuse component being at least comparable to the thermal and even as large or larger than the turbulent energy. Such strong magnetic fields prevent the formation of stars because patches of the interstellar medium are magnetically subcritical. Here we present the results from global numerical simulations of strongly magnetised and self-gravitating galactic discs, which show that the buoyancy of the magnetic field due to the Parker instability leads at first to the formation of giant filamentary regions. These filamentary structures become gravitationally unstable and fragment into $\sim10^5 M_{\odot}$ clouds that attract kpc long, coherent filamentary flows that build them into GMCs. Our results thus provide a solution to the long-standing problem of how the transition from sub- to supercritical regions in the interstellar medium proceeds.
△ Less
Submitted 22 May, 2018;
originally announced May 2018.
-
Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE
Authors:
H. Beuther,
J. C. Mottram,
A. Ahmadi,
F. Bosco,
H. Linz,
Th. Henning,
P. Klaassen,
J. M. Winters,
L. T. Maud,
R. Kuiper,
D. Semenov,
C. Gieser,
T. Peters,
J. S. Urquhart,
R. Pudritz,
S. E. Ragan,
S. Feng,
E. Keto,
S. Leurini,
R. Cesaroni,
M. Beltran,
A. Palau,
A. Sanchez-Monge,
R. Galvan-Madrid,
Q. Zhang
, et al. (8 additional authors not shown)
Abstract:
Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores.
Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high…
▽ More
Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores.
Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun.
Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths.
Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.
△ Less
Submitted 3 May, 2018;
originally announced May 2018.
-
Connecting Planetary Composition with Formation
Authors:
Ralph E. Pudritz,
Alex J. Cridland,
Matthew Alessi
Abstract:
The rapid advances in observations of the different populations of exoplanets, the characterization of their host stars and the links to the properties of their planetary systems, the detailed studies of protoplanetary disks, and the experimental study of the interiors and composition of the massive planets in our solar system provide a firm basis for the next big question in planet formation theo…
▽ More
The rapid advances in observations of the different populations of exoplanets, the characterization of their host stars and the links to the properties of their planetary systems, the detailed studies of protoplanetary disks, and the experimental study of the interiors and composition of the massive planets in our solar system provide a firm basis for the next big question in planet formation theory. How do the elemental and chemical compositions of planets connect with their formation? The answer to this requires that the various pieces of planet formation theory be linked together in an end-to-end picture that is capable of addressing these large data sets. In this review, we discuss the critical elements of such a picture and how they affect the chemical and elemental make up of forming planets. Important issues here include the initial state of forming and evolving disks, chemical and dust processes within them, the migration of planets and the importance of planet traps, the nature of angular momentum transport processes involving turbulence and/or MHD disk winds, planet formation theory, and advanced treatments of disk astrochemistry. All of these issues affect, and are affected by the chemistry of disks which is driven by X-ray ionization of the host stars. We discuss how these processes lead to a coherent end-to-end model and how this may address the basic question.
△ Less
Submitted 11 April, 2018;
originally announced April 2018.
-
Formation of Planetary Populations I: Metallicity & Envelope Opacity Effects
Authors:
Matthew Alessi,
Ralph E. Pudritz
Abstract:
We present a comprehensive body of simulations of the formation of exoplanetary populations that incorporate the role of planet traps in slowing planetary migration. The traps we include in our model are the water ice line, the disk heat transition, and the dead zone outer edge. We reduce our model parameter set to two physical parameters: the opacity of the accreting planetary atmospheres (…
▽ More
We present a comprehensive body of simulations of the formation of exoplanetary populations that incorporate the role of planet traps in slowing planetary migration. The traps we include in our model are the water ice line, the disk heat transition, and the dead zone outer edge. We reduce our model parameter set to two physical parameters: the opacity of the accreting planetary atmospheres ($κ_{\rm{env}}$) and a measure of the efficiency of planetary accretion after gap opening ($f_{\rm{max}}$). We perform planet population synthesis calculations based on the initial observed distributions of host star and disk properties - their disk masses, lifetimes, and stellar metallicities. We find the frequency of giant planet formation scales with disk metallicity, in agreement with the observed Jovian planet frequency-metallicity relation. We consider both X-ray and cosmic ray disk ionization models, whose differing ionization rates lead to different dead zone trap locations. In both cases, Jovian planets form in our model out to 2-3 AU, with a distribution at smaller radii dependent on the disk ionization source and the setting of envelope opacity. We find that low values of $κ_{\rm{env}}$ (0.001-0.002 cm$^2$ g$^{-1}$) and X-ray disk ionization are necessary to obtain a separation between hot Jupiters near 0.1 AU, and warm Jupiters outside 0.6 AU, a feature present in the data. Our model also produces a large number of super Earths, but the majority are outside of 2 AU. As our model assumes a constant dust to gas ratio, we suggest that radial dust evolution must be taken into account to reproduce the observed super Earth population.
△ Less
Submitted 2 May, 2018; v1 submitted 3 April, 2018;
originally announced April 2018.
-
Simulating the UV Escape Fractions from Molecular Cloud Populations in Star-forming Dwarf and Spiral Galaxies
Authors:
Corey S. Howard,
Ralph E. Pudritz,
Bill E. Harris,
Ralf S. Klessen
Abstract:
The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) --- Giant Molecular Clouds (GMCs) --- is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, f$_{\text{esc,GMC}}$, from a suite of individual GMC simulations with masses in the range 10$^{4-6}$ M…
▽ More
The escape of ultraviolet photons from the densest regions of the interstellar medium (ISM) --- Giant Molecular Clouds (GMCs) --- is a poorly constrained parameter which is vital to understanding the ionization of the ISM and the intergalactic medium. We characterize the escape fraction, f$_{\text{esc,GMC}}$, from a suite of individual GMC simulations with masses in the range 10$^{4-6}$ M$_{\odot}$ using the adaptive-mesh refinement code FLASH. We find significantly different f$_{\text{esc,GMC}}$ depending on the GMC mass which can reach $>$90% in the evolution of 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ clouds or remain low at $\sim$5% for most of the lifetime of more massive GMCs. All clouds show fluctuations over short, sub-Myr timescales produced by flickering HII regions. We combine our results to calculate the total escape fraction (f$_{\text{esc,tot}}$) from GMC populations in dwarf starburst and spiral galaxies by randomly drawing clouds from a GMC mass distribution (dN/dM$\propto$M$^α$, where $α$ is either -1.5 or -2.5) over fixed time intervals. We find typical f$_{\text{esc,tot}}$ values of 8% for both the dwarf and spiral models. The fluctuations of f$_{\text{esc,tot}}$, however, are much larger for the dwarf models with values as high as 90%. The photons escaping from the 5$\times$10$^4$ and 10$^{5}$ M$_{\odot}$ GMCs are the dominant contributors to f$_{\text{esc,tot}}$ in all cases. We also show that the accompanying star formation rates (SFRs) of our model ($\sim$2$\times$10$^{-2}$ and 0.73 M$_{\odot}$yr$^{-1}$) are consistent with observations of SFRs in dwarf starburst and spiral galaxies, respectively.
△ Less
Submitted 18 December, 2017; v1 submitted 11 October, 2017;
originally announced October 2017.
-
Origin of the RNA World: The Fate of Nucleobases in Warm Little Ponds
Authors:
Ben K. D. Pearce,
Ralph E. Pudritz,
Dmitry A. Semenov,
Thomas K. Henning
Abstract:
Prior to the origin of simple cellular life, the building blocks of RNA (nucleotides) had to form and polymerize in favourable environments on the early Earth. At this time, meteorites and interplanetary dust particles delivered organics such as nucleobases (the characteristic molecules of nucleotides) to warm little ponds whose wet-dry cycles promoted rapid polymerization. We build a comprehensiv…
▽ More
Prior to the origin of simple cellular life, the building blocks of RNA (nucleotides) had to form and polymerize in favourable environments on the early Earth. At this time, meteorites and interplanetary dust particles delivered organics such as nucleobases (the characteristic molecules of nucleotides) to warm little ponds whose wet-dry cycles promoted rapid polymerization. We build a comprehensive numerical model for the evolution of nucleobases in warm little ponds leading to the emergence of the first nucleotides and RNA. We couple Earth's early evolution with complex prebiotic chemistry in these environments. We find that RNA polymers must have emerged very quickly after the deposition of meteorites (< a few years). Their constituent nucleobases were primarily meteoritic in origin and not from interplanetary dust particles. Ponds appeared as continents rose out of the early global ocean but this increasing availability of "targets" for meteorites was offset by declining meteorite bombardment rates. Moreover, the rapid losses of nucleobases to pond seepage during wet periods, and to UV photodissociation during dry periods means that the synthesis of nucleotides and their polymerization into RNA occurred in just one to a few wet-dry cycles. Under these conditions, RNA polymers likely appeared prior to 4.17 billion years ago.
Significance: There are two competing hypotheses for the site at which an RNA world emerged: hydrothermal vents in the deep ocean and warm little ponds. Because the former lacks wet and dry cycles, which are well known to promote polymerization (in this case, of nucleotides into RNA), we construct a comprehensive model for the origin of RNA in the latter sites. Our model advances the story and timeline of the RNA world by constraining the source of biomolecules, the environmental conditions, the timescales of reaction, and the emergence of first RNA polymers.
△ Less
Submitted 1 October, 2017;
originally announced October 2017.
-
Simulating radiative feedback and star cluster formation in GMCs: II. Mass dependence of cloud destruction and cluster properties
Authors:
Corey Howard,
Ralph Pudritz,
William Harris
Abstract:
The process of radiative feedback in Giant Molecular Clouds (GMCs) is an important mechanism for limiting star cluster formation through the heating and ionization of the surrounding gas. We explore the degree to which radiative feedback affects early ($\lesssim$5 Myr) cluster formation in GMCs having masses that range from 10$^{4-6}$ M$_{\odot}$ using the FLASH code. The inclusion of radiative fe…
▽ More
The process of radiative feedback in Giant Molecular Clouds (GMCs) is an important mechanism for limiting star cluster formation through the heating and ionization of the surrounding gas. We explore the degree to which radiative feedback affects early ($\lesssim$5 Myr) cluster formation in GMCs having masses that range from 10$^{4-6}$ M$_{\odot}$ using the FLASH code. The inclusion of radiative feedback lowers the efficiency of cluster formation by 20-50\% relative to hydrodynamic simulations. Two models in particular --- 5$\times$10$^4$ and 10$^5$ M$_{\odot}$ --- show the largest suppression of the cluster formation efficiency, corresponding to a factor of $\sim$2. For these clouds only, the internal energy, a measure of the energy injected by radiative feedback, exceeds the gravitational potential for a significant amount of time. We find a clear relation between the maximum cluster mass, M$_{cl,max}$, formed in a GMC of mass M$_{GMC}$; M$_{cl,max}\propto$ M$_{GMC}^{0.81}$. This scaling result suggests that young globular clusters at the necessary scale of $10^6 M_{\odot}$ form within host GMCs of masses near $\sim 5 \times 10^7 M_{\odot}$. We compare simulated cluster mass distributions to the observed embedded cluster mass function ($dlog(N)/dlog(M) \propto M^β$ where $β$ = -1) and find good agreement ($β$ = -0.99$\pm$0.14) only for simulations including radiative feedback, indicating this process is important in controlling the growth of young clusters. However, the high star formation efficiencies, which range from 16-21\%, and high star formation rates compared to locally observed regions suggest other feedback mechanisms are also important during the formation and growth of stellar clusters.
△ Less
Submitted 1 June, 2017;
originally announced June 2017.
-
Composition of Early Planetary Atmospheres II: Coupled Dust and Chemical Evolution in Protoplanetary Disks
Authors:
A. J. Cridland,
Ralph E. Pudritz,
Tilman Birnstiel,
L. Ilsedore Cleeves,
Edwin A. Bergin
Abstract:
We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter…
▽ More
We present the next step in a series of papers devoted to connecting the composition of the atmospheres of forming planets with the chemistry of their natal evolving protoplanetary disks. The model presented here computes the coupled chemical and dust evolution of the disk and the formation of three planets per disk model. Our three canonical planet traps produce a Jupiter near 1 AU, a Hot Jupiter and a Super-Earth. We study the dependency of the final orbital radius, mass, and atmospheric chemistry of planets forming in disk models with initial disk masses that vary by 0.02 $M_\odot$ above and below our fiducial model ($M_{disk,0} = 0.1 ~M_\odot$). We compute C/O and C/N for the atmospheres formed in our 3 models and find that C/O$_{\rm planet}\sim$ C/O$_{\rm disk}$, which does not vary strongly between different planets formed in our model. The nitrogen content of atmospheres can vary in planets that grow in different disk models. These differences are related to the formation history of the planet, the time and location that the planet accretes its atmosphere, and are encoded in the bulk abundance of NH$_3$. These results suggest that future observations of atmospheric NH$_3$ and an estimation of the planetary C/O and C/N can inform the formation history of particular planetary systems.
△ Less
Submitted 5 May, 2017;
originally announced May 2017.
-
Meteorites and the RNA World: A Thermodynamic Model of Nucleobase Synthesis within Planetesimals
Authors:
Ben K. D. Pearce,
Ralph E. Pudritz
Abstract:
The possible meteorite parent body origin of Earth's pregenetic nucleobases is substantiated by the guanine (G), adenine (A) and uracil (U) measured in various meteorites. Cytosine (C) and thymine (T) however are absent in meteorites, making the emergence of a RNA and later RNA/DNA/protein world problematic. We investigate the meteorite parent body (planetesimal) origin of all nucleobases by compu…
▽ More
The possible meteorite parent body origin of Earth's pregenetic nucleobases is substantiated by the guanine (G), adenine (A) and uracil (U) measured in various meteorites. Cytosine (C) and thymine (T) however are absent in meteorites, making the emergence of a RNA and later RNA/DNA/protein world problematic. We investigate the meteorite parent body (planetesimal) origin of all nucleobases by computationally modeling 18 reactions that potentially contribute to nucleobase formation in such environments. Out of this list, we identify the two most important reactions for each nucleobase and find that these involve small molecules such as HCN, CO, NH3, and water that ultimately arise from the protoplanetary disks in which planetesimals are built. The primary result of this study is that cytosine is unlikely to persist within meteorite parent bodies due to aqueous deamination. Thymine has a thermodynamically favourable reaction pathway from uracil, formaldehyde and formic acid, but likely did not persist within planetesimals containing H2O2 due to an oxidation reaction with this molecule. Finally, while FT synthesis is found to be the dominant source of nucleobases within our model planetesimal, NC synthesis may still be significant under certain chemical conditions (e.g. within CR2 parent bodies). We discuss several major consequences of our results for the origin of the RNA world.
△ Less
Submitted 22 November, 2016;
originally announced November 2016.
-
Radial Drift of Dust in Protoplanetary Disks: The Evolution of Ice lines and Dead zones
Authors:
Alex J. Cridland,
Ralph E. Pudritz,
Tilman Birnstiel
Abstract:
We have developed a new model for the astrochemical structure of a viscously evolving protoplanetary disk that couples an analytic description of the disk's temperature and density profile, chemical evolution, and an evolving dust distribution. We compute evolving radial distributions for a range of dust grain sizes, which depend on coagulation, fragmentation and radial drift processes. In particu…
▽ More
We have developed a new model for the astrochemical structure of a viscously evolving protoplanetary disk that couples an analytic description of the disk's temperature and density profile, chemical evolution, and an evolving dust distribution. We compute evolving radial distributions for a range of dust grain sizes, which depend on coagulation, fragmentation and radial drift processes. In particular we find that the water ice line plays an important role in shaping the radial distribution of the maximum grain size because ice coated grains are significantly less susceptible to fragmentation than their dry counterparts. This in turn has important effects on disk ionization and therefore on the location of dead zones. In comparison to a simple constant gas-to-dust ratio model for the dust as an example, we find that the new model predicts an outer dead zone edge that moves in by a factor of about 3 at 1 Myr (to 5 AU) and by a factor of about 14 by 3 Myr (to 0.5 AU). We show that the changing position of the dead zone and heat transition traps have important implications for the formation and trapping of planets in protoplanetary disks. Finally, we consider our results in light of recent ALMA observations of HL Tau and TW Hya.
△ Less
Submitted 8 November, 2016;
originally announced November 2016.
-
Ultra Violet Escape Fractions from Giant Molecular Clouds During Early Cluster Formation
Authors:
Corey S. Howard,
Ralph E. Pudritz,
Ralf S. Klessen
Abstract:
The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the Interstellar Medium (ISM), and extragalactic processes, such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction (f$_{esc}$) arising as a consequence of star cluster formation in a turbulent, 10$^6$ M$_{\odot}$ GMC, simulated using the…
▽ More
The UV photon escape fraction from molecular clouds is a key parameter for understanding the ionization of the Interstellar Medium (ISM), and extragalactic processes, such as cosmic reionization. We present the ionizing photon flux and the corresponding photon escape fraction (f$_{esc}$) arising as a consequence of star cluster formation in a turbulent, 10$^6$ M$_{\odot}$ GMC, simulated using the code FLASH. We make use of sink particles to represent young, star-forming clusters coupled with a radiative transfer scheme to calculate the emergent UV flux. We find that the ionizing photon flux across the cloud boundary is highly variable in time and space due to the turbulent nature of the intervening gas. The escaping photon fraction remains at $\sim$5% for the first 2.5 Myr, followed by two pronounced peaks at 3.25 and 3.8 Myr with a maximum f$_{esc}$ of 30% and 37%, respectively. These peaks are due to the formation of large HII regions, that expand into regions of lower density and some of which reach the cloud surface. However, these phases are short lived and f$_{esc}$ drops sharply as the HII regions are quenched by the central cluster passing through high-density material due to the turbulent nature of the cloud. We find an average f$_{esc}$ of 15% with factor of two variations over 1 Myr timescales. Our results suggest that assuming a single value for f$_{esc}$ from a molecular cloud is in general a poor approximation, and that the dynamical evolution of the system leads to large temporal variation.
△ Less
Submitted 10 January, 2017; v1 submitted 8 November, 2016;
originally announced November 2016.
-
On the Formation and Chemical Composition of Super Earths
Authors:
Matthew Alessi,
Ralph E. Pudritz,
Alex J. Cridland
Abstract:
Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they a…
▽ More
Super Earths are the largest population of exoplanets and are seen to exhibit a rich diversity of compositions as inferred through their mean densities. Here we present a model that combines equilibrium chemistry in evolving disks with core accretion that tracks materials accreted onto planets during their formation. In doing so, we aim to explain why super Earths form so frequently and how they acquire such a diverse range of compositions. A key feature of our model is disk inhomogeneities, or planet traps, that act as barriers to rapid type-I migration. The traps we include are the dead zone, which can be caused by either cosmic ray or X-ray ionization, the ice line, and the heat transition. We find that in disks with sufficiently long lifetimes ($\gtrsim$ 4 Myr), all traps produce Jovian planets. In these disks, planet formation in the heat transition and X-ray dead zone produces hot Jupiters while the ice line and cosmic ray dead zones produce Jupiters at roughly 1 AU. Super Earth formation takes place within short-lived disks ($\lesssim$ 2 Myr), whereby the disks are photoevaporated while planets are in a slow phase of gas accretion. We find that super Earth compositions range from dry and rocky ($<$ 6 % ice by mass) to those with substantial water contents ($>$ 30 % ice by mass). The traps play a crucial role in our results, as they dictate where in the disk particular planets can accrete from, and what compositions they are able to acquire.
△ Less
Submitted 13 September, 2016; v1 submitted 29 June, 2016;
originally announced June 2016.
-
Simulating radiative feedback and star cluster formation in GMCs: I. Dependence on gravitational boundedness
Authors:
Corey Howard,
Ralph Pudritz,
William Harris
Abstract:
Radiative feedback is an important consequence of cluster formation in Giant Molecular Clouds (GMCs) in which newly formed clusters heat and ionize their surrounding gas. The process of cluster formation, and the role of radiative feedback, has not been fully explored in different GMC environments. We present a suite of simulations which explore how the initial gravitational boundedness, and radia…
▽ More
Radiative feedback is an important consequence of cluster formation in Giant Molecular Clouds (GMCs) in which newly formed clusters heat and ionize their surrounding gas. The process of cluster formation, and the role of radiative feedback, has not been fully explored in different GMC environments. We present a suite of simulations which explore how the initial gravitational boundedness, and radiative feedback, affect cluster formation. We model the early evolution ($<$ 5 Myr) of turbulent, 10$^6$ M$_{\odot}$ clouds with virial parameters ranging from 0.5 to 5. To model cluster formation, we use cluster sink particles, coupled to a raytracing scheme, and a custom subgrid model which populates a cluster via sampling an IMF with an efficiency of 20\% per freefall time. We find that radiative feedback only decreases the cluster particle formation efficiency by a few percent. The initial virial parameter plays a much stronger role in limiting cluster formation, with a spread of cluster formation efficiencies of 37\% to 71\% for the most unbound to the most bound model. The total number of clusters increases while the maximum mass cluster decreases with an increasing initial virial parameter, resulting in steeper mass distributions. The star formation rates in our cluster particles are initially consistent with observations but rise to higher values at late times. This suggests that radiative feedback alone is not responsible for dispersing a GMC over the first 5 Myr of cluster formation.
△ Less
Submitted 16 June, 2016; v1 submitted 31 May, 2016;
originally announced June 2016.
-
Composition of Early Planetary Atmospheres I: Connecting Disk Astrochemistry to the Formation of Planetary Atmospheres
Authors:
Alex J. Cridland,
Ralph E. Pudritz,
Matthew Alessi
Abstract:
We present a model of the early chemical composition and elemental abundances of planetary atmospheres based on the cumulative gaseous chemical species that are accreted onto planets forming by core accretion from evolving protoplanetary disks. The astrochemistry of the host disk is computed using an ionization driven, non-equilibrium chemistry network within viscously evolving disk models. We acc…
▽ More
We present a model of the early chemical composition and elemental abundances of planetary atmospheres based on the cumulative gaseous chemical species that are accreted onto planets forming by core accretion from evolving protoplanetary disks. The astrochemistry of the host disk is computed using an ionization driven, non-equilibrium chemistry network within viscously evolving disk models. We accrete gas giant planets whose orbital evolution is controlled by planet traps using the standard core accretion model and track the chemical composition of the material that is accreted onto the protoplanet. We choose a fiducial disk model and evolve planets in 3 traps - water ice line, dead zone and heat transition. For a disk with a lifetime of 4.1 Myr we produce two Hot Jupiters (M = 1.43, 2.67 M$_{\rm Jupiter}$, r = 0.15, 0.11 AU) in the heat transition and ice line trap and one failed core (M = 0.003 M$_{\rm Jupiter}$, r =3.7 AU) in the dead zone. These planets are found with mixing ratios for CO and H$_2$O of $1.99\times 10^{-4}$, $5.0\times 10^{-4}$ respectively for both Hot Jupiters. Additionally for these planets we find CO$_2$ and CH$_4$, with mixing ratios of $1.8\times 10^{-6}\rightarrow 9.8\times 10^{-10}$ and $1.1\times 10^{-8}\rightarrow 2.3\times 10^{-10}$ respectively. These ranges correspond well with the mixing ratio ranges that have been inferred through the detection of emission spectra from Hot Jupiters by multiple authors. We compute a carbon-to-oxygen ratio of 0.227 for the ice line planet and 0.279 for the heat transition planet. These planets accreted their gas inside the ice line, hence the sub-solar C/O.
△ Less
Submitted 21 June, 2016; v1 submitted 30 May, 2016;
originally announced May 2016.
-
Filamentary flow and magnetic geometry in evolving cluster-forming molecular cloud clumps
Authors:
Mikhail Klassen,
Ralph E. Pudritz,
Helen Kirk
Abstract:
We present an analysis of the relationship between the orientation of magnetic fields and filaments that form in 3D magnetohydrodynamic simulations of cluster-forming, turbulent molecular cloud clumps. We examine simulated cloud clumps with size scales of L ~ 2-4 pc and densities of n ~ 400-1000 cm^-3 with Alfven Mach numbers near unity. We simulated two cloud clumps of different masses, one in vi…
▽ More
We present an analysis of the relationship between the orientation of magnetic fields and filaments that form in 3D magnetohydrodynamic simulations of cluster-forming, turbulent molecular cloud clumps. We examine simulated cloud clumps with size scales of L ~ 2-4 pc and densities of n ~ 400-1000 cm^-3 with Alfven Mach numbers near unity. We simulated two cloud clumps of different masses, one in virial equilibrium, the other strongly gravitationally bound, but with the same initial turbulent velocity field and similar mass-to-flux ratio. We apply various techniques to analyze the filamentary and magnetic structure of the resulting cloud, including the DisPerSE filament-finding algorithm in 3D. The largest structure that forms is a 1-2 parsec-long filament, with smaller connecting sub-filaments. We find that our simulated clouds, wherein magnetic forces and turbulence are comparable, coherent orientation of the magnetic field depends on the virial parameter. Subvirial clumps undergo strong gravitational collapse and magnetic field lines are dragged with the accretion flow. We see evidence of filament-aligned flow and accretion flow onto the filament in the subvirial cloud. Magnetic fields oriented more parallel in the subvirial cloud and more perpendicular in the denser, marginally bound cloud. Radiative feedback from a 16 Msun star forming in a cluster in one of our simulations ultimately results in the destruction of the main filament, the formation of an HII region, and the sweeping up of magnetic fields within an expanding shell at the edges of the HII region.
△ Less
Submitted 7 November, 2016; v1 submitted 27 May, 2016;
originally announced May 2016.