New advancements in radio data post-processing are underway within the Square Kilometre Array (SKA) precursor community, aiming to facilitate the extraction of scientific results from survey images through a semi-automated approach. Several of these developments leverage deep learning methodologies for diverse tasks, including source detection, object or morphology classification, and anomaly detection. Despite substantial progress, the full potential of these methods often remains untapped due to challenges associated with training large supervised models, particularly in the presence of small and class-unbalanced labelled datasets.

Self-supervised learning has recently established itself as a powerful methodology to deal with some of the aforementioned challenges, by directly learning a lower-dimensional representation from large samples of unlabelled data. The resulting model and data representation can then be used for data inspection and various downstream tasks if a small subset of labelled data is available.

In this work, we explored contrastive learning methods to learn suitable radio data representations by training the SimCLR model on large collections of unlabelled radio images taken from the ASKAP EMU and SARAO MeerKAT GPS surveys. The resulting models were fine-tuned over smaller labelled datasets, including annotated images from various radio surveys, and evaluated on radio source detection and classification tasks. Additionally, we employed the trained self-supervised models to extract features from radio images, which were used in an unsupervised search for objects with peculiar morphology in the ASKAP EMU pilot survey data. For all considered downstream tasks, we reported the model performance metrics and discussed the benefits brought by self-supervised pre-training, paving the way for building radio foundational models in the SKA era.

We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3 405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor and pathfinder instruments; however, we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We demonstrate that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the Spectrum Roentgen Gamma spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org.

There are different classes of pulsating stars in the H-R diagram. While many of those classes are undisputed, some remain a mystery such as the objects historically called ‘Maia variables’. Whereas the presence of such a class was suggested seven decades ago, no pulsational driving mechanism is known that could excite short-period oscillations in these late B to early A-type stars. Alternative hypotheses that would render the reports of variability of those stars erroneous have been proposed such as incorrect effective temperatures, binarity or rapid rotation, but no certain conclusions have been reached yet. Therefore, the existence of these variables as a homogeneous class of pulsating star is still under discussion. Meanwhile, many new candidates of these variables have been claimed especially by using photometric observations of space telescopes. In this study, we examined 31 objects that are alleged members of this hypothetical group and carried out detailed spectroscopic and photometric analyses to test the proposed hypotheses for their cause of variability. The $T_\textrm{eff}$, $\log g$, $v \sin i$, and chemical abundances of the targets were determined and the TESS photometric data were examined. As a result, we found that most of these targets are located inside the $\delta$ Scuti, $\beta$ Cephei, or SPB star instability strips, a few show evidence for binarity and others for rapid rotation. We give arguments that none of the apparently rapid pulsations in our targets is caused by a star outside any known instability strip. By extrapolation, we argue that most stars proposed as pulsators outside well-established instability domains are misclassified. Hence there is no sufficient evidence justifying the existence of a class of pulsating stars formerly known as the ‘Maia variables’.

New time series photometry of WISE J152614.95-111326.4, an eclipsing binary candidate, has been obtained. Full cycles of variation were covered in five filters, ranging from B to z. Archival time series photometry is also available from several sources. The phased light curve shape changes from a double wave form in the red, to a single wave at shorter wavelengths. Analysis of the spectral energy distribution and SALT spectra shows the presence of a cool ($\sim$7 250–7 900 K) white dwarf and an M6 star. The light curves can be explained by a hot spot on the opposing hemisphere of the white dwarf. The star may be in a pre-cataclysmic variable phase with a very low rate of mass flow from the red dwarf to the white dwarf, such that no flickering is evident. Evidence in favour of this hypothesis is that the period of the system (2.25 h) is in the cataclysmic variable period gap. It is speculated that a weak magnetic field associated with the white dwarf funnels accreted material onto a magnetic pole. Amplitudes of the W1 and W2 WISE light curves are anomalously large. The possibility is discussed that variability in this spectral region is primarily driven by electron cyclotron radiation.

