We performed a detailed spectral and timing analysis of a Seyfert 1 galaxy Mrk 509 using data from the Neil Gehrels Swift observatory that spanned over $\sim$ 13 years between 2006 and 2019. To study the variability properties from the optical/UV to X-ray emission, we used a total of 275 pointed observations in this work. The average spectrum over the entire duration exhibits a strong soft X-ray excess above the power law continuum. The soft X-ray excess is well described by two thermal components with temperatures of kT $_{\rm BB1}\sim$ 120 eV and kT $_{\rm BB2}\sim$ 460 eV. The warm thermal component is likely due to the presence of an optically thick and warm Comptonizing plasma in the inner accretion disk. The fractional variability amplitude is found to be decreasing with increasing wavelength, i.e., from the soft X-ray to UV/optical emission. However, the hard X-ray (2–8 keV) emission shows very low variability. The strength of the correlation within the UV and the optical bands (0.95–0.99) is found to be stronger than the correlation between the UV/optical and X-ray bands (0.40–0.53). These results clearly suggest that the emitting regions of the X-ray and UV/optical emission are likely distinct or partly interacting. Having removed the slow variations in the light curves, we find that the lag spectrum is well described by the 4/3 rule for the standard Shakura–Sunyaev accretion disk when we omit X-ray lags. All these results suggest that the real disk is complex, and the UV emission is likely reprocessed in the accretion disk to give X-ray and optical emission.

Equilibrium configurations of the internal magnetic field of a pulsar play a key role in modelling astrophysical phenomena from glitches to gravitational wave emission. In this paper, we present a numerical scheme for solving the Grad–Shafranov equation and calculating equilibrium configurations of pulsars, accounting for superconductivity in the core of the neutron star, and for the Hall effect in the crust of the star. Our numerical code uses a finite difference method in which the source term appearing in the Grad–Shafranov equation, which is used to model the magnetic equilibrium is non-linear. We obtain solutions by linearising the source and applying an under-relaxation scheme at each step of computation to improve the solver’s convergence. We have developed our code in both C++ and Python, and our numerical algorithm can further be adapted to solve any non-linear PDEs appearing in other areas of computational astrophysics. We produce mixed toroidal–poloidal field configurations, and extend the portion of parameter space that can be investigated with respect to previous studies. We find that in even in the more extreme cases, the magnetic energy in the toroidal component does not exceed approximately 5% of the total. We also find that if the core of the star is superconducting, the toroidal component is entirely confined to the crust of the star, which has important implications for pulsar glitch models which rely on the presence of a strong toroidal field region in the core of the star, where superfluid vortices pin to superconducting fluxtubes.

We present a collection of double-lobed sources towards a $20\,\mathrm{deg}^2$ area of the Cygnus region at the northern sky, observed at 325 and 610 MHz with the Giant Metrewave Radio Telescope. The 10 $^{\prime\prime}$ resolution achieved at 325 MHz is 5.5 times better than previous studies, while at 610 MHz, these are the first results ever of such a large area, mapped with 6 $^{\prime\prime}$ angular resolution. After a thorough visual inspection of the images at the two bands, we found 43 double-lobed source candidates, proposed as such due to the presence of 2 bright peaks, within a few arcminutes apart, joined by a bridge or a central nucleus. All but two are presented here as a double-lobed candidates for the first time. Thirty nine of the candidates were covered at both bands, and we provide the spectral index information for them. We have searched for positional coincidences between the detected sources/components and other objects from the literature, along the electromagnetic spectrum. Twenty-three candidates possess radio counterpart(s), 12 present infrared counterparts, and 1 showed an overlapping X-ray source. We analysed each candidate considering morphology, counterparts, and spectral indices. Out of the 43 candidates, 37 show characteristics compatible with an extragalactic nature, 2 of probably Galactic origin, 3 remain as dubious cases, though with feature(s) compatible with an extragalactic nature, and the remaining one, evidence of physically unrelated components. The median spectral index of the 40 putative extragalactic sources is $-1.0$ . Their celestial surface density at 610 MHz resulted in $1.9\,\mathrm{per\ deg}^2$ , across a region lying at the Galactic plane.

