We perform simulations of giant impacts onto the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100--1000$\times$ improvement in particle number reveals that simulations with below 10^7 particles fail to converge on even bulk properties like the post-impact rotation period, or on the detailed erosion of the atmosphere. Higher resolutions appear to determine these large-scale results reliably, but even 10^8 particles may not be sufficient to study the detailed composition of the debris -- finding that almost an order of magnitude more rock is ejected beyond the Roche radius than with 10^5 particles. We present two software developments that enable this increase in the feasible number of particles. First, we present an algorithm to place any number of particles in a spherical shell such that they all have an SPH density within 1% of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1% of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Second, we develop the hydrodynamics code SWIFT for planetary simulations. SWIFT uses task-based parallelism and other modern algorithmic approaches to take full advantage of contemporary supercomputer architectures. Both the particle placement code and SWIFT are publicly released.

We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors.

We advance the modeling of rubble-pile solid bodies by re-examining the tidal breakup of comet Shoemaker-Levy 9, an event that occurred during a 1.33 Jupiter radii encounter with that planet in July 1992. Tidal disruption of the comet nucleus led to a chain of sub-nuclei about 100-1000 m in diameter; these went on to collide with the planet two years later (Chodas & Yeomans 1996). They were intensively studied prior to and during the collisions, making SL9 the best natural benchmark for physical models of small body disruption. For the first time in the study of this event, we use numerical codes treating rubble-piles as collections of polyhedra (Korycansky & Asphaug 2009). This introduces forces of dilatation and friction, and inelastic response. As in our previous studies (Asphaug & Benz 1994,1996) we conclude that the progenitor must have been a rubble-pile, and we obtain approximately the same pre-breakup diameter (about 1.5 km) in our best fits to the data. We find that the inclusion of realistic fragment shapes leads to grain locking and dilatancy, so that even in the absence of friction or other dissipation we find that disruption is overall more difficult than in our spheres-based simulations. We constrain the comet's bulk density at about 300-400 kg/m^3, half that of our spheres-based predictions and consistent with recent estimates derived from spacecraft observations.

We have investigated the 2009 July impact event on Jupiter using the ZEUS-MP 2 three-dimensional hydrodynamics code. We studied the impact itself and the following plume development. Eight impactors were considered: 0.5 km and 1 km porous (\rho = 1.760 g cm^{-3}) and non-porous (\rho = 2.700 g cm^{-3}) basalt impactors, and 0.5 km and 1 km porous (\rho = 0.600 g cm^{-3}) and non-porous \rho = 0.917 g cm^{-3}) ice impactors. The simulations consisted of these bolides colliding with Jupiter at an incident angle of \theta = 69 degrees from the vertical and with an impact velocity of v = 61.4 km s^{-1}. Our simulations show the development of relatively larger, faster plumes created after impacts involving 1 km diameter bodies. Comparing simulations of the 2009 event with simulations of the Shoemaker-Levy 9 events reveals a difference in plume development, with the higher incident angle of the 2009 impact leading to a shallower terminal depth and a smaller and slower plume. We also studied the amount of dynamical chaos present in the simulations conducted at the 2009 incident angle. Compared to the chaos of the SL9 simulations, where \theta is approximately 45 degrees, we find no significant difference in chaos at the higher 2009 incident angle.

We have studied plume formation after a Jovian comet impact using the ZEUS-MP 2 hydrodynamics code. The three-dimensional models followed objects with 500, 750, and 1000 meter diameters. Our simulations show the development of a fast, upward-moving component of the plume in the wake of the impacting comet that "pinches off" from the bulk of the cometary material ~50 km below the 1 bar pressure level, ~100 km above the depth of greatest mass and energy deposition. The fast-moving component contains about twice the mass of the initial comet, but consists almost entirely (>99.9%) of Jovian atmosphere rather than cometary material. The ejecta rise mainly along the impact trajectory, but an additional vertical velocity component due to buoyancy establishes itself within seconds of impact, leading to an asymmetry in the ejecta with respect to the entry trajectory. The mass of the upward-moving component follows a velocity distribution M(>v) approximately proportional to v^-1.4 (v^-1.6 for the 750 m and 500 m cases) in the velocity range 0.1 < v < 10 km/s.

We have run high-resolution, three-dimensional, hydrodynamic simulations of the impact of comet Shoemaker-Levy 9 into the atmosphere of Jupiter. We find that the energy deposition profile is largely similar to the previous two-dimensional calculations of Mac Low and Zahnle (1994), though perhaps somewhat broader in the range of height over which the energy is deposited. As with similar calculations for impacts into the Venusian atmosphere, there is considerable sensitivity in the results to small changes in the initial conditions, indicating dynamical chaos. We calculated the median depth of energy deposition (the height z at which 50% of the bolide's energy has been released) per run. The mean value among runs is approx 70 km below the 1-bar level, for a 1-km diameter impactor of porous ice of density rho=0.6 g cm^{-3}. The standard deviation among these runs is 14 km. We find little evidence of a trend in these results with the resolution of the calculations (up to 57 cells across the impactor radius, or 8.8-m resolution), suggesting that resolutions as low as 16 grid cells across the radius of the bolide may yield good results for this particular quantity.

The Sun's gradual brightening will seriously compromise the Earth's biosphere within ~ 1E9 years. If Earth's orbit migrates outward, however, the biosphere could remain intact over the entire main-sequence lifetime of the Sun. In this paper, we explore the feasibility of engineering such a migration over a long time period. The basic mechanism uses gravitational assists to (in effect) transfer orbital energy from Jupiter to the Earth, and thereby enlarges the orbital radius of Earth. This transfer is accomplished by a suitable intermediate body, either a Kuiper Belt object or a main belt asteroid. The object first encounters Earth during an inward pass on its initial highly elliptical orbit of large (~ 300 AU) semimajor axis. The encounter transfers energy from the object to the Earth in standard gravity-assist fashion by passing close to the leading limb of the planet. The resulting outbound trajectory of the object must cross the orbit of Jupiter; with proper timing, the outbound object encounters Jupiter and picks up the energy it lost to Earth. With small corrections to the trajectory, or additional planetary encounters (e.g., with Saturn), the object can repeat this process over many encounters. To maintain its present flux of solar energy, the Earth must experience roughly one encounter every 6000 years (for an object mass of 1E22 g). We develop the details of this scheme and discuss its ramifications.