BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Morris Podolak

Morris Podolak 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

Catalog footprint

What is connected

17works
1topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

On the structure and long-term evolution of ice-rich bodies

The interest in the structure of ice-rich planetary bodies, in particular the differentiation between ice and rock, has grown due to the discovery of Kuiper belt objects and exoplanets. We thus carry out a parameter study for a range of planetary masses $M$, yielding radii $50 \aplt R \aplt 3000$~km, and for rock/ice mass ratios between 0.25 and 4, evolving them for 4.5~Gyr in a cold environment, to obtain the present structure. We use a thermal evolution model that allows for liquid and vapor flow in a porous medium, solving mass and energy conservation equations under hydrostatic equilibrium for a spherical body in orbit around a central star. The model includes the effect of pressure on porosity and on the melting temperature, heating by long-lived radioactive isotopes, and temperature-dependent serpentinization and dehydration. We obtain the boundary in parameter space [size, rock-content] between bodies that differentiate, forming a rocky core, and those which remain undifferentiated: small bodies, bodies with a low rock content, and the largest bodies considered, which develop high internal pressures and barely attain the melting temperature. The final differentiated structure comprises a rocky core, an ice-rich mantle, and a thin dense crust below the surface. We obtain and discuss the bulk density-radius relationship. The effect of a very cold environment is investigated and we find that at an ambient temperature of $\sim$20~K, small bodies preserve the ice in amorphous form to the present.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Random Models for Exploring Planet Compositions I: Uranus as an Example

Modeling the interior of a planet is difficult because the small number of measured parameters is insufficient to constrain the many variables involved in describing the interior structure and composition. One solution is to invoke additional constraints based on arguments about how the planet formed. However, a planet's actual structure and composition may hold clues to its formation which would be lost if this structure were not allowed by the initial assumptions. It is therefore interesting to explore the space of allowable compositions and structures in order to better understand which cosmogonic constraints are absolutely necessary. To this end, we describe a code for generating random, monotonic, density distributions, rho(r), that fit a given mass, radius, and moment of inertia. Integrating the equation of hydrostatic equilibrium gives the pressure, P(r), at each point in the body. We then provide three algorithms for generating a monotonic temperature distribution, T(r), and an associated composition that is consistent with the rho-P relation, and realistic equations of state. We apply this code to Uranus as a proof of concept, and show that the ratio of rock to water cannot be much larger than 2.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Effect of Accretion Rate and Composition on the Structure of Ice-rich Super-Earths

It is reasonable to assume that the structure of a planet and the interior distribution of its components are determined by its formation history. We thus follow the growth of a planet from a small embryo through its subsequent evolution. We estimate the accretion rate range based on a protoplanetary disk model at a large enough distance from the central star, for water ice to be a major component. We assume the accreted material to be a mixture of silicate rock and ice, with no H-He envelope, as the accretion timescale is much longer than the time required for the nebular gas to dissipate. We adopt a thermal evolution model that includes accretional heating, radioactive energy release, and separation of ice and rock. Taking the Safronov parameter and the ice-to-rock ratio as free parameters, we compute growth and evolutionary sequences for different parameter combinations, for 4.6 Gyr. We find the final structure to depend significantly on both parameters. Low initial ice to rock ratios and high accretion rates, each resulting in increased heating rate, lead to the formation of extended rocky cores, while the opposite conditions leave the composition almost unchanged and result in relatively low internal temperatures. When rocky cores form, the ice-rich outer mantles still contain rock mixed with the ice. We find that a considerable fraction of the ice evaporates upon accretion, depending on parameters, and assume it is lost, thus the final surface composition and bulk density of the planet do not necessarily reflect the protoplanetary disk composition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detailed Calculations of the Efficiency of Planetesimal Accretion in the Core-Accretion Model

