Source author record

Ralph E. Pudritz

Ralph E. Pudritz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph.GA astro-ph physics.chem-ph physics.ao-ph physics.comp-ph astro-ph.CO Biomolecules physics.geo-ph

Catalog footprint

What is connected

50works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Combined Effects of Disk Winds and Turbulence-Driven Accretion on Planet Populations

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.

preprint2022arXiv

Meteorites and the RNA World: Synthesis of Nucleobases in Carbonaceous Planetesimals and the Role of Initial Volatile Content

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.

preprint2022arXiv

Possible Ribose Synthesis in Carbonaceous Planetesimals

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.

preprint2022arXiv

Towards RNA life on Early Earth: From atmospheric HCN to biomolecule production in warm little ponds

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.

preprint2022arXiv

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Formation of Planetary Populations II: Effects of Initial Disk Size & Radial Dust Drift

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.

preprint2020arXiv

Formation of Planetary Populations III: Core Composition & Atmospheric Evaporation

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.

preprint2020arXiv

HCN production in Titan's Atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling

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.

preprint2018arXiv

Meteoritic Abundances of Fatty Acids and Potential Reaction Pathways in Planetesimals

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.

preprint2017arXiv

Composition of Early Planetary Atmospheres II: Coupled Dust and Chemical Evolution in Protoplanetary Disks

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.

preprint2016arXiv

Filamentary flow and magnetic geometry in evolving cluster-forming molecular cloud clumps

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.

preprint2016arXiv

Meteorites and the RNA World: A Thermodynamic Model of Nucleobase Synthesis within Planetesimals

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.

preprint2016arXiv

On the Formation and Chemical Composition of Super Earths

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.

preprint2016arXiv

Radial Drift of Dust in Protoplanetary Disks: The Evolution of Ice lines and Dead zones

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.

preprint2016arXiv

The Search for Extraterrestrial Intelligence in Earth's Solar Transit Zone

Over the past few years, astronomers have detected thousands of planets and candidate planets by observing their periodic transits in front of their host stars. Related methods might soon allow studies of the chemical imprints of life in extrasolar planetary atmospheres. Here, we address the reciprocal question, namely, from where is Earth detectable by extrasolar observers using similar methods. We explore Earth's transit zone (ETZ), the projection of a band around Earth's ecliptic onto the celestial plane, where observers can detect Earth transits across the Sun. The ETZ is between $0.520^\circ$ and $0.537^\circ$ wide due to the non-circular Earth orbit. The restricted ETZ (rETZ), where Earth transits the Sun less than 0.5 solar radii from its center, is ~$0.262^\circ$ wide. We first compile a target list of 45 K and 37 G dwarf stars inside the rETZ and within 1 kiloparsec (3260 lightyears) using the Hipparcos catalog. We then greatly enlarge the number of potential targets by constructing an analytic galactic disk model and find that ~$10^5$ K and G dwarf stars should reside within the rETZ. The ongoing GAIA space mission can potentially discover all G dwarfs among them (several $10^4$) within the next five years. Many more potentially habitable planets orbit dim, unknown M stars in the ETZ and other stars that traversed the ETZ long time ago. If any of these planets host intelligent observers, they could have identified Earth as a habitable, or even as a living, world, and we could be receiving their broadcasts today. The K2 mission, the Allen Telescope Array, the upcoming Square Kilometer Array, or the Green Bank Telescope might detect such deliberate extraterrestrial messages. Ultimately, the ETZ would be an ideal region to monitor by the Breakthrough Listen Initiatives, an upcoming survey that will constitute the most comprehensive search for extraterrestrial intelligence so far.

preprint2015arXiv

Conditions for water ice lines and Mars-mass exomoons around accreting super-Jovian planets at 1 - 20 AU from Sun-like stars

Exomoon detections might be feasible with NASA's Kepler or ESA's upcoming PLATO mission or the ground-based E-ELT. To use observational resources most efficiently we need to know where the largest, most easily detected moons can form. We explore the possibility of large exomoons by following the movement of water (H2O) ice lines in the accretion disks around young super-Jovian planets. We want to know how different heating sources in those disks affect the H2O ice lines. We simulate 2D rotationally symmetric accretion disks in hydrostatic equilibrium around super-Jovian exoplanets. The energy terms in our semi-analytical model -- (1) viscous heating, (2) planetary illumination, (3) accretional heating, and (4) stellar illumination -- are fed by precomputed planet evolution tracks. We consider planets accreting 1 to 12 Jupiter masses at distances between 1 and 20 AU to a Sun-like star. Accretion disks around Jupiter-mass planets closer than ~4.5 AU to Sun-like stars do not feature H2O ice lines, but the most massive super-Jovians can form icy satellites as close as ~3 AU to Sun-like stars. Super-Jovian planets forming beyond ~5 AU can host Mars-mass moons. We study a broad range of disk parameters for planets at 5.2 AU and find that the H2O ice lines are universally between ~15 and 30 Jupiter radii when the last generation of moons is forming. If the abundant population of super-Jovian planets at ~1 AU formed in situ, then they should lack giant icy moons because their disks did not host H2O ice in the final stages of accretion. In the more likely case that these planets migrated to their current locations from beyond a few AU, they might be orbited by large, H2O-rich moons. In this case, Mars-mass ocean moons might be common in the stellar habitable zones. Future exomoon searches can provide powerful constraints on the formation and migration history of giant exoplanets.

preprint2015arXiv

Nature's Starships II: Simulating the Synthesis of Amino Acids in Meteorite Parent Bodies