The Vela pulsar (J0835$-$4510) is known to exhibit variations in Faraday rotation and dispersion on multi-decade timescales due to the changing sightline through the surrounding Vela supernova remnant and the Gum Nebula. Until now, variations in Faraday rotation towards Vela have not been studied on timescales less than around a decade. We present the results of a high-cadence observing campaign carried out with the Aperture Array Verification System 2 (AAVS2), a prototype SKA-Low station, which received a significant bandwidth upgrade in 2022. We collected observations of the Vela pulsar and PSR J0630$-$2834 (a nearby pulsar located outside the Gum Nebula), spanning $\sim$1 and $\sim$0.3 yr, respectively, and searched for linear trends in the rotation measure (RM) as a function of time. We do not detect any significant trends on this timescale ($\sim$months) for either pulsar, but the constraints could be greatly improved with more accurate ionospheric models. For the Vela pulsar, the combination of our data and historical data from the published literature have enabled us to model long-term correlated trends in RM and dispersion measure (DM) over the past two decades. We detect a change in DM of $\sim$0.3 $\mathrm{cm}^{-3}\,\mathrm{pc}$ which corresponds to a change in electron density of $\sim$$10^5\,\mathrm{cm}^{-3}$ on a transverse length scale of $\sim$1–2 au. The apparent magnetic field strength in the time-varying region changes from $240^{+30}_{-20}\,\mu\mathrm{G}$ to $-6.2^{+0.7}_{-0.9}\,\mu\mathrm{G}$ over the time span of the dataset. As well as providing an important validation of polarimetry, this work highlights the pulsar monitoring capabilities of SKA-Low stations, and the niche science opportunities they offer for high-precision polarimetry and probing the microstructure of the magneto-ionic interstellar medium.

We present a model-independent way to characterise properties of the magnetic-field turbulence in the emitting regions of Gamma-Ray Burst afterglows. Our only assumption is that afterglows’ synchrotron radiation is efficient. It turns out that the gyroradius of plasma particles must be smaller (with a good margin) than the correlation length of the magnetic-field fluctuations. Such turbulence is essentially non-linear and therefore must be produced by some kind of magnetohydrodynamical instability, likely acting on top of kinetic Weibel instability. We also find that the emitting particles are loosely confined to local magnetic-field structures and diffusion allows them to sample the entire distribution of local magnetisation values. This means that one-zone approach to modelling the afterglow spectra is still valid despite the non-linear nature of the magnetic turbulence. However, the non-linear turbulence may (and likely will) change the synchrotron spectrum of individual electrons.

Drift scan observations provide the broad sky coverage and instrumental stability needed to measure the Epoch of Reionization (EoR) 21-cm signal. In such observations, the telescope’s pointing centre (PC) moves continuously on the sky. The Tracking Tapered Gridded Estimator (TTGE) combines observations from different PC to estimate $P(k_{\perp}, k_{\parallel})$ the 21-cm power spectrum, centred on a tracking centre (TC) which remains fixed on the sky. The tapering further restricts the sky response to a small angular region around TC, thereby mitigating wide-field foregrounds. Here we consider $154.2\,\mathrm{MHz}$ ($z = 8.2$) Murchison Widefield Array (MWA) drift scan observations. The periodic pattern of flagged channels, present in MWA data, is known to introduce artefacts which pose a challenge for estimating $P(k_{\perp}, k_{\parallel})$. Here we have validated the TTGE using simulated MWA drift scan observations which incorporate the flagged channels same as the data. We demonstrate that the TTGE is able to recover $P(k_{\perp}, k_{\parallel})$ without any artefacts and estimate $P(k)$ within $5 \%$ accuracy over a large $k$-range. We also present preliminary results for a single PC, combining 9 nights of observation $(17 \, \mathrm{min}$ total). We find that $P(k_{\perp}, k_{\parallel})$ exhibits streaks at a fixed interval of $k_{\parallel}=0.29 \, \mathrm{Mpc}^{-1}$, which matches $\Delta \nu_\mathrm{per}=1.28 \, \mathrm{MHz}$ that is the period of the flagged channels. Since the simulations demonstrate that the TTGE is impervious to the flagged channels, the streaks seen for the actual data are possibly caused by some systematic that has the same period as the flagged channels. These streaks are more than 3–4 orders of magnitude smaller than the peak foreground power $\mid P(k_{\perp}, k_{\parallel}) \mid \approx 10^{16} \, \mathrm{mK^2}\, \mathrm{Mpc^3}$ at $k_{\parallel}=0$. The streaks are not as pronounced at larger $k_{\parallel}$, and in some cases they do not appear to extend across the entire $k_{\perp}$ range. The rectangular region $0.05 \leq k_{\perp} \leq 0.16 \, \mathrm{Mpc^{-1}}$ and $0.9 \leq k_{\parallel}\leq 4.6 \, \mathrm{Mpc^{-1}}$ is found to be relatively free of foreground contamination and artefacts, and we have used this to place the $2\unicode{x03C3}$ upper limit $\Delta^2(k) < (1.85\times10^4)^2\, \mathrm{mK^2}$ on the EoR 21-cm mean squared brightness temperature fluctuations at $k=1 \,\mathrm{Mpc}^{-1}$.

