We demonstrate the importance of radio selection in probing heavily obscured galaxy populations. We combine Evolutionary Map of the Universe (EMU) Early Science data in the Galaxy and Mass Assembly (GAMA) G23 field with the GAMA data, providing optical photometry and spectral line measurements, together with Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry, providing IR luminosities and colours. We investigate the degree of obscuration in star-forming galaxies, based on the Balmer decrement (BD), and explore how this trend varies, over a redshift range of $0<z<0.345$. We demonstrate that the radio-detected population has on average higher levels of obscuration than the parent optical sample, arising through missing the lowest BD and lowest mass galaxies, which are also the lower star formation rate (SFR) and metallicity systems. We discuss possible explanations for this result, including speculation around whether it might arise from steeper stellar initial mass functions in low mass, low SFR galaxies.

We investigate the connection between a cluster’s structural configuration and observable measures of its gas emission that can be obtained in X-ray and Sunyaev–Zeldovich (SZ) surveys. We present an analytic model for the intracluster gas density profile: parameterised by the dark matter halo’s inner logarithmic density slope, $\alpha$, the concentration, c, the gas profile’s inner logarithmic density slope, $\varepsilon$, the dilution, d, and the gas fraction, $\eta$, normalised to cosmological content. We predict four probes of the gas emission: the emission-weighted, $T_\mathrm{X}$, and mean gas mass-weighted, $T_\mathrm{m_g}$, temperatures, and the spherically, $Y_\mathrm{sph}$, and cylindrically, $Y_\mathrm{cyl}$, integrated Compton parameters. Over a parameter space of clusters, we constrain the X-ray temperature scaling relations, $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$, within $57.3\%$ and $41.6\%$, and $M_{200} - T_\mathrm{m_g}$ and $M_{500} - T_\mathrm{m_g}$, within $25.7\%$ and $7.0\%$, all respectively. When excising the cluster’s core, the $M_{200} - T_\mathrm{X}$ and $M_{500} - T_\mathrm{X}$ relations are further constrained, to within $31.3\%$ and $17.1\%$, respectively. Similarly, we constrain the SZ scaling relations, $M_{200} - Y_\mathrm{sph}$ and $M_{500} - Y_\mathrm{sph}$, within $31.1\%$ and $17.7\%$, and $M_{200} - Y_\mathrm{cyl}$ and $M_{500} - Y_\mathrm{cyl}$, within $25.2\%$ and $22.0\%$, all respectively. The temperature observable $T_\mathrm{m_g}$ places the strongest constraint on the halo mass, whilst $T_\mathrm{X}$ is more sensitive to the parameter space. The SZ constraints are sensitive to the gas fraction, whilst insensitive to the form of the gas profile itself. In all cases, the halo mass is recovered with an uncertainty that suggests the cluster’s structural profiles only contribute a minor uncertainty in its scaling relations.

This study builds upon our prior work to further explore and unravel the effects of saturated thermal conduction within a viscous resistive MHD framework on the intricate transport mechanisms of angular momentum and energy in disc-jet systems. We conducted a series of 2.5-dimensional non-relativistic time-dependent numerical simulations using the PLUTO code. Employing a saturation parameter spanning [0.002-0.01], our results are consistent with previous investigations that omitted consideration of thermal conduction, affirming the established understanding that kinetic torque plays a predominant role in governing the total accretion angular momentum, surpassing the magnetic contribution within the disc. At the initial time steps of our calculations, we find that thermal conduction enhances this kinetic contribution, while concurrently diminishing the effect of magnetic contribution. In contrast to the prevailing influence of kinetic torque within the disc, we also assert the magnetic torque as the primary contributor to the total ejection angular momentum. We further unveil that doubling the saturation parameter leads to bolstering of approximately $23.7\%$ in the integral dominance of magnetic torque compared to kinetic torque within the jet. Our findings reveal that doubling the effect of thermal conduction improves the integral total accretion power by approximately 2%, thereby slightly amplifying the energy content within the system and increasing overall energy output. We underscore that as the local energy dissipation within the disc intensifies, the significance of the enthalpy accretion flux increases at the expense of the jet power. We reveal that increasing the saturation parameter mitigates enthalpy accumulation within the disc, and further restricts the jet’s energy extraction from the disc. This limitation is determined in our analysis through the decrease in the integral ratio between the bipolar jet and liberated power of approximately $13.8\%$, for twice the strength of the saturation parameter. We identify the Poynting flux as the primary contributor to total jet power, with thermal conduction exerting minimal influence on magnetic contributions. Additionally, we emphasise the integration of jet enthalpy as another significant factor in determining overall jet power, highlighting a distinct correlation between the rise in saturation parameter and heightened enthalpy contribution. Moreover, we observe the promotion of Poynting flux over kinetic flux at advanced time steps of our simulations, a trend supported by the presence of thermal conduction, which demonstrates an integral increase of approximately $11.2\%$ when considering a doubling of the saturation parameter.

