• We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the S\'{e}rsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.
  • Galaxy clusters are known to host a variety of extended radio sources: tailed radio galaxies whose shape is modelled by the interaction with the intra-cluster medium (ICM); radio bubbles filling cavities in the ICM distribution and rising buoyantly through the thermal gas; diffuse giant radio sources ("halos" and "relics") revealing the presence of relativistic electrons and magnetic fields in the intra-cluster volume. It is currently the subject of an active debate how the non-thermal components that we observe at radio wavelengths affect the physical properties of the ICM and depend on the dynamical state of galaxy clusters. In this work we start our SKA1 feasibility study of the "radio cluster zoo" through simulations of a typical radio-loud cluster, hosting several bright tailed radio galaxies and a diffuse radio halo. Realistic simulations of SKA1 observations are obtained through the MeqTrees software. A new deconvolution algorithm, based on sparse representations and optimised for the detection of faint diffuse astronomical sources, is tested and compared to the classical CLEAN method.
  • GPU-based beamforming is a relatively unexplored area in radio astronomy, possibly due to the assumption that any such system will be severely limited by the PCIe bandwidth required to transfer data to the GPU. We have developed a CUDA-based GPU implementation of a coherent beamformer, specifically designed and optimised for deployment at the BEST-2 array which can generate an arbitrary number of synthesized beams for a wide range of parameters. It achieves $\sim$1.3 TFLOPs on an NVIDIA Tesla K20, approximately 10x faster than an optimised, multithreaded CPU implementation. This kernel has been integrated into two real-time, GPU-based time-domain software pipelines deployed at the BEST-2 array in Medicina: a standalone beamforming pipeline and a transient detection pipeline. We present performance benchmarks for the beamforming kernel as well as the transient detection pipeline with beamforming capabilities as well as results of test observation.
  • We report on the MIT Epoch of Reionization (MITEoR) experiment, a pathfinder low-frequency radio interferometer whose goal is to test technologies that improve the calibration precision and reduce the cost of the high-sensitivity 3D mapping required for 21 cm cosmology. MITEoR accomplishes this by using massive baseline redundancy, which enables both automated precision calibration and correlator cost reduction. We demonstrate and quantify the power and robustness of redundancy for scalability and precision. We find that the calibration parameters precisely describe the effect of the instrument upon our measurements, allowing us to form a model that is consistent with $\chi^2$ per degree of freedom < 1.2 for as much as 80% of the observations. We use these results to develop an optimal estimator of calibration parameters using Wiener filtering, and explore the question of how often and how finely in frequency visibilities must be reliably measured to solve for calibration coefficients. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious Hydrogen Epoch of Reionization Array (HERA) project and other next-generation instruments, which would incorporate many identical or similar technologies.
  • In this memo we investigate the applicability of NVIDIA Graphics Processing Units (GPUs) for SKA1-Low station and Central Signal Processing (CSP)-level processing. Station-level processing primarily involves generating a single station beam which will then be correlated with other beams in CSP. Fine channelisation can be performed either at the station of CSP-level, while coarse channelisation is assumed to be performed on FPGA-based Tile Processors, together with A/D conversion, equilisation and other processes. Rough estimates for number of GPUs required and power requirements will also be provided.
  • Conformal Weyl Gravity (CWG) predicts galactic rotation curves without invoking dark matter. This makes CWG an interesting candidate theory as an alternative to GR. This removal of the necessity of dark matter arises because the exact exterior solution in CWG for a static, spherically symmetric source provides a linear potential $\gamma r$, which can explain the observed galactic rotation curves with a value for $\gamma$ given by $\sim + 10^{-26} \mathrm{m}^{-1}$. Previous work has also shown that CWG can explain lensing observations, but with $\gamma\sim - 10^{-26} \mathrm{m}^{-1}$ in order to ensure converging rays of light rather than diverging ones. Even though different expressions for the angle of deflection have been derived in CWG in order to remove this inconsistency, this paper shows that the $\gamma$ required to fit the lensing data is several orders of magnitude higher than that required to fit the galactic rotation curves.
  • A new digital backend has been developed for the BEST-2 array at Radiotelescopi di Medicina, INAF-IRA, Italy which allows concurrent operation of an FX correlator, and a direct-imaging correlator and beamformer. This backend serves as a platform for testing some of the spatial Fourier transform concepts which have been proposed for use in computing correlations on regularly gridded arrays. While spatial Fourier transform-based beamformers have been implemented previously, this is to our knowledge, the first time a direct-imaging correlator has been deployed on a radio astronomy array. Concurrent observations with the FX and direct-imaging correlator allows for direct comparison between the two architectures. Additionally, we show the potential of the direct-imaging correlator for time-domain astronomy, by passing a subset of beams though a pulsar and transient detection pipeline. These results provide a timely verification for spatial Fourier transform-based instruments that are currently in commissioning. These instruments aim to detect highly-redshifted hydrogen from the Epoch of Reionization and/or to perform widefield surveys for time-domain studies of the radio sky. We experimentally show the direct-imaging correlator architecture to be a viable solution for correlation and beamforming.
  • Gravitomagnetic effects are characterized by two phenomena: first, the geodetic effect which describes the precession of the spin of a gyroscope in a free orbit around a massive object, second, the Lense-Thirring effect which describes the precession of the orbital plane about a rotating source mass. We calculate both these effects in the fourth-order theory of conformal Weyl gravity for the test case of circular orbits. We show that for the geodetic effect a linear term arises which may be interesting for high radial orbits, whereas for the Lense-Thirring effect the additional term has a diminishing effect for most orbits. Circular orbits are also considered in general leading up to a generalization of Kepler's third law.
  • Mapping our universe in 3D by imaging the redshifted 21 cm line from neutral hydrogen has the potential to overtake the cosmic microwave background as our most powerful cosmological probe, because it can map a much larger volume of our Universe, shedding new light on the epoch of reionization, inflation, dark matter, dark energy, and neutrino masses. We report on MITEoR, a pathfinder low-frequency radio interferometer whose goal is to test technologies that greatly reduce the cost of such 3D mapping for a given sensitivity. MITEoR accomplishes this by using massive baseline redundancy both to enable automated precision calibration and to cut the correlator cost scaling from N^2 to NlogN, where N is the number of antennas. The success of MITEoR with its 64 dual-polarization elements bodes well for the more ambitious HERA project, which would incorporate many identical or similar technologies using an order of magnitude more antennas, each with dramatically larger collecting area.
  • Radio transient discovery using next generation radio telescopes will pose several digital signal processing and data transfer challenges, requiring specialized high-performance backends. Several accelerator technologies are being considered as prototyping platforms, including Graphics Processing Units (GPUs). In this paper we present a real-time pipeline prototype capable of processing multiple beams concurrently, performing Radio Frequency Interference (RFI) rejection through thresholding, correcting for the delay in signal arrival times across the frequency band using brute-force dedispersion, event detection and clustering, and finally candidate filtering, with the capability of persisting data buffers containing interesting signals to disk. This setup was deployed at the BEST-2 SKA pathfinder in Medicina, Italy, where several benchmarks and test observations of astrophysical transients were conducted. These tests show that on the deployed hardware eight 20MHz beams can be processed simultaneously for $\sim$640 Dispersion Measure (DM) values. Furthermore, the clustering and candidate filtering algorithms employed prove to be good candidates for online event detection techniques. The number of beams which can be processed increases proportionally to the number of servers deployed and number of GPUs, making it a viable architecture for current and future radio telescopes.
  • Considering cylindrical topology we present the static solution for a charged black hole in conformal gravity. We show that unlike the general relativistic case there are two different solutions, both including a factor that when set to zero recovers the familiar static charged black string solution in Einstein's theory. This factor gives rise to a linear term in the potential that also features in the neutral case and may have significant ramifications for particle trajectories.
  • We present the exact exterior solution for a static and neutral cylindrically symmetric source in locally conformal invariant Weyl gravity. As a special case the general relativity analogue still can be attained, however only as a sub-family of solutions. Our solution contains a linear term that would thus result in a potential that grows linearly over large distances. This may have implications for exotic astrophysical structures as well as matter fields on the extremely small scale.
  • The formation of cosmological black holes is investigated using the functional Schroedinger equation as observed by an asymptotic observer, assuming a spherical domain wall collapse process. The mass formula of the Sultana-Dyer black hole is derived using Israel's domain wall mathematical framework. This is used to examine the semi-classical and quantum nature of the collapsing domain wall in the general scale of an arbitrary scale factor, while ignoring evaporation and back reaction mass losses. Particular FLRW scale factors are then explored for black hole horizon formation times, all yielding the same over-all classical result, namely that an arbitrary amount of time is required for an asymptotic observer to register the formation of this surface.
  • The Square Kilometre Array (SKA) is a radio telescope designed to operate between 70MHz and 10GHz. Due to this large bandwidth, the SKA will be built out of different collectors, namely antennas and dishes to cover the frequency range adequately. In order to deal with this bandwidth, innovative feeds and detectors must be designed and introduced in the initial phases of development. Moreover, the required level of resolution may only be achieved through a groundbreaking configuration of dishes and antennas. Due to the large collecting area and the specifications required for the SKA to deliver the promised science, the configuration of the dishes and the antennas within stations is an important question. This research builds on the work done before by Cohanim et al. (2004), Hassan et al. (2005) and Grigorescu et al. (2009) to further investigate the applicability of machine learning techniques to determine the optimum configurations for the collecting elements within the SKA. This work primarily uses genetic algorithms to search a large space of optimum layouts. Every genetic step provides a population with candidate individuals each of which encodes a possible solution. These are randomly generated or created through the combination of previous encodings. In this study, a number of fitness functions that rank individuals within a population of dish configurations are investigated. The UV density, connecting wire length and power spectra are considered to determine a good dish layout.
  • The identification, and subsequent discovery, of fast radio transients through blind-search surveys requires a large amount of processing power, in worst cases scaling as $\mathcal{O}(N^3)$. For this reason, survey data are generally processed offline, using high-performance computing architectures or hardware-based designs. In recent years, graphics processing units have been extensively used for numerical analysis and scientific simulations, especially after the introduction of new high-level application programming interfaces. Here we show how GPUs can be used for fast transient discovery in real-time. We present a solution to the problem of de-dispersion, providing performance comparisons with a typical computing machine and traditional pulsar processing software. We describe the architecture of a real-time, GPU-based transient search machine. In terms of performance, our GPU solution provides a speed-up factor of between 50 and 200, depending on the parameters of the search.
  • It has recently been pointed out that arbitrary center-of-mass energies may be obtained for particle collisions near the horizon of an extremal Kerr black hole. We investigate this mechanism in cylindrical topology. In particular we consider the center-of-mass energies of a cylindrical black hole with an extremal rotation and charge parameter. The geodesics are first derived with a rotating charged cylindrical black hole producing the background gravitational field. Finally the center-of-mass is determined for this background and its extremal limit is taken.
  • Using f (R) gravity in the Palatini formularism, the metric for a charged spherically symmetric black hole is derived, taking the Ricci scalar curvature to be constant. The generalized uncertainty principle is then used to calculate the temperature of the resulting black hole, through this the entropy is found correcting the Bekenstein-Hawking entropy in this case. Using the entropy the tunneling probability and heat capacity are calculated up to the order of the Planck length, which produces an extra factor that becomes important as black holes become small, such as in the case of mini black holes.
  • The authors apologise for the preliminary submission of this paper.
  • The aim of this paper is to present a new fast-convergent numerically stable space-time adaptive processing (STAP) algorithm derived using a novel technique of feedback orthogonalization. The main advantages of this approach lie in its perfected stability to computational errors and faults which makes its real-time implementation on substantially faster and cheaper regular fixed-point processors possible.
  • The EPR correlation has become an integral part of quantum communications as has general relativity in classical communication theory, however when combined an apparent deterioration is observed for spin states. We consider appropriate changes in directions of measurement to exploit full EPR entanglement for a pair of particles and show that it can be deduced only up to the outer even horizon of a Kerr-Newman black hole, even in the case of freely falling observer.
  • In this work, decision tree learning algorithms and fuzzy inferencing systems are applied for galaxy morphology classification. In particular, the CART, the C4.5, the Random Forest and fuzzy logic algorithms are studied and reliable classifiers are developed to distinguish between spiral galaxies, elliptical galaxies or star/unknown galactic objects. Morphology information for the training and testing datasets is obtained from the Galaxy Zoo project while the corresponding photometric and spectra parameters are downloaded from the SDSS DR7 catalogue.
  • This paper focuses on utilizing two different Bayesian methods to deal with a variety of toy problems which occur in data analysis. In particular we implement the Variational Bayesian and Nested Sampling methods to tackle the problems of polynomial selection and Gaussian Mixture Models, comparing the algorithms in terms of processing speed and accuracy. In the problems tackled here it is the Variational Bayesian algorithms which are the faster though both results give similar results.
  • Densely-packed, all-digital aperture arrays form a key area of technology development required for the Square Kilometre Array (SKA) radio telescope. The design of real-time signal processing systems for digital aperture arrays is currently a central challenge in pathfinder projects worldwide. We describe interim results of such work; an heirarchical, frequency-domain beamforming architecture for synthesising a sky beam from the wideband antenna feeds of digital aperture arrays.