We describe 14 years of public data from the Parkes Pulsar Timing Array (PPTA), an ongoing project that is producing precise measurements of pulse times of arrival from 26 millisecond pulsars using the 64-m Parkes radio telescope with a cadence of approximately three weeks in three observing bands. A comprehensive description of the pulsar observing systems employed at the telescope since 2004 is provided, including the calibration methodology and an analysis of the stability of system components. We attempt to provide full accounting of the reduction from the raw measured Stokes parameters to pulse times of arrival to aid third parties in reproducing our results. This conversion is encapsulated in a processing pipeline designed to track provenance. Our data products include pulse times of arrival for each of the pulsars along with an initial set of pulsar parameters and noise models. The calibrated pulse profiles and timing template profiles are also available. These data represent almost 21,000 hrs of recorded data spanning over 14 years. After accounting for processes that induce time-correlated noise, 22 of the pulsars have weighted root-mean-square timing residuals of < 1 ${\mu}$s in at least one radio band. The data should allow end users to quickly undertake their own gravitational-wave analyses (for example) without having to understand the intricacies of pulsar polarisation calibration or attain a mastery of radio-frequency interference mitigation as is required when analysing raw data files.

The stronger than expected 21-cm absorption was observed by EDGES recently, and another anomaly of $^8$Be transitions would be signatures of new interactions. These two issues may be related to each other, e.g., pseudoscalar $A$ mediated fermionic millicharged dark matter (DM), and the 21-cm absorption could be induced by photon mediated scattering between MeV millicharged DM and hydrogen. This will be explored in this paper. For fermionic millicharged DM $\bar{\chi} \chi$ with masses in a range of $2 m_A < 2 m_{\chi} < 3 m_A$, the p-wave annihilation $\bar{\chi} \chi \to A A$ would be dominant during DM freeze-out. The s-wave annihilation $\bar{\chi} \chi$ $\to A, \gamma $ $\to e^+ e^-$ is tolerant by constraints from CMB and the 21-cm absorption. The millicharged DM can evade constraints from direct detection experiments. The process of $K^+ \to \pi^+ \pi^0$ with the invisible decay $\pi^0 \to \bar{\chi} \chi$ could be employed to search for the millicharged DM, and future high intensity $K^+$ sources, such as NA62, will do the job.

We present first results from pilot observations using a phased array feed (PAF) mounted on the Parkes 64-m radio telescope. The observations presented here cover a frequency range from 1150 to 1480 MHz and are used to show the ability of PAFs to suppress standing wave problems by a factor of $\sim10$ which afflict normal feeds. We also compare our results with previous HIPASS observations and with previous HI images of the Large Magellanic Cloud. Drift scan observations of the GAMA G23 field resulted in direct HI detections at $z=0.0043$ and $z=0.0055$ of HIPASS galaxies J2242-30 and J2309-30. Our new measurements generally agree with archival data in spectral shape and flux density, with small differences being due to differing beam patterns. We also detect signal in the stacked HI data of 1094 individually undetected galaxies in the GAMA G23 field in the redshift range $0.05 \leq z \leq 0.075$. Finally, we use the low standing wave ripple and wide bandwidth of the PAF to set a $3\sigma$ upper limit to any positronium recombination line emission from the Galactic Centre of $<0.09$ K, corresponding to a recombination rate of $<3.0\times10^{45}\,\mathrm{s}^{-1}$.

During February 2016, CSIRO Astronomy and Space Science and the Max-Planck-Institute for Radio Astronomy installed, commissioned and carried out science observations with a phased array feed (PAF) receiver system on the 64m diameter Parkes radio telescope. Here we demonstrate that the PAF can be used for pulsar observations and we highlight some unique capabilities. We demonstrate that the pulse profiles obtained using the PAF can be calibrated and that multiple pulsars can be simultaneously observed. Significantly, we find that an intrinsic polarisation leakage of -31dB can be achieved with a PAF beam offset from the centre of the field of view. We discuss the possibilities for using a PAF for future pulsar observations and for searching for fast radio bursts with the Parkes and Effelsberg telescopes.

