After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass loading factor for every galaxy. Building on the model’s resolved angular momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus–age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z = 0$ stellar and H i mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and python, is publicly available at https://github.com/arhstevens/DarkSage.

This study presents an analysis of the optical variability of the blazar 1E 1458.8+2249 on diverse timescales using multi-band observations, including observations in the optical BVRI bands carried out with the T60 and T100 telescopes from 2020 to 2023 and ZTF gri data from 2018 to 2023. On seven nights, we searched for intraday variability using the power-enhanced F-test and the nested ANOVA test, but no significant variability was found. The long-term light curve shows a variability behaviour in the optical BVRI bands with amplitudes of $\sim$100% and in the gri bands with amplitudes of $\sim$120%, including short-term variability of up to $\sim$1.1 mag. Correlation analysis revealed a strong correlation between the optical multi-band emissions without any time lag. From 62 nightly spectral energy distributions, we obtained spectral indices between 0.826 and 1.360, with an average of $1.128\pm0.063$. The relationships of both spectral indices and colour with respect to brightness indicate a mild BWB trend throughout the observation period, both intraday and long-term. We also performed a periodicity search using the weighted wavelet Z-transform and Lomb–Scargle methods. A recurrent optical emission pattern with a quasi-periodicity of $\sim$340 days is detected in the combined V- and R-band light curves. The observational results indicate that the blazar 1E 1458.8+2249 has a complex variability, while emphasising the need for future observations to unravel its underlying mechanisms.

We present radio observations of the galaxy cluster Abell S1136 at 888 MHz, using the Australian Square Kilometre Array Pathfinder radio telescope, as part of the Evolutionary Map of the Universe Early Science program. We compare these findings with data from the Murchison Widefield Array, XMM-Newton, the Wide-field Infrared Survey Explorer, the Digitised Sky Survey, and the Australia Telescope Compact Array. Our analysis shows the X-ray and radio emission in Abell S1136 are closely aligned and centered on the Brightest Cluster Galaxy, while the X-ray temperature profile shows a relaxed cluster with no evidence of a cool core. We find that the diffuse radio emission in the centre of the cluster shows more structure than seen in previous low-resolution observations of this source, which appeared formerly as an amorphous radio blob, similar in appearance to a radio halo; our observations show the diffuse emission in the Abell S1136 galaxy cluster contains three narrow filamentary structures visible at 888 MHz, between $\sim$80 and 140 kpc in length; however, the properties of the diffuse emission do not fully match that of a radio (mini-)halo or (fossil) tailed radio source.

The baryon mass content (i.e. stellar and gas mass) of dark matter halos in the early Universe depends on both global factors – for example, ionising ultraviolet (UV) radiation background – and local factors – for example, star formation efficiency and assembly history. We use a lightweight semi-analytical model to investigate how both local and global factors impact the halo baryon mass content at redshifts of $z\geq 5$. Our model incorporates a time delay between when stars form and when they produce feedback of $0\leq t^d/\mathrm{Myr} \leq 30$, which can drive bursts of star formation, and a mass and redshift-dependent UV background, which captures the influence of cosmological reionisation on gas accretion onto halos. We use statistically representative halo assembly histories and assume that the cosmological gas accretion rate is proportional to the halo mass accretion rate. Delayed ($t^d$>0) feedback leads to oscillations in gas mass with cosmic time, behaviour that cannot be captured with instantaneous feedback ($t^d$=0). Highly efficient star formation drives stronger oscillations, while strong feedback impacts when oscillations occur; in contrast, inefficient star formation and weak feedback produce similar long-term behaviour to that observed in instantaneous feedback models. If the delayed feedback timescale is too long, a halo retains its gas reservoir but the feedback suppresses star formation. Our model predicts that lower mass systems (halo masses $m_\mathrm{h} \leq 10^7 \mathrm{M}_\odot$) at $z \leq 10$ should be strongly gas deficient ($m_\mathrm{g}\rightarrow 0$), whereas higher mass systems retain their gas reservoirs because they are sufficiently massive to continue accreting gas through cosmological reionisation. Interestingly, in higher mass halos, the median $m_\star/(m_\star+m_\mathrm{g}) \simeq 0.01-0.05$, but is a factor of 3–5 smaller when feedback is delayed. Our model does not include seed supermassive black hole feedback, which is necessary to explain massive quenched galaxies in the early Universe.

