Final water inventories of newly formed terrestrial planets are shaped by their collision history. A setting where volatiles are transported from beyond the snowline to habitable-zone planets suggests collisions of very dry with water-rich bodies. By means of smooth particle hydrodynamics (SPH) simulations we study water delivery in scenarios where a dry target is hit by a water-rich projectile, focusing on hit-and-run encounters with two large surviving bodies, which probably comprise about half of all similar-sized collisions.

Today there are numerous studies on asteroid mining. They elaborate on selecting the right objects, prospecting missions, potential asteroid redirection, and the mining process itself. For economic reasons, most studies focus on mining candidates in the 100-500m size-range. Also, suggestions regarding the design and implementation of space stations or even colonies inside the caverns of mined asteroids exist. Caverns provide the advantages of confined material in near-zero gravity during mining and later the hull will shield the inside from radiation. Existing studies focus on creating the necessary artificial gravity by rotating structures that are built inside the asteroid. Here, we assume the entire mined asteroid to rotate at a sufficient rate for artificial gravity and investigate its use for housing a habitat inside. In this study we present how to estimate the necessary spin rate assuming a cylindrical space station inside a mined asteroid and discuss the implications arising from substantial material stress given the required rotation rate. We estimate the required material strength using two relatively simple analytical models and apply them to fictitious, yet realistic rocky near-Earth asteroids.

The scenario and efficiency of water transport by icy asteroids and comets are still amongst the most important unresolved questions of planetary systems. A better understanding of cometary dynamics in extrasolar systems shall provide information about cometary reservoirs and give an insight into water transport especially to planets in the habitable zone. The detection of Proxima Centauri-b (PCb), which moves in the habitable zone of this system, triggered a debate whether or not this planet can be habitable. In this work, we focus on the stability of an additional planet in the system and on water transport by minor bodies. We perform numerous N-body simulations with PCb and an outer Oort-cloud like reservoir of comets. We investigate close encounters and collisions with the planet, which are important for the transport of water. Observers found hints for a second planet with a period longer than 60 days. Our dynamical studies show that two planets in this system are stable even for a more massive second planet (~12 Earth masses). Furthermore, we perform simulations including exocomets, a second planet, and the influence of the binary Alpha Centauri. The studies on the dynamics of exocomets reveal that the outer limit for water transport is around 200 au. In addition we show that water transport would be possible from a close-in planetesimal cloud (1-4 au). From our simulations, based on typical M-star protoplanetary disks, we estimate the water mass delivered to the planets up to 51 Earth oceans.

In this work we describe a genetic algorithm which is used in order to study orbits of minor bodies in the frames of close encounters. We find that the algorithm in combination with standard orbital numerical integrators can be used as a good proxy for finding typical orbits of minor bodies in close encounters with planets and even their moons, saving a lot of computational time compared to long-term orbital numerical integrations. Here, we study close encounters of Centaurs with Callisto and Ganymede in particular. We also perform n-body numerical simulations for comparison. We find typical impact velocities to be between $v_{rel}=20[v_{esc}]$ and $v_{rel}=30[v_{esc}]$ for Ganymede and between $v_{rel}=25[v_{esc}]$ and $v_{rel}=35[v_{esc}]$ for Callisto.

It has been suggested that the comet-like activity of Main Belt Comets is due to the sublimation of sub-surface water-ice that is exposed when these objects are impacted by meter-sized bodies. We recently examined this scenario and showed that such impacts can in fact excavate ice and present a plausible mechanism for triggering the activation of MBCs (Haghighipour et al. 2016). However, because the purpose of that study was to prove the concept and identify the most viable ice-longevity model, the porosity of the object and the loss of ice due to the heat of impact were ignored. In this paper, we extend our impact simulations to porous materials and account for the loss of ice due to an impact. We show that for a porous MBC, impact craters are deeper, reaching to approximately 15 m implying that if the activation of MBCs is due to the sublimation of sub-surface ice, this ice has to be within the top 15 m of the object. Results also indicate that the loss of ice due to the heat of impact is negligible, and the re-accretion of ejected ice is small. The latter suggests that the activities of current MBCs are most probably from multiple impact sites. Our study also indicates that in order for sublimation from multiple sites to account for the observed activity of the currently known MBCs, the water content of MBCs (and their parent asteroids) needs to be larger than the values traditionally considered in models of terrestrial planet formation.

