Astrophysics

In this work, we have searched quasi-periodic oscillations (QPOs) in the 15 GHz light curve of the FSRQ PKS J0805-0111 monitored by the Owens Valley Radio Observatory (OVRO) 40 m telescope during the period from January 9,2008 to May 9,2019, using the weighted wavelet Z-transform (WWZ) and the Lomb-Scargle Periodogram (LSP) techniques. This is the first time to search for periodic radio signal in the FSRQ PKS J0805-0111 by these two methods. All two methods consistently reveal a repeating signal with a periodicity of 3.38 ± 0.8 years ( > 99.7\% confidence level). In order to determining the significance of the periods, the false alarm probability method was used, and a large number of Monte Carlo simulations were performed. As possible explanations, we discuss a number of scenarios including the thermal instability of thin disks scenario, the spiral jet scenario and the binary supermassive black hole scenario, we expected that the binary black hole scenario, where the QPO is caused by the precession of the binary black holes, is the most likely explanation. FSRQ PKS J0805-0111 thus could be a good binary black hole candidate. In the binary black hole scenario, the distance between the primary black hole and the secondary black hole is about 1.71× 10 16 cm.

Gravitational waves from neutron star mergers have long been considered a promising way to measure the Hubble constant, H 0 , which describes the local expansion rate of the Universe. While black hole mergers are more abundantly observed, their expected lack of electromagnetic emission and poor gravitational-wave localization makes them less suited for measuring H 0 . Black hole mergers within the disks of Active Galactic Nuclei (AGN) could be an exception. Accretion from the AGN disk may produce an electromagnetic signal, pointing observers to the host galaxy. Alternatively, the low number density of AGNs could help identify the host galaxy of 1−5% of mergers. Here we show that black hole mergers in AGN disks may be the most sensitive way to determine H 0 with gravitational waves. If 1% of LIGO/Virgo's observations occur in AGN disks with identified host galaxies, we could measure H 0 with 1% uncertainty within five years, likely beyond the sensitivity of neutrons star mergers.

Very long base interferometry (VLBI) radio images recently proved to be essential in breaking the degeneracy in the ejecta model for the neutron star merger event GW170817. We discuss the properties of synthetic radio images of merger jet afterglow by using semi-analytic models of laterally spreading or non-spreading jets. The image centroid initially moves away from the explosion point in the sky with an apparent superlumianal velocity. After reaching a maximum displacement its motion is reversed. This behavior is in line with that found in full hydrodynamics simulations. Since the evolution of the centroid shift and jet image size are significantly different in the two jet models, observations of these characteristics for very bright events might be able to confirm or constrain the lateral expansion law of merger jets. We explicitly demonstrate how θ obs is obtained by the centroid shift of radio images or its apparent velocity provided the ratio of the jet core size θ c and the viewing angle θ obs is determined by afterglow light curves. We show that a simple method based on a point-source approximation provides reasonable angular estimates ( 10??0% errors at most). By taking a sample of structured Gaussian jet results, we find that the model with θ obs ??.32 rad can explain the main features of the GW170817 afterglow light curves and the radio images.