We present results of a detailed study of the rate of the accretion of planetesimals by a growing proto-Jupiter in the core-accretion model. Using a newly developed code, we accurately combine a detailed three-body trajectory calculation with gas drag experienced during the passage of planetesimals in the protoplanet's envelope. We find that the motion of planetesimals is excited to the extent that encounters with the proto-planetary envelope become so fast that ram pressure breaks up the planetesimals in most encounters. As a result, the accretion rate is largely independent of the planetesimal size and composition. For the case we explored of a planet forming at 5.2 AU from the Sun in a disk with a solid surface density of 6 g/cm^2 (Lozovsky et al. 2017) the accretion rate we compute differs in several respects from that assumed by those authors. We find that only 4-5 M_Earth is accreted in the first 1.5x10^6 years before the onset of rapid gas accretion. Most of the mass, some 10 M_Earth, is accreted simultaneously with this rapid gas accretion. In addition, we find that the mass accretion rate remains small, but non-zero for at least a million years after this point, and an additional 0.3-0.4 M_Earth is accreted during that time. This late accretion, together with a rapid infall of gas could lead to the accreted material being mixed throughout the outer regions, and may account for the enhancement of high-Z material in Jupiter's envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Simple Model for Understanding the DIM Dust Measurement at Comet 67P/Churyumov-Gerasimenko

We present a simple model for gas and dust flow from 67P/Churyumov-Gerasimenko that can be used to understand the grain impact observed by the DIM instrument on Philae (Krueger et al. 2015) We show how model results when applied to the GIADA measurements (Rotundi et al. 2015; Della Corte et al. 2015) can be used, in conjunction with the results found by the MIRO (Schloerb et al. 2015) and VIRTIS (De Sanctis et al. 2015) instruments to infer surface properties such as surface temperature and surface ice fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Capture and Evolution of Planetesimals in Circumjovian Disks

We study the evolution of planetesimals in evolved gaseous disks, which orbit a solar-mass star and harbor a Jupiter-mass planet at a_p~5AU. The gas dynamics is modeled with a three-dimensional hydrodynamics code that employes nested-grids and achieves a resolution of one Jupiter's radius in the circumplanetary disk. The code models solids as individual particles. Planetesimals are subjected to gravitational forces by the star and the planet, drag force by the gas, disruption via ram pressure, and mass loss through ablation. The mass evolution of solids is calculated self-consistently with their temperature, velocity, and position. We consider icy and icy/rocky bodies of radius 0.1-100km, initially deployed on orbits around the star within a few Hill radii (Rhill) of the planet's orbit. Planetesimals are scattered inward, outward, and toward disk regions of radius r>>a_p. Scattering can relocate significant amounts of solids, provided that regions |r-a_p|~ 3Rhill are replenished with planetesimals. Scattered bodies can be temporarily captured on planetocentric orbits. Ablation consumes nearly all solids at gas temperatures > ~220K. Super-keplerian rotation around and beyond the outer edge of the gas gap can segregate < ~0.1km bodies, producing solid gap edges at size-dependent radial locations. Capture, break-up, and ablation of solids result in a dust-laden circumplanetary disk with low surface densities of km-size planetesimals, implying relatively long timescales for satellite formation. After a giant planet acquires most of its mass, accretion of solids is unlikely to alter significantly its heavy-element content. The luminosity generated by solids' accretion can be of a similar order of magnitude to the contraction luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convection and Mixing in Giant Planet Evolution

The primordial internal structures of gas giant planets are unknown. Often giant planets are modeled under the assumption that they are adiabatic, convective, and homogeneously mixed, but this is not necessarily correct. In this work, we present the first self-consistent calculation of convective transport of both heat and material as the planets evolve. We examine how planetary evolution depends on the initial composition and its distribution, whether the internal structure changes with time, and if so, how it affects the evolution. We consider various primordial distributions, different compositions, and different mixing efficiencies and follow the distribution of heavy elements in a Jupiter-mass planet as it evolves. We show that a heavy-element core cannot be eroded by convection if there is a sharp compositional change at the core-envelope boundary. If the heavy elements are initially distributed within the planet according to some compositional gradient, mixing occurs in the outer regions resulting in a compositionally homogeneous outer envelope. Mixing of heavy materials that are injected in a convective gaseous envelope are found to mix efficiently. Our work demonstrates that the primordial internal structure of a giant planet plays a substantial role in determining its long-term evolution and that giant planets can have non-adiabatic interiors. These results emphasize the importance of coupling formation, evolution, and internal structure models of giant planets self-consistently.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dust Impact Monitor (SESAME-DIM) Measurements at Comet 67P/Churyumov-Gerasimenko

