While observations of many high-precision radio pulsars of order $\lesssim$ 1 $\mu$s across the sky are needed for the detection and characterization of a stochastic background of low-frequency gravitational waves (GWs), sensitivity to single sources of GWs requires even higher timing precision. The Argentine Institute of Radio Astronomy (IAR; Instituto Argentino de Radioastronom\'ia) has begun observations of the brightest-known millisecond pulsar, J0437-4715. Even though the two antennas are smaller than other single-dish telescopes previously used for pulsar timing array (PTA) science, the IAR's capability to monitor this pulsar daily coupled with the pulsar's brightness allows for high-precision pulse arrival-time measurements. While upgrades of the facility are currently underway, we show that modest improvements beyond current plans will provide IAR with unparalleled sensitivity to this pulsar. The most stringent upper limits on single GW sources come from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). Observations of PSR J0437-4715 will provide a significant sensitivity increase in NANOGrav's "blind spot" in the sky where fewer pulsars are currently being observed. With state-of-the-art instrumentation installed, we estimate the array's sensitivity will improve by a factor of $\approx$2-7 over 10 years for 20% of the sky with the inclusion of this pulsar as compared to a static version of the PTA used in NANOGrav's most recent limits. More modest instrumentation result in factors of $\approx$1.4-4. We identify four other candidate pulsars as suitable for inclusion in PTA efforts. International PTA efforts will also benefit from inclusion of these data given the potential achievable sensitivity.

Blue Large-Amplitude Pulsators (BLAPs) are a new class of pulsating variable star. They are located close to the hot subdwarf branch in the Hertzsprung-Russell diagram and have spectral classes of late O or early B. Stellar evolution models indicate that these stars are likely radially pulsating, driven by iron group opacity in their interiors. A number of variable stars with a similar driving mechanism exist near the hot subdwarf branch with multi-periodic oscillations caused by either pressure (p) or gravity (g) modes. No multi-periodic signals were detected in the OGLE discovery light curves since it would be difficult to detect short period signals associated with higher-order p modes with the OGLE cadence. Using the RISE instrument on the Liverpool Telescope, we produced high cadence light curves of two BLAPs, OGLE-BLAP-009 ($m_{\mathrm{v}}=15.65$ mag) and OGLE-BLAP-014 ($m_{\mathrm{v}}=16.79$ mag) using a $720$ nm longpass filter. Frequency analysis of these light curves identify a primary oscillation with a period of $31.935\pm0.0098$ mins and an amplitude from a Fourier series fit of $0.236$ mag for BLAP-009. The analysis of BLAP-014 identifies a period of $33.625\pm0.0214$ mins and an amplitude of $0.225$ mag. Analysis of the residual light curves reveals no additional short period variability down to an amplitude of $15.20\pm0.26$ mmag for BLAP-009 and $58.60\pm3.44$ mmag for BLAP-014 for minimum periods of $20$ s and $60$ s respectively. These results further confirm that the BLAPs are monoperiodic.

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.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

Ultra-cool white dwarfs are among the oldest stellar remnants in the Universe. Their efficient gravitational settling and low effective temperatures are responsible for the smooth spectra they exhibit. For that reason, it is not possible to derive their radial velocities or to find the chemistry of the progenitors. The best that can be done is to infer such properties from associated sources, which are coeval. The simplest form of such a system is a common proper motion pair where one star is an evolved degenerate and the other a main-sequence star. In this work, we present the discovery of the first of such a system, the M dwarf LHS 6328 and the ultra-cool white dwarf PSO J1801+625, from the Pan-STARRS 1 3$\pi$ survey and the $Gaia$ Data Release 2. Follow-up spectra were collected covering a usable wavelength range of 3750-24500 {\AA}. Their spectra show that the white dwarf has an effective temperature of 3550K and surface gravity of $\log g = 7.45 \pm 0.13$ or $\log g = 7.49 \pm 0.13$ for a CO or He core, respectively, when compared against synthetic spectra of ultra-cool white dwarf atmosphere models. The system has slightly subsolar metallicity with -0.25 < [Fe/H] < 0.0, and a spatial velocity of (U, V, W) = ($-114.26 \pm 0.24$, $222.94 \pm 0.60$, $10.25 \pm 0.34$) km s $^{-1}$ the first radial velocity and metallicity measurements of an ultra-cool white dwarf. This makes it the first and only benchmark of its kind to date.

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).

Radio emission propagating over an Earth-pulsar line of sight provides a unique probe of the intervening ionized interstellar medium (ISM). Variations in the integrated electron column density along this line of sight, or dispersion measure (DM), have been observed since shortly after the discovery of pulsars. As early as 2006, frequency-dependent dispersion measures have been observed and attributed to several possible causes. Ray-path averaging over different effective light-cone volumes through the turbulent ISM contributes to this effect as will DM misestimation due to radio propagation across compact lensing structures such as those caused by "extreme scattering events". We present methods to assess the variations in frequency-dependent dispersion measures due to the turbulent ISM versus these compact lensing structures along the line of sight. We analyze recent Low-Frequency Array (LOFAR) observations of PSR J2219+4754 to test the underlying physical mechanism of the observed frequency-dependent DM. Previous analyses have indicated the presence of strong lensing due to compact overdensities halfway between the Earth and pulsar. Instead we find the frequency dependence of the DM timeseries for PSR J2219+4754 is consistent with being due solely to ISM turbulence and there is no evidence for any extreme scattering event or small-scale lensing structure. The data show possible deviations from a uniform turbulent medium, suggesting that there may be an enhanced scattering screen near one of the two ends of the line of sight. We present this analysis as an example of the power of low-frequency observations to distinguish the underlying mechanisms in frequency-dependent propagation effects.

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.