Odd Radio Circles (ORCs) are a class of low surface brightness, circular objects approximately one arcminute in diameter. ORCs were recently discovered in the Australian Square Kilometre Array Pathfinder (ASKAP) data and subsequently confirmed with follow-up observations on other instruments, yet their origins remain uncertain. In this paper, we suggest that ORCs could be remnant lobes of powerful radio galaxies, re-energised by the passage of a shock. Using relativistic hydrodynamic simulations with synchrotron emission calculated in post-processing, we show that buoyant evolution of remnant radio lobes is alone too slow to produce the observed ORC morphology. However, the passage of a shock can produce both filled and edge-brightnened ORC-like morphologies for a wide variety of shock and observing orientations. Circular ORCs are predicted to have host galaxies near the geometric centre of the radio emission, consistent with observations of these objects. Significantly offset hosts are possible for elliptical ORCs, potentially causing challenges for accurate host galaxy identification. Observed ORC number counts are broadly consistent with a paradigm in which moderately powerful radio galaxies are their progenitors.

Post-asymptotic giant branch stars (post-AGB) in binary systems, with typical orbital periods between $\sim\!100$ to $\sim$1 000 days, result from a poorly understood interaction that terminates their precursory AGB phase. The majority of these binaries display a photospheric anomaly called ‘chemical depletion’, thought to arise from an interaction between the circumbinary disc and the post-AGB star, leading to the reaccretion of pure gas onto the star, devoid of refractory elements due to dust formation. In this paper, we focus on a subset of chemically peculiar binary post-AGBs in the Galaxy and the Magellanic Clouds (MCs). Our detailed stellar parameter and chemical abundance analysis utilising high-resolution optical spectra from VLT+UVES revealed that our targets span a $T_{\rm eff}$ of 4 900–7 250 K and [Fe/H] of −0.5 - −1.57 dex. Interestingly, these targets exhibit a carbon ([C/Fe] ranging from 0.5 - 1.0 dex, dependant on metallicity) and s-process enrichment ($\textrm{[s/Fe]}\,\geq\!1$dex) contrary to the commonly observed chemical depletion pattern. Using spectral energy distribution (SED) fitting and period–luminosity–colour (PLC) relation methods, we determine the luminosity of the targets (2 700–8 300 $\rm L_{\odot}$), which enables confirmation of their evolutionary phase and estimation of initial masses (as a function of metallicity) (1–2.5 $\textrm{M}_{\odot}$). In conjunction with predictions from dedicated ATON stellar evolutionary models, our results indicate a predominant intrinsic enrichment of carbon and s-process elements in our binary post-AGB targets. We qualitatively rule out extrinsic enrichment and inherited s-process enrichment from the host galaxy as plausible explanations for the observed overabundances. Our chemically peculiar subset of intrinsic carbon and s-process enriched binary post-AGBs also hints at potential variation in the efficiency of chemical depletion between stars with C-rich and O-rich circumbinary disc chemistries. However, critical observational studies of circumbinary disc chemistry, along with specific condensation temperature estimates in C-rich environments, are necessary to address gaps in our current understanding of disc-binary interactions inducing chemical depletion in binary post-AGB systems.