Carbonaceous chondrite meteorites are known for having high water and organic material contents, including amino acids. Here we address the origin of amino acids in the warm interiors of their parent bodies (planetesimals) within a few million years of their formation, and connect this with the astrochemistry of their natal protostellar disks. We compute both the total amino acid abundance pattern as well as the relative frequencies of amino acids within the CM2 (e.g. Murchison) and CR2 chondrite subclasses based on Strecker reactions within these bodies. We match the relative frequencies to well within an order of magnitude among both CM2 and CR2 meteorites for parent body temperatures $<$ 200$^{\circ}$C. These temperatures agree with 3D models of young planetesimal interiors. We find theoretical abundances of approximately 7x10$^5$ parts-per-billion (ppb), which is in agreement with the average observed abundance in CR2 meteorites of 4$\pm$7x10$^5$, but an order of magnitude higher than the average observed abundance in CM2 meteorites of 2$\pm$2x10$^4$. We find that the production of hydroxy acids could be favoured over the production of amino acids within certain meteorite parent bodies (e.g. CI1, CM2) but not others (e.g. CR2). This could be due to the relatively lower NH$_3$ abundances within CI1 and CM2 meteorite parent bodies, which leads to less amino acid synthesis. We also find that the water content in planetesimals is likely to be the main cause of variance between carbonaceous chondrites of the same subclass. We propose that amino acid abundances are primarily dependent on the ammonia and water content of planetesimals that are formed in chemically distinct regions within their natal protostellar disks.

preprint2015arXiv

Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations

We present an algorithm for solving the radiative transfer problem on massively parallel computers using adaptive mesh refinement and domain decomposition. The solver is based on the method of characteristics which requires an adaptive raytracer that integrates the equation of radiative transfer. The radiation field is split into local and global components which are handled separately to overcome the non-locality problem. The solver is implemented in the framework of the magneto-hydrodynamics code FLASH and is coupled by an operator splitting step. The goal is the study of radiation in the context of star formation simulations with a focus on early disc formation and evolution. This requires a proper treatment of radiation physics that covers both the optically thin as well as the optically thick regimes and the transition region in particular. We successfully show the accuracy and feasibility of our method in a series of standard radiative transfer problems and two 3D collapse simulations resembling the early stages of protostar and disc formation.

preprint2015arXiv

Seeding the Pregenetic Earth: Meteoritic Abundances of Nucleobases and Potential Reaction Pathways

Carbonaceous chondrites are a class of meteorite known for having a high content of water and organics. In this study, abundances of the nucleobases, i.e., the building blocks of RNA and DNA, found in carbonaceous chondrites are collated from a variety of published data and compared across various meteorite classes. An extensive review of abiotic chemical reactions producing nucleobases is then performed. These reactions are then reduced to a list of 15 individual reaction pathways that could potentially occur within meteorite parent bodies. The nucleobases guanine, adenine and uracil are found in carbonaceous chondrites in the amounts of 1$-$500 ppb. It is currently unknown which reaction is responsible for their synthesis within the meteorite parent bodies. One class of carbonaceous meteorites dominate the abundances of both amino acids and nucleobases$-$the so-called CM2 (e.g. Murchison meteorite). CR2 meteorites (e.g. Graves Nunataks) also dominate the abundances of amino acids, but are the least abundant in nucleobases. The abundances of total nucleobases in these two classes are $330 \pm250$ ppb and $16 \pm13$ ppb respectively. Guanine most often has the greatest abundances in carbonaceous chondrites with respect to the other nucleobases, but is 1$-$2 orders of magnitude less abundant in CM2 meteorites than glycine (the most abundant amino acid). Our survey of the reaction mechanisms for nucleobase formation suggests that Fischer-Tropsch synthesis (i.e. CO, H$_2$ and NH$_3$ gases reacting in the presence of a catalyst such as alumina or silica) is the most likely for conditions that characterize early states of planetesimals.

preprint2014arXiv

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodynamics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion (FLD) solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

preprint2014arXiv

Dust in Protoplanetary Disks: A Clue as to the Critical Mass of Planetary Cores

Dust in protoplanetary disks is recognized as the building blocks of planets. In the core accretion scenario, the abundance of dust in disks (or metallicity) is crucial for forming cores of gas giants. We present our recent progress on the relationship between the metallicity and planet formation, wherein planet formation frequencies (PFFs) as well as the critical mass of planetary cores ($M_{c,crit}$) that can initiate gas accretion are statistically examined. We focus on three different planetary populations that are prominent for observed exoplanets in the mass-semimajor axis diagram: hot Jupiters, exo-Jupiters that are densely populated around 1 AU, and low-mass planets in tight orbits. We show that the resultant PFFs for both Jovian planets are correlated positively with the metallicity whereas low-mass planets form efficiently for a wide range of metallicities. This is consistent with the observed Planet-Metallicity correlation. Examining the statistically averaged value of $M_{c,crit}$ (defined as $<M_{c,crit}>$), we find that the correlation originates from the behavior of $<M_{c,crit}>$ that increases steadily with metallicity for two kinds of the Jovian planets while the low-mass planets obtain a rather constant value for $<M_{c,crit}>$. Such a difference in $<M_{c,crit}>$ can define transition metallicities (TMs) above which the Jovian planets gain a larger value of $<M_{c,crit}>$ than the low-mass planets, and hence gas giant formation takes place more efficiently. We find that TMs are sensitive to the important parameter that involves $M_{c,crit}$. We show, by comparing with the observations, that a most likely value of $M_{c,crit}$ is $\simeq 5M_{\oplus}$, which is smaller than the conventional value in the literature ($\simeq 10M_{\oplus}$). Our results suggest that opacities in planetary atmospheres play an important role for lowering $M_{c,crit}$.

preprint2014arXiv

Nature's Starships. I. Observed Abundances and Relative Frequencies of Amino Acids in Meteorites

The class of meteorites called carbonaceous chondrites are examples of material from the solar system which have been relatively unchanged from the time of their initial formation. These meteorites have been classified according to the temperatures and physical conditions of their parent planetesimals. We collate available data on amino acid abundance in these meteorites and plot the concentrations of different amino acids for each meteorite within various meteorite subclasses. We plot average concentrations for various amino acids across meteorites separated by subclass and petrologic type. We see a predominance in the abundance and variety of amino acids in CM2 and CR2 meteorites. The range in temperature corresponding to these subclasses indicates high degrees of aqueous alteration, suggesting aqueous synthesis of amino acids. Within the CM2 and CR2 subclasses, we identify trends in relative frequencies of amino acids to investigate how common amino acids are as a function of their chemical complexity. These two trends (total abundance and relative frequencies) can be used to constrain formation parameters of amino acids within planetesimals. Our organization of the data supports an onion shell model for the temperature structure of planetesimals. The least altered meteorites (type 3) and their amino acids originated near cooler surface regions. The most active amino acid synthesis likely took place at intermediate depths (type 2). The most altered materials (type 1) originated furthest toward parent body cores. This region is likely too hot to either favor amino acid synthesis or for amino acids to be retained after synthesis.

