[abridged] Accurately accounting for the AGN phase in galaxy evolution requires a large, clean AGN sample. This is now possible with SRG/eROSITA. The public Data Release 1 (DR1, Jan 31, 2024) includes 930,203 sources from the Western Galactic Hemisphere. The data enable the selection of a large AGN sample and the discovery of rare sources. However, scientific return depends on accurate characterisation of the X-ray emitters, requiring high-quality multiwavelength data. This paper presents the identification and classification of optical and infrared counterparts to eRASS1 sources using Gaia DR3, CatWISE2020, and Legacy Survey DR10 (LS10) with the Bayesian NWAY algorithm and trained priors. Sources were classified as Galactic or extragalactic via a Machine Learning model combining optical/IR and X-ray properties, trained on a reference sample. For extragalactic LS10 sources, photometric redshifts were computed using Circlez. Within the LS10 footprint, all 656,614 eROSITA/DR1 sources have at least one possible optical counterpart; about 570,000 are extragalactic and likely AGN. Half are new detections compared to AllWISE, Gaia, and Quaia AGN catalogues. Gaia and CatWISE2020 counterparts are less reliable, due to the surveys shallowness and the limited amount of features available to assess the probability of being an X-ray emitter. In the Galactic Plane, where the overdensity of stellar sources also increases the chance of associations, using conservative reliability cuts, we identify approximately 18,000 Gaia and 55,000 CatWISE2020 extragalactic sources. We release three high-quality counterpart catalogues, plus the training and validation sets, as a benchmark for the field. These datasets have many applications, but in particular empower researchers to build AGN samples tailored for completeness and purity, accelerating the hunt for the Universe most energetic engines.

We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $\mu$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance.

Blazars are often observed to flare across multiple wavelengths. Orphan flares from blazars have been only detected a few times, providing an opportunity to understand the structure of the jet in the accreting system. We report a remarkable orphan X-ray flare from a blazar candidate EP240709a, detected by Einstein Probe (EP) in July 2024. The multi-band spectral properties and variability support EP240709a as a high-energy peaked BL Lacertae-type object. The flux in 0.5-10 keV increases by at least 28 times to the value of low state in 2020, with non-detection of remarkable flaring in other bands during the same period. EP240709a exhibits the harder-when-brighter tendency in the X-ray band during the orphan flare, while its infrared-optical spectra are featureless. We employ one-zone and two-zone leptonic synchrotron self-Compton models to perform the spectral energy distribution fitting. Detecting this rare orphan flare shows the potential of EP in discovering peculiar activities from AGN in high-cadence X-ray sky surveys.

Long-period radio transients (LPTs) are an emerging group of radio transients that show periodic polarised radio bursts with periods varying from a few minutes to a few hours. Fewer than a dozen LPTs have been detected so far, and their origin (source and emission mechanism) remains unclear. Here, we report the discovery of a 1.5 hr LPT, ASKAP J144834-685644, adding to the current sample of sources. ASKAP J144834-685644 is one of the very few LPTs that has been detected from X-rays to radio. It shows a steep radio spectrum and polarised radio bursts, which resemble the radio emission in known LPTs. In addition, it also shows highly structured and periodic narrow-band radio emission. Multi-wavelength properties suggest that the spectral energy distribution (SED) peaks at near ultraviolet wavelengths, indicating the presence of a hot magnetic source. Combining multi-wavelength information, we infer that ASKAP J144834-685644 may be a near edge-on magnetic white dwarf binary (MWD), although we can not fully rule out ASKAP J144834-685644 being an isolated white dwarf pulsar or even a transitional millisecond pulsar (despite the lack of radio pulsations). If ASKAP J144834-685644 is an MWD binary, the observed broadband spectral energy distribution can be explained by emission from an accretion disk. This hints that some fraction of optically bright LPTs may be accreting binaries with the radio period being the orbital period. It might further suggest a connection between optically bright non-accreting synchronized LPTs, such as polars, and non-accreting asynchronous WD pulsars, such as AR Sco and J1912-4410.

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources.