The Rosetta lander Philae successfully landed on the nucleus of comet 67P/Churyumov-Gerasimenko on 12 November 2014. Philae carries the Dust Impact Monitor (DIM) on board, which is part of the Surface Electric Sounding and Acoustic Monitoring Experiment (SESAME). DIM employs piezoelectric PZT sensors to detect impacts by sub-millimeter and millimeter-sized ice and dust particles that are emitted from the nucleus and transported into the cometary coma. The DIM sensor measures dynamical data like flux and the directionality of the impacting particles. Mass and speed of the particles can be constrained assuming density and elastic particle properties. DIM was operated during three mission phases of Philae at the comet: (1) Before Philae's separation from Rosetta at distances of about 9.6 km, 11.8 km, and 25.3 km from the nucleus barycenter. In this mission phase particles released from the nucleus on radial trajectories remained undetectable because of significant obscuration by the structures of Rosetta, and no dust particles were indeed detected. (2) During Philae's descent to its nominal landing site Agilkia, DIM detected one approximately millimeter-sized particle at a distance of 5.0 km from the nucleus' barycenter, corresponding to an altitude of 2.4 km from the surface. This is the closest ever dust detection at a cometary nucleus by a dedicated in-situ dust detector. (3) At Philae's final landing site, Abydos, DIM detected no dust impact which may be due to low cometary activity in the vicinity of Philae, or due to shading by obstacles close to Philae, or both. Laboratory calibration experiments showed that the material properties of the detected particle are compatible with a porous particle having a bulk density of approximately $250\, \mathrm{kg\,m^{-3}}$. The particle could have been lifted off from the comet's surface by sublimating water ice.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Methane Planets and their Mass-Radius Relation

Knowledge of both the mass and radius of an exoplanet allows us to estimate its mean density, and therefore its composition. Exoplanets seem to fill a very large parameter space in terms of mass and composition, and unlike the solar-system's planets, exoplanets also have intermediate masses (~ 5 - 50 M_Earth) with various densities. In this letter, we investigate the behavior of the Mass-Radius relation for methane (CH_4) planets and show that when methane planets are massive enough (Mp >~ 15 M_Earth), the methane can dissociate and lead to a differentiated planet with a carbon core, a methane envelope, and a hydrogen atmosphere. The contribution of a rocky core to the behavior of CH_4 planet is considered as well. We also develop interior models for several detected intermediate-mass planets that could, in principle, be methane/methane-rich planets. The example of methane planets emphasizes the complexity of the Mass-Radius relation and the challenge involved in uniquely inferring the planetary composition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structure and Dynamics of Cold Water Super-Earths: The Case of Occluded CH4 and its Outgassing

We study the transport of methane in the external water envelopes surrounding water-rich super-Earths and estimate its outgassing into the atmosphere. We investigate the influence of methane on the thermodynamics and mechanics of the water mantle. We find that including methane in the water matrix introduces a new phase (filled ice) resulting in hotter planetary interiors. This effect renders the super-ionic and reticulating phases accessible to relatively low mass planets lacking a H/He atmosphere. We model the thermal and structural profile of the planetary crust and discuss five possible crustal regimes. The formation of methane clathrate in the subsurface is shown to inhibit the formation of a subterranean ocean. This effect results in increased stresses on the lithosphere making modes of ice plate tectonics possible. The dynamics of the tectonic plates are analysed. We derive overturn and resurfacing time scales as well as the melt fraction underneath spreading centers. Ice mantle dynamics is found to be important for assessing the composition of the atmosphere. We formulate the relation between the outgassing flux of methane and the tectonic mode dynamics. We give numerical estimates for the global outgassing rate of methane into the atmosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Volatile Transport inside Super-Earths by Entrapment in the Water Ice Matrix

Whether volatiles can be entrapped in a background matrix composing planetary envelopes and be dragged via convection to the surface is a key question in understanding atmospheric fluxes, cycles and composition. In this paper we consider super-Earths with an extensive water mantle (i.e. water planets), and the possibility of entrapment of methane in their extensive water ice envelopes. We adopt the theory developed by van der Waals & Platteeuw (1959) for modelling solid solutions, often used for modelling clathrate hydrates, and modify it in order to estimate the thermodynamic stability field of a new phase, called methane filled ice Ih. We find that in comparison to water ice VII the filled ice Ih structure may be stable not only at the high pressures but also at the high temperatures expected at the core-water mantle transition boundary of water planets.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmochemical Consequences of Particle Trajectories During FU Orionis Outbursts by the Early Sun

