We describe an improved technique for using the backscattered phase from meteor radar echo measurements just prior to the specular point ($t_{0}$) to calculate meteor speeds and their uncertainty. Our method, which builds on earlier work of Cervera et al (1997), scans possible speeds in the Fresnel distance - time domain with a dynamic, sliding window and derives a best-speed estimate from the resultant speed distribution. We test the performance of our method, called pre-$t_{0}$ speeds by sliding-slopes technique (PSSST), on transverse scattered meteor echoes observed by the Middle Atmosphere Alomar Radar System (MAARSY) and the Canadian Meteor Orbit Radar (CMOR), and compare the results to time-of-flight and Fresnel transform speed estimates. Our novel technique is shown to produce good results when compared to both model and speed measurements using other techniques. We show that our speed precision is $\pm$5$\%$ at speeds less than 40 km/s and we find that more than 90$\%$ of all CMOR multi-station echoes have PSSST solutions. For CMOR data, PSSST is robust against the selection of critical phase value and poor phase unwrapping. Pick errors of up to $\pm$6 pulses for meteor speeds less than about 50 km/s produce errors of less than $\pm$5$\%$ of the meteoroid speed. In addition, the width of the PSSST speed Kernel density estimate (KDE) is used as a natural measure of uncertainty that captures both noise and $t_0$ pick uncertainties.

We perform a self-calibration of the richness-to-mass ($N$-$M$) relation of CAMIRA galaxy clusters with richness $N\geq15$ at redshift $0.2\leq z<1.1$ by modeling redshift-space two-point correlation functions, namely, the auto-correlation function $\xi_{\mathrm{cc}}$ of CAMIRA clusters, the auto-correlation function $\xi_{\mathrm{gg}}$ of the CMASS galaxies spectroscopically observed in the BOSS survey, and the cross-correlation function $\xi_{\mathrm{cg}}$ between these two samples. We focus on constraining the normalization $A_{\mathrm{N}}$ of the $N$-$M$ relation in a forward-modeling approach, carefully accounting for the redshift-space distortion, the Finger-of-God effect, and the uncertainty in photometric redshifts of CAMIRA clusters. The modeling also takes into account the projection effect on the halo bias of CAMIRA clusters. The parameter constraints are shown to be unbiased according to validation tests using a large set of mock catalogs constructed from N-body simulations. At the pivotal mass $M_{500}=10^{14}h^{-1}M_{\odot}$ and the pivotal redshift $z_{\mathrm{piv}} = 0.6$, the resulting normalization $A_{\mathrm{N}}$ is constrained as $13.8^{+5.8}_{-4.2}$, $13.2^{+3.4}_{-2.7}$, and $11.9^{+3.0}_{-1.9}$ by modeling $\xi_{\mathrm{cc}}$, $\xi_{\mathrm{cc}}+\xi_{\mathrm{cg}}$, and $\xi_{\mathrm{cc}} + \xi_{\mathrm{cg}} + \xi_{\mathrm{gg}}$, with average uncertainties at levels of $36\%$, $23\%$, and $21\%$, respectively. We find that the resulting $A_{\mathrm{N}}$ is statistically consistent with those independently obtained from weak-lensing magnification and from a joint analysis of shear and cluster abundance, with a preference for a lower value at a level of $\lesssim1.9\sigma$. This implies that the absolute mass scale of CAMIRA clusters inferred from clustering is mildly higher than those from the independent methods. [abridged]

This paper reports on a cavity haloscope search for dark matter axions in the galactic halo in the mass range $2.81$-$3.31$ ${\mu}eV$. This search excludes the full range of axion-photon coupling values predicted in benchmark models of the invisible axion that solve the strong CP problem of quantum chromodynamics, and marks the first time a haloscope search has been able to search for axions at mode crossings using an alternate cavity configuration. Unprecedented sensitivity in this higher mass range is achieved by deploying an ultra low-noise Josephson parametric amplifier as the first stage signal amplifier.

We show that peculiar velocities of Type Ia supernovae can be used to derive constraints on the sum of neutrino masses, $\Sigma m_{\nu}$, and dark energy equation of state, $w = w_0+w_a(1-a)$, from measurements of the magnitude-redshift relation, complementary to galaxy redshift and weak lensing surveys. Light from a supernova propagates through a perturbed Universe so the luminosity distance is modified from its homogeneous prediction. This modification is proportional to the matter density fluctuation and its time derivative due to gravitational lensing and peculiar velocity respectively. At low redshifts, the peculiar velocity signal dominates while at high redshifts lensing does. We show that using lensing and peculiar velocity of supernovae from the upcoming surveys WFIRST and ZTF, without other observations, we can constrain $\Sigma m_{\nu} \lesssim 0.31$ eV, $\sigma(w_0) \lesssim 0.02$, and ${\sigma(w_a)} \lesssim 0.18$ ($1-\sigma$ CL) in the $\Sigma m_{\nu}$-$w_0$-$w_a$ parameter space, where all the other cosmological parameters are fixed. We find that adding peculiar velocity information from low redshifts shrinks the volume of the parameter ellipsoid in this space by $\sim 33$%. We also allow $\Omega_{\text{CDM}}$ to vary as well as $\Sigma m_{\nu}$, $w_0$ and $w_a$, and demonstrate how these constraints degrade as a consequence.

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others.

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