preprint2014arXiv

Planet Traps and Planetary Cores: Origins of the Planet-Metallicity Correlation

Massive exoplanets are observed preferentially around high metallicity ([Fe/H]) stars while low-mass exoplanets do not show such an effect. This so-called planet-metallicity correlation generally favors the idea that most observed gas giants at $r<10$ AU are formed via a core accretion process. We investigate the origin of this phenomenon using a semi-analystical model, wherein the standard core accretion takes place at planet traps in protostellar disks where rapid type I migrators are halted. We focus on the three major exoplanetary populations - hot-Jupiters, exo-Jupiters located at $r \simeq 1$ AU, and the low-mass planets. We show using a statistical approach that the planet-metallicity correlations are well reproduced in these models. We find that there are specific transition metallicities with values [Fe/H]$=-0.2$ to $-0.4$, below which the low-mass population dominates, and above which the Jovian populations take over. The exo-Jupiters significantly exceed the hot-Jupiter population at all observed metallicities. The low-mass planets formed via the core accretion are insensitive to metallicity, which may account for a large fraction of the observed super-Earths and hot-Neptunes. Finally, a controlling factor in building massive planets is the critical mass of planetary cores ($M_{c,crit}$) that regulates the onset of runaway gas accretion. Assuming the current data is roughly complete at [Fe/H]$>-0.6$, our models predict that the most likely value of the "mean" critical core mass of Jovian planets is $\braket{M_{c,crit}} \simeq 5 M_{\oplus}$ rather than $10 M_{\oplus}$. This implies that grain opacities in accreting envelopes should play an important role in lowering $M_{c,crit}$.

preprint2014arXiv

The Earliest Stages of Star and Planet Formation: Core Collapse, and the Formation of Disks and Outflows

(Abridged) In this review we focus on the observations and theory of the formation of early disks and outflows, and their connections with the first phases of planet formation. Large rotationally supported circumstellar disks, although common around more evolved young stellar objects, are rarely detected during the earliest, "Class 0" phase; however, a few excellent candidates have been discovered recently around both low and high mass protostars. In this early phase, prominent outflows are ubiquitously observed; they are expected to be associated with at least small magnetized disks. Disk formation - once thought to be a simple consequence of the conservation of angular momentum during hydrodynamic core collapse - is far more subtle in magnetized gas. In this case, the rotation can be strongly magnetically braked. Indeed, both analytic arguments and numerical simulations have shown that disk formation is suppressed in the strict ideal magnetohydrodynamic (MHD) limit for the observed level of core magnetization. We review what is known about this "magnetic braking catastrophe", possible ways to resolve it, and the current status of early disk observations. Outflows are also intimately linked to disk formation; they are a natural product of magnetic fields and rotation and are important signposts of star formation. We review new developments on early outflow generation since PPV. The properties of early disks and outflows are a key component of planet formation in its early stages and we review these major connections.

preprint2013arXiv

Cluster formation in molecular clouds: I. stellar populations, star formation rates, and ionizing radiation

We present a model for the radiative output of star clusters in the process of star formation suitable for use in hydrodynamical simulations of radiative feedback. Gas in a clump, defined as a region whose density exceeds 10^4 cm^-3, is converted to stars via the random sampling of the Chabrier IMF. A star formation efficiency controls the rate of star formation. We have completed a suite of simulations which follow the evolution of accretion-fed clumps with initial masses ranging from 0 to 10^5 M_sol and accretion rates ranging from 10^-5 to 10^-1 M_sol yr^-1. The stellar content is tracked over time which allows the aggregate luminosity, ionizing photon rate, number of stars, and star formation rate (SFR) to be determined. For a fiducial clump of 10^4 M_sol, the luminosity is ~4x10^6 L_sol with a SFR of roughly 3x10^-3 M_sol yr^-1. We identify two regimes in our model. The accretion-dominated regime obtains the majority of its gas through accretion and is characterized by an increasing SFR while the reservoir-dominated regime has the majority of its mass present in the initial clump with a decreasing SFR. We show that our model can reproduce the expected number of O stars, which dominate the radiative output of the cluster. We find a nearly linear relationship between SFR and mass as seen in observations. We conclude that our model is an accurate and straightforward way to represent the output of clusters in hydrodynamical simulations with radiative feedback.

preprint2013arXiv

Disk Inhomogeneities and the Origins of Planetary System Architectures and Observational Properties

Recent high-resolution observations show that protoplanetary disks have various kinds of structural properties or inhomogeneities. These are the consequence of a mixture of a number of physical and chemical processes taking place in the disks. Here, we discuss the results of our comprehensive investigations on how disk inhomogeneities affect planetary migration. We demonstrate that disk inhomogeneities give rise to planet traps - specific sites in protoplanetary disks at which rapid type I migration is halted. We show that up to three types of traps (heat transitions, ice lines and dead zones) can exist in a single disk, and that they move differently as the disk accretion rate decreases with time. We also demonstrate that the position of planet traps strongly depends on stellar masses and disk accretion rates. This indicates that host stars establish preferred (initial) scales of their planetary systems. Finally, we discuss the possible observational signatures of disk inhomogeneities.

preprint2013arXiv

From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation

Recent studies of the nearest star-forming clouds of the Galaxy at submillimeter wavelengths with the Herschel Space Observatory have provided us with unprecedented images of the initial and boundary conditions of the star formation process. The Herschel results emphasize the role of interstellar filaments in the star formation process and connect remarkably well with nearly a decade's worth of numerical simulations and theory that have consistently shown that the ISM should be highly filamentary on all scales and star formation is intimately related to self-gravitating filaments. In this review, we trace how the apparent complexity of cloud structure and star formation is governed by relatively simple universal processes - from filamentary clumps to galactic scales. We emphasize two crucial and complementary aspects: (i) the key observational results obtained with Herschel over the past three years, along with relevant new results obtained from the ground on the kinematics of interstellar structures, and (ii) the key existing theoretical models and the many numerical simulations of interstellar cloud structure and star formation. We then synthesize a comprehensive physical picture that arises from the confrontation of these observations and simulations.

