Gravitational microlensing can detect isolated stellar-mass black holes (BHs), which are believed to be the dominant form of Galactic BHs according to population synthesis models. Previous searches for BH events in microlensing data focused on long-timescale events with significant microlensing parallax detections. Here we show that, although BH events preferentially have long timescales, the microlensing parallax amplitudes are so small that in most cases the parallax signals cannot be detected statistically significantly. We then identify OGLE-2006-BLG-044 to be a candidate BH event because of its long timescale and small microlensing parallax. Our findings have implications to future BH searches in microlensing data.

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269, with planet-host mass ratio $q \sim 6\times10^{-4}$, i.e., $0.6$ times smaller than the Jupiter/Sun mass ratio. Combined with the $Gaia$ parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant $(M_{2} = 0.69_{-0.22}^{+0.44}\,M_{\rm J})$ orbiting a Sun-like star $(M_{1} = 1.13_{-0.35}^{+0.72}\,M_{\odot})$ at a distance of $D_{\rm L} = 2.56_{-0.62}^{+0.92}\,{\rm kpc}$. The projected planet-host separation is $a_{\perp} = 4.61_{-1.17}^{+1.70}\,{\rm au}$. Using {\it Gaia} astrometry, we show that the blended light lies $\lesssim 12\,$mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at $>99.6\%$. The host is therefore a subgiant. For host metallicities in the range of $0.0 \leq {\rm [Fe/H]} \leq +0.3$, the host and planet masses are then in the range of $1.16 \leq M_{1}/M_{\odot} \leq 1.38$ and $0.74 \leq M_{2}/M_{\rm J} \leq 0.89$, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic followup both to test the microlensing solution and to further characterize the system.

We report the discovery and analysis of the planetary microlensing event OGLE-2017-BLG-0406, which was observed both from the ground and by the ${\it Spitzer}$ satellite in a solar orbit. At high magnification, the anomaly in the light curve was densely observed by ground-based-survey and follow-up groups, and it was found to be explained by a planetary lens with a planet/host mass ratio of $q=7.0 \times 10^{-4}$ from the light-curve modeling. The ground-only and ${\it Spitzer}$-"only" data each provide very strong one-dimensional (1-D) constraints on the 2-D microlens parallax vector $\bf{\pi_{\rm E}}$. When combined, these yield a precise measurement of $\bf{\pi_{\rm E}}$, and so of the masses of the host $M_{\rm host}=0.56\pm0.07\,M_\odot$ and planet $M_{\rm planet} = 0.41 \pm 0.05\,M_{\rm Jup}$. The system lies at a distance $D_{\rm L}=5.2 \pm 0.5 \ {\rm kpc}$ from the Sun toward the Galactic bulge, and the host is more likely to be a disk population star according to the kinematics of the lens. The projected separation of the planet from the host is $a_{\perp} = 3.5 \pm 0.3 \ {\rm au}$, i.e., just over twice the snow line. The Galactic-disk kinematics are established in part from a precise measurement of the source proper motion based on OGLE-IV data. By contrast, the ${\it Gaia}$ proper-motion measurement of the source suffers from a catastrophic $10\,\sigma$ error.

The first measurements from the CIBER experiment of extragalactic background light (EBL) in near-infrared (NIR) band exhibit a higher intensity than those inferred through $\gamma$-ray observations. Recent theoretical-EBL intensities are typically consistent with the very high energy (VHE) $\gamma$-ray observations. Yet, it is possible that the excess NIR radiation is a new component of EBL and not in tension with the TeV spectra of distant blazars, since the hypothetical axion-like particle (ALP) may lead to a reduced opacity of the Universe for VHE $\gamma$-rays. In order to probe whether the excess component arises mainly from EBL, thirteen observed spectra in high energy and VHE ranges from ten distant TeV BL Lac objects are fitted by four theoretical spectra which involve theoretical EBL (Gilmore ), Gilmore's EBL model including photon/ALP coupling, Gilmore's EBL with CIBER excess and the latter including photon/ALP coupling respectively. We find the goodness of fit for the model with CIBER excess can be improved with a significance of $7.6~\sigma$ after including the photon/ALP coupling; Thus, the ALP/photon mixing mechanism can effectively alleviate the tension; However, the Gilmore EBL model, on the whole, is more compatible with the observed spectra compared to those with ALP, although individual blazars such as PKS 1424+240 and 1ES 1101-232 prefer the ALP-model. Our results suggest that the recent EBL models can solely explain the VHE $\gamma$-ray observations, and assuming the existence of the ALP to alleviate the tension is not required in a statistical sense, thus the excess over the EBL models is less likely to be a new EBL component.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