We present high-resolution, phase-resolved spectroscopic observations of the polar EF Eri, obtained with SALT and the SAAO 1.9-m telescope during its recent emergence from a three-decade-long low state. The average spectrum shows strong emission from the Balmer lines (H$\alpha$ and H$\beta$) and He~\textsc{ii} 4686 Å, along with weaker emission from the He~\textsc{i} lines and the Bowen fluorescence (C~\textsc{iii}/N~\textsc{iii}) blend at 4650 Å. The emission lines redward of 5500 Å transition to pure absorption at orbital phases $\sim$0.75--0.95, which we attribute to obscuration of the line-emitting region by the accretion stream. Trailed spectra of the emission lines reveal multicomponent structures consistent with other polars. In this first Doppler study of EF Eri, tomograms of the strongest lines (He~\textsc{ii} 4686 Å and the Balmer lines), using both the standard and inside-out projections, identify three key emission regions: the irradiated face of the secondary star, the ballistic and threading regions of the accretion stream, and the magnetically confined flow. Our Doppler maps show not only the ballistic stream but also two unambiguous magnetic accretion flows, which is consistent with the presence of multiple magnetic accretion regions.

After its discovery in 2016, the white dwarf binary AR Scorpii (AR Sco) remained for several years the only white dwarf system to show pulsed radio emission associated with a fast-spinning white dwarf. The evolutionary origin and the emission mechanism for AR Sco are not completely understood, with different models proposed. Testing and improving these models requires observational input. Here we report the results of a targeted search for other binary white dwarf pulsars like AR Sco. Using data from Gaia and WISE, we identified 56 candidate systems with similar properties to AR Sco, of which 26 were previously uncharacterised. These were subject to spectroscopic and photometric follow-up observations. Aside from one new binary white dwarf pulsar found, J191213.72-441045.1, which was reported in a separate work, we find no other systems whose characteristics are akin to AR Sco. The newly characterised systems are primarily young stellar objects (with 10 found) or cataclysmic variables (7 identifications), with the remaining being either blended or non-variable on short timescales.

The critical lack of structured terminological data for South Africa's official languages hampers progress in multilingual NLP, despite the existence of numerous government and academic terminology lists. These valuable assets remain fragmented and locked in non-machine-readable formats, rendering them unusable for computational research and development. \emph{Marito} addresses this challenge by systematically aggregating, cleaning, and standardising these scattered resources into open, interoperable datasets. We introduce the foundational \emph{Marito} dataset, released under the equitable, Africa-centered NOODL framework. To demonstrate its immediate utility, we integrate the terminology into a Retrieval-Augmented Generation (RAG) pipeline. Experiments show substantial improvements in the accuracy and domain-specific consistency of English-to-Tshivenda machine translation for large language models. \emph{Marito} provides a scalable foundation for developing robust and equitable NLP technologies, ensuring South Africa's rich linguistic diversity is represented in the digital age.

We present a comprehensive photometric, spectroscopic, and polarimetric study of the intermediate polar (IP) 1RXS J080114.6-462324, using observations from the South African Astronomical Observatory (SAAO) 1.0-m and 1.9-m telescopes and the Southern African Large Telescope (SALT), complemented by archival TESS photometry. Photometric and photo-polarimetric data reveal a coherent modulation at the white dwarf (WD) spin period. TESS confirms periodicities of 1307.517 s (spin) and 11.803 h (orbital). Photopolarimetric and circular spectropolarimetric measurements show circular polarisation reaching ~+5%, modulated with the WD spin period, consistent with cyclotron emission from an accreting magnetic pole. Time-resolved optical spectra display prominent Balmer (H$\gamma$, H$\beta$, and H$\alpha$) and He\textsc{ii} $\lambda$4686 emission features and additional He\textsc{i} emissions, all exhibiting minimal radial-velocity variations. We detect red-shifted absorption dips adjacent to the He\textsc{ii} and H$\beta$ lines, modulated at the WD spin period; periodogram analysis of the emission lines also yields spin modulation. These observations indicate that the system is a disc-fed, low-inclination IP. The combination of circular polarisation and spin-modulated absorption by infalling accretion curtain material supports this classification. Its comparatively long orbital period among IPs and the detection of polarised emission render 1RXS J080114.6-462324 an appealing candidate for evolutionary studies, potentially offering insight into how magnetic accretion systems evolve toward synchronism at longer orbital periods.