In the solar wind, power spectral density (PSD) of the magnetic field fluctuations generally follow the so-called Kolmogorov spectrum f^-5/3 in the inertial range, where the dynamics is thought to be dominated by nonlinear interactions between counter-propagating incompressible Alfv\'en wave parquets. These features are thought to be ubiquitous in space plasmas. The present study gives a new and more complex picture of magnetohydrodynamics (MHD) turbulence as observed in the terrestrial magnetosheath. The study uses three years of in-situ data from the Cluster mission to explore the nature of the magnetic fluctuations at MHD scales in different locations within the magnetosheath, including flanks and subsolar regions. It is found that the magnetic field fluctuations at MHD scales generally have a PSD close to f^-1 (shallower than the Kolmogorov one f^-5/3) down to the ion characteristic scale, which recalls the energy containing scales of solar wind turbulence. The Kolmogorov spectrum is observed only away from the bow shock toward the flank and the magnetopause regions in 17% of the analyzed time intervals. Measuring the magnetic compressibility, it is shown that only a fraction (35%) of the observed Kolmogorov spectra were populated by shear Alfv\'enic fluctuations, whereas the majority of the events (65%) was found to be dominated by compressible magnetosonic-like fluctuations, which contrasts with well-known turbulence properties in the solar wind. This study gives a first comprehensive view of the origin of the f^-1 and the transition to the Kolmogorov inertial range; both questions remain controversial in solar wind turbulence.

We report the observations of an electron vortex magnetic hole corresponding to a new type of coherent structures in the magnetosheath turbulent plasma using the Magnetospheric Multiscale (MMS) mission data. The magnetic hole is characterized by a magnetic depression, a density peak, a total electron temperature increase (with a parallel temperature decrease but a perpendicular temperature increase), and strong currents carried by the electrons. The current has a dip in the center of the magnetic hole and a peak in the outer region of the magnetic hole. The estimated size of the magnetic hole is about 0.23 \r{ho}i (~ 30 \r{ho}e) in the circular cross-section perpendicular to its axis, where \r{ho}i and \r{ho}e are respectively the proton and electron gyroradius. There are no clear enhancement seen in high energy electron fluxes, but an enhancement in the perpendicular electron fluxes at ~ 90{\deg} pitch angles inside the magnetic hole is seen, implying that the electron are trapped within it. The variations of the electron velocity components Vem and Ven suggest that an electron vortex is formed by trapping electrons inside the magnetic hole in the circular cross-section (in the M-N plane). These observations demonstrate the existence of a new type of coherent structures behaving as an electron vortex magnetic hole in turbulent space plasmas as predicted by recent kinetic simulations.

We present the first installation and characterization of a phased array feed (PAF) on the 64 m Parkes radio telescope. The combined system operates best between 0.8 GHz and 1.74 GHz where the beamformed noise temperature is between 45 K and 60 K, the aperture efficiency ranges from 70% to 80%, and the effective field of view is 1.4 deg$^2$ at 1310 MHz. After a 6-month trial observing program at Parkes, the PAF will be installed on the 100 m antenna at Effelsberg. This is the first time a PAF has been installed on a large single-antenna radio telescope and made available to astronomers.

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.

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

Gravitational wave background results from the superposition of gravitational waves generated from all sources across the Universe. Previous efforts on detecting such a background with pulsar timing arrays assume it is an isotropic Gaussian background with a power law spectrum. However, when the number of sources is limited, the background might be non-Gaussian or the spectrum might not be a power law. Correspondingly previous analysis may not work effectively. Here we use a method --- Bayesian Nonparametric Analysis --- to try to detect a generic gravitational wave background, which directly sets constraints on the feasible shapes of the pulsar timing signals induced by a gravitational wave background and allows more flexible forms of the background. Our Bayesian nonparametric analysis will infer if a gravitational wave background is present in the data, and also estimate the parameters that characterize the background. This method will be much more effective than the conventional one assuming the background spectrum follows a power law in general cases. While the context of our discussion focuses on pulsar timing arrays, the analysis itself is directly applicable to detect and characterize any signals that arise from the superposition of a large number of astrophysical events.

Using two volume-limited Main galaxy samples of the Sloan Digital Sky Survey Data Release 10 (SDSS DR10), I investigate the environmental dependence of galaxy age, and get the same conclusions in two volume-limited Main galaxy samples: old galaxies exist preferentially in the densest regions of the universe, while young galaxies are located preferentially in low density regions. Such an age-density relation is likely a combination of a strong age-stellar mass relation and the stellar mass-density relation.

Gravitational wave burst is a catch-all category for signals whose durations are shorter than the observation period. We apply a method new to gravitational wave data analysis --- Bayesian non-parameterics --- to the problem of gravitational wave detection, with an emphasis on pulsar timing array observations. In Bayesian non-parametrics, constraints are set on the function space that may be reasonably thought to characterize the range of gravitational-wave signals. This differs from the approaches currently employed or proposed, which focus on introducing parametric signal models or looking for excess power as evidence of the presence of a gravitational wave signal. Our Bayesian nonparametrics analysis method addresses two issues: (1) investigate if a gravitational wave burst is present in the data; (2) infer the sky location of the source and the duration of the burst. Compared with the popular method proposed by Finn \& Lommen, our method improves in two aspects: (1) we can estimate the burst duration by adding the prior that the gravitational wave signals are smooth, while Finn \& Lommen ignored this important point; (2) we perform a full Bayesian analysis by marginalizing over all possible parameters and provide robust inference on the presence of gravitational waves, while Finn \& Lommen chose to optimize over parameters, which would increase false alarm risk and also underestimate the parameter uncertainties.

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.