We present a broadband radio study of the transient jets ejected from the black hole candidate X-ray binary MAXI J1535–571, which underwent a prolonged outburst beginning on 2017 September 2. We monitored MAXI J1535–571 with the Murchison Widefield Array (MWA) at frequencies from 119 to 186 MHz over six epochs from 2017 September 20 to 2017 October 14. The source was quasi-simultaneously observed over the frequency range 0.84–19 GHz by UTMOST (the Upgraded Molonglo Observatory Synthesis Telescope) the Australian Square Kilometre Array Pathfinder (ASKAP), the Australia Telescope Compact Array (ATCA), and the Australian Long Baseline Array (LBA). Using the LBA observations from 2017 September 23, we measured the source size to be $34\pm1$ mas. During the brightest radio flare on 2017 September 21, the source was detected down to 119 MHz by the MWA, and the radio spectrum indicates a turnover between 250 and 500 MHz, which is most likely due to synchrotron self-absorption (SSA). By fitting the radio spectrum with a SSA model and using the LBA size measurement, we determined various physical parameters of the jet knot (identified in ATCA data), including the jet opening angle ( $\phi_{\rm op} = 4.5\pm1.2^{\circ}$ ) and the magnetic field strength ( $B_{\rm s} = 104^{+80}_{-78}$ mG). Our fitted magnetic field strength agrees reasonably well with that inferred from the standard equipartition approach, suggesting the jet knot to be close to equipartition. Our study highlights the capabilities of the Australian suite of radio telescopes to jointly probe radio jets in black hole X-ray binaries via simultaneous observations over a broad frequency range, and with differing angular resolutions. This suite allows us to determine the physical properties of X-ray binary jets. Finally, our study emphasises the potential contributions that can be made by the low-frequency part of the Square Kilometre Array (SKA-Low) in the study of black hole X-ray binaries.

We present the data and initial results from the first pilot survey of the Evolutionary Map of the Universe (EMU), observed at 944 MHz with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The survey covers $270 \,\mathrm{deg}^2$ of an area covered by the Dark Energy Survey, reaching a depth of 25–30 $\mu\mathrm{Jy\ beam}^{-1}$ rms at a spatial resolution of $\sim$ 11–18 arcsec, resulting in a catalogue of $\sim$ 220 000 sources, of which $\sim$ 180 000 are single-component sources. Here we present the catalogue of single-component sources, together with (where available) optical and infrared cross-identifications, classifications, and redshifts. This survey explores a new region of parameter space compared to previous surveys. Specifically, the EMU Pilot Survey has a high density of sources, and also a high sensitivity to low surface brightness emission. These properties result in the detection of types of sources that were rarely seen in or absent from previous surveys. We present some of these new results here.

We produce the first low to mid-frequency radio simulation that incorporates both traditional extragalactic radio sources as well as synchrotron cosmic web emission. The FIlaments & GAlactic RadiO (FIGARO) simulation includes 10 unique 4° × 4° fields, incorporating active galactic nucleii (AGNs), star-forming galaxies (SFGs), and synchrotron cosmic web emission out to a redshift of z = 0.8 and over the frequency range 100–1 400 MHz. To do this, the simulation brings together a recent 1003 Mpc3 magnetohydrodynamic simulation (Vazza et al. 2019, A&A, 627, A5), calibrated to match observed radio relic population statistics, alongside updated ‘T-RECS’ code for simulating extragalactic radio sources (Bonaldi et al. 2019, MNRAS, 482, 2). Uniquely, the AGNs and SFGs are populated and positioned in accordance with the underlying matter density of the cosmological simulation. In this way, the simulation provides an accurate understanding of the apparent morphology, angular scales, and brightness of the cosmic web as well as—crucially—the clustering properties of the cosmic web with respect to the embedded extragalactic radio population. We find that the synchrotron cosmic web does not closely trace the underlying mass distribution of the cosmic web, but is instead dominated by shocked shells of emission surrounding dark matter halos and resembles a large, undetected population of radio relics. We also show that, with accurate kernels, the cosmic web radio emission is clearly detectable by cross-correlation techniques and this signal is separable from the embedded extragalactic radio population. We offer the simulation as a public resource towards the development of techniques for detecting and measuring the synchrotron cosmic web.

