The Fermi Large Area Telescope gamma-ray source 3FGL J2039.6−5618 contains a periodic optical and X-ray source that was predicted to be a “redback” millisecond pulsar (MSP) binary system. However, the conclusive identification required the detection of pulsations from the putative MSP. To better constrain the orbital parameters for a directed search for gamma-ray pulsations, we obtained new optical light curves in 2017 and 2018, which revealed long-term variability from the companion star. The resulting orbital parameter constraints were used to perform a targeted gamma-ray pulsation search using the Einstein@Home distributed volunteer computing system. This search discovered pulsations with a period of 2.65 ms, confirming the source as a binary MSP now known as PSR J2039−5617. Optical light curve modelling is complicated, and likely biased, by asymmetric heating on the companion star and long-term variability, but we find an inclination i ≳ 60 ○, for a low pulsar mass between 1.1 M⊙ < Mpsr < 1.6 M⊙, and a companion mass of 0.15–0.22 M⊙, confirming the redback classification. Timing the gamma-ray pulsations also revealed significant variability in the orbital period, which we find to be consistent with quadrupole moment variations in the companion star, suggestive of convective activity. We also find that the pulsed flux is modulated at the orbital period, potentially due to inverse Compton scattering between high-energy leptons in the pulsar wind and the companion star’s optical photon field.

Context. The gravitational strong equivalence principle (SEP) is a cornerstone of the general theory of relativity (GR). Hence, testing the validity of SEP is of great importance when confronting GR, or its alternatives, with experimental data. Pulsars that are orbited by white dwarf companions provide an excellent laboratory, where the extreme difference in binding energy between neutron stars and white dwarfs allows for precision tests of the SEP via the technique of radio pulsar timing. Aims. To date, the best limit on the validity of SEP under strong-field conditions was obtained with a unique pulsar in a triple stellar system, PSR J0337+1715. We report here on an improvement of this test using an independent data set acquired over a period of 6 years with the Nançay radio telescope. The improvements arise from a uniformly sampled data set, a theoretical analysis, and a treatment that fixes some short-comings in the previously published results, leading to better precision and reliability of the test. Methods. In contrast to the previously published test, we use a different long-term timing data set, developed a new timing model and an independent numerical integration of the motion of the system, and determined the masses and orbital parameters with a different methodology that treats the parameter Δ, describing a possible strong-field SEP violation, identically to all other parameters. Results. We obtain a violation parameter Δ = ( + 0.5 ± 1.8) × 10−6 at 95% confidence level, which is compatible with and improves upon the previous study by 30%. This result is statistics-limited and avoids limitation by systematics as previously encountered. We find evidence for red noise in the pulsar spin frequency, which is responsible for up to 10% of the reported uncertainty. We use the improved limit on SEP violation to place constraints on a class of well-studied scalar-tensor theories, in particular we find ωBD >  140 000 for the Brans-Dicke parameter. The conservative limits presented here fully take into account current uncertainties in the equation for state of neutron-star matter.

For 80 days in 2017, the Kepler Space Telescope continuously observed the transitional millisecond pulsar system PSR J1023+0038 in its accreting state. We present analyses of the 59-second cadence data, focusing on investigations of the orbital light curve of the irradiated companion star, and of flaring activity in the neutron star's accretion disc. The underlying orbital modulation from the companion star retains a similar amplitude and asymmetric heating profile as seen in previous photometric observations of the system in its radio pulsar state, suggesting that the heating mechanism has not been affected by the state change. We also find tentative evidence that this asymmetry may vary with time. The light curve also exhibits "flickering" activity, evident as short time-scale flux correlations throughout the observations, and periods of rapid modeswitching activity on time scales shorter than the observation cadence. Finally, the system spent ∼ 20% of the observations in a flaring state, with the length of these flares varying from < 2 minutes up to several hours. The flaring behaviour is consistent with a self-organised criticality mechanism, most likely related to the build up and release of mass at the inner edge of the accretion disc.

Kelvin probe force microscopy at normal pressure was performed by two different groups on the same Au-coated planar sample used to measure the Casimir interaction in a sphere-plane geometry. The obtained voltage distribution was used to calculate the separation dependence of the electrostatic pressure Pres(D) in the configuration of the Casimir experiments. In the calculation it was assumed that the potential distribution in the sphere has the same statistical properties as the measured one, and that there are no correlation effects on the potential distributions due to the presence of the other surface. The result of this calculation, using the currently available knowledge, is that Pres(D) does not explain the magnitude or the separation dependence of the difference ∆P (D) between the measured Casimir pressure and the one calculated using a Drude model for the electromagnetic response of Au. We discuss in the conclusions the points which have to be checked out by future work, including the influence of pressure and a more accurate determination of the patch distribution, in order to confirm these results.