Ground-based and space-borne observatories studying exoplanetary transits now and in the future will considerably increase the number of known exoplanets and the precision of the measured times of transit minima. Variations in the transit times can not only be used to infer the presence of additional planets, but might also provide opportunities for testing new physics in the places beyond the Solar system. In this work, we take deviation from the inverse-square law of gravity as an example, focus on the fifth-force-like Yukawa-type correction to the Newtonian gravitational force which parameterizes this deviation, investigate its effects on the secular transit timing variations and analyze their observability in exoplanetary systems. It is found that the most optimistic values of Yukawa-type secular transit timing variations are at the level of $\sim 0.1$ seconds per year. Those values unfortunately appear only in rarely unique cases and, most importantly, they are still at least two orders of magnitude below the current capabilities of observations. Such a deviation from the inverse-square law of gravity is likely too small to detect for the foreseeable future. Meanwhile, systematic uncertainties, such as the presence of additional and unknown planets, will likely be exceptionally difficult to remove from a signal that should be seen.

Using the LOWZ and CMASS samples of the ninth data release (DR9) from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), I investigate properties of star forming galaxies and AGNs. The CMASS sample seriously suffers from the radial selection effect, even within the redshift 0.44<z<0.6, which will lead to statistical conclusions in the CMASS sample being less likely robust. In the LOWZ sample, the fraction of star-forming galaxies is nearly constant from the lowest density regime to the densest regime; the AGN fraction is also insensitive to the local environment. In addition, I note that in the LOWZ sample, stellar mass and stellar velocity dispersion distributions of star forming galaxies and AGNs are nearly same.

In this study, I investigate the environmental dependence of galaxy colors in the CMASS sample of the Sloan Digital Sky Survey Data Release 9 (SDSS DR9). To decrease the radial selection effect, I divide the CMASS sample into subsamples with a redshift binning size of 0.01 and analyze the environmental dependence of the u-r, u-g, g-r, r-i and i-z colors for these subsamples in each redshift bin. Statistical analysis shows that all the five colors weakly correlate with the local environment, which may mean that the environmental processes responsible for a galaxy's properties proceed slowly over cosmic time.

We present a first statistical study of subproton and electron scales turbulence in the terrestrial magnetosheath using the Cluster Search Coil Magnetometer (SCM) waveforms of the STAFF instrument measured in the frequency range [1,180] Hz. It is found that clear spectral breaks exist near the electron scale, which separate two power-law like frequency bands referred to as the dispersive and the electron dissipation ranges. The frequencies of the breaks f_b are shown to be well correlated with the electron gyroscale \rho_e rather than with the electron inertial length de. The distribution of the slopes below fb was found to be narrow and peaks near -2.9, while that of the slopes above fb was found broader, peaks near -5.2 and has values as low as -7.5. This is the first time that such steep power-law spectra are reported in space plasma turbulence. These observations provide strong constraints on theoretical modeling of kinetic turbulence and dissipation in collisionless magnetized plasmas.

We demonstrate how observations of pulsars can be used to help navigate a spacecraft travelling in the solar system. We make use of archival observations of millisecond pulsars from the Parkes radio telescope in order to demonstrate the effectiveness of the method and highlight issues, such as pulsar spin irregularities, which need to be accounted for. We show that observations of four millisecond pulsars every seven days using a realistic X-ray telescope on the spacecraft throughout a journey from Earth to Mars can lead to position determinations better than approx. 20km and velocity measurements with a precision of approx. 0.1m/s.

The Yukawa correction to the Newtonian gravitational force is accepted as a parameterization of deviations from the inverse-square law of gravity which might be caused by new physics beyond the standard model of particles and the general theory of relativity. We investigate these effects on the clock onboard a drag-free satellite: dynamics of the satellite and influence on the time transfer link. We find the Yukawa signal in the time transfer is much harder to detect with current state of clocks than those effects on the dynamics, especially the secular change of periastron, by laser ranging in the case of an artificial Earth satellite carrying a frequency standard with an orbit of $a=10^7$\,m and $e=0.01$.