As TeV gamma-ray astronomy progresses into the era of the Cherenkov Telescope Array (CTA), instantaneously following up on gamma-ray transients is becoming more important than ever. To this end, a worldwide network of Imaging Atmospheric Cherenkov Telescopes has been proposed. Australia is ideally suited to provide coverage of part of the Southern Hemisphere sky inaccessible to H.E.S.S. in Namibia and the upcoming CTA-South in Chile. This study assesses the sources detectable by a small, transient-focused array in Australia based on CTA telescope designs. The TeV emission of extragalactic sources (including the majority of gamma-ray transients) can suffer significant absorption by the extragalactic background light. As such, we explored the improvements possible by implementing stereoscopic and topological triggers, as well as lowered image cleaning thresholds, to access lower energies. We modelled flaring gamma-ray sources based on past measurements from the satellite-based gamma-ray telescope Fermi-LAT. We estimate that an array of four Medium-Sized Telescopes (MSTs) would detect $\sim$24 active galactic nucleus flares >5$\sigma$ per year, up to a redshift of $z\approx1.5$. Two MSTs achieved $\sim$80–90% of the detections of four MSTs. The modelled Galactic transients were detectable within the observation time of one night, 11 of the 21 modelled gamma-ray bursts were detectable, as were $\sim$10% of unidentified transients. An array of MST-class telescopes would thus be a valuable complementary telescope array for transient TeV gamma-ray astronomy.

Environment has long been known to have an impact on the evolution of galaxies, but disentangling its impact from mass evolution requires the careful analysis of statistically significant samples. By implementing cutting-edge visualisation methods to test and validate group-finding algorithms, we utilise a mass-complete sample of galaxies to $z \lt 0.1$ comprised of spectroscopic redshifts from prominent surveys such as the 2-degree Field Galaxy Redshift Survey and the Galaxy and Mass Assembly Survey. Utilising our group finding methods, we find 1 413 galaxy groups made up of 8 990 galaxies corresponding to 36% of galaxies associated with group environments. We also search for close pairs, with separations of $r_\mathrm{sep} \lt 50$ $\text{h}^{-1}\text{kpc}$ and $v_\mathrm{sep} \lt 500 \: \text{km s}^{-1}$ within our sample and further classified them into major ($M_{sec}/M_{prim} \leq$ 0.25) and minor ($M_{sec}/M_{prim} \gt $ 0.25) pairs. To examine the impact of environmental factors, we employ bespoke WISE photometry, which facilitates accurate measurements of stellar mass and star formation rates and hence the best possible description of the variation of galaxy properties as a function of the local environment. Our analysis, employing a derived star-forming main sequence relation, reveals that star-formation (SF) within galaxies are pre-processed as a function of group membership. This is evident from the evolution of the star-forming and quenched population of galaxies. We see an increase in the fraction of quiescent galaxies relative to the field as group membership increases, and this excess of quenched galaxies relative to the field is later quantified through the use of the environmental quenching efficiency ($\varepsilon_{env}$) metric. Within the star-forming population, we observe SF pre-processing with the relative difference in specific star formation rates ($\Delta sSFR$), where we see a net decrease in SF as group membership increases, particularly at larger stellar masses. We again quantify this change within the SF population with our star formation deficiency ($\varepsilon_{SFD}$) metric. Our sample of close pairs at low stellar masses exhibit enhanced star formation efficiencies compared to the field, and at larger stellar mass ranges show large deficiencies. Separating the close pairs into major/minors and primary/secondaries reveals SF enhancements projected separation decreases within the minor pairs, this effect is even more pronounced within minor primaries. This research emphasises the importance of carefully studying the properties of galaxies within group environments to better understand the pre-processing of SF within galaxies. Our results show that the small-scale environments of galaxies influence star-forming properties even when stellar masses are kept constant. This demonstrates that galaxies do not evolve in isolation over cosmic time but are shaped by a complex interaction between their internal dynamics and external influences.