In gamma-ray binaries neutron star is orbiting a companion that produces a strong stellar wind. We demonstrate that observed properties of ‘stellar wind’–‘pulsar wind’ interaction depend both on the overall wind thrust ratio, as well as more subtle geometrical factors: the relative direction of the pulsar’s spin, the plane of the orbit, the direction of motion, and the instantaneous line of sight. Using fully 3D relativistic magnetohydrodynamical simulations we find that the resulting intrinsic morphologies can be significantly orbital phase-dependent: a given system may change from tailward-open to tailward-closed shapes. As a result, the region of unshocked pulsar wind can change by an order of magnitude over a quarter of the orbit. We calculate radiation maps and synthetic light curves for synchrotron (X-ray) and inverse-Compton emission (GeV-TeV), taking into account $\gamma $–$\gamma $ absorption. Our modelled light curves are in agreement with the phase-dependent observed light curves of LS5039.

Recent discoveries of multiple long-period pulsars (periods ${\sim}10\,$s or larger) are starting to challenge the conventional notion that coherent radio emission cannot be produced by objects that are below the many theorised death lines. Many of the past pulsar surveys and software have been prone to selection effects that restricted their sensitivities towards long-period and sporadically emitting objects. Pulsar surveys using new-generation low-frequency facilities are starting to employ longer dwell times, which makes them significantly more sensitive in detecting long-period or nulling pulsars. There have also been software advancements to aid more sensitive searches towards long-period objects. Furthermore, recent discoveries suggest that nulling may be a key aspect of the long-period pulsar population. We simulate both long-period and nulling pulsar signals, using the Southern-sky MWA Rapid Two-meter (SMART) survey data as reference and explore the detection efficacy of popular search methods such as the fast Fourier transform (FFT), fast-folding algorithm (FFA) and single pulse search (SPS). For FFT-based search and SPS, we make use of the PRESTO implementation, and for FFA we use RIPTIDE. We find RIPTIDE’s FFA to be more sensitive; however, it is also the slowest algorithm. PRESTO’s FFT, although faster than others, also shows some unexpected inaccuracies in detection properties. SPS is highly sensitive to long-period and nulling signals, but only for pulses with high intrinsic signal-to-noise ratios. We use these findings to inform current and future pulsar surveys that aim to uncover a large population of long-period or nulling objects and comment on how to make optimal use of these methods in unison.

Traditional pulsar surveys have primarily employed time-domain periodicity searches. However, these methods are susceptible to effects like scattering, eclipses, and orbital motion. At lower radio frequencies ($\lesssim$300 MHz), factors such as dispersion measure and pulse broadening become more prominent, reducing the detection sensitivity. On the other hand, image domain searches for pulsars are not limited by these effects and can extend the parameter space to regions inaccessible to traditional search techniques. Therefore, we have developed a pipeline to form 1-second full Stokes images from offline correlated high time-resolution data from the Murchison Widefield Array (MWA). This led to the development of image-based methodologies to identify new pulsar candidates. In this paper, we applied these methodologies to perform a low-frequency image-based pulsar census of the galactic plane (12 MWA observations, covering $\sim$6 000 $\textrm{deg}^\textrm{2}$ sky). This work focuses on the detection of the known pulsar population which were present in the observed region of the sky using both image-based and beamformed methods. This resulted in the detection of 83 known pulsars, with 16 pulsars found only in Stokes I images but not in periodicity searches applied in beamformed data. Notably, for 14 pulsars these are the first reported low-frequency detections. This underscores the importance of image-based searches for pulsars that may be undetectable in time-series data, due to scattering and/or dispersive smearing at low frequencies. This highlights the importance of low-frequency flux density measurements in refining pulsar spectral models and investigating the spectral turnover of pulsars at low frequencies.

