Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

We present results from the FUNK experiment in the search for hidden-photon dark matter. Near the surface of a mirror, hidden photons may be converted into ordinary photons. These photons are emitted perpendicular to the surface and have an energy equal to the mass of the dark matter hidden photon. Our experimental setup consists of a large, spherical mirror with an area of more than 14 m$^2$, which concentrates the emitted photons into its central point. Using a detector sensitive to visible and near-UV photons, we can exclude a kinetic-mixing coupling of stronger than $\chi \approx 10^{-12}$ in the mass range of 2.5 to 7 eV, assuming hidden photons comprise all of the dark matter. The experimental setup and analysis used to obtain this limit are discussed in detail.

We present strong-lensing analyses of three galaxy clusters, RXJ2129.4+0009 (z=0.235), MS0451.6-0305 (z=0.55), and MACSJ2129.4-0741 (z=0.589), using the powerful combination of Hubble Space Telescope (HST) multi-band observations, and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. We confirm the spectroscopic identification of 4 new multiple image systems, and 36 cluster members in RXJ2129. Our resulting mass model is constrained by 8 multiple image systems, and composed of 71 halos including one dark matter cluster-scale halo and 2 galaxy-scale halos optimized individually. For MS0451, we report the spectroscopic identification of 2 new systems of multiple images in the poorly-constrained Northern region, and 44 cluster members. Our mass model uses 16 multiple image systems, and 146 halos, including 2 large-scale halos, and 7 galaxy-scale halos independently optimized. For MACSJ2129, we confirm the spectroscopic identification of 4 new systems of multiple images, and 52 cluster members. We also report the spectroscopic identification of one new multiple image system at $z=4.41$. Our mass model is optimized by 14 multiple image systems, and composed of 151 halos, including 2 large-scale halos and 4 galaxy-scale halos independently optimized. Our best models have rms of 0.33",0.6", 0.73" in the image plane for RXJ2129, MS0451, and MACSJ2129 respectively. This analysis presents a detailed comparison with the existing literature showing excellent agreements, and discuss specific studies of lensed galaxies, e.g. a group of submilimeter galaxies at z=2.9 in MS0451, and a bright z=2.1472 red singly imaged galaxy in MACSJ2129.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

With the upcoming generation of telescopes, cluster scale strong gravitational lenses will act as an increasingly relevant probe of cosmology and dark matter. The better resolved data produced by current and future facilities requires faster and more efficient lens modeling software. Consequently, we present Lenstool-HPC, a strong gravitational lens modeling and map generation tool based on High Performance Computing (HPC) techniques and the renowned Lenstool software. We also showcase the HPC concepts needed for astronomers to increase computation speed through massively parallel execution on supercomputers. Lenstool-HPC was developed using lens modelling algorithms with high amounts of parallelism. Each algorithm was implemented as a highly optimised CPU, GPU and Hybrid CPU-GPU version. The software was deployed and tested on the Piz Daint cluster of the Swiss National Supercomputing Centre (CSCS). Lenstool-HPC perfectly parallel lens map generation and derivative computation achieves a factor 30 speed-up using only 1 GPUs compared to Lenstool. Lenstool-HPC hybrid Lens-model fit generation tested at Hubble Space Telescope precision is scalable up to 200 CPU-GPU nodes and is faster than Lenstool using only 4 CPU-GPU nodes.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

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.

According to current evidence the water inventory of Earth (and perhaps similar exoplanets) was transported inwards via (giant) collisions during the chaotic final phase of planet formation. In dynamical simulations water delivery is still studied almost exclusively by assuming oversimplified perfect merging (PM), even though it is particularly prone to collisional transfer and loss. To close this gap we have developed a framework to model collisional water transport by direct combination of long-term N-body computations with dedicated 3D SPH simulations for each collision. Post-collision water inventories are self-consistently traced further, in accretionary or erosive as well as hit-and-run encounters. The latter are frequent outcomes among protoplanets, where besides collisional losses, water transfer between the encountering bodies has to be considered. This hybrid approach enables us for the first time to trace the full dynamical and collisional evolution of ~200 bodies throughout the whole late-stage accretion phase (several 100 Myrs). As a first application we choose a Solar System-like architecture, with already formed giant planets on either circular or eccentric orbits and a debris disk spanning from 0.5 - 4 au. Realistic collision treatment leads to considerably different results than PM, with lower mass planets and water inventories reduced by a factor of two or more. Due to a combination of collisional losses and considerably lengthened accretion, final water contents especially with giant planets on circular orbits are strongly reduced to more Earth-like values. Water delivery to potentially habitable planets is dominated by few decisive collisions, mostly with embryo-sized or larger bodies. The high frequency of hit-and-run, with generally low water (and mass) transfer efficiencies, are a crucial part of this process, and of system-wide evolution in general.