Collisions between large, similar-sized bodies are believed to shape the final characteristics and composition of terrestrial planets. Their inventories of volatiles such as water, are either delivered or at least significantly modified by such events. Besides the transition from accretion to erosion with increasing impact velocity, similar-sized collisions can also result in hit-and-run outcomes for sufficiently oblique impact angles and large enough projectile-to-target mass ratios. We study volatile transfer and loss focusing on hit-and-run encounters by means of Smooth Particle Hydrodynamics simulations, including all main parameters: impact velocity, impact angle, mass ratio, and also the total colliding mass. We find a broad range of overall water losses, up to 75% in the most energetic hit-and-run events, and confirm the much more severe consequences for the smaller body also for stripping of volatile layers. Transfer of water between projectile and target inventories is found to be mostly rather inefficient, and final water contents are dominated by pre-collision inventories reduced by impact losses, for similar pre-collision water mass fractions. Comparison with our numerical results shows that current collision outcome models are not accurate enough to reliably predict these composition changes in hit-and-run events. To also account for non-mechanical losses we estimate the amount of collisionally vaporized water over a broad range of masses, and find that these contributions are particularly important in collisions of ~Mars-sized bodies, with sufficiently high impact energies, but still relatively low gravity. Our results clearly indicate that the cumulative effect of several (hit-and-run) collisions can efficiently strip protoplanets of their volatile layers, especially the smaller body, as it might be common e.g. for Earth-mass planets in systems with Super-Earths.

We investigate the role of secular and mean motion resonances on the water transport from a belt of icy asteroids onto planets or embryos orbiting inside the circumprimary habitable zone (HZ) of a binary star system. In addition, the host-star has an accompanying gas giant planet. For a comparison, we perform two case studies where a secular resonance (SR) is located either inside the HZ close to 1.0 au (causing eccentric motion of a planet or embryos therein) or in the asteroid belt, beyond the snow line. In the latter case, a higher flux of icy objects moving towards the HZ is expected. Collisions between asteroids and objects in the HZ are treated analytically. Our purely dynamical study shows that the SR in the HZ boosts the water transport; however, collisions can occur at very high impact speeds. In this paper, we treat for the first time, realistic collisions using a GPU 3D-SPH code to assess the water loss in the projectile. Including the water loss into the dynamical results, we get more realistic values for the water mass fraction of the asteroid during an impact. We highlight that collisions occurring at high velocities greatly reduce the water content of the projectile and thus the amount of water transported to planets or embryos orbiting inside the HZ. Moreover, we discuss other effects that could modify our results, namely the asteroid's surface rate recession due to ice sublimation and the atmospheric drag contribution on the asteroids' mass loss.

We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation and failure flows for dynamical behaviour of regolith. In the context of projected small body (asteroid or small moons) sample return missions, we investigate the efficiency and feasibility of a particular material sampling method: Brushes sweep material from the asteroid's surface into a collecting tray. We analyze the influence of different material parameters of regolith such as cohesion and angle of internal friction on the sampling rate. Furthermore, we study the sampling process in two environments by varying the surface gravity (Earth's and Phobos') and we apply different rotation rates for the brushes. We find good agreement of our sampling simulations on Earth with experiments and provide estimations for the influence of the material properties on the collecting rate.