Accretion discs around black holes power some of the most luminous objects in the Universe. Discs that are misaligned to the black hole spin can become warped over time by Lense-Thirring precession. Recent work has shown that strongly warped discs can become unstable, causing the disc to break into discrete rings producing a more dynamic and variable accretion flow. In a companion paper, we present numerical simulations of this instability and the resulting dynamics. In this paper, we discuss the implications of this dynamics for accreting black hole systems, with particular focus on the variability of Active Galactic Nuclei (AGN). We discuss the timescales on which variability might manifest, and the impact of the observer orientation with respect to the black hole spin axis. When the disc warp is unstable near the inner edge of the disc, we find quasi periodic behaviour of the inner disc which may explain the recent quasi periodic eruptions observed in, for example, the Seyfert 2 galaxy GSN 069 and in the galactic nucleus of RX J1301.9+2747. These eruptions are thought to be similar to the `heartbeat' modes observed in some X-ray binaries (e.g. GRS 1915+105 and IGR J17091-3624). When the instability manifests at larger radii in the disc, we find that the central accretion rate can vary on timescales that may be commensurate with, e.g., changing-look AGN. We therefore suggest that some of the variability properties of accreting black hole systems may be explained by the disc being significantly warped, leading to disc tearing.

We present numerical simulations of misaligned discs around a spinning black hole covering a range of parameters. Previous simulations have shown that discs that are strongly warped by a forced precession -- in this case the Lense-Thirring effect from the spinning black hole -- can break apart into discrete discs or rings that can behave quasi-independently for short timescales. With the simulations we present here, we confirm that thin and highly inclined discs are more susceptible to disc tearing than thicker or low inclination discs, and we show that lower values of the disc viscosity parameter lead to instability at lower warp amplitudes. This is consistent with detailed stability analysis of the warped disc equations. We find that the growth rates of the instability seen in the numerical simulations are similar across a broad range of parameters, and are of the same order as the predicted growth rates. However, we did not find the expected trend of growth rates with viscosity parameter. This may indicate that the growth rates are affected by numerical resolution, or that the wavelength of the fastest growing mode is a function of local disc parameters. Finally, we also find that disc tearing can occur for discs with a viscosity parameter that is higher than predicted by a local stability analysis of the warped disc equations. In this case, the instability manifests differently producing large changes in the disc tilt locally in the disc, rather than the large changes in disc twist that typically occur in lower viscosity discs.

We performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope (LAT) gamma-ray sources using the Giant Metrewave Radio Telescope (GMRT) at 322 and 607 MHz. In this paper we report the discovery of three millisecond pulsars (MSPs), PSR J0248+4230, PSR J1207 ??5050 and PSR J1536 ??4948. We conducted follow up timing observations for around 5 years with the GMRT and derived phase coherent timing models for these MSPs. PSR J0248 + 4230 and J1207 ??5050 are isolated MSPs having periodicities of 2.60 ms and 4.84 ms. PSR J1536-4948 is a 3.07 ms pulsar in a binary system with orbital period of around 62 days about a companion of minimum mass 0.32 solar mass. We also present multi-frequency pulse profiles of these MSPs from the GMRT observations. PSR J1536-4948 is an MSP with an extremely wide pulse profile having multiple components. Using the radio timing ephemeris we subsequently detected gamma-ray pulsations from these three MSPs, confirming them as the sources powering the gamma-ray emission. For PSR J1536-4948 we performed combined radio-gamma-ray timing using around 11.6 years of gamma-ray pulse times of arrivals (TOAs) along with the radio TOAs. PSR J1536-4948 also shows evidence for pulsed gamma-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a >10 GeV point source. The multi-wavelength pulse profiles of all three MSPs offer challenges to models of radio and gamma-ray emission in pulsar magnetospheres.

We present Chandra X-ray observations of 14 radio-loud quasars at redshifts 3<z<4 , selected from a well-defined sample. All quasars are detected in the 0.5-7.0 keV energy band, and resolved X-ray features are detected in five of the objects at distances of 1-12" from the quasar core. The X-ray features are spatially coincident with known radio features for four of the five quasars. This indicates that these systems contain X-ray jets. X-ray fluxes and luminosities are measured, and jet-to-core X-ray flux ratios are estimated. The flux ratios are consistent with those observed for nearby jet systems, suggesting that the observed X-ray emission mechanism is independent of redshift. For quasars with undetected jets, an upper limit on the average X-ray jet intensity is estimated using a stacked image analysis. Emission spectra of the quasar cores are extracted and modeled to obtain best-fit photon indices, and an Fe K emission line is detected from one quasar in our sample. We compare X-ray spectral properties with optical and radio emission in the context of both our sample and other quasar surveys.

PSR J1023+0038 is a rapidly-spinning neutron star with a low-mass-binary companion that switches between a radio pulsar and low-luminosity disk state. In 2013, it transitioned to its current disk state accompanied by brightening at all observed wavelengths. In this state, PSR J1023+0038 now shows optical and X-ray pulsations and abrupt X-ray luminosity switches between discrete 'low' and 'high' modes. Continuum radio emission, denoting an outflow, is also present and brightens during the X-ray low modes. Here, we present a simultaneous optical, ultraviolet (UV) and X-ray campaign comprising Kepler ( 400??00 nm), Hubble Space Telescope ( 180??80 nm), XMM-Newton ( 0.3??0 keV) and NuSTAR ( 3??9 keV). We demonstrate that low and high luminosity modes in the UV band are strictly simultaneous with the X-ray modes and change the UV brightness by a factor of ??5 \% on top of a much brighter persistent UV component. We find strong evidence for UV pulsations (pulse fraction of 0.82±0.19 \%) in the high-mode, with a similar waveform as the X-ray pulsations making it the first known UV millisecond pulsar. Lastly, we find that the optical mode changes occur synchronously with the UV/X-ray mode changes, but optical modes are inverted compared to the higher frequencies. There appear to be two broad-band emission components: one from radio to near-infrared/optical that is brighter when the second component from optical to hard X-rays is dimmer (and vice-versa). We suggest that these components trace switches between accretion into the neutron star magnetosphere (high-energy high-mode) versus ejection of material (low-energy high-mode). Lastly, we propose that optical/UV/X-ray pulsations can arise from a shocked accretion flow channeled by the neutron star's magnetic field.

We report the discovery of 1.97 ms period gamma-ray pulsations from the 75 minute orbital-period binary pulsar now named PSR J1653-0158. The associated Fermi Large Area Telescope gamma-ray source 4FGL J1653.6-0158 has long been expected to harbor a binary millisecond pulsar. Despite the pulsar-like gamma-ray spectrum and candidate optical/X-ray associations -- whose periodic brightness modulations suggested an orbit -- no radio pulsations had been found in many searches. The pulsar was discovered by directly searching the gamma-ray data using the GPU-accelerated Einstein@Home distributed volunteer computing system. The multi-dimensional parameter space was bounded by positional and orbital constraints obtained from the optical counterpart. More sensitive analyses of archival and new radio data using knowledge of the pulsar timing solution yield very stringent upper limits on radio emission. Any radio emission is thus either exceptionally weak, or eclipsed for a large fraction of the time. The pulsar has one of the three lowest inferred surface magnetic-field strengths of any known pulsar with B surf ≈4× 10 7 G. The resulting mass function, combined with models of the companion star's optical light curve and spectra, suggests a pulsar mass ≳2 M ⊙ . The companion is light-weight with mass ∼0.01 M ⊙ , and the orbital period is the shortest known for any rotation-powered binary pulsar. This discovery demonstrates the Fermi Large Area Telescope's potential to discover extreme pulsars that would otherwise remain undetected.

We have discovered a new candidate redback millisecond pulsar binary near the center of the error ellipse of the bright unassociated Fermi-LAT γ -ray source 4FGL J0940.3-7610. The candidate counterpart is a variable optical source that also shows faint X-ray emission. Optical photometric and spectroscopic monitoring with the SOAR telescope indicates the companion is a low-mass star in a 6.5-hr orbit around an invisible primary, showing both ellipsoidal variations and irradiation and consistent with the properties of known redback millisecond pulsar binaries. Given the orbital parameters, preliminary modeling of the optical light curves suggests an edge-on inclination and a low-mass ( ??.2 - 1.4 M ??) neutron star, along with a secondary mass somewhat more massive than typical ??.4 M ??. This combination of inclination and secondary properties could make radio eclipses more likely for this system, explaining its previous non-discovery in radio pulsation searches. Hence 4FGL J0940.3-7610 may be a strong candidate for a focused search for γ -ray pulsations to enable the future detection of a millisecond pulsar.