Light curves of photospheric radius expansion (PRE) bursts, a subset of type I X-ray bursts, have been used as standard candles to estimate the ‘nominal PRE distances’ for 63% of PRE bursters (bursters), assuming PRE burst emission is spherically symmetric. Model-independent geometric parallaxes of bursters provide a valuable chance to test models of PRE bursts (PRE models) and can be provided in some cases by Gaia astrometry of the donor stars in bursters. We searched for counterparts to 115 known bursters in the Gaia Early Data Release 3 and confirmed 4 bursters with Gaia counterparts that have detected ( $\gt\!3\,\sigma$ , prior to zero-point correction) parallaxes. We describe a generic approach to the Gaia parallax zero point as well as its uncertainty using an ensemble of Gaia quasars individually determined for each target. Assuming the spherically symmetric PRE model is correct, we refined the resultant nominal PRE distances of three bursters (i.e. Cen $\textrm{X}-4$ , Cyg $\textrm{X}-2$ , and $4\textrm{U}\,0919-54$ ) and put constraints on their compositions of the nuclear fuel powering the bursts. Finally, we describe a method for testing the correctness of the spherically symmetric PRE model using parallax measurements and provide preliminary results.

We previously reported a putative detection of a radio galaxy at $z=10.15$ , selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The redshift of this source, GLEAM J0917–0012, was based on three weakly detected molecular emission lines observed with the Atacama Large Millimetre Array (ALMA). In order to confirm this result, we conducted deep spectroscopic follow-up observations with ALMA and the Karl Jansky Very Large Array (VLA). The ALMA observations targeted the same CO lines previously reported in Band 3 (84–115 GHz) and the VLA targeted the CO(4-3) and [CI(1-0)] lines for an independent confirmation in Q-band (41 and 44 GHz). Neither observation detected any emission lines, removing support for our original interpretation. Adding publicly available optical data from the Hyper Suprime-Cam survey, Widefield Infrared Survey Explorer (WISE), and Herschel Space Observatory in the infrared, as well as $<$ 10 GHz polarisation and 162 MHz inter-planetary scintillation observations, we model the physical and observational characteristics of GLEAM J0917–0012 as a function of redshift. Comparing these predictions and observational relations to the data, we are able to constrain its nature and distance. We argue that if GLEAM J0917–0012 is at $z<3,$ then it has an extremely unusual nature, and that the more likely solution is that the source lies above $z=7$ .

The recent increase in well-localised fast radio bursts (FRBs) has facilitated in-depth studies of global FRB host properties, the source circumburst medium, and the potential impacts of these environments on the burst properties. The Australian Square Kilometre Array Pathfinder (ASKAP) has localised 11 FRBs with sub-arcsecond to arcsecond precision, leading to sub-galaxy localisation regions in some cases and those covering much of the host galaxy in others. The method used to astrometrically register the FRB image frame for ASKAP, in order to align it with images taken at other wavelengths, is currently limited by the brightness of continuum sources detected in the short-duration (‘snapshot’) voltage data captured by the Commensal Real-Time ASKAP Fast Transients (CRAFT) software correlator, which are used to correct for any frame offsets due to imperfect calibration solutions and estimate the accuracy of any required correction. In this paper, we use dedicated observations of bright, compact radio sources in the low- and mid-frequency bands observable by ASKAP to investigate the typical astrometric accuracy of the positions obtained using this so-called ‘snapshot’ technique. Having captured these data with both the CRAFT software and ASKAP hardware correlators, we also compare the offset distributions obtained from both data products to estimate a typical offset between the image frames resulting from the differing processing paths, laying the groundwork for future use of the longer duration, higher signal-to-noise ratio (S/N) data recorded by the hardware correlator. We find typical offsets between the two frames of ${\sim}0.6$ and ${\sim}0.3$ arcsec in the low- and mid-band data, respectively, for both RA and Dec. We also find reasonable agreement between our offset distributions and those of the published FRBs. We detect only a weak dependence in positional offset on the relative separation in time and elevation between target and calibrator scans, with the trends being more pronounced in the low-band data and in Dec. Conversely, the offsets show a clear dependence on frequency in the low band, which we compare to the frequency-dependent Dec. offsets found in FRB 200430. In addition, we present a refined methodology for estimating the overall astrometric accuracy of CRAFT FRBs.