Precision pulsar timing can be used for a variety of astrophysical tests from the detection of gravitational waves to probing the properties of the interstellar medium (ISM). Here we present analyses of the noise contributions to pulsar timing residuals from continuous multi-hour observations of seven millisecond pulsars (MSPs). We present scintillation bandwidth measurements for all MSPs in the sample, some for the first time, and scintillation timescale measurements and lower limits for all MSPs for the first time. In addition, we present upper limits on the contribution of pulse phase jitter to the timing residual error for all MSPs. These long observations also allow us to constrain variations in dispersion measures (DMs) on hour-long timescales for several millisecond pulsars. We find that there are no apparent DM variations in any of the MSPs studied on these timescales as expected. In light of new radio telescopes such as the Canadian Hydrogen Intensity Mapping Experiment (CHIME), which will be able to time many pulsars for a short time each day, we search for differences in timing precisions from continuous TOAs and from equivalent length time-discontinuous TOAs. We find no differences in the precision for any MSP in our sample, as expected. We conclude that the TOA variations are consistent with the expected breakdown into template-fitting, jitter, and scintillation errors.

We present a detailed analysis of the International Pulsar Timing Array (IPTA) Second Mock Data Challenge. We tested our analysis methods using the open datasets, and then analyzed the closed datasets. In both the open and the closed datasets, we were able to detect some, but not all, of the injected gravitational wave signals. This work presents two search cases that are not well explored in the pulsar timing array (PTA) literature: a simultaneous search for a stochastic GW background and an individual loud super-massive black hole binary (SMBHB) and a simultaneous search for two SMBHB sources. While we have constructed a cohesive framework for performing these GW searches, our analyses required fine-tuning of the sampling method used in order to appropriately converge. Given the nature of real PTA data in which multiple sources will be present in data, improved techniques will be required in the future to accurately detect and characterize these GW signals.

Wavelength calibration is a routine and critical part of any spectral work-flow, but many astronomers still resort to matching detected peaks and emission lines by hand. We present RASCAL (RANSAC Assisted Spectral CALibration), a python library for automated wavelength calibration of astronomical spectrographs. RASCAL implements recent state-of-the-art methods for wavelength calibration and requires minimal input from a user. In this paper we discuss the implementation of the library and apply it to real-world calibration spectra.

We aim to provide a suite of publicly available spectral data reduction software to facilitate rapid scientific outcomes from time-domain observations. For time resolved observations, an automated pipeline frees astronomers from performance of the routine data analysis tasks to concentrate on interpretation, planning future observations and communication with international collaborators. The project consists of two parts: data processing (ASPIRED) and a graphical user interface (gASPIRED).ASPIREDis written in python 3, and was intentionally developed as a self-consistent reduction pipeline with its own diagnostics and error handling. The pipeline can reduce2D spectral data from raw image to wavelength and flux calibrated 1D spectrum automatically without any user input.gASPIREDis a cross-platform software developedwithElectronon a single code base. It brings interactivity to the software with a well-maintained and user-friendly environment.

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 have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

This paper presents the result of $H\alpha$ imaging for luminous infrared galaxies (LIRGs) and ultra-luminous infrared galaxies (ULIRGs). \textbf{It is } a \textbf{complete subsample of Great Observatories All-sky LIRG Survery (GOALS) with $Dec. \ge -30^{\circ}$}, \textbf{and} consists 148 galaxies with $log(L_{IR}/L_{\odot}) \ge 11.0$. All the $H\alpha$ images were carried out using the 2.16-m telescope \textbf{at the Xinglong Station of the} National Astronomy Observatories, Chinese Academy of Sciences\textbf{ (NAOC),} during the year from 2006 to 2009. We obtained pure $H\alpha$ luminosity for each galaxy and corrected the luminosity for $[NII]$ emission, filter transmission and extinction. We also classified these galaxies based on their morphology and interaction. We found that the distribution of star-forming \textbf{regions} in these galaxies is related to this classification. As the merging process advanced, these galaxies tend to have a more compact distribution of star-forming region, higher \textbf{$L_{IR}$} and warmer IR-color ($f_{60}/f_{100}$). \textbf{These} results imply that the degree of dynamical disturbance plays an important role in determining the distribution of star-forming region.

