We present a new interatomic potential for solids and liquids called Spectral Neighbor Analysis Potential (SNAP). The SNAP potential has a very general form and uses machine-learning techniques to ...reproduce the energies, forces, and stress tensors of a large set of small configurations of atoms, which are obtained using high-accuracy quantum electronic structure (QM) calculations. The local environment of each atom is characterized by a set of bispectrum components of the local neighbor density projected onto a basis of hyperspherical harmonics in four dimensions. The bispectrum components are the same bond-orientational order parameters employed by the GAP potential 1. The SNAP potential, unlike GAP, assumes a linear relationship between atom energy and bispectrum components. The linear SNAP coefficients are determined using weighted least-squares linear regression against the full QM training set. This allows the SNAP potential to be fit in a robust, automated manner to large QM data sets using many bispectrum components. The calculation of the bispectrum components and the SNAP potential are implemented in the LAMMPS parallel molecular dynamics code. We demonstrate that a previously unnoticed symmetry property can be exploited to reduce the computational cost of the force calculations by more than one order of magnitude. We present results for a SNAP potential for tantalum, showing that it accurately reproduces a range of commonly calculated properties of both the crystalline solid and the liquid phases. In addition, unlike simpler existing potentials, SNAP correctly predicts the energy barrier for screw dislocation migration in BCC tantalum.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
ABSTRACT
Detecting a signal from the Epoch of Reionization (EoR) requires an exquisite understanding of Galactic and extragalactic foregrounds, low-frequency radio instruments, instrumental ...calibration, and data analysis pipelines. In this work, we build upon existing work that aims to understand the impact of calibration errors on 21-cm power spectrum (PS) measurements. It is well established that calibration errors have the potential to inhibit EoR detections by introducing additional spectral features that mimic the structure of EoR signals. We present a straightforward way to estimate the impact of a wide variety of modelling residuals in EoR PS estimation. We apply this framework to the specific case of broken dipoles in Murchison Widefield Array (MWA) to understand its effect and estimate its impact on PS estimation. Combining an estimate of the percentage of MWA tiles that have at least one broken dipole (15–40 per cent) with an analytic description of beam errors induced by such dipoles, we compute the residuals of the foregrounds after calibration and source subtraction. We find that that incorrect beam modelling introduces bias in the 2D-PS on the order of $\sim 10^3\, \mathrm{mK}^2 \, h^{-3}\, \mathrm{Mpc}^{3}$. Although this is three orders of magnitude lower than current lowest limits, it is two orders of magnitude higher than the expected signal. Determining the accuracy of both current beam models and direction-dependent calibration pipelines is therefore crucial in our search for an EoR signal.
We present a survey for transient and variable sources, on time-scales from 28 s to ∼1 yr, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no ...transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of <4.1 × 10−7 deg−2 for the shortest time-scale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest time-scale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of <82 FRBs per sky per day for dispersion measures <700 pc cm−3. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities ≳0.5 Jy and identify no sources with significant variability in their light curves. However, we note that substantial further work is required to fully characterize both the short-term and low-level variability within this field.
The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and ...redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hr of data targeting one of the MWA's epoch of reionization (EoR) fields observed in 2016. We use both improved analysis techniques presented in Barry et al. and several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8, and 7.1 based on our deep analysis on this 40 hr data set. These limits span a range in k-space of 0.18 h Mpc−1 < k < 1.6 h Mpc−1, with a lowest measurement of Δ2 ≤ 2.39 × 103 mK2 at k = 0.59 h Mpc−1 and z = 6.5.
Abstract
We detail new techniques for analysing ionospheric activity, using Epoch of Reionization data sets obtained with the Murchison Widefield Array, calibrated by the ‘real-time system’ (RTS). ...Using the high spatial- and temporal-resolution information of the ionosphere provided by the RTS calibration solutions over 19 nights of observing, we find four distinct types of ionospheric activity, and have developed a metric to provide an ‘at a glance’ value for data quality under differing ionospheric conditions. For each ionospheric type, we analyse variations of this metric as we reduce the number of pierce points, revealing that a modest number of pierce points is required to identify the intensity of ionospheric activity; it is possible to calibrate in real-time, providing continuous information of the phase screen. We also analyse temporal correlations, determine diffractive scales, examine the relative fractions of time occupied by various types of ionospheric activity and detail a method to reconstruct the total electron content responsible for the ionospheric data we observe. These techniques have been developed to be instrument agnostic, useful for application on LOw Frequency ARray and Square Kilometre Array-Low.
We present a sophisticated statistical point-source foreground model for low-frequency radio Epoch of Reionization (EoR) experiments using the 21 cm neutral hydrogen emission line. Motivated by our ...understanding of the low-frequency radio sky, we enhance the realism of two model components compared with existing models: the source count distributions as a function of flux density and spatial position (source clustering), extending current formalisms for the foreground covariance of 2D power-spectral modes in 21 cm EoR experiments. The former we generalize to an arbitrarily broken power law, and the latter to an arbitrary isotropically correlated field. This paper presents expressions for the modified covariance under these extensions, and shows that for a more realistic source spatial distribution, extra covariance arises in the EoR window that was previously unaccounted for. Failure to include this contribution can yield bias in the final power-spectrum and under-estimate uncertainties, potentially leading to a false detection of signal. The extent of this effect is uncertain, owing to ignorance of physical model parameters, but we show that it is dependent on the relative abundance of faint sources, to the effect that our extension will become more important for future deep surveys. Finally, we show that under some parameter choices, ignoring source clustering can lead to false detections on large scales, due to both the induced bias and an artificial reduction in the estimated measurement uncertainty.
ABSTRACT Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization (EoR) and estimation of its basic physical parameters are principal scientific aims of many current ...low-frequency radio telescopes. Here we describe the Cosmological H i Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array, to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximizing the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hr of data, we set a 2 upper limit on the EoR dimensionless power at Mpc−1 of mK2 in the redshift range z = 6.2-6.6, consistent with previous estimates.
ABSTRACT
Observations in the lowest Murchison Widefield Array (MWA) band between 75 and 100 MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at ...redshift ∼13–17. Using 15 h of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction-independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction-dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 h of the cleanest data. Using these data, we obtain 2σ upper limits on the 21 cm power spectrum in the range of $0.1~ h~{\mathrm{ Mpc}}^{-1}\lessapprox k \lessapprox 1 ~ ~h~{\mathrm{ Mpc}}^{-1}$ and at z = 14.2, 15.2, and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $k=0.14 ~h~{\mathrm{ Mpc}}^{-1}$ and at z = 15.2 with a possibility of a few per cent of signal loss due to direction-independent calibration.
The 2D power spectrum is a cornerstone of the modern toolkit for analysis of the low-frequency radio interferometric observations of the 21 cm signal arising from the early universe. Its familiar ...form disentangles a great deal of systematic information concerning both the sky and telescope and is displayed as a foreground-dominated "brick" and "wedge" on large line-of-sight scales and a complementary "window" on smaller scales. This paper builds on many previous works in the literature that seek to elucidate the varied instrumental and foreground factors that contribute to these familiar structures in the 2D power spectrum. In particular, we consider the effects of uv sampling on the emergence of the wedge. Our results verify the expectation that arbitrarily dense instrument layouts in principle restore the missing information that leads to mode mixing and can therefore mitigate the wedge feature. We derive rule-of-thumb estimates for the required baseline density for complete wedge mitigation, showing that these will be unachievable in practice. We also discuss the optimal shape of the layout, showing that logarithmic regularity in the radial separation of baselines is favorable. While complete suppression of foreground leakage into the wedge is practically unachievable, we find that designing layouts that promote radial density and regularity is able to reduce the amplitude of foreground power by one to three orders of magnitude.
Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect ...of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the Murchison Widefield Array (MWA) EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by point spread function sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source.