Source author record

Jeffrey N. Cuzzi

Jeffrey N. Cuzzi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

17 published item(s)

preprint2022arXiv

Depletion of Moderately Volatile Elements by Open-System loss in the Early Solar Nebula

Rocky bodies of the inner solar system display a systematic depletion of the "Moderately Volatile Elements" (MVEs) that correlates with the expected condensation temperature of their likely host materials under protoplanetary nebula conditions. In this paper, we present and test a new hypothesis in which open system loss processes irreversibly remove vaporized MVEs from high nebula altitudes, leaving behind the more refractory solids residing much closer to the midplane. The MVEs irreversibly lost from the nebula through these open system loss processes are then simply unavailable for condensation onto planetesimals forming even much later, after the nebula cooled, overcoming a critical difficulty encountered by previous models of this type. We model open system loss processes operating at high nebula altitudes, such as resulting from disk winds flowing out of the system entirely, or layered accretion directly onto the young sun. We find that mass loss rates higher than found in typical T-Tauri disk winds, lasting short periods of time, are most satisfactory, pointing to multiple intense early outburst stages. Using our global nebula model, incorporating realistic particle growth and inward drift for solids, we constrain how much the MVE depletion signature in the inner region is diluted by the drift of undepleted material from the outer nebula. We also find that a significant irreversible loss of the common rock-forming elements (Fe, Mg, Si) can occur, leading to a new explanation of another longstanding puzzle of the apparent "enhancement" in the relative abundance of highly refractory elements in chondrites.

preprint2022arXiv

Global Modeling of Nebulae With Particle Growth, Drift, and Evaporation Fronts. II. The Influence of Porosity on Solids Evolution

Incremental particle growth in turbulent protoplanetary nebulae is limited by a combination of barriers that can slow or stall growth. Moreover, particles that grow massive enough to decouple from the gas are subject to inward radial drift which could lead to the depletion of most disk solids before planetesimals can form. Compact particle growth is probably not realistic. Rather, it is more likely that grains grow as fractal aggregates which may overcome this so-called radial drift barrier because they remain more coupled to the gas than compact particles of equal mass. We model fractal aggregate growth and compaction in a viscously evolving solar-like nebula for a range of turbulent intensities $α_{\rm{t}} = 10^{-5}-10^{-2}$. We do find that radial drift is less influential for porous aggregates over much of their growth phase; however, outside the water snowline fractal aggregates can grow to much larger masses with larger Stokes numbers more quickly than compact particles, leading to rapid inward radial drift. As a result, disk solids outside the snowline out to $\sim 10-20$ AU are depleted earlier than in compact growth models, but outside $\sim 20$ AU material is retained much longer because aggregate Stokes numbers there remain lower initially. Nevertheless, we conclude even fractal models will lose most disk solids without the intervention of some leap-frog planetesimal forming mechanism such as the Streaming Instability (SI), though conditions for the SI are generally never satisfied, except for a brief period %for a brief stage around $\sim 0.2$ Myr at the snowline for $α_{\rm{t}}=10^{-5}$.

preprint2022arXiv

Global Modeling of Nebulae With Particle Growth, Drift, and Evaporation Fronts. III. Redistribution of Refractories and Volatiles

Formation of the first planetesimals remains an unsolved problem. Growth by sticking must initiate the process, but multiple studies have revealed a series of barriers that can slow or stall growth, most of them due to nebula turbulence. In a companion paper, we study the influence of these barriers on models of fractal aggregate and solid, compact particle growth in a viscously evolving solar-like nebula for a range of turbulent intensities $α_{\rm{t}} = 10^{-5}-10^{-2}$. Here, we examine how disk composition in these same models changes with time. We find that advection and diffusion of small grains and vapor, and radial inward drift for larger compact particles and fractal aggregates, naturally lead to diverse outcomes for planetesimal composition. Larger particles can undergo substantial inward radial migration due to gas drag before being collisionally fragmented or partially evaporating at various temperatures. This leads to enhancement of the associated volatile in both vapor inside, and solids outside, their respective evaporation fronts, or ``snowlines''. In cases of lower $α_{\rm{t}}$, we see narrow belts of volatile or supervolatile material develop in the outer nebula, which could be connected to the bands of ``pebbles" seen by ALMA. Volatile bands, which migrate inwards as the disk cools, can persist over long timescales as their gas phase continues to advect or diffuse outward across its evaporation front. These belts could be sites where supervolatile-rich planetesimals form, such as the rare CO-rich and water-poor comets; giant planets formed just outside the H$_2$O snowline may be enhanced in water.