Using the galaxy data of the Sloan Digital Sky Survey Data Release 8 (SDSS DR8), I explore whether the concentration index is a good morphological classification tool and find that a reasonably pure late-type galaxy sample can be constructed with the choice of r-band concentration index ci=2.85. The opposite is not true, however, due to the fairly high contamination of an early-type sample by late-type galaxies. In such an analysis, the influence of selection effects is less important. To disentangle correlations of the morphology and concentration index with stellar mass, star formation rate (SFR), specific star formation rate (SSFR) and active galactic nucleus (AGN) activity, I investigate correlations of the concentration index with these properties at a fixed morphology and correlations of the morphology with these properties at a fixed concentration index. It is found that at a fixed morphology, high-concentration galaxies are preferentially more massive and have a lower SFR and SSFR than low-concentration galaxies, while at a fixed concentration index, elliptical galaxies are preferentially more massive and have a lower SFR and SSFR than spiral galaxies. This result shows that the stellar mass, SFR and SSFR of a galaxy are correlated with its concentration index as well as its morphology. In addition, I note that AGNs are preferentially found in more concentrated galaxies only in the spiral galaxy sample.

Using two volume-limited Main galaxy samples of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we explore influences of galaxy interactions on AGN activity. It is found that in the faint volume-limited sample, paired galaxies have a slightly higher AGN fraction than isolated galaxies, whereas in the luminous volume-limited sample, an opposite trend can be observed. The significance is <1 sigma. Thus, we do not observe strong evidence that interactions or mergers likely trigger the AGN activity.

For pulsar projects it is often necessary to predict the pulse phase in advance, for example, when preparing for new observations. Interpolation of the pulse phase between existing measurements is also often required, for example, when folding X-ray or gamma-ray observations according to the radio pulse phase. Until now these procedures have been done using various ad hoc methods. The purpose of this paper is to show how to interpolate or predict the pulse phase optimally using statistical models of the various noise processes and the phase measurement uncertainty.

Pulsar timing arrays act to detect gravitational waves by observing the small, correlated effect the waves have on pulse arrival times at Earth. This effect has conventionally been evaluated assuming the gravitational wave phasefronts are planar across the array, an assumption that is valid only for sources at distances $R\gg2\pi{}L^2/\lambda$, where $L$ is physical extent of the array and $\lambda$ the radiation wavelength. In the case of pulsar timing arrays (PTAs) the array size is of order the pulsar-Earth distance (kpc) and $\lambda$ is of order pc. Correspondingly, for point gravitational wave sources closer than $\sim100$~Mpc the PTA response is sensitive to the source parallax across the pulsar-Earth baseline. Here we evaluate the PTA response to gravitational wave point sources including the important wavefront curvature effects. Taking the wavefront curvature into account the relative amplitude and phase of the timing residuals associated with a collection of pulsars allows us to measure the distance to, and sky position of, the source.

From the Main galaxy sample of the Sloan Digital Sky Survey Data Release 6 (SDSS DR6), we construct two volume-limited samples above and below the value of M, to explore the difference of the environmental dependence of galaxy properties between galaxies above and below the value of M . We measure the local three-dimensional galaxy density in a comoving sphere with a radius of the distance to the 5th nearest galaxy for each galaxy, and compare basic properties of galaxies in the lowest density regime with those of galaxies in the densest regime. It is found that the galaxy luminosity strongly depend on local environments only for galaxies above M, but this dependence is very weak for galaxies below M . It is worth noting that g-r color, concentration index ci and galaxy morphologies strongly depend on local environments for all galaxies with different luminosities. This shows that M is an characteristic parameter only for the environmental dependence of galaxy luminosity.

In order to explore the statistical properties of galaxy morphological types in compact groups (CGs), we construct a random group sample which has the same distributions of redshift and number of member galaxies as those of the CG sample. It turns out that the proportion of early-type galaxies in different redshift bins for the CG sample is statistically higher than that for random group sample, and with growing redshift z this kind of difference becomes more significant. This may be due to the existence of interactions and mergers within a significant fraction of SDSS CGs. We also compare statistical results of CGs with those of more compact groups and pairs, but do not observe as large statistical difference as Hickson (1982)'results.

Using a volume-limited sample of Main Galaxies from the SDSS Data Release 5, we investigate the dependence of galaxy properties on local environment. For each galaxy, the local three-dimensional density is calculated. We find that galaxy morphologies strongly depend on local environment: galaxies in dense environments have predominantly early type morphologies, but other galaxy properties do not present significant dependence on local environment. Clearly, this puts a important constraint on proposed physical mechanisms.