In order to study the structure and temperature distribution within high-mass star-forming clumps, we employed the Australia Telescope Compact Array to image the $\mathrm{NH}_3$ (J,K) = (1,1) through (6,6) and the (2,1) inversion transitions, the $\mathrm{H}_2\mathrm{O}$ $6_{16}$-$5_{23}$ maser line at 22.23508 GHz, several $\mathrm{CH}_3\mathrm{OH}$ lines and hydrogen and helium recombination lines. In addition, 22- and 24-GHz radio continuum emission was also imaged.

The $\mathrm{NH}_3$ lines probe the optical depth and gas temperature of compact structures within the clumps. The $\mathrm{H}_2\mathrm{O}$ maser pinpoints the location of shocked gas associated with star formation. The recombination lines and the continuum emission trace the ionised gas associated with hot OB stars. The paper describes the data and presents sample images and spectra towards select clumps. The technique for estimating gas temperature from $\mathrm{NH}_3$ line ratios is described. The data show widespread hyperfine intensity anomalies in the $\mathrm{NH}_3$ (1,1) images, an indicator of non-LTE $\mathrm{NH}_3$ excitation. We also identify several new $\mathrm{NH}_3$ (3,3) masers associated with shocked gas. Towards AGAL328.809+00.632, the $\mathrm{H}_2\mathrm{O}$ $6_{16}$-$5_{23}$ line, normally seen as a maser, is instead seen as a thermally excited absorption feature against a strong background continuum. The data products are described in detail.

Simulations suggest that slow rotating galaxies are the result of galaxy-galaxy mergers that have a tendency to randomise stellar orbits. The exact pathway for slow rotator formation, however, is still unclear. Our aim is to see whether there is a relationship between fossil groups - whose central galaxies are thought to have undergone more major merging than other central galaxies – and the stellar kinematic properties of those central galaxies. We classify all galaxy groups in the GAMA redshift survey whose central galaxies were observed with SAMI as: (i) fossil groups, (ii) mass gap groups (fossil-like groups), and (iii) groups that are not dynamically evolved (NDEGs, i.e. controls). We compare the following properties of centrals across the three different group types: spin ($\lambda_{Re}$), the fraction of slow rotators ($f_{SR}$), and age. We also repeat our analysis on data from the EAGLE and Magneticum hydrodynamical cosmological simulations. In SAMI, we find that the spin parameter, slow rotator fraction, and age are broadly consistent across our three group types, i.e. the fossil groups, mass gap groups and NDEGs. We do find a weak indication that $f_{SR}$ is slightly lower for fossil group centrals as compared to NDEG centrals. In contrast, in EAGLE and Magneticum, fossil and mass gap group centrals typically have a significantly lower $\lambda_{Re}$ than NDEG centrals. Our results for SAMI suggest that the types of mergers that form fossil groups are not the types of mergers that form slow rotators. Merger count may be less important for slow rotator formation than specific merger conditions, such as the gas content of progenitors. When and where the merging occurs are also suspected to play an important role in slow rotator formation, and these conditions may differ for fossil group formation.

We conducted a study on the X-ray polarisation properties of MCG-5-23-16 by analysing long-term monitoring data from NuSTAR jointly with IXPE observations made in May and November 2022. The re-analysis of IXPE data gives model-dependent polarisation degree, PD (%) = $1.08\pm0.66$ in the energy band 2–8 keV, which agrees with previous studies within error bars. The model-independent analysis of PD poses an upper limit of $\leq3.8$ ($1\sigma$ level) for the same energy band. The observed upper limit of PD, along with broadband spectral analysis (2–79 keV) using an accretion-ejection based model, allowed us to derive the corona geometry (i.e. radius and height) and the accretion disc inclination ($\sim33^\circ$). Additional NuSTAR observations were also analysed to gain insights into the accretion flow properties of the source and to estimate the expected polarisation during those epochs with PD $\sim 4.3\%$. The radius and height of the corona varies between $28.2\pm3.1 - 39.8\pm4.6$ r$_s$ and $14.3\pm1.7-21.4\pm1.9$ r$_s$ respectively, with a mass outflow rate from the corona measuring $0.14\pm0.03-0.2\pm0.03$ Eddington rate ($\dot m_{\mathrm{Edd}}$). The estimated PD values were nearly constant up to a certain radial distance and height of the corona and then decreased for increasing corona geometry. The spectral analysis further provided an estimate for the mass of the central black hole $\sim2\times 10^7$ M$_\odot$ and the velocity of the outflowing gas $\sim0.16-0.19c$. A comparative broadband spectral study using reflection-based models estimates the disc inclination between $\sim 31^\circ\pm8^\circ-45^\circ\pm7^\circ$, and yields an expected PD of 3.4–6.0%. We also found a weak reflection fraction and a less ionised distant reflecting medium. The expected PD measured using accretion-ejection and reflection models is less compared to the expected PD measured for a given disc inclination of $45^\circ$. Our modelling of the disc-corona-outflows and polarisation connection can be extended and validated with data from the recently launched XPoSat, India’s first X-ray Polarimeter Satellite, offering potential applications to other sources.