preprint2020arXiv

Cascade Model for Planetesimal Formation by Turbulent Clustering

We use a newly developed cascade model of turbulent concentration of particles in protoplanetary nebulae to calculate several properties of interest to the formation of primitive planetesimals and to the meteorite record. The model follows, and corrects, calculations of the primary planetesimal Initial Mass Function (IMF) by Cuzzi et al. (2010), in which an incorrect cascade model was used. Here we use the model of Hartlep et al. (2017), which has been validated against several published numerical simulations of particle concentration in turbulence. We find that, for a range of nebula and particle properties, planetesimals may be "born big", formed as sandpiles with diameters in the 10-100 km range, directly from freely floating particles. The IMFs have a modal nature, with a well-defined peak rather than a powerlaw size dependence. Predictions for the inner and outer nebula behave similarly in these regards, and observations of inner and outer nebula primitive bodies support such modal IMFs. Also, we present predictions of local particle concentrations on several lengthscales in which particles "commonly" find themselves, which have significance for meteoritical observations of the redox state and isotopic fractionation in regions of chondrule formation. An important difference between these results, and those of Cuzzi et al. (2010), is that particle growth-by-sticking must proceed to at least the 1-few cm radius range for the IMF and meteoritical properties to be most plausibly satisfied. That is, as far as the inner nebula goes, the predominant "particles" must be aggregates of chondrules (or chondrule-size precursors) rather than individual chondrules themselves.

preprint2020arXiv

Streaming Instability in Turbulent Protoplanetary Disks

The streaming instability for solid particles in protoplanetary disks is re-examined assuming the familiar alpha ($α$) model for isotropic turbulence. Turbulence always reduces the growth rates of the streaming instability relative to values calculated for globally laminar disks. While for small values of the turbulence parameter, $α< 10^{-5}$, the wavelengths of the fastest-growing disturbances are small fractions of the local gas vertical scale height $H$, we find that for moderate values of the turbulence parameter, i.e., $α\sim 10^{-5}-10^{-3}$, the lengthscales of maximally growing disturbances shift toward larger scales, approaching $H$. At these moderate turbulent intensities and for local particle to gas mass density ratios $ε< 0.5$, the vertical scales of the most unstable modes begin to exceed the corresponding radial scales so that the instability appears in the form of vertically oriented sheets extending well beyond the particle scale height. We find that for hydrodynamical turbulent disk models reported in the literature, with $α= 4\times 10^{-5} - 5\times 10^{-4}$, together with state of the art global evolution models of particle growth, the streaming instability is predicted to be viable within a narrow triangular patch of $α$--$τ_s$ parameter space centered on Stokes numbers, $τ_s \sim 0.01$ and $α\sim 4\times 10^{-5}$ and, further, exhibits growth rates on the order of several hundred to thousands of orbit times for disks with 1 percent ($Z= 0.01$) cosmic solids abundance or metallicity. Our results are consistent with, and place in context, published numerical studies of streaming instabilities.

preprint2020arXiv

The Saturn Ring Skimmer Mission Concept: The next step to explore Saturn's rings, atmosphere, interior, and inner magnetosphere

The innovative Saturn Ring Skimmer mission concept enables a wide range of investigations that address fundamental questions about Saturn and its rings, as well as giant planets and astrophysical disk systems in general. This mission would provide new insights into the dynamical processes that operate in astrophysical disk systems by observing individual particles in Saturn's rings for the first time. The Ring Skimmer would also constrain the origin, history, and fate of Saturn's rings by determining their compositional evolution and material transport rates. In addition, the Ring Skimmer would reveal how the rings, magnetosphere, and planet operate as an inter-connected system by making direct measurements of the ring's atmosphere, Saturn's inner magnetosphere and the material owing from the rings into the planet. At the same time, this mission would clarify the dynamical processes operating in the planet's visible atmosphere and deep interior by making extensive high-resolution observations of cloud features and repeated measurements of the planet's extremely dynamic gravitational field. Given the scientific potential of this basic mission concept, we advocate that it be studied in depth as a potential option for the New Frontiers program.

preprint2016arXiv

Critical layers and protoplanetary disk turbulence