preprint2013arXiv

Giant Molecular Cloud Formation in Disk Galaxies: Characterizing Simulated versus Observed Cloud Catalogues

We present the results of a study of simulated Giant Molecular Clouds (GMCs) formed in a Milky Way-type galactic disk with a flat rotation curve. This simulation, which does not include star formation or feedback, produces clouds with masses ranging between 10^4 Msun and 10^7 Msun. We compare our simulated cloud population to two observational surveys; The Boston University- Five College Radio Astronomy Observatory Galactic Ring Survey and the BIMA All-Disk Survey of M33. An analysis of the global cloud properties as well as a comparison of Larson's scaling relations is carried out. We find that simulated cloud properties agree well with the observed cloud properties, with the closest agreement occurring between the clouds at comparable resolution in M33. Our clouds are highly filamentary - a property that derives both from their formation due to gravitational instability in the sheared galactic environment, as well as to cloud- cloud gravitational encounters. We also find that the rate at which potentially star forming gas accumulates within dense regions - wherein n_{thresh} > 10^4 cm^{-3} - is 3% per 10 Myr, in clouds of roughly 10^6 Msun. This suggests that star formation rates in observed clouds are related to the rates at which gas can be accumulated into dense subregions within GMCs via filamentary flows. The most internally well-resolved clouds are chosen for listing in a catalogue of simulated GMCs; the first of its kind. The catalogued clouds are available as an extracted data set from the global simulation.

preprint2013arXiv

Planetary Populations in the Mass-Period Diagram: A Statistical Treatment of Exoplanet Formation and the Role of Planet Traps

The rapid growth in the number of known exoplanets has revealed the existence of several distinct planetary populations in the observed mass-period diagram. Two of the most surprising are, (1) the concentration of gas giants around 1AU and (2) the accumulation of a large number of low-mass planets with tight orbits, also known as super-Earths and hot Neptunes. We have recently shown that protoplanetary disks have multiple planet traps that are characterized by orbital radii in the disks and halt rapid type I planetary migration. By coupling planet traps with the standard core accretion scenario, we showed that one can account for the positions of planets in the mass-period diagram. In this paper, we demonstrate quantitatively that most gas giants formed at planet traps tend to end up around 1 AU with most of these being contributed by dead zones and ice lines. In addition, we show that a large fraction of super-Earths and hot Neptunes are formed as "failed" cores of gas giants - this population being constituted by comparable contributions from dead zone and heat transition traps. Our results are based on the evolution of forming planets in an ensemble of disks where we vary only the lifetimes of disks as well as their mass accretion rates onto the host star. We show that a statistical treatment of the evolution of a large population of planetary cores initially caught in planet traps accounts for the existence of three distinct exoplantary populations - the hot Jupiters, the more massive planets at roughly orbital radii around 1 AU orbital, and the short period SuperEarths and hot Neptunes. There are very few evolutionary tracks that feed into the large orbital radii characteristic of the imaged Jovian planet, which agrees with recent surveys. Finally, we find that low-mass planets in tight orbits become the dominant planetary population for low mass stars ($M_* \le 0.7 M_{\odot}$).

preprint2012arXiv

Evolutionary Tracks of Trapped, Accreting Protoplanets: the Origin of the Observed Mass-Period Relation

The large number of observed exoplanets ($\gtrsim $ 700) provides important constraints on their origin as deduced from the mass-period diagram of planets. The most surprising features in the diagram are 1) the (apparent) pile up of gas giants at a period of $\sim 500$ days ($\sim1$ AU) and 2) the so-called mass-period relation which indicates that planetary mass is an increasing function of orbital period. We construct the evolutionary tracks of growing planets at planet traps in evolving protoplanetary disks and show that they provide a good physical understanding of how these observational properties arise. The fundamental feature of our model is that inhomogeneities in protoplanetary disks give rise to multiple (up to 3) trapping sites for rapid (type I) planetary migration of planetary cores. The viscous evolution of disks results in the slow radial movement of the traps and their cores from large to small orbital periods. In our model, the slow inward motion of planet traps is coupled with the standard core accretion scenario for planetary growth. As planets grow, type II migration takes over. Planet growth and radial movement are ultimately stalled by the dispersal of gas disks via photoevaporation. Our model makes a number of important predictions: that distinct sub-populations of planets that reflect the properties of planet traps where they have grown result in the mass-period relation; that the presence of these sub-populations naturally explains a pile-up of planets at $\sim 1$ AU; and that evolutionary tracks from the ice line do put planets at short periods and fill an earlier claimed "planet desert" - sparse population of planets in the mass-semi-major axis diagram.

preprint2012arXiv

HII region variability and pre-main-sequence evolution

Recent observations and simulations have suggested that HII regions around massive stars may vary in their size and emitted flux on timescales short enough to be observed. This variability can have a number of causes, ranging from environmental causes to variability of the ionizing source itself. We explore the latter possibility by considering the pre-main-sequence evolution of massive protostars and conducting numerical simulations with ionizing radiation feedback using the FLASH AMR hydrodynamics code. We investigate three different models: a simple ZAMS model, a self-consistent one-zone model by Offner et al. (2009), and a model fit to the tracks computed by Hosokawa & Omukai (2009). The protostellar models show that hypercompact HII regions around massive, isolated protostars collapse or shrink from diameters of 80 or 300 AU, depending on the model choice, down to near absence during the swelling of stellar radius that accompanies the protostar's transition from a convective to a radiative internal structure. This occurs on timescales as short as ~3000 years.

preprint2012arXiv

Shock interactions, turbulence, and the origin of the stellar mass spectrum