A planetary microlensing signal is generally characterized by a short-term perturbation to the standard single lensing light curve. A subset of binary-source events can produce perturbations that mimic planetary signals, thereby introducing an ambiguity between the planetary and binary-source interpretations. In this paper, we present analysis of the microlensing event MOA-2012-BLG-486, for which the light curve exhibits a short-lived perturbation. Routine modeling not considering data taken in different passbands yields a best-fit planetary model that is slightly preferred over the best-fit binary-source model. However, when allowed for a change in the color during the perturbation, we find that the binary-source model yields a significantly better fit and thus the degeneracy is clearly resolved. This event not only signifies the importance of considering various interpretations of short-term anomalies, but also demonstrates the importance of multi-band data for checking the possibility of false-positive planetary signals.

We present the first NuSTAR X-ray observation of EF Eri, a well-known polar system. The NuSTAR observation was conducted in conjunction with NICER, shortly after EF Eri entered a high accretion state following an unprecedented period of low activity lasting 26 years since 1997. NuSTAR detected hard X-ray emission up to 50 keV with an X-ray flux of $1.2\times10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ ($3\rm{-}50$ keV). Folded X-ray lightcurves exhibit a single peak with $\sim65\%$ spin modulation throughout the $3\rm{-}50$ keV band. We found no evidence of QPO signals at $\nu = 0.1\rm{-}100$ Hz with an upper limit on the QPO amplitude below $5\%$ ($90\%$ CL) at $\nu \sim 0.5$ Hz where the optical QPO was previously detected. Our 1-D accretion column model, called $\texttt{MCVSPEC}$, was fitted to the NuSTAR spectral data, yielding an accurate WD mass measurement of $M = (0.55\rm{-}0.63) M_\odot$. ${\tt MCVSPEC}$ accounts for radiative cooling by thermal bremsstrahlung and cyclotron emission, X-ray reflection off the WD surface, and a previously constrained range of the accretion column area. The derived WD mass range is in excellent agreement with the previous measurement of $M = (0.55\rm{-}0.65) M_\odot$ in the optical band. This demonstrates a combination of broadband X-ray spectral analysis and the ${\tt MCVSPEC}$ model that can be employed in our ongoing NuSTAR observation campaign of other polars to determine their WD masses accurately.

NGC 3105 is a young open cluster hosting blue, yellow and red supergiants. This rare combination makes it an excellent laboratory to constrain evolutionary models of high-mass stars. It is poorly studied and fundamental parameters such as its age or distance are not well defined. We intend to characterize in an accurate way the cluster as well as its evolved stars, for which we derive for the first time atmospheric parameters and chemical abundances. We identify 126 B-type likely members within a radius of 2.7$\pm$0.6 arcmin, which implies an initial mass, $M_{cl}\approx$4100 M$_{\odot}$. We find a distance of 7.2$\pm$0.7 kpc for NGC 3105, placing it at $R_{GC}$=10.0$\pm$1.2 kpc. Isochrone fitting supports an age of 28$\pm$6 Ma, implying masses around 9.5 M$_{\odot}$ for the supergiants. A high fraction of Be stars ($\approx$25 %) is found at the top of the main sequence down to spectral type b3. From the spectral analysis we estimate for the cluster a $v_{rad}$=+46.9$\pm$0.9 km s$^{-1}$ and a low metallicity, [Fe/H]=-0.29$\pm$0.22. We also have determined, for the first time, chemical abundances for Li, O, Na, Mg, Si, Ca, Ti, Ni, Rb, Y, and Ba for the evolved stars. The chemical composition of the cluster is consistent with that of the Galactic thin disc. An overabundance of Ba is found, supporting the enhanced $s$-process. NGC 3105 has a low metallicity for its Galactocentric distance, comparable to typical LMC stars. It is a valuable spiral tracer in a very distant region of the Carina-Sagittarius spiral arm, a poorly known part of the Galaxy. As one of the few Galactic clusters containing blue, yellow and red supergiants, it is massive enough to serve as a testbed for theoretical evolutionary models close to the boundary between intermediate and high-mass stars.

