Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “protoplanetary disks”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Initial Conditions of Planet Formation: Time Constraints from Small Bodies and the Lifetime of Reservoirs in the Solar Protoplanetary Disk

This review explores the timescales of the initial phase of planet formation, from nebular dust (CAIs and chondrules) to planetesimal accretion and differentiation, using evidence from meteorite research. Aluminium-Mg systematics of CAIs indicate either an extended period of CAI formation (∼0.3 Ma) or an initial 26 Al heterogeneity, with evidence supporting a homogeneous 26 Al abundance in the protoplanetary disk. Thermal and aqueous alteration on the parent body can disturb the U-Pb and Al-Mg chronometers in chondrules. Focusing on relatively robust isochron data from plagioclase of pristine (types ≤3.05) chondrites indicates a shift in chondrule formation locations, moving from the inner to the outer disk over time. Ages of basaltic achondrites show that silicate differentiation on small bodies was well underway within the first few million years (Ma) of our solar system. Their age record, however, reveals inconsistencies between different chronometers, partly caused by (i) secondary disturbances due to thermal metamorphism, aqueous alteration, or impacts, (ii) the presence of xenolithic minerals, and (iii) potentially variable initial 26 Al abundances due to disturbances at the mineral scale. Nucleosynthetic isotope data indicate that parent bodies of iron and stony meteorites formed in two distinct regions within the protoplanetary disk: the inner, non-carbonaceous (NC) and the outer, carbonaceous (CC) region. Based on Hf-W chronometry it has been demonstrated that NC and CC parent bodies of magmatic iron meteorites segregated their cores within ∼1–3 Ma after CAI formation, implying that parent body accretion occurred within <1 Ma in both reservoirs. Combining accretion ages with nucleosynthetic data further reveals that, at first order, NC and CC reservoirs in the solar protoplanetary disk were established within 1 Ma and existed over several Ma with limited exchange between them. In the CR chondrite accretion region of the disk, planetary bodies formed over at least 3 Ma, while in most other regions, formation spanned at least 1 Ma, with minimal changes in nucleosynthetic isotope compositions. Aerodynamical size sorting of dust likely introduced or amplified some of these variations.

26Al distribution↗
DOE OSTI2026

The Extent of Solar Energetic Particle Irradiation in the Sun’s Protoplanetary Disk

Solar flares emit X-rays and high-energy (MeV to GeV) ions called solar energetic particles (SEPs). Astronomical observations show solar-mass protostellar fluxes are a factor Φ ≈ 3 × 10 2 –3 × 10 3 times higher than the present-day Sun. Constraining Φ in the early solar system is important for modeling ionization in the Sun’s protoplanetary disk, the extent of magnetorotational instability or magnetocentrifugal outflows, or even production of short-lived radionuclides. Recent interpretations of meteoritic data—including cosmogenic Ne in hibonite grains, initial ( 10 Be/ 9 Be) 0 ratios in Ca-rich, Al-rich inclusions (CAIs), or even inferences of live 7 Be in CAIs—suggest values Φ > 10 5 , reaching as high as Φ ≈ 6 × 10 6 , which would make the young Sun extraordinarily active, even for a protostar. We constrain Φ by reexamining these data. We conclude that cosmogenic Ne was produced in hibonite grains as they resided in the disk; 36 Cl was created in Cl-poor grains after the disk dissipated; 10 Be was inherited from the molecular cloud, with almost no (<1%) 10 Be created in the disk; and there is no evidence whatsoever for any live 7 Be in CAIs. We show these data are consistent with a value Φ ≈ 3 × 10 3 for the first >5 Myr of the solar nebula. The early Sun evidently emitted a flux of X-rays and SEPs not atypical for a protostar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Zonal Flow and Self-regulating Mechanism in a Hydrodynamic Disk

This study addresses key aspects of momentum transport in hydrodynamic disks, which is critical for understanding zonal flow generation and turbulence in compressible hydrodynamic disks. We find that nonlinear momentum/density transport leads to the formation of zonal flows from the Rossby wave instability in disks. We analytically derive the generation and location of zonal flows and describe a modified Taylor identity applicable to compressible disk flows. We further present a self-regulation model, revealing a dynamic interplay between zonal flow and fluctuations driven by Rossby wave instability that regulates the nonlinear saturation state. This theoretical framework contributes insights into the dynamics of disks such as protoplanetary disks, shedding light on the intricate processes governing momentum/density transport and the emergence of zonal flows in the saturation of protoplanetary disks.

Astrophysical fluid dynamics↗
DOE OSTI2025

The variety and origin of materials accreted by Bennu’s parent asteroid