Supersonic turbulence is an essential element in understanding how structure within interstellar gas is created and shaped. In the context of star formation, many computational studies show that the mass spectrum of density and velocity fluctuations within dense clouds, as well as the distribution of their angular momenta, trace their origin to the statistical and physical properties of gas that is lashed with shock waves. In this article, we review the observations, simulations, and theories of how turbulent-like processes can account for structures we see in molecular clouds. We then compare traditional ideas of supersonic turbulence with a simpler physical model involving the effects of multiple shock waves and their interaction in the interstellar medium. Planar intersecting shock waves produce dense filaments, and generate vortex sheets that are essential to create the broad range of density and velocity structure in clouds. As an example, the lower mass behaviour of the stellar initial mass function can be traced to the tendency of a collection of shock waves to build-up a log-normal density distribution (or column density). Vorticity - which is essential to produce velocity structure over a very broad range of length scales in shocked clouds - can also be generated by the passage of curved shocks or intersecting planar shocks through such media. Two major additional physical forces affect the structure of star forming gas - gravity and feedback processes from young stars. Both of these can produce power-law tails at the high mass end of the IMF.

preprint2012arXiv

Simulating protostellar evolution and radiative feedback in the cluster environment

Radiative feedback is among the most important consequences of clustered star formation inside molecular clouds. At the onset of star formation, radiation from massive stars heats the surrounding gas, which suppresses the formation of many low-mass stars. When simulating pre-main-sequence stars, their stellar properties must be defined by a prestellar model. Different approaches to prestellar modeling may yield quantitatively different results. In this paper, we compare two existing prestellar models under identical initial conditions to gauge whether the choice of model has any significant effects on the final population of stars. The first model treats stellar radii and luminosities with a ZAMS model, while separately estimating the accretion luminosity by interpolating to published prestellar tracks. The second, more accurate prestellar model self-consistently evolves the radius and luminosity of each star under highly variable accretion conditions. Each is coupled to a raytracing-based radiative feedback code that also treats ionization. The impact of the self-consistent model is less ionizing radiation and less heating during the early stages of star formation. This may affect final mass distributions. We noted a peak stellar mass reduced by 8% from 47.3 Msun to 43.5 Msun in the evolutionary model, relative to the track-fit model. Also, the difference in mass between the two largest stars in each case is reduced from 14 Msun to 7.5 Msun. The HII regions produced by these massive stars were also seen to flicker on timescales down to the limit imposed by our timestep (< 560 years), rapidly changing in size and shape, confirming previous cluster simulations using ZAMS-based estimates for prestellar ionizing flux.

preprint2012arXiv

The long-term evolution of warped, magnetised discs, and precessing outflows in collapsing pre-stellar cores

(abridged) The nature of early Class 0/I protostellar discs is not clearly understood. Early protostellar discs are needed to drive molecular outflows and jets observed in star forming regions, but there has been some debate to how they form. From a theoretical perspective, the consequences of disc and outflow generation are crucial to understanding the very nature of how stars are assembled. We have performed 3D ideal magnetohydrodynamic (MHD) simulations of collapsing Bonnor-Ebert spheres, employing sink particles with a radius of 3.2 AU alongside an AMR grid and using a cooling function to model radiative cooling of the gas. This has allowed us to explore 2-8x10^4 yr further into the evolution of an early Class 0 disc-outflow system than previous simulations. Our outflow is precessing, kinked, turbulent, contains internal shocks and has a scale of 0.1 pc end-to-end. We form a rotationally dominated disc with a radius of 100 AU embedded inside a transient, unstable, flattened, rotating core extending out to 2000 AU. The larger flattened structure launches a low speed wind (v_r < 1.5 km/s) dominated by B_phi, while the inner disc launches a centrifugally driven jet dominated by B_p with speeds up to 20$ km/s. From the inner disk, the value of dM/dt_out/dM/dt_in ~ 0.1, wheras in the outer core dM/dt_out/dM/dt_in ~ 1.0. The inner disc becomes unstable to a warping instability due to the magnetic structure of the outflow and warps to 30 deg with respect to the z-axis by the end of the simulation. The envelope is cleared out and is less massive than the disc. We measure star formation efficiencies of eta_core=0.63 (and growing), higher than theoretical predictions. This indicates that outflows are not as efficient at expelling envelope mass as some current models estimate. We discuss the relevance of our disc misalignment concerning the formation of mis-aligned hot Jupiters.

preprint2011arXiv

Dead Zones and the Diversity of Exoplanetary Systems

Planetary migration provides a theoretical basis for the observed diversity of exoplanetary systems. We demonstrate that dust settling - an inescapable feature of disk evolution - gives even more rapid type I migration by up to a factor of about 2 than occurs in disks with fully mixed dust. On the other hand, type II migration becomes slower by a factor of 2 due to dust settling. This even more problematic type I migration can be resolved by the presence of a dead zone; the inner, high density region of a disk which features a low level of turbulence. We show that enhanced dust settling in the dead zone leaves a dusty wall at its outer edge. Back-heating of the dead zone by this wall produces a positive radial gradient for the disk temperature, which acts as a barrier for type I migration.

preprint2011arXiv

Spin Evolution of Accreting Young Stars. II. Effect of Accretion-Powered Stellar Winds

We present a model for the rotational evolution of a young, solar-mass star interacting magnetically with an accretion disk. As in a previous paper (Paper I), the model includes changes in the star's mass and radius as it descends the Hayashi track, a decreasing accretion rate, and a prescription for the angular momentum transfer between the star and disk. Paper I concluded that, for the relatively strong magnetic coupling expected in real systems, additional processes are necessary to explain the existence of slowly rotating pre-main-sequence stars. In the present paper, we extend the stellar spin model to include the effect of a spin-down torque that arises from an accretion-powered stellar wind. For a range of magnetic field strengths, accretion rates, initial spin rates, and mass outflow rates, the modeled stars exhibit rotation periods within the range of 1--10 days in the age range of 1--3 Myr. This range coincides with the bulk of the observed rotation periods, with the slow rotators corresponding to stars with the lowest accretion rates, strongest magnetic fields, and/or highest stellar wind mass outflow rates. We also make a direct, quantitative comparison between the accretion-powered stellar wind scenario and the two types of disk-locking models (namely the X-wind and Ghosh & Lamb type models) and identify some remaining theoretical issues for understanding young star spins.