Some short gamma-ray bursts (SGRBs) are thought to be caused by the mergers of binary neutron stars which may sometimes produce massive neutron star remnants capable of producing extragalactic fast radio bursts (FRBs). We conducted a deep search for FRBs from the sites of six low-redshift SGRBs. We collected high time- and frequency-resolution data from each of the sites for 10 hours using the 2 GHz receiver of the Green Bank Telescope. Two of the SGRB sites we targeted were visible with the Arecibo Radio Telescope with which we conducted an additional 10 hours of 1.4 GHz observations for each. We searched our data for FRBs using the GPU-optimized dedispersion algorithm $\texttt{heimdall}$ and the machine-learning-based package $\texttt{FETCH}$ (Fast Extragalactic Transient Candidate Hunter). We did not discover any FRBs, but would have detected any with peak flux densities in excess of 87 mJy at the Green Bank Telescope or 21 mJy at Arecibo with a signal-to-noise ratio of at least 10. The isotropic-equivalent energy of any FRBs emitted from these sites in our bands during our observations must not have exceeded a few times $10^{38}$ erg, comparable to some of the lowest energy bursts yet seen from the first known repeating FRB 121102.

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.

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 calculate the dispersion measure (DM) contributed by the intergalactic medium (IGM) to the total measured DM for fast radio bursts (FRBs). We use the MareNostrum Instituto de Ciencias del Espacio (MICE) Onion Universe simulation (Fosalba et al. 2008) to track the evolution of the dark matter particle density over a large range of redshifts. We convert this dark matter particle number density to the corresponding free electron density and then integrate it to find the DM as a function of redshift. This approach yields an intergalactic DM of $\rm DM_{\rm IGM}(z = 1) = 800^{+7000}_{-170}$ pc cm$^{-3}$, with the large errors representative of the structure in the IGM. We place limits on the redshifts of the current population of observed FRBs. We also use our results to estimate the host galaxy contribution to the DM for the first repeater, FRB 121102, and show that the most probable host galaxy DM contribution, $\rm DM_{\rm host} \approx 310$ pc cm$^{-3}$, is consistent with the estimate made using the Balmer emission lines in the spectrum of the host galaxy, $\rm DM_{\rm Balmer} = 324$ pc cm$^{-3}$ (Tendulkar et al. 2017). We also compare our predictions for the host galaxy contribution to the DM for the observations of FRB 180924 (Bannister et al. 2019) and FRB 190523 (Ravi et al. 2019), both of which have been localized to a host galaxy.

We describe the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) and its efforts to directly detect and study gravitational waves and other synergistic physics and astrophysics using radio timing observations of millisecond pulsars.

The new field of gravitational wave astrophysics requires a growing pool of students and researchers with unique, interdisciplinary skill sets. It also offers an opportunity to build a diverse, inclusive astronomy community from the ground up. We describe the efforts used by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) NSF Physics Frontiers Center to foster such growth by involving students at all levels in low-frequency gravitational wave astrophysics with pulsar timing arrays (PTAs) and establishing collaboration policies that ensure broad participation by diverse groups. We describe and illustrate the impact of these techniques on our collaboration as a case study for other distributed collaborations.

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.

The current status of and motivation for the 4MOST Survey Strategy, as developed by the Consortium science team, are presented here. Key elements of the strategy are described, such as sky coverage, number of visits and total exposure times in different parts of the sky, and how to deal with different observing conditions. The task of organising the strategy is not simple, with many different surveys that have vastly different target brightnesses and densities, sample completeness levels, and signal-to-noise requirements. We introduce here a number of concepts that we will use to ensure all surveys are optimised. Astronomers who are planning to submit a Participating Survey proposal are strongly encouraged to read this article and any relevant 4MOST Survey articles in this issue of The Messenger such that they can optimally complement and benefit from the planned surveys of the 4MOST Consortium.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We present optical luminosity functions (LFs) of galaxies in the $^{0.1}g$, $^{0.1}r$, and $^{0.1}i$ bands, calculated using data in $\sim 40$ $deg^{2}$ sky area of LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in Southern Galactic Cap. Redshifts for galaxies brighter $r = 18.1$ were obtained mainly with LAMOST. In each band, LFs derived using both parametric and non-parametric maximum likelihood methods agree well with each other. In the $^{0.1}r$ band, our fitting parameters of the Schechter function are $\phi_{*}=(1.65\pm0.36)\times10^{-2}h^{3}Mpc^{-3}$, $M_{*}=-20.69\pm0.06$ mag, and $\alpha=-1.12\pm0.08$, in agreements with previous studies. Separate LFs are also derived for emission line galaxies and absorption line galaxies, respectively. The LFs of absorption line galaxies show a dip at $^{0.1}r \sim 18.5$ and can be well fitted by a double-Gaussian function, suggesting a bi-modality in passive galaxies.

After a decade of great progress in understanding gas flow into, out of, and through the Milky Way, we are poised to merge observations with simulations to build a comprehensive picture of the multi-scale magnetized interstellar medium (ISM). These insights will also be crucial to four bold initiatives in the 2020s: detecting nanohertz gravitational waves with pulsar timing arrays (PTAs), decoding fast radio bursts (FRBs), cosmic B-mode detection, and imaging the Milky Way's black hole with the Event Horizon Telescope (EHT).

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

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 investigate the resolved relation between local extinction and star formation surface density within nearby star-forming galaxies selected from the MaNGA survey. Balmer decrement measurements imply an extinction of the H{\alpha} line emission which scales approximately linearly with the logarithm of the star formation surface density: $ A_{H{\alpha}} = 0.46 \log(\Sigma_{SFR}) + 1.53$. Secondary dependencies are observed such that, at a given $\Sigma_{SFR}$, regions of lower metallicity and/or enhanced H{\alpha} equivalent width (EW) suffer less obscuration than regions of higher metallicity and/or lower H{\alpha} EW. Spaxels lying above the mean relation also tend to belong to galaxies that are more massive, larger and viewed under higher inclination than average. We present a simple model in which the observed trends can be accounted for by a metallicity-dependent scaling between $\Sigma_{SFR}$ and $\Sigma_{dust}$ via a super-linear Kennicutt-Schmidt relation ($n_{KS} \sim 1.47$) and a dust-to-gas ratio which scales linearly with metallicity (DGR($Z_{\odot}$) = 0.01). The relation between the resulting total dust column and observed effective extinction towards nebular regions requires a geometry for the relative distribution of H{\alpha} emitting regions and dust that deviates from a uniform foreground screen and also from an entirely homogeneous mixture of dust and emitting sources. The best-fit model features an H{\alpha} EW and galactocentric distance dependent fraction of the dust mass in a clumpy foreground screen in front of a homogeneous mixture.

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.