Spider pulsars are binary systems containing an energetic millisecond pulsar that intensely irradiates a closely orbiting low-mass companion. Modelling their companion’s optical light curves is essential to the study of the orbital properties of the binary, including the determination of the pulsar mass, characterising the pulsar wind and the star itself. We aim to generalise the traditional direct heating model of irradiation, whereby energy deposited by the pulsar wind into the stellar envelope is locally re-emitted, by introducing heat redistribution via diffusion and convection within the outer stellar envelope. We approximate the irradiated stellar envelope as a two-dimensional shell. This allows us to propose an effective equation of energy conservation that can be solved at a reduced computational cost. We then implement this model in the Icarus software and use evidence sampling to determine the most likely convection and diffusion laws for the light curve of the redback companion of PSR J2215+5135. Redistribution effects concentrate near the terminator line of pulsar irradiation, and can create apparent hot and cold spots. Among the models tested for PSR J2215+5135, we find that all models with heat redistribution are more likely than symmetric direct heating. The best-fitting redistribution model involves diffusion together with a uniformly rotating envelope. However, we caution that all models still present serious systematic effects, and that prior knowledge from pulsar timing, spectroscopy and distance are key to determine with certainty the most accurate redistribution law. We propose an extension of the direct heating framework that allows for exploring a variety of heat redistribution effects. Future work is necessary to determine the relevant laws from first principles and empirically using complementary observations.

Kelvin probe force microscopy at normal pressure was performed by two different groups on the same Au-coated planar sample used to measure the Casimir interaction in a sphere-plane geometry. The obtained voltage distribution was used to calculate the separation dependence of the electrostatic pressure Pres(D) in the configuration of the Casimir experiments. In the calculation it was assumed that the potential distribution in the sphere has the same statistical properties as the measured one, and that there are no correlation effects on the potential distributions due to the presence of the other surface. The result of this calculation, using the currently available knowledge, is that Pres(D) does not explain the magnitude or the separation dependence of the difference ∆P (D) between the measured Casimir pressure and the one calculated using a Drude model for the electromagnetic response of Au. We discuss in the conclusions the points which have to be checked out by future work, including the influence of pressure and a more accurate determination of the patch distribution, in order to confirm these results.

Aims. We investigate the evaporation of close-by pulsar companions, such as planets, asteroids, and white dwarfs, by induction heating. Methods. Assuming that the outflow energy is dominated by a Poynting flux (or pulsar wave) at the location of the companions, we calculate their evaporation timescales, by applying the Mie theory. Results. Depending on the size of the companion compared to the incident electromagnetic wavelength, the heating regime varies and can lead to a total evaporation of the companion. In particular, we find that inductive heating is mostly inefficient for small pulsar companions, although it is generally considered the dominant process. Conclusions. Small objects like asteroids can survive induction heating for 10 4 yr at distances as small as 1 R from the neutron star. For degenerate companions, induction heating cannot lead to evaporation and another source of heating (likely by kinetic energy of the pulsar wind) has to be considered. It was recently proposed that bodies orbiting pulsars are the cause of fast radio bursts; the present results explain how those bodies can survive in the pulsar's highly energetic environment.

Spider millisecond pulsars are, along with some eclipsing post-common envelope systems and cataclysmic variables, part of an expanding category of compact binaries with low-mass companions for which puzzling timing anomalies have been observed. The most prominent type of irregularities seen in them are orbital period variations, a phenomenon which has been proposed to originate from changes in the gravitational quadrupole moment of the companion star. A physically sound modelling of the timing of these systems is key to understanding their structure and evolution. In this paper we argue that a complete timing model must account for relativistic corrections as well as rotationally and tidally induced quadrupole distortions. We solve for the resulting orbital dynamics using perturbation theory and derive the corresponding timing model in the low eccentricity limit. We find that the expected strong quadrupole deformation of the companion star results in an effective minimum orbital eccentricity. It is accompanied by a fast periastron precession which, if not taken into account, averages out any measurement of the said eccentricity. We show that, with our model, detection of both eccentricity and precession is likely to be made in many if not all spider pulsar systems. Combined with optical light curves, this will allow us to measure the apsidal motion constant, connecting the quadrupole deformation to the internal structure, and thus opening a new window into probing the nature of their exotic stellar interiors. Moreover, more accurate timing may eventually lead spider pulsars to be used for high-precision timing experiments such as pulsar timing arrays.