preprint2011arXiv

The Origin of Planetary System Architectures. I. Multiple Planet Traps in Gaseous Discs

The structure of planetary systems around their host stars depends on their initial formation conditions. Massive planets will likely be formed as a consequence of rapid migration of planetesimals and low mass cores into specific trapping sites in protoplanetary discs. We present analytical modeling of inhomogeneities in protoplanetary discs around a variety of young stars, - from Herbig Ae/Be to classical T Tauri and down to M stars, - and show how they give rise to planet traps. The positions of these traps define the initial orbital distribution of multiple protoplanets. We investigate both corotation and Lindblad torques, and show that a new trap arises from the (entropy-related) corotation torque. This arises at that disc radius where disc heating changes from viscous to stellar irradiation dominated processes. We demonstrate that up to three traps (heat transitions, ice lines and dead zones) can exist in a single disc, and that they move differently as the disc accretion rate $\dot{M}$ decreases with time. The interaction between the giant planets which grow in such traps may be a crucial ingredient for establishing planetary systems. We also demonstrate that the position of planet traps strongly depends on stellar masses and disc accretion rates. This indicates that host stars establish preferred scales of planetary systems formed around them. We discuss the potential of planet traps induced by ice lines of various molecules such as water and CO, and estimate the maximum and minimum mass of planets which undergo type I migration. We finally apply our analyses to accounting for the initial conditions proposed in the Nice model for the origin of our Solar system.

preprint2010arXiv

Dead Zones as Thermal Barriers to Rapid Planetary Migration in Protoplanetary Disks

Planetary migration in standard models of gaseous protoplanetary disks is known to be very rapid ($\sim 10^5$ years) jeopardizing the existence of planetary systems. We present a new mechanism for significantly slowing rapid planetary migration, discovered by means of radiative transfer calculations of the thermal structure of protoplanetary disks irradiated by their central stars. Rapid dust settling in a disk's dead zone - a region with very little turbulence - leaves a dusty wall at its outer edge. We show that the back-heating of the dead zone by this irradiated wall produces a positive gradient of the disk temperature which acts as a thermal barrier to planetary migration which persists for the disk lifetime. Although we analyze in detail the migration of a Super-Earth in a low mass disk around an M star, our findings can apply to wide variety of young planetary systems. We compare our findings with other potentially important stopping mechanisms and show that there are large parameter spaces for which dead zones are likely to play the most important role for reproducing the observed mass-period relation in longer planetary periods.

preprint2010arXiv

Discs, outflows, and feedback in collapsing magnetized cores

The pre-stellar cores in which low mass stars form are generally well magnetized. Our simulations show that early protostellar discs are massive and experience strong magnetic torques in the form of magnetic braking and protostellar outflows. Simulations of protostellar disk formation suggest that these torques are strong enough to suppress a rotationally supported structure from forming for near critical values of mass-to-flux. We demonstrate through the use of a 3D adaptive mesh refinement code -- including cooling, sink particles and magnetic fields -- that one produces transient 1000 AU discs while simultaneously generating large outflows which leave the core region, carrying away mass and angular momentum. Early inflow/outflow rates suggest that only a small fraction of the mass is lost in the initial magnetic tower/jet event.

preprint2010arXiv

Dust Settling and Rapid Planetary Migration

Planetary migration is essential to explain the observed mass-period relation for exoplanets. Without some stopping mechanism, the tidal, resonant interaction between planets and their gaseous disc generally causes the planets to migrate inward so efficiently that they plunge into the host star within the gaseous disc lifetime ($\sim $ 1-3 Myrs). We investigate planetary migration by analytically calculating the migration rate and time within self-consistently computed, radiatively heated discs around M stars in which the effects of dust settling are included. We show that dust settling lowers the disc temperature and raises the gas density in the mid-plane. This inescapable evolution of disc structure speeds up type I planetary migration for lower mass bodies by up to a factor of about 2. We also examine the effects of dust settling on the gap-opening mass and type II migration, and find that the gap-opening mass is reduced by a factor of 2 and type II migration becomes slower by a factor of 2. While dust settling can somewhat alleviate the problem of planetary migration for more massive planets, the more rapid migration of low mass planets and planetary cores requires a robust slowing mechanism.

preprint2010arXiv

Simulations of the IMF in Clusters

We review computational approaches to understanding the origin of the Initial Mass Function (IMF) during the formation of star clusters. We examine the role of turbulence, gravity and accretion, equations of state, and magnetic fields in producing the distribution of core masses - the Core Mass Function (CMF). Observations show that the CMF is similar in form to the IMF. We focus on feedback processes such as stellar dynamics, radiation, and outflows can reduce the accreted mass to give rise to the IMF. Numerical work suggests that filamentary accretion may play a key role in the origin of the IMF.

preprint2009arXiv

A thermodynamic basis for prebiotic amino acid synthesis and the nature of the first genetic code

Of the twenty amino acids used in proteins, ten were formed in Miller's atmospheric discharge experiments. The two other major proposed sources of prebiotic amino acid synthesis include formation in hydrothermal vents and delivery to Earth via meteorites. We combine observational and experimental data of amino acid frequencies formed by these diverse mechanisms and show that, regardless of the source, these ten early amino acids can be ranked in order of decreasing abundance in prebiotic contexts. This order can be predicted by thermodynamics. The relative abundances of the early amino acids were most likely reflected in the composition of the first proteins at the time the genetic code originated. The remaining amino acids were incorporated into proteins after pathways for their biochemical synthesis evolved. This is consistent with theories of the evolution of the genetic code by stepwise addition of new amino acids. These are hints that key aspects of early biochemistry may be universal.

preprint2009arXiv

Radiatively heated, protoplanetary discs with dead zones. I. Dust settling and thermal structure of discs around M stars

