Normalized to: Stewart, S.
[1]
oai:arXiv.org:2006.01992 [pdf] - 2107392
Giant impacts stochastically change the internal pressures of
terrestrial planets
Submitted: 2020-06-02
Pressure is a key parameter in the physics and chemistry of planet formation
and evolution. Previous studies have erroneously assumed that internal
pressures monotonically increase with the mass of a body. Using smoothed
particle hydrodynamics and potential field method calculations, we demonstrate
that the hot, rapidly-rotating bodies produced by giant impacts can have much
lower internal pressures than cool, slowly-rotating planets of the same mass.
Pressures subsequently increase due to thermal and rotational evolution of the
body. Using the Moon-forming impact as an example, we show that the internal
pressures after the collision could have been less than half that in
present-day Earth. The current pressure profile was not established until Earth
cooled and the Moon receded, a process that may take up to 10s Myr after the
last giant impact. Our work defines a new paradigm for pressure evolution
during accretion of terrestrial planets: stochastic changes driven by impacts.
[2]
oai:arXiv.org:2002.00998 [pdf] - 2057908
Silicate Melting and Vaporization during Rocky Planet Formation
Submitted: 2020-02-03
Collisions that induce melting and vaporization can have a substantial effect
on the thermal and geochemical evolution of planets. However, the
thermodynamics of major minerals are not well known at the extreme conditions
attained during planet formation. We obtained new data at the Sandia Z Machine
and use published thermodynamic data for the major mineral forsterite
(Mg$_2$SiO$_4$) to calculate the specific entropy in the liquid region of the
principal Hugoniot. We use our calculated specific entropy of shocked
forsterite, and revised entropies for shocked silica, to determine the critical
impact velocities for melting or vaporization upon decompression from the
shocked state to 1 bar and the triple points, which are near the pressures of
the solar nebula. We also demonstrate the importance of the initial temperature
on the criteria for vaporization. Applying these results to $N$-body
simulations of terrestrial planet formation, we find that up to 20 to 40% of
the total system mass is processed through collisions with velocities that
exceed the criteria for incipient vaporization at the triple point. Vaporizing
collisions between small bodies are an important component of terrestrial
planet formation.
[3]
oai:arXiv.org:1910.04687 [pdf] - 2032844
The Shock Physics of Giant Impacts: Key Requirements for the Equations
of State
Submitted: 2019-10-10, last modified: 2020-01-16
We discuss major challenges in modeling giant impacts between planetary
bodies, focusing on the equations of state (EOS). During the giant impact stage
of planet formation, rocky planets are melted and partially vaporized. However,
most EOS models fail to reproduce experimental constraints on the thermodynamic
properties of the major minerals over the required phase space. Here, we
present an updated version of the widely-used ANEOS model that includes a
user-defined heat capacity limit in the thermal free energy term. Our revised
model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky
planets, provides a better fit to material data over most of the phase space of
giant impacts. We discuss the limitations of this model and the Tillotson
equation of state, a commonly used alternative model.
[4]
oai:arXiv.org:1912.04936 [pdf] - 2042313
The energy budgets of giant impacts
Submitted: 2019-12-10
Giant impacts dominate the final stages of terrestrial planet formation and
set the configuration and compositions of the final system of planets. A giant
impact is believed to be responsible for the formation of Earth's Moon, but the
specific impact parameters are under debate. Because the canonical Moon-forming
impact is the most intensely studied scenario, it is often considered the
archetypal giant impact. However, a wide range of impacts with different
outcomes are possible. Here we examine the total energy budgets of giant
impacts that form Earth-mass bodies and find that they differ substantially
across the wide range of possible Moon-forming events. We show that
gravitational potential energy exchange is important, and we determine the
regime in which potential energy has a significant effect on the collision
outcome. Energy is deposited heterogeneously within the colliding planets,
increasing their internal energies, and portions of each body attain sufficient
entropy for vaporization. After gravitational re-equilibration, post-impact
bodies are strongly thermally stratified, with varying amounts of vaporized and
supercritical mantle. The canonical Moon-forming impact is a relatively low
energy event and should not be considered the archetype of accretionary giant
impacts that form Earth-mass planets. After a giant impact, bodies are
significantly inflated in size compared to condensed planets of the same mass,
and there are substantial differences in the magnitudes of their potential,
kinetic and internal energy components. As a result, the conditions for
metal-silicate equilibration and the subsequent evolution of the planet may
vary widely between different impact scenarios.
[5]
oai:arXiv.org:1910.11293 [pdf] - 2026131
Early Dynamics of the Lunar Core
Submitted: 2019-10-24
The Moon is known to have a small liquid core, and it is thought that in the
distant past the core may have produced strong magnetic fields recorded in
lunar samples. Here we implement a numerical model of lunar orbital and
rotational dynamics that includes the effects of a liquid core. In agreement
with previous work, we find that the lunar core is dynamically decoupled from
the lunar mantle, and that this decoupling happened very early in lunar
history. Our model predicts that the lunar core rotates sub-synchronously, and
the difference between the core and the mantle rotational rates was significant
when the Moon had a high forced obliquity during and after the Cassini State
transition. We find that the presence of the lunar liquid core further
destabilizes synchronous rotation of the mantle for a wide range of semimajor
axes centered around the Cassini State transition. CMB torques make it even
more likely that the Moon experienced large-scale inclination damping during
the Cassini State transition. We present estimates for the mutual core-mantle
obliquity as a function of Earth-Moon distance, and we discuss plausible
absolute time-lines for this evolution. We conclude that our results are
consistent with the hypothesis of a precession-driven early lunar dynamo and
may explain the variability of the inferred orientation of the past lunar
dynamo.
[6]
oai:arXiv.org:1910.00619 [pdf] - 2025919
The energy budget and figure of Earth during recovery from the
Moon-forming giant impact
Submitted: 2019-10-01
Quantifying the energy budget of Earth in the first few million years
following the Moon-forming giant impact is vital to understanding Earth's
initial thermal state and the dynamics of lunar tidal evolution. After the
impact, the body was substantially vaporized and rotating rapidly, very
different from the planet we know today. The subsequent evolution of Earth's
energy budget, as the body cooled and angular momentum was transferred during
lunar tidal recession, has not been accurately calculated with all relevant
energy components included. Here, we use giant impact simulations and planetary
structure models to calculate the energy budget at stages in Earth's evolution.
We show that the figure and internal structure of Earth changed substantially
during its post-impact evolution and that changes in kinetic, potential, and
internal energy were all significant. These changes have important implications
for the dynamics of tidal recession and the thermal structure of early Earth.
[7]
oai:arXiv.org:1802.10223 [pdf] - 1642759
The origin of the Moon within a terrestrial synestia
Submitted: 2018-02-27
The giant impact hypothesis remains the leading theory for lunar origin.
However, current models struggle to explain the Moon's composition and isotopic
similarity with Earth. Here we present a new lunar origin model. High-energy,
high-angular momentum giant impacts can create a post-impact structure that
exceeds the corotation limit (CoRoL), which defines the hottest thermal state
and angular momentum possible for a corotating body. In a typical super-CoRoL
body, traditional definitions of mantle, atmosphere and disk are not
appropriate, and the body forms a new type of planetary structure, named a
synestia. Using simulations of cooling synestias combined with dynamic,
thermodynamic and geochemical calculations, we show that satellite formation
from a synestia can produce the main features of our Moon. We find that cooling
drives mixing of the structure, and condensation generates moonlets that orbit
within the synestia, surrounded by tens of bars of bulk silicate Earth (BSE)
vapor. The moonlets and growing moon are heated by the vapor until the first
major element (Si) begins to vaporize and buffer the temperature. Moonlets
equilibrate with BSE vapor at the temperature of silicate vaporization and the
pressure of the structure, establishing the lunar isotopic composition and
pattern of moderately volatile elements. Eventually, the cooling synestia
recedes within the lunar orbit, terminating the main stage of lunar accretion.
Our model shifts the paradigm for lunar origin from specifying a certain impact
scenario to achieving a Moon-forming synestia. Giant impacts that produce
potential Moon-forming synestias were common at the end of terrestrial planet
formation.
[8]
oai:arXiv.org:1802.03356 [pdf] - 1632542
Tidal evolution of the Moon from a high-obliquity, high-angular-momentum
Earth
Submitted: 2018-02-09
In the giant impact hypothesis for lunar origin, the Moon accreted from an
equatorial circum-terrestrial disk; however the current lunar orbital
inclination of 5 degrees requires a subsequent dynamical process that is still
debated. In addition, the giant impact theory has been challenged by the Moon's
unexpectedly Earth-like isotopic composition. Here, we show that tidal
dissipation due to lunar obliquity was an important effect during the Moon's
tidal evolution, and the past lunar inclination must have been very large,
defying theoretical explanations. We present a new tidal evolution model
starting with the Moon in an equatorial orbit around an initially
fast-spinning, high-obliquity Earth, which is a probable outcome of giant
impacts. Using numerical modeling, we show that the solar perturbations on the
Moon's orbit naturally induce a large lunar inclination and remove angular
momentum from the Earth-Moon system. Our tidal evolution model supports recent
high-angular momentum giant impact scenarios to explain the Moon's isotopic
composition and provides a new pathway to reach Earth's climatically favorable
low obliquity.
[9]
oai:arXiv.org:1712.08234 [pdf] - 1607127
The Elusive Origin of Mercury
Submitted: 2017-12-21
The MESSENGER mission sought to discover what physical processes determined
Mercury's high metal to silicate ratio. Instead, the mission has discovered
multiple anomalous characteristics about our innermost planet. The lack of FeO
and the reduced oxidation state of Mercury's crust and mantle are more extreme
than nearly all other known materials in the solar system. In contrast,
moderately volatile elements are present in abundances comparable to the other
terrestrial planets. No single process during Mercury's formation is able to
explain all of these observations. Here, we review the current ideas for the
origin of Mercury's unique features. Gaps in understanding the innermost
regions of the solar nebula limit testing different hypotheses. Even so, all
proposed models are incomplete and need further development in order to unravel
Mercury's remaining secrets.
[10]
oai:arXiv.org:1712.02790 [pdf] - 1622381
Collisional stripping of planetary crusts
Submitted: 2017-12-07
Geochemical studies of planetary accretion and evolution have invoked various
degrees of collisional erosion to explain differences in bulk composition
between planets and chondrites. Here we undertake a full, dynamical evaluation
of 'crustal stripping' during accretion and its key geochemical consequences.
We present smoothed particle hydrodynamics simulations of collisions between
differentiated rocky planetesimals and planetary embryos. We find that the
crust is preferentially lost relative to the mantle during impacts, and we have
developed a scaling law that approximates the mass of crust that remains in the
largest remnant. Using this scaling law and a recent set of N-body simulations,
we have estimated the maximum effect of crustal stripping on incompatible
element abundances during the accretion of planetary embryos. We find that on
average one third of the initial crust is stripped from embryos as they
accrete, which leads to a reduction of ~20% in the budgets of the heat
producing elements if the stripped crust does not reaccrete. Erosion of crusts
can lead to non-chondritic ratios of incompatible elements, but the magnitude
of this effect depends sensitively on the details of the crust-forming melting
process. The Lu/Hf system is fractionated for a wide range of crustal formation
scenarios. Using eucrites (the products of planetesimal silicate melting,
thought to represent the crust of Vesta) as a guide to the Lu/Hf of
planetesimal crust partially lost during accretion, we predict the Earth could
evolve to a superchondritic 176-Hf/177-Hf (3-5 parts per ten thousand) at
present day. Such values are in keeping with compositional estimates of the
bulk Earth. Stripping of planetary crusts during accretion can lead to
detectable changes in bulk composition of lithophile elements, but the
fractionation is relatively subtle, and sensitive to the efficiency of
reaccretion.
[11]
oai:arXiv.org:1711.06214 [pdf] - 1608655
Col-OSSOS: Colors of the Interstellar Planetesimal 1I/`Oumuamua
Bannister, Michele T.;
Schwamb, Megan E.;
Fraser, Wesley C.;
Marsset, Michael;
Fitzsimmons, Alan;
Benecchi, Susan D.;
Lacerda, Pedro;
Pike, Rosemary E.;
Kavelaars, J. J.;
Smith, Adam B.;
Stewart, Sunny O.;
Wang, Shiang-Yu;
Lehner, Matthew J.
Submitted: 2017-11-16, last modified: 2017-12-06
The recent discovery by Pan-STARRS1 of 1I/2017 U1 (`Oumuamua), on an unbound
and hyperbolic orbit, offers a rare opportunity to explore the planetary
formation processes of other stars, and the effect of the interstellar
environment on a planetesimal surface. 1I/`Oumuamua's close encounter with the
inner Solar System in 2017 October was a unique chance to make observations
matching those used to characterize the small-body populations of our own Solar
System. We present near-simultaneous g$^\prime$, r$^\prime$, and J photometry
and colors of 1I/`Oumuamua from the 8.1-m Frederick C. Gillett Gemini North
Telescope, and $gri$ photometry from the 4.2 m William Herschel Telescope. Our
g$^\prime$r$^\prime$J observations are directly comparable to those from the
high-precision Colours of the Outer Solar System Origins Survey (Col-OSSOS),
which offer unique diagnostic information for distinguishing between outer
Solar System surfaces. The J-band data also provide the highest signal-to-noise
measurements made of 1I/`Oumuamua in the near-infrared. Substantial, correlated
near-infrared and optical variability is present, with the same trend in both
near-infrared and optical. Our observations are consistent with 1I/`Oumuamua
rotating with a double-peaked period of $8.10 \pm 0.42$ hours and being a
highly elongated body with an axial ratio of at least 5.3:1, implying that it
has significant internal cohesion. The color of the first interstellar
planetesimal is at the neutral end of the range of Solar System $g-r$ and $r-J$
solar-reflectance colors: it is like that of some dynamically excited objects
in the Kuiper belt and the less-red Jupiter Trojans.
[12]
oai:arXiv.org:1705.07858 [pdf] - 1583719
The structure of terrestrial bodies: Impact heating, corotation limits
and synestias
Submitted: 2017-05-22
During accretion, terrestrial bodies attain a wide range of thermal and
rotational states, which are accompanied by significant changes in physical
structure (size, shape, pressure and temperature profile, etc.). However,
variations in structure have been neglected in most studies of rocky planet
formation and evolution. Here, we present a new code, HERCULES, that solves for
the equilibrium structure of planets as a series of overlapping
constant-density spheroids. Using HERCULES and a smoothed particle
hydrodynamics code, we show that Earth-like bodies display a dramatic range of
morphologies. For any rotating planetary body, there is a thermal limit beyond
which the rotational velocity at the equator intersects the Keplerian orbital
velocity. Beyond this corotation limit (CoRoL), a hot planetary body forms a
structure, which we name a synestia, with a corotating inner region connected
to a disk-like outer region. By analyzing calculations of giant impacts and
models of planet formation, we show that typical rocky planets are
substantially vaporized multiple times during accretion. For the expected
angular momentum of growing planets, a large fraction of post-impact bodies
will exceed the CoRoL and form synestias. The common occurrence of hot,
rotating states during accretion has major implications for planet formation
and the properties of the final planets. In particular, the structure of
post-impact bodies influences the physical processes that control accretion,
core formation and internal evolution. Synestias also lead to new mechanisms
for satellite formation. Finally, the wide variety of possible structures for
terrestrial bodies also expands the mass-radius range for rocky exoplanets.
[13]
oai:arXiv.org:1511.04830 [pdf] - 1505522
Erosive Hit-and-Run Impact Events: Debris Unbound
Submitted: 2015-11-16, last modified: 2016-10-31
Erosive collisions among planetary embryos in the inner solar system can lead
to multiple remnant bodies, varied in mass, composition and residual velocity.
Some of the smaller, unbound debris may become available to seed the main
asteroid belt. The makeup of these collisionally produced bodies is different
from the canonical chondritic composition, in terms of rock/iron ratio and may
contain further shock-processed material. Having some of the material in the
asteroid belt owe its origin from collisions of larger planetary bodies may
help in explaining some of the diversity and oddities in composition of
different asteroid groups.
[14]
oai:arXiv.org:1605.09788 [pdf] - 1457358
Nominal values for selected solar and planetary quantities: IAU 2015
Resolution B3
Prsa, Andrej;
Harmanec, Petr;
Torres, Guillermo;
Mamajek, Eric;
Asplund, Martin;
Capitaine, Nicole;
Christensen-Dalsgaard, Joergen;
Depagne, Eric;
Haberreiter, Margit;
Hekker, Saskia;
Hilton, James;
Kopp, Greg;
Kostov, Veselin;
Kurtz, Donald W.;
Laskar, Jacques;
Mason, Brian D.;
Milone, Eugene F.;
Montgomery, Michele;
Richards, Mercedes;
Schmutz, Werner;
Schou, Jesper;
Stewart, Susan G.
Submitted: 2016-05-31
In this brief communication we provide the rationale for, and the outcome of
the International Astronomical Union (IAU) resolution vote at the XXIX-th
General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal
conversion constants for selected solar and planetary properties. The problem
addressed by the resolution is a lack of established conversion constants
between solar and planetary values and SI units: a missing standard has caused
a proliferation of solar values (e.g., solar radius, solar irradiance, solar
luminosity, solar effective temperature and solar mass parameter) in the
literature, with cited solar values typically based on best estimates at the
time of paper writing. As precision of observations increases, a set of
consistent values becomes increasingly important. To address this, an IAU
Working Group on Nominal Units for Stellar and Planetary Astronomy formed in
2011, uniting experts from the solar, stellar, planetary, exoplanetary and
fundamental astronomy, as well as from general standards fields to converge on
optimal values for nominal conversion constants. The effort resulted in the IAU
2015 Resolution B3, passed at the IAU General Assembly by a large majority. The
resolution recommends the use of nominal solar and planetary values, which are
by definition exact and are expressed in SI units. These nominal values should
be understood as conversion factors only, not as the true solar/planetary
properties or current best estimates. Authors and journal editors are urged to
join in using the standard values set forth by this resolution in future work
and publications to help minimize further confusion.
[15]
oai:arXiv.org:1510.06262 [pdf] - 1300416
IAU 2015 Resolution B2 on Recommended Zero Points for the Absolute and
Apparent Bolometric Magnitude Scales
Mamajek, E. E.;
Torres, G.;
Prsa, A.;
Harmanec, P.;
Asplund, M.;
Bennett, P. D.;
Capitaine, N.;
Christensen-Dalsgaard, J.;
Depagne, E.;
Folkner, W. M.;
Haberreiter, M.;
Hekker, S.;
Hilton, J. L.;
Kostov, V.;
Kurtz, D. W.;
Laskar, J.;
Mason, B. D.;
Milone, E. F.;
Montgomery, M. M.;
Richards, M. T.;
Schou, J.;
Stewart, S. G.
Submitted: 2015-10-21, last modified: 2015-10-26
The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on
recommended zero points for the absolute and apparent bolometric magnitude
scales. The resolution was proposed by the IAU Inter-Division A-G Working Group
on Nominal Units for Stellar and Planetary Astronomy after consulting with a
broad spectrum of researchers from the astronomical community. Resolution B2
resolves the long-standing absence of an internationally-adopted zero point for
the absolute and apparent bolometric magnitude scales. Resolution B2 defines
the zero point of the absolute bolometric magnitude scale such that a radiation
source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The
zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0)
corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points
were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU
2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the
value most commonly adopted in recent literature. The nominal total solar
irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds
approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832.
Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric
magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based
on the IAU 2012 Resolution B2 definition of the astronomical unit.
[16]
oai:arXiv.org:1510.07674 [pdf] - 1300455
IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for
Selected Solar and Planetary Properties
Mamajek, E. E.;
Prsa, A.;
Torres, G.;
Harmanec, P.;
Asplund, M.;
Bennett, P. D.;
Capitaine, N.;
Christensen-Dalsgaard, J.;
Depagne, E.;
Folkner, W. M.;
Haberreiter, M.;
Hekker, S.;
Hilton, J. L.;
Kostov, V.;
Kurtz, D. W.;
Laskar, J.;
Mason, B. D.;
Milone, E. F.;
Montgomery, M. M.;
Richards, M. T.;
Schou, J.;
Stewart, S. G.
Submitted: 2015-10-26
Astronomers commonly quote the properties of celestial objects in units of
parameters for the Sun, Jupiter, or the Earth. The resolution presented here
was proposed by the IAU Inter-Division Working Group on Nominal Units for
Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly
in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial,
and jovian conversion constants for stellar and (exo)planetary astronomy which
are defined to be exact SI values. While the nominal constants are based on
current best estimates (CBEs; which have uncertainties, are not secularly
constant, and are updated regularly using new observations), they should be
interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts
five solar conversion constants (nominal solar radius, nominal total solar
irradiance, nominal solar luminosity, nominal solar effective temperature, and
nominal solar mass parameter) and six planetary conversion constants (nominal
terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian
equatorial radius, nominal jovian polar radius, nominal terrestrial mass
parameter, and nominal jovian mass parameter).
[17]
oai:arXiv.org:1509.07504 [pdf] - 1304165
Compositional evolution during rocky protoplanet accretion
Submitted: 2015-09-24
The Earth appears non-chondritic in its abundances of refractory lithophile
elements, posing a significant problem for our understanding of its formation
and evolution. It has been suggested that this non-chondritic composition may
be explained by collisional erosion of differentiated planetesimals of
originally chondritic composition. In this work, we present N-body simulations
of terrestrial planet formation that track the growth of planetary embryos from
planetesimals. We simulate evolution through the runaway and oligarchic growth
phases under the Grand Tack model and in the absence of giant planets. These
simulations include a state-of-the-art collision model which allows multiple
collision outcomes, such as accretion, erosion, and bouncing events, that
enables tracking of the evolving core mass fraction of accreting planetesimals.
We show that the embryos grown during this intermediate stage of planet
formation exhibit a range of core mass fractions, and that with significant
dynamical excitation, enough mantle can be stripped from growing embryos to
account for the Earth's non-chondritic Fe/Mg ratio. We also find that there is
a large diversity in the composition of remnant planetesimals, with both
iron-rich and silicate-rich fragments produced via collisions.
[18]
oai:arXiv.org:1410.3421 [pdf] - 1222543
A Collisional Origin to Earth's Non-chondritic Composition?
Submitted: 2014-10-13
Several lines of evidence indicate a non-chondritic composition for Bulk
Earth. If Earth formed from the accretion of chondritic material, its
non-chondritic composition, in particular the super-chondritic 142Nd/144Nd and
low Mg/Fe ratios, might be explained by the collisional erosion of
differentiated planetesimals during its formation. In this work we use an
N-body code, that includes a state-of-the-art collision model, to follow the
formation of protoplanets, similar to proto-Earth, from differentiated
planetesimals (> 100 km) up to isolation mass (> 0.16 M_Earth). Collisions
between differentiated bodies have the potential to change the core-mantle
ratio of the accreted protoplanets. We show that sufficient mantle material can
be stripped from the colliding bodies during runaway and oligarchic growth,
such that the final protoplanets could have Mg/Fe and Si/Fe ratios similar to
that of bulk Earth, but only if Earth is an extreme case and the core is
assumed to contain 10% silicon by mass. This may indicate an important role for
collisional differentiation during the giant impact phase if Earth formed from
chondritic material.
[19]
oai:arXiv.org:1109.4588 [pdf] - 507864
Collisions between Gravity-Dominated Bodies: 2. The Diversity of Impact
Outcomes during the End Stage of Planet Formation
Submitted: 2011-09-21, last modified: 2012-05-02
Numerical simulations of the stochastic end stage of planet formation
typically begin with a population of embryos and planetesimals that grow into
planets by merging. We analyzed the impact parameters of collisions leading to
the growth of terrestrial planets from recent $N$-body simulations that assumed
perfect merging and calculated more realistic outcomes using a new analytic
collision physics model. We find that collision outcomes are diverse and span
all possible regimes: hit-and-run, merging, partial accretion, partial erosion,
and catastrophic disruption. The primary outcomes of giant impacts between
planetary embryos are approximately evenly split between partial accretion,
graze-and-merge, and hit-and-run events. To explore the cumulative effects of
more realistic collision outcomes, we modeled the growth of individual planets
with a Monte Carlo technique using the distribution of impact parameters from
$N$-body simulations. We find that fewer planets reached masses $>0.7 M_{\rm
Earth}$ using the collision physics model compared to simulations that assumed
every collision results in perfect merging. For final planets with masses $>0.7
M_{\rm Earth}$, 60% are enriched in their core-to-mantle mass fraction by >10%
compared to the initial embryo composition. Fragmentation during planet
formation produces significant debris ($\sim15$% of the final mass) and occurs
primarily by erosion of the smaller body in partial accretion and hit-and-run
events. In partial accretion events, the target body grows by preferentially
accreting the iron core of the projectile and the escaping fragments are
derived primarily from the silicate mantles of both bodies. Thus, the bulk
composition of a planet can evolve via stochastic giant impacts.
[20]
oai:arXiv.org:1106.6084 [pdf] - 966951
Collisions Between Gravity-Dominated Bodies: 1. Outcome Regimes and
Scaling Laws
Submitted: 2011-06-29, last modified: 2012-01-05
Collisions are the core agent of planet formation. In this work, we derive an
analytic description of the dynamical outcome for any collision between
gravity-dominated bodies. We conduct high-resolution simulations of collisions
between planetesimals; the results are used to isolate the effects of different
impact parameters on collision outcome. During growth from planetesimals to
planets, collision outcomes span multiple regimes: cratering, merging,
disruption, super-catastrophic disruption, and hit-and-run events. We derive
equations (scaling laws) to demarcate the transition between collision regimes
and to describe the size and velocity distributions of the post-collision
bodies. The scaling laws are used to calculate maps of collision outcomes as a
function of mass ratio, impact angle, and impact velocity, and we discuss the
implications of the probability of each collision regime during planet
formation.
The analytic collision model presented in this work will significantly
improve the physics of collisions in numerical simulations of planet formation
and collisional evolution. (abstract abridged)
[21]
oai:arXiv.org:1011.0387 [pdf] - 1041661
Whole Earth Telescope Observations of the subdwarf B star KPD 1930+2752:
A rich, short period pulsator in a close binary
Reed, M. D.;
Harms, S. L.;
Poindexter, S.;
Zhou, A. -Y.;
Eggen, J. R.;
Morris, M. A.;
Quint, A. C.;
McDaniel, S.;
Baran, A.;
Dolez, N.;
Kawaler, S. D.;
Kurtz, D. W.;
Moskalik, P.;
Riddle, R.;
Zola, S.;
Ostensen, R. H.;
Solheim, J. -E.;
Kepler, S. O.;
Costa, A. F. M.;
Provencal, J. L.;
Mullally, F.;
Winget, D. W.;
Vuckovic, M.;
Crowe, R.;
Terry, D.;
Avila, R.;
Berkey, B.;
Stewart, S.;
Bodnarik, J.;
Bolton, D.;
Binder, P. -M.;
Sekiguchi, K.;
Sullivan, D. J.;
Kim, S. -L.;
Chen, W. -P.;
Chen, C. -W.;
Lin, H. -C.;
Jian, X. -J.;
Wu, H.;
Gou, J. -P.;
Liu, Z.;
Leibowitz, E.;
Lipkin, Y.;
Akan, C.;
Cakirli, O.;
Janulis, R.;
Pretorius, R.;
Ogloza, W.;
Stachowski, G.;
Paparo, M.;
Szabo, R.;
Csubry, Z.;
Zsuffa, D.;
Silvotti, R.;
Marinoni, S.;
Bruni, I.;
Vauclair, G.;
Chevreton, M.;
Matthews, J. M.;
Cameron, C.;
Pablo, H.
Submitted: 2010-11-01, last modified: 2010-11-10
KPD 1930+2752 is a short-period pulsating subdwarf B (sdB) star. It is also
an ellipsoidal variable with a known binary period just over two hours. The
companion is most likely a white dwarf and the total mass of the system is
close to the Chandresakhar limit. In this paper we report the results of Whole
Earth Telescope (WET) photometric observations during 2003 and a smaller
multisite campaign from 2002. From 355 hours of WET data, we detect 68
pulsation frequencies and suggest an additional 13 frequencies within a crowded
and complex temporal spectrum between 3065 and 6343 $\mu$Hz (periods between
326 and 157 s). We examine pulsation properties including phase and amplitude
stability in an attempt to understand the nature of the pulsation mechanism. We
examine a stochastic mechanism by comparing amplitude variations with simulated
stochastic data. We also use the binary nature of KPD 1930+2752 for identifying
pulsation modes via multiplet structure and a tidally-induced pulsation
geometry. Our results indicate a complicated pulsation structure that includes
short-period ($\approx 16$ h) amplitude variability, rotationally split modes,
tidally-induced modes, and some pulsations which are geometrically limited on
the sdB star.
[22]
oai:arXiv.org:1007.3212 [pdf] - 1033764
Water/Icy Super-Earths: Giant Impacts and Maximum Water Content
Submitted: 2010-07-19
Water-rich super-Earth exoplanets are expected to be common. We explore the
effect of late giant impacts on the final bulk abundance of water in such
planets. We present the results from smoothed particle hydrodynamics
simulations of impacts between differentiated water(ice)-rock planets with
masses between 0.5 and 5 M_Earth and projectile to target mass ratios from 1:1
to 1:4. We find that giant impacts between bodies of similar composition never
decrease the bulk density of the target planet. If the commonly assumed maximum
water fraction of 75wt% for bodies forming beyond the snow line is correct,
giant impacts between similar composition bodies cannot serve as a mechanism
for increasing the water fraction. Target planets either accrete materials in
the same proportion, leaving the water fraction unchanged, or lose material
from the water mantle, decreasing the water fraction. The criteria for
catastrophic disruption of water-rock planets are similar to those found in
previous work on super-Earths of terrestrial composition. Changes in bulk
composition for giant impacts onto differentiated bodies of any composition
(water-rock or rock-iron) are described by the same equations. These general
laws can be incorporated into future N-body calculations of planet formation to
track changes in composition from giant impacts.
[23]
oai:arXiv.org:1003.5822 [pdf] - 1026020
The Formation of the Collisional Family around the Dwarf Planet Haumea
Submitted: 2010-03-30
Haumea, a rapidly rotating elongated dwarf planet (~ 1500 km in diameter),
has two satellites and is associated with a "family" of several smaller Kuiper
Belt objects (KBOs) in similar orbits. All members of the Haumea system share a
water ice spectral feature that is distinct from all other KBOs. The relative
velocities between the Haumea family members are too small to have formed by
catastrophic disruption of a large precursor body, which is the process that
formed families around much smaller asteroids in the Main Belt. Here we show
that all of the unusual characteristics of the Haumea system are explained by a
novel type of giant collision: a graze-and-merge impact between two comparably
sized bodies. The grazing encounter imparted the high angular momentum that
spun off fragments from the icy crust of the elongated merged body. The
fragments became satellites and family members. Giant collision outcomes are
extremely sensitive to the impact parameters. Compared to the Main Belt, the
largest bodies in the Kuiper Belt are more massive and experience slower
velocity collisions; hence, outcomes of giant collisions are dramatically
different between the inner and outer solar system. The dwarf planets in the
Kuiper Belt record an unexpectedly large number of giant collisions, requiring
a special dynamical event at the end of solar system formation.
[24]
oai:arXiv.org:1003.0451 [pdf] - 1025464
Minimum Radii of Super-Earths: Constraints from Giant Impacts
Submitted: 2010-03-01
The detailed interior structure models of super-Earth planets show that there
is degeneracy in the possible bulk compositions of a super-Earth at a given
mass and radius, determined via radial velocity and transit measurements,
respectively. In addition, the upper and lower envelopes in the mass--radius
relationship, corresponding to pure ice planets and pure iron planets,
respectively, are not astrophysically well motivated with regard to the
physical processes involved in planet formation. Here we apply the results of
numerical simulations of giant impacts to constrain the lower bound in the
mass--radius diagram that could arise from collisional mantle stripping of
differentiated rocky/iron planets. We provide a very conservative estimate for
the minimum radius boundary for the entire mass range of large terrestrial
planets. This envelope is a readily testable prediction for the population of
planets to be discovered by the Kepler mission.
[25]
oai:arXiv.org:0912.1847 [pdf] - 298957
Constraints on the Source of Lunar Cataclysm Impactors
Submitted: 2009-12-09
Multiple impact basins formed on the Moon about 3.8 Gyr ago in what is known
as the lunar cataclysm or late heavy bombardment. Many workers currently
interpret the lunar cataclysm as an impact spike primarily caused by main-belt
asteroids destabilized by delayed planetary migration. We show that
morphologically fresh (class 1) craters on the lunar highlands were mostly
formed during the brief tail of the cataclysm, as they have absolute crater
number density similar to that of the Orientale basin and ejecta blanket. The
connection between class 1 craters and the cataclysm is supported by the
similarity of their size-frequency distribution to that of
stratigraphically-identified Imbrian craters. Majority of lunar craters younger
than the Imbrium basin (including class 1 craters) thus record the
size-frequency distribution of the lunar cataclysm impactors. This distribution
is much steeper than that of main-belt asteroids. We argue that the projectiles
bombarding the Moon at the time of the cataclysm could not have been main-belt
asteroids ejected by purely gravitational means.
[26]
oai:arXiv.org:0907.0234 [pdf] - 25880
Collisional Stripping and Disruption of Super-Earths
Submitted: 2009-07-02
The final stage of planet formation is dominated by collisions between
planetary embryos. The dynamics of this stage determine the orbital
configuration and the mass and composition of planets in the system. In the
solar system, late giant impacts have been proposed for Mercury, Earth, Mars,
and Pluto. In the case of Mercury, this giant impact may have significantly
altered the bulk composition of the planet. Here we present the results of
smoothed particle hydrodynamics simulations of high-velocity (up to ~5 v_esc)
collisions between 1 and 10 M_Earth planets of initially terrestrial
composition to investigate the end stages of formation of extrasolar
super-Earths. As found in previous simulations of collisions between smaller
bodies, when collision energies exceed simple merging, giant impacts are
divided into two regimes: (1) disruption and (2) hit-and-run (a grazing
inelastic collision and projectile escape). Disruption occurs when the impact
parameter is near zero, when the projectile mass is small compared to the
target, or at extremely high velocities. In the disruption regime, we derive
the criteria for catastrophic disruption (when half the total colliding mass is
lost), the transition energy between accretion and erosion, and a scaling law
for the change in bulk composition (iron-to-silicate ratio) resulting from
collisional stripping of a mantle.
[27]
oai:arXiv.org:0811.0175 [pdf] - 18106
Full Numerical Simulations of Catastrophic Small Body Collisions
Submitted: 2008-11-02
The outcome of collisions between small icy bodies, such as Kuiper belt
objects, is poorly understood and yet a critical component of the evolution of
the trans-Neptunian region. The expected physical properties of outer solar
system materials (high porosity, mixed ice-rock composition, and low material
strength) pose significant computational challenges. We present results from
catastrophic small body collisions using a new hybrid hydrocode to $N$-body
code computational technique. This method allows detailed modeling of shock
propagation and material modification as well as gravitational reaccumulation.
Here, we consider a wide range of material strengths to span the possible range
of Kuiper belt objects. We find that the shear strength of the target is
important in determining the collision outcome for 2 to 50-km radius bodies,
which are traditionally thought to be in a pure gravity regime. The
catastrophic disruption and dispersal criteria, $Q_D^*$, can vary by up to a
factor of three between strong crystalline and weak aggregate materials. The
material within the largest reaccumulated remnants experiences a wide range of
shock pressures. The dispersal and reaccumulation process results in the
material on the surfaces of the largest remnants having experienced a wider
range of shock pressures compared to material in the interior. Hence, depending
on the initial structure and composition, the surface materials on large,
reaccumulated bodies in the outer solar system may exhibit complex spectral and
albedo variations. Finally, we present revised catastrophic disruption criteria
for a range of impact velocities and material strengths for outer solar system
bodies.
[28]
oai:arXiv.org:0810.4010 [pdf] - 17720
The pulsating hot subdwarf Balloon 090100001: results of the 2005
multisite campaign
Baran, A.;
Oreiro, R.;
Pigulski, A.;
Hernandez, F. Perez;
Ulla, A.;
Reed, M. D.;
Rodriguez-Lopez, C.;
Moskalik, P.;
Kim, S. -L.;
Chen, W. -P.;
Crowe, R.;
Siwak, M.;
Armendarez, L.;
Binder, P. M.;
Choo, K. -J.;
Dye, A.;
Eggen, J. R.;
Garrido, R.;
Perez, J. M. Gonzalez;
Harms, S. L.;
Huang, F. -Y.;
Koziel, D.;
Lee, H. -T.;
MacDonald, J.;
Machado, L. Fox;
Monserrat, T.;
Stevick, J.;
Stewart, S.;
Terry, D.;
Zhou, A. -Y.;
Zola, S.
Submitted: 2008-10-22
We present the results of a multisite photometric campaign on the pulsating
sdB star Balloon 090100001. The star is one of the two known hybrid hot
subdwarfs with both long- and short-period oscillations. The campaign involved
eight telescopes with three obtaining UBVR data, four B-band data, and one
Stromgren uvby photometry. The campaign covered 48 nights, providing a temporal
resolution of 0.36microHz with a detection threshold of about 0.2mmag in
B-filter data.
Balloon 090100001 has the richest pulsation spectrum of any known pulsating
subdwarf B star and our analysis detected 114 frequencies including 97
independent and 17 combination ones. The strongest mode (f_1) in the 2.8mHz
region is most likely radial while the remaining ones in this region form two
nearly symmetric multiplets: a triplet and quintuplet, attributed to
rotationally split \ell=1 and 2 modes, respectively. We find clear increases of
splitting in both multiplets between the 2004 and 2005 observing campaigns,
amounting to 15% on average. The observed splittings imply that the rotational
rate in Bal09 depends on stellar latitude and is the fastest on the equator. We
use a small grid of models to constrain the main mode (f_1), which most likely
represents the radial fundamental pulsation. The groups of p-mode frequencies
appear to lie in the vicinity of consecutive radial overtones, up to the third
one. Despite the large number of g-mode frequencies observed, we failed to
identify them, most likely because of the disruption of asymptotic behaviour by
mode trapping. The observed frequencies were not, however, fully exploited in
terms of seismic analysis which should be done in the future with a larger grid
of reliable evolutionary models of hot subdwarfs.
[29]
oai:arXiv.org:0705.3943 [pdf] - 1643
Physical effects of collisions in the Kuiper belt
Submitted: 2007-05-27
Collisions are a major modification process over the history of the Kuiper
Belt. Recent work illuminates the complex array of possible outcomes of
individual collisions onto porous, volatile bodies. The cumulative effects of
such collisions on the surface features, composition, and internal structure of
Kuiper Belt Objects are not yet known. In this chapter, we present the current
state of knowledge of the physics of cratering and disruptive collisions in KBO
analog materials. We summarize the evidence for a rich collisional history in
the Kuiper Belt and present the range possible physical modifications on
individual objects. The question of how well present day bodies represent
primordial planetesimals can be addressed through future studies of the coupled
physical and collisional evolution of Kuiper Belt Objects.
[30]
oai:arXiv.org:astro-ph/0510326 [pdf] - 76822
Asteroseismology of the Beta Cephei star 12 (DD) Lacertae: photometric
observations, pulsational frequency analysis and mode identification
Handler, G.;
Jerzykiewicz, M.;
Rodriguez, E.;
Uytterhoeven, K.;
Amado, P. J.;
Dorokhova, T. N.;
Dorokhov, N. I.;
Poretti, E.;
Sareyan, J. -P.;
Parrao, L.;
Lorenz, D.;
Zsuffa, D.;
Drummond, R.;
Daszynska-Daszkiewicz, J.;
Verhoelst, T.;
De Ridder, J.;
Acke, B.;
Bourge, P. -O.;
Movchan, A. I.;
Garrido, R.;
Paparo, M.;
Sahin, T.;
Antoci, V.;
Udovichenko, S. N.;
Csorba, K.;
Crowe, R.;
Berkey, B.;
Stewart, S.;
Terry, D.;
Mkrtichian, D. E.;
Aerts, C.
Submitted: 2005-10-11
We report a multisite photometric campaign for the Beta Cephei star 12
Lacertae. 750 hours of high-quality differential photoelectric Stromgren,
Johnson and Geneva time-series photometry were obtained with 9 telescopes
during 190 nights. Our frequency analysis results in the detection of 23
sinusoidal signals in the light curves. Eleven of those correspond to
independent pulsation modes, and the remainder are combination frequencies. We
find some slow aperiodic variability such as that seemingly present in several
Beta Cephei stars. We perform mode identification from our colour photometry,
derive the spherical degree l for the five strongest modes unambiguously and
provide constraints on l for the weaker modes. We find a mixture of modes of 0
<= l <= 4. In particular, we prove that the previously suspected rotationally
split triplet within the modes of 12 Lac consists of modes of different l;
their equal frequency splitting must thus be accidental.
One of the periodic signals we detected in the light curves is argued to be a
linearly stable mode excited to visible amplitude by nonlinear mode coupling
via a 2:1 resonance. We also find a low-frequency signal in the light
variations whose physical nature is unclear; it could be a parent or daughter
mode resonantly coupled. The remaining combination frequencies are consistent
with simple light-curve distortions.
The range of excited pulsation frequencies of 12 Lac may be sufficiently
large that it cannot be reproduced by standard models. We suspect that the star
has a larger metal abundance in the pulsational driving zone, a hypothesis also
capable of explaining the presence of Beta Cephei stars in the LMC.
[31]
oai:arXiv.org:astro-ph/0007077 [pdf] - 36915
UV Observations of the Powering Source of the Supergiant Shell in IC2574
Submitted: 2000-07-06
A multi-band analysis of the region containing the supergiant HI shell in the
nearby dwarf irregular galaxy IC2574 presents evidence of a causal relationship
between a central star cluster, the surrounding expanding HI shell, and
secondary star formation sites on the rim of the HI shell. Comparisons of the
far-UV (FUV, 1521 A), optical broad-band, H-alpha, X-ray, and HI morphologies
suggest that the region is in an auspicious moment of star formation triggered
by the central stellar cluster. The derived properties of the HI shell, the
central stellar cluster, and the star forming regions on the rim support this
scenario: The kinematic age of the HI shell is <14 Myr and in agreement with
the age of the central stellar cluster derived from the FUV observations (sim
11 Myr). An estimate for the mechanical energy input from SN and stellar winds
of the central stellar cluster made from FUV photometry and the derived cluster
age is 4.1 x 10^52 erg, roughly a few times higher than the kinetic energy of
the HI shell. The requisite energy input needed to create the HI shell, derived
in the `standard' fashion from the HI observations (using the numerical models
of Chevalier), is 2.6 x 10^53 erg which is almost an order of magnitude higher
than the estimated energy input as derived from the FUV data. Given the
overwhelming observational evidence that the central cluster is responsible for
the expanding HI shell, this discrepancy suggests that the required energy
input is overestimated using the `standard' method. This may explain why some
other searches for remnant stellar clusters in giant HI holes have been
unsuccessful so far. Our observations also show that stellar clusters are
indeed able to create supergiant HI shells, even at large galactocentric radii,
a scenario which has recently been questioned by a number of authors.
[32]
oai:arXiv.org:astro-ph/9910098 [pdf] - 108652
Star Formation Triggering Mechanisms in Dwarf Galaxies: The
Far-Ultraviolet, H-alpha, and HI Morphology of Holmberg II
Stewart, S. G.;
Fanelli, M. N.;
Byrd, G. G.;
Hill, J. K.;
Westpfahl, D. J.;
Cheng, K. P.;
O'Connell, R. W.;
Roberts, M. S.;
Neff, S. G.;
Smith, A. M.;
Stecher, T. P.
Submitted: 1999-10-05
Far-UV (FUV), H-alpha, and HI observations of dwarf galaxy Holmberg II are
used to trace the interaction between sites of massive star formation and the
neutral and ionized components of the surrounding ISM. The data emphasize the
importance of local conditions in regulating star formation from evidence such
as massive stars inside ionized shells, compact HII regions surrounding aging
clusters, and stars formed in chains of progressing age. Large-scale triggering
by HI shells is supported by observations of progenitor populations as well as
secondary sites of star formation associated with their dense rims. Analysis of
the energy available from massive stars inside HI shells indicates energy
deposited into the ISM from supernovae and stellar winds is sufficient to
account for the HI morphology. Ages of individual star forming regions are
derived using B, H-alpha, and FUV photometry. The distribution of ages is
reconciled with the HI morphology, showing a clear preference of young regions
for areas of dense HI and old regions for HI voids. Global kinematical
properties may also play a role in the star formation process since differences
in the rotation characteristics of the neutral gas disk correlate with
differences in triggering mechanisms.