sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Stewart, Sarah T.

Normalized to: Stewart, S.

32 article(s) in total. 173 co-authors, from 1 to 10 common article(s). Median position in authors list is 3,0.

[1]  oai:arXiv.org:2006.01992  [pdf] - 2107392
Giant impacts stochastically change the internal pressures of terrestrial planets
Comments: Main text: 23 pages, 6 figures. Supplementary material: 10 pages, 8 figures. Published in Science Advances
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
Comments: 15 pages, 8 figures, Accepted for publication in JGR: Planets
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
Comments: M-ANEOS source code public release (https://github.com/isale-code/M-ANEOS); forsterite EOS public release (https://github.com/ststewart/aneos-forsterite-2019); manuscript accepted (no changes): Stewart, S., et al. (accepted). In J. Lane, T. Germann, and M. Armstrong (Eds.), 21st Biennial APS Conference on Shock Compression of Condensed Matter (SCCM19). AIP Publishing
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
Comments: 16 pages, 11 figures. Accepted for publication in JGR: Planets. Supplementary material is available from http://philipjcarter.com/energybudgets/EBSI.pdf . Accompanying animations are available from http://philipjcarter.com/energybudgets/
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
Comments: Accepted for JGR:Planets
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
Comments: 23 pages, 13 figures. Accepted in Earth and Planetary Science Letters
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
Comments: Accepted for publication in Journal of Geophysical Research: Planets. Main text: 44 pages, 24 figures. Supplement: 16 pages, 5 figures, 3 tables
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
Comments: Preprint version of Nature vol. 539, pp-402-406 (2016)
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
Comments: To appear in "Mercury: The View after MESSENGER" edited by Solomon, Nittler & Anderson (www.cambridge.org/9781107154452). This version is free to view and download for personal use only. Not for re-distribution, re-sale or use in derivative works. 37 pages, 5 figures
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
Comments: 15 pages, 9 figures. Accepted for publication in EPSL. Abstract shortened. Accompanying animations can be found at http://www.star.bris.ac.uk/pcarter/crust_strip/
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
Comments: Accepted to ApJL
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
Comments: Main text: 31 pages, 15 figures. Supporting information: 15 pages, 8 figures
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
Comments: 7 pages, 3 figures
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
Comments: 15 pages, 2 tables, accepted to AJ
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
Comments: 4 pages, IAU 2015 Resolution B2, passed by the XXIXth IAU General Assembly in Honolulu, 13 August 2015. All IAU resolutions can be viewed at https://www.iau.org/administration/resolutions/general_assemblies/
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
Comments: 6 pages, IAU 2015 Resolution B3, passed by the XXIXth IAU General Assembly in Honolulu, 13 August 2015. All IAU resolutions can be viewed at https://www.iau.org/administration/resolutions/general_assemblies/
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
Comments: 19 pages, 15 figures. Accepted for publication in ApJ. Accompanying animations can be found at http://www.star.bris.ac.uk/pcarter/comp_evo_15
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?
Comments: 14 pages, 9 figures, accepted for publication in Icarus
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
Comments: Accepted version of manuscript
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
Comments: Version 3, accepted to ApJ in Nov. 2011 published online Dec. 2011. Abstract abridged
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
Comments: 30 pages, 17 figures, 11 tables. Accepted for publication by MNRAS
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
Comments: 9 pages, 4 figures, Accepted for publication in ApJ Letters
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
Comments: Accepted for publication in ApJ, 12 pages, 4 figures, 2 tables.
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
Comments: 8 pages, 4 figures, Accepted for publication in ApJ Letters
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
Comments: Accepted for Icarus
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
Comments: 10 pages, 1 table, 4 figures. Accepted for publication in ApJ Letters
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
Comments: 46 pages, 11 figures, 5 tables, accepted for publication in Icarus. Figures 1 and 8 are low resolution in astro-ph version. Contact Z. M. Leinhardt for full resolution version of paper
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
Comments: accepted for publication in MNRAS
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
Comments: To appear in The Solar System Beyond Neptune, published by Arizona Press
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
Comments: 12 pages, 7 figures, MNRAS, in press
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
Comments: AJ, accepted, 16 pages, 6 figures
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
Comments: 24 Pages, 20 Figures, to appear in ApJ
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.