Clusters of galaxies have been found to host Mpc-scale diffuse, non-thermal radio emission in the form of central radio halos and peripheral relics. Turbulence and shock-related processes in the intra-cluster medium are generally considered responsible for the emission, though details of these processes are still not clear. The low surface brightness makes detection of the emission a challenge, but with recent surveys with high-sensitivity radio telescopes we are beginning to build large samples of these sources. The Evolutionary Map of the Universe (EMU) is a Southern Sky survey being performed by the Australian SKA Pathfinder (ASKAP) over the next few years and is well-suited to detect and characterise such emission. To assess prospects of the full survey, we have performed a pilot search of diffuse sources in 71 clusters from the Planck Sunyaev–Zeldovich (SZ) cluster catalogue (PSZ2) found in archival ASKAP observations. After re-imaging the archival data and performing both (u, v)-plane and image-plane angular scale filtering, we detect 21 radio halos (12 for the first time, excluding an additional six candidates), 11 relics (in seven clusters, and six for the first time, excluding a further five candidate relics), along with 12 other, unclassified diffuse radio sources. From these detections, we predict the full EMU survey will uncover up to $\approx 254$ radio halos and $\approx 85$ radio relics in the 858 PSZ2 clusters that will be covered by EMU. The percentage of clusters found to host diffuse emission in this work is similar to the number reported in recent cluster surveys with the LOw Frequency ARray (LOFAR) Two-metre Sky Survey [Botteon, et al. 2022a, A&A, 660, A78], suggesting EMU will complement similar searches being performed in the Northern Sky and provide us with statistically significant samples of halos and relics at the completion of the full survey. This work presents the first step towards large samples of the diffuse radio sources in Southern Sky clusters with ASKAP and eventually the SKA.

We present source detection and catalogue construction pipelines to build the first catalogue of radio galaxies from the 270 $\rm deg^2$ pilot survey of the Evolutionary Map of the Universe (EMU-PS) conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The detection pipeline uses Gal-DINO computer vision networks (Gupta et al. 2024, PASA, 41, e001) to predict the categories of radio morphology and bounding boxes for radio sources, as well as their potential infrared host positions. The Gal-DINO network is trained and evaluated on approximately 5 000 visually inspected radio galaxies and their infrared hosts, encompassing both compact and extended radio morphologies. We find that the Intersection over Union (IoU) for the predicted and ground-truth bounding boxes is larger than 0.5 for 99% of the radio sources, and 98% of predicted host positions are within $3^{\prime \prime}$ of the ground-truth infrared host in the evaluation set. The catalogue construction pipeline uses the predictions of the trained network on the radio and infrared image cutouts based on the catalogue of radio components identified using the Selavy source finder algorithm. Confidence scores of the predictions are then used to prioritise Selavy components with higher scores and incorporate them first into the catalogue. This results in identifications for a total of 211 625 radio sources, with 201 211 classified as compact and unresolved. The remaining 10 414 are categorised as extended radio morphologies, including 582 FR-I, 5 602 FR-II, 1 494 FR-x (uncertain whether FR-I or FR-II), 2 375 R (single-peak resolved) radio galaxies, and 361 with peculiar and other rare morphologies. Each source in the catalogue includes a confidence score. We cross-match the radio sources in the catalogue with the infrared and optical catalogues, finding infrared cross-matches for 73% and photometric redshifts for 36% of the radio galaxies. The EMU-PS catalogue and the detection pipelines presented here will be used towards constructing catalogues for the main EMU survey covering the full southern sky.