A number of nearby dwarf galaxies have globular cluster (GC) candidates that require spectroscopic confirmation. Here, we present Keck telescope spectra for 15 known GCs and GC candidates that may be associated with a host dwarf galaxy and an additional 3 GCs in the halo of M31 that are candidates for accretion from a now-disrupted dwarf galaxy. We confirm six star clusters (of intermediate-to-old age) to be associated with NGC 247. The vast bulk of its GC system remains to be studied spectroscopically. We also confirm the GC candidates in F8D1 and DDO190, finding both to be young star clusters. The three M31 halo GCs all have radial velocities consistent with M31 and are old and very metal-poor. Their ages and metallicities are consistent with accretion from a low-mass satellite galaxy. Finally, three objects are found to be background galaxies – two are projected near NGC 247 and one (candidate GCC7) is near the IKN dwarf. The IKN dwarf thus has only five confirmed GCs but still a remarkable specific frequency of 124.

In order to study exoplanets, a comprehensive characterisation of the fundamental properties of the host stars – such as angular diameter, temperature, luminosity, and age, is essential, as the formation and evolution of exoplanets are directly influenced by the host stars at various points in time. In this paper, we present interferometric observations taken of directly imaged planet host 51 Eridani at the CHARA Array. We measure the limb-darkened angular diameter of 51 Eridani to be $\theta_\mathrm{LD} = 0.450\pm 0.006$ mas and combining with the Gaia zero-point corrected parallax, we get a stellar radius of $1.45 \pm 0.02$ R$_{\odot}$. We use the PARSEC isochrones to estimate an age of $23.2^{+1.7}_{-2.0}$ Myr and a mass of $1.550^{+0.006}_{-0.005}$ M$_{\odot}$. The age and mass agree well with values in the literature, determined through a variety of methods ranging from dynamical age trace-backs to lithium depletion boundary methods. We derive a mass of $4.1\pm0.4$ M$_\mathrm{Jup}$ for 51 Eri b using the Sonora Bobcat models, which further supports the possibility of 51 Eri b forming under either the hot-start formation model or the warm-start formation model.

Fast Radio Bursts (FRBs) are millisecond dispersed radio pulses of predominately extra-galactic origin. Although originally discovered at GHz frequencies, most FRBs have been detected between 400 and 800 MHz. Nevertheless, only a handful of FRBs were detected at radio frequencies $\le$400 MHz. Searching for FRBs at low frequencies is computationally challenging due to increased dispersive delay that must be accounted for. Nevertheless, the wide field of view (FoV) of low-frequency telescopes – such as the the Murchison Widefield Array (MWA), and prototype stations of the low-frequency Square Kilometre Array (SKA-Low) – makes them promising instruments to open a low-frequency window on FRB event rates, and constrain FRB emission models. The standard approach, inherited from high-frequencies, is to form multiple tied-array beams to tessellate the entire FoV and perform the search on the resulting time series. This approach, however, may not be optimal for low-frequency interferometers due to their large FoVs and high spatial resolutions leading to a large number of beams. Consequently, there are regions of parameter space in terms of number of antennas and resolution elements (pixels) where interferometric imaging is computationally more efficient. Here we present a new high-time resolution imager BLINK implemented on modern graphical processing units (GPUs) and intended for radio astronomy data. The main goal for this imager is to become part of a fully GPU-accelerated FRB search pipeline. We describe the imager and present its verification on real and simulated data processed to form all-sky and widefield images from the MWA and prototype SKA-Low stations. We also present and compare benchmarks of the GPU and CPU code executed on laptops, desktop computers, and Australian supercomputers. The code is publicly available at https://github.com/PaCER-BLINK-Project/imager and can be applied to data from any radio telescope.

We report the results of our analysis of six gravity-mode pulsating hot subdwarf stars observed in the short cadence mode by Transiting Exoplanet Survey Satellite. We detected at least 10 pulsation periods in each star, searched for multiplets, and used an asymptotic period spacing to identify modes. We used a grid of evolutionary and pulsation models calculated with the MESA and GYRE, along with spectroscopic parameters and modal degree identification, to derive the physical properties of the stars. We checked the relation between the helium content and pulsations and found that no pulsator exists among the extremely helium-rich hot subdwarfs, while the number of detected pulsators in other helium groups increases as the helium content decreases. We found p- and g-mode hot subdwarfs pulsators in all Galactic populations.

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.

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.

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.

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 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.

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.

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.

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.