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

I review what properties of the white-dwarf companions of recycled pulsars can be inferred from optical observations, and discuss how these can help us understand the characteristics and evolution of these binaries. I focus on spectroscopic observations, describing results obtained recently, and looking forward to what may come.

The White Dwarf Companions of Recycled Pulsars

We measured the radial-velocity curve of HD77581, the B-supergiant companion of the X-ray pulsar Vela X-1, using 183 high-resolution optical spectra obtained in a nine-month campaign. We derive radial-velocity amplitudes for different lines and wavelength regions, and find all are consistent with each other, as well as with values found in previous analyses. We show that one apparent exception, an anomalously low value derived from ultra-violet spectra obtained with the IUE, was due to an error in the analysis procedures. We re-analyse all IUE spectra, and combine the resulting velocities with the ones derived from the new optical spectra presented here, as well as those derived from optical spectra published earlier. As in previous analyses, the radial velocities show strong deviations from those expected for a pure Keplerian orbit. The deviations likely are related to the pronounced line-profile variations seen in our spectra. It turns out that systematic deviations as a function of orbital phase are present. Our best estimate of the radial-velocity amplitude, Kopt = 21.7 +/- 1.6 km/s, has an uncertainty not much reduced to that found in previous analyses, in which the systematic deviations had not been taken into account. Combining our velocity amplitude with the accurate orbital elements of the X-ray pulsar, we infer M_ns sin^3i = 1.78 +/- 0.15 Msun.

/pdf/the-mass-of-the-neutron-star-in-vela-x-1-4pxh6i76r5.pdf

The mass of the neutron star in Vela X-1

We present K band photometry of the region including the radio error box of the soft X-ray transient 4U 1630-47 during its 1998 outburst. We detect a variable source at K=16.1 mag located inside the radio error circle which we identify as the counterpart to the X-ray source. We discuss the properties of the source, and conclude that it is most likely a black-hole X-ray binary similar to 4U 1543-47, GRO J1655-40 or SAX J1819.3-2525, containing a relatively early-type secondary.

/pdf/the-ir-counterpart-of-the-black-hole-candidate-4u-1630-47-sl4as5wawc.pdf

The IR counterpart of the black-hole candidate 4U 1630-47

We present phase-resolved spectroscopy of the op- tical counterparts V 801 Ara and V 926 Sco of the low-mass X-ray binaries 1636-536 and 1735-444. To assist in the in- terpretation of the spectroscopic observations we derive new ephemerides for the photometric variations of both sources. Su- perior conjunction of the line emission region occurs at photo- metric phase 0.7 in both sources, which indicates that these variations are dominated by a component originating from the point where the mass stream from the donor star intersects the outer disk. We nd that the properties of 1636-536/V 801 Ara and 1735-444/V 926 Sco are very similar.

/pdf/phase-resolved-spectroscopy-of-the-low-mass-x-ray-binaries-50lo03c0ho.pdf

Phase-resolved spectroscopy of the low-mass X-ray binaries 1636-536/V801 Arae and 1735-444/V926 Scorpii

/pdf/hst-observations-of-psr-j2051-0827-asymmetry-variability-and-geruxbh2bc.pdf

M. H. van Kerkwijk

Papers

The White Dwarf Companions of Recycled Pulsars

The mass of the neutron star in Vela X-1

The IR counterpart of the black-hole candidate 4U 1630-47

Phase-resolved spectroscopy of the low-mass X-ray binaries 1636-536/V801 Arae and 1735-444/V926 Scorpii

HST Observations of PSR J2051-0827: Asymmetry, Variability and Modelling