Credit scorecards are models used for the modelling of the probability of default of clients. The decision to extend credit to an applicant, as well as the price of the credit, is often based on these models. In order to ensure that scorecards remain accurate over time, the hypothesis of population stability is tested periodically; that is, the hypothesis that the distributions of the attributes of clients at the time when the scorecard was developed is still representative of these distributions at review is tested. A number of measures of population stability are used in practice, with several being proposed in the recent literature. This paper provides a critical review of several testing procedures for the mentioned hypothesis. The widely used population stability index is discussed alongside two recently proposed techniques. Additionally, the use of classical goodness-of-fit techniques is considered and the problems associated with large samples are investigated. In addition to the existing testing procedures, we propose two new techniques which can be used to test population stability. The first is based on the calculation of effect sizes which does not suffer the same problems as classical goodness-of-fit techniques when faced with large samples. The second proposed procedure is the so-called overlapping statistic. We argue that this simple measure can be useful due to its intuitive interpretation. In order to demonstrate the use of the various measures, as well as to highlight their strengths and weaknesses, several numerical examples are included.

We report the discovery of a peculiar X-ray transient, EP240408a, by Einstein Probe (EP) and follow-up studies made with EP, Swift, NICER, GROND, ATCA and other ground-based multi-wavelength telescopes. The new transient was first detected with Wide-field X-ray Telescope (WXT) on board EP on April 8th, 2024, manifested in an intense yet brief X-ray flare lasting for 12 seconds. The flare reached a peak flux of 3.9x10^(-9) erg/cm2/s in 0.5-4 keV, about 300 times brighter than the underlying X-ray emission detected throughout the observation. Rapid and more precise follow-up observations by EP/FXT, Swift and NICER confirmed the finding of this new transient. Its X-ray spectrum is non-thermal in 0.5-10 keV, with a power-law photon index varying within 1.8-2.5. The X-ray light curve shows a plateau lasting for about 4 days, followed by a steep decay till becoming undetectable about 10 days after the initial detection. Based on its temporal property and constraints from previous EP observations, an unusual timescale in the range of 7-23 days is found for EP240408a, which is intermediate between the commonly found fast and long-term transients. No counterparts have been found in optical and near-infrared, with the earliest observation at 17 hours after the initial X-ray detection, suggestive of intrinsically weak emission in these bands. We demonstrate that the remarkable properties of EP240408a are inconsistent with any of the transient types known so far, by comparison with, in particular, jetted tidal disruption events, gamma-ray bursts, X-ray binaries and fast blue optical transients. The nature of EP240408a thus remains an enigma. We suggest that EP240408a may represent a new type of transients with intermediate timescales of the order of about 10 days. The detection and follow-ups of more of such objects are essential for revealing their origin.

Symbiotic stars are binaries in which a cool and evolved star of luminosity class I-III accretes onto a smaller companion. However, direct accretion signatures like disc flickering and boundary layer X-rays are typically outshone or suppressed by the luminous giant, shell burning on the accreting white dwarf, and the illuminated wind nebula. We present a new way to find symbiotics that is less biased against directly-detectable accretion discs than methods based on narrow-band H$\alpha$ photometry or objective prism plate surveys. We identified outliers in SkyMapper survey photometry, using reconstructed uvg snapshot colours and rapid variability among the three exposures of each 20-minute SkyMapper Main Survey filter sequence, from a sample of 366,721 luminous red objects. We found that SkyMapper catalog colours of large-amplitude pulsating giants must be corrected for variability, and that flickering is detectable with only three data points. Our methods probed a different region of parameter space than a recent search for accreting-only symbiotics in the GALAH survey, while being surprisingly concordant with another survey's infrared detection algorithm. We discovered 12 new symbiotics, including four with optical accretion disc flickering. Two of the optical flickerers exhibited boundary-layer hard X-rays. We also identified 10 symbiotic candidates, and discovered likely optical flickering in the known symbiotic V1044 Cen (CD-36 8436). We conclude that at least 20% of the true population of symbiotics exhibit detectable optical flickering from the inner accretion disc, the majority of which do not meet the H$\alpha$ detection thresholds used to find symbiotics in typical narrow-band surveys.

