Array programming with NumPy.

Array programming provides a powerful, compact and expressive syntax for accessing, manipulating and operating on data in vectors, matrices and higher-dimensional arrays. NumPy is the primary array programming library for the Python language. It has an essential role in research analysis pipelines in fields as diverse as physics, chemistry, astronomy, geoscience, biology, psychology, materials science, engineering, finance and economics. For example, in astronomy, NumPy was an important part of the software stack used in the discovery of gravitational waves1 and in the first imaging of a black hole2. Here we review how a few fundamental array concepts lead to a simple and powerful programming paradigm for organizing, exploring and analysing scientific data. NumPy is the foundation upon which the scientific Python ecosystem is constructed. It is so pervasive that several projects, targeting audiences with specialized needs, have developed their own NumPy-like interfaces and array objects. Owing to its central position in the ecosystem, NumPy increasingly acts as an interoperability layer between such array computation libraries and, together with its application programming interface (API), provides a flexible framework to support the next decade of scientific and industrial analysis.

Pulsar emission amplified and resolved by plasma lensing in an eclipsing binary.

Radio pulsars scintillate because their emission travels through the ionized interstellar medium along multiple paths, which interfere with each other. It has long been realized that, independent of their nature, the regions responsible for the scintillation could be used as 'interstellar lenses' to localize pulsar emission regions1,2. Most such lenses, however, resolve emission components only marginally, limiting results to statistical inferences and detections of small positional shifts3-5. As lenses situated close to their source offer better resolution, it should be easier to resolve emission regions of pulsars located in high-density environments such as supernova remnants 6 or binaries in which the pulsar's companion has an ionized outflow. Here we report observations of extreme plasma lensing in the 'black widow' pulsar, B1957+20, near the phase in its 9.2-hour orbit at which its emission is eclipsed by its companion's outflow7-9. During the lensing events, the observed radio flux is enhanced by factors of up to 70-80 at specific frequencies. The strongest events clearly resolve the emission regions: they affect the narrow main pulse and parts of the wider interpulse differently. We show that the events arise naturally from density fluctuations in the outer regions of the outflow, and we infer a resolution of our lenses that is comparable to the pulsar's radius, about 10 kilometres. Furthermore, the distinct frequency structures imparted by the lensing are reminiscent of what is observed for the repeating fast radio burst FRB 121102, providing observational support for the idea that this source is observed through, and thus at times strongly magnified by, plasma lenses 10 .

A massive pulsar in a compact relativistic binary.

Many physically motivated extensions to general relativity (GR) predict substantial deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 ± 0.04 solar mass (M⊙) pulsar in a 2.46-hour orbit with a 0.172 ± 0.003 M⊙ white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.