The inversion of the white dwarf luminosity function provides an independent way to prove the past star formation history of the Milky Way independent of any cosmological models. In Rowell & Hambly (2011), the effective volume method uses the average properties of all the objects in a given bin, so a significant amount of information is lost in the early stage of the analysis, in this work, I explore the possibility of assigning objects individually in a probabilistic way using the generalised Schmidt density estimator (1/Vmax).

A large sample of white dwarfs is selected by both proper motion and colours from the Pan-STARRS 1 3{\pi} Steradian Survey Processing Version 2 to construct the White Dwarf Luminosity Functions of the discs and halo in the solar neighbourhood. Four-parameter astrometric solutions were recomputed from the epoch data. The generalised maximum volume method is then used to calculate the density of the populations. After removal of crowded areas near the Galactic plane and centre, the final sky area used by this work is 7.833 sr, which is 83% of the 3{\pi} sky and 62% of the whole sky. By dividing the sky using Voronoi tessellation, photometric and astrometric uncertainties are recomputed at each step of the integration to improve the accuracy of the maximum volume. Interstellar reddening is considered throughout the work. We find a disc-to-halo white dwarf ratio of about 100.

The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope which will be online in the early-2020s; and with the development of a new data archive facility for a range of telescopes at the National Astronomical Research Institute of Thailand. The Newton Fund enabled us to collaborate in designing a new versatile generic system that serves all purposes. In the end, we conclude that a single system would not meet the needs of all parties and only adopt similar front-ends while the back-ends are bespoke to our respective systems and data-flows.

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.

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.

Observations of low-frequency gravitational waves will require the highest possible timing precision from an array of the most spin-stable pulsars. We can improve the sensitivity of a pulsar timing array (PTA) to different gravitational-wave sources by observing pulsars with low timing noise over years to decades and distributed across the sky. We discuss observing strategies for a PTA focused on a stochastic gravitational-wave background such as from unresolved supermassive black hole binaries as well as focused on single continuous-wave sources. First we describe the method to calculate a PTA's sensitivity to different gravitational-wave-source classes. We then apply our method to the 45 pulsars presented in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set. For expected amplitudes of the stochastic background, we find that all pulsars contribute significantly over the timescale of decades; the exception is for very pessimistic values of the stochastic background amplitude. For individual single sources, we find that a number of pulsars contribute to the sensitivity of a given source but that which pulsars contribute are different depending on the source, or versus an all-sky metric. Our results seem robust to the presence of temporally-correlated red noise in pulsar arrival times. It is critical to obtain more robust pulsar-noise parameters as they heavily affect our results. Our results show that it is also imperative to locate and time as many high-precision pulsars as possible, as quickly as possible, to maximize the sensitivity of next-generation PTA detectors.

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 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 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.

Precision pulsar timing requires optimization against measurement errors and astrophysical variance from the neutron stars themselves and the interstellar medium. We investigate optimization of arrival time precision as a function of radio frequency and bandwidth. We find that increases in bandwidth that reduce the contribution from receiver noise are countered by the strong chromatic dependence of interstellar effects and intrinsic pulse-profile evolution. The resulting optimal frequency range is therefore telescope and pulsar dependent. We demonstrate the results for five pulsars included in current pulsar timing arrays and determine that they are not optimally observed at current center frequencies. For those objects, we find that better choices of total bandwidth as well as center frequency can improve the arrival-time precision. Wideband receivers centered at somewhat higher frequencies with respect to the currently adopted receivers can reduce required overall integration times and provide significant improvements in arrival time uncertainty by a factor of ~sqrt(2) in most cases, assuming a fixed integration time. We also discuss how timing programs can be extended to pulsars with larger dispersion measures through the use of higher-frequency observations.

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 characteristics of IR properties and MIR variability of RSGs in the LMC are analyzed based on 12 bands of NIR to MIR co-added data from 2MASS, Spitzer and WISE, and $\sim$6.6 years of MIR time-series data collected by the ALLWISE and NEOWISE-R projects. 773 RSGs candidates were compiled from the literature and verified by using the CMD, SED and MIR variability. About 15\% of valid targets in the $IRAC1-IRAC2$/$IRAC2-IRAC3$ diagram may show PAH emission. We show that arbitrary dereddening Q parameters related to the IRAC4, S9W, WISE3, WISE4, and MIPS24 bands could be constructed based on a precise measurement of MIR interstellar extinction law. Several peculiar outliers in our sample are discussed, in which one outlier might be a RSG right before the explosion or an x-AGB star in the very late evolutionary stage based on the MIR spectrum and photometry. There are 744 identified RSGs in the final sample having both the WISE1- and WISE2-band time-series data. The results show that the MIR variability is increasing along with the increasing of brightness. There is a relatively tight correlation between the MIR variability, MLR, and the warm dust or continuum, where the MIR variability is evident for the targets with $K_S-WISE3>1.0~mag$ and $WISE4<6.5~mag$, while the rest of the targets show much smaller MIR variability. The MIR variability is also correlated with the MLR for which targets with larger variability also show larger MLR with an approximate upper limit of $-6.1~M_\odot/yr^{-1}$. Both the variability and the luminosity may be important for the MLR since the WISE4-band flux is increasing exponentially along with the degeneracy of luminosity and variability. The identified RSG sample has been compared with the theoretical evolutionary models and shown that the discrepancy between observation and evolutionary models can be mitigated by considering both variability and extinction.

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.

