The giant impact hypothesis for Moon formation successfully explains the dynamic properties of the Earth-Moon system but remains challenged by the similarity of isotopic fingerprints of the terrestrial and lunar mantles. Moreover, recent geochemical evidence suggests that the Earth's mantle preserves ancient (or "primordial") heterogeneity that predates the Moon-forming giant impact. Using a new hydrodynamical method, we here show that Moon-forming giant impacts lead to a stratified starting condition for the evolution of the terrestrial mantle. The upper layer of the Earth is compositionally similar to the disk, out of which the Moon evolves, whereas the lower layer preserves proto-Earth characteristics. As long as this predicted compositional stratification can at least partially be preserved over the subsequent billions of years of Earth mantle convection, the compositional similarity between the Moon and the accessible Earth's mantle is a natural outcome of realistic and high-probability Moon-forming impact scenarios. The preservation of primordial heterogeneity in the modern Earth not only reconciles geochemical constraints but is also consistent with recent geophysical observations. Furthermore, for significant preservation of a proto-Earth reservoir, the bulk composition of the Earth-Moon system may be systematically shifted towards chondritic values.

Giant impacts (GIs) are common in the late stage of planet formation. The Smoothed Particle Hydrodynamics (SPH) method is widely used for simulating the outcome of such violent collisions, one prominent example being the formation of the Moon. However, a decade of numerical studies in various areas of computational astrophysics has shown that the standard formulation of SPH suffers from several shortcomings such as artificial surface tension and its tendency to promptly damp turbulent motions on scales much larger than the physical dissipation scale, both resulting in the suppression of mixing. In order to quantify how severe these limitations are when modeling GIs we carried out a comparison of simulations with identical initial conditions performed with the standard SPH as well as with the novel Lagrangian Meshless Finite Mass (MFM) method in the GIZMO code. We confirm the lack of mixing between the impactor and target when SPH is employed, while MFM is capable of driving vigorous sub-sonic turbulence and leads to significant mixing between the two bodies. Modern SPH variants with artificial conductivity, a different formulation of the hydro force or reduced artificial viscosity, do not improve mixing as significantly. Angular momentum is conserved similarly well in both methods, but MFM does not suffer from spurious transport induced by artificial viscosity, resulting in a slightly higher angular momentum of the proto-lunar disk. Furthermore, SPH initial conditions exhibit an unphysical density discontinuity at the core-mantle boundary which is easily removed in MFM.

There is a growing population of relativistically relevant minor bodies in the Solar System and a growing population of massive extrasolar planets with orbits very close to the central star where relativistic effects should have some signature. Our purpose is to review how general relativity affects the orbital dynamics of the planetary systems and to define a suitable relativistic correction for Solar System orbital studies when only point masses are considered. Using relativistic formulae for the N body problem suited for a planetary system given in the literature we present a series of numerical orbital integrations designed to test the relevance of the effects due to the general theory of relativity in the case of our Solar System. Comparison between different algorithms for accounting for the relativistic corrections are performed. Relativistic effects generated by the Sun or by the central star are the most relevant ones and produce evident modifications in the secular dynamics of the inner Solar System. The Kozai mechanism, for example, is modified due to the relativistic effects on the argument of the perihelion. Relativistic effects generated by planets instead are of very low relevance but detectable in numerical simulations.