The solar nebula is thought to have undergone a number of episodes of FU Orionis outbursts during its early evolution. We present here the first calculations of the trajectories of particles in a marginally gravitationally unstable solar nebula during an FU Orionis outburst, which show that 0.1 to 10 cm-sized particles traverse radial distances of 10 AU or more, inward and outward, in less than 200 yrs, exposing the particles to temperatures from $\sim$ 60 K to $\sim$ 1500 K. Such trajectories can thus account for the discovery of refractory particles in comets. Refractory particles should acquire Wark-Lovering-like rims as they leave the highest temperature regions of the disk, and these rims should have significant variations in their stable oxygen isotope ratios. Particles are likely to be heavily modified or destroyed if they pass within 1 AU of the Sun, and so are only likely to survive if they formed in the final few FU Orionis outbursts, or were transported to the outer reaches of the solar system. Calcium, aluminum-rich inclusions (CAIs) from primitive meteorites are the oldest known solar system objects and have a very narrow age range. Most CAIs may have formed at the end of the FU Orionis outbursts phase, with an age range reflecting the period between the last few outbursts.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

What do we Really Know about Uranus and Neptune?

The internal structures and compositions of Uranus and Neptune are not well constrained due to the uncertainty in rotation period and flattening, as well as the relatively large error bars on the gravitational coefficients. While Uranus and Neptune are similar in mass and radius, they differ in other physical properties such as thermal emission, obliquity, and inferred atmospheric enrichment. In this letter we consider the uncertainty in the planetary rotation periods, show that rotation periods more consistent with the measured oblateness imply that Uranus and Neptune have different internal structures, and speculate on the source of that difference. We conclude that Uranus and Neptune might have very different structures and/or compositions despite their similar masses and radii. We point out that understanding these differences can have important implications for our view of the formation and evolution of Uranus and Neptune as well as intermediate-mass extra-solar planets in general.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Estimating the Density of Intermediate Size KBOs from Considerations of Volatile Retention

By using a hydrodynamic atmospheric escape mechanism \citep{amit09} we show how the unusually high mass density of Quaoar could have been predicted (constrained), without any knowledge of a binary companion. We suggest an explanation of the recent spectroscopic observations of Orcus and Charon \citep{delsanti10,cook07}. We present a simple relation between the detection of certain volatile ices and the body mass density and diameter. As a test case we implement the relations on the KBO 2003 AZ$_{84}$ and give constraints on its mass density. We also present a method of relating the latitude-dependence of hydrodynamic gas escape to the internal structure of a rapidly rotating body and apply it to Haumea.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 179, 415--431; Lissauer, J. J., Hubickyj, O., D'Angelo, G., Bodenheimer, P., 2009. Icarus 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 Earth masses, for solid surface densities, sigma, of 10, 6, and 4 grams per squared centimeter, respectively. For sigma=10 and sigma=6 g/cm^2, respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interior Models of Uranus and Neptune

'Empirical' models (pressure vs. density) of Uranus and Neptune interiors constrained by the gravitational coefficients J_2, J_4, the planetary radii and masses, and Voyager solid-body rotation periods are presented. The empirical pressure-density profiles are then interpreted in terms of physical equations of state of hydrogen, helium, ice (H_2O), and rock (SiO_2) to test the physical plausibility of the models. The compositions of Uranus and Neptune are found to be similar with somewhat different distributions of the high-Z material. The big difference between the two planets is that Neptune requires a non-solar envelope while Uranus is best matched with a solar composition envelope. Our analysis suggests that the heavier elements in both Uranus' and Neptune's interior might increase gradually towards the planetary centers. Indeed it is possible to fit the gravitational moments without sharp compositional transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Location of the Snow Line in Protostellar Disks

The snow line in a gas disk is defined as the distance from the star beyond which the water ice is stable against evaporation. Since oxygen is the most abundant element after hydrogen and helium, the presence of ice grains can have important consequences for disk evolution. However, determining the position of the snow line is not simple. I discuss some of the important processes that affect the position of the snow line.

0 citations0 reviewsNo code linkedOpen access