Detecting and mitigating radio frequency interference (RFI) is critical for enabling and maximising the scientific output of radio telescopes. The emergence of machine learning (ML) methods capable of handling large datasets has led to their application in radio astronomy, particularly in RFI detection. Spiking neural networks (SNNs), inspired by biological systems, are well suited for processing spatio-temporal data. This study introduces the first exploratory application of SNNs to an astronomical data processing task, specifically RFI detection. We adapt the nearest latent neighbours (NLNs) algorithm and auto-encoder architecture proposed by previous authors to SNN execution by direct ANN2SNN conversion, enabling simplified downstream RFI detection by sampling the naturally varying latent space from the internal spiking neurons. Our subsequent evaluation aims to determine whether SNNs are viable for future RFI detection schemes. We evaluate detection performance with the simulated HERA telescope and hand-labelled LOFAR observation dataset the original authors provided. We additionally evaluate detection performance with a new MeerKAT-inspired simulation dataset that provides a technical challenge for machine-learnt RFI detection methods. This dataset focuses on satellite-based RFI, an increasingly important class of RFI and is an additional contribution. Our SNN approach remains competitive with the original NLN algorithm and AOFlagger in AUROC, AUPRC, and F1-scores for the HERA dataset but exhibits difficulty in the LOFAR and Tabascal datasets. However, our method maintains this accuracy while completely removing the compute and memory-intense latent sampling step found in NLN. This work demonstrates the viability of SNNs as a promising avenue for ML-based RFI detection in radio telescopes by establishing a minimal performance baseline on traditional and nascent satellite-based RFI sources and is the first work to our knowledge to apply SNNs in astronomy.

Generation of science-ready data from processed data products is one of the major challenges in next-generation radio continuum surveys with the Square Kilometre Array (SKA) and its precursors, due to the expected data volume and the need to achieve a high degree of automated processing. Source extraction, characterization, and classification are the major stages involved in this process. In this work we focus on the classification of compact radio sources in the Galactic plane using both radio and infrared images as inputs. To this aim, we produced a curated dataset of $\sim$20 000 images of compact sources of different astronomical classes, obtained from past radio and infrared surveys, and novel radio data from pilot surveys carried out with the Australian SKA Pathfinder. Radio spectral index information was also obtained for a subset of the data. We then trained two different classifiers on the produced dataset. The first model uses gradient-boosted decision trees and is trained on a set of pre-computed features derived from the data, which include radio-infrared colour indices and the radio spectral index. The second model is trained directly on multi-channel images, employing convolutional neural networks. Using a completely supervised procedure, we obtained a high classification accuracy (F1-score > 90%) for separating Galactic objects from the extragalactic background. Individual class discrimination performances, ranging from 60% to 75%, increased by 10% when adding far-infrared and spectral index information, with extragalactic objects, PNe and Hii regions identified with higher accuracies. The implemented tools and trained models were publicly released and made available to the radioastronomical community for future application on new radio data.