Classification of intermediate redshift ($z$ = 0.3--0.8) emission line galaxies as star-forming galaxies, composite galaxies, active galactic nuclei (AGN), or low-ionization nuclear emission regions (LINERs) using optical spectra alone was impossible because the lines used for standard optical diagnostic diagrams: [NII], H$\alpha$, and [SII] are redshifted out of the observed wavelength range. In this work, we address this problem using four supervised machine learning classification algorithms: $k$-nearest neighbors (KNN), support vector classifier (SVC), random forest (RF), and a multi-layer perceptron (MLP) neural network. For input features, we use properties that can be measured from optical galaxy spectra out to $z < 0.8$---[OIII]/H$\beta$, [OII]/H$\beta$, [OIII] line width, and stellar velocity dispersion---and four colors ($u-g$, $g-r$, $r-i$, and $i-z$) corrected to $z=0.1$. The labels for the low redshift emission line galaxy training set are determined using standard optical diagnostic diagrams. RF has the best area under curve (AUC) score for classifying all four galaxy types, meaning highest distinguishing power. Both the AUC scores and accuracies of the other algorithms are ordered as MLP$>$SVC$>$KNN. The classification accuracies with all eight features (and the four spectroscopically-determined features only) are 93.4% (92.3%) for star-forming galaxies, 69.4% (63.7%) for composite galaxies, 71.8% (67.3%) for AGNs, and 65.7% (60.8%) for LINERs. The stacked spectrum of galaxies of the same type as determined by optical diagnostic diagrams at low redshift and RF at intermediate redshift are broadly consistent. Our publicly available code (https://github.com/zkdtc/MLC_ELGs) and trained models will be instrumental for classifying emission line galaxies in upcoming wide-field spectroscopic surveys.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

Strong gravitational lensing is a powerful probe of cosmology and the dark matter distribution. Efficient lensing software is already a necessity to fully use its potential and the performance demands will only increase with the upcoming generation of telescopes. In this paper, we study the possible impact of High Performance Computing techniques on a performance-critical part of the widely used lens modeling software LENSTOOL. We implement the algorithm once as a highly optimized CPU version and once with graphics card acceleration for a simple parametric lens model. In addition, we study the impact of finite machine precision on the lensing algorithm. While double precision is the default choice for scientific applications, we find that single precision can be sufficiently accurate for our purposes and lead to a big speedup. Therefore we develop and present a mixed precision algorithm which only uses double precision when necessary. We measure the performance of the different implementations and find that the use of High Performance Computing Techniques dramatically improves the code performance both on CPUs and GPUs. Compared to the current LENSTOOL implementation on 12 CPU cores, we obtain speedup factors of up to 170. We achieve this optimal performance by using our mixed precision algorithm on a high-end GPU which is common in modern supercomputers. We also show that these techniques reduce the energy consumption by up to 98%. Furthermore, we demonstrate that a highly competitive speedup can be reached with consumer GPUs. While they are an order of magnitude cheaper than the high-end graphics cards, they are rarely used for scientific computations due to their low double precision performance. Our mixed precision algorithm unlocks their full potential. The consumer GPU delivers a speedup which is only a factor of four lower than the best speedup achieved by a high-end GPU.

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.

We investigate the outcome of collisions in very different mass regimes, but an otherwise identical parameter setup, comprising the impact velocity ($v/v_\mathrm{esc}$), impact angle, mass ratio, and initial composition, w.r.t. simple hydrodynamic scaling. The colliding bodies' masses range from $\simeq 10^{16}$ to $10^{24}$ kg, which includes km-sized planetesimals up to planetary-sized objects. Our analysis of the results comprises the time evolution of fragment masses, the fragments' water contents and fragment dynamics, where we start with bodies consisting of basalt and water ice. The usual assumption of hydrodynamic scaling over a wider range of masses is based on material behavior similar to a fluid, or a rubble pile, respectively. All our simulations are carried out once including full solid-body physics, and once for strengthless - but otherwise identical - bodies, to test for the influence of material strength. We find that scale-invariance over a wider range of masses is mostly only a very crude approximation at best, but can be applied to constrained mass ranges if tested carefully. For the chosen scenarios the outcomes of solid-body objects compared to strengthless fluid bodies differ most for our intermediate masses, but are similar for the lowest and highest masses. The most energetic, planet-sized collisions produce considerably faster and more fragments, which is also reflected in high water losses - above 50 percent in a single collision.

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ~4 lenses per square degree down to i=24.8. This is equivalent to ~ 60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

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.

If dark matter consists of hidden-sector photons which kinetically mix with regular photons, a tiny oscillating electric-field component is present wherever we have dark matter. In the surface of conducting materials this induces a small probability to emit single photons almost perpendicular to the surface, with the corresponding photon frequency matching the mass of the hidden photons. We report on a construction of an experimental setup with a large ~14 m2 spherical metallic mirror that will allow for searches of hidden-photon dark matter in the eV and sub-eV range by application of different electromagnetic radiation detectors. We discuss sensitivity and accessible regions in the dark matter parameter space.

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

Following the recent insight in the material structure of comets, protoplanetesimals are assumed to have low densities and to be highly porous agglomerates. It is still unclear if planetesimals can be formed from these objects by collisional growth. Therefore, it is important to study numerically the collisional outcome from low velocity impacts of equal sized porous agglomerates which are too large to be examined in a laboratory experiment. We use the Lagrangian particle method Smooth Particle Hydrodynamics to solve the equations that describe the dynamics of elastic and plastic bodies. Additionally, to account for the influence of porosity, we follow a previous developed equation of state and certain relations between the material strength and the relative density. Collisional growth seems possible for rather low collision velocities and particular material strengths. The remnants of collisions with impact parameters that are larger than 50% of the radius of the colliding objects tend to rotate. For small impact parameters, the colliding objects are effectively slowed down without a prominent compaction of the porous structure, which probably increases the possibility for growth. The protoplanetesimals, however, do not stick together for the most part of the employed material strengths. An important issue in subsequent studies has to be the influence of rotation to collisional growth. Moreover, for realistic simulations of protoplanetesimals it is crucial to know the correct material parameters in more detail.