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 R⊕. 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 M⊕ 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 M⊕ impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (∼3 × 1036 kg m2 s−1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 M⊕ impactors than 3 M⊕ ones, because it requires tidal disruption that is suppressed by the more massive impactors.
ABSTRACT
We perform simulations of giant impacts on to the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100–1000 × improvement in particle number ...reveals that simulations with below 107 particles fail to converge on even bulk properties such as 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 108 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 105 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 per cent of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1 per cent of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Secondly, we develop the hydrodynamics code sph with interdependent fine-grained tasking(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 examine the mechanisms by which the atmosphere can be eroded by giant impacts onto Earth-like planets with thin atmospheres, using 3D smoothed particle hydrodynamics simulations with sufficient ...resolution to directly model the fate of low-mass atmospheres. We present a simple scaling law to estimate the fraction lost for any impact angle and speed in this regime. In the canonical Moon-forming impact, only around 10% of the atmosphere would have been lost from the immediate effects of the collision. There is a gradual transition from removing almost none to almost all of the atmosphere for a grazing impact as it becomes more head-on or increases in speed, including complex, nonmonotonic behavior at low impact angles. In contrast, for head-on impacts, a slightly greater speed can suddenly remove much more atmosphere. Our results broadly agree with the application of 1D models of local atmosphere loss to the ground speeds measured directly from our simulations. However, previous analytical models of shock-wave propagation from an idealized point-mass impact significantly underestimate the ground speeds and hence the total erosion. The strong dependence on impact angle and the interplay of multiple nonlinear and asymmetrical loss mechanisms highlight the need for 3D simulations in order to make realistic predictions.
Abstract
We simulate the collision of precursor icy moons analogous to Dione and Rhea as a possible origin for Saturn’s remarkably young rings. Such an event could have been triggered a few hundred ...million years ago by resonant instabilities in a previous satellite system. Using high-resolution smoothed particle hydrodynamics simulations, we find that this kind of impact can produce a wide distribution of massive objects and scatter material throughout the system. This includes the direct placement of pure-ice ejecta onto orbits that enter Saturn’s Roche limit, which could form or rejuvenate rings. In addition, fragments and debris of rock and ice totaling more than the mass of Enceladus can be placed onto highly eccentric orbits that would intersect with any precursor moons orbiting in the vicinity of Mimas, Enceladus, or Tethys. This could prompt further disruption and facilitate a collisional cascade to distribute more debris for potential ring formation, the re-formation of the present-day moons, and evolution into an eventual cratering population of planetocentric impactors.
We present a new scaling law to predict the loss of atmosphere from planetary collisions for any speed, angle, impactor mass, target mass, and body composition, in the regime of giant impacts onto ...broadly terrestrial planets with relatively thin atmospheres. To this end, we examine the erosion caused by a wide range of impacts, using 3D smoothed particle hydrodynamics simulations with sufficiently high resolution to directly model the fate of low-mass atmospheres around 1% of the target's mass. Different collision scenarios lead to extremely different behaviors and consequences for the planets. In spite of this complexity, the fraction of lost atmosphere is fitted well by a power law. Scaling is independent of the system mass for a constant impactor mass ratio. Slow atmosphere-hosting impactors can also deliver a significant mass of atmosphere, but always accompanied by larger proportions of their mantle and core. Different Moon-forming impact hypotheses suggest that around 10%-60% of a primordial atmosphere could have been removed directly, depending on the scenario. We find no evident departure from the scaling trends at the extremes of the parameters explored. The scaling law can be incorporated readily into models of planet formation.
Abstract
The Moon is traditionally thought to have coalesced from the debris ejected by a giant impact onto the early Earth. However, such models struggle to explain the similar isotopic compositions ...of Earth and lunar rocks at the same time as the system’s angular momentum, and the details of potential impact scenarios are hotly debated. Above a high resolution threshold for simulations, we find that giant impacts can immediately place a satellite with similar mass and iron content to the Moon into orbit far outside Earth’s Roche limit. Even satellites that initially pass within the Roche limit can reliably and predictably survive, by being partially stripped and then torqued onto wider, stable orbits. Furthermore, the outer layers of these directly formed satellites are molten over cooler interiors and are composed of around 60% proto-Earth material. This could alleviate the tension between the Moon’s Earth-like isotopic composition and the different signature expected for the impactor. Immediate formation opens up new options for the Moon’s early orbit and evolution, including the possibility of a highly tilted orbit to explain the lunar inclination, and offers a simpler, single-stage scenario for the origin of the Moon.
ABSTRACT
We simulate the hypothesized collision between the proto-Earth and a Mars-sized impactor that created the Moon. Among the resulting debris disc in some impacts, we find a self-gravitating ...clump of material. It is roughly the mass of the Moon, contains $\sim 1{{\ \rm per\ cent}}$ iron like the Moon, and has its internal composition resolved for the first time. The clump contains mainly impactor material near its core but becomes increasingly enriched in proto-Earth material near its surface. The formation of this Moon-sized clump depends sensitively on the spin of the impactor. To explore this, we develop a fast method to construct models of multilayered, rotating bodies and their conversion into initial conditions for smoothed particle hydrodynamical (SPH) simulations. We use our publicly available code to calculate density and pressure profiles in hydrostatic equilibrium and then generate configurations of over a billion particles with SPH densities within 1 per cent of the desired values. This algorithm runs in a few minutes on a desktop computer, for 107 particles, and allows direct control over the properties of the spinning body. In comparison, alternative relaxation or spin-up techniques take hours on a supercomputer and the structure of the rotating body cannot be known beforehand. Collisions that differ only in the impactor’s initial spin reveal a wide variety of outcomes: a merger, a grazing hit-and-run, or the creation of an orbiting proto-Moon.
ABSTRACT
Density discontinuities cannot be precisely modelled in standard formulations of smoothed particles hydrodynamics (SPH) because the density field is defined smoothly as a kernel-weighted sum ...of neighbouring particle masses. This is a problem when performing simulations of giant impacts between proto-planets, for example, because planets typically do have density discontinuities both at their surfaces and at any internal boundaries between different materials. The inappropriate densities in these regions create artificial forces that effectively suppress mixing between particles of different material and, as a consequence, this problem introduces a key unknown systematic error into studies that rely on SPH simulations. In this work, we present a novel, computationally cheap method that deals simultaneously with both of these types of density discontinuity in SPH simulations. We perform standard hydrodynamical tests and several example giant impact simulations, and compare the results with standard SPH. In a simulated Moon-forming impact using 107 particles, the improved treatment at boundaries affects at least 30${{\ \rm per\ cent}}$ of the particles at some point during the simulation.
ABSTRACT
Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass MBH and host galaxy velocity dispersion σe. ...Outside the local Universe, measurements of MBH are usually only possible for SMBHs in an active state: limiting sample size and introducing selection biases. Gravitational lensing makes it possible to measure the mass of non-active SMBHs. We present models of the $z$ = 0.169 galaxy-scale strong lens Abell 1201. A cD galaxy in a galaxy cluster, it has sufficient ‘external shear’ that a magnified image of a $z$ = 0.451 background galaxy is projected just ∼1 kpc from the galaxy centre. Using multiband Hubble Space Telescope imaging and the lens modelling software PYAUTOLENS, we reconstruct the distribution of mass along this line of sight. Bayesian model comparison favours a point mass with MBH = 3.27 ± 2.12 × 1010 M⊙ (3σ confidence limit); an ultramassive black hole. One model gives a comparable Bayesian evidence without an SMBH; however, we argue this model is nonphysical given its base assumptions. This model still provides an upper limit of MBH ≤ 5.3 × 1010 M⊙, because an SMBH above this mass deforms the lensed image ∼1 kpc from Abell 1201’s centre. This builds on previous work using central images to place upper limits on MBH, but is the first to also place a lower limit and without a central image being observed. The success of this method suggests that surveys during the next decade could measure thousands more SMBH masses, and any redshift evolution of the MBH−σe relation. Results are available at https://github.com/Jammy2211/autolens_abell_1201.
Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by the presence of low-density granulocytes (LDGs) with a heightened capacity for spontaneous NETosis, but the contribution ...of LDGs to SLE pathogenesis remains unclear. To characterize LDGs in human SLE, gene expression profiles derived from isolated LDGs were characterized by weighted gene coexpression network analysis, and a 92-gene module was identified. The LDG gene signature was enriched in genes related to neutrophil degranulation and cell cycle regulation. This signature was assessed in gene expression datasets from two large-scale SLE clinical trials to study associations between LDG enrichment, SLE manifestations, and treatment regimens. LDG enrichment in the blood was associated with corticosteroid treatment as well as anti-dsDNA, low serum complement, renal manifestations, and vasculitis, but the latter two of these associations were dependent on concomitant corticosteroid treatment. In addition, LDG enrichment was associated with enrichment of gene signatures induced by type I IFN and TNF irrespective of corticosteroid treatment. Notably, LDG enrichment was not found in numerous tissues affected by SLE. Comparison with relevant reference datasets indicated that LDG enrichment is likely reflective of increased granulopoiesis in the bone marrow and not peripheral neutrophil activation. The results have uncovered important determinants of the appearance of LDGs in SLE and have emphasized the likely role of LDGs in specific aspects of lupus pathogenesis.