The $\mu$eV axion is a well-motivated extension to the standard model. The Axion Dark Matter eXperiment (ADMX) collaboration seeks to discover this particle by looking for the resonant conversion of dark-matter axions to microwave photons in a strong magnetic field. In this Letter, we report results from a pathfinder experiment, the ADMX "Sidecar," which is designed to pave the way for future, higher mass, searches. This testbed experiment lives inside of and operates in tandem with the main ADMX experiment. The Sidecar experiment excludes masses in three widely spaced frequency ranges (4202-4249, 5086-5799, and 7173-7203 MHz). In addition, Sidecar demonstrates the successful use of a piezoelectric actuator for cavity tuning. Finally, this publication is the first to report data measured using both the TM$_{010}$ and TM$_{020}$ modes.

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method.

We calculate the bispectrum of scale-invariant tensor modes sourced by spectator SU(2) gauge fields during inflation in a model containing a scalar inflaton, a pseudoscalar axion and SU(2) gauge fields. A large bispectrum is generated in this model at tree-level as the gauge fields contain a tensor degree of freedom, and its production is dominated by self-coupling of the gauge fields. This is a unique feature of non-Abelian gauge theory. The shape of the tensor bispectrum is approximately an equilateral shape for $3\lesssim m_Q\lesssim 4$, where $m_Q$ is an effective dimensionless mass of the SU(2) field normalised by the Hubble expansion rate during inflation. The amplitude of non-Gaussianity of the tensor modes, characterised by the ratio $B_h/P^2_h$, is inversely proportional to the energy density fraction of the gauge field. This ratio can be much greater than unity, whereas the ratio from the vacuum fluctuation of the metric is of order unity. The bispectrum is effective at constraining large $m_Q$ regions of the parameter space, whereas the power spectrum constrains small $m_Q$ regions.

We show that an inflation model in which a spectator axion field is coupled to an SU(2) gauge field produces a large three-point function (bispectrum) of primordial gravitational waves, $B_{h}$, on the scales relevant to the cosmic microwave background experiments. The amplitude of the bispectrum at the equilateral configuration is characterized by $B_{h}/P_h^2=\mathcal{O}(10)\times \Omega_A^{-1}$, where $\Omega_A$ is a fraction of the energy density in the gauge field and $P_h$ is the power spectrum of gravitational waves produced by the gauge field.

We demonstrate equivalence of the in-in formalism and Green's function method for calculating the bispectrum of primordial gravitational waves generated by vacuum fluctuations of the metric. The tree-level bispectrum from the field equation, $B_h$, agrees with the results obtained previously using the in-in formalism exactly. Characterising non-Gaussianity of the fluctuations using the ratio $B_h/P^2_h$ in the equilateral configuration, where $P_h$ is the power spectrum of scale-invariant gravitational waves, we show that it is much weaker than in models with spectator gauge fields. We also calculate the tree-level bispectrum of two right-handed and one left-handed gravitational wave using Green's function, reproducing the results from in-in formalism, and show that it can be as large as the bispectrum of three right-handed gravitational waves.

We present a public code to generate a mock galaxy catalog in redshift space assuming a log-normal probability density function (PDF) of galaxy and matter density fields. We draw galaxies by Poisson-sampling the log-normal field, and calculate the velocity field from the linearised continuity equation of matter fields, assuming zero vorticity. This procedure yields a PDF of the pairwise velocity fields that is qualitatively similar to that of N-body simulations. We check fidelity of the catalog, showing that the measured two-point correlation function and power spectrum in real space agree with the input precisely. We find that a linear bias relation in the power spectrum does not guarantee a linear bias relation in the density contrasts, leading to a cross-correlation coefficient of matter and galaxies deviating from unity on small scales. We also find that linearising the Jacobian of the real-to-redshift space mapping provides a poor model for the two-point statistics in redshift space. That is, non-linear redshift-space distortion is dominated by non-linearity in the Jacobian. The power spectrum in redshift space shows a damping on small scales that is qualitatively similar to that of the well-known Fingers-of-God (FoG) effect due to random velocities, except that the log-normal mock does not include random velocities. This damping is a consequence of non-linearity in the Jacobian, and thus attributing the damping of the power spectrum solely to FoG, as commonly done in the literature, is misleading.

This paper presents SunPy (version 0.5), a community-developed Python package for solar physics. Python, a free, cross-platform, general-purpose, high-level programming language, has seen widespread adoption among the scientific community, resulting in the availability of a large number of software packages, from numerical computation (NumPy, SciPy) and machine learning (scikit-learn) to visualisation and plotting (matplotlib). SunPy is a data-analysis environment specialising in providing the software necessary to analyse solar and heliospheric data in Python. SunPy is open-source software (BSD licence) and has an open and transparent development workflow that anyone can contribute to. SunPy provides access to solar data through integration with the Virtual Solar Observatory (VSO), the Heliophysics Event Knowledgebase (HEK), and the HELiophysics Integrated Observatory (HELIO) webservices. It currently supports image data from major solar missions (e.g., SDO, SOHO, STEREO, and IRIS), time-series data from missions such as GOES, SDO/EVE, and PROBA2/LYRA, and radio spectra from e-Callisto and STEREO/SWAVES. We describe SunPy's functionality, provide examples of solar data analysis in SunPy, and show how Python-based solar data-analysis can leverage the many existing tools already available in Python. We discuss the future goals of the project and encourage interested users to become involved in the planning and development of SunPy.