We study in this paper bright late-type giants found in the First Byurakan Survey (FBS) data base. Phase dependent light-curves from large sky area variability data bases such as Catalina Sky Survey (CSS) and All-Sky Automated Survey for Supernovae (ASAS-SN), and the early installment of the third Gaia data release (Gaia EDR3) photometric and astrometric data have been used to characterize our sample of 1 100 M-type giants and 130 C-type stars found at high latitudes. Gaia radial velocities (RV) are available for 134 and luminosities for 158 stars out of 1 100. We show the behaviour of our sample stars in a Gaia color–absolute magnitude diagram (CaMD), the Gaia-2MASS-diagram from Lebzelter et al. with some alternative versions. In this way we explore the potential of these diagrams and their combination for the analysis and interpretation of datasets of LPVs. We show the possibility to classify stars into M- and C-types and to identify the mass of the bulk of the sample stars.

Galaxy clusters have been found to host a range of diffuse, non-thermal emission components, generally with steep, power law spectra. In this work we report on the detection and follow-up of radio halos, relics, remnant radio galaxies, and other fossil radio plasmas in Southern Sky galaxy clusters using the Murchison Widefield Array and the Australian Square Kilometre Array Pathfinder. We make use of the frequency coverage between the two radio interferometers—from 88 to $\sim\!900$ MHz—to characterise the integrated spectra of these sources within this frequency range. Highlights from the sample include the detection of a double relic system in Abell 3186, a mini-halo in RXC J0137.2–0912, a candidate halo and relic in Abell 3399, and a complex multi-episodic head-tail radio galaxy in Abell 3164. We compare this selection of sources and candidates to the literature sample, finding sources consistent with established radio power–cluster mass scaling relations. Finally, we use the low-frequency integrated spectral index, $\alpha$ ( $S_v \propto v^\alpha$ ), of the detected sample of cluster remnants and fossil sources to compare with samples of known halos, relics, remnants and fossils to investigate a possible link between their electron populations. We find the distributions of $\alpha$ to be consistent with relic and halo emission generated by seed electrons that originated in fossil or remnant sources. However, the present sample sizes are insufficient to rule out other scenarios.

The Variables and Slow Transients Survey (VAST) on the Australian Square Kilometre Array Pathfinder (ASKAP) is designed to detect highly variable and transient radio sources on timescales from 5 s to $\sim\!5$ yr. In this paper, we present the survey description, observation strategy and initial results from the VAST Phase I Pilot Survey. This pilot survey consists of $\sim\!162$ h of observations conducted at a central frequency of 888 MHz between 2019 August and 2020 August, with a typical rms sensitivity of $0.24\ \mathrm{mJy\ beam}^{-1}$ and angular resolution of $12-20$ arcseconds. There are 113 fields, each of which was observed for 12 min integration time, with between 5 and 13 repeats, with cadences between 1 day and 8 months. The total area of the pilot survey footprint is 5 131 square degrees, covering six distinct regions of the sky. An initial search of two of these regions, totalling 1 646 square degrees, revealed 28 highly variable and/or transient sources. Seven of these are known pulsars, including the millisecond pulsar J2039–5617. Another seven are stars, four of which have no previously reported radio detection (SCR J0533–4257, LEHPM 2-783, UCAC3 89–412162 and 2MASS J22414436–6119311). Of the remaining 14 sources, two are active galactic nuclei, six are associated with galaxies and the other six have no multi-wavelength counterparts and are yet to be identified.

