Researcher profile

Morris Podolak

Morris Podolak contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 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.

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.

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.

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.

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.