The angular correlation is a method for measuring the distribution of structure in the Universe, through the statistical properties of the angular distribution of galaxies on the sky. We measure the angular correlation of galaxies from the second data release of the GaLactic and Extragalactic All-sky Murchison Widefield Array eXtended survey (GLEAM-X) survey, a low-frequency radio survey covering declinations below $+30^\circ$. We find an angular distribution consistent with the $\Lambda$CDM cosmological model assuming the best fitting cosmological parameters from Planck Collaboration et al. (2020, A&A, 641, A6). We fit a bias function to the discrete tracers of the underlying matter distribution, finding a bias that evolves with redshift in either a linear or exponential fashion to be a better fit to the data than a constant bias. We perform a covariance analysis to obtain an estimation of the properties of the errors, by analytic, jackknife, and sample variance means. Our results are consistent with previous studies on the topic, and also the predictions of the $\Lambda$CDM cosmological model.

We present deep near-infrared $K_\textrm{s}$-band imaging for 35 of the 53 sources from the high-redshift ($z \gt 2$) radio galaxy candidate sample defined in Broderick et al. (2022, PASA, 39, e061). These images were obtained using the High-Acuity Widefield K-band Imager (HAWK-I) on the Very Large Telescope. Host galaxies are detected for 27 of the sources, with $K_\textrm{s} \approx 21.6$–23.0 mag (2$^{\prime\prime}$ diameter apertures; AB). The remaining eight targets are not detected to a median $3\unicode{x03C3}$ depth of $K_\textrm{s} \approx 23.3$ mag (2$^{\prime\prime}$ diameter apertures). We examine the radio and near-infrared flux densities of the 35 sources, comparing them to the known $z \gt 3$ powerful radio galaxies with 500-MHz radio luminosities $L_{500\,\textrm{MHz}} \gt 10^{27}$ W Hz$^{-1}$. By plotting 150-MHz flux density versus $K_\textrm{s}$-band flux density, we find that, similar to the sources from the literature, these new targets have large radio to near-infrared flux density ratios, but extending the distribution to fainter flux densities. Five of the eight HAWK-I deep non-detections have a median $3\unicode{x03C3}$ lower limit of $K_\textrm{s} \gtrsim 23.8$ mag (1$.\!^{\prime\prime}$5 diameter apertures); these five targets, along with a further source from Broderick et al. (2022, PASA, 39, e061) with a deep non-detection ($K_\textrm{s} \gtrsim 23.7$ mag; $3\unicode{x03C3}$; 2$^{\prime\prime}$ diameter aperture) in the Southern H-ATLAS Regions $K_\textrm{s}$-band Survey, are considered candidates to be ultra-high-redshift ($z \gt 5$) radio galaxies. The extreme radio to near-infrared flux density ratios ($\gt 10^5$) for these six sources are comparable to TN J0924$-$2201, GLEAM J0856$+$0223 and TGSS J1530$+$1049, the three known powerful radio galaxies at $z \gt 5$. For a selection of galaxy templates with different stellar masses, we show that $z \gtrsim 4.2$ is a plausible scenario for our ultra-high-redshift candidates if the stellar mass $M_\textrm{*} \gtrsim 10^{10.5}$ M$_\odot$. In general, the 35 targets studied have properties consistent with the previously known class of infrared-faint radio sources. We also discuss the prospects for finding more UHzRG candidates from wide and deep near-infrared surveys.