In this era of spatially resolved observations of planet-forming disks with Atacama Large Millimeter Array (ALMA) and large ground-based telescopes such as the Very Large Telescope (VLT), Keck, and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPace Infrared telescope for Cosmology and Astrophysics (SPICA) is an infrared space mission concept developed jointly by Japan Aerospace Exploration Agency (JAXA) and European Space Agency (ESA) to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage $10{-}220\,\mu\mathrm{m}$ , (2) the high line detection sensitivity of $(1{-}2) \times 10^{-19}\,\mathrm{W\,m}^{-2}$ with $R \sim 2\,000{-}5\,000$ in the far-IR (SAFARI), and $10^{-20}\,\mathrm{W\,m}^{-2}$ with $R \sim 29\,000$ in the mid-IR (SPICA Mid-infrared Instrument (SMI), spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet-forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid-state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. We demonstrate that the SPICA mission concept would allow us to achieve the above ambitious science goals through large surveys of several hundred disks within $\sim\!2.5$ months of observing time.

We explore the effect of anisotropic wind driving on the properties of accretion onto black holes (BHs) in close binaries. We specifically focus on line-driven winds, which are common in high-mass X-ray binaries (HMXBs). In close binary systems, the tidal force from the companion star can modify the wind structure in two different ways. One is the reduction of wind terminal velocity due to the weaker effective surface gravity. The other is the reduction in mass flux due to gravity darkening (GD). We incorporate these effects into the so-called CAK theory in a simple way and investigate the wind flow around the accretor on the orbital scale. We find that a focused accretion stream is naturally formed when the Roche lobe filling factor is ${\gtrsim}0.8$ –0.9, analogous to that of wind Roche lobe overflow, but only when the velocity reduction is taken into account. The formation of a stream is necessary to bring in sufficient angular momentum to form an accretion disc around the BH. GD effects reduce the amount of accreted angular momentum, but not enough to prevent the formation of a disc. Based on these results, we expect there to be a discrete step in the observability of HMXBs depending on whether the donor Roche lobe filling factor is below or above ${\sim}$ 0.8–0.9.

One of the principal systematic constraints on the Epoch of Reionisation (EoR) experiment is the accuracy of the foreground calibration model. Recent results have shown that highly accurate models of extended foreground sources, and including models for sources in both the primary beam and its sidelobes, are necessary for reducing foreground power. To improve the accuracy of the source models for the EoR fields observed by the Murchison Widefield Array (MWA), we conducted the MWA Long Baseline Epoch of Reionisation Survey (LoBES). This survey consists of multi-frequency observations of the main MWA EoR fields and their eight neighbouring fields using the MWA Phase II extended array. We present the results of the first half of this survey centred on the MWA EoR0 observing field (centred at RA (J2000) $0^\mathrm{h}$ , Dec (J2000) $-27^{\circ}$ ). This half of the survey covers an area of 3 069 degrees $^2$ , with an average rms of 2.1 mJy beam–1. The resulting catalogue contains a total of 80 824 sources, with 16 separate spectral measurements between 100 and 230 MHz, and spectral modelling for 78 $\%$ of these sources. Over this region we estimate that the catalogue is 90 $\%$ complete at 32 mJy, and 70 $\%$ complete at 10.5 mJy. The overall normalised source counts are found to be in good agreement with previous low-frequency surveys at similar sensitivities. Testing the performance of the new source models we measure lower residual rms values for peeled sources, particularly for extended sources, in a set of MWA Phase I data. The 2-dimensional power spectrum of these data residuals also show improvement on small angular scales—consistent with the better angular resolution of the LoBES catalogue. It is clear that the LoBES sky models improve upon the current sky model used by the Australian MWA EoR group for the EoR0 field.

The Rapid ASKAP Continuum Survey (RACS) is the first large sky survey using the Australian Square Kilometre Array Pathfinder (ASKAP), covering the sky south of $+41^\circ$ declination. With ASKAP’s large, instantaneous field of view, ${\sim}31\,\mathrm{deg}^2$ , RACS observed the entire sky at a central frequency of 887.5 MHz using 903 individual pointings with 15 minute observations. This has resulted in the deepest radio survey of the full Southern sky to date at these frequencies. In this paper, we present the first Stokes I catalogue derived from the RACS survey. This catalogue was assembled from 799 tiles that could be convolved to a common resolution of $25^{\prime\prime}$ , covering a large contiguous region in the declination range $\delta=-80^{\circ}$ to $+30^\circ$ . The catalogue provides an important tool for both the preparation of future ASKAP surveys and for scientific research. It consists of $\sim$ 2.1 million sources and excludes the $|b|<5^{\circ}$ region around the Galactic plane. This provides a first extragalactic catalogue with ASKAP covering the majority of the sky ( $\delta<+30^{\circ}$ ). We describe the methods to obtain this catalogue from the initial RACS observations and discuss the verification of the data, to highlight its quality. Using simulations, we find this catalogue detects 95% of point sources at an integrated flux density of $\sim$ 5 mJy. Assuming a typical sky source distribution model, this suggests an overall 95% point source completeness at an integrated flux density $\sim$ 3 mJy. The catalogue will be available through the CSIRO ASKAP Science Data Archive (CASDA).

We re-evaluate the outer edge of orbital stability for possible exomoons orbiting the radial velocity planet discovered in the HD 23079 system. In this system, a solar-type star hosts a Jupiter-mass planet in a nearly circular orbit in the outer stellar habitable zone. The outer stability limit of exomoons is deduced using N-body and tidal migration simulations considering a large range of initial conditions, encompassing both prograde and retrograde orbits. In particular, we extend previous works by evaluating many values in the satellite mean anomaly to identify and exclude regions of quasi-stability. Future observations of this system can make use of our results through a scale factor relative to the currently measured minimum mass. Using a constant time lag tidal model (Hut 1981), we find that plausible tidal interactions within the system are insufficient to induce significant outward migration toward the theoretical stability limit. While current technologies are incapable of detecting exomoons in this system, we comment on the detectability of putative moons through Doppler monitoring within direct imaging observations in view of future research capacities.

Next-generation spectro-polarimetric broadband surveys will probe cosmic magnetic fields in unprecedented detail, using the magneto-optical effect known as Faraday rotation. However, non-parametric methods such as RMCLEAN can introduce non-observable linearly polarised flux into a fitted model at negative wavelengths squared. This leads to Faraday rotation structures that are consistent with the observed data, but would be impossible or difficult to measure. We construct a convex non-parametric QU-fitting algorithm to constrain the flux at negative wavelengths squared to be zero. This allows the algorithm to recover structures that are limited in complexity to the observable region in wavelength squared. We verify this approach on simulated broadband data sets where we show that it has a lower root mean square error and that it can change the scientific conclusions for real observations. We advise using this prior in next-generation broadband surveys that aim to uncover complex Faraday depth structures. We provide a public Python implementation of the algorithm at https://github.com/Luke-Pratley/Faraday-Dreams.

Curiously, our Universe was born in a low entropy state, with abundant free energy to power stars and life. The form that this free energy takes is usually thought to be gravitational: the Universe is almost perfectly smooth, and so can produce sources of energy as matter collapses under gravity. It has recently been argued that a more important source of low-entropy energy is nuclear: the Universe expands too fast to remain in nuclear statistical equilibrium, effectively shutting off nucleosynthesis in the first few minutes, providing leftover hydrogen as fuel for stars. Here, we fill in the astrophysical details of this scenario and seek the conditions under which a Universe will emerge from early nucleosynthesis as almost-purely iron. In so doing, we identify a hitherto-overlooked character in the story of the origin of the second law: matter–antimatter asymmetry.

The cosmic evolution of the chemical elements from the Big Bang to the present time is driven by nuclear fusion reactions inside stars and stellar explosions. A cycle of matter recurrently re-processes metal-enriched stellar ejecta into the next generation of stars. The study of cosmic nucleosynthesis and this matter cycle requires the understanding of the physics of nuclear reactions, of the conditions at which the nuclear reactions are activated inside the stars and stellar explosions, of the stellar ejection mechanisms through winds and explosions, and of the transport of the ejecta towards the next cycle, from hot plasma to cold, star-forming gas. Due to the long timescales of stellar evolution, and because of the infrequent occurrence of stellar explosions, observational studies are challenging, as they have biases in time and space as well as different sensitivities related to the various astronomical methods. Here, we describe in detail the astrophysical and nuclear-physical processes involved in creating two radioactive isotopes useful in such studies, $^{26}\mathrm{Al}$ and $^{60}\mathrm{Fe}$ . Due to their radioactive lifetime of the order of a million years, these isotopes are suitable to characterise simultaneously the processes of nuclear fusion reactions and of interstellar transport. We describe and discuss the nuclear reactions involved in the production and destruction of $^{26}\mathrm{Al}$ and $^{60}\mathrm{Fe}$ , the key characteristics of the stellar sites of their nucleosynthesis and their interstellar journey after ejection from the nucleosynthesis sites. This allows us to connect the theoretical astrophysical aspects to the variety of astronomical messengers presented here, from stardust and cosmic-ray composition measurements, through observation of $\gamma$ rays produced by radioactivity, to material deposited in deep-sea ocean crusts and to the inferred composition of the first solids that have formed in the Solar System. We show that considering measurements of the isotopic ratio of $^{26}\mathrm{Al}$ to $^{60}\mathrm{Fe}$ eliminate some of the unknowns when interpreting astronomical results, and discuss the lessons learned from these two isotopes on cosmic chemical evolution. This review paper has emerged from an ISSI-BJ Team project in 2017–2019, bringing together nuclear physicists, astronomers, and astrophysicists in this inter-disciplinary discussion.