A linear analysis of the zombie vortex instability is performed in a stratified shearing sheet setting for three model barotropic shear flows. The linear analysis is done by utilizing a Green's function formulation to resolve the critical layers of the associated normal-mode problem. The instability is the result of a resonant interaction between a Rossby wave and a gravity wave which we refer to as Z-modes. The associated critical layer is the location where the Doppler shifted frequency of a distant Rossby wave equals the local Brunt-Vaisala frequency. The minimum required Rossby number for instability, Ro= 0.2, is confirmed for parameter values reported in the literature. It is also found that the shear layer supports the instability in the limit where stratification vanishes. The zombie vortex instability is examined in a jet model, finding that the instability can occur for Ro= 0.05. Nonlinear vorticity forcing due to unstable Z-modes is shown to result in the creation of a jet flow at the critical layer emerging as the result of the competition between the vertical lifting of perturbation radial vorticity and the radial transport of perturbation vertical vorticity. We find that the picture of this instability leading to a form of nonlinearly driven self-replicating pattern of creation and destruction is warranted: a parent jet spawns a growing child jet at associated critical layers. A mature child jet creates a next generation of child jets at associated critical layers while simultaneously destroying its parent jet via the Rossby wave instability.

preprint2015arXiv

Global Modeling of Nebulae with Particle Growth, Drift and Evaporation Fronts. I: Methodology and Typical Results

We model particle growth in a turbulent, viscously evolving protoplanetary nebula, incorporating sticking, bouncing, fragmentation, and mass transfer at high speeds. We treat small particles using a moments method and large particles using a traditional histogram binning, including a probability distribution function of collisional velocities. The fragmentation strength of the particles depends on their composition (icy aggregates are stronger than silicate aggregates). The particle opacity, which controls the nebula thermal structure, evolves as particles grow and mass redistributes. While growing, particles drift radially due to nebula headwind drag. Particles of different compositions evaporate at "evaporation fronts" (EFs) where the midplane temperature exceeds their respective evaporation temperatures. We track the vapor and solid phases of each component, accounting for advection and radial and vertical diffusion. We present characteristic results in evolutions lasting $2 \times 10^5$ years. In general, (a) mass is transferred from the outer to inner nebula in significant amounts, creating radial concentrations of solids at EFs, (b) particle sizes are limited by a combination of fragmentation, bouncing, and drift, (c) "lucky" large particles never represent a significant amount of mass, and (d) restricted radial zones just outside each EF become compositionally enriched in the associated volatiles. We point out implications for mm-submm SEDs and inference of nebula mass, radial banding, the role of opacity on new mechanisms for generating turbulence, enrichment of meteorites in heavy oxygen isotopes, variable and nonsolar redox conditions, primary accretion of silicate and icy planetesimals, and the makeup of Jupiter's core.

preprint2015arXiv

New Paradigms For Asteroid Formation

Asteroids and meteorites provide key evidence on the formation of planetesimals in the Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation within a population of initially km-scale planetesimals. However, new models challenge this idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the gravitational collapse of small particles which have organised themselves in dense filaments and clusters in the turbulent gas. Particles concentrate passively between eddies down to the smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The streaming instability causes particles to take an active role in the concentration, by piling up in dense filaments whose friction on the gas reduces the radial drift compared to that of isolated particles. In this chapter we review new paradigms for asteroid formation and compare critically against the observed properties of asteroids as well as constraints from meteorites. Chondrules of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. We review recent progress on understanding the incorporation of chondrules into the asteroids, including layered accretion models where chondrules are accreted onto asteroids over millions of years. We highlight in the end ten unsolved questions in asteroid formation where we expect that progress will be made over the next decade.

preprint2015arXiv

The Gravitational instability of solids assisted by gas drag: slowing by turbulent mass diffusivity

The Goldreich and Ward (1973) (axisymmetric) gravitational instability of a razor thin particle layer occurs when the Toomre parameter $Q_T \equiv c_p Ω_0 / πG Σ_p < 1$ ($c_p$ being the particle dispersion velocity). Ward(1976,2000) extended this analysis by adding the effect of gas drag upon particles and found that even when $Q_T > 1$, sufficiently long waves were always unstable. Youdin (2005a,b) carried out a detailed analysis and showed that the instability allows chondrule-sized ($\sim 1 $ mm) particles to undergo radial clumping with reasonable growth times even in the presence of a moderate amount of turbulent stirring. The analysis of Youdin includes the role of turbulence in setting the thickness of the dust layer and in creating a turbulent particle pressure in the momentum equation. However, he ignores the effect of turbulent mass diffusivity on the disturbance wave. Here we show that including this effect reduces the growth-rate significantly, by an amount that depends on the level of turbulence, and reduces the maximum intensity of turbulence the instability can withstand by 1 to 3 orders of magnitude. The instability is viable only when turbulence is extremely weak and the solid to gas surface density of the particle layer is considerably enhanced over minimum-mass-nebula values. A simple mechanistic explanation of the instability shows how the azimuthal component of drag promotes instability while the radial component hinders it. A gravito-diffusive overstability is also possible but never realized in the nebula models.

