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

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.

We present cosmological parameter results from the ﬁnal full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to signiﬁcant gains in the precision of other correlated parameters Improved modelling of the small-scale polarization leads to more robust constraints on manyparameters,withresidualmodellinguncertaintiesestimatedtoaﬀectthemonlyatthe05σlevelWeﬁndgoodconsistencywiththestandard spatially-ﬂat6-parameter ΛCDMcosmologyhavingapower-lawspectrumofadiabaticscalarperturbations(denoted“base ΛCDM”inthispaper), from polarization, temperature, and lensing, separately and in combination A combined analysis gives dark matter density Ωch2 = 0120±0001, baryon density Ωbh2 = 00224±00001, scalar spectral index ns = 0965±0004, and optical depth τ = 0054±0007 (in this abstract we quote 68% conﬁdence regions on measured parameters and 95% on upper limits) The angular acoustic scale is measured to 003% precision, with 100θ∗ = 10411±00003Theseresultsareonlyweaklydependentonthecosmologicalmodelandremainstable,withsomewhatincreasederrors, in many commonly considered extensions Assuming the base-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: HubbleconstantH0 = (674±05)kms−1Mpc−1;matterdensityparameterΩm = 0315±0007;andmatterﬂuctuationamplitudeσ8 = 0811±0006 We ﬁnd no compelling evidence for extensions to the base-ΛCDM model Combining with baryon acoustic oscillation (BAO) measurements (and consideringsingle-parameterextensions)weconstraintheeﬀectiveextrarelativisticdegreesoffreedomtobe Neﬀ = 299±017,inagreementwith the Standard Model prediction Neﬀ = 3046, and ﬁnd that the neutrino mass is tightly constrained toPmν < 012 eV The CMB spectra continue to prefer higher lensing amplitudesthan predicted in base ΛCDM at over 2σ, which pulls some parameters that aﬀect thelensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAOdataThejointconstraintwithBAOmeasurementsonspatialcurvatureisconsistentwithaﬂatuniverse, ΩK = 0001±0002Alsocombining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w0 = −103±003, consistent with a cosmological constant We ﬁnd no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r0002 < 006 Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower ﬂuctuation amplitudes or matter density parameters), and in signiﬁcant, 36σ, tension with local measurements of the Hubble constant (which prefer a higher value) Simple model extensions that can partially resolve these tensions are not favoured by the Planck data

/pdf/planck-2018-results-vi-cosmological-parameters-3d8tp95x87.pdf

Planck 2018 results. VI. Cosmological parameters

On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160)  Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

/pdf/the-observation-of-gravitational-waves-from-a-binary-black-58srpupgg9.pdf

The Observation of Gravitational Waves from a Binary Black Hole Merger

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

/pdf/planck-2018-results-vi-cosmological-parameters-gjusa8fy3b.pdf

We present X-ray and radio monitoring observations of the gamma-ray binary PSR J2032+4127/MT91 213 during its periastron passage in late 2017. Dedicated Chandra, XMM-Newton,NuSTAR X-ray observations and VLA radio observations of this long orbit (50 years), 143 ms pulsar/Be star system clearly revealed flux and spectral variability during the passage. The X-ray spectrum hardened near periastron, with a significant decrease in the power-law photon index from \Gamma ~ 2 to 1.2 and evidence of an increased absorption column density. We identified a possible spectral break at a few keV in the spectrum that suggests synchrotron cooling. A coincident radio and X-ray flare occurred one week after periastron, which is possibly the result of the pulsar wind interacting with the Be stellar disk and generating synchrotron radiation. However, a multi-wavelength comparison indicate that the X-ray and radio spectra cannot be simply connected by a single power-law component. Hence, the emission in these two energy bands must originate from different particle populations.

/pdf/x-ray-and-radio-variabilities-of-psr-j2032-4127-near-127jrp0ayl.pdf

X-ray and Radio Variabilities of PSR J2032+4127 near Periastron

Analysis of Chandra observations of the neutron star (NS) in the centre of the Cassiopeia A supernova remnant taken in the subarray (FAINT) mode of the ACIS detector performed by Posselt and collaborators revealed, after inclusion of the most recent (May 2020) observations, a signiﬁcant decrease of the source surface temperature from 2006 to 2020. The obtained cooling rate is consistent with those obtained from analysis of the 2000–2019 data taken in the GRADED mode of the ACIS detector, which is potentially more strongly aﬀected by instrumental eﬀects. We performed a joint spectral analysis using all ACIS data to constrain the NS parameters and cooling rate. We constrain the mass of the Cassiopeia A NS at 𝑀 = 1 . 55 ± 0 . 25 𝑀 (cid:12) , and its radius at 𝑅 = 13 . 5 ± 1 . 5 km. The surface temperature cooling rate is found to be 2 . 2 ± 0 . 3 per cent in 10 years if the absorbing hydrogen column density is allowed to vary and 1 . 6 ± 0 . 2 per cent in 10 years if it is ﬁxed. The observed cooling can be explained by enhanced neutrino emission from the superﬂuid NS interior due to Cooper Pair Formation (CPF) process. Based on analysis of all ACIS data, we constrain the maximal critical temperature of triplet neutron pairing within the NS core at ( 4 − 9 . 5 ) × 10 8 K. In accordance with previous studies, the required eﬀective strength of the CPF neutrino emission is at least a factor of 2 higher than existing microscopic calculations suggest.

Constraints on neutron star superfluidity from the cooling neutron star in Cassiopeia A using all Chandra ACIS-S observations

Supernovae and Their Consequences: Studies with the Generation-X Mission

Recent observations by XMM-Newton detected rotational pulsations in the total brightness and spectrum of several neutron stars. To properly interpret the data, accurate modeling of neutron star emission is necessary. Detailed analysis of the shape and strength of the rotational variations allows a measurement of the surface composition and magnetic field, as well as constrains the nuclear equation of state. We discuss our models of the spectra and light curves of two of the most observed neutron stars, RX J1856.5-3754 and 1E 1207.4-5209, and discuss some implications of our results and the direction of future work.

Modeling Phase-resolved Observations of the Surfaces of Magnetic Neutron Stars

We study several effects that influence the strength of the proton cyclotron and atomic features in the thermal spectra of magnetic neutron stars. Magnetic field variations over the neutron star surface leads to broadening of features. Vacuum polarization can strongly suppress spectral lines when B ≳ 1014 G. The surface spectrum is unaffected by vacuum polarization when B ≲ 7 × 1013 G; thus the proton cyclotron absorption line (and atomic lines) can have a large equivalent width, possibly explaining the features seen in some isolated neutron stars.

/pdf/modeling-spectral-features-from-isolated-neutron-stars-30fvdb5w5m.pdf

Wynn C. G. Ho

Papers

X-ray and Radio Variabilities of PSR J2032+4127 near Periastron

Constraints on neutron star superfluidity from the cooling neutron star in Cassiopeia A using all Chandra ACIS-S observations

Supernovae and Their Consequences: Studies with the Generation-X Mission

Modeling Phase-resolved Observations of the Surfaces of Magnetic Neutron Stars

Modeling Spectral Features from Isolated Neutron Stars