By precisely monitoring the "ticks" of Nature's most precise clocks (millisecond pulsars), scientists are trying to detect the "ripples in spacetime" (gravitational waves) produced by the inspirals of supermassive black holes in the centers of distant merging galaxies. Here we describe a relatively simple demonstration that uses two metronomes and a microphone to illustrate several techniques used by pulsar astronomers to search for gravitational waves. An adapted version of this demonstration could be used as an instructional laboratory investigation at the undergraduate level.

We present the observed H$\alpha$ flux and derived star formation rates (SFRs) for a fall sample of low$-$surface$-$brightness galaxies (LSBGs). The sample is selected from the fall sky region of the 40$\%$ ALFALFA {\ion{H}{1}} survey $-$ SDSS DR7 photometric data, and all the $H\alpha$ images were obtained using the 2.16 m telescope, operated by the National Astronomy Observatories, Chinese Academy of Sciences. A total of 111 LSBGs were observed and $H\alpha$ flux was measured in 92 of them. Though almost all the LSBGs in our sample are {\ion{H}{1}}$-$rich, their SFRs derived from the extinction and filter$-$transmission$-$corrected $H\alpha$ flux, are less than 1$M_{\sun}$$yr^{-1}$. LSBGs and star forming galaxies have similar {\ion{H}{1}} surface densities, but LSBGs have much lower SFRs and SFR surface densities than star$-$forming galaxies. Our results show that LSBGs deviate from the Kennicutt-Schmidt law significantly, which indicate that they have low star formation efficiency. The SFRs of LSBGs are close to average SFRs in Hubble time and support the previous arguments that most of the LSBGs are stable systems and they tend to seldom contain strong interactions or major mergers during their star formation histories.

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.

Astronomical catalogs derived from wide-field imaging surveys are an important tool for understanding the Universe. We construct an astronomical catalog from 55 TB of imaging data using Celeste, a Bayesian variational inference code written entirely in the high-productivity programming language Julia. Using over 1.3 million threads on 650,000 Intel Xeon Phi cores of the Cori Phase II supercomputer, Celeste achieves a peak rate of 1.54 DP PFLOP/s. Celeste is able to jointly optimize parameters for 188M stars and galaxies, loading and processing 178 TB across 8192 nodes in 14.6 minutes. To achieve this, Celeste exploits parallelism at multiple levels (cluster, node, and thread) and accelerates I/O through Cori's Burst Buffer. Julia's native performance enables Celeste to employ high-level constructs without resorting to hand-written or generated low-level code (C/C++/Fortran), and yet achieve petascale performance.

Gravitational wave (GW) detection with pulsar timing arrays (PTAs) requires accurate noise characterization. The noise of our Galactic-scale GW detector has been systematically evaluated by the Noise Budget and Interstellar Medium Mitigation working groups within the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Intrinsically, individual radio millisecond pulsars (MSPs) used by NANOGrav can have some degree of achromatic red spin noise, as well as white noise due to pulse phase jitter. Along any given line-of-sight, the ionized interstellar medium contributes chromatic noise through dispersion measure (DM) variations, interstellar scintillation, and scattering. These effects contain both red and white components. In the future, with wideband receivers, the effects of frequency-dependent DM will become important. Having anticipated and measured these diverse sources of detector noise, the NANOGrav PTA remains well-poised to detect low-frequency GWs.

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

