Swift: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications

Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code Swift. The software package exploits hybrid task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. Swift also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarize the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with Swift.Comment: 39 pages, 18 figures, submitted to MNRAS. Code, documentation, and examples available at www.swiftsim.co

Swift: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications

Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code Swift. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. Swiftalso evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with ≈300 billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with Swift

Measurement of the free neutron lifetime using the neutron spectrometer on NASA's Lunar Prospector mission

We use data from the Lunar Prospector Neutron Spectrometer to make the second space-based measurement of the free neutron lifetime finding τ n = 887 ± 14 stat + 7 − 3 syst s, which is within 1 σ of the accepted value. This measurement expands the range of planetary bodies where the neutron lifetime has been quantified from space, and by extending the modeling to account for nonuniform elemental composition, we mitigated a significant source of systematic uncertainty on the previous space-based lifetime measurement. This modeling moves space-based neutron lifetime measurement towards the ultimate goal of reducing the magnitude of the systematics on a future space-measurement to the level of those seen in laboratory-based experiments

Space-based measurement of the neutron lifetime using data from the neutron spectrometer on NASA's MESSENGER mission

We establish the feasibility of measuring the neutron lifetime via an alternative, space-based class of methods, which use neutrons generated by galactic cosmic ray spallation of planets surfaces and atmospheres. Free neutrons decay via the weak interaction with a mean lifetime of around 880 s. This lifetime constrains the unitarity of the CKM matrix and is a key parameter for studies of Big-Bang nucleosynthesis. However, current laboratory measurements, using two independent approaches, differ by over 4σ. Using data acquired in 2007 and 2008 during flybys of Venus and Mercury by NASA’s MESSENGER spacecraft, which was not designed to make this measurement, we estimate the neutron lifetime to be 780 ± 60stat ± 70syst s, thereby demonstrating the viability of this new approach

Evidence for a Localized Source of the Argon in the Lunar Exosphere

We perform the first tests of various proposed explanations for observed features of the Moon's argon exosphere, including models of the following: spatially varying surface interactions; a source that reflects the lunar near-surface potassium distribution; and temporally varying cold trap areas. Measurements from the Lunar Atmosphere and Dust Environment Explorer (LADEE) and the Lunar Atmosphere Composition Experiment (LACE) are used to test whether these models can reproduce the data. The spatially varying surface interactions hypothesized in previous work cannot reproduce the persistent argon enhancement observed over the western maria. They also fail to match the observed local time of the near-sunrise peak in argon density, which is the same for the highland and mare regions and is well reproduced by simple surface interactions with a ubiquitous desorption energy of 28 kJ mol−1. A localized source can explain the observations, with a trade-off between an unexpectedly localized source or an unexpectedly brief lifetime of argon atoms in the exosphere. To match the observations, a point-like source requires source and loss rates of ∼1.9 × 1021 atoms s−1. A more diffuse source, weighted by the near-surface potassium, requires much higher rates of ∼1.1 × 1022 atoms s−1, corresponding to a mean lifetime of just 1.4 lunar days. We do not address the mechanism for producing a localized source, but demonstrate that this appears to be the only model that can reproduce the observations. Large, seasonally varying cold traps could explain the long-term fluctuation in the global argon density observed by LADEE, but not that by LACE

Abell 1201: detection of an ultramassive black hole in a strong gravitational lens

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