The irradiation of protoplanetary discs by central stars is the main heating mechanism for discs, resulting in their flared geometric structure. In a series of papers, we investigate the deep links between 2D self-consistent disc structure and planetary migration in irradiated discs, focusing particularly on those around M stars. In this first paper, we analyse the thermal structure of discs that are irradiated by an M star by solving the radiative transfer equation by means of a Monte Carlo code. Our simulations of irradiated hydrostatic discs are realistic and self-consistent in that they include dust settling with multiple grain sizes (N=15), the gravitational force of an embedded planet on the disc, and the presence of a dead zone (a region with very low levels of turbulence) within it. We show that dust settling drives the temperature of the mid-plane from an $r^{-3/5}$ distribution (well mixed dust models) toward an $r^{-3/4}$. The dead zone, meanwhile, leaves a dusty wall at its outer edge because dust settling in this region is enhanced compared to the active turbulent disc at larger disc radii. The disc heating produced by this irradiated wall provides a positive gradient region of the temperature in the dead zone in front of the wall. This is crucially important for slowing planetary migration because Lindblad torques are inversely proportional to the disc temperature. Furthermore, we show that low turbulence of the dead zone is self-consistently induced by dust settling, resulting in the Kelvin-Helmholtz instability (KHI). We show that the strength of turbulence arising from the KHI in the dead zone is $α=10^{-5}$.

preprint2009arXiv

Shock generated vorticity in the interstellar medium and the origin of the stellar initial mass function

The main observational evidence for turbulence in the interstellar medium (ISM) and molecular clouds is the power-law energy spectrum for velocity fluctuations, E(k) \propto k^α. The Kolmogorov scaling exponent, α=-5/3, is typical. At the same time, the observed probability distribution function (PDF) of gas densities in both the ISM as well as in molecular clouds is a log-normal distribution, which is similar to the initial mass function (IMF) that describes the distribution of stellar masses. We examine the density and velocity structure of interstellar gas traversed by curved shock waves in the kinematic limit. We demonstrate mathematically that just a few passages of curved shock waves generically produces a log-normal density PDF. This explains the ubiquity of the log-normal PDF in many different numerical simulations. We also show that subsequent interaction with a spherical blast wave generates a power-law density distribution at high densities, qualitatively similar to the Salpeter power-law for the IMF. Finally, we show that a focused shock produces a downstream flow with energy spectrum exponent α=-2. Subsequent shock passages reduce this slope, achieving α\approx -5/3 after a few passages. Subsequent dissipation of energy piled up at the small scales will act to maintain the spectrum very near to the Kolomogorov value. Therefore, fully-developed turbulence may not be required to explain the observed energy spectrum and density PDF. We argue that the self-similar spherical blast wave arising from expanding HII regions or stellar winds from massive stars may ultimately be responsible for creating the high mass, power-law, tail in the IMF.

preprint2009arXiv

The early history of protostellar disks, outflows, and binary stars

In star formation, magnetic fields act as a cosmic angular momentum extractor that increases mass accretion rates onto protostars and in the process, creates spectacular outflows. However, recently it has been argued that this magnetic brake is so strong that early protostellar disks -- the cradles of planet formation -- cannot form. Our three-dimensional numerical simulations of the early stages of collapse (\lesssim 10^5 yr) of overdense star--forming clouds form early outflows and have magnetically regulated and rotationally dominated disks (inside 10 AU) with high accretion rates, despite the slip of the field through the mostly neutral gas. We find that in three dimensions, magnetic fields suppress gravitationally driven instabilities which would otherwise prevent young, well ordered disks from forming. Our simulations have surprising consequences for the early formation of disks, their density and temperature structure, the mechanism and structure of early outflows, the flash heating of dust grains through ambipolar diffusion, and the origin of planets and binary stars.

preprint2008arXiv

The Growth & Migration of Jovian Planets in Evolving Protostellar Disks with Dead Zones

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the most important problems that needs to be solved in light of observations of the exosolar planets. Studies of the migration of planets in standard gas disk models routinely show that migration is too fast to form Jovian planets, and that such migrating planetary cores generally plunge into the central stars in less than a Myr. In previous work, we have shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to include dead zone evolution along with the disk, as well as planet formation via accretion of rocky and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, even in evolving disks wherein migrating planets grow by accretion, still play a fundamental role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks that are ten times more massive than a minimum mass solar nebula (MMSN) with an opacity reduction and without slowing down migration artificially. Our simulations indicate that protoplanetary disks with an initial mass comparable to the MMSN only produce Neptunian mass planets. We also find that planet migration does not help core accretion as much in the oligarchic planetesimal accretion scenario as it was expected in the runaway accretion scenario. Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting gas disks.

preprint2007arXiv

Saving Planetary Systems: Dead Zones & Planetary Migration

The tidal interaction between a disk and a planet leads to the planet's migration. A long-standing question regarding this mechanism is how to stop the migration before planets plunge into their central stars. In this paper, we propose a new, simple mechanism to significantly slow down planet migration, and test the possibility by using a hybrid numerical integrator to simulate the disk-planet interaction. The key component of the scenario is the role of low viscosity regions in protostellar disks known as dead zones, which affect planetary migration in two ways. First of all, it allows a smaller-mass planet to open a gap, and hence switch the faster type I migration to the slower type II migration. Secondly, a low viscosity slows down type II migration itself, because type II migration is directly proportional to the viscosity. We present numerical simulations of planetary migration by using a hybrid symplectic integrator-gas dynamics code. Assuming that the disk viscosity parameter inside the dead zone is (alpha=1e-4-1e-5), we find that, when a low-mass planet (e.g. 1-10 Earth masses) migrates from outside the dead zone, its migration is stopped due to the mass accumulation inside the dead zone. When a low-mass planet migrates from inside the dead zone, it opens a gap and slows down its migration. A massive planet like Jupiter, on the other hand, opens a gap and slows down inside the dead zone, independent of its initial orbital radius. The final orbital radius of a Jupiter mass planet depends on the dead zone's viscosity. For the range of alpha's noted above, this can vary anywhere from 7 AU, to an orbital radius of 0.1 AU that is characteristic of the hot Jupiters.

preprint2006arXiv

Massive star formation via high accretion rates and early disk-driven outflows