The radio interferometric closure phases can be a valuable tool for studying cosmological HI from the early Universe. Closure phases have the advantage of being immune to element-based gains and associated calibration errors. Thus, calibration and errors therein, which are often sources of systematics limiting standard visibility-based approaches, can be avoided altogether in closure phase analysis. In this work, we present the first results of the closure phase power spectrum of HI 21-cm fluctuations using the Murchison Widefield Array (MWA), with $\sim12$ h of MWA phase II observations centred around redshift, $z\approx 6.79$, during the Epoch of Reionisation. On analysing three redundant classes of baselines – 14, 24, and 28 m equilateral triads, our estimates of the $2\sigma$ (95% confidence interval) 21-cm power spectra are $\lesssim(184)^2 pseudo\,\mathrm{mK}^2$ at ${k}_{||} = 0.36 pseudo\ h \mathrm{Mpc}^{-1}$ in the EoR1 field for the 14 m baseline triads, and $\lesssim(188)^2 pseudo\,\mathrm{mK}^2$ at $k_{||} = 0.18 \,pseudo\ h \mathrm{Mpc}^{-1}$ in the EoR0 field for the 24 m baseline triads. The ‘pseudo’ units denote that the length scale and brightness temperature should be interpreted as close approximations. Our best estimates are still 3-4 orders high compared to the fiducial 21-cm power spectrum; however, our approach provides promising estimates of the power spectra even with a small amount of data. These data-limited estimates can be further improved if more datasets are included into the analysis. The evidence for excess noise has a possible origin in baseline-dependent systematics in the MWA data that will require careful baseline-based strategies to mitigate, even in standard visibility-based approaches.

We present high-resolution observations of nearby ($z\lesssim0.1$) galaxies that have hosted Type Ia supernovae to measure systemic spectroscopic redshifts using the wide field spectrograph (WiFeS) instrument on the Australian National University 2.3 m telescope at Siding Spring Observatory. While most of the galaxies targeted have previous spectroscopic redshifts, we provide demonstrably more accurate and precise redshifts with competitive uncertainties, motivated by potential systematic errors that could bias estimates of the Hubble constant ($H_0$). The WiFeS instrument is remarkably stable; after calibration, the wavelength solution varies by $\lesssim$0.5 Å in red and blue with no evidence of a trend over the course of several years. By virtue of the $25\times 38$ arcsec field of view, we are always able to measure the redshift of the galactic core, or the entire galaxy in the cases where its angular extent is smaller than the field of view, reducing any errors due to galaxy rotation. We observed 185 southern SN Ia host galaxies and measured the redshift of each via at least one spatial region of (a) the core and (b) the average over the full-field/entire galaxy. Overall, we find stochastic differences between historical redshifts and our measured redshifts on the order of $\lesssim10^{-3}$ with a mean offset of 4.3${\times 10^{-5}}$ and normalised median absolute deviation of 1.2${\times 10^{-4}}$. We show that a systematic redshift offset at this level is not enough to bias cosmology, as $H_0$ shifts by $+0.1$ km s$^{-1}$ Mpc$^{-1}$ when we replace Pantheon+ redshifts with our own, but the occasional large differences are interesting to note.

The proto-Milky Way epoch forms the earliest stars in our galaxy and sets the initial conditions for the subsequent disk formation. Recent observations showed that the [$\alpha$/Fe] ratio among in situ metal-poor stars declined between [Fe/H] $=-3$ and $-1.3$ until it reached the lowest value ($\sim$0.25) and rose to the traditional value associated with the high-[$\alpha$/Fe] thick disk ($\sim$0.3) at [Fe/H] $\approx$ -1.0. It was suggested that the rise in [$\alpha$/Fe] could be caused by an increase in the star formation efficiency (SFE), known as the ‘simmering’ phase scenario. However, gas inflow also plays a vital role in shaping the star formation history and chemical evolution of galaxies, especially during the earliest epoch of the universe. We investigate this unexpected [$\alpha$/Fe]-rise with an experiment involving a galactic chemical evolution model. Our model has five free parameters: the mass of the initial reservoir of the cold interstellar medium (ISM) at birth, the frequency of Type Ia supernovae (SNe Ia), the cooling timescale of the warm ISM, the SFE, and the inflow rate of fresh gas. The last two free parameters were allowed to change after [$\alpha$/Fe] reached its lowest value, dividing the proto-Galaxy epoch into two phases. The models that reproduced the observed [Fe/H]-[$\alpha$/Fe]-track provided estimates for these fundamental parameters of the proto-Milky Way. We find that the rise in [$\alpha$/Fe] could also be caused by a large inflow of high-[$\alpha$/Fe] gas and conclude that the [$\alpha$/Fe]-rise could be a signature of the gas accretion that fuelled the formation of the Milky Way disk.