The first bodies to form in the Solar System acquired their materials from stars, the presolar molecular cloud and the protoplanetary disk. Asteroids that have not undergone planetary differentiation retain evidence of these primary accreted materials. However, geologic processes such as hydrothermal alteration can dramatically change their bulk mineralogy, isotopic compositions and chemistry. Here we analyse the elemental and isotopic compositions of samples from asteroid Bennu to uncover the sources and types of material accreted by its parent body. We show that some primary accreted materials escaped the extensive aqueous alteration that occurred on the parent asteroid, including presolar grains from ancient stars, organic matter from the outer Solar System or molecular cloud, refractory solids that formed close to the Sun, and dust enriched in neutron-rich Ti isotopes. We find Bennu to be richer in isotopically anomalous organic matter, anhydrous silicates, and light isotopes of K and Zn than its closest compositional counterparts, asteroid Ryugu and Ivuna-type (CI) carbonaceous chondrite meteorites. We propose that the parent bodies of Bennu, Ryugu and CI chondrites formed from a common but spatially and/or temporally heterogeneous reservoir of materials in the outer protoplanetary disk.

Kuiper belt↗
DOE OSTI2025

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Does the Sun have a dark disk?

The Sun is not quite a perfect sphere, and its oblateness, thought to be induced through its rotation, has been measured using optical observations of its radius. Its gravitational quadrupole moment can then be deduced using solar models, or through helioseismology, and it can also be determined from measurements of its gravitational effects on Mercury’s orbit. The various assessments do not appear to agree, with the most complete and precise orbital assessments being in slight excess of other determinations. This may speak to the existence of a nonluminous disk or ring, where we also note evidence for a circumsolar dust ring within Mercury’s orbit from the Solar TErrestrial RElations Observatory (STEREO) mission. Historically, too, a protoplanetary disk may have been key to reconciling the Sun’s metallicity with its neutrino yield. The distribution of the nonluminous mass within Mercury’s orbit can modify the relative size of the optical and orbital quadrupole moments in different ways. We develop how we can use these findings to limit the mass of a dark disk, ring, or halo in the immediate vicinity of the Sun, and we note how future observational studies of the inner Solar System can not only refine these constraints but can also help to identify and to assess the mass of its dark-matter component. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A collective trigger for widespread planetesimal formation revealed by accretion ages

The formation of planetesimals was an integral part of the cascading series of processes that built the terrestrial planets. To illuminate planetesimal formation, here we develop a refined thermal evolution model to calculate the formation ages of meteorite parent planetesimals. This model includes chemical reactions and phase changes during heating, as well as natural variations in the proportions of the constituent phases of these planetesimals. We find that the parent bodies of non-carbonaceous (NC) and carbonaceous (CC) iron meteorites start forming at very similar times ( ∼ 0.95 Myr after calcium-aluminium-rich inclusion [CAI] formation) and occupy overlapping time windows. NC and CC chondrite parent bodies formed later during non-overlapping periods. We combine these ages with proportions of isotopic end-members we recover from mixing models to construct records of motion throughout the protoplanetary disk. These records argue that NC and CC material traversed the barrier in the disk after ∼ 0.95 Myr after CAI formation. The onset of this motion coincided with planetesimal formation, indicating that the phenomenon that drove motion also triggered planetesimal formation. We argue that this feature also served as the semi-permeable barrier in the disk. Although its identity is uncertain, the effects this phenomenon had on the timing of planetesimal formation and motion through the disk can now serve as constraints on models of disk evolution. Models that reproduce these effects would elucidate the nature and implications of this phenomenon, which is key to unlocking a holistic model of terrestrial planet building.

Astronomy and AstroPhysics↗
DOE OSTI2025

Abundant ammonia and nitrogen-rich soluble organic matter in samples from asteroid (101955) Bennu

Organic matter in meteorites reveals clues about early Solar System chemistry and the origin of molecules important to life, but terrestrial exposure complicates interpretation. Samples returned from the B-type asteroid Bennu by the Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer mission enabled us to study pristine carbonaceous astromaterial without uncontrolled exposure to Earth’s biosphere. Here we show that Bennu samples are volatile rich, with more carbon, nitrogen and ammonia than samples from asteroid Ryugu and most meteorites. Nitrogen-15 isotopic enrichments indicate that ammonia and other N-containing soluble molecules formed in a cold molecular cloud or the outer protoplanetary disk. We detected amino acids (including 14 of the 20 used in terrestrial biology), amines, formaldehyde, carboxylic acids, polycyclic aromatic hydrocarbons and N-heterocycles (including all five nucleobases found in DNA and RNA), along with ~10,000 N-bearing chemical species. All chiral non-protein amino acids were racemic or nearly so, implying that terrestrial life’s left-handed chirality may not be due to bias in prebiotic molecules delivered by impacts. The relative abundances of amino acids and other soluble organics suggest formation and alteration by low-temperature reactions, possibly in NH 3 -rich fluids. Bennu’s parent asteroid developed in or accreted ices from a reservoir in the outer Solar System where ammonia ice was stable.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Building wet planets through high-pressure magma–hydrogen reactions