NASA’s all-sky survey mission, the Transiting Exoplanet Survey Satellite (TESS), is specifically engineered to detect exoplanets that transit bright stars. Thus far, TESS has successfully identified approximately 400 transiting exoplanets, in addition to roughly 6 000 candidate exoplanets pending confirmation. In this study, we present the results of our ongoing project, the Validation of Transiting Exoplanets using Statistical Tools (VaTEST). Our dedicated effort is focused on the confirmation and characterisation of new exoplanets through the application of statistical validation tools. Through a combination of ground-based telescope data, high-resolution imaging, and the utilisation of the statistical validation tool known as TRICERATOPS, we have successfully discovered eight potential super-Earths. These planets bear the designations: TOI-238b (1.61$^{+0.09} _{-0.10}$ R$_\oplus$), TOI-771b (1.42$^{+0.11} _{-0.09}$ R$_\oplus$), TOI-871b (1.66$^{+0.11} _{-0.11}$ R$_\oplus$), TOI-1467b (1.83$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-1739b (1.69$^{+0.10} _{-0.08}$ R$_\oplus$), TOI-2068b (1.82$^{+0.16} _{-0.15}$ R$_\oplus$), TOI-4559b (1.42$^{+0.13} _{-0.11}$ R$_\oplus$), and TOI-5799b (1.62$^{+0.19} _{-0.13}$ R$_\oplus$). Among all these planets, six of them fall within the region known as ‘keystone planets’, which makes them particularly interesting for study. Based on the location of TOI-771b and TOI-4559b below the radius valley we characterised them as likely super-Earths, though radial velocity mass measurements for these planets will provide more details about their characterisation. It is noteworthy that planets within the size range investigated herein are absent from our own solar system, making their study crucial for gaining insights into the evolutionary stages between Earth and Neptune.

We present results regarding the longitudinal migrations of cool stellar spots that exhibit remarkable oscillations and explore their possible causes. We conducted analyses using high-quality data from nine target systems of various spectral types, spanning from F to M, which were observed by the Kepler Satellite. The systems in which the behaviour of the spots was examined are as follows: KIC 4357272, KIC 6025466, KIC 6058875, KIC 6962018, KIC 7798259, KIC 9210828, KIC 11706658, KIC 12599700, and KIC 8669092. Basic stellar parameters were calculated from light curve analysis using the PHOEBE V.0.32 software, and light curves were modelled to obtain sinusoidal variations occurring out-of-eclipses phases, induced by rotational modulation. Subsequently, we calculated the minimum times of the obtained sinusoidal variations using the Fourier transform. The distributions of ${\theta}_{\min}$ corresponding to these minimum times over time were computed using linear fits to determine the longitudinal migrations of the spotted areas. We then compared the longitudinal migration periods with the stellar parameters found in the literature. In addition, we also found a secondary variation in the spot migrations apart from the linear models. Our results revealed that the longitudinal migration periods vary in relation to the $B-V$ colour index of the stars.

We report the discovery of a bow-shock pulsar wind nebula (PWN), named Potoroo, and the detection of a young pulsar J1638$-$4713 that powers the nebula. We present a radio continuum study of the PWN based on 20-cm observations obtained from the Australian Square Kilometre Array Pathfinder (ASKAP) and MeerKAT. PSR J1638$-$4713 was identified using Parkes radio telescope observations at frequencies above 3 GHz. The pulsar has the second-highest dispersion measure of all known radio pulsars (1 553 pc cm$^{-3}$), a spin period of 65.74 ms and a spin-down luminosity of $\dot{E}=6.1\times10^{36}$ erg s$^{-1}$. The PWN has a cometary morphology and one of the greatest projected lengths among all the observed pulsar radio tails, measuring over 21 pc for an assumed distance of 10 kpc. The remarkably long tail and atypically steep radio spectral index are attributed to the interplay of a supernova reverse shock and the PWN. The originating supernova remnant is not known so far. We estimated the pulsar kick velocity to be in the range of 1 000–2 000 km s$^{-1}$ for ages between 23 and 10 kyr. The X-ray counterpart found in Chandra data, CXOU J163802.6$-$471358, shows the same tail morphology as the radio source but is shorter by a factor of 10. The peak of the X-ray emission is offset from the peak of the radio total intensity (Stokes $\rm I$) emission by approximately 4.7$^{\prime\prime}$, but coincides well with circularly polarised (Stokes $\rm V$) emission. No infrared counterpart was found.

Since the turn of the century, astronomers have been exploiting the rich information afforded by combining stellar kinematic maps and imaging in an attempt to recover the intrinsic, three-dimensional (3D) shape of a galaxy. A common intrinsic shape recovery method relies on an expected monotonic relationship between the intrinsic misalignment of the kinematic and morphological axes and the triaxiality parameter. Recent studies have, however, cast doubt about underlying assumptions relating shape and intrinsic kinematic misalignment. In this work, we aim to recover the 3D shape of individual galaxies using their projected stellar kinematic and flux distributions using a supervised machine learning approach with mixture density network (MDN). Using a mock dataset of the EAGLE hydrodynamical cosmological simulation, we train the MDN model for a carefully selected set of common kinematic and photometric parameters. Compared to previous methods, we demonstrate potential improvements achieved with the MDN model to retrieve the 3D galaxy shape along with the uncertainties, especially for prolate and triaxial systems. We make specific recommendations for recovering galaxy intrinsic shapes relevant for current and future integral field spectroscopic galaxy surveys.

Exploration of the 21cm signal during the Cosmic Dawn and the Epoch of Reionisation (EoR) can unravel the mysteries of the early Universe when the first stars and galaxies were born and ionised, respectively. However, the 21 cm signal is exceptionally weak, and thus, the detection amidst the bright foregrounds is extremely challenging. The Murchison Widefield Array (MWA) aims to measure the brightness temperature fluctuations of neutral hydrogen from the early Universe. The MWA telescope observes the radio sky with a large field of view (FoV) that causes the bright galaxies, especially near the horizon, to contaminate the measurements. These foregrounds contaminating the EoR datasets must be meticulously removed or treated to detect the signal successfully. The Central Redundant Array Mega-tile (CRAM) is a zenith-pointing new instrument, installed at the centre of the MWA Phase II southern hexagonal configuration, comprising of 64 dipoles in an $8 \times 8$ configuration with a FoV half the width of the MWA’s at every frequency under consideration. The primary objective of this new instrument is to mitigate the impact of bright radio sources near the field centre in accordance with the reduced primary beamshape and to reduce the contamination of foreground sources near the horizon with the reduced sidelobe response of the larger array configuration. In this paper, we introduce the new instrument to the community and present the system architecture and characteristics of the instrument. Using the first light observations, we determine the CRAM system temperature and system performance.

In large-scale galaxy surveys, particularly deep ground-based photometric studies, galaxy blending was inevitable. Such blending posed a potential primary systematic uncertainty for upcoming surveys. Current deblenders predominantly depended on analytical modelling of galaxy profiles, facing limitations due to inflexible and imprecise models. We presented a novel approach, using a U-net structured transformer-based network for deblending astronomical images, which we term the CAT-deblender. It was trained using both RGB and the grz-band images, spanning two distinct data formats present in the Dark Energy Camera Legacy Survey (DECaLS) database, including galaxies with diverse morphologies in the training dataset. Our method necessitated only the approximate central coordinates of each target galaxy, sourced from galaxy detection, bypassing assumptions on neighbouring source counts. Post-deblending, our RGB images retained a high signal-to-noise peak, consistently showing superior structural similarity against ground truth. For multi-band images, the ellipticity of central galaxies and median reconstruction error for r-band consistently lie within $\pm$0.025 to $\pm$0.25, revealing minimal pixel residuals. In our comparison of deblending capabilities focused on flux recovery, our model showed a mere 1% error in magnitude recovery for quadruply blended galaxies, significantly outperforming SExtractor’s higher error rate of 4.8%. Furthermore, by cross-matching with the publicly accessible overlapping galaxy catalogs from the DECaLS database, we successfully deblended 433 overlapping galaxies. Moreover, we have demonstrated effective deblending of 63 733 blended galaxy images, randomly chosen from the DECaLS database.

High-precision pulsar timing observations are limited in their accuracy by the jitter noise that appears in the arrival time of pulses. Therefore, it is important to systematically characterise the amplitude of the jitter noise and its variation with frequency. In this paper, we provide jitter measurements from low-frequency wideband observations of PSR J0437$-$4715 using data obtained as part of the Indian Pulsar Timing Array experiment. We were able to detect jitter in both the 300–500 MHz and 1 260–1 460 MHz observations of the upgraded Giant Metrewave Radio Telescope (uGMRT). The former is the first jitter measurement for this pulsar below 700 MHz, and the latter is in good agreement with results from previous studies. In addition, at 300–500 MHz, we investigated the frequency dependence of the jitter by calculating the jitter for each sub-banded arrival time of pulses. We found that the jitter amplitude increases with frequency. This trend is opposite as compared to previous studies, indicating that there is a turnover at intermediate frequencies. It will be possible to investigate this in more detail with uGMRT observations at 550–750 MHz and future high-sensitive wideband observations from next generation telescopes, such as the Square Kilometre Array. We also explored the effect of jitter on the high precision dispersion measure (DM) measurements derived from short duration observations. We find that even though the DM precision will be better at lower frequencies due to the smaller amplitude of jitter noise, it will limit the DM precision for high signal-to-noise observations, which are of short durations. This limitation can be overcome by integrating for a long enough duration optimised for a given pulsar.

In radio astronomy, the science output of a telescope is often limited by computational resources. This is especially true for transient and technosignature surveys that need to search high-resolution data across a large parameter space. The tremendous data volumes produced by modern radio array telescopes exacerbate these processing challenges. Here, we introduce a ‘reduced-resolution’ beamforming approach to alleviate downstream processing requirements. Our approach, based on post-correlation beamforming, allows sensitivity to be traded against the number of beams needed to cover a given survey area. Using the MeerKAT and Murchison Widefield Array telescopes as examples, we show that survey speed can be vastly increased, and downstream signal processing requirements vastly decreased, if a moderate sacrifice to sensitivity is allowed. We show the reduced-resolution beamforming technique is intimately related to standard techniques used in synthesis imaging. We suggest that reduced-resolution beamforming should be considered to ease data processing challenges in current and planned searches; further, reduced-resolution beamforming may provide a path to computationally expensive search strategies previously considered infeasible.

The unprecedented imaging power of James Webb Space Telescope (JWST) provides new abilities to observe the shapes of objects in the early Universe in a way that has not been possible before. Recently, JWST acquired a deep field image inside the same field imaged in the past as the Hubble Space Telescope (HST) Ultra Deep Field. Computer-based quantitative analysis of spiral galaxies in that field shows that among 34 galaxies for which their rotation of direction can be determined by the shapes of the arms, 24 rotate clockwise, and just 10 rotate counterclockwise. The one-tailed binomial distribution probability to have asymmetry equal or stronger than the observed asymmetry by chance is $\sim$0.012. While the analysis is limited by the small size of the data, the observed asymmetry is aligned with all relevant previous large-scale analyses from all premier digital sky surveys, all show a higher number of galaxies rotating clockwise in that part of the sky, and the magnitude of the asymmetry increases as the redshift gets higher. This paper also provides data and analysis to reproduce previous experiments suggesting that the distribution of galaxy rotation in the Universe is random, to show that the exact same data used in these studies in fact show non-random distribution, and in excellent agreement with the results shown here. These findings reinforce consideration of the possibility that the directions of rotation of spiral galaxies as observed from Earth are not necessarily randomly distributed. The explanation can be related to the large-scale structure of the Universe, but can also be related to a possible anomaly in the physics of galaxy rotation.

The nature of the first galaxies that reionised the universe during the Epoch of Reionisation (EoR) remains unclear. Attempts to directly determine spectral properties of these early galaxies are affected by both limited photometric constraints across the spectrum and by the opacity of the intergalactic medium to the Lyman Continuum (LyC) at high redshift. We approach this by analysing properties of analogous extreme emission line galaxies (EELGs, [OIII]+Hbeta EW $\gt 400$) at $2.5\lt z\lt 4$ from the ZFOURGE survey using the Multi-wavelength Analysis of Galaxy Physical Properties (MAGPHYS) SED fitting code. We compare these to galaxies at $z \gt 5.5$ observed with the James Webb Space Telesope with self-consistent spectral energy distribution fitting methodology. This work focuses on the comparison of their UV slopes (${\unicode{x03B2}}_P$), ionising photon production efficiencies $\xi_{ion}$, star formation rates and dust properties to determine the effectiveness of this analogue selection technique. We report the median ionising photon production efficiencies as log$_{10}(\xi_{ion}/(Hz\ {\rm erg}^{-1}))=$$25.14^{+0.06}_{-0.04}$, $25.16^{+0.06}_{-0.05}$, $25.16^{+0.04}_{-0.05}$, $25.18^{+0.06}_{-0.07}$ for our ZFOURGE control, ZFOURGE EELG, JADES, and CEERS samples, respectively. ZFOURGE EELGs are 0.57 dex lower in stellar mass and have half the dust extinction, compared to their ZFOURGE control counterparts. They also have a similar specific star formation rates and ${\unicode{x03B2}}_P$ to the $z\gt 5.5$ samples. We find that EELGs at low redshift ($2.5\lt z\lt 4$) are analogous to EoR galaxies in their dust attenuation and specific star formation rates. Their extensive photometric coverage and the accessibility of their LyC region opens pathways to infer stellar population properties in the EoR.