We present the fourth data release (DR4) of the SkyMapper Southern Survey (SMSS), the last major step in our hemispheric survey with six optical filters: u, v, g, r, i, z. SMSS DR4 covers 26 000 deg$^{2}$ from over 400 000 images acquired by the 1.3 m SkyMapper telescope between 2014-03 and 2021-09. The 6-band sky coverage extends from the South Celestial Pole to $\delta=+16^{\circ}$, with some images reaching $\delta\sim +28^{\circ}$. In contrast to previous DRs, we include all good-quality images from the facility taken during that time span, not only those explicitly taken for the public Survey. From the image dataset, we produce a catalogue of over 15 billion detections made from $\sim$700 million unique astrophysical objects. The typical 10$\sigma$ depths for each field range between 18.5 and 20.5 mag, depending on the filter, but certain sky regions include longer exposures that reach as deep as 22 mag in some filters. As with previous SMSS catalogues, we have cross-matched with a host of other imaging and spectroscopic datasets to facilitate additional science outcomes. SMSS DR4 is now available to the worldwide astronomical community.

PSR J0837$-$2454 is a young 629 ms radio pulsar whose uncertain distance has important implications. A large distance would place the pulsar far out of the Galactic plane and suggest it is the result of a runaway star, while a short distance would mean the pulsar is extraordinarily cold. Here we present further radio observations and the first deep X-ray observation of PSR J0837$-$2454. Data from the Parkes Murriyang telescope show flux variations over short and long timescales and also yield an updated timing model, while the position and proper motion (and, less strongly, parallax) of the pulsar are constrained by a number of low-significance detections with the Very Long Baseline Array. XMM-Newton data enable detection of X-ray pulsations for the first time from this pulsar and yield a spectrum that is thermal and blackbody-like, with a cool blackbody temperature $\approx$$70\ \mbox{eV}$ or atmosphere temperature $\approx$$50\ \mbox{eV}$, as well as a small hotspot. The spectrum also indicates the pulsar is at a small distance of $\lesssim$$1\ \mbox{kpc}$, which is compatible with the marginal VLBA parallax constraint that favours a distance of $\gtrsim$330 pc. The low implied luminosity ($\sim7.6\times10^{31}\mbox{erg\, s}^{-1}$ at 0.9 kpc) suggests PSR J0837$-$2454 has a mass high enough that fast neutrino emission from direct Urca reactions operates in this young star and points to a nuclear equation of state that allows for direct Urca reactions at the highest densities present in neutron star cores.

The study of very short-period contact binaries provides an important laboratory in which the most important and problematic astrophysical processes of stellar evolution take place. Short-period contact systems, such as CC Com, are particularly important for binary evolution. Close binary systems, especially those with multiple system members, have significant period variations, angular momentum loss mechanisms predominance, and pre-merger stellar evolution, making them valuable astrophysical laboratories. In this study, observations of CC Com, previously reported as a binary system, and new observations from the TÜBİTAK National Observatory (TUG) and the space-based telescope TESS have revealed that there is a third object with a period of about eight years and a fourth object with a period of about a century orbiting the binary system. From simultaneous analysis of all available light curves and radial velocities, the sensitive orbital and physical parameters of the system components are derived. The orbital parameters of the components are P$_\mathrm{A}=0.2206868 \pm 0.0000002$ days, P$_\mathrm{B}=7.9\pm0.1$ yr, P$_\mathrm{C}=98\pm5$ yr, $e_3$ = 0.06, $e_4$ = 0.44 and the physical parameters as M$_\mathrm{A1}=0.712\pm0.009$ M$_{\odot}$, M$_\mathrm{ A2}=0.372\pm0.005$ M$_{\odot}$, $m_{B;i^{\prime}=90^\circ}=0.074$ M$_{\odot}$, $m_{C;i^{\prime}=90^\circ}=0.18$ M$_{\odot}$, R$_\mathrm{A1}=0.693\pm0.006$ R$_{\odot}$, R$_\mathrm{A2}=0.514\pm0.005$ R$_{\odot}$, L$_\mathrm{A1}$ = 0.103 L$_\odot$, L$_\mathrm{A2}$ = 0.081 L$_\odot$. Finally, the evolutionary status of the multiple system CC Com and its component stars is discussed.