So far, more than 130 extrasolar planets have been found in multiple stellar systems. Dynamical simulations show that the outcome of the planetary formation process can lead to different planetary architectures (i.e. location, size, mass, and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (HZ). In this study, we make a comparison of several binary star systems and aim to show how efficient they are at moving icy asteroids from beyond the snow line into orbits crossing the HZ. We also analyze the influence of secular and mean motion resonances on the water transport towards the HZ. Our study shows that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result in a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star not only favors a faster depletion of our disk of planetesimals, but can also bring 4-5 times more water into the whole HZ. However, due to the secular resonance located either inside the HZ or inside the asteroid belt, impacts between icy planetesimals from the disk and big objects in the HZ can occur at high impact speed. Therefore, real collision modeling using a GPU 3D-SPH code show that in reality, the water content of the projectile is greatly reduced and therefore, also the water transported to planets or embryos initially inside the HZ.

It has been suggested that the comet-like activity of Main Belt Comets are due to the sublimation of sub-surface water-ice that has been exposed as a result of their surfaces being impacted by m-sized bodies. We have examined the viability of this scenario by simulating impacts between m-sized and km-sized objects using a smooth particle hydrodynamics approach. Simulations have been carried out for different values of the impact velocity and impact angle as well as different target material and water-mass fraction. Results indicate that for the range of impact velocities corresponding to those in the asteroid belt, the depth of an impact crater is slightly larger than 10 m suggesting that if the activation of MBCs is due to the sublimation of sub-surface water-ice, this ice has to exist no deeper than a few meters from the surface. Results also show that ice-exposure occurs in the bottom and on the interior surface of impact craters as well as the surface of the target where some of the ejected icy inclusions are re-accreted. While our results demonstrate that the impact scenario is indeed a viable mechanism to expose ice and trigger the activity of MBCs, they also indicate that the activity of the current MBCs is likely due to ice sublimation from multiple impact sites and/or the water contents of these objects (and other asteroids in the outer asteroid belt) is larger than the 5% that is traditionally considered in models of terrestrial planet formation providing more ice for sublimation. We present details of our simulations and discuss their results and implications.

Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request.

The problem of the formation of the Moon is still not explained satisfactorily. While it is a generally accepted scenario that the last giant impact on Earth between some 50 to 100 million years after the starting of the formation of the terrestrial planets formed our natural satellite, there are still many open questions like the isotopic composition which is identical for these two bodies. In our investigation we will not deal with these problems of chemical composition but rather undertake a purely dynamical study to find out the probability of a Mars-sized body to collide with the Earth shortly after the formation of the Earth-like planets. For that we assume an additional massive body between Venus and Earth, respectively Earth and Mars which formed there at the same time as the other terrestrial planets. We have undertaken massive n-body integrations of such a planetary system with 4 inner planets (we excluded Mercury but assumed one additional body as mentioned before) for up to tens of millions of years. Our results led to a statistical estimation of the collision velocities as well as the collision angles which will then serve as the basis of further investigation with detailed SPH computations. We find a most probable origin of the Earth impactor at a semi-major axis of approx. 1.16 AU.

So far, multiple stellar systems harbor more than 130 extra solar planets. Dynamical simulations show that the outcome of planetary formation process can lead to various planetary architecture (i.e. location, size, mass and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (hereafter HZ). In this study, we make a comparison of several binary star systems and their efficiency to move icy asteroids from beyond the snow-line into orbits crossing the HZ. We modeled a belt of 10000 asteroids (remnants from the late phase of planetary formation process) beyond the snow-line. The planetesimals are placed randomly around the primary star and move under the gravitational influence of the two stars and a gas giant. As the planetesimals do not interact with each other, we divided the belt into 100 subrings which were separately integrated. In this statistical study, several double star configurations with a G-type star as primary are investigated. Our results show that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result into a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star can not only favor a faster depletion of our disk of planetesimals but can also bring 4 -- 5 times more water into the whole HZ.

We present preliminary results of terrestrial planet formation using on the one hand classical numerical integration of hundreds of small bodies on CPUs and on the other hand -- for comparison reasons -- the results of our GPU code with thousands of small bodies which then merge to larger ones. To be able to determine the outcome of collision events we use our smooth particle hydrodynamics (SPH) code which tracks how water is lost during such events.

As part of a national scientific network 'Pathways to Habitability' the formation of planets and the delivery of water onto these planets is a key question as water is essential for the development of life. In the first part of the paper we summarize the state of the art of planet formation - which is still under debate in the astronomical community - before we show our results on this topic. The outcome of our numerical simulations depends a lot on the choice of the initial distribution of planetesimals and planetary embryos after gas disappeared in the protoplanetary disk. We also take into account that some of these planetesimals of sizes in the order of the mass of the Moon already contained water; the quantity depends on the distance from the Sun - close-by bodies are dry, but starting from a distance of about 2 AU they can contain substantial amounts of water. We assume that the gas giants and terrestrial planets are already formed when we check the collisions of the small bodies containing water (in the order of a few percent) with the terrestrial planets. We thus are able to give an estimate of the respective contribution to the actual water content (of some Earth-oceans) in the mantle, in the crust and on the surface of Earth. In the second part we discuss in more detail how the formation of larger bodies after a collision may happen as the outcome depends on parameters like collision velocity, impact angle, and the materials involved. We present results obtained by SPH (Smooth Particle Hydrodynamics) simulations. We briefly describe this method and show different scenarios with respect to the formed bodies, possible fragmentation and the water content before and after the collision. In an appendix we discuss detection methods for extrasolar planets (close to 2000 such objects have been discovered so far).

We investigate the outcome of collisions of Ceres-sized planetesimals composed of a rocky core and a shell of water ice. These collisions are not only relevant for explaining the formation of planetary embryos in early planetary systems, but also provide insight into the formation of asteroid families and possible water transport via colliding small bodies. Earlier studies show characteristic collision velocities exceeding the bodies' mutual escape velocity which - along with the distribution of the impact angles - cover the collision outcome regimes 'partial accretion', 'erosion', and 'hit-and-run' leading to different expected fragmentation scenarios. Existing collision simulations use bodies composed of strengthless material; we study the distribution of fragments and their water contents considering the full elasto-plastic continuum mechanics equations also including brittle failure and fragmentation.

Context. We investigate the effects of including material strength in multi-material planetesimal collisions. Aims. The differences between strengthless material models and including the full elasto-plastic model for solid bodies with brittle failure and fragmentation when treating collisions of asteroid-sized bodies as they occur frequently in early planetary systems are demonstrated. Methods. We study impacts of bodies of Ceres-mass with a solid rock impactor and a target with 30 weight-% water content as surface ice.} The initial impact velocities and impact parameters are varied between the escape velocity $v_\mathrm{esc}$ to about 6 $v_\mathrm{esc}$ and from head-on collisions to close fly-bys, respectively. We simulate the collisions using our own SPH code utilizing both strengthless material and the full elasto-plastic material model including brittle failure. Results. One of the most prominent differences is the higher degree of fragmentation and shattered debris clouds in the solid model. In most collision scenarios however, the final outcomes are very similar and differ primarily by the about one order of magnitude higher degree of fragmentation in the solid case. Also, the survivors tend to be of less mass in the solid case which also predicts a higher water loss than the strengthless hydro model. This may be an effect of the relatively low-energy impacts that cannot destroy the solid material instantly. As opposed to giant impacts we also observe an indication that some water ice gets transferred between the bodies.

We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.

We investigate the distribution of encounter velocities and impact angles describing collisions in the habitable zone of the early planetary system. Here we present a catalogue of collision characteristics for a particular mass ratio of the colliding bodies and seven different planetesimal masses ranging from a tenth of Ceres' mass to 10 times the mass of the Moon. We show that there are virtually no collisions with impact speeds lower than the surface escape velocity and a similar velocity-impact angle distribution for different planetesimal masses if velocities are normalized using the escape velocity. An additional perturbing Jupiter-like object distorts the collision velocity and impact picture in the sense that grazing impacts at higher velocities are promoted if the perturber's orbit is close to the habitable zone whereas a more distant perturber has more the effect of a mere widening of the velocity dispersion.

We give a brief introduction to smoothed particle hydrodynamics methods for continuum mechanics. Specifically, we present our 3D SPH code to simulate and analyze collisions of asteroids consisting of two types of material: basaltic rock and ice. We consider effects like brittle failure, fragmentation, and merging in different impact scenarios. After validating our code against previously published results we present first collision results based on measured values for the Weibull flaw distribution parameters of basalt.

This paper presents new and synthetic narrow band photometry of ellipticals, spirals, Seyferts and interacting galaxies in an attempt to identify the cause of the unusually high fraction of blue cluster galaxies in distant clusters (Butcher-Oemler effect). The properties and distribution of the low redshift sample specifically points to starbursts as the origin of the blue narrow band colors in interacting Arp galaxies. Comparison between theoretical models and multicolor diagrams, particularly 4000A break colors, indicates a photometric signature which differs from both normal disk galaxy star formation and nonthermal components. This photometric signature is absent for the Butcher-Oemler galaxies whose general color distribution, compared to present-day clusters, is consistent with a majority of the blue population involved in normal star formation rates (spiral-like) with the addition of a small fraction of bright, blue interacting/merger systems. This photometric picture of the Butcher- Oemler galaxies is in agreement with the morophological evidence from HST imaging.