preprint2015arXiv

The Spherically Symmetric Gravitational Collapse of a Clump of Solids in a Gas

Several mechanisms have been identified that create dense particle clumps in the solar nebula. The present work is concerned with the gravitational collapse of such clumps, idealized as being spherically symmetric. Calculations using the two-fluid model are performed (almost) up to the time when a central density singularity forms. The end result of the study is a parametrization for this time, in order that it may be compared with timescales for various disruptive effects to which clumps may be subject. An important effect is that as the clump compresses, it also compresses the gas due to drag. This increases gas pressure which retards particle collapse and leads to oscillation in the size and density of the clump. The ratio of gravitational force to gas pressure gives a two-phase Jeans parameter, $J_t$, which is the classical Jeans parameter with the sound speed replaced by an the wave speed in a coupled two-fluid medium. Its use makes the results insensitive to the initial density ratio of particles to gas as a separate parameter. An ordinary differential equation model is developed which takes the form of two coupled non-linear oscillators and reproduces key features of the simulations. Finally, a parametric study of the time to collapse is performed and a formula (fit to the simulations) is developed. In the incompressible limit $J_t \to 0$, collapse time equals sedimentation time. As $J_t$ increases, the collapse time decreases roughly linearly with $J_t$ until $J_t \gtrsim 0.4$ when it becomes approximately equal to the dynamical time.

preprint2014arXiv

Cloud Base Signature in Transmission Spectra of Exoplanet Atmospheres

We present an analytical model for the transmission spectrum of a transiting exoplanet, showing that a cloud base can produce an observable inflection point in the spectrum. The wavelength and magnitude of the inflection can be used to break the degeneracy between the atmospheric pressure and the abundance of the main cloud material, however, the abundance still depends on cloud particle size. An observed inflection also provides a specific point on the atmospheric P-T profile, giving us a "thermometer" to directly validate or rule out postulated cloud species. We apply the model to the transit spectrum of HD 189733b.

preprint2013arXiv

Clouds and Hazes in Exoplanet Atmospheres

Clouds and hazes are commonplace in the atmospheres of solar system planets and are likely ubiquitous in the atmospheres of extrasolar planets as well. Clouds affect every aspect of a planetary atmosphere, from the transport of radiation, to atmospheric chemistry, to dynamics and they influence - if not control - aspects such as surface temperature and habitability. In this review we aim to provide an introduction to the role and properties of clouds in exoplanetary atmospheres. We consider the role clouds play in influencing the spectra of planets as well as their habitability and detectability. We briefly summarize how clouds are treated in terrestrial climate models and consider the far simpler approaches that have been taken so far to model exoplanet clouds, the evidence for which we also review. Since clouds play a major role in the atmospheres of certain classes of brown dwarfs we briefly discuss brown dwarf cloud modeling as well. We also review how the scattering and extinction efficiencies of cloud particles may be approximated in certain limiting cases of small and large particles in order to facilitate physical understanding. Since clouds play such important roles in planetary atmospheres, cloud modeling may well prove to be the limiting factor in our ability to interpret future observations of extrasolar planets.

preprint2013arXiv

Utilitarian Opacity Model for Aggregate Particles in Protoplanetary Nebulae and Exoplanet Atmospheres

As small solid grains grow into larger ones in protoplanetary nebulae, or in the cloudy atmospheres of exoplanets, they generally form porous aggregates rather than solid spheres. A number of previous studies have used highly sophisticated schemes to calculate opacity models for irregular, porous particles with size much smaller than a wavelength. However, mere growth itself can affect the opacity of the medium in far more significant ways than the detailed compositional and/or structural differences between grain constituents once aggregate particle sizes exceed the relevant wavelengths. This physics is not new; our goal here is to provide a model that provides physical insight and is simple to use in the increasing number of protoplanetary nebula evolution, and exoplanet atmosphere, models appearing in recent years, yet quantitatively captures the main radiative properties of mixtures of particles of arbitrary size, porosity, and composition. The model is a simple combination of effective medium theory with small-particle closed-form expressions, combined with suitably chosen transitions to geometric optics behavior. Calculations of wavelength-dependent emission and Rosseland mean opacity are shown and compared with Mie theory. The model's fidelity is very good in all comparisons we have made, except in cases involving pure metal particles or monochromatic opacities for solid particles with size comparable to the wavelength.