The survey volume of a proper motion-limited sample is typically much smaller than a magnitude-limited sample. This is because of the noisy astrometric measurements from detectors that are not dedicated for astrometric missions. In order to apply an empirical completeness correction, existing works limit the survey depth to the shallower parts of the sky that hamper the maximum potential of a survey. The number of epoch of measurement is a discrete quantity that cannot be interpolated across the projected plane of observation, so that the survey properties change in discrete steps across the sky. This work proposes a method to dissect the survey into small parts with Voronoi tessellation using candidate objects as generating points, such that each part defines a `mini-survey' that has its own properties. Coupling with a maximum volume density estimator, the new method is demonstrated to be unbiased and recovered {\sim}20% more objects than the existing method in a mock catalogue of a white dwarf-only solar neighbourhood with Pan--STARRS 1-like characteristics. Towards the end of this work, we demonstrate one way to increase the tessellation resolution with artificial generating points, which would be useful for analysis of rare objects with small number counts.

We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases, a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models, the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein time-scale of 7.5/10.8 yr, respectively. In one case, the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects, there is spectroscopic evidence for an intervening absorber but the extra structure seen in the light curves requires a more complex lensing scenario to adequately explain.

White dwarfs are among the most common objects in the stellar halo; however, due to their low luminosity and low number density compared to the stars in the discs of the Milky Way, they are scarce in the observable volume. Hence, they are still poorly understood one hundred years after their discovery as relatively few have been observed. They are crucial to the understanding of several fundamental properties of the Galaxy {\mdash} the geometry, kinematics and star formation history, as well as to the study of the end-stage of stellar evolution for low- and intermediate-mass stars. White dwarfs were traditionally identified by their ultraviolet (UV) excess, however, if they have cooled for a long time, they become so faint in that part of the spectrum that they cannot be seen by the most sensitive modern detectors. Proper motion was then used as a means to identify white dwarf candidates, due to their relatively large space motions compared to other objects with the same colour. The use of proper motion as a selection criterion has proven effective and has yielded large samples of candidates with the SuperCOSMOS Sky Survey and Sloan Digital Sky Survey. In this work I will further increase the sample size with the Panchromatic Synoptic Telescope And Rapid Response System 1 (Pan{\ndash}STARRS1). (cont.)

We discuss the Gaia Data Release 1 (September 2016) and preliminary work on maximising the benefit for cool white dwarf (WD) science in advance of the full parallax catalogue which will appear around one year later in DR2. The Tycho catalogue is used in conjunction with the all-sky ground based astrometric/photometric SuperCOSMOS Sky Survey in order to identify candidate faint common proper motion objects to the Tycho stars. Gaia DR1 is supplemented by the Tycho-Gaia Astrometric Solution catalogue containing some 2 million parallaxes with Hipparcos-like precision for Tycho stars. While hotter, brighter WDs are present in Tycho, cooler examples are much rarer (if present at all) and CPM offers one method to infer precision distances for a statistically useful sample of these very faint WDs.

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.

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 pulsar timing, timing residuals are the differences between the observed times of arrival and the predictions from the timing model. A comprehensive timing model will produce featureless residuals, which are presumably composed of dominating noise and weak physical effects excluded from the timing model (e.g. gravitational waves). In order to apply the optimal statistical methods for detecting the weak gravitational wave signals, we need to know the statistical properties of the noise components in the residuals. In this paper we utilize a variety of non-parametric statistical tests to analyze the whiteness and Gaussianity of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 5-year timing data which are obtained from the Arecibo Observatory and the Green Bank Telescope from 2005 to 2010 (Demorest et al. 2013). We find that most of the data are consistent with white noise; Many data deviate from Gaussianity at different levels, nevertheless, removing outliers in some pulsars will mitigate the deviations.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

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.

We analyze deterministic and random temporal variations in dispersion measure (DM) from the full three-dimensional velocities of pulsars with respect to the solar system, combined with electron-density variations on a wide range of length scales. Previous treatments have largely ignored the pulsar's changing distance while favoring interpretations involving the change in sky position from transverse motion. Linear trends in pulsar DMs seen over 5-10~year timescales may signify sizable DM gradients in the interstellar medium (ISM) sampled by the changing direction of the line of sight to the pulsar. We show that motions parallel to the line of sight can also account for linear trends, for the apparent excess of DM variance over that extrapolated from scintillation measurements, and for the apparent non-Kolmogorov scalings of DM structure functions inferred in some cases. Pulsar motions through atomic gas may produce bow-shock ionized gas that also contributes to DM variations. We discuss possible causes of periodic or quasi-periodic changes in DM, including seasonal changes in the ionosphere, annual variation of the solar elongation angle, structure in the heliosphere-ISM boundary, and substructure in the ISM. We assess the solar cycle's role on the amplitude of ionospheric and solar-wind variations. Interstellar refraction can produce cyclic timing variations from the error in transforming arrival times to the solar system barycenter. We apply our methods to DM time series and DM gradient measurements in the literature and assess consistency with a Kolmogorov medium. Finally, we discuss the implications of DM modeling in precision pulsar timing experiments.

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.

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.

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.

Since the release of LAMOST (the Large Sky Area Multi-Object Fiber Spectroscopic Telescope) catalog, we have the opportunity to use the LAMOST DR2 stellar catalog and \emph{WISE All-sky catalog} to search for 22 $\mu$m excess candidates. In this paper, we present 10 FGK candidates which show an excess in the infrared (IR) at 22 $\mu$m. The ten sources are all the newly identified 22 $\mu$m excess candidates. Of these 10 stars, 5 stars are F type and 5 stars are G type. The criterion for selecting candidates is $K_s-[22]_{\mu m}\geq0.387$. In addition, we present the spectral energy distributions (SEDs) covering wavelength from optical to mid-infrared band. Most of them show an obvious excess from 12 $\mu$m band and three candidates even show excess from 3.4 $\mu$m. To characterize the amount of dust, we also estimate the fractional luminosity of ten 22 $\mu$m excess candidates.

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.

We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains only a negligible part of their noise budget.

We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.

Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.

We report the discovery of an isolated compact galaxy triplet SDSS J084843.45+164417.3, which is first detected by the LAMOST spectral survey and then confirmed by the spectroscopic observation of the BFOSC of the 2.16 meter telescope. It is found that this triplet is an isolated and extremely compact system, which has an aligned configuration and very small radial velocity dispersion. The member galaxies have similar colors and show marginal star formation activities. These results enhance the opinion that the compact triplets are well-evolved systems rather than the hierarchically forming structures. This occasional discovery reveals the limitations of the fiber spectral redshift surveys in studying such compact system, and declares the necessity of additional observations to complete the current redshift sample.

We report 21-yr timing of one of the most precise pulsars: PSR J1713+0747. Its pulse times of arrival are well modeled by a comprehensive pulsar binary model including its three-dimensional orbit and a noise model that incorporates correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square $\sim 92$ ns. The new dataset allows us to update and improve previous measurements of the system properties, including the masses of the neutron star ($1.31\pm0.11$ $M_{\odot}$) and the companion white dwarf ($0.286\pm0.012$ $M_{\odot}$) and the parallax distance $1.15\pm0.03$ kpc. We measured the intrinsic change in orbital period, $\dot{P}^{\rm Int}_{\rm b}$, is $-0.20\pm0.17$ ps s$^{-1}$, which is not distinguishable from zero. This result, combined with the measured $\dot{P}^{\rm Int}_{\rm b}$ of other pulsars, can place a generic limit on potential changes in the gravitational constant $G$. We found that $\dot{G}/G$ is consistent with zero [$(-0.6\pm1.1)\times10^{-12}$ yr$^{-1}$, 95\% confidence] and changes at least a factor of $31$ (99.7\% confidence) more slowly than the average expansion rate of the Universe. This is the best $\dot{G}/G$ limit from pulsar binary systems. The $\dot{P}^{\rm Int}_{\rm b}$ of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation $\kappa_D=(-0.9\pm3.3)\times10^{-4}$ (95\% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters $\Delta$, which describes the violation of strong equivalence principle, and $\hat{\alpha}_3$, which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95\% confidence, $\Delta<0.01$ and $\hat{\alpha}_3<2\times10^{-20}$ based on PSR J1713+0747.

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.

The traditional Schmidt density estimator has been proven to be unbiased and effective in a magnitude-limited sample. Previously, efforts have been made to generalize it for populations with non-uniform density and proper motion-limited cases. This work shows that the then-good assumptions for a proper motion-limited sample are no longer sufficient to cope with modern data. Populations with larger differences in the kinematics as compared to the local standard of rest are most severely affected. We show that this systematic bias can be removed by treating the discovery fraction inseparable from the generalized maximum volume integrand. The treatment can be applied to any proper motion-limited sample with good knowledge of the kinematics. This work demonstrates the method through application to a mock catalogue of a white dwarf-only solar neighbourhood for various scenarios and compared against the traditional treatment using a survey with Pan-STARRS-like characteristics.

The population of Low Surface Brightness (LSB) galaxies is crucial for understanding the extremes of galaxy formation and evolution of the universe. As LSB galaxies are mostly rich in gas (HI), the alpha.40-SDSS DR7 sample is absolutely one of the best survey combinations to select a sample of them in the local Universe. Since the sky backgrounds are systematically overestimated for galaxy images by the SDSS photometric pipeline, particularly for luminous galaxies or galaxies with extended low surface brightness outskirts, in this paper, we above all estimated the sky backgrounds of SDSS images in the alpha.40-SDSS DR7 sample, using a precise method of sky subtraction. Once subtracting the sky background, we did surface photometry with the Kron elliptical aperture and fitted geometric parameters with an exponential profile model for each galaxy image. Basing on the photometric and geometric results, we further calculated the B-band central surface brightness, mu_{0}(B), for each galaxy and ultimately defined a sample of LSB galaxies consisting of 1129 galaxies with mu_{0}(B) > 22.5 mag arcsec^{-2} and the minor-to-major axis ratio b/a > 0.3. This HI-selected LSB galaxy sample from the alpha.40-SDSS DR7 is a relatively unbiased sample of gas-rich and disk-dominated LSB galaxies, which is complete both in HI observation and the optical magnitude within the limit of SDSS DR7 photometric survey. We made analysis on optical and radio 21cm HI properties of this LSB galaxy sample. Additionally, we statistically investigated the environment of our LSB galaxies, and found that up to 92% of the total LSB galaxies have less than 8 neighbouring galaxies, which strongly evidenced that LSB galaxies prefer to reside in the low-density environment.

The use of Reduced Proper Motion in identifying isolated white dwarfs has long been used as a proxy for the absolute magnitude in a population with known kinematics. This, however, introduces a proper motion detection limit on top of the existing photometric limit. How the survey volume is hampered by this extra parameter is discussed in Hambly et al. 2012. In this work, we discuss some robust outlier rejection methods in order to minimise the proper motion limit and hence maximise the survey volume. The generalised volume, corrected for the distance of the Sun from the Galactic Plane, is integrated explicitly.

Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.

Free electrons in the interstellar medium cause frequency-dependent delays in pulse arrival times due to both scattering and dispersion. Multi-frequency measurements are used to estimate and remove dispersion delays. In this paper, we focus on the effect of any non-simultaneity of multi-frequency observations on dispersive delay estimation and removal. Interstellar density variations combined with changes in the line-of-sight from pulsar and observer motions cause dispersion measure variations with an approximately power-law power spectrum, augmented in some cases by linear trends. We simulate time series, estimate the magnitude and statistical properties of timing errors that result from non-simultaneous observations, and derive prescriptions for data acquisition that are needed in order to achieve a specified timing precision. For nearby, highly stable pulsars, measurements need to be simultaneous to within about one day in order that the timing error from asynchronous DM correction is less than about 10 ns. We discuss how timing precision improves when increasing the number of dual-frequency observations used in dispersion measure estimation for a given epoch. For a Kolmogorov wavenumber spectrum, we find about a factor of two improvement in precision timing when increasing from two to three observations but diminishing returns thereafter.

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

We present new narrow-band (H$\alpha$ and [OIII]) imagings and optical spectrophotometry of HII regions for a gas-rich low surface brightness irregular galaxy, KKR 17. The central surface brightness of the galaxy is $\mu_0(B)$ = 24.15 $\pm$0.03 mag~sec$^{-2}$. The galaxy was detected by \emph{Arecibo Legacy Fast ALFA survey} (ALFALFA), and its mass is dominated by neutral hydrogen (HI) gas. In contrast, both the stellar masses of the bright HII and diffuse stellar regions are small. In addition, the fit to the spectral energy distribution to each region shows the stellar populations of HII and diffuse regions are different. The bright HII region contains a large fraction of O-type stars, revealing the recent strong star formation, whereas the diffuse region is dominated by median age stars, which has a typical age of $\sim$ 600 Myrs. Using the McGaugh's abundance model, we found that the average metallicity of KKR 17 is 12 + (O/H) = 8.0 $\pm$ 0.1. The star formation rate of KKR 17 is 0.21$\pm$0.04 M$_{\odot}$/yr, which is $\sim$1/5 of our Milky Way's. Based on the analysis results to young stellar clusters in HII region, it is found that the bright HII region showed two sub-components with different velocities and metallicities. This may be caused by the outflow of massive stars or merging events. However, the mechanism of triggering star formation in the HII region is still uncertain.

With the newly available SPIRE images at 250 and 500 micron from Herschel Space Observatory, we study quantitative correlations over a sub-kpc scale among three distinct emission components in the interstellar medium of the nearby spiral galaxy M 81 (NGC 3031): (a) $I_{8}$ or $I_{24}$, the surface brightness of the mid-infrared emission observed in the Spitzer IRAC 8 or MIPS 24 micron band, with $I_8$ and $I_{24}$ being dominated by the emissions from Polycyclic Aromatic Hydrocarbons (PAHs) and very small grains (VSGs) of dust, respectively; (b) $I_{500}$, that of the cold dust continuum emission in the Herschel SPIRE 500 micron band, dominated by the emission from large dust grains heated by evolved stars, and (c) $I_{{\rm H}\alpha}$, a nominal surface brightness of the H$\alpha$ line emission, from gas ionized by newly formed massive stars. The results from our correlation study, free from any assumption on or modeling of dust emissivity law or dust temperatures, present solid evidence for significant heating of PAHs and VSGs by evolved stars. In the case of M 81, about 67% (48%) of the 8 micron (24 micron) emission derives its heating from evolved stars, with the remainder attributed to radiation heating associated with ionizing stars.

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.

We present an analysis of the metallicity and star formation activities of HII regions in the interacting system Arp 86, based on the first scientific observation of the multi-object spectroscopy on the 2.16m Telescope at Xinglong Observatory. We find that the oxygen abundance gradient in Arp 86 is flatter than that in normal disk galaxies, which confirms that gas inflows caused by tidal forces during encounters can flatten the metallicity distributions in galaxies. The companion galaxy NGC 7752 is currently experiencing a galaxy-wide starburst with higher surface density of star formation rate than the main galaxy NGC 7753, which can be explained that the companion galaxy is more susceptible to the effects of interaction than the primary. We also find that the galaxy 2MASX J23470758+2926531 has similar abundance and star formation properties to NGC 7753, and may be a part of the Arp 86 system.

In this paper we present a catalog which includes 141 bright candidates ($\leq10.27$ mag, V band) showing the infrared (IR) excess at 22 $\mu$m. Of which, 38 stars are known IR excess stars or disk, 23 stars are double or multiple stars and 4 are Be stars. While the remaining more than 70 stars are identified as the 22 $\mu$m excess candidates in our work. The criterion of selecting candidates is $K_s-[22]_{\mu m}$. All these candidates are selected from \emph{WISE} All-sky data cross-correlated with \emph{Hipparcos} Main Catalog and the likelihood-ratio technique is employed. Considering the effect of background, we introduce the \emph{IRAS} 100 $\mu$m level to exclude the high background. We also estimated the coincidence probability of these sources. In addition, we presented the optical to mid-infrared SEDs and optical images of all the candidates, and gave the observed optical spectra of 6 stars with NAOC's 2.16-m telescope. To measure for the dust amount around each star, the fractional luminosity is also provided. We also test whether our method of selecting IR excess stars can be used to search for extra-solar planets, we cross-matched our catalog with known IR-excess stars having planets but none is matched. Finally, we give the fraction of stars showing IR-excess for different spectral type of main-sequence stars.

Using the Herschel ATLAS science demonstration phase data crossidentified with SDSS DR7 spectra, we select 297 galaxies with F250{\mu}m > 5{\sigma}. The sample galaxies are classified into five morphological types, and more than 40% of the galaxies are peculiar/compact galaxies. The peculiar galaxies show higher far-infrared/submillimeter luminosity-to-mass ratios than the other types. We perform and analyze the correlations of far-infrared/submillimeter and H{\alpha} luminosities for different morphological types and different spectral types. The Spearman rank coefficient decreases and the scatter increases with the wavelength increasing from 100 {\mu}m to 500 {\mu}m. We conclude that a single Herschel SPIRE band is not good for tracing star formation activities in galaxies. AGNs contribute less to the far-infrared/submillimeter luminosities and do not show a difference from star-forming galaxies. However, the earlier type galaxies present significant deviations from the best fit of star-forming galaxies.