SciPy is an open source scientific computing library for the Python programming language. SciPy 1.0 was released in late 2017, about 16 years after the original version 0.1 release. SciPy has become a de facto standard for leveraging scientific algorithms in the Python programming language, with more than 600 unique code contributors, thousands of dependent packages, over 100,000 dependent repositories, and millions of downloads per year. This includes usage of SciPy in almost half of all machine learning projects on GitHub, and usage by high profile projects including LIGO gravitational wave analysis and creation of the first-ever image of a black hole (M87). The library includes functionality spanning clustering, Fourier transforms, integration, interpolation, file I/O, linear algebra, image processing, orthogonal distance regression, minimization algorithms, signal processing, sparse matrix handling, computational geometry, and statistics. In this work, we provide an overview of the capabilities and development practices of the SciPy library and highlight some recent technical developments.

/pdf/scipy-1-0-fundamental-algorithms-for-scientific-computing-in-111lr790ol.pdf

SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python

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. NumPy is the primary array programming library for Python; here its fundamental concepts are reviewed and its evolution into a flexible interoperability layer between increasingly specialized computational libraries is discussed.

/pdf/array-programming-with-numpy-31s3jl24ra.pdf

Array programming with NumPy

SciPy is an open-source scientific computing library for the Python programming language. Since its initial release in 2001, SciPy has become a de facto standard for leveraging scientific algorithms in Python, with over 600 unique code contributors, thousands of dependent packages, over 100,000 dependent repositories and millions of downloads per year. In this work, we provide an overview of the capabilities and development practices of SciPy 1.0 and highlight some recent technical developments.

/pdf/scipy-1-0-fundamental-algorithms-for-scientific-computing-in-45othhotbu.pdf

SciPy 1.0: fundamental algorithms for scientific computing in Python.

The status of experimental tests of general relativity and of theoretical frameworks for analyzing them is reviewed and updated. Einstein’s equivalence principle (EEP) is well supported by experiments such as the Eotvos experiment, tests of local Lorentz invariance and clock experiments. Ongoing tests of EEP and of the inverse square law are searching for new interactions arising from unification or quantum gravity. Tests of general relativity at the post-Newtonian level have reached high precision, including the light deflection, the Shapiro time delay, the perihelion advance of Mercury, the Nordtvedt effect in lunar motion, and frame-dragging. Gravitational wave damping has been detected in an amount that agrees with general relativity to better than half a percent using the Hulse-Taylor binary pulsar, and a growing family of other binary pulsar systems is yielding new tests, especially of strong-field effects. Current and future tests of relativity will center on strong gravity and gravitational waves.

/pdf/the-confrontation-between-general-relativity-and-experiment-5dqdube5p9.pdf

The Confrontation Between General Relativity and Experiment

Many physically motivated extensions to general relativity (GR) predict significant deviations at energies present in massive neutron stars. We report the measurement of a 2.01 \(\pm \) 0.04 solar mass (M\(_\odot \)) pulsar in a 2.46-h orbit around a 0.172 \(\pm \) 0.003 M\(_\odot \) 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 detection experiments. Additionally, the system strengthens recent constraints on the properties of dense matter and provides novel insight to binary stellar astrophysics and pulsar recycling.

/pdf/a-massive-pulsar-in-a-compact-relativistic-binary-4fslj44n2n.pdf

A massive pulsar in a compact relativistic binary

We present optical and infrared observations taken with the Very Large Telescope of the eclipsing binary pulsar system PSR B1718-19. The candidate companion of the pulsar, identified earlier in Hubble Space Telescope observations, has been detected in all three bands, R , I , and J . These detections allowed us to derive constraints on temperature, radius, and mass, pointing to a companion that has expanded to a radius between one of a main sequence star and one at the Roche-limit. We focus on the role of tidal circularisation in the system, which will have transformed the initially eccentric orbit expected from formation scenarios into the nearly circular orbit presently observed. Based on simple energy balance arguments, we are able to draw a picture of the companion's evolution resulting from the energy deposition in the star due to circularisation. In this picture, our measurement of the companion's parameters is consistent with the expected initial eccentricity. However, with the present understanding of tidal dissipation, it remains difficult to account for the short time in which the system was circularised.

/pdf/observations-of-the-companion-to-the-pulsar-psr-b1718-19-the-4k9b0bc0db.pdf

Observations of the companion to the pulsar PSR B1718-19 - The role of tidal circularisation

Xte j1814-338

I describe recent high-resolution X-ray spectroscopy of surface emission from nearby, thermally emitting neutron stars. I focus on RX J0720.4-3125, RX J1308.6+2127, and RX J1605.3+3249, all of which have similar temperature, but differ in the presence and strength of absorption features in their spectra. I discuss possible causes for the absorption we see in two sources, and conclude that it may be proton cyclotron line absorption, but weakened due to the strong-field quantum electrodynamics effect of vacuum resonance mode conversion.

X-ray Spectroscopy of Thermally Emitting Neutron Stars

A NLTE-analysis is presented of high S/N spectra of the optical component of the standard massive X-ray binary Vela X-1. In combination with the orbital parameters we conclude that the optical star is highly helium enriched and is significantly overluminous compared to standard evolutionary tracks of massive accretion stars. We then propose a new accretion model able to explain these features.

The mass and helium discrepancy in massive stars: the case Vela X-1.

We report the discovery of PSR J1909-3744, a 2.95 millisecond pulsar in a nearly circular 1.53 day orbit. Its narrow pulse width of 43 microseconds allows pulse arrival times to be determined with great accuracy. We have spectroscopically identified the companion as a moderately hot (T = 8500 K) white dwarf with strong absorption lines. Radial velocity measurements of the companion will yield the mass ratio of the system. Our timing data suggest the presence of Shapiro delay; we expect that further timing observations, combined with the mass ratio, will allow the first accurate determination of a millisecond pulsar mass. We have measured the timing parallax and proper motion for this pulsar which indicate a transverse velocity of 140 (+80/-40) km/s. This pulsar's stunningly narrow pulse profile makes it an excellent candidate for precision timing experiments that attempt to detect low frequency gravitational waves from coalescing supermassive black hole binaries.

M. H. van Kerkwijk

Papers

Observations of the companion to the pulsar PSR B1718-19 - The role of tidal circularisation

Xte j1814-338

X-ray Spectroscopy of Thermally Emitting Neutron Stars

The mass and helium discrepancy in massive stars: the case Vela X-1.

PSR J1909-3744, a Binary Millisecond Pulsar with a Very Small Duty Cycle