We present timing solutions and analyses of 11 pulsars discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). These pulsars were discovered using an ultra-wide bandwidth receiver in drift-scan observations made during the commissioning phase of FAST, and were then confirmed and timed using the 64-m Parkes Radio Telescope. Each pulsar has been observed over a span of at least one year. Highlighted discoveries include PSR J0344-0901, which displays mode-changing behaviour and may belong to the class of so-called `swooshing' pulsars (alongside PSRs B0919+06 and B1859+07); PSR J0803-0942, whose emission is almost completely linearly polarised; and PSRs J1900-0134 and J1945+1211, whose well defined polarisation angle curves place stringent constraints on their emission geometry. We further discuss the detectability of these pulsars by earlier surveys, and highlight lessons learned from our work in carrying out confirmation and monitoring observations of pulsars discovered by a highly sensitive telescope, many of which may be applicable to next-generation pulsar surveys. This paper marks one of the first major releases of FAST-discovered pulsars, and paves the way for future discoveries anticipated from the Commensal Radio Astronomy FAST Survey (CRAFTS).

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

To investigate the use of saliency-map analysis to aid in searches for transient signals, such as fast radio bursts and individual pulses from radio pulsars. We aim to demonstrate that saliency maps provide the means to understand predictions from machine learning algorithms and can be implemented in piplines used to search for transient events. We have implemented a new deep learning methodology to predict whether or not any segment of the data contains a transient event. The algorithm has been trained using real and simulated data sets. We demonstrate that the algorithm is able to identify such events. The output results are visually analysed via the use of saliency maps. We find that saliency maps can produce an enhanced image of any transient feature without the need for de-dispersion or removal of radio frequency interference. Such maps can be used to understand which features in the image were used in making the machine learning decision and to visualise the transient event. Even though the algorithm reported here was developed to demonstrate saliency-map analysis, we have detected, in archival data, a single burst event with dispersion measure of $41$\,cm$^{-3}$pc that is not associated with any currently known pulsar.

We reveal that the energy spectra of electrons-positrons in primary cosmic rays measured at atmosphere top have double structures: an excess component $\Phi^s_{e^+}(E)=\Phi^s_{e^-}(E)$ around $400 GeV$, which origins from a strong $e^+e^-$-source and the distorted background $\Phi^0_{e^-}(E)$. We supposed that the difference between AMS-CALET and Fermi-LAT-DAMPE data origins from the energy loss of the fluxes due to the anomalous bremsstrahlung effect at a special window. The evolution of spectra under anomalous bremsstrahlung effect satisfies an improved electromagnetic cascade equation. The above spectra are parameterized and they can be regarded as the subjects exploring new physics. We suggest to check the previous applications of the Bethe-Heitler formula in the study of the propagation of high energy electrons and photons.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.

We carry out a comparison study on the bar structure in the Illustris-1 and TNG100 simulations. At $z=0$, 8.9\% of 1232 disc galaxies with stellar mass $>10^{10.5}M_{\odot}$ in Illustris-1 are barred, while the numbers are 55\% of 1269 in TNG100. The bar fraction as a function of stellar mass in TNG100 agrees well with the survey $S^4G$. The median redshift of bar formation are $\sim 0.4-0.5$ and $\sim 0.25$ in TNG100 and Illustris-1 respectively. Bar fraction generally increases with stellar mass and decreases with gas fraction in both simulations. Barred galaxy had higher gas fraction at high redshift tend to form bar later. When the bars were formed, the disc gas fractions were mostly lower than 0.4. The much higher bar fraction in TNG100 probably have benefit from much lower gas fraction in massive disc galaxies since $z\sim3$, which may result from the combination of more effective stellar and AGN feedback. The latter may be the primary factor at $z<2$. Meanwhile, in both simulations, barred galaxies have higher star formation rate before bar formation, and stronger AGN feedback all the time than unbarred galaxies. The properties of dark matter halos hosting massive disc galaxies are similar between two simulations, and should have minor effect on the different bar frequency. For individual galaxies under similar halo environment cross two simulations, different baryonic physics can lead to striking discrepancy on morphology. The morphology of individual galaxies is subject to combined effects of environment and internal baryonic physics, and is often not predictable.

We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon.

We report the discovery of an eclipsing binary millisecond pulsar in the globular cluster M92 (NGC6341) with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). PSR J1717+4308A, or M92A, has a pulse frequency of 316.5~Hz (3.16~ms) and a dispersion measure of 35.45 pc cm$^{-3}$. The pulsar is a member of a binary system with an orbital period of 0.20~days around a low-mass companion which has a median mass of $\sim$0.18~\Ms. From observations so far, at least two eclipsing events have been observed in each orbit. The longer one lasted for ~5000~s in the orbital phase range 0.1--0.5. The other lasted for ~500~s and occurred between 1000--2000~s before or after the longer eclipsing event. The lengths of these two eclipsing events also change. These properties suggest that J1717+4308A is a ``red-back'' system with a low-mass main sequence or sub-giant companion. Timing observations of the pulsar and further searches of the data for additional pulsars are ongoing.

Hardening of cosmic proton- and nuclei-spectra is explained by using the gluon condensation (GC) model, which states that a large amount of gluons in proton may condense near the high energy threshold. The results present the universality of the GC-effect, which may be a common origin of a series of anomalous astrophysical phenomena including the broken power-law in gamma-ray spectra, the excess in positron- and electron-spectra and hardening of proton- and nuclei-spectra.

The Bethe-Heitler formula describes bremsstrahlung and it's a typical and important example in quantum electromagnetic dynamics (QED). This formula is widely applied in many branches of physics and astrophysics. We find that the integrated bremsstrahlung cross section at the static approximation and high energy limit has an unexpected big increment, which is missed by previous bremsstrahlung theory. This anomalous effect also exists in electron-positron pair creation. We derive the relating formulas and point out that electromagnetic cascades at the top of atmosphere can test this effect.

Taking a large amount of images, the Cassini Imaging Science Subsystem (ISS) has been routinely used in astrometry. In ISS images, disk-resolved objects often lead to false detection of stars that disturb the camera pointing correction. The aim of this study was to develop an automated processing method to remove the false image stars in disk-resolved objects in ISS images. The method included the following steps: extracting edges, segmenting boundary arcs, fitting circles and excluding false image stars. The proposed method was tested using 200 ISS images. Preliminary experimental results show that it can remove the false image stars in more than 95% of ISS images with disk-resolved objects in a fully automatic manner, i.e. outperforming the traditional circle detection based on Circular Hough Transform (CHT) by 17%. In addition, its speed is more than twice as fast as that of the CHT method. It is also more robust (no manual parameter tuning is needed) when compared with CHT. The proposed method was also applied to a set of ISS images of Rhea to eliminate the mismatch in pointing correction in automatic procedure. Experiment results showed that the precision of final astrometry results can be improve by roughly 2 times than that of automatic procedure without the method. It proved that the proposed method is helpful in the astrometry of ISS images in fully automatic manner.

Wave effects are often neglected in microlensing studies; however, for coherent point-like sources, such as pulsars and fast radio bursts (FRBs), wave effects will become important in their gravitational lensing. In this paper, we describe the wave optics formalism, its various limits, and the conditions for which these limits hold. Using the simple point lens as an example, we will show that the frequency dependence of wave effects breaks degeneracies that are present in the usual geometric optics limit, and constructive interference results in larger magnifications further from the lens. This latter fact leads to a generic increase in cross section for microlensing events in the wave-optics regime compared to the geometric optics regime. For realistic percent-level spectral sensitivities, this leads to a relative boost in lensing cross section of more than an order of magnitude. We apply the point-lens model to the lensing of FRBs and pulsars and find that these radio sources will be lensed in the full wave-optics regime by isolated masses in the range of $0.1-100\,M_\oplus$, which includes free-floating planets (FFPs), whose Einstein radius is smaller than the Fresnel scale. More generally, the interference pattern allows an instantaneous determination of lens masses, unlike traditional microlensing techniques which only yield a mass inference from the event timescale.

Recent studies claimed that planets around the same star have similar sizes and masses and regular spacings, and that planet pairs usually show ordered sizes such that the outer planet is usually the larger one. Here I show that these patterns can be largely explained by detection biases. The \emph{Kepler} planet detections are set by the transit signal-to-noise ratio (S/N). For different stellar properties and orbital period values, the same S/N corresponds to different planetary sizes. This variation in the detection threshold naturally leads to apparent correlations in planet sizes and the observed size ordering. The apparently correlated spacings, measured in period ratios, between adjacent planet pairs in systems with at least three detected planets are partially due to the arbitrary upper limit that the earlier study imposed on the period ratio, and partially due to the varying stability threshold for different planets. After these detection biases are taken into account, we do not find strong evidence for the so-called "intra-system uniformity" or the size ordering effect. Instead, the physical properties of \emph{Kepler} planets are largely independent of the properties of their siblings and the parent star. It is likely that the dynamical evolution has erased the memory of \emph{Kepler} planets about their initial formation conditions. In other words, it will be difficult to infer the initial conditions from the observed properties and the architecture of \emph{Kepler} planets.

The thermal history of cosmic gas in the Dark Ages remains largely unknown. It is important to quantify the impact of relevant physics on the IGM temperature between $z=10$ and $z \sim 30$, in order to interpret recent and oncoming observations, including results reported by EDGES. We revisit the gas heating due to structure formation shocks in this era, using a set of fixed grid cosmological hydrodynamical simulations performed by three different codes. In all our simulations, the cosmic gas is predicted to be in multiphase state since $z>30$. The gas surrounding high density peaks gradually develops a relation more sharp than $T \propto \rho^{2/3}$, approximately $T \propto \rho^{2}$, from $z=30$ to $z=11$, might due to shock heating. Meanwhile, the gas in void region tends to have a large local mach number, and their thermal state varies significantly from code to code. In the redshift range $11-20$, the mass fraction of gas shock heated above the CMB temperature in our simulations is larger than previous semi-analytical results by a factor of 2 to 8. At $z=15$, the fraction varies from $\sim 19\%$ to $52 \%$ among different codes. Between $z=11$ and $z=20$, the gas temperature $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be $\sim 10-20$ K by two codes, much higher than the adiabatic cooling model and some previous works. However, in our simulations performed by RAMSES, $<1/T_{\rm{K}}>_M^{-1}$ is predicted to be even below the temperature required to explain result of the EDGES. Given the fact that different codes give different predictions, currently, it seems a challenge to make solid prediction on the temperature of gas at $z \sim 17$ in simulations.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

The gluon condensation effect in cosmic proton-proton collisions at high energy is used to explain an excess in the positron spectrum observed by AMS. We find that this excess may originate from the GC-effect in Tycho's supernova remnant.

At $q=1.81\pm 0.20 \times 10^{-5}$, KMT-2018-BLG-0029Lb has the lowest planet-host mass ratio $q$ of any microlensing planet to date by more than a factor of two. Hence, it is the first planet that probes below the apparent "pile-up" at $q=5$--10 $\times 10^{-5}$. The event was observed by {\it Spitzer}, yielding a microlens-parallax $\pi_{\rm E}$ measurement. Combined with a measurement of the Einstein radius $\theta_{\rm E}$ from finite-source effects during the caustic crossings, these measurements imply masses of the host $M_{\rm host}=1.14^{+0.10}_{-0.12}\, M_\odot$ and planet $M_{\rm planet} = 7.59^{+0.75}_{-0.69}\,M_\oplus$, system distance $D_L = 3.38^{+0.22}_{-0.26}\,\,{\rm kpc}$ and projected separation $a_\perp = 4.27^{+0.21}_{-0.23}\,{\rm au}$. The blended light, which is substantially brighter than the microlensed source, is plausibly due to the lens and could be observed at high resolution immediately.

We report the discovery of a super-Jovian planet in the microlensing event KMT-2016-BLG-1836, which was found by the Korea Microlensing Telescope Network's high-cadence observations (\Gamma ~ 4~{hr}^{-1}). The planet-host mass ratio q ~ 0.004. A Bayesian analysis indicates that the planetary system is composed of a super-Jovian M_{planet} = 2.2_{-1.1}^{+1.9} M_{J} planet orbiting an M or K dwarf M_{\rm host} = 0.49_{-0.25}^{+0.38} M_{Sun}, at a distance of D_{L} = 7.1_{-2.4}^{+0.8} kpc. The projected planet-host separation is 3.5^{+1.1}_{-0.9} AU, implying that the planet is located beyond the snowline of the host star. Future high-resolution images can potentially strongly constrain the lens brightness and thus the mass and distance of the planetary system. Without considering detailed detection efficiency, selection or publication biases, we find a potential "mass ratio desert" at -3.7 \lesssim \log q \lesssim -3.0 for the 31 published KMTNet planets.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

We report the lens mass and distance measurements of the nearby microlensing event TCP J05074264+2447555. We measure the microlens parallax vector ${\pi}_{\rm E}$ using Spitzer and ground-based light curves with constraints on the direction of lens-source relative proper motion derived from Very Large Telescope Interferometer (VLTI) GRAVITY observations. Combining this ${\pi}_{\rm E}$ determination with the angular Einstein radius $\theta_{\rm E}$ measured by VLTI GRAVITY observations, we find that the lens is a star with mass $M_{\rm L} = 0.495 \pm 0.063~M_{\odot}$ at a distance $D_{\rm L} = 429 \pm 21~{\rm pc}$. We find that the blended light basically all comes from the lens. The lens-source proper motion is $\mu_{\rm rel,hel} = 26.55 \pm 0.36~{\rm mas\,yr^{-1}}$, so with currently available adaptive-optics (AO) instruments, the lens and source can be resolved in 2021. This is the first microlensing event whose lens mass is unambiguously measured by interferometry + satellite parallax observations, which opens a new window for mass measurements of isolated objects such as stellar-mass black holes.

In order to search for further observational evidence of cloud-cloud collisions in one of the promising candidates, L1188, we carried out observations of multiple molecular lines toward the intersection region of the two nearly orthogonal filamentary molecular clouds in L1188. Based on these observations, we find two parallel filamentary structures, both of which have at least two velocity components being connected with broad bridging features. We also found a spatially complementary distribution between the two molecular clouds, as well as enhanced $^{13}$CO emission and $^{12}$CO self-absorption toward their abutting regions. At the most blueshifted velocities, we unveil a 1~pc-long arc ubiquitously showing $^{12}$CO line wings. We discover two 22 GHz water masers, which are the first maser detections in L1188. An analysis of line ratios at a linear resolution of 0.2 pc suggests that L1188 is characterised by kinetic temperatures of 13--23~K and H$_{2}$ number densities of 10$^{3}$--10$^{3.6}$ cm$^{-3}$. On the basis of previous theoretical predictions and simulations, we suggest that these observational features can be naturally explained by the scenario of a cloud-cloud collision in L1188, although an additional contribution of stellar feedback from low-mass young stellar objects cannot be ruled out.

We report the discovery and the analysis of the short (tE < 5 days) planetary microlensing event, OGLE-2015-BLG-1771. The event was discovered by the Optical Gravitational Lensing Experiment (OGLE), and the planetary anomaly (at I ~ 19) was captured by The Korea Microlensing Telescope Network (KMTNet). The event has three surviving planetary models that explain the observed light curves, with planet-host mass ratio q \~ 5.4 * 10^{-3}, 4.5 * 10^{-3} and 4.5 * 10^{-2}, respectively. The first model is the best-fit model, while the second model is disfavored by \Delta\chi^2 ~ 3. The last model is strongly disfavored by \Delta\chi^2 ~ 15 but not ruled out. A Bayesian analysis using a Galactic model indicates that the first two models are probably composed of a Saturn-mass planet orbiting a late M dwarf, while the third one could consist of a super-Jovian planet and a mid-mass brown dwarf. The source-lens relative proper motion is mu_rel ~ 9 mas/yr, so the source and lens could be resolved by current adaptive-optics (AO) instruments in 2021 if the lens is luminous.

The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 ($\mu_{\rm cl}$(N,E) = (+0.36+-0.10, +1.42+-0.10) mas/yr) as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using OGLE photometry is consistent with the value found based on the light curve displacement between Earth and Spitzer.

We report the first wideband observations of pulsars C, D and J in the globular cluster 47Tucanae (NGC 104) using the Ultra-Wideband Low (UWL) receiver system recently installed on the Parkes 64 m radio telescope. The wide frequency range of the UWL receiver (704-4032 MHz), along with the well-calibrated system, allowed us to obtain flux density measurements and polarization pulse profiles. The mean pulse profiles have significant linear and circular polarization, allowing for determination of the Faraday rotation measure for each pulsar. Precise measurements of the dispersion measures show a significant deviation in the value for pulsar D compared to earlier results. Searches for new pulsars in the cluster are on-going and we have determined optimal bands for such searches using the Parkes UWL receiver system.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.

The rate of fast radio bursts (FRBs) in the direction of nearby galaxy clusters is expected to be higher than the mean cosmological rate if intrinsically faint FRBs are numerous. In this paper, we describe a targeted search for faint FRBs near the core of the Virgo cluster using the Australian Square Kilometer Array Pathfinder telescope. During 300 hr of observations, we discovered one burst, FRB 180417, with dispersion measure DM $=474.8$ cm$^{-3}$pc. The FRB was promptly followed up by several radio telescopes for 27 h, but no repeat bursts were detected. An optical follow-up of FRB 180417 using the PROMPT5 telescope revealed no new sources down to an $R$-band magnitude of 20.1. We argue that FRB 180417 is likely behind the Virgo cluster as the Galactic and intracluster DM contribution are small compared to the DM of the FRB, and there are no galaxies in the line of sight. The non-detection of FRBs from Virgo constrains the faint-end slope, $\alpha<1.52$ (at 68\% confidence limit), and the minimum luminosity, $L_{\rm min}\gtrsim 2\times 10^{40}$ erg s$^{-1}$ (at 68\% confidence limit), of the FRB luminosity function assuming cosmic FRB rate of $10^4$ FRBs sky$^{-1}$ day$^{-1}$ with flux above 1 Jy located out to redshift of 1. Further FRB surveys of galaxy clusters with high-sensitivity instruments will tighten the constraints on the faint end of the luminosity function and, thus, are strongly encouraged.

Nuclear Coulomb potential at completely ionized and extremely thin atmosphere can leaks to a macroscopic spatial scale. This effect is used to explain the difference between the energy spectra of cosmic-ray electrons around $1 TeV$ measured by different experimental groups. The result inspires us to review the traditional electromagnetic shower theory at the extreme conditions. It is also reminds us that the energy spectrum of cosmic-ray electrons, which are measured by Fermi-LAT and DAMPE at a higher altitude is more closer to a true signature.

In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.

We report on eight millisecond pulsars (MSPs) in binary systems discovered with the Arecibo PALFA survey. Phase-coherent timing solutions derived from 2.5 to 5 years of observations carried out at Arecibo and Jodrell Bank observatories are provided. PSR J1921+1929 is a 2.65-ms pulsar in a 39.6-day orbit for which we detect $\gamma$-ray pulsations in archival Fermi data. PSR J1928+1245 is a very low-mass-function system with an orbital period of 3.3 hours that belongs to the non-eclipsing black widow population. We also present PSR J1932+1756, the longest-orbital-period (41.5 days) intermediate-mass binary pulsar known to date. In light of the numerous discoveries of binary MSPs over the past years, we characterize the Galactic distribution of known MSP binaries in terms of binary class. Our results support and strengthen previous claims that the scatter in the Galactic scale height distribution correlates inversely with the binary mass function. We provide evidence of observational biases against detecting the most recycled pulsars near the Galactic plane, which overestimates the scale height of lighter systems. A possible bimodality in the mass function of MSPs with massive white dwarfs is also reported.

We present the analysis of the event OGLE-2017-BLG-1186 from the 2017 Spitzer microlensing campaign. This is a remarkable microlensing event because its source is photometrically bright and variable, which makes it possible to perform an asteroseismic analysis using ground-based data. We find that the source star is an oscillating red giant with average timescale of $\sim 9$ d. The asteroseismic analysis also provides us source properties including the source angular size ($\sim 27~\mu{\rm as}$) and distance ($\sim 11.5$ kpc), which are essential for inferring the properties of the lens. When fitting the light curve, we test the feasibility of Gaussian Processes (GPs) in handling the correlated noise caused by the variable source. We find that the parameters from the GP model are generally more loosely constrained than those from the traditional $\chi^2$ minimization method. We note that this event is the first microlensing system for which asteroseismology and GPs have been used to account for the variable source. With both finite-source effect and microlens parallax measured, we find that the lens is likely a $\sim 0.045~M_{\odot}$ brown dwarf at distance $\sim 9.0$ kpc, or a $\sim 0.073~M_{\odot}$ ultracool dwarf at distance $\sim 9.8$ kpc. Combining the estimated lens properties with a Bayesian analysis using a Galactic model, we find a $\sim 35$ per cent probability for the lens to be a bulge object and $\sim 65$ per cent to be a background disc object.

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state.

In this paper, we present the results of timing observations of PSRs J1949+3106 and J1950+2414, two binary millisecond pulsars discovered in data from the Arecibo ALFA pulsar survey (PALFA). The timing parameters include precise measurements of the proper motions of both pulsars, which show that PSR J1949+3106 has a transversal motion very similar to that of an object in the local standard of rest. The timing also includes measurements of the Shapiro delay and the rate of advance of periastron for both systems. Assuming general relativity, these allow estimates of the masses of the components of the two systems; for PSR J1949+3106, the pulsar mass is $M_p \, = \, 1.34^{+0.17}_{-0.15} \, M_{\odot}$ and the companion mass $M_c \, = \, 0.81^{+0.06}_{-0.05}\, M_{\odot}$; for PSR J1950+2414 $M_p \, = \, 1.496 \, \pm \, 0.023\, M_{\odot}$ and $M_c \, = \, 0.280^{+0.005}_{-0.004}\, M_{\odot}$ (all values 68.3 % confidence limits). We use these masses and proper motions to investigate the evolutionary history of both systems: PSR J1949+3106 is likely the product of a low-kick supernova; PSR J1950+2414 is a member of a new class of eccentric millisecond pulsar binaries with an unknown formation mechanism. We discuss the proposed hypotheses for the formations of these systems in light of our new mass measurements.

We report a multiplanetary system found from the analysis of microlensing event OGLE-2018-BLG-1011, for which the light curve exhibits a double-bump anomaly around the peak. We find that the anomaly cannot be fully explained by the binary-lens or binary-source interpretations and its description requires the introduction of an additional lens component. The 3L1S (3 lens components and a single source) modeling yields three sets of solutions, in which one set of solutions indicates that the lens is a planetary system in a binary, while the other two sets imply that the lens is a multiplanetary system. By investigating the fits of the individual models to the detailed light curve structure, we find that the multiple-planet solution with planet-to-host mass ratios $\sim 9.5\times 10^{-3}$ and $\sim 15\times 10^{-3}$ are favored over the other solutions. From the Bayesian analysis, we find that the lens is composed of two planets with masses $1.8^{+3..4}_{-1.1}~M_{\rm J}$ and $2.8^{+5.1}_{-1.7}~M_{\rm J}$ around a host with a mass $0.18^{+0.33}_{-0.10}~M_\odot$ and located at a distance $7.1^{+1.1}_{-1.5}~{\rm kpc}$. The estimated distance indicates that the lens is the farthest system among the known multiplanetary systems. The projected planet-host separations are $a_{\perp,2}=1.8^{+2.1}_{-1.5}~{\rm au}$ ($0.8^{+0.9}_{-0.6}~{\rm au}$) and $a_{\perp,3}=0.8^{+0.9}_{-0.6}~{\rm au}$, where the values of $a_{\perp,2}$ in and out the parenthesis are the separations corresponding to the two degenerate solutions, indicating that both planets are located beyond the snow line of the host, as with the other four multiplanetary systems previously found by microlensing.

We present the analysis of the microlensing event OGLE-2015-BLG-1670, detected in a high-extinction field, very close to the Galactic plane. Due to the dust extinction along the line of sight, this event was too faint to be detected before it reached the peak of magnification. The microlensing light-curve models indicate a high-magnification event with a maximum of $A_\mathrm{max}\gtrsim200$, very sensitive to planetary deviations. An anomaly in the light curve has been densely observed by the microlensing surveys MOA, KMTNet, and OGLE. From the light-curve modeling, we find a planetary anomaly characterized by a planet-to-host mass ratio, $q=\left(1.00^{+0.18}_{-0.16}\right)\times 10^{-4}$, at the peak recently identified in the mass-ratio function of microlensing planets. Thus, this event is interesting to include in future statistical studies about planet demography. We have explored the possible degeneracies and find two competing planetary models resulting from the $s\leftrightarrow1/s$ degeneracy. However, because the projected separation is very close to $s=1$, the physical implications for the planet for the two solutions are quite similar, except for the value of $s$. By combining the light-curve parameters with a Galactic model, we have estimated the planet mass $M_2=17.9^{+9.6}_{-8.8}\,\mathrm{M}_\oplus$ and the lens distance $D_\mathrm{L}=6.7^{+1.0}_{-1.3}\,\mathrm{kpc}$, corresponding to a Neptune-mass planet close to the Galactic bulge. Such events with a low absolute latitude ($|b|\approx 1.1\,\mathrm{deg}$) are subject to both high extinction and more uncertain source distances, two factors that may affect the mass measurements in the provisional Wide Field Infrared Survey Telescope fields. More events are needed to investigate the potential trade-off between the higher lensing rate and the difficulty in measuring masses in these low-latitude fields.

We combine statistical arguments and dynamical analysis to study the orbital configuration of the microlensing planetary system OGLE-2012-BLG-0026L. This system is known to host two massive planets, with both projected close to the Einstein ring at the time of the detection. Assuming an isotropic distribution for the orbital orientation, we find that the two planets should also be closely spaced in three-dimensional space and that the ratio of their orbital periods is almost certainly less than two. With N-body numerical integrations, we then identify two types of stable configurations: the two planets can be in first-order mean-motion resonances (MMRs) and have significant ($\gtrsim0.1$) eccentricities, or they stay out of MMRs and have nearly circular orbits. The latter is disfavored, given the absence of similar planet pairs in radial velocity (RV) observations as well as the theoretical difficulties in forming such a configuration. Therefore, the two massive planets in OGLE-2012-BLG-0026L are likely in a resonance configuration. Our work shows that the microlensing technique, which usually only measures the projected configurations, can also probe the detailed dynamical state of multi-planet systems. We also discuss theoretical implications of measuring the multiplicity and the orbital architecture of cold planets.

Currently, we have only limited means to probe the presence of planets at large orbital separations. Foreman-Mackey et al. searched for long-period transiting planets in the Kepler light curves using an automated pipeline. Here, we apply their pipeline, with minor modifications, to a larger sample and use updated stellar parameters from Gaia DR2. The latter boosts the stellar radii for most of the planet candidates found by FM16, invalidating a number of them as false positives. We identify 15 candidates, including two new ones. All have sizes from 0.3 to 1 $R_{\rm J}$, and all but two have periods from 2 to 10 yr. We report two main findings based on this sample. First, the planet occurrence rate for the above size and period ranges is $0.70^{+0.40}_{-0.20}$ planets per Sun-like star, with the frequency of cold Jupiters agreeing with that from radial velocity surveys. Planet occurrence rises with decreasing planet size, roughly describable as $dN/d\log R \propto R^{\alpha}$ with $\alpha = -1.6^{+1.0}_{-0.9}$, i.e., Neptune-sized planets are some four times more common than Jupiter-sized ones. Second, five out of our 15 candidates orbit stars with known transiting planets at shorter periods, including one with five inner planets. We interpret this high incidence rate to mean: (1) almost all our candidates should be genuine; (2) across a large orbital range (from $\sim 0.05$ to a few astronomical units), mutual inclinations in these systems are at most a few degrees; and (3) large outer planets exist almost exclusively in systems with small inner planets.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

We report the discovery of a $Spitzer$ microlensing planet OGLE-2018-BLG-0596Lb, with preferred planet-host mass ratio $q \sim 2\times10^{-4}$. The planetary signal, which is characterized by a short $(\sim 1~{\rm day})$ "bump" on the rising side of the lensing light curve, was densely covered by ground-based surveys. We find that the signal can be explained by a bright source that fully envelops the planetary caustic, i.e., a "Hollywood" geometry. Combined with the source proper motion measured from $Gaia$, the $Spitzer$ satellite parallax measurement makes it possible to precisely constrain the lens physical parameters. The preferred solution, in which the planet perturbs the minor image due to lensing by the host, yields a Uranus-mass planet with a mass of $M_{\rm p} = 13.9\pm1.6~M_{\oplus}$ orbiting a mid M-dwarf with a mass of $M_{\rm h} = 0.23\pm0.03~M_{\odot}$. There is also a second possible solution that is substantially disfavored but cannot be ruled out, for which the planet perturbs the major image. The latter solution yields $M_{\rm p} = 1.2\pm0.2~M_{\oplus}$ and $M_{\rm h} = 0.15\pm0.02~M_{\odot}$. By combining the microlensing and $Gaia$ data together with a Galactic model, we find in either case that the lens lies on the near side of the Galactic bulge at a distance $D_{\rm L} \sim 6\pm1~{\rm kpc}$. Future adaptive optics observations may decisively resolve the major image/minor image degeneracy.

We report the discovery of the planet OGLE-2018-BLG-0532Lb, with very obvious signatures in the light curve that lead to an estimate of the planet-host mass ratio $q=M_{\rm planet}/M_{\rm host}\simeq 1\times10^{-4}$. Although there are no obvious systematic residuals to this double-lens/single-source (2L1S) fit, we find that $\chi^2$ can be significantly improved by adding either a third lens (3L1S, $\Delta\chi^2=81$) or second source (2L2S, $\Delta\chi^2=65$) to the lens-source geometry. After thorough investigation, we conclude that we cannot decisively distinguish between these two scenarios and therefore focus on the robustly-detected planet. However, given the possible presence of a second planet, we investigate to what degree and with what probability such additional planets may affect seemingly single-planet light curves. Our best estimates for the properties of the lens star and the secure planet are: a host mass $M\sim 0.25\,M_\odot$, system distance $D_L\sim 1\,$kpc and planet mass $m_{p,1}= 8\,M_\oplus$ with projected separation $a_{1,\perp}=1.4\,$au. However, there is a relatively bright $I=18.6$ (and also relatively blue) star projected within $<50\,$mas of the lens, and if future high-resolution images show that this is coincident with the lens, then it is possible that it is the lens, in which case, the lens would be both more massive and more distant than the best-estimated values above.

Using the Ly$\alpha$ mass assignment scheme (LyMAS), we make theoretical predictions for the 3-dimensional 3-point correlation function (3PCF) of the Ly$\alpha$ forest at redshift $z=2.3$. We bootstrap results from the (100 $h^{-1} \mbox{ Mpc}$)$^3$ Horizon hydrodynamic simulation to a (1 $h^{-1}$ Gpc)$^3$ $N$-body simulation, considering both a uniform UV background (UVB) and a fluctuating UVB sourced by quasars with a comoving $n_q \approx 10^{-5}$ $h^3$ Mpc$^{-3}$ placed either in massive halos or randomly. On scales of $10-30$ $h^{-1} \mbox{ Mpc}$, the flux 3PCF displays hierarchical scaling with the square of the 2PCF, but with an unusual value of $Q \equiv \zeta_{123}/(\xi_{12} \xi_{13} + \xi_{12} \xi_{23} + \xi_{13} \xi_{23}) \approx -4.5$ that reflects the low bias of the Ly$\alpha$ forest and the anti-correlation between mass density and transmitted flux. For halo-based quasars and an ionizing photon mean free path of $\lambda = 300$ $h^{-1} \mbox{ Mpc}$ comoving, UVB fluctuations moderately depress the 2PCF and 3PCF, with cancelling effects on $Q$. For $\lambda = 100$ $h^{-1} \mbox{ Mpc}$ or 50 $h^{-1} \mbox{ Mpc}$, UVB fluctuations substantially boost the 2PCF and 3PCF on large scales, shifting the hierarchical ratio to $Q \approx -3$. We scale our simulation results to derive rough estimate of the 3PCF detectability in observational data sets for the redshift range $z=2.1 - 2.6$. At $r = 10$ $h^{-1} \mbox{ Mpc}$ and 20 $h^{-1} \mbox{ Mpc}$, we predict a signal-to-noise (SNR) of $\sim$ 9 and $\sim$ 7, respectively, for both BOSS and eBOSS, and $\sim$ 37 and $\sim$ 25 for DESI. At $r = 40$ $h^{-1} \mbox{ Mpc}$ the predicted SNR is lower by $\sim$ 3$-$5 times. Measuring the flux 3PCF would be a novel test of the conventional paradigm of the Ly$\alpha$ forest and help separate the contributions of UVB fluctuations and density fluctuations to Ly$\alpha$ forest clustering.

We report the mass and distance measurements of two single-lens events from the 2017 Spitzer microlensing campaign. The ground-based observations yield the detection of finite-source effects, and the microlens parallaxes are derived from the joint analysis of ground-based observations and Spitzer observations. We find that the lens of OGLE-2017-BLG-1254 is a $0.60 \pm 0.03 M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.11~\text{kpc}$, where $D_{\rm LS}$ is the distance between the lens and the source. The second event, OGLE-2017-BLG-1161, is subject to the known satellite parallax degeneracy, and thus is either a $0.51^{+0.12}_{-0.10} M_{\odot}$ star with $D_{\rm LS} = 0.40 \pm 0.12~\text{kpc}$ or a $0.38^{+0.13}_{-0.12} M_{\odot}$ star with $D_{\rm LS} = 0.53 \pm 0.19~\text{kpc}$. Both of the lenses are therefore isolated stars in the Galactic bulge. By comparing the mass and distance distributions of the eight published Spitzer finite-source events with the expectations from a Galactic model, we find that the Spitzer sample is in agreement with the probability of finite-source effects occurrence in single lens events.

We study the dynamical evolution of the stellar mass binary black holes (BBHs) in a galactic nucleus that contains a massive black hole (MBH). For a comprehensive study of their merging events, we consider simultaneously the non-resonant and resonant relaxations of the BBHs, the binary-single encounters of the BBHs with the field stars, the Kozai-Lidov (KL) oscillation and the close encounters between the BBHs and the central MBH, which usually lead to binaries' tidal disruptions. As the BBHs are usually heavier than the background stars, they sink to the center by mass segregation, making the KL oscillation an important effect in merging BBHs. The binary-single encounters can not only lead to softening and ionization of the BBHs, it can also make them hardening, that increases the merging rates significantly. The mergers of BBHs are mainly contributed by galaxies containing MBHs less massive than $10^8 M_\odot$ and the total event rates are likely in orders of $1$--$10$ Gpc$^{-3}$ yr$^{-1}$, depending on the detailed assumptions of the nucleus clusters. About $3-10\%$ of these BBH mergers are with eccentricity $\ge 0.01$ when their gravitational wave oscillating frequencies enter the LIGO band ($10$\,Hz). Our results show that merging the BBHs within galactic nuclei can be an important source of the merging events detected by the Advanced LIGO/Virgo detectors, and they can be distinguished from BBH mergers from the galactic fields and globular clusters when enough events are accumulated.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.

With the largest dish Five-hundred-meter Aperture Spherical radio Telescope (FAST), both the mean and single pulses of PSR B2016$+$28, especially including the single-pulse structure, are investigated in detail in this study. The mean pulse profiles at different frequencies can be well fitted in a conal model, and the peak separation of intensity-dependent pulse profiles increases with intensity. The integrated pulses are obviously frequency dependent (pulse width decreases by $\sim\,20\%$ as frequency increases from 300 MHz to 750 MHz), but the structure of single pulses changes slightly (the corresponding correlation scale decreases by only $\sim\,1\%$). This disparity between mean and single pulses provides independent evidence for the existence of the RS-type vacuum inner gap, indicating a strong bond between particles on the pulsar surface. Diffused drifting sub-pulses are analyzed. The results show that the modulation period along pulse series ($P_3$) is positively correlated to the separation between two adjacent sub-pulses ($P_2$). This correlation may hint a rough surface on the pulsar, eventually resulting in the irregular drift of sparks. All the observational results may have significant implications in the dynamics of pulsar magnetosphere and are discussed extensively in this paper.

The Commensal Radio Astronomy Five-hundred-meter Aperture Spherical radio Telescope (FAST) Survey (CRAFTS) utilizes the novel drift-scan commensal survey mode of FAST and can generate billions of pulsar candidate signals. The human experts are not likely to thoroughly examine these signals, and various machine sorting methods are used to aid the classification of the FAST candidates. In this study, we propose a new ensemble classification system for pulsar candidates. This system denotes the further development of the pulsar image-based classification system (PICS), which was used in the Arecibo Telescope pulsar survey, and has been retrained and customized for the FAST drift-scan survey. In this study, we designed a residual network model comprising 15 layers to replace the convolutional neural networks (CNNs) in PICS. The results of this study demonstrate that the new model can sort >96% of real pulsars to belong the top 1% of all candidates and classify >1.6 million candidates per day using a dual--GPU and 24--core computer. This increased speed and efficiency can help to facilitate real-time or quasi-real-time processing of the pulsar-search data stream obtained from CRAFTS. In addition, we have published the labeled FAST data used in this study online, which can aid in the development of new deep learning techniques for performing pulsar searches.

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$.

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.

We analyze the binary microlensing event MOA-2016-BLG-231, which was observed from the ground and from Spitzer. The lens is composed of very low-mass brown dwarfs (BDs) with $M_1 = 21^{+12}_{-5} \ M_J$ and $M_2 = 9^{+5}_{-2}\ M_J$, and it is located in the Galactic disk $D_{\rm L} = 2.85^{+0.88}_{-0.50}\ {\rm kpc}$. This is the fifth binary brown dwarf discovered by microlensing, and the BD binary is moving counter to the orbital motion of disk stars. Constraints on the lens physical properties come from late time, non-caustic-crossing features of the Spitzer light curve. Thus, MOA-2016-BLG-231 shows how Spitzer plays a crucial role in resolving the nature of BDs in binary BD events with short timescale ($\lesssim 10$ days).

We study the influence of stellar metallicity on the fraction of stars with planets (i.e., the occurrence rate of planetary systems) and the average number of planets per star (i.e., the occurrence rate of planets). The former directly reveals the planet formation efficiency, whereas the latter reveals the final product of formation and evolution. We show that these two occurrence rates have different dependences on stellar metallicity. Specifically, the fraction of stars with planets rises gradually with metallicity, from $\sim$25% to $\sim$36% for 0.4 dex of [Fe/H] for all Kepler-like planets (period $P<400$ days and radius $R_{\rm p}\gtrsim R_\oplus$). The average number of planets per star reaches a plateau (or possibly starts declining) at [Fe/H]$\gtrsim0.1$. This is plausibly caused by the emergence of distant giant planets at high metallicities, given that the close-in small planets and the distant giants preferentially co-exist in the same system.

The kinematics of isolated brown dwarfs in the Galaxy, beyond the solar neighborhood, is virtually unknown. Microlensing has the potential to probe this hidden population, as it can measure both the mass and five of the six phase-space coordinates (all except the radial velocity) even of a dark isolated lens. However, the measurements of both the microlens parallax and finite-source effects are needed in order to recover the full information. Here, we combine $Spitzer$ satellite parallax measurement with the ground-based light curve, which exhibits strong finite-source effects, of event OGLE-2017-BLG-0896. We find two degenerate solutions for the lens (due to the known satellite-parallax degeneracy), which are consistent with each other except for their proper motion. The lens is an isolated brown dwarf with a mass of either $18\pm1M_J$ or $20\pm1M_J$. This is the lowest isolated-object mass measurement to date, only $\sim$45\% more massive than the theoretical deuterium-fusion boundary at solar metallicity, which is the common definition of a free-floating planet. The brown dwarf is located at either $3.9\pm0.1$ kpc or $4.1\pm0.1$ kpc toward the Galactic bulge, but with proper motion in the opposite direction of disk stars, with one solution suggesting it is moving within the Galactic plane. While it is possibly a halo brown dwarf, it might also represent a different, unknown population.

We report the discovery of a low mass-ratio planet $(q = 1.3\times10^{-4})$, i.e., 2.5 times higher than the Neptune/Sun ratio. The planetary system was discovered from the analysis of the KMT-2017-BLG-0165 microlensing event, which has an obvious short-term deviation from the underlying light curve produced by the host of the planet. Although the fit improvement with the microlens parallax effect is relatively low, one component of the parallax vector is strongly constrained from the light curve, making it possible to narrow down the uncertainties of the lens physical properties. A Bayesian analysis yields that the planet has a super-Neptune mass $(M_{2}=34_{-12}^{+15}~M_{\oplus})$ orbiting a Sun-like star $(M_{1}=0.76_{-0.27}^{+0.34}~M_{\odot})$ located at $4.5~{\rm kpc}$. The blended light is consistent with these host properties. The projected planet-host separation is $a_{\bot}={3.45_{-0.95}^{+0.98}}~{\rm AU}$, implying that the planet is located outside the snowline of the host, i.e., $a_{sl}\sim2.1~{\rm AU}$. KMT-2017-BLG-0165Lb is the sixteenth microlensing planet with mass ratio $q<3\times10^{-4}$. Using the fifteen of these planets with unambiguous mass-ratio measurements, we apply a likelihood analysis to investigate the form of the mass-ratio function in this regime. If we adopt a broken power law for the form of this function, then the break is at $q_{\rm br}\simeq0.55\times10^{-4}$, which is much lower than previously estimated. Moreover, the change of the power law slope, $\zeta>3.3$ is quite severe. Alternatively, the distribution is also suggestive of a "pile-up" of planets at Neptune-like mass ratios, below which there is a dramatic drop in frequency.

We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.

Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.

We have studied 4265 giant pulses (GPs) from the millisecond pulsar B1937+21; the largest-ever sample gathered for this pulsar, in observations made with the Large European Array for Pulsars. The pulse energy distribution of GPs associated with the interpulse are well-described by a power law, with index $\alpha = -3.99 \pm 0.04$, while those associated with the main pulse are best-described by a broken power law, with the break occurring at $\sim7$ Jy $\mu$s, with power law indices $\alpha_{\text{low}} = -3.48 \pm 0.04$ and $\alpha_{\text{high}} = -2.10 \pm 0.09$. The modulation indices of the GP emission are measured, which are found to vary by $\sim0.5$ at pulse phases close to the centre of the GP phase distributions. We find the frequency-resolved structure of GPs to vary significantly, and in a manner that cannot be attributed to the interstellar medium influence on the observed pulses. We examine the distribution of polarisation fractions of the GPs and find no correlation between GP emission phase and fractional polarisation. We use the GPs to time PSR B1937+21 and although the achievable time of arrival precision of the GPs is approximately a factor of two greater than that of the average pulse profile, there is a negligible difference in the precision of the overall timing solution when using the GPs.

We report the discovery of a giant planet in the KMT-2016-BLG-1397 microlensing event, which was found by The Korea Microlensing Telescope Network (KMTNet) alone. The time scale of this event is t_E = 40.0 +- 0.5 days and the mass ratio between the lens star and its companion is q = 0.016 +- 0.002. The planetary perturbation in the light curve is a smooth bump, resulting in the classical binary-lens/binary-source (2L1S/1L2S) degeneracy. We measure the V - I color of the (putative) two sources in the 1L2S model, and then effectively rule out the binary source solution. The finite-source effect is marginally detected. Combined with the limits on the blend flux and the probability distribution of the source size normalized by the Einstein radius rho, a Bayesian analysis yields the lens mass M_L = 0.45+0.33-0.28 M_sun, at distance of D_L = 6.60+1.10-1.30 kpc. Thus the companion is a super-Jupiter of a mass m_p = 7.0+5.2-4.3 M_J , at a projected separation r = 5.1+1.5-1.7 AU, indicating that the planet is well beyond the snow line of the host star.

A new type of gamma ray spectrum is predicted in a general hadronic framework if considering the gluon condensation effects in proton. The result presents the power-law with a sharp break in the gamma ray spectra at the TeV-band. We suggest to probe this GC-signature in the Earth limb gamma-ray spectra using the DArk Matter Particle Explorer and the CALorimetric Electron Telescope on orbit.

Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.

The gluon condensation (GC)-effects in high energy cosmic rays are investigated. After a brief review of the GC, several examples including gamma-, electron-, and positron-spectra in a broad GeV$\sim $TeV region can be explained by the GC-effects. We find that the GC may break the power-law of the cosmic ray spectra if the energy of accelerated protons exceeds the GC-threshold. The GC is a new phenomenon that is not yet known, it provides a new window to understand the high energy cosmic ray spectra.

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.

Here we draw attention to a candidate system with one hot Jupiter and one small, nearby companion, revealed by the recent Kepler data release (DR25). The hot Jupiter, Kepler-730b, has radius $R_{\rm p}=11.36^{+1.14}_{-0.98}~R_\oplus$ and orbital period $P=6.492$ d, and the newly discovered companion, KOI-929.02, has $R_{\rm p}=1.45^{+0.15}_{-0.20}~R_\oplus$ and $P=2.852$ d. No transit timing variation (TTV) was detected because of the marginal detection of the small companion transit, but this small companion passed all the validation tests. This system is probably another planetary system with hot Jupiter and small, nearby companion, after the outstanding WASP-47 system, and so far the only such system (out of 46) in the prime Kepler mission. The nature of this system, if confirmed, would suggest that hot Jupiters with small, nearby companions are probably more common than we used to believe. There remains a possibility that the small companion actually transits a star that is different from the hot Jupiter host, and follow-up observations with 10-m telescopes can potentially resolve this issue.

We develop stellar population models to predict the number of binaries with a single luminous member is Gaia and Hipparcos. Our models yield dozens of detections of black hole - luminous companion binaries (BHLC) and hundreds to thousands of neutron star - luminous companion binaries (NSLC) with Gaia. Interestingly, our models also yield a single detection of BHLC binary with Hipparcos, and a few NSLC binaries. We also show how the statistical distribution of detected binaries with a single luminous companion can be used to constrain the formation process of neutron stars and black holes.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.

The scattering of fast radio bursts (FRBs) by the intergalactic medium (IGM) is explored using cosmological hydrodynamical simulations. We confirm that the scattering by the clumpy IGM has significant line-of-sight variations. We demonstrate that the scattering by the IGM in the voids and walls of the cosmic web is weak, but it can be significantly enhanced by the gas in clusters and filaments. The observed non-monotonic dependence of the FRB widths on the dispersion measures (DM) cannot determine whether the IGM is an important scattering matter or not. The IGM may dominate the scattering of some FRBs, and the host galaxy dominates others. For the former case, the scattering should be primarily caused by the medium in clusters. A mock sample of 500 sources shows that $\tau_{\rm{IGM}} \propto \rm{DM_{IGM}}^{1.6-2.1}$ at $z<1.5$. Assuming that the turbulence follows Kolmogorov scaling, we find that an outer scale of $L_0\sim 5\,$pc is required to make $\tau_{\rm{IGM}} \sim 1-10\,$ms at $\nu=1\, $ GHz. The required $L_0\sim 5\, $pc can alleviate the tension in the timescales of turbulent heating and cooling but is still $\sim 4$ orders of magnitude lower than the presumed injection scale of turbulence in the IGM. The gap is expected to be effectively shortened if the simulation resolution is further increased. The mechanisms that may further reduce the gap are shortly discussed. If future observations can justify the role of the IGM in the broadening of FRBs, it can help to probe the gas in clusters and filaments.

We present a newly implemented single-pulse pipeline for the PALFA survey to efficiently identify single radio pulses from pulsars, Rotating Radio Transients (RRATs) and Fast Radio Bursts (FRBs). We have conducted a sensitivity analysis of this new pipeline in which multiple single pulses with a wide range of parameters were injected into PALFA data sets and run through the pipeline. Based on the recovered pulses, we find that for pulse widths $\rm < 5\ ms$ the sensitivity of the PALFA pipeline is at most a factor of $\rm \sim 2$ less sensitive to single pulses than our theoretical predictions. For pulse widths $\rm > 10\ ms$, as the $\rm DM$ decreases, the degradation in sensitivity gets worse and can increase up to a factor of $\rm \sim 4.5$. Using this pipeline, we have thus far discovered 7 pulsars and 2 RRATs and identified 3 candidate RRATs and 1 candidate FRB. The confirmed pulsars and RRATs have DMs ranging from 133 to 386 pc cm$^{-3}$ and flux densities ranging from 20 to 160 mJy. The pulsar periods range from 0.4 to 2.1 s. We report on candidate FRB 141113, which we argue is likely astrophysical and extragalactic, having $\rm DM \simeq 400\ pc~cm^{-3}$, which represents an excess over the Galactic maximum along this line of sight of $\rm \sim$ 100 - 200 pc cm$^{-3}$. We consider implications for the FRB population and show via simulations that if FRB 141113 is real and extragalactic, the slope $\alpha$ of the distribution of integral source counts as a function of flux density ($N (>S) \propto S^{-\alpha}$) is $1.4 \pm 0.5$ (95% confidence range). However this conclusion is dependent on several assumptions that require verification.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

We report the discovery and the analysis of the short timescale binary-lens microlensing event, MOA-2015-BLG-337. The lens system could be a planetary system with a very low mass host, around the brown dwarf/planetary mass boundary, or a brown dwarf binary. We found two competing models that explain the observed light curves with companion/host mass ratios of q~0.01 and ~0.17, respectively. From the measurement of finite source effects in the best-fit planetary model, we find a relatively small angular Einstein radius of theta_E ~ 0.03 mas which favors a low mass lens. We conduct a Bayesian analysis to obtain the probability distribution of the lens properties. The results for the planetary models strongly depend on the minimum mass, M_min, in the assumed mass function. In summary, there are two solutions of the lens system: (1) a brown dwarf/planetary mass boundary object orbited by a super-Neptune (the planetary model with M_min=0.001 M_sun) and (2) a brown dwarf binary (the binary model). If the planetary models is correct, this system can be one of a new class of planetary system, having a low host mass and also a planetary mass ratio (q <0.03) between the host and its companion. The discovery of the event is important for the study of planetary formation in very low mass objects. In addition, it is important to consider all viable solutions in these kinds of ambiguous events in order for the future comprehensive statistical analyses of planetary/binary microlensing events.

The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.

We report relations between inner ($<1$ au) super Earths (planets with mass/radius between Earth and Neptune) and outer ($>1$ au) giant planets (mass $>0.3~M_{\rm J}$, or cold Jupiters) around Sun-like stars, based on data from both ground-based radial velocity (RV) observations and the Kepler mission. We find that cold Jupiters appear three times more often around hosts of super Earths than they do around field stars. Given the prevalence of the super Earth systems, their cold Jupiters can account for nearly all cold Jupiters. In other words, cold Jupiters are almost certainly ($\sim90\%$) companied by super Earths. A few corollaries follow: (1) around metal-rich ([Fe/H]$>0.1$) stars, the fraction of super Earths with cold Jupiters can rise to $60\%$ or higher; (2) the inner architecture can be strongly impacted by the outer giant and we report some observational evidence for this; (3) planetary systems like our own, with cold Jupiters but no super Earths, should be rare ($\sim1\%$). The strong correlation between super Earths and cold Jupiters establish that super Earths and cold Jupiters do not compete for solid material, rather, they share similar origins, with the cold Jupiter formation requiring a somewhat more stringent condition. Lastly, we propose a few immediate observational tests of our results, using ground-based RV observations and ongoing/planned space missions.

We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.

We present the analysis of the binary-lens microlensing event OGLE-2017-BLG-0537. The light curve of the event exhibits two strong caustic-crossing spikes among which the second caustic crossing was resolved by high-cadence surveys. It is found that the lens components with a mass ratio $\sim 0.5$ are separated in projection by $\sim 1.3\thetae$, where $\thetae$ is the angular Einstein radius. Analysis of the caustic-crossing part yields $\thetae=1.77\pm 0.16$~mas and a lens-source relative proper motion of $\mu =12.4\pm 1.1~{\rm mas}~{\rm yr}^{-1}$. The measured $\mu$ is the third highest value among the events with measured proper motions and $\sim 3$ times higher than the value of typical Galactic bulge events, making the event a strong candidate for follow-up observations to directly image the lens by separating it from the source. From the angular Einstein radius combined with the microlens parallax, it is estimated that the lens is composed of two main-sequence stars with masses $M_1\sim 0.4~M_\odot$ and $M_2\sim 0.2~M_\odot$ located at a distance of $D_{\rm L}\sim 1.2$~kpc. However, the physical lens parameters are not very secure due to the weak microlens-parallax signal, and thus we cross check the parameters by conducting a Bayesian analysis based on the measured Einstein radius and event timescale combined with the blending constraint. From this, we find that the physical parameters estimated from the Bayesian analysis are consistent with those based on the measured microlens parallax. Resolving the lens from the source can be done in about 5 years from high-resolution follow-up observations and this will provide a rare opportunity to test and refine the microlensing model.

Microlensing is a powerful and unique technique to probe isolated objects in the Galaxy. To study the characteristics of these interesting objects based on the microlensing method, measurement of the microlens parallax is required to determine the properties of the lens. Of the various methods to measure microlens parallax, the most robust way is to make simultaneous ground- and space-based observations, i.e., by measuring the space-based microlens parallax. However, space-based campaigns usually require "expensive" resources. Gould & Yee (2012) proposed an idea called the "cheap space-based microlens parallax" that can measure the lens-parallax using only $2$ or $3$ space-based observations of high-magnification events. This cost-effective observation strategy to measure microlens parallaxes could be used by space-borne telescopes to build a complete sample for studying isolated objects. This would enable a direct measurement of the mass function including both extremely low-mass objects and high-mass stellar remnants. However, to adopt this idea requires a test to check how it would work in actual situations. Thus, we present the first practical test of the idea using the high-magnification microlensing event OGLE-2016-BLG-1045, for which a subset of Spitzer observations fortuitously duplicate the prescription of Gould & Yee (2012). From the test, we confirm that the measurement of the lens-parallax adopting this idea has sufficient accuracy to determine the physical properties of the isolated lens.

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, cosmic strings, and a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using recent simulations, yielding an SSE-marginalized $95\%$ upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.

The Pulsar Arecibo L-Band Feed Array (PALFA) survey, the most sensitive blind search for radio pulsars yet conducted, is ongoing at the Arecibo Observatory in Puerto Rico. The vast majority of the 180 pulsars discovered by PALFA have spin periods shorter than 2 seconds. Pulsar surveys may miss long-period radio pulsars due to the summing of a finite number of harmonic components in conventional Fourier analyses (typically $\sim$16), or due to the strong effect of red noise at low modulation frequencies. We address this reduction in sensitivity by using a time-domain search technique: the Fast-Folding Algorithm (FFA). We designed a program that implements a FFA-based search in the PALFA processing pipeline, and tested the efficiency of the algorithm by performing tests under both ideal, white noise conditions, as well as with real PALFA observational data. In the two scenarios, we show that the time-domain algorithm has the ability to outperform the FFT-based periodicity search implemented in the survey. We perform simulations to compare the previously reported PALFA sensitivity with that obtained using our new FFA implementation. These simulations show that for a pulsar having a pulse duty cycle of roughly 3%, the performance of our FFA pipeline exceeds that of our FFT pipeline for pulses with DM $\lesssim$ 40 pc cm$^{-3}$ and for periods as short as $\sim$500 ms, and that the survey sensitivity is improved by at least a factor of two for periods $\gtrsim$ 6 sec. Discoveries from the implementation of the algorithm in PALFA are also presented in this paper.

We constrain the intrinsic architecture of Kepler planetary systems by modeling the observed multiplicities of the transiting planets (tranets) and their transit timing variations (TTVs). We robustly determine that the fraction of Sun-like stars with Kepler-like planets, $\eta_{\rm Kepler}$, is $30\pm3\%$. Here Kepler-like planets are planets that have radii $R_{\rm p} \gtrsim R_\oplus$ and orbital periods $P<400$~days. Our result thus significantly revises previous claims that more than 50\% of Sun-like stars have such planets. Combining with the average number of Kepler planets per star ($\sim0.9$), we obtain that on average each planetary system has $3.0\pm0.3$ planets within 400 days. We also find that the dispersion in orbital inclinations of planets within a given planetary system, $\sigma_{i,k}$, is a steep function of its number of planets, $k$. This can be parameterized as $\sigma_{i,k}\propto k^\alpha$ and we find that $-4<\alpha<-2$ at 2-$\sigma$ level. Such a distribution well describes the observed multiplicities of both tranets and TTVs with no excess of single tranets. Therefore we do not find evidence supporting the so-called "Kepler dichotomy." Together with a previous study on orbital eccentricities, we now have a consistent picture: the fewer planets in a system, the hotter it is dynamically. We discuss briefly possible scenarios that lead to such a trend. Despite our Solar system not belonging to the Kepler club, it is interesting to notice that the Solar system also has three planets within 400 days and that the inclination dispersion is similar to Kepler systems of the same multiplicity.

The gas is the dominant component of baryonic matter in most galaxy groups and clusters. The spatial offsets of gas centre from the halo centre could be an indicator of the dynamical state of cluster. Knowledge of such offsets is important for estimate the uncertainties when using clusters as cosmological probes. In this paper, we study the centre offsets $r_{\rm off}$ between the gas and that of all the matter within halo systems in $\Lambda$CDM cosmological hydrodynamic simulations. We focus on two kinds of centre offsets: one is the three-dimensional PB offsets between the gravitational potential minimum of the entire halo and the barycentre of the ICM, and the other is the two-dimensional PX offsets between the potential minimum of the halo and the iterative centroid of the projected synthetic X-ray emission of the halo. Halos at higher redshifts tend to have larger values of rescaled offsets $r_{\rm off}/r_{200}$ and larger gas velocity dispersion $\sigma_{v}^{\rm gas}/\sigma_{200}$. For both types of offsets, we find that the correlation between the rescaled centre offsets $r_{\rm off}/r_{200}$ and the rescaled 3D gas velocity dispersion, $\sigma_v^{\rm gas}/\sigma_{200}$ can be approximately described by a quadratic function as $r_{{\rm off}}/r_{200} \propto (\sigma_{v}^{\rm gas}/\sigma_{200} - k_2)^{2}$. A Bayesian analysis with MCMC method is employed to estimate the model parameters. Dependence of the correlation relation on redshifts and the gas mass fraction are also investigated.

We show that dense OGLE and KMTNet $I$-band survey data require four bodies (sources plus lenses) to explain the microlensing light curve of OGLE-2015-BLG-1459. However, these can equally well consist of three lenses and one source (3L1S), two lenses and two sources (2L2S) or one lens and three sources (1L3S). In the 3L1S and 2L2S interpretations, the host is a brown dwarf and the dominant companion is a Neptune-class planet, with the third body (in the 3L1S case) being a Mars-class object that could have been a moon of the planet. In the 1L3S solution, the light curve anomalies are explained by a tight (five stellar radii) low-luminosity binary source that is offset from the principal source of the event by $\sim 0.17\,\au$. These degeneracies are resolved in favor of the 1L3S solution by color effects derived from comparison to MOA data, which are taken in a slightly different ($R/I$) passband. To enable current and future ($WFIRST$) surveys to routinely characterize exomoons and distinguish among such exotic systems requires an observing strategy that includes both a cadence faster than 9 min$^{-1}$ and observations in a second band on a similar timescale.

We analyze the combined Spitzer and ground-based data for OGLE-2017-BLG-1140 and show that the event was generated by a Jupiter-class $(m_p\simeq 1.6\,M_{\rm jup})$ planet orbiting a mid-late M dwarf $(M\simeq 0.2\,M_\odot)$ that lies $D_{LS}\simeq 1.0\,\mathrm{kpc}$ in the foreground of the microlensed, Galactic-bar, source star. The planet-host projected separation is $a_\perp \simeq 1.0\,\mathrm{au}$, i.e., well-beyond the snow line. By measuring the source proper motion ${\mathbf{\mu}}_s$ from ongoing, long-term OGLE imaging, and combining this with the lens-source relative proper motion ${\mathbf{\mu}}_\mathrm{rel}$ derived from the microlensing solution, we show that the lens proper motion ${\mathbf{\mu}}_l={\mathbf{\mu}}_\mathrm{rel} + {\mathbf{\mu}}_s$ is consistent with the lens lying in the Galactic disk, although a bulge lens is not ruled out. We show that while the Spitzer and ground-based data are comparably well fitted by planetary (i.e., binary-lens, 2L1S) models and by binary-source (1L2S) models, the combination of Spitzer and ground-based data decisively favor the planetary model. This is a new channel to resolve the 2L1S/1L2S degeneracy, which can be difficult to break in some cases.

We present Korea Microlensing Telescope Network (KMTNet) light curves for microlensing-event candidates for the 2016 season, which covers an area of 97 sq.deg observed at cadences ranging from Gamma=0.2/hr to Gamma=8/hr from three southern sites in Chile, South Africa, and Australia. These 2163 light curves are comprised of 1856 "clear microlensing" and 307 "possible microlensing" events (including 265 previously released from the K2 C9 field). The data policy is very similar to the one governing the 2015 release. The changes relative to 2015 in the algorithms to find and vet microlensing events are comprehensively described.

We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.

We report the discovery of microlensing planet OGLE-2017-BLG-0373Lb. We show that while the planet-host system has an unambiguous microlens topology, there are two geometries within this topology that fit the data equally well, which leads to a factor 2.5 difference in planet-host mass ratio, i.e., $q=1.5\times 10^{-3}$ vs. $q=0.6\times 10^{-3}$. We show that this is an "accidental degeneracy" in the sense that it is due to a gap in the data. We dub it "the caustic-chirality degeneracy". We trace the mathematical origins of this degeneracy, which should enable similar degenerate solutions to be easily located in the future. A Bayesian estimate, based on a Galactic model, yields a host mass $M=0.25^{+0.30}_{-0.15} M_\odot$ at a distance $D_L=5.9^{+1.3}_{-1.95}$ kpc. The lens-source relative proper motion is relatively fast, $\mu=9$ mas/yr, which implies that the host mass and distance can be determined by high-resolution imaging after about 10 years. The same observations could in principle resolve the discrete degeneracy in $q$, but this will be more challenging.

We report the discovery of a planetary system in which a super-earth orbits a late M-dwarf host. The planetary system was found from the analysis of the microlensing event OGLE-2017-BLG-0482, wherein the planet signal appears as a short-term anomaly to the smooth lensing light curve produced by the host. Despite its weak signal and short duration, the planetary signal was firmly detected from the dense and continuous coverage by three microlensing surveys. We find a planet/host mass ratio of $q\sim 1.4\times 10^{-4}$. We measure the microlens parallax $\pi_{\rm E}$ from the long-term deviation in the observed lensing light curve, but the angular Einstein radius $\theta_{\rm E}$ cannot be measured because the source trajectory did not cross the planet-induced caustic. Using the measured event timescale and the microlens parallax, we find that the masses of the planet and the host are $M_{\rm p}=9.0_{-4.5}^{+9.0}\ M_\oplus$ and $M_{\rm host}=0.20_{-0.10}^{+0.20}\ M_\odot$, respectively, and the projected separation between them is $a_\perp=1.8_{-0.7}^{+0.6}$ au. The estimated distance to the lens is $D_{\rm L}=5.8_{-2.1}^{+1.8}$ kpc. The discovery of the planetary system demonstrates that microlensing provides an important method to detect low-mass planets orbiting low-mass stars.

K2 Campaign 9 (K2C9) was the first space-based microlensing parallax survey capable of measuring microlensing parallaxes of free-floating planet candidate microlensing events. Simultaneous to K2C9 observations we conducted the K2C9 Canada-France-Hawaii Telescope Multi-Color Microlensing Survey (K2C9-CFHT MCMS) in order to measure the colors of microlensing source stars to improve the accuracy of K2C9's parallax measurements. We describe the difference imaging photometry analysis of the K2C9-CFHT MCMS observations, and present the project's first data release. This includes instrumental difference flux lightcurves of 217 microlensing events identified by other microlensing surveys, reference image photometry calibrated to PanSTARRS data release 1 photometry, and tools to convert between instrumental and calibrated flux scales. We derive accurate analytic transformations between the PanSTARRS bandpasses and the Kepler bandpass, as well as angular diameter-color relations in the PanSTARRS bandpasses. To demonstrate the use of our data set, we analyze ground-based and K2 data of a short timescale microlensing event, OGLE-2016-BLG-0795. We find the event has a timescale $t_{\rm E}=4.5 \pm 0.1$~days and microlens parallax $\pi_{\rm E}=0.12 \pm 0.03$ or $0.97 \pm 0.04$, subject to the standard satellite parallax degeneracy. We argue that the smaller value of the parallax is more likely, which implies that the lens is likely a stellar-mass object in the Galactic bulge as opposed to a super-Jupiter mass object in the Galactic disk.

We report the discovery of a giant planet in the OGLE-2017-BLG-1522 microlensing event. The planetary perturbations were clearly identified by high-cadence survey experiments despite the relatively short event timescale of $t_{\rm E} \sim 7.5$ days. The Einstein radius is unusually small, $\theta_{\rm E} = 0.065\,$mas, implying that the lens system either has very low mass or lies much closer to the microlensed source than the Sun, or both. A Bayesian analysis yields component masses $(M_{\rm host}, M_{\rm planet})=(46_{-25}^{+79}, 0.75_{-0.40}^{+1.26})~M_{\rm J}$ and source-lens distance $D_{\rm LS} = 0.99_{-0.54}^{+0.91}~{\rm kpc}$, implying that this is a brown-dwarf/Jupiter system that probably lies in the Galactic bulge, a location that is also consistent with the relatively low lens-source relative proper motion $\mu = 3.2 \pm 0.5~{\rm mas}~{\rm yr^{-1}}$. The projected companion-host separation is $0.59_{-0.11}^{+0.12}~{\rm AU}$, indicating that the planet is placed beyond the snow line of the host, i.e., $a_{sl} \sim 0.12~{\rm AU}$. Planet formation scenarios combined with the small companion-host mass ratio $q \sim 0.016$ and separation suggest that the companion could be the first discovery of a giant planet that formed in a protoplanetary disk around a brown dwarf host.

We present the analysis of KMT-2016-BLG-0212, a low flux-variation $(I_{\rm flux-var}\sim 20$) microlensing event, which is well-covered by high-cadence data from the three Korea Microlensing Telescope Network (KMTNet) telescopes. The event shows a short anomaly that is incompletely covered due to the brief visibility intervals that characterize the early microlensing season when the anomaly occurred. We show that the data are consistent with two classes of solutions, characterized respectively by low-mass brown-dwarf $(q=0.037)$ and sub-Neptune $(q<10^{-4})$ companions, respectively. Future high-resolution imaging should easily distinguish between these solutions.

We analyze the binary gravitational microlensing event OGLE-2017-BLG-1130 (mass ratio q~0.45), the first published case in which the binary anomaly was only detected by the Spitzer Space Telescope. This event provides strong evidence that some binary signals can be missed by observations from the ground alone but detected by Spitzer. We therefore invert the normal procedure, first finding the lens parameters by fitting the space-based data and then measuring the microlensing parallax using ground-based observations. We also show that the normal four-fold space-based degeneracy in the single-lens case can become a weak eight-fold degeneracy in binary-lens events. Although this degeneracy is resolved in event OGLE-2017-BLG-1130, it might persist in other events.

Mass measurements of gravitational microlenses require one to determine the microlens parallax $\pie$, but precise $\pie$ measurement, in many cases, is hampered due to the subtlety of the microlens-parallax signal combined with the difficulty of distinguishing the signal from those induced by other higher-order effects. In this work, we present the analysis of the binary-lens event OGLE-2017-BLG-0329, for which $\pie$ is measured with a dramatically improved precision using additional data from space-based $Spitzer$ observations. We find that while the parallax model based on the ground-based data cannot be distinguished from a zero-$\pie$ model at 2$\sigma$ level, the addition of the $Spitzer$ data enables us to identify 2 classes of solutions, each composed of a pair of solutions according to the well-known ecliptic degeneracy. It is found that the space-based data reduce the measurement uncertainties of the north and east components of the microlens-parallax vector $\pivec_{\rm E}$ by factors $\sim 18$ and $\sim 4$, respectively. With the measured microlens parallax combined with the angular Einstein radius measured from the resolved caustic crossings, we find that the lens is composed of a binary with components masses of either $(M_1,M_2)\sim (1.1,0.8)\ M_\odot$ or $\sim (0.4,0.3)\ M_\odot$ according to the two solution classes. The first solution is significantly favored but the second cannot be securely ruled out based on the microlensing data alone. However, the degeneracy can be resolved from adaptive optics observations taken $\sim 10$ years after the event.

Symmetries play an important role in modern theories of gravity. The strong equivalence principle (SEP) constitutes a collection of gravitational symmetries which are all implemented by general relativity. Alternative theories, however, are generally expected to violate some aspects of SEP. We test three aspects of SEP using observed change rates in the orbital period and eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's constancy as part of locational invariance of gravitation; 2. the post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance; 3. the universality of free fall (UFF) for strongly self-gravitating bodies. Based on the pulsar timing result of the combined dataset from the North American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|< 0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 < 4\times10^{-20}$ at 95% confidence. These are the first direct UFF and $\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various limitations of previous tests.

We report the discovery, via the microlensing method, of a new very-low-mass binary system. By combining measurements from Earth and from the Spitzer telescope in Earth-trailing orbit, we are able to measure the microlensing parallax of the event, and find that the lens likely consists of an $(12.0 \pm 0.6) M_{\rm J}$ + $(15.7 \pm 1.5) M_{\rm J}$ super-Jupiter / brown-dwarf pair. The binary is located at a distance of $(3.08 \pm 0.18)$ kpc in the Galactic Plane, and the components have a projected separation of $(0.43 \pm 0.03)$ AU. Two alternative solutions with much lower likelihoods are also discussed, an 8- and 6-$M_{\rm J}$ model and a 90- and 70-$M_{\rm J}$ model. Although disfavored at the 3-$\sigma$ and 5-$\sigma$ levels, these alternatives cannot be rejected entirely. We show how the more-massive of these models could be tested with future direct imaging.

Current microlensing surveys are sensitive to free-floating planets down to Earth-mass objects. All published microlensing events attributed to unbound planets were identified based on their short timescale (below two days), but lacked an angular Einstein radius measurement (and hence lacked a significant constraint on the lens mass). Here, we present the discovery of a Neptune-mass free-floating planet candidate in the ultrashort ($t_{\rm E}=0.320\pm0.003$ days) microlensing event OGLE-2016-BLG-1540. The event exhibited strong finite-source effects, which allowed us to measure its angular Einstein radius of $\theta_{\rm E}=9.2\pm0.5\,\mu$as. There remains, however, a degeneracy between the lens mass and distance. The combination of the source proper motion and source-lens relative proper motion measurements favors a Neptune-mass lens located in the Galactic disk. However, we cannot rule out that the lens is a Saturn-mass object belonging to the bulge population. We exclude stellar companions up to 15 au.

We report the discovery of a cold Super-Earth planet (m_p=4.4 +/- 0.5 M_Earth) orbiting a low-mass (M=0.23 +/- 0.03 M_Sun) M dwarf at projected separation a_perp = 1.18 +/- 0.10 AU, i.e., about 1.9 times the snow line. The system is quite nearby for a microlensing planet, D_Lens = 0.86 +/- 0.09 kpc. Indeed, it was the large lens-source relative parallax pi_rel=1.0 mas (combined with the low mass M) that gave rise to the large, and thus well-measured, "microlens parallax" that enabled these precise measurements. OGLE-2017-BLG-1434Lb is the eighth microlensing planet with planet-host mass ratio q < 1 * 10^-4. We apply a new planet-detection sensitivity method, which is a variant of "V/V_max", to seven of these eight planets to derive the mass-ratio function in this regime. We find dN/d(ln q) ~ q^p, with p = 1.05 (+0.78,-0.68), which confirms the "turnover" in the mass function found by Suzuki et al. relative to the power law of opposite sign n = -0.93 +/- 0.13 at higher mass ratios q >~ 2 * 10^-4. We combine our result with that of Suzuki et al. to obtain p = 0.73 (+0.42,-0.34).

We report the discovery and initial follow-up of a double neutron star (DNS) system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA) survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\, 0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds using a new phase binning mode. We have searched multiwavelength catalogs for coincident sources but did not find any counterparts. The improved position enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \, = \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates that this pulsar has been recycled. This and the orbital eccentricity lead to the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and the shortest estimated merger timescale, 46 Myr; at that time it will display the largest spin effects on gravitational wave waveforms of any such system discovered to date. We have measured the advance of periastron passage for this system, $\dot{\omega} \, = \, 25.6 \, \pm \, 0.3\, \deg \rm yr^{-1}$, implying a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the lowest mass DNS systems. This total mass measurement combined with the minimum companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$.

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models.

We present Korea Microlensing Telescope Network (KMTNet) light curves for microlensing-event candidates in the Kepler K2 C9 field having peaks within 3 effective timescales of the Kepler observations. These include 181 "clear microlensing" and 84 "possible microlensing" events found by the KMTNet event finder, plus 56 other events found by OGLE and/or MOA that were not found by KMTNet. All data for the first two classes are immediately available for public use without restriction.

We report the discovery of a sub-Jupiter mass planet orbiting beyond the snow line of an M-dwarf most likely in the Galactic disk as part of the joint Spitzer and ground-based monitoring of microlensing planetary anomalies toward the Galactic bulge. The microlensing parameters are strongly constrained by the light curve modeling and in particular by the Spitzer-based measurement of the microlens parallax, $\pi_\mathrm{E}$. However, in contrast to many planetary microlensing events, there are no caustic crossings, so the angular Einstein radius, $\theta_\mathrm{E}$ has only an upper limit based on the light curve modeling alone. Additionally, the analysis leads us to identify 8 degenerate configurations: the four-fold microlensing parallax degeneracy being doubled by a degeneracy in the caustic structure present at the level of the ground-based solutions. To pinpoint the physical parameters, and at the same time to break the parallax degeneracy, we make use of a series of arguments: the $\chi^2$ hierarchy, the Rich argument, and a prior Galactic model. The preferred configuration is for a host at $D_L=3.73_{-0.67}^{+0.66}~\mathrm{kpc}$ with mass $M_\mathrm{L}=0.30_{-0.12}^{+0.15}~\mathrm{M_\odot}$, orbited by a Saturn-like planet with $M_\mathrm{planet}=0.43_{-0.17}^{+0.21}~\mathrm{M_\mathrm{Jup}}$ at projected separation $a_\perp = 1.70_{-0.39}^{+0.38}~\mathrm{au}$, about 2.1 times beyond the system snow line. Therefore, it adds to the growing population of sub-Jupiter planets orbiting near or beyond the snow line of M-dwarfs discovered by microlensing. Based on the rules of the real-time protocol for the selection of events to be followed up with Spitzer, this planet will not enter the sample for measuring the Galactic distribution of planets.

We present the analysis of the binary-microlensing event OGLE-2014-BLG-0289. The event light curve exhibits very unusual five peaks where four peaks were produced by caustic crossings and the other peak was produced by a cusp approach. It is found that the quintuple-peak features of the light curve provide tight constraints on the source trajectory, enabling us to precisely and accurately measure the microlensing parallax $\pi_{\rm E}$. Furthermore, the three resolved caustics allow us to measure the angular Einstein radius $\thetae$. From the combination of $\pi_{\rm E}$ and $\thetae$, the physical lens parameters are uniquely determined. It is found that the lens is a binary composed of two M dwarfs with masses $M_1 = 0.52 \pm 0.04\ M_\odot$ and $M_2=0.42 \pm 0.03\ M_\odot$ separated in projection by $a_\perp = 6.4 \pm 0.5$ au. The lens is located in the disk with a distance of $D_{\rm L} = 3.3 \pm 0.3$~kpc. It turns out that the reason for the absence of a lensing signal in the {\it Spitzer} data is that the time of observation corresponds to the flat region of the light curve.

We present the first results from an optical reverberation mapping campaign executed in 2014, targeting the active galactic nuclei (AGN) MCG+08-11-011, NGC 2617, NGC 4051, 3C 382, and Mrk 374. Our targets have diverse and interesting observational properties, including a "changing look" AGN and a broad-line radio galaxy. Based on continuum-H$\beta$ lags, we measure black hole masses for all five targets. We also obtain H$\gamma$ and He{\sc ii}\,$\lambda 4686$ lags for all objects except 3C 382. The He{\sc ii}\,$\lambda 4686$ lags indicate radial stratification of the BLR, and the masses derived from different emission lines are in general agreement. The relative responsivities of these lines are also in qualitative agreement with photoionization models. These spectra have extremely high signal-to-noise ratios (100--300 per pixel) and there are excellent prospects for obtaining velocity-resolved reverberation signatures.

We present microlensing events in the 2015 Korea Microlensing Telescope Network (KMTNet) data and our procedure for identifying these events. In particular, candidates were detected with a novel "completed event" microlensing event-finder algorithm. The algorithm works by making linear fits to a (t0,teff,u0) grid of point-lens microlensing models. This approach is rendered computationally efficient by restricting u0 to just two values (0 and 1), which we show is quite adequate. The implementation presented here is specifically tailored to the commission-year character of the 2015 data, but the algorithm is quite general and has already been applied to a completely different (non-KMTNet) data set. We outline expected improvements for 2016 and future KMTNet data. The light curves of the 660 "clear microlensing" and 182 "possible microlensing" events that were found in 2015 are presented along with our policy for their public release.

The Large European Array for Pulsars (LEAP) is a European Pulsar Timing Array project that combines the Lovell, Effelsberg, Nan\c{c}ay, Sardinia, and Westerbork radio telescopes into a single tied-array, and makes monthly observations of a set of millisecond pulsars (MSPs). The overview of our experiment is presented in Bassa et al. (2016). Baseband data are recorded at a central frequency of 1396 MHz and a bandwidth of 128 MHz at each telescope, and are correlated offline on a cluster at Jodrell Bank Observatory using a purpose-built correlator, detailed in Smits et al. (2017). LEAP offers a substantial increase in sensitivity over that of the individual telescopes, and can operate in timing and imaging modes (notably in observations of the galactic centre radio magnetar; Wucknitz 2015). To date, 4 years of observations have been reduced. Here, we report on the scientific projects which have made use of LEAP data.

We report the discovery of OGLE-2016-BLG-1190Lb, which is likely to be the first Spitzer microlensing planet in the Galactic bulge/bar, an assignation that can be confirmed by two epochs of high-resolution imaging of the combined source-lens baseline object. The planet's mass M_p= 13.4+-0.9 M_J places it right at the deuterium burning limit, i.e., the conventional boundary between "planets" and "brown dwarfs". Its existence raises the question of whether such objects are really "planets" (formed within the disks of their hosts) or "failed stars" (low mass objects formed by gas fragmentation). This question may ultimately be addressed by comparing disk and bulge/bar planets, which is a goal of the Spitzer microlens program. The host is a G dwarf M_host = 0.89+-0.07 M_sun and the planet has a semi-major axis a~2.0 AU. We use Kepler K2 Campaign 9 microlensing data to break the lens-mass degeneracy that generically impacts parallax solutions from Earth-Spitzer observations alone, which is the first successful application of this approach. The microlensing data, derived primarily from near-continuous, ultra-dense survey observations from OGLE, MOA, and three KMTNet telescopes, contain more orbital information than for any previous microlensing planet, but not quite enough to accurately specify the full orbit. However, these data do permit the first rigorous test of microlensing orbital-motion measurements, which are typically derived from data taken over <1% of an orbital period.

We present microlensing planet OGLE-2017-BLG-0173Lb, with planet-host mass ratio either $q\simeq 2.5\times 10^{-5}$ or $q\simeq 6.5\times 10^{-5}$, the lowest or among the lowest ever detected. The planetary perturbation is strongly detected, $\Delta\chi^2\sim 10,000$, because it arises from a bright (therefore, large) source passing over and enveloping the planetary caustic: a so-called "Hollywood" event. The factor $\sim 2.5$ offset in $q$ arises because of a previously unrecognized discrete degeneracy between Hollywood events in which the caustic is fully enveloped and those in which only one flank is enveloped, which we dub "Cannae" and "von Schlieffen", respectively. This degeneracy is "accidental" in that it arises from gaps in the data. Nevertheless, the fact that it appears in a $\Delta\chi^2=10,000$ planetary anomaly is striking. We present a simple formalism to estimate the sensitivity of other Hollywood events to planets and show that they can lead to detections close to, but perhaps not quite reaching, the Earth/Sun mass ratio of $3\times 10^{-6}$. This formalism also enables an analytic understanding of the factor $\sim 2.5$ offset in $q$ between the Cannae and von Schlieffen solutions. The Bayesian estimates for the host-mass, system distance, and planet-host projected separation are $M=0.39^{+0.40}_{-0.24}\,M_\odot$, $D_L=4.8^{+1.5}_{-1.8}\,\kpc$, and $a_\perp=3.8\pm 1.6\,\au$. The two estimates of the planet mass are $m_p=3.3^{+3.8}_{-2.1}\,M_\oplus$ and $m_p=8^{+11}_{-6}\,M_\oplus$. The measured lens-source relative proper motion $\mu=6\,\masyr$ will permit imaging of the lens in about 15 years or at first light on adaptive-optics imagers on next-generation telescopes. These will allow to measure the host mass but probably cannot resolve the planet-host mass-ratio degeneracy.

We analyze an ensemble of microlensing events from the 2015 Spitzer microlensing campaign, all of which were densely monitored by ground-based high-cadence survey teams. The simultaneous observations from Spitzer and the ground yield measurements of the microlensing parallax vector $\pi_{\rm E}$, from which compact constraints on the microlens properties are derived, including $\lesssim$25\% uncertainties on the lens mass and distance. With the current sample, we demonstrate that the majority of microlenses are indeed in the mass range of M dwarfs. The planet sensitivities of all 41 events in the sample are calculated, from which we provide constraints on the planet distribution function. In particular, assuming a planet distribution function that is uniform in $\log{q}$, where $q$ is the planet-to-star mass ratio, we find a $95\%$ upper limit on the fraction of stars that host typical microlensing planets of 49\%, which is consistent with previous studies. Based on this planet-free sample, we develop the methodology to statistically study the Galactic distribution of planets using microlensing parallax measurements. Under the assumption that the planet distributions are the same in the bulge as in the disk, we predict that $\sim$1/3 of all planet detections from the microlensing campaigns with Spitzer should be in the bulge. This prediction will be tested with a much larger sample, and deviations from it can be used to constrain the abundance of planets in the bulge relative to the disk.

We present the analysis of OGLE-2016-BLG-0613, for which the lensing light curve appears to be that of a typical binary-lens event with two caustic spikes but with a discontinuous feature on the trough between the spikes. We find that the discontinuous feature was produced by a planetary companion to the binary lens. We find 4 degenerate triple-lens solution classes, each composed of a pair of solutions according to the well-known wide/close planetary degeneracy. One of these solution classes is excluded due to its relatively poor fit. For the remaining three pairs of solutions, the most-likely primary mass is about $M_1\sim 0.7\,M_\odot$ while the planet is a super-Jupiter. In all cases the system lies in the Galactic disk, about half-way toward the Galactic bulge. However, in one of these three solution classes, the secondary of the binary system is a low-mass brown dwarf, with relative mass ratios (1 : 0.03 : 0.003), while in the two others the masses of the binary components are comparable. These two possibilities can be distinguished in about 2024 when the measured lens-source relative proper motion will permit separate resolution of the lens and source.

We present the result of microlensing event MOA-2016-BLG-290, which received observations from the two-wheel Kepler (K2), Spitzer, as well as ground-based observatories. A joint analysis of data from K2 and the ground leads to two degenerate solutions of the lens mass and distance. This degeneracy is effectively broken once the (partial) Spitzer light curve is included. Altogether, the lens is found to be an extremely low-mass star located in the Galactic bulge. MOA-2016-BLG-290 is the first microlensing event for which we have signals from three well-separated ($\sim1$ AU) locations. It demonstrates the power of two-satellite microlensing experiment in reducing the ambiguity of lens properties, as pointed out independently by S. Refsdal and A. Gould several decades ago.

We present an analysis of microlensing event OGLE-2016-BLG-0693, based on the survey-only microlensing observations by the OGLE and KMTNet groups. In order to analyze the light curve, we consider the effects of parallax, orbital motion, and baseline slope, and also refine the result using a Galactic model prior. From the microlensing analysis, we find that the event is a binary composed of a low-mass brown dwarf 49+-20 M_J companion and a K- or G-dwarf host, which lies at a distance 5.0+-0.6 kpc toward the Galactic bulge. The projected separation between the brown dwarf and its host star is less than 5 AU, and thus it is likely that the brown dwarf companion is located in the brown dwarf desert.

The preheating of intergalactic medium(IGM) by structure collapsing and ultraviolet background(UVB) are investigated in cosmological hydrodynamical simulations. When gravitational collapsing is the sole heating mechanism, we find that (1) $60\%, 45\%$ of the IGM are heated up to $S>8, 17$ kev cm$^2$ respectively at $z=0$, but the fractions drop rapidly to a few percents at $z=2$, (2) the entropy of the circum-halo gas $S_{\rm{cir}}$ is higher than the virial entropy for more than $75 \%$ of the halos with masses $M<10^{11.5}$ $M_{\odot}$ since $z=2$, but the fraction higher than the entropy, $S_{\rm{pr}}$, required in preventive model of galaxies formation is only $15-20 \%$ for halos with $M<10^{10.5} M_{\odot}$ at $z=0$, and decreases as redshift increases, (3)assuming a metallicity of $Z \leq 0.03 Z_{\odot}$, the fraction of halos whose circum-halo gas having a cooling time longer than the Hubble time $t_{cool,cir}>t_{\rm{H}}$ is merely $5-10 \%$ at $z \lesssim 0.5$, and even less at $z \geq 1$ for halos with $M<10^{10.5} M_{\odot}$. (4) gas in the filaments undergoes the strongest preheating. Furthermore, we show that the UVB can not enhance the fraction of IGM with $S>17$ kev cm$^2$, but can increase the fraction of low mass halos($<10^{10.5} M_{\odot}$) that having $S_{\rm{cir}}>S_{\rm{pr}}$ to $\sim 70 \%$ at $z=0$, and that having $t_{\rm{cool, cir}}>t_{\rm{H}}$ to $15-30 \%$ at $z \lesssim 0.5$. Our results indicate that preheating due to gravitational collapsing and UVB are inadequate to fulfil the needs of preventative model, especially for halos with $10^{10.5}<M<10^{11.5} M_{\odot}$. Nevertheless, these two mechanisms might cause large scale galactic conformity.

The reduction of the \emph{K2}'s Campaign 9 (K2C9) microlensing data is challenging mostly because of the very crowded field and the unstable pointing of the spacecraft. In this work, we present the first method that can extract microlensing signals from this K2C9 data product. The raw light curves and the astrometric solutions are first derived, using the techniques from Soares-Furtado et al. and Huang et al. for \emph{K2} dense field photometry. We then minimize and remove the systematic effect by performing simultaneous modeling with the microlensing signal. We also derive precise $(K_p-I)$ vs. $(V-I)$ color-color relations that can predict the microlensing source flux in the \emph{Kepler} bandpass. By implementing the color-color relation in the light curve modeling, we show that the microlensing parameters can be better constrained. In the end, we use two example microlensing events, OGLE-2016-BLG-0980 and OGLE-2016-BLG-0940, to test our method.

We report the discovery of a planet-mass companion to the microlens OGLE-2016-BLG-0263L. Unlike most low-mass companions that were detected through perturbations to the smooth and symmetric light curves produced by the primary, the companion was discovered through the channel of a repeating event, in which the companion itself produced its own single-mass light curve after the event produced by the primary had ended. Thanks to the continuous coverage of the second peak by high-cadence surveys, the possibility of the repeating nature due to source binarity is excluded with a $96\%$ confidence level. The mass of the companion estimated by a Bayesian analysis is $M_{\rm p}=4.1_{-2.5}^{+6.5}\ M_{\rm J}$. The projected primary-companion separation is $a_\perp = 6.5^{+1.3}_{-1.9}$ au. The ratio of the separation to the snow-line distance of $a_\perp/a_{\rm sl}\sim 15.4$ corresponds to the region beyond Neptune, the outermost planet of the solar system. We discuss the importance of high-cadence surveys in expanding the range of microlensing detections of low-mass companions and future space-based microlensing surveys.

During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly.

We present the analysis of the binary gravitational microlensing event MOA-2015-BLG-020. The event has a fairly long timescale (about 63 days) and thus the light curve deviates significantly from the lensing model that is based on the rectilinear lens-source relative motion. This enables us to measure the microlensing parallax through the annual parallax effect. The microlensing parallax parameters constrained by the ground-based data are confirmed by the Spitzer observations through the satellite parallax method. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined. It is found that the binary lens is composed of two dwarf stars with masses $M_1 = 0.606 \pm 0.028M_\odot$ and $M_2 = 0.125 \pm 0.006M_\odot$ in the Galactic disk. Assuming the source star is at the same distance as the bulge red clump stars, we find the lens is at a distance $D_L = 2.44 \pm 0.10 kpc$. In the end, we provide a summary and short discussion of all published microlensing events in which the annual parallax effect is confirmed by other independent observations.

The microlens parallax is a crucial observable for conclusively identifying the nature of lens systems in microlensing events containing or composed of faint (even dark) astronomical objects such as planets, neutron stars, brown dwarfs, and black holes. With the commencement of a new era of microlensing in collaboration with space-based observations, the microlens parallax can be routinely measured. In addition, space-based observations can provide opportunities to verify the microlens parallax measured from ground-only observations and to find a unique solution of the lensing lightcurve analysis. However, since most space-based observations cannot cover the full lightcurves of lensing events, it is also necessary to verify the reliability of the information extracted from fragmentary space-based lightcurves. We conduct a test based on the microlensing event OGLE-2016-BLG-0168 created by a binary lens system consisting of almost equal mass M-dwarf stars to demonstrate that it is possible to verify the microlens parallax and to resolve degeneracies by using the space-based lightcurve even though the observations are fragmentary. Since space-based observatories will frequently produce fragmentary lightcurves due to their short observing windows, the methodology of this test will be useful for next-generation microlensing experiments that combine space-based and ground-based collaboration.

We report the discovery of a binary composed of two brown dwarfs, based on the analysis of the microlensing event OGLE-2016-BLG-1469. Thanks to detection of both finite-source and microlens-parallax effects, we are able to measure both the masses $M_1\sim 0.05\ M_\odot$, $M_2\sim 0.01\ M_\odot$, and distance $D_{\rm L} \sim 4.5$ kpc, as well as the projected separation $a_\perp \sim 0.33$ au. This is the third brown-dwarf binary detected using the microlensing method, demonstrating the usefulness of microlensing in detecting field brown-dwarf binaries with separations less than 1 au.

We report the analysis of the first resolved caustic-crossing binary-source microlensing event OGLE-2016-BLG-1003. The event is densely covered by the round-the-clock observations of three surveys. The light curve is characterized by two nested caustic-crossing features, which is unusual for typical caustic-crossing perturbations. From the modeling of the light curve, we find that the anomaly is produced by a binary source passing over a caustic formed by a binary lens. The result proves the importance of high-cadence and continuous observations, and the capability of second-generation microlensing experiments to identify such complex perturbations that are previously unknown. However, the result also raises the issues of the limitations of current analysis techniques for understanding lens systems beyond two masses and of determining the appropriate multiband observing strategy of survey experiments.

We combine $Spitzer$ and ground-based KMTNet microlensing observations to identify and precisely measure an Earth-mass ($1.43^{+0.45}_{-0.32} M_\oplus$) planet OGLE-2016-BLG-1195Lb at $1.16^{+0.16}_{-0.13}$ AU orbiting a $0.078^{+0.016}_{-0.012} M_\odot$ ultracool dwarf. This is the lowest-mass microlensing planet to date. At $3.91^{+0.42}_{-0.46}$ kpc, it is the third consecutive case among the $Spitzer$ "Galactic distribution" planets toward the Galactic bulge that lies in the Galactic disk as opposed to the bulge itself, hinting at a skewed distribution of planets. Together with previous microlensing discoveries, the seven Earth-size planets orbiting the ultracool dwarf TRAPPIST-1, and the detection of disks around young brown dwarfs, OGLE-2016-BLG-1195Lb suggests that such planets might be common around ultracool dwarfs. It therefore sheds light on the formation of both ultracool dwarfs and planetary systems at the limit of low-mass protoplanetary disks.

We analyze the single microlensing event OGLE-2015-BLG-1482 simultaneously observed from two ground-based surveys and from \textit{Spitzer}. The \textit{Spitzer} data exhibit finite-source effects due to the passage of the lens close to or directly over the surface of the source star as seen from \textit{Spitzer}. Such finite-source effects generally yield measurements of the angular Einstein radius, which when combined with the microlens parallax derived from a comparison between the ground-based and the \textit{Spitzer} light curves, yields the lens mass and lens-source relative parallax. From this analysis, we find that the lens of OGLE-2015-BLG-1482 is a very low-mass star with the mass $0.10 \pm 0.02 \ M_\odot$ or a brown dwarf with the mass $55\pm 9 \ M_{J}$, which are respectively located at $D_{\rm LS} = 0.80 \pm 0.19\ \textrm{kpc}$ and $ D_{\rm LS} = 0.54 \pm 0.08\ \textrm{kpc}$, and thus it is the first isolated low-mass microlens that has been decisively located in the Galactic bulge. The fundamental reason for the degeneracy is that the finite-source effect is seen only in a single data point from \textit{Spitzer} and this single data point gives rise to two solutions for $\rho$. Because the $\rho$ degeneracy can be resolved only by relatively high cadence observations around the peak, while the \textit{Spitzer} cadence is typically $\sim 1\,{\rm day}^{-1}$, we expect that events for which the finite-source effect is seen only in the \textit{Spitzer} data may frequently exhibit this $\rho$ degeneracy. For OGLE-2015-BLG-1482, the relative proper motion of the lens and source for the low-mass star is $\mu_{\rm rel} = 9.0 \pm 1.9\ \textrm{mas yr$^{-1}$}$, while for the brown dwarf it is $5.5 \pm 0.5\ \textrm{mas yr$^{-1}$}$. Hence, the degeneracy can be resolved within $\sim 10\ \rm yrs$ from direct lens imaging by using next-generation instruments with high spatial resolution.

The Large European Array for Pulsars combines Europe's largest radio telescopes to form a tied-array telescope that provides high signal-to-noise observations of millisecond pulsars (MSPs) with the objective to increase the sensitivity of detecting low-frequency gravitational waves. As part of this endeavor we have developed a software correlator and beamformer which enables the formation of a tied-array beam from the raw voltages from each of telescopes. We explain the concepts and techniques involved in the process of adding the raw voltages coherently. We further present the software processing pipeline that is specifically designed to deal with data from widely spaced, inhomogeneous radio telescopes and describe the steps involved in preparing, correlating and creating the tied-array beam. This includes polarization calibration, bandpass correction, frequency dependent phase correction, interference mitigation and pulsar gating. A link is provided where the software can be obtained.

Here we present an in-depth study of the behaviour of the Fast Folding Algorithm, an alternative pulsar searching technique to the Fast Fourier Transform. Weaknesses in the Fast Fourier Transform, including a susceptibility to red noise, leave it insensitive to pulsars with long rotational periods (P > 1 s). This sensitivity gap has the potential to bias our understanding of the period distribution of the pulsar population. The Fast Folding Algorithm, a time-domain based pulsar searching technique, has the potential to overcome some of these biases. Modern distributed-computing frameworks now allow for the application of this algorithm to all-sky blind pulsar surveys for the first time. However, many aspects of the behaviour of this search technique remain poorly understood, including its responsiveness to variations in pulse shape and the presence of red noise. Using a custom CPU-based implementation of the Fast Folding Algorithm, ffancy, we have conducted an in-depth study into the behaviour of the Fast Folding Algorithm in both an ideal, white noise regime as well as a trial on observational data from the HTRU-S Low Latitude pulsar survey, including a comparison to the behaviour of the Fast Fourier Transform. We are able to both confirm and expand upon earlier studies that demonstrate the ability of the Fast Folding Algorithm to outperform the Fast Fourier Transform under ideal white noise conditions, and demonstrate a significant improvement in sensitivity to long-period pulsars in real observational data through the use of the Fast Folding Algorithm.

We investigate the evolution of cosmic web since $z=5$ in grid based cosmological hydrodynamical simulations, focusing on the mass and velocity field of both baryonic and cold dark matter. The tidal tensor of density is used as the main method for web identification, with $\lambda_{th}=0.2-1.2$. The evolution trends in baryonic and dark matter are similar, while moderate differences are observed. Sheets appeared early and their large scale pattern may have been set up by $z=3$. In term of mass, filaments superseded sheets as the primary collapsing structures at $z\sim2-3$. Tenuous filaments assembled with each other to form prominent ones at $z<~2$. In accordance with the construction of the frame of sheets, the cosmic divergence velocity, $v_{div}$, had been well developed above 2-3 Mpc by z=3. Afterwards, curl velocity, $v_{curl}$, grown dramatically along with the rising of filaments, become comparable to $v_{div}$, for $<2-3 Mpc$ at $z=0$. The scaling of $v_{curl}$ can be described by the hierarchical turbulence model. The alignment between vorticity and eigenvectors of shear tensor in baryonic matter field resembles dark matter, and is even moderately stronger between $\vec{\omega}$ and $\vec{e}_1$, and $\vec{e}_3$. Compared with dark matter, mildly less baryonic matter is found residing in filaments and clusters, and its vorticity has been developed more significantly below $2-3 Mpc$. These differences may be underestimated due to the limited resolution and lack of star formation in our simulation. The impact of the change of dominant structures in over-dense regions at $z\sim2-3$ on galaxy formation and evolution is shortly discussed.

We report the discovery of a high mass-ratio planet $q=0.012$, i.e., 13 times higher than the Jupiter/Sun ratio. The host mass is not presently measured but can be determined or strongly constrained from adaptive optics imaging. The planet was discovered in a small archival study of high-magnification events in pure-survey microlensing data, which was unbiased by the presence of anomalies. The fact that it was previously unnoticed may indicate that more such planets lie in archival data and could be discovered by similar systematic study. In order to understand the transition from predominantly survey+followup to predominately survey-only planet detections, we conduct the first analysis of these detections in the observational $(s,q)$ plane. Here $s$ is projected separation in units of the Einstein radius. We find some evidence that survey+followup is relatively more sensitive to planets near the Einstein ring, but that there is no statistical difference in sensitivity by mass ratio.

In the process of analyzing an observed light curve, one often confronts various scenarios that can mimic the planetary signals causing difficulties in the accurate interpretation of the lens system. In this paper, we present the analysis of the microlensing event OGLE-2016-BLG-0733. The light curve of the event shows a long-term asymmetric perturbation that would appear to be due to a planet. From the detailed modeling of the lensing light curve, however, we find that the perturbation originates from the binarity of the source rather than the lens. This result demonstrates that binary sources with roughly equal-luminosity components can mimic long-term perturbations induced by planets with projected separations near the Einstein ring. The result also represents the importance of the consideration of various interpretations in planet-like perturbations and of high-cadence observations for ensuring the unambiguous detection of the planet.

We report that the saturation/CGC model of gluon distribution is unstable under action of the chaotic solution in a nonlinear QCD evolution equation, and it evolves to the distribution with a sharp peak at the critical momentum. We find that this gluon condensation is caused by a new kind of shadowing-antishadowing effects, and it leads to a series of unexpected effects in high energy hadron collisions including astrophysical events. For example, the extremely intense fluctuations in the transverse-momentum and rapidity distributions of the gluon jets present the gluon-jet bursts; a sudden increase of the proton-proton cross sections may fill the GZK suppression; the blocking QCD evolution will restrict the maximum available energy of the hadron-hadron colliders.

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548.

We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199-day binary orbit with a companion having a minimum mass of 0.33 solar masses. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches.

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 13 pulsars show both effects. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized ISM features and, in one case, an upper limit to the size of the dispersing region of 13.2 AU. Finally, five pulsars show very nearly quadratic structure functions, which could be indicative of an underlying Kolmogorov medium. Four pulsars show roughly Kolmogorov structure functions and another four show structure functions less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and show that the presence of other trends in the data is the most likely cause.

We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognised, underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars.

Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.

In this paper, we present the analysis of the binary gravitational microlensing event OGLE-2015-BLG-0196. The event lasted for almost a year and the light curve exhibited significant deviations from the lensing model based on the rectilinear lens-source relative motion, enabling us to measure the microlens parallax. The ground-based microlens parallax is confirmed by the data obtained from space-based microlens observations using the {\it Spitzer} telescope. By additionally measuring the angular Einstein radius from the analysis of the resolved caustic crossing, the physical parameters of the lens are determined up to the two-fold degeneracy: $u_0<0$ and $u_0>0$ solutions caused by the well-known "ecliptic" degeneracy. It is found that the binary lens is composed of two M dwarf stars with similar masses $M_1=0.38\pm 0.04\ M_\odot$ ($0.50\pm 0.05\ M_\odot)$ and $M_2=0.38\pm 0.04\ M_\odot$ ($0.55\pm 0.06\ M_\odot$) and the distance to the lens is $D_{\rm L}=2.77\pm 0.23$ kpc ($3.30\pm 0.29$ kpc). Here the physical parameters out and in the parenthesis are for the $u_0<0$ and $u_0>0$ solutions, respectively.

We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613$-$0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as $m_{\rm p} = 1.18^{+0.10}_{-0.09}\text{M}_{\odot}$ for PSR J1918$-$0642 and as high as $m_{\rm p} = 1.928^{+0.017}_{-0.017}\text{M}_{\odot}$ for PSR J1614$-$2230 (both 68.3\% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918$-$0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be $m_{\rm p} = 1.41^{+0.21}_{-0.18}\text{M}_{\odot}$ (68.3\% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of $m_{\rm p} = 1.65^{+0.02}_{-0.02}\text{M}_{\odot}$ (68.3\% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer.

The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for Gravitational Waves to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of RMS jitter in our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve gravitational wave detection pipelines.

The dependence of gas giant planet occurrence rate on stellar metallicity has been firmly established. We extend this so-called planet-metallicity correlation to broader ranges of metallicities and planet masses/radii. In particular, we assume that the planet-metallicity correlation is a power law below some critical saturation threshold, and that the probability of hosting at least one planet is unity for stars with metallicity above the threshold. We then are able to explain the discrepancy between the tentative detection and null detection in previous studies regarding the planet-metallicity correlation for small planets. In particular, we find that the null detection of this correlation can be attributed to the combination of high planet occurrence rate and low detection efficiency. Therefore, a planet-metallicity correlation for small planets cannot be ruled out. We propose that stars with metallicities lower than the Solar value are better targets for testing the planet-metallicity correlation for small planets.

Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of 7 and 3 pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.

We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1 $M_\odot$ companion and possible eclipses that last ~10% of the orbit. The position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.

We report here the Einstein@Home discovery of PSR J1913+1102, a 27.3-ms pulsar found in data from the ongoing Arecibo PALFA pulsar survey. The pulsar is in a 4.95-hr double neutron star (DNS) system with an eccentricity of 0.089. From radio timing with the Arecibo 305-m telescope, we measure the rate of advance of periastron to be 5.632(18) deg/yr. Assuming general relativity accurately models the orbital motion, this corresponds to a total system mass of 2.875(14) solar masses, similar to the mass of the most massive DNS known to date, B1913+16, but with a much smaller eccentricity. The small eccentricity indicates that the second-formed neutron star (the companion of PSR J1913+1102) was born in a supernova with a very small associated kick and mass loss. In that case this companion is likely, by analogy with other systems, to be a light (1.2 solar mass) neutron star; the system would then be highly asymmetric. A search for radio pulsations from the companion yielded no plausible detections, so we can't yet confirm this mass asymmetry. By the end of 2016, timing observations should permit the detection of two additional post-Keplerian parameters: the Einstein delay, which will enable precise mass measurements and a verification of the possible mass asymmetry of the system, and the orbital decay due to the emission of gravitational waves, which will allow another test of the radiative properties of gravity. The latter effect will cause the system to coalesce in ~0.5 Gyr.

Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.

A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.

PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.

We present a combined analysis of the observations of the gravitational microlensing event OGLE-2015-BLG-0479 taken both from the ground and by the {\it Spitzer Space Telescope}. The light curves seen from the ground and from space exhibit a time offset of $\sim 13$ days between the caustic spikes, indicating that the relative lens-source positions seen from the two places are displaced by parallax effects. From modeling the light curves, we measure the space-based microlens parallax. Combined with the angular Einstein radius measured by analyzing the caustic crossings, we determine the mass and distance of the lens. We find that the lens is a binary composed of two G-type stars with masses $\sim 1.0\ M_\odot$ and $\sim 0.9\ M_\odot$ located at a distance $\sim 3$ kpc. In addition, we are able to constrain the complete orbital parameters of the lens thanks to the precise measurement of the microlens parallax derived from the joint analysis. In contrast to the binary event OGLE-2014-BLG-1050, which was also observed by {\it Spitzer}, we find that the interpretation of OGLE-2015-BLG-0479 does not suffer from the degeneracy between $(\pm,\pm)$ and $(\pm,\mp)$ solutions, confirming that the four-fold parallax degeneracy in single-lens events collapses into the two-fold degeneracy for the general case of binary-lens events. The location of the blend in the color-magnitude diagram is consistent with the lens properties, suggesting that the blend is the lens itself. The blend is bright enough for spectroscopy and thus this possibility can be checked from future follow-up observations.

Simultaneous observations of microlensing events from multiple locations allow for the breaking of degeneracies between the physical properties of the lensing system, specifically by exploring different regions of the lens plane and by directly measuring the "microlens parallax". We report the discovery of a 30-55$M_J$ brown dwarf orbiting a K dwarf in microlensing event OGLE-2015-BLG-1319. The system is located at a distance of $\sim$5 kpc toward the Galactic bulge. The event was observed by several ground-based groups as well as by $Spitzer$ and $Swift$, allowing the measurement of the physical properties. However, the event is still subject to an 8-fold degeneracy, in particular the well-known close-wide degeneracy, and thus the projected separation between the two lens components is either $\sim$0.25 AU or $\sim$45 AU. This is the first microlensing event observed by $Swift$, with the UVOT camera. We study the region of microlensing parameter space to which $Swift$ is sensitive, finding that while for this event $Swift$ could not measure the microlens parallax with respect to ground-based observations, it can be important for other events. Specifically, for detecting nearby brown dwarfs and free-floating planets in high magnification events.

We calculate the microlensing event rate and typical time-scales for the free-floating planet (FFP) population that is predicted by the core accretion theory of planet formation. The event rate is found to be ~$1.8\times 10^{-3}$ of that for the stellar population. While the stellar microlensing event time-scale peaks at around 20 days, the median time-scale for FFP events (~0.1 day) is much shorter. Our values for the event rate and the median time-scale are significantly smaller than those required to explain the \cite{Sum+11} result, by factors of ~13 and ~16, respectively. The inclusion of planets at wide separations does not change the results significantly. This discrepancy may be too significant for standard versions of both the core accretion theory and the gravitational instability model to explain satisfactorily. Therefore, either a modification to the planet formation theory is required, or other explanations to the excess of short-time-scale microlensing events are needed. Our predictions can be tested by ongoing microlensing experiment such as KMTNet, and by future satellite missions such as WFIRST and Euclid.

We report on the mass and distance measurements of two single-lens events from the 2015 \emph{Spitzer} microlensing campaign. With both finite-source effect and microlens parallax measurements, we find that the lens of OGLE-2015-BLG-1268 is very likely a brown dwarf. Assuming that the source star lies behind the same amount of dust as the Bulge red clump, we find the lens is a $45\pm7$ $M_{\rm J}$ brown dwarf at $5.9\pm1.0$ kpc. The lens of of the second event, OGLE-2015-BLG-0763, is a $0.50\pm0.04$ $M_\odot$ star at $6.9\pm1.0$ kpc. We show that the probability to definitively measure the mass of isolated microlenses is dramatically increased once simultaneous ground- and space-based observations are conducted.

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

$K2$'s Campaign 9 ($K2$C9) will conduct a $\sim$3.7 deg$^{2}$ survey toward the Galactic bulge from 7/April through 1/July of 2016 that will leverage the spatial separation between $K2$ and the Earth to facilitate measurement of the microlens parallax $\pi_{\rm E}$ for $\gtrsim$127 microlensing events. These will include several that are planetary in nature as well as many short-timescale microlensing events, which are potentially indicative of free-floating planets (FFPs). These satellite parallax measurements will in turn allow for the direct measurement of the masses of and distances to the lensing systems. In this white paper we provide an overview of the $K2$C9 space- and ground-based microlensing survey. Specifically, we detail the demographic questions that can be addressed by this program, including the frequency of FFPs and the Galactic distribution of exoplanets, the observational parameters of $K2$C9, and the array of resources dedicated to concurrent observations. Finally, we outline the avenues through which the larger community can become involved, and generally encourage participation in $K2$C9, which constitutes an important pathfinding mission and community exercise in anticipation of $WFIRST$.

Fast Radio Bursts are millisecond-duration astronomical radio pulses of unknown physical origin that appear to come from extragalactic distances. Previous follow-up observations have failed to find additional bursts at the same dispersion measures (i.e. integrated column density of free electrons between source and telescope) and sky position as the original detections. The apparent non-repeating nature of the fast radio bursts has led several authors to hypothesise that they originate in cataclysmic astrophysical events. Here we report the detection of ten additional bursts from the direction of FRB121102, using the 305-m Arecibo telescope. These new bursts have dispersion measures and sky positions consistent with the original burst. This unambiguously identifies FRB121102 as repeating and demonstrates that its source survives the energetic events that cause the bursts. Additionally, the bursts from FRB121102 show a wide range of spectral shapes that appear to be predominantly intrinsic to the source and which vary on timescales of minutes or shorter. While there may be multiple physical origins for the population of fast radio bursts, the repeat bursts with high dispersion measure and variable spectra specifically seen from FRB121102 support models that propose an origin in a young, highly magnetised, extragalactic neutron star.

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.

The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.