NGC 1566 is known for exhibiting recurrent outbursts, which are accompanied by changes in spectral type. The most recent transient event occurred from 2017 to 2019 and was reported to be accompanied by a change in Seyfert classification from Seyfert 1.8 to Seyfert 1.2. We analyze data from an optical spectroscopic variability campaign of NGC 1566 taken with the 9.2m SALT between July 2018 and October 2019 and supplement our data set with optical to near-infrared spectroscopic archival data taken by VLT/MUSE in September 2015 and October 2017. We observe the emergence and fading of a strong power-law-like blue continuum as well as strong variations in the Balmer, HeI, HeII lines and the coronal lines [FeVII], [FeX] and [FeXI]. Moreover, we detect broad double-peaked emission line profiles of OI 8446 and the CaII 8498,8542,8662 triplet. This is the first time that genuine double-peaked OI 8446 and CaII 8498,8542,8662 emission in AGN is reported in the literature. All broad lines show a clear redward asymmetry with respect to their central wavelength and we find indications for a significant blueward drift of the total line profiles during the transient event. We show that the double-peaked emission line profiles are well approximated by emission from a low-inclination, relativistic eccentric accretion disk, and that single-peaked profiles can be obtained by broadening due to scale-height dependent turbulence. Small-scale features in the OI and CaII lines suggest the presence of inhomogeneities in the broad-line region. We conclude that the broad-line region in NGC 1566 is dominated by the kinematics of a relativistic eccentric accretion disk. The broad-line region can be modeled to be vertically stratified with respect to scale-height turbulence. The observed blueward drift might be attributed to a low-optical-depth wind launched during the transient event.

We present multi-wavelength temporal and spectral characteristics of a magnetic cataclysmic variable (MCV) Swift J0503.7-2819, using far ultraviolet (FUV) and X-ray data from AstroSat, supplemented with optical data from the Southern African Large Telescope and X-ray data from the XMM-Newton and Swift observatories. The X-ray modulations at 4897.6657 s and 3932.0355 s are interpreted as the orbital ($P_{\Omega}$) and spin ($P_{\omega}$) period, respectively, and are consistent with prior reports. With a spin-orbit period ratio of 0.8 and $P_{\Omega}$ falling below the period gap (2-3 hrs) of CVs, Swift J0503.7-2819 would be the newest addition to the growing population of nearly synchronous MCVs, which we call EX Hya-like systems. Hard X-ray luminosity of < $2.5\times10^{32} erg s^{-1}$, as measured with the Swift Burst Alert Telescope, identifies it to be a low-luminosity intermediate polar, similar to other EX Hya-like systems. The phenomenology of the light curves and the spectral characteristics rule out a purely disc-fed/stream-fed model and instead reveal the presence of complex accretion structures around the white dwarf. We propose a ring-like accretion flow, akin to EX Hya, using period ratio, stability arguments, and observational features. An attempt is made to differentiate between the asynchronous polar/nearly-synchronous intermediate polar nature of Swift J0503.7-2819. Further, we note that with the advent of sensitive surveys, a growing population of MCVs that exhibit characteristics of both polars and intermediate polars is beginning to be identified, likely forming a genealogical link between the two conventional classes of MCVs.

White dwarf stars are the most common stellar fossils. When in binaries, they make up the dominant form of compact object binary within the Galaxy and can offer insight into different aspects of binary formation and evolution. One of the most remarkable white dwarf binary systems identified to date is AR Scorpii (henceforth AR Sco). AR Sco is composed of an M-dwarf star and a rapidly-spinning white dwarf in a 3.56-hour orbit. It shows pulsed emission with a period of 1.97 minutes over a broad range of wavelengths, which led to it being known as a white dwarf pulsar. Both the pulse mechanism and the evolutionary origin of AR Sco provide challenges to theoretical models. Here we report the discovery of the first sibling of AR Sco, J191213.72-441045.1 (henceforth J1912-4410), which harbours a white dwarf in a 4.03-hour orbit with an M-dwarf and exhibits pulsed emission with a period of 5.30 minutes. This discovery establishes binary white dwarf pulsars as a class and provides support for proposed formation models for white dwarf pulsars.

In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$\mu$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$\mu$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data.

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $\gamma$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $\gamma$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $\gamma$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $\gamma$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $\gamma$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $\gamma$-ray cosmology.