We present an investigation of massive star formation that results from the gravitational collapse of massive, magnetized molecular cloud cores. We investigate this by means of highly resolved, numerical simulations of initial magnetized Bonnor-Ebert-Spheres that undergo collapse and cooling. By comparing three different cases - an isothermal collapse, a collapse with radiative cooling, and a magnetized collapse - we show that massive stars assemble quickly with mass accretion rates exceeding 10^-3 Msol/yr. We confirm that the mass accretion during the collapsing phase is much more efficient than predicted by selfsimilar collapse solutions, i.e. dM/dt ~ c^3/G. We find that during protostellar assembly the mass accretion reaches 20 - 100 c^3/G. Furthermore, we determined the self-consistent structure of bipolar outflows that are produced in our three dimensional magnetized collapse simulations. These outflows produce cavities out of which radiation pressure can be released, thereby reducing the limitations on the final mass of massive stars formed by gravitational collapse. Moreover, we argue that the extraction of angular momentum by disk-threaded magnetic fields and/or by the appearance of bars with spiral arms significantly enhance the mass accretion rate, thereby helping the massive protostar to assemble more quickly.

preprint2005arXiv

Accretion-powered Stellar Winds as a Solution to the Stellar Angular Momentum Problem

We compare the angular momentum extracted by a wind from a pre-main-sequence star to the torques arising from the interaction between the star and its Keplerian accretion disk. We find that the wind alone can counteract the spin-up torque from mass accretion, solving the mystery of why accreting pre-main-sequence stars are observed to spin at less than 10% of break-up speed, provided that the mass outflow rate in the stellar winds is ~10% of the accretion rate. We suggest that such massive winds will be driven by some fraction $ε$ of the accretion power. For observationally constrained typical parameters of classical T-Tauri stars, $ε$ needs to be between a few and a few tens of percent. In this scenario, efficient braking of the star will terminate simultaneously with accretion, as is usually assumed to explain the rotation velocities of stars in young clusters.

Ralph E. Pudritz

What is connected

Connect this record

See the researcher in context

Building this map preview

50 published item(s)

Combined Effects of Disk Winds and Turbulence-Driven Accretion on Planet Populations

Meteorites and the RNA World: Synthesis of Nucleobases in Carbonaceous Planetesimals and the Role of Initial Volatile Content

Possible Ribose Synthesis in Carbonaceous Planetesimals

Towards RNA life on Early Earth: From atmospheric HCN to biomolecule production in warm little ponds

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

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

Formation of Planetary Populations II: Effects of Initial Disk Size & Radial Dust Drift

Formation of Planetary Populations III: Core Composition & Atmospheric Evaporation

HCN production in Titan's Atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling

Meteoritic Abundances of Fatty Acids and Potential Reaction Pathways in Planetesimals

Composition of Early Planetary Atmospheres II: Coupled Dust and Chemical Evolution in Protoplanetary Disks

Filamentary flow and magnetic geometry in evolving cluster-forming molecular cloud clumps

Meteorites and the RNA World: A Thermodynamic Model of Nucleobase Synthesis within Planetesimals

On the Formation and Chemical Composition of Super Earths

Radial Drift of Dust in Protoplanetary Disks: The Evolution of Ice lines and Dead zones

The Search for Extraterrestrial Intelligence in Earth's Solar Transit Zone

Conditions for water ice lines and Mars-mass exomoons around accreting super-Jovian planets at 1 - 20 AU from Sun-like stars

Nature's Starships II: Simulating the Synthesis of Amino Acids in Meteorite Parent Bodies

Radiation Hydrodynamics using Characteristics on Adaptive Decomposed Domains for Massively Parallel Star Formation Simulations

Seeding the Pregenetic Earth: Meteoritic Abundances of Nucleobases and Potential Reaction Pathways

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

Dust in Protoplanetary Disks: A Clue as to the Critical Mass of Planetary Cores

Nature's Starships. I. Observed Abundances and Relative Frequencies of Amino Acids in Meteorites

Planet Traps and Planetary Cores: Origins of the Planet-Metallicity Correlation

The Earliest Stages of Star and Planet Formation: Core Collapse, and the Formation of Disks and Outflows

Cluster formation in molecular clouds: I. stellar populations, star formation rates, and ionizing radiation

Disk Inhomogeneities and the Origins of Planetary System Architectures and Observational Properties

From Filamentary Networks to Dense Cores in Molecular Clouds: Toward a New Paradigm for Star Formation

Giant Molecular Cloud Formation in Disk Galaxies: Characterizing Simulated versus Observed Cloud Catalogues

Planetary Populations in the Mass-Period Diagram: A Statistical Treatment of Exoplanet Formation and the Role of Planet Traps

Evolutionary Tracks of Trapped, Accreting Protoplanets: the Origin of the Observed Mass-Period Relation

HII region variability and pre-main-sequence evolution

Shock interactions, turbulence, and the origin of the stellar mass spectrum

Simulating protostellar evolution and radiative feedback in the cluster environment

The long-term evolution of warped, magnetised discs, and precessing outflows in collapsing pre-stellar cores

Dead Zones and the Diversity of Exoplanetary Systems

Spin Evolution of Accreting Young Stars. II. Effect of Accretion-Powered Stellar Winds

The Origin of Planetary System Architectures. I. Multiple Planet Traps in Gaseous Discs

Dead Zones as Thermal Barriers to Rapid Planetary Migration in Protoplanetary Disks

Discs, outflows, and feedback in collapsing magnetized cores

Dust Settling and Rapid Planetary Migration

Simulations of the IMF in Clusters

A thermodynamic basis for prebiotic amino acid synthesis and the nature of the first genetic code

Radiatively heated, protoplanetary discs with dead zones. I. Dust settling and thermal structure of discs around M stars

Shock generated vorticity in the interstellar medium and the origin of the stellar initial mass function

The early history of protostellar disks, outflows, and binary stars

The Growth & Migration of Jovian Planets in Evolving Protostellar Disks with Dead Zones

Saving Planetary Systems: Dead Zones & Planetary Migration

Massive star formation via high accretion rates and early disk-driven outflows

Accretion-powered Stellar Winds as a Solution to the Stellar Angular Momentum Problem