A 150 MHz all sky survey with the Precision Array to Probe the Epoch of Reionization
- Authors: Chege, James Kariuki
- Date: 2020
- Subjects: Epoch of reionization -- Research , Astronomy -- Observations , Radio interferometers
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/117733 , vital:34556
- Description: The Precision Array to Probe the Epoch of Reionization (PAPER) was built to measure the redshifted 21 cm line of hydrogen from cosmic reionization. Such low frequency observations promise to be the best means of understanding the cosmic dawn; when the first galaxies in the universe formed, and also the Epoch of Reionization; when the intergalactic medium changed from neutral to ionized. The major challenges to these observations is the presence of astrophysical foregrounds that are much brighter than the cosmological signal. Here, I present an all-sky survey at 150 MHz obtained from the analysis of 300 hours of PAPER observations. Particular focus is given to the calibration and imaging techniques that need to deal with the wide field of view of a non-tracking instrument. The survey covers ~ 7000 square degrees of the southern sky. From a sky area of 4400 square degrees out of the total survey area, I extract a catalogue of sources brighter than 4 Jy whose accuracy was tested against the published GLEAM catalogue, leading to a fractional difference rms better than 20%. The catalogue provides an all-sky accurate model of the extragalactic foreground to be used for the calibration of future Epoch of Reionization observations and to be subtracted from the PAPER observations themselves in order to mitigate the foreground contamination.
- Full Text:
- Date Issued: 2020
A Bayesian approach to tilted-ring modelling of galaxies
- Authors: Maina, Eric Kamau
- Date: 2020
- Subjects: Bayesian statistical decision theory , Galaxies , Radio astronomy , TiRiFiC (Tilted Ring Fitting Code) , Neutral hydrogen , Spectroscopic data cubes , Galaxy parametrisation
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/145783 , vital:38466
- Description: The orbits of neutral hydrogen (H I) gas found in most disk galaxies are circular and also exhibit long-lived warps at large radii where the restoring gravitational forces of the inner disk become weak (Spekkens and Giovanelli 2006). These warps make the tilted-ring model an ideal choice for galaxy parametrisation. Analysis software utilizing the tilted-ring-model can be grouped into two and three-dimensional based software. Józsa et al. (2007b) demonstrated that three dimensional based software is better suited for galaxy parametrisation because it is affected by the effect of beam smearing only by increasing the uncertainty of parameters but not with the notorious systematic effects observed for two-dimensional fitting techniques. TiRiFiC, The Tilted Ring Fitting Code (Józsa et al. 2007b), is a software to construct parameterised models of high-resolution data cubes of rotating galaxies. It uses the tilted-ring model, and with that, a combination of some parameters such as surface brightness, position angle, rotation velocity and inclination, to describe galaxies. TiRiFiC works by directly fitting tilted-ring models to spectroscopic data cubes and hence is not affected by beam smearing or line-of-site-effects, e.g. strong warps. Because of that, the method is unavoidable as an analytic method in future Hi surveys. In the current implementation, though, there are several drawbacks. The implemented optimisers search for local solutions in parameter space only, do not quantify correlations between parameters and cannot find errors of single parameters. In theory, these drawbacks can be overcome by using Bayesian statistics, implemented in Multinest (Feroz et al. 2008), as it allows for sampling a posterior distribution irrespective of its multimodal nature resulting in parameter samples that correspond to the maximum in the posterior distribution. These parameter samples can be used as well to quantify correlations and find errors of single parameters. Since this method employs Bayesian statistics, it also allows the user to leverage any prior information they may have on parameter values.
- Full Text:
- Date Issued: 2020
A pilot wide-field VLBI survey of the GOODS-North field
- Authors: Akoto-Danso, Alexander
- Date: 2019
- Subjects: Radio astronomy , Very long baseline interferometry , Radio interometers , Imaging systems in astronomy , Hubble Space Telescope (Spacecraft) -- Observations
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/72296 , vital:30027
- Description: Very Long Baseline Interferometry (VLBI) has significant advantages in disentangling active galactic nuclei (AGN) from star formation, particularly at intermediate to high-redshift due to its high angular resolution and insensitivity to dust. Surveys using VLBI arrays are only just becoming practical over wide areas with numerous developments and innovations (such as multi-phase centre techniques) in observation and data analysis techniques. However, fully automated pipelines for VLBI data analysis are based on old software packages and are unable to incorporate new calibration and imaging algorithms. In this work, the researcher developed a pipeline for VLBI data analysis which integrates a recent wide-field imaging algorithm, RFI excision, and a purpose-built source finding algorithm specifically developed for the 64kx64k wide-field VLBI images. The researcher used this novel pipeline to process 6% (~ 9 arcmin2 of the total 160 arcmin2) of the data from the CANDELS GOODS- North extragalactic field at 1.6 GHz. The milli-arcsec scale images have an average rms of a ~ 10 uJy/beam. Forty four (44) candidate sources were detected, most of which are at sub-mJy flux densities, having brightness temperatures and luminosities of >5x105 K and >6x1021 W Hz-1 respectively. This work demonstrates that automated post-processing pipelines for wide-field, uniform sensitivity VLBI surveys are feasible and indeed made more efficient with new software, wide-field imaging algorithms and more purpose-built source- finders. This broadens the discovery space for future wide-field surveys with upcoming arrays such as the African VLBI Network (AVN), MeerKAT and the Square Kilometre Array (SKA).
- Full Text:
- Date Issued: 2019
A statistical study of travelling ionospheric disturbances over the African-European and American sectors
- Authors: Thaganyana, Golekamang Piet
- Date: 2023-03-31
- Subjects: Uncatalogued
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/422541 , vital:71956 , DOI 10.21504/10962/422543
- Description: This research presents a long-term statistical study of travelling ionospheric disturbances (TIDs) of low- and high-latitude origin over the American and African-European sectors between 2010 and 2018. The TIDs of low latitude origin (hereafter known as poleward TIDs) were studied in both quiet and disturbed conditions, whereas the equatorward TIDs were only studied during quiet conditions. The Kp > 4 and Dst_ -50 nT was used as a criterion for geomagnetic disturbed conditions, while the four geomagnetically quiet days were selected each month based on Kp < 3. Observations of TIDs are made using Global Navigational Satellite Systems (GNSS) total electron content derived data. During quiet conditions, seven and two transhemispheric TIDs were identified over the African-European and American sectors, respectively. The observed TIDs originated from the wintertime hemisphere and propagated into the summertime hemisphere. The horizontal velocity, periods, and horizontal wavelengths of TIDs are in range of cH = 120-274 m/s, 48-80 min and _H = 379-1104 km, respectively. These quiet-time equatorward TIDs have been associated with tertiary gravity waves (GWs) from the dissipation of secondary GWs which are in turn generated from the dissipation of mountain waves (MWs) as a result of excited orographic forcing. The poleward TIDs during geomagnetically quiet conditions over the African and American sectors occur mainly during local daytime. Poleward TIDs were observed mostly in the African-European sector than the American sector. Their horizontal propagation velocities and periods range between 129-280 m/s and 39-70 min over African-European and American sectors. Although the mechanisms responsible for launching quiet-time poleward TIDs have not been established in this study, lower atmospheric processes such as convection systems, sudden stratospheric warming and cold weather fronts may have a role in their generation. During geomagnetic storms in the African sector, almost all poleward TIDs (with the exception of two cases) during the main phase were large-scale with horizontal velocities and periods ranging from 250-503 m/s and 30 min to 2 hours. During recovery phase, poleward TIDs fall under the category of medium scale. In the American sector, the majority of poleward TIDs occurred during the storm's main phase, as opposed to the African-European sector, which experienced a significant number of poleward TIDs during the recovery phase. The periods and horizontal velocities of TIDs range from 45 min-1.5 h and 180-296 m/s during main phase. During the recovery phase, the horizontal velocity and period range from 177-271 m/s and 40-1.5 h, respectively. Overall, it has been shown that statistically, changes in equatorial electrodynamics related to enhanced eastward electric _eld and hence increased equatorial electrojet (vertical E_B drift) correlates highly with the reported poleward TIDs. , Thesis (PhD) -- Faculty of Science, Physics and Electronics, 2023
- Full Text:
- Date Issued: 2023-03-31
A study of why some physic concepts in the South African Physical Science curriculum are poorly understood in order to develop a targeted action-research intervention for Newton’s second law
- Authors: Cobbing, Kathleen Margaret
- Date: 2020
- Subjects: Physics -- Study and teaching (Secondary) -- South Africa , Physics -- Examinations, questions, etc. -- South Africa , Motion -- Study and teaching (Secondary) -- South Africa
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/146903 , vital:38575
- Description: Globally, many students show a poor understanding of concepts in high school physics and lack the necessary problem-solving skills that the course demands. The application of Newton’s second law was found to be particularly problematic through document analysis of South African examination feedback reports, as well as from an analysis of the physics examinations at a pair of well-resourced South African independent schools that follow the Independent Examination Board curriculum. Through an action-research approach, a resource for use by students was designed and modified to improve students’ understanding of this concept, while modelling problemsolving methods. The resource consisted of brief revision notes, worked examples and scaffolded exercises. The design of the resource was influenced by the theory of cognitive apprenticeship, cognitive load theory and conceptual change theory. One of the aims of the resource was to encourage students to translate between the different representations of a problem situation: symbolic, abstract, model and concrete. The impact of this resource was evaluated at a pair of schools using a mixed methods approach. This incorporated pre- and post-tests for a quantitative assessment, qualitative student evaluations and the analysis of examination scripts. There was an improvement from pre- to post-test for all four iterations of the intervention and these improvements were shown to be significant. The use of the resource led to an increase in the quality and quantity of diagrams drawn by students in subsequent assessments.
- Full Text:
- Date Issued: 2020
Accelerated implementations of the RIME for DDE calibration and source modelling
- Authors: Van Staden, Joshua
- Date: 2021
- Subjects: Radio astronomy , Radio inferometers , Radio inferometers -- Calibration , Radio astronomy -- Data processing , Radio inferometers -- Data processing , Radio inferometers -- Calibration -- Data processing
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/172422 , vital:42199
- Description: Second- and third-generation calibration methods filter out subtle effects in interferometer data, and therefore yield significantly higher dynamic ranges. The basis of these calibration techniques relies on building a model of the sky and corrupting it with models of the effects acting on the sources. The sensitivities of modern instruments call for more elaborate models to capture the level of detail that is required to achieve accurate calibration. This thesis implements two types of models to be used in for second- and third-generation calibration. The first model implemented is shapelets, which can be used to model radio source morphologies directly in uv space. The second model implemented is Zernike polynomials, which can be used to represent the primary beam of the antenna. We implement these models in the CODEX-AFRICANUS package and provide a set of unit tests for each model. Additionally, we compare our implementations against other methods of representing these objects and instrumental effects, namely NIFTY-GRIDDER against shapelets and a FITS-interpolation method against the Zernike polynomials. We find that to achieve sufficient accuracy, our implementation of the shapelet model has a higher runtime to that of the NIFTY-GRIDDER. However, the NIFTY-GRIDDER cannot simulate a component-based sky model while the shapelet model can. Additionally, the shapelet model is fully parametric, which allows for integration into a parameterised solver. We find that, while having a smaller memory footprint, our Zernike model has a greater computational complexity than that of the FITS-interpolated method. However, we find that the Zernike implementation has floating-point accuracy in its modelling, while the FITS-interpolated model loses some accuracy through the discretisation of the beam.
- Full Text:
- Date Issued: 2021
Advanced radio interferometric simulation and data reduction techniques
- Authors: Makhathini, Sphesihle
- Date: 2018
- Subjects: Interferometry , Radio interferometers , Algorithms , Radio telescopes , Square Kilometre Array (Project) , Very Large Array (Observatory : N.M.) , Radio astronomy
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/57348 , vital:26875
- Description: This work shows how legacy and novel radio Interferometry software packages and algorithms can be combined to produce high-quality reductions from modern telescopes, as well as end-to-end simulations for upcoming instruments such as the Square Kilometre Array (SKA) and its pathfinders. We first use a MeqTrees based simulations framework to quantify how artefacts due to direction-dependent effects accumulate with time, and the consequences of this accumulation when observing the same field multiple times in order to reach the survey depth. Our simulations suggest that a survey like LADUMA (Looking at the Distant Universe with MeerKAT Array), which aims to achieve its survey depth of 16 µJy/beam in a 72 kHz at 1.42 GHz by observing the same field for 1000 hours, will be able to reach its target depth in the presence of these artefacts. We also present stimela, a system agnostic scripting framework for simulating, processing and imaging radio interferometric data. This framework is then used to write an end-to-end simulation pipeline in order to quantify the resolution and sensitivity of the SKA1-MID telescope (the first phase of the SKA mid-frequency telescope) as a function of frequency, as well as the scale-dependent sensitivity of the telescope. Finally, a stimela-based reduction pipeline is used to process data of the field around the source 3C147, taken by the Karl G. Jansky Very Large Array (VLA). The reconstructed image from this reduction has a typical 1a noise level of 2.87 µJy/beam, and consequently a dynamic range of 8x106:1, given the 22.58 Jy/beam flux Density of the source 3C147.
- Full Text:
- Date Issued: 2018
An investigation of traveling ionospheric disturbances (TIDs) in the SANAE HF radar data
- Authors: Atilaw, Tsige Yared
- Date: 2022-04-07
- Subjects: Ionospheric storms Antarctica , Radar Antarctica , Range time-intensity (RTI) , South African National Antarctic Expedition (SANAE) , Super Dual Auroral Radar Network (SuperDARN)
- Language: English
- Type: Doctoral thesis , text
- Identifier: http://hdl.handle.net/10962/232377 , vital:49986 , DOI 10.21504/10962/232377
- Description: This thesis aims to study the characteristics of traveling ionospheric disturbances (TIDs) as identified in the radar data of the South African National Antarctic Expedition (SANAE) Super Dual Auroral Radar Network (SuperDARN) radar located in Antarctica. For this project, 22 TIDs were identified from visual inspection of range time-intensity (RTI) plots of backscattered power and Doppler velocity parameters of the SANAE radar between 2005âAS2015. These events were studied to determine their characteristics and driving mechanisms. Where good quality data were available, the SANAE HF radar data were supplemented by Halley radar data, which has large area of overlapping field of view (FOV) with the SANAE radar, and also by GPS TEC data. This provided a multi-instrument data analysis of some TID events. Different spectral analysis methods, namely the multitaper method (MTM), Fast Fourier transform (FFT) and the Lomb-Scargle periodogram were used to obtain spectral information of the observed waves. The advantage of using multiple windowing in MTM over the traditional windowing method was illustrated using one of the TID events. In addition, the analytic signal of the wave from the MTM method was used to estimate the instantaneous phase velocity and propagation azimuth of the wave, which was able to track the change in the characteristics of the medium-scale TID (MSTID) efficiently throughout the duration of the event. This is a clear advantage over other windowing techniques. The energy contribution by this MSTID through Joule heating was estimated over the region where spectral analysis of both SANAE and Halley data showed it to be present. The majority of the TIDs (65.4%) could be classified as MSTIDs with periods of 20–60 minutes, velocities of 50–333 ms1 and wavelengths of 129–833 km. The TID occurrence rate was high around the March equinox with 12 out of the 16 event days being during March–May. March had a particularly high number of occurrences of TIDs (46%). The majority of the TIDs observed during this month propagated northward or southeastward. In terms of prevailing geomagnetic conditions, 6 out of 16 event days were geomagnetically quiet, while 10 occurred during geomagnetic storms and substorms. During quiet conditions, TIDs could be linked to Es and polarised electric fields in 2 of these events. The other quiet time events could not be related to Es instability and polarised electric field either because their exact propagation direction could not be determined or data quality from the Es region scatter was too poor to perform spectral analysis. The storm-/substorm-related TIDs are possibly generated through Joule heating, the Lorentz force and energetic particle precipitation. , Thesis (PhD) -- Faculty of Science, Physics and Electronics, 2022
- Full Text:
- Date Issued: 2022-04-07
Behaviour of quiet time ionospheric disturbances at African equatorial and midlatitude regions
- Authors: Orford, Nicola Diane
- Date: 2018
- Subjects: Ionospheric storms , Ionospheric storms -- Africa , Ionosphere , Plasmasphere , Q-disturbances , Total electron content (TEC)
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/62672 , vital:28228
- Description: Extreme ionospheric and geomagnetic disturbances affect technology adversely. Prestorm enhancements, considered a potential predictor of geomagnetic storms, occur during quiet conditions prior to geomagnetic disturbances. The ionosphere experiences general disturbances during quiet geomagnetic conditions and these Q- disturbances remain unexplored over Africa. This study used TEC data to characterize the morphology of Q-disturbances over Africa, exploring variations with solar cycle, season, time of occurrence and latitude. Observations from 10 African GPS stations in the equatorial and midlatitude regions show that Q-disturbances in the equatorial region are predominantly driven by E x B variations, while multiple mechanisms affect the midlatitude region. Q- disturbances occur more frequently during nighttime than during daytime and no seasonal trend is observed. Midlatitude Q-disturbance mechanisms are explored in depth, considering substorm activity, the plasmaspheric contribution to GPS TEC and plasma transfer between conjugate points. Substorm activity is not a dominant mechanism, although Q-disturbances occurring under elevated substorm conditions tend to have longer duration and larger amplitude than general Q-disturbances. Many observed Q-disturbances become non-significant once the plasmaspheric contribution to the TEC measurements is removed, indicating that these disturbances occur within the plasmasphere, and not the ionosphere. Transfer of plasma between conjugate points does not seem to be a mechanism driving Q-disturbances, as the corresponding nighttime behaviour expected between depletions in the summer hemisphere and enhancements in the winter hemisphere is not observed. Pre-storm enhancements occur infrequently, rendering them a poor predictor of geomagnetic disturbances. Pre-storm enhancement morphology does not differ significantly from general quiet time enhancement morphology, suggesting pre-storms are not a special case of Q-disturbances.
- Full Text:
- Date Issued: 2018
Combined spectral and stimulated luminescence study of charge trapping and recombination processes in α-Al2O3:C
- Authors: Nyirenda, Angel Newton
- Date: 2018
- Subjects: Luminescence , Thermoluminescence , Luminescence spectroscopy , Carbon-doped aluminium oxide , Radioluminescence , Time-resolved X-ray excited optical luminescence
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/62683 , vital:28235
- Description: The main objective of this project was to gain a deeper and better understanding of the luminescence processes in a-Al₂O₃:C, a highly-sensitive dosimetric material, using a combined spectral and stimulated luminescence study. The spectral studies concentrated on the emission spectra obtained using X-ray induced radioluminescence (XERL), thermoluminescence (XETL) and time-resolved X-ray excited optical luminescence (TR-XEOL) techniques. The stimulated luminescence studies were based on thermoluminescence (TL), optically stimulated luminescence (OSL) and phototransferred TL (PTTL) methods that were used in the study of the radiation-induced defects at high beta-doses and the deep traps, that is, traps with thermal depths beyond 500°C. The spectral and stimulated luminescence measurements were carried out using a high sensitivity luminescence spectrometer and a Ris0 TL/OSL Model DA-20 Reader, respectively. The XERL emission spectrum measured at room temperature shows seven gaussian peaks associated with F-centres (420 nm), F+-centres (334 nm), F2+-centres (559 nm), Stoke’s vibronic band of Cr3+ (671 nm), Cr3+ R-line emission (694 nm), anti-Stokes vibronic band of Cr3+ (710 nm) and an unidentified emission band (260-300 nm) which we associate with hole recombinations at a luminescence centre. The 694-nm R-line emission from Cr3+ impurity ions is most likely due to recombination of holes at Cr2+ during stimulated luminescence and as a result of an intracentre excitation of Cr3+ in photoluminescence (PL) due to photon absorption. The Cr3+ emission decreases in intensity, whereas the intensity of F-centre emission band is almost constant with repeated XERL measurements. Depending on the amount of X-ray irradiation dose, both holes and/or electrons may take place in the emission processes of peaks I (30-80°C), II (90-250°C) and III (250-320°C) during a TL readout, albeit, electron recombination is dominant regardless of dose. At higher doses, the XETL emission spectra indicate that the dominant band associated with TL peak III (250-320°C) in the material, shifts from F-centre to Cr3+. Using the deep-traps OSL, it has been confirmed that the main TL trap is also the main OSL trap whereas the TL traps lying in the temperature range of 400-550°C constitute the secondary OSL traps. There is evidence of strong retrapping at the main trap during optical stimulation of charges from the secondary OSL traps and the deep traps and that the retrapping occurs via the delocalized bands. At high-irradiation beta-doses, aggregate defect centres which significantly alter the TL and OSL properties, are induced in the material. The induced aggregate centres get completely obliterated by heating a sample to 700°C. The radiation-induced defects cause the main TL peak to shift towards higher temperatures, increase its FWHM, reduce its maximum intensity and cause an underestimation of both the activation energy and order of kinetics of the peak. On the other hand, the OSL response of the material is enhanced following a high-irradiation dose. During sample storage in the dark at ambient temperature, charges do migrate from the deep traps (donors) to the main and intermediate traps (acceptors) and that the major donor traps during this charge transfer phenomenon lie between 500-600°C.
- Full Text:
- Date Issued: 2018
CubiCal: a fast radio interferometric calibration suite exploiting complex optimisation
- Authors: Kenyon, Jonathan
- Date: 2019
- Subjects: Interferometry , Radio astronomy , Python (Computer program language) , Square Kilometre Array (Project)
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/92341 , vital:30711
- Description: The advent of the Square Kilometre Array and its precursors marks the start of an exciting era for radio interferometry. However, with new instruments producing unprecedented quantities of data, many existing calibration algorithms and implementations will be hard-pressed to keep up. Fortunately, it has recently been shown that the radio interferometric calibration problem can be expressed concisely using the ideas of complex optimisation. The resulting framework exposes properties of the calibration problem which can be exploited to accelerate traditional non-linear least squares algorithms. We extend the existing work on the topic by considering the more general problem of calibrating a Jones chain: the product of several unknown gain terms. We also derive specialised solvers for performing phase-only, delay and pointing error calibration. In doing so, we devise a method for determining update rules for arbitrary, real-valued parametrisations of a complex gain. The solvers are implemented in an optimised Python package called CubiCal. CubiCal makes use of Cython to generate fast C and C++ routines for performing computationally demanding tasks whilst leveraging multiprocessing and shared memory to take advantage of modern, parallel hardware. The package is fully compatible with the measurement set, the most common format for interferometer data, and is well integrated with Montblanc - a third party package which implements optimised model visibility prediction. CubiCal's calibration routines are applied successfully to both simulated and real data for the field surrounding source 3C147. These tests include direction-independent and direction dependent calibration, as well as tests of the specialised solvers. Finally, we conduct extensive performance benchmarks and verify that CubiCal convincingly outperforms its most comparable competitor.
- Full Text:
- Date Issued: 2019
Design patterns and software techniques for large-scale, open and reproducible data reduction
- Authors: Molenaar, Gijs Jan
- Date: 2021
- Subjects: Radio astronomy -- Data processing , Radio astronomy -- Data processing -- Software , Radio astronomy -- South Africa , ASTRODECONV2019 dataset , Radio telescopes -- South Africa , KERN (omputer software)
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/172169 , vital:42172 , 10.21504/10962/172169
- Description: The preparation for the construction of the Square Kilometre Array, and the introduction of its operational precursors, such as LOFAR and MeerKAT, mark the beginning of an exciting era for astronomy. Impressive new data containing valuable science just waiting for discovery is already being generated, and these devices will produce far more data than has ever been collected before. However, with every new data instrument, the data rates grow to unprecedented quantities of data, requiring novel new data-processing tools. In addition, creating science grade data from the raw data still requires significant expert knowledge for processing this data. The software used is often developed by a scientist who lacks proper training in software development skills, resulting in the software not progressing beyond a prototype stage in quality. In the first chapter, we explore various organisational and technical approaches to address these issues by providing a historical overview of the development of radioastronomy pipelines since the inception of the field in the 1940s. In that, the steps required to create a radio image are investigated. We used the lessons-learned to identify patterns in the challenges experienced, and the solutions created to address these over the years. The second chapter describes the mathematical foundations that are essential for radio imaging. In the third chapter, we discuss the production of the KERN Linux distribution, which is a set of software packages containing most radio astronomy software currently in use. Considerable effort was put into making sure that the contained software installs appropriately, all items next to one other on the same system. Where required and possible, bugs and portability fixes were solved and reported with the upstream maintainers. The KERN project also has a website, and issue tracker, where users can report bugs and maintainers can coordinate the packaging effort and new releases. The software packages can be used inside Docker and Singularity containers, enabling the installation of these packages on a wide variety of platforms. In the fourth and fifth chapters, we discuss methods and frameworks for combining the available data reduction tools into recomposable pipelines and introduce the Kliko specification and software. This framework was created to enable end-user astronomers to chain and containerise operations of software in KERN packages. Next, we discuss the Common Workflow Language (CommonWL), a similar but more advanced and mature pipeline framework invented by bio-informatics scientists. CommonWL is supported by a wide range of tools already; among other schedulers, visualisers and editors. Consequently, when a pipeline is made with CommonWL, it can be deployed and manipulated with a wide range of tools. In the final chapter, we attempt something unconventional, applying a generative adversarial network based on deep learning techniques to perform the task of sky brightness reconstruction. Since deep learning methods often require a large number of training samples, we constructed a CommonWL simulation pipeline for creating dirty images and corresponding sky models. This simulated dataset has been made publicly available as the ASTRODECONV2019 dataset. It is shown that this method is useful to perform the restoration and matches the performance of a single clean cycle. In addition, we incorporated domain knowledge by adding the point spread function to the network and by utilising a custom loss function during training. Although it was not possible to improve the cleaning performance of commonly used existing tools, the computational time performance of the approach looks very promising. We suggest that a smaller scope should be the starting point for further studies and optimising of the training of the neural network could produce the desired results.
- Full Text:
- Date Issued: 2021
Dynamics of charge movement in ∞-Al2O3:C,Mg using thermoluminescence phototransferred and optically stimulated luminescence
- Authors: Lontsi Sob, Aaron Joel
- Date: 2022-04-08
- Subjects: Thermoluminescence , Optically stimulated luminescence , Phototransfer , Deep traps , Phototransferred thermoluminescence (PTTL)
- Language: English
- Type: Academic theses , Doctoral theses , text
- Identifier: http://hdl.handle.net/10962/294607 , vital:57237 , DOI 10.21504/10962/294607
- Description: The dosimetric features of ∞-Al2O3:C,Mg have been investigated for unannealed and annealed samples. The unannealed sample is referred to as sample A whereas the samples annealed at 700, 900 and 1200°C for 15 minutes each are referred to as samples B, C and D respectively. A glow curve of unannealed ∞-Al2O3:C,Mg measured at 1°C/s after irradiation to 2.0 Gy consists of peaks at 43, 73, 164, 195, 246, 284, 336 and 374°C respectively. For sample B (annealed at 700°C), a glow curve measured at 1°C/s after irradiation to 3.0 Gy has peaks at 46, 76, 100, 170, 199, 290, 330 and 375°C whereas the glow curve of sample C (annealed at 900°C) recorded under the same conditions consists of peaks at 49, 80, 100, 174, 206, 235, 290, 335 and 375°C respectively. Sample D (annealed at 1200°C) is the most sensitive of the four samples. A glow curve of sample D measured at 1°C/s after irradiation to 0.2 Gy has peaks at 52, 82, 102, 174, 234, 288 and 384°C respectively. The peaks are labelled I-VIII in order of appearance. The 100°C peak, labelled IIa, is induced by annealing at or above 700°C. The dose response of these peaks was studied for doses within 0.1-8.2 Gy. The reported peaks follow first-order kinetics irrespective of annealing temperature. Peaks I-III of each sample are reproduced under phototransfer for preheating up to 400°C. For the unannealed sample, the reproduced peaks are labelled A1-A3 whereas for the annealed samples, they are labelled B1-B3, C1-C3 and D1-D3 respectively. The annealing-induced peak at 100°C is reproduced as B2a, C2a and D2a for samples B, C and D respectively. A PTTL peak labelled C2b or D2b is also observed near 140°C in samples C and D. In addition to these PTTL peaks, a PTTL peak corresponding to peak IV is also found for sample D and for the unannealed sample. As the corresponding conventional peaks, the PTTL peaks of each sample follow first-order kinetics. Peak I and its corresponding PTTL peak for each sample are unstable and fade to a minimal level after 300 s of storage time. On the other hand, peak II of each sample and its corresponding PTTL peak could still be observed with delay up to 5000 s. Peak III of the unannealed sample remains stable with storage time up to 48 hours. Irrespective of annealing, the trap corresponding to peak III is the most sensitive to optical stimulation. Time-dependent profiles of PTTL from unannealed and annealed ∞-Al2O3:C,Mg were also studied. The mathematical analysis of the PTTL time-response profiles is based on experimental results. The role of various electron traps in PTTL was determined by using pulse annealing and by monitoring the dependence of peak intensity on duration of illumination for peaks not removed by preheating. The presence and role of deep traps were further demonstrated with thermally assisted optically stimulated luminescence. For the unannealed sample, the activation energy for thermal assistance is 0.033 ± 0.001 eV and the activation energy for thermal i quenching is 1.043 ± 0.001 eV. For sample C, the activation energy for thermal assistance is 0.044 ± 0.003 eV whereas that for thermal quenching is 1.110 ± 0.006 eV. The values for the activation energy for thermal assistance are lower than those reported in literature. Only the values for the activation energy for thermal quenching are somewhat comparable to values reported elsewhere. , Thesis (PhD) -- Faculty of Science, Physics and Electronics, 2022
- Full Text:
- Date Issued: 2022-04-08
Dynamics of stimulated luminescence in natural quartz: Thermoluminescence and phototransferred thermoluminescence
- Authors: Folley, Damilola Esther
- Date: 2020
- Subjects: Thermoluminescence , Quartz
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/146255 , vital:38509
- Description: Natural quartz has remained an important mineral that is of topical interest in luminescence and dosimetry-related research. We investigate the dynamics of stimulated luminescence on this material through thermoluminescence (TL) and phototransferred thermoluminescence (PTTL). Measurements were made on unannealed natural quartz as well as quartz annealed at 800 and 1000̊C. The samples were annealed for 10 minutes and for 1 hour. The material, in its un- and annealed state has its main peak between 68 and 72̊C when measured at 1Cs ̃1 after a dose of 50 Gy. A study of dosimetric features and kinetic analysis was carried out on two prominent peaks, peak I and III for all the samples. The peaks show a sublinear dose response for irradiation doses between 10 and 300 Gy. Kinetic analysis shows that peak I is a first-order peak and peak III a general-order peak. Interestingly, we observe for peak I for the sample annealed at 800̊C for 1 hour an inverse thermal quenching behaviour. We demonstrate that a peak affected with an inverse thermal quenching-like behaviour can still show effect of thermal quenching when the dose the sample is irradiated to is significantly reduced. We ascribe the apparent dependence of thermal quenching on dose to competition between radiative and non-radiative transitions at the recombination centre. Peaks I, II, and III for all the samples were reproduced under phototransfer when the peaks, initially removed by preheating to a certain temperature are exposed to 470 and 525 nm light. The infuence of duration of illumination on the PTTL intensity of these peaks corresponding to various preheating temperatures is modelled using coupled first-order dfferential equations. The model is based on systems of acceptors and donors whose number and role depends on preheating temperature
- Full Text:
- Date Issued: 2020
Finite precision arithmetic in Polyphase Filterbank implementations
- Authors: Myburgh, Talon
- Date: 2020
- Subjects: Radio interferometers , Interferometry , Radio telescopes , Gate array circuits , Floating-point arithmetic , Python (Computer program language) , Polyphase Filterbank , Finite precision arithmetic , MeerKAT
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/146187 , vital:38503
- Description: The MeerKAT is the most sensitive radio telescope in its class, and it is important that systematic effects do not limit the dynamic range of the instrument, preventing this sensitivity from being harnessed for deep integrations. During commissioning, spurious artefacts were noted in the MeerKAT passband and the root cause was attributed to systematic errors in the digital signal path. Finite precision arithmetic used by the Polyphase Filterbank (PFB) was one of the main factors contributing to the spurious responses, together with bugs in the firmware. This thesis describes a software PFB simulator that was built to mimic the MeerKAT PFB and allow investigation into the origin and mitigation of the effects seen on the telescope. This simulator was used to investigate the effects in signal integrity of various rounding techniques, overflow strategies and dual polarisation processing in the PFB. Using the simulator to investigate a number of different signal levels, bit-width and algorithmic scenarios, it gave insight into how the periodic dips occurring in the MeerKAT passband were the result of the implementation using an inappropriate rounding strategy. It further indicated how to select the best strategy for preventing overflow while maintaining high quantization effciency in the FFT. This practice of simulating the design behaviour in the PFB independently of the tools used to design the DSP firmware, is a step towards an end-to-end simulation of the MeerKAT system (or any radio telescope using nite precision digital signal processing systems). This would be useful for design, diagnostics, signal analysis and prototyping of the overall instrument.
- Full Text:
- Date Issued: 2020
Foreground simulations for observations of the global 21-cm signal
- Authors: Klutse, Diana
- Date: 2019
- Subjects: Cosmic background radiation , Astronomy -- Observations , Electromagnetic waves , Radiation, Background
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/76398 , vital:30557
- Description: The sky-averaged (global) spectrum of the redshifted 21-cm line promises to be a direct probe of the Dark Ages, the period before the first luminous sources formed and the Epoch of Reionization during which these sources produced enough ionizing photons to ionize the neutral intergalactic medium. However, observations of this signal are contaminated by both astrophysical foregrounds which are orders of magnitude brighter than the cosmological signal and by non-astrophysical and non-ideal instrumental effects. It is therefore crucial to understand all these data components and their impacts on the cosmological signal, for successful signal extraction. In this view, we investigated the impact that small scale spatial structures of diffuse Galactic foreground has on the foreground spectrum as observed by a global 21-cm observation. We simulated two different sets of observations using a realistic dipole beam model of two synchotron foreground templates that differ from each other in the small scale structure: the original 408 MHz all-sky map by Haslam et al. (1982) and a version where the calibration was improved to remove artifcats and point sources (Remazeilles et al., 2015). We generated simulated foreground spectra and modeled them using a polynomial expansion in frequency. We found that the different foreground templates have a modest impact on the simulated spectra, generate differences up to 2% in the root mean square of residual spectra after the log-polynomial best fit was subtracted out.
- Full Text:
- Date Issued: 2019
Influence of argon ion implantation on the thermoluminescence properties of aluminium oxide
- Authors: Khabo, Bokang
- Date: 2022-04-06
- Subjects: Aluminum oxide , Thermoluminescence , Ion implantation , Kinetic analysis , Oxygen vacancies , Argon , Irradiation
- Language: English
- Type: Master's thesis , text
- Identifier: http://hdl.handle.net/10962/234220 , vital:50173
- Description: The influence of argon ion implantation on the thermoluminescence properties (TL) of aluminium oxide (alumina) was investigated. Aluminium oxide (Al2O3) samples were implanted with 80 keV Ar ions. An unimplanted sample and samples implanted at fluences of 1×1014, 5×1014, 1×1015, 5×1015, 1×1016 Ar+/cm2 were irradiated at a dose of 40 Gy and heated at a rate of 1°C/s using a Risø reader model TL/OSL-DA-20 equipped with a Hoya U-340 filter. The thermoluminescence glow curves showed five distinct peaks with main peaks at 178°C, 188°C, 176°C, 208°C, 216°C and 204°C for the unimplanted sample as well as implanted samples. The peak positions of the samples were independent of the irradiation dose suggesting that the samples were characterised by first order kinetics. This was also confirmed by the TM-TSTOP analysis. It was observed that the TL intensity decreases with fluence of implantation. This observation suggests that the concentration of electron traps responsible for thermoluminescence decreases with ion implantation. The decrease in electron concentration might be due to the formation of non-radiative transition bands or the creation of defect clusters and extended defects following the ion implantation and ion fluence increases. The stopping and range of atoms in matter (SRIM) program was used to correlate the TL response of Al2O3 with defects under ion implantation. Subsequent to ion implantation, it was found that the number of oxygen vacancies which are related to electron traps are higher than the number of aluminium vacancies. Kinetic analysis was carried out using the initial rise, Chens peak shape, various heating rate, the whole glow curve, glow curve fitting and the isothermal decay methods. The activation energy was found to be around 0.8 eV and the frequency factor to be of the order 108 𝑠−1 regardless of the implantation fluence. This means that argon ion implantation did not affect the nature of electron traps. The dosimetric features of samples were also investigated at doses in the range of 40 – 200 Gy. Samples generally showed a superlinear response at doses less than 140 Gy and sublinear response at doses higher than 160 Gy. , Thesis (MSc) -- Faculty of Science, Physics and Electronics, 2022
- Full Text:
- Date Issued: 2022-04-06
Ionospheric disturbances during magnetic storms at SANAE
- Authors: Hiyadutuje, Alicreance
- Date: 2017
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/54956 , vital:26639
- Description: The coronal mass ejections (CMEs) and solar flares associated with extreme solar activity may strike the Earth's magnetosphere and give rise to geomagnetic storms. During geomagnetic storms, the polar plasma dynamics may influence the middle and low-latitude ionosphere via travelling ionospheric disturbances (TIDs). These are wave-like electron density disturbances caused by atmospheric gravity waves propagating in the ionosphere. TIDs focus and defocus SuperDARN signals producing a characteristic pattern of ground backscattered power (Samson et al., 1989). Geomagnetic storms may cause a decrease of total electron content (TEC), i.e. a negative storm effect, or/and an increase of TEC, i.e. a positive storm effect. The aim of this project was to investigate the ionospheric response to strong storms (Dst < -100 nT) between 2011 and 2015, using TEC and scintillation measurements derived from GPS receivers as well as SuperDARN power, Doppler velocity and convection maps. In this study the ionosphere's response to geomagnetic storms is determined by the magnitude and time of occurrence of the geomagnetic storm. The ionospheric TEC results of this study show that most of the storm effects observed were a combination of both negative and positive per storm per station (77.8%), and only 8.9% and 13.3% of effects on TEC were negative and positive respectively. The highest number of storm effects occurred in autumn (36.4%), while 31.6%, 28.4% and 3.6% occurred in winter, spring and summer respectively. During the storms studied, 71.4% had phase scintillation in the range of 0.7 - 1 radians, and only 14.3% of the storms had amplitude scintillations near 0.4. The storms studied at SANAE station generated TIDs with periods of less than an hour and amplitudes in the range 0.2 - 5 TECU. These TIDs were found to originate from the high-velocity plasma flows, some of which are visible in SuperDARN convection maps. Early studies concluded that likely sources of these disturbances correspond to ionospheric current surges (Bristow et al., 1994) in the dayside auroral zone (Huang et al., 1998).
- Full Text:
- Date Issued: 2017
Long-term analysis of ionospheric response during geomagnetic storms in mid, low and equatorial latitudes
- Authors: Matamba, Tshimangadzo Merline
- Date: 2018
- Subjects: Ionospheric storms , Coronal mass ejections , Corotating interaction regions , Solar flares , Global Positioning System , Ionospheric critical frequencies , Equatorial Ionization Anomaly (EIA)
- Language: English
- Type: text , Thesis , Doctoral , PhD
- Identifier: http://hdl.handle.net/10962/63991 , vital:28517
- Description: Understanding changes in the ionosphere is important for High Frequency (HF) communications and navigation systems. Ionospheric storms are the disturbances in the Earth’s upper atmosphere due to solar activities such as Coronal Mass Ejections (CMEs), Corotating interaction Regions (CIRs) and solar flares. This thesis reports for the first time on an investigation of ionospheric response to great geomagnetic storms (Disturbance storm time, Dst ≤ −350 nT) that occurred during solar cycle 23. The storm periods analysed were 29 March - 02 April 2001, 27 - 31 October 2003, 18 - 23 November 2003 and 06 - 11 November 2004. Global Navigation Satellite System (GNSS), Total Electron Content (TEC) and ionosonde critical frequency of F2 layer (foF2) data over northern hemisphere (European sector) and southern hemisphere (African sector) mid-latitudes were used to study the ionospheric responses within 15E° - 40°E longitude and ±31°- ±46° geomagnetic latitude. Mid-latitude regions within the same longitude sector in both hemispheres were selected in order to assess the contribution of the low latitude changes especially the expansion of Equatorial Ionization Anomaly (EIA) also known as the dayside ionospheric super-fountain effect during these storms. In all storm periods, both negative and positive ionospheric responses were observed in both hemispheres. Negative ionospheric responses were mainly due to changes in neutral composition, while the expansion of the EIA led to pronounced positive ionospheric storm effect at mid-latitudes for some storm periods. In other cases (e.g 29 October 2003), Prompt Penetration Electric Fields (PPEF), EIA expansion and large scale Traveling Ionospheric Disturbances (TIDs) were found to be present during the positive storm effect at mid-latitudes in both hemispheres. An increase in TEC on the 28 October 2003 was because of the large solar flare with previously determined intensity of X45± 5. A further report on statistical analysis of ionospheric storm effects due to Corotating Interaction Region (CIR)- and Coronal Mass Ejection (CME)-driven storms was performed. The storm periods analyzed occurred during the period 2001 - 2015 which covers part of solar cycles 23 and 24. Dst≤ -30 nT and Kp≥ 3 indices were used to identify the storm periods considered. Ionospheric TEC derived from IGS stations that lie within 30°E - 40°E geographic longitude in mid, low and equatorial latitude over the African sector were used. The statistical analysis of ionospheric storm effects were compared over mid, low and equatorial latitudes in the African sector for the first time. Positive ionospheric storm effects were more prevalent during CME-driven and CIR-driven over all stations considered in this study. Negative ionospheric storm effects occurred only during CME-driven storms over mid-latitude stations and were more prevalent in summer. The other interesting finding is that for the stations considered over mid-, low, and equatorial latitudes, negative-positive ionospheric responses were only observed over low and equatorial latitudes. A significant number of cases where the electron density changes remained within the background variability during storm conditions were observed over the low latitude stations compared to other latitude regions.
- Full Text:
- Date Issued: 2018
Machine learning methods for calibrating radio interferometric data
- Authors: Zitha, Simphiwe Nhlanhla
- Date: 2019
- Subjects: Calibration , Radio astronomy -- Data processing , Radio astronomy -- South Africa , Karoo Array Telescope (South Africa) , Radio telescopes -- South Africa , Common Astronomy Software Application (Computer software)
- Language: English
- Type: text , Thesis , Masters , MSc
- Identifier: http://hdl.handle.net/10962/97096 , vital:31398
- Description: The applications of machine learning have created an opportunity to deal with complex problems currently encountered in radio astronomy data processing. Calibration is one of the most important data processing steps required to produce high dynamic range images. This process involves the determination of calibration parameters, both instrumental and astronomical, to correct the collected data. Typically, astronomers use a package such as Common Astronomy Software Applications (CASA) to compute the gain solutions based on regular observations of a known calibrator source. In this work we present applications of machine learning to first generation calibration (1GC), using the KAT-7 telescope environmental and pointing sensor data recorded during observations. Applying machine learning to 1GC, as opposed to calculating the gain solutions in CASA, has shown evidence of reducing computation, as well as accurately predict the 1GC gain solutions representing the behaviour of the antenna during an observation. These methods are computationally less expensive, however they have not fully learned to generalise in predicting accurate 1GC solutions by looking at environmental and pointing sensors. We call this multi-output regression model ZCal, which is based on random forest, decision trees, extremely randomized trees and K-nearest neighbor algorithms. The prediction error obtained during the testing of our model on testing data is ≈ 0.01 < rmse < 0.09 for gain amplitude per antenna, and 0.2 rad < rmse <0.5 rad for gain phase. This shows that the instrumental parameters used to train our model more strongly correlate with gain amplitude effects than phase.
- Full Text:
- Date Issued: 2019