preprint2010arXiv

Cassini imaging search rules out rings around Rhea

We have conducted an intensive search for any material that may orbit Rhea, using images obtained by the Cassini ISS narrow-angle camera. We find no evidence for any such material, contradicting an earlier and surprising inference that Rhea, the second-largest moon of Saturn, possesses a system of narrow rings embedded in a broad circum-satellite disk or cloud (Jones et al. 2008, Science). If one accepts the calculations of the previous authors, our results would require any narrow rings around Rhea to be composed of objects no smaller than 10 meters in radius, in order to be massive enough to account for the observed charged-particle absorptions but diffuse enough to evade detection in our images. Such a particle-size distribution is unrealistic given that erosion processes will break down 10-meter objects to smaller sizes, which would then have been seen in our images. Furthermore, Jones et al. (2008) assumed that the absorption of electrons by hypothetical circum-Rhea particles is proportional to the volume of the particles, even when they are much larger than the electron penetration depth (van Allen 1983, 1987, JGR). If we assume instead that large particles can only absorb electrons efficiently near their surfaces, then our results rule out any solid material orbiting Rhea as the explanation for the absorptions detected by Jones et al. (2008), regardless of particle size, by a margin of four (for narrow rings) or two (for a broad diffuse cloud) orders of magnitude. We conclude that the anomalous electron absorptions reported by Jones et al. (2008) are the signature of a new and hitherto unknown magnetospheric process.

preprint2010arXiv

Towards Initial Mass Functions for Asteroids and Kuiper Belt Objects

Our goal is to understand primary accretion of the first planetesimals. The primitive meteorite record suggests that sizeable planetesimals formed in the asteroid belt over a period longer than a million years, each composed entirely of an unusual, but homogeneous, mixture of mm-size particles. We sketch a scenario in which primary accretion of 10-100km size planetesimals proceeds directly, if sporadically, from aerodynamically-sorted mm-size particles (generically "chondrules"). These planetesimal sizes are in general agreement with the currently observed asteroid mass peak near 100km diameter, which has been identified as a "fossil" property of the pre-erosion, pre-depletion population. We extend our primary accretion theory to make predictions for outer solar system planetesimals, which may also have a preferred size in the 100km diameter range. We estimate formation rates of planetesimals and assess the conditions needed to match estimates of both asteroid and Kuiper Belt Object (KBO) formation rates. For nebula parameters that satisfy observed mass accretion rates of Myr-old protoplanetary nebulae, the scenario is roughly consistent with not only the "fossil" sizes of the asteroids, and their estimated production rates, but also with the observed spread in formation ages of chondrules in a given chondrite, and with a tolerably small radial diffusive mixing during this time between formation and accretion (the model naturally helps explain the peculiar size distribution of chondrules within such objects). The scenario also produces 10-100km diameter primary KBOs. The optimum range of parameters, however, represents a higher gas density and fractional abundance of solids, and a smaller difference between keplerian and pressure-supported orbital velocities, than "canonical" models of the solar nebula. We discuss several potential explanations for these differences.

preprint2007arXiv

A Cascade Model for Particle Concentration and Enstrophy in Fully Developed Turbulence with Mass Loading Feedback

A cascade model is described based on multiplier distributions determined from 3D direct numerical simulations (DNS) of turbulent particle laden flows, which include two-way coupling between the phases at global mass loadings equal to unity. The governing Eulerian equations are solved using pseudo-spectral methods on up to 512**3 computional grid points. DNS results for particle concentration and enstrophy at Taylor microscale Reynolds numbers in the range 34 - 170 were used to directly determine multiplier distributions (PDFs) on spatial scales 3 times the Kolmogorov length scale. The width of the PDFs, which is a measure of intermittency, decreases with increasing mass loading within the local region where the multipliers are measured. The functional form of this dependence is not sensitive to Reynolds numbers in the range considered. A partition correlation probability is included in the cascade model to account for the observed spatial anticorrelation between particle concentration and enstrophy. Joint probability distribution functions of concentration and enstrophy generated using the cascade model are shown to be in excellent agreement with those derived directly from our 3D simulations. Probabilities predicted by the cascade model are presented at Reynolds numbers well beyond what is achievable by direct simulation. These results clearly indicate that particle mass loading significantly reduces the probabilities of high particle concentration and enstrophy relative to those resulting from unloaded runs. Particle mass density appears to reach a limit at around 100 times the gas density. This approach has promise for significant computational savings in certain applications.