Close-in transiting sub-Neptunes are abundant in our Galaxy. Planetary interior models based on their observed radius–mass relationship suggest that sub-Neptunes contain a discernible amount of either hydrogen (dry planets) or water (wet planets) blanketing a core composed of rocks and metal. Water-rich sub-Neptunes have been believed to form farther from the star and then migrate inwards to their present orbits. Here we report experimental evidence of reactions between warm, dense hydrogen fluid and silicate melt that release silicon from the magma to form alloys and hydrides at high pressures. We found that oxygen liberated from the silicate melt reacts with hydrogen, producing an appreciable amount of water up to a few tens of weight per cent, which is much greater than previously predicted based on low-pressure ideal gas extrapolation. Consequently, these reactions can generate a spectrum of water contents in hydrogen-rich planets, with the potential to reach water-rich compositions for some sub-Neptunes, implying an evolutionary relationship between hydrogen-rich and water-rich planets. Therefore, detection of a large amount of water in exoplanet atmospheres may not be the optimal evidence for planet migration in the protoplanetary disk, calling into question the assumed link between composition and planet formation location.

Horn, H. W. [Arizona State University, Tempe, AZ (↗
DOE OSTI2025

The formation and structure of iron-dominated planetesimals

Metal-rich asteroids and iron meteorites are considered core remnants of differentiated planetesimals and/or products of oxygen-depleted accretion. Investigating the origins of iron-rich planetesimals could provide key insights into planet formation mechanisms. Using differentiation models, we evaluate the interior structure and composition of representative-sized planetesimals (~200 km diameter), while varying oxygen fugacity and initial bulk meteoritic composition. Under the oxygen-poor conditions that likely existed early in the inner regions of the Solar System and other protoplanetary disks, core fractions remain relatively consistent across a range of bulk compositions (CI, H, EH, and CBa). Some of these cores could incorporate significant amounts of silicon (10–30 weight%) and explain the metal fractions of Fe-rich bodies in the absence of mantle stripping. Conversely, planetesimals forming under more oxidizing conditions, such as beyond snow lines, could exhibit smaller cores, enriched in carbon, sulfur (>1 wt%), and oxides. Sulfur-rich cores, like those formed from EH and H bulk compositions, could remain partly molten, sustain dynamos, and even drive sulfur-rich volcanism. Additionally, bodies with high carbon contents, such as CI compositions, can form graphitic outer layers. These variations highlight the importance of initial formation conditions in shaping planetesimal structures. Future missions, such as NASA’s Psyche mission, offer an opportunity to measure the relative abundances of key elements (Fe, Ni, Si, and S) necessary to distinguish among formation scenarios and structure models for Fe-rich and reduced planetesimals.

meteorites↗
DOE OSTI2025

exoALMA. V. Gaseous Emission Surfaces and Temperature Structures

An analysis of the gaseous component in protoplanetary disks can inform us about their thermal and physical structure, chemical composition, and kinematic properties, all of which are crucial for understanding various processes within the disks. By exploiting the asymmetry of the line emission, or via line profile analysis, we can locate the emitting surfaces. Here, we present the emission surfaces of the exoALMA sources in 12 CO J = 3–2, 13 CO J = 3–2, and CS J = 7–6. We find that 12 CO traces the upper disk atmosphere, with mean $\langle$ z/r $\rangle$ values of ≈0.28, while 13 CO and CS trace lower regions of the disk with mean $\langle$ z/r $\rangle$ values of ≈0.16 and ≈0.18, respectively. We find that 12 CO $\langle$ z/r $\rangle$ and the disk mass are positively correlated with each other; this relationship offers a straightforward way to infer the disk mass. We derive 2D r – z temperature distributions of the disks. Additionally, we search for substructure in the surfaces and radial intensity profiles; we find evidence of localized substructure in the emission surfaces and peak intensity profiles of nearly every disk, with this substructure often being coincident between molecular tracers, intensity profiles, and kinematic perturbations. Four disks display evidence of potential photodesorption, implying that this effect may be common even in low far-ultraviolet star-forming regions. For most disks, we find that the physical and thermal structure is more complex than analytical models can account for, highlighting a need for more theoretical work and a better understanding of the role of projection effects on our observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2026

Isotopic Evidence for A Cold and Distant Origin of 3I/Atlas

Interstellar objects provide the only directly observable samples of icy planetesimals formed around other stars, and can therefore provide insight into the diversity of physical and chemical conditions occurring during exoplanet formation1−3. Here we report isotopic measurements of the interstellar comet 3I/ATLAS, which reveal an elemental composition unlike any Solar System body. The water in 3I/ATLAS is enriched in deuterium, at a level of D/H = (0.98 ± 0.06)%, which is more than an order of magnitude higher than in known comets, while its range of 12 C/ 13 C ratios (141–191 for CO 2 and 123–172 for CO) exceeds typical values found in the Solar System, as well as nearby interstellar clouds and protoplanetary disks. Such extreme isotopic signatures indicate formation at temperatures ≲ 30 K in a relatively metal-poor environment. When interpreted with respect to models for Galactic chemical evolution, the carbon isotopic composition implies that 3I/ATLAS may have accreted as long ago as 12 billion years, following a period of intense, early star formation. 3I/ATLAS thus represents a preserved fragment of an ancient planetary system.

Martin A Cordiner↗
NASA NTRS2026

Isotopic Signature of Organic Molecules from Beyond the Solar System: An Enriched Methane D/H Ratio in the Interstellar Object 3I/ATLAS

Interstellar objects are interlopers from other planetary systems, and their volatile compositions provide a glimpse into planet formation around their host star. We present near-infrared spectra of the coma of interstellar object 3I/ATLAS measured with the James Webb Space Telescope. Our results demonstrate an unexpectedly high D/H = (3.31 ± 0.34)% for methane and represent an exceedingly rare detection of deuterated organic molecules in an interstellar object. This D/H ratio is a factor of 14 ± 2 higher than that measured in comet 67P/Churyumov-Gerasimenko by the Rosetta spacecraft, the only other comet for which CH3D has been detected, yet the ratio of deuteration in methane compared with water is consistent for both comets within 1.2σ. The D/H ratio in methane is observationally unconstrained in extrasolar sources to date, but the enriched ratio in 3I/ATLAS is most similar to those measured in other organic molecules toward primitive environments. The high D/H ratios of water and methane in 3I/ATLAS are a natural consequence of formation in a high D/H elemental ratio environment as a result of locally cold conditions in the protoplanetary disk and prior interstellar cloud. Thus, 3I/ATLAS formed in an environment very different from that in which our Sun and planets originated.

Interstellar Object 3I/ATLAS↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dust supply to close binary systems

Binary systems can be born surrounded by circumbinary discs. The gaseous disc around either of the two stellar companions can have its life extended by the supply of mass arriving from the circumbinary disc. The objective of this study is to investigate the gravitational interactions exerted by a compact and eccentric binary system on the circumbinary and circumprimary discs, and the resulting transport of gas and solids between the disc components. We assume that the gas in the system behaves as a fluid, and we model its evolution by means of high-resolution hydrodynamical simulations. Dust grains are modelled as Lagrangian particles that interact with the gas and the stars. Our models indicate that significant fluxes of gas and dust proceed from the circumbinary disc towards the circumprimary disc. For the applied system parameters, grains of certain sizes are segregated outside the tidal gap generated by the stars. Consequently, the size distribution of the transported dust is not continuous, but presents a gap in the millimetre size range. In close binaries, the lifetime of an isolated circumprimary disc is found to be short, ~10 5 years, because of its low mass. However, because of the influx of gas from beyond the tidal gap, the disc around the primary star can survive much longer, ~10 6 years, as long as gas accretion from the circumbinary disc continues. The supply of solids and the extended lifetime of a circumbinary disc also aids in the possible formation of giant planets. Compared to close binary systems without a circumbinary disc, we expect a higher frequency of singleplanet or multiple-planet systems. Additionally, a planetesimal or debris belt can form in the proximity of the truncation radius of the circumprimary disc and/or around the location of the exterior edge of the tidal gap.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. II. Locality of Gravitoturbulence, Clumpy Spirals, and Implications for Observable Substructure

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local α viscosity to zeroth order in aspect ratio (H/R), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal{O}$(α 1/2 ) deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in H/R and α 1/2 . The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different m. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with m-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

Gravitational instability↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. I. Stochastic Fragmentation with Optical-depth-dependent Rate and Universal Fragment Mass

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area p frag showing an exponential dependence on the parameter β ≡ Ω K t cool , where Ω K is the Keplerian rotation frequency, and t cool is the average cooling timescale. Compared to a prescribed constant β, radiative cooling in the optically thin/thick regime makes p frag decrease slower/faster in β; the critical β corresponding to ∼1 fragment per orbit is ≈3, 5, and 2 for constant β, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is m frag ≈ 40M tot h 3 , with M tot being the total (star+disk) mass and h = H/R being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find m frag ∼ 20M J , suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the m frag distribution, forming massive planets is also possible.

gravitational instability↗