Bayesian inference is a powerful tool in gravitational-wave astronomy. It enables us to deduce the properties of merging compact-object binaries and to determine how these mergers are distributed as a population according to mass, spin, and redshift. As key results are increasingly derived using Bayesian inference, there is increasing scrutiny on Bayesian methods. In this review, we discuss the phenomenon of model misspecification, in which results obtained with Bayesian inference are misleading because of deficiencies in the assumed model(s). Such deficiencies can impede our inferences of the true parameters describing physical systems. They can also reduce our ability to distinguish the ‘best fitting’ model: it can be misleading to say that Model A is preferred over Model B if both models are manifestly poor descriptions of reality. Broadly speaking, there are two ways in which models fail. Firstly, models that fail to adequately describe the data (either the signal or the noise) have misspecified likelihoods. Secondly, population models—designed, for example, to describe the distribution of black hole masses—may fail to adequately describe the true population due to a misspecified prior. We recommend tests and checks that are useful for spotting misspecified models using examples inspired by gravitational-wave astronomy. We include companion python notebooks to illustrate essential concepts.

A plausible formation scenario for the Galactic globular clusters 47 Tucanae (47 Tuc) and Omega Centauri $(\omega$ Cen) is that they are tidally stripped remnants of dwarf galaxies, in which case they are likely to have retained a fraction of their dark matter cores. In this study, we have used the ultra-wide band receiver on the Parkes telescope (Murriyang) to place upper limits on the annihilation rate of exotic Light Dark Matter particles $(\chi)$ via the $\chi\chi\rightarrow e^+e^-$ channel using measurements of the recombination rate of positronium (Ps). This is an extension of a technique previously used to search for Ps in the Galactic Centre. However, by stacking of spectral data at multiple line frequencies, we have been able to improve sensitivity. Our measurements have resulted in $3-\sigma$ flux density (recombination rate) upper limits of 1.7 mJy $\left(1.4\times 10^{43}\, \mathrm{s}^{-1}\right)$ and 0.8 mJy $\left(1.1 \times 10^{43} \mathrm{s}^{-1}\right)$ for 47 Tuc and $\omega$ Cen, respectively. Within the Parkes beam at the cluster distances, which varies from 10–23 pc depending on the frequency of the recombination line, and for an assumed annihilation cross-section $\langle\sigma v\rangle = 3\times 10^{-29} \mathrm{cm}^3\, \mathrm{s}^{-1}$ , we calculate upper limits to the dark matter mass and rms dark matter density of ${\lesssim} 1.2-1.3\times 10^5 f_n^{-0.5}$ $\left(m_\chi/\mathrm{MeV\, c}^{-2}\right)$ $\mathrm{M}_{\odot}$ and ${\lesssim} 48-54 f_n^{-0.5}$ $\left(m_\chi/\mathrm{MeV\, c}^{-2}\right)$ $\mathrm{M}_{\odot} \mathrm{pc}^{-3}$ for the clusters, where $f_n=R_n/R_p$ is the ratio of Ps recombination transitions to annihilations, estimated to be ${\sim}0.01$ . The radio limits for $\omega$ Cen suggest that, for a fiducial dark/luminous mass ratio of ${\sim}0.05$ , any contribution from Light Dark Matter is small unless $\langle\sigma v\rangle < 7.9\times 10^{-28}\ \left(m_\chi/\mathrm{MeV\, c}^{-2}\right)^2 \mathrm{cm}^3 \mathrm{s}^{-1}$ . Owing to the compactness and proximity of the clusters, archival 511-keV measurements suggest even tighter limits than permitted by CMB anisotropies, $\langle\sigma v\rangle < 8.6\times 10^{-31}\ (m_\chi/\mathrm{MeV\, c}^{-2})^2 \mathrm{cm}^3 \mathrm{s}^{-1}$ . Due to the very low synchrotron radiation background, our recombination rate limits substantially improve on previous radio limits for the Milky Way.

Ultra-compact H ii (UC HII) regions are an important phase in the formation and early evolution of massive stars and a key component of the interstellar medium (ISM). The main objectives of this work are to study the young stellar population associated with the G45.07+0.13 and G45.12+0.13 UC HII regions, as well as the ISM in which they are embedded. We determined the distribution of the hydrogen column density (N( $\mathrm{H}_2$ )) and dust temperature ( $T_d$ ) in the molecular cloud using Modified blackbody fitting on Herschel images obtained in four bands: 160, 250, 350, and $500\,\unicode{x03BC}\mathrm{m}$ . We used near-, mid-, and far-infrared photometric data to identify and classify the young stellar objects (YSOs). Their main parameters were determined by the radiation transfer models. We also constructed a colour-magnitude diagram and K luminosity functions (KLFs) to compare the parameters of stellar objects with the results of the radiative transfer models. We found that N( $\mathrm{H}_2$ ) varies from ${\sim}3.0 \times 10^{23}$ to $5.5 \times 10^{23}\,\mathrm{cm}^{-2}$ within the G45.07+0.13 and G45.12+0.13 regions, respectively. The maximum $T_d$ value is 35 K in G45.12+0.13 and 42 K in G45.07+0.13. $T_d$ then drops significantly from the centre to the periphery, reaching about 18–20 K at distances of ${\sim}2.6$ and ${\sim}3.7\,\mathrm{pc}$ from InfraRed Astronomical Satellite (IRAS) 19110+1045 (G45.07+0.13) and IRAS 19111+1048 (G45.12+0.13), respectively. The gas plus dust mass value included in G45.12+0.13 is ${\sim}3.4 \times 10^5\,\mathrm{M}_\odot$ and ${\sim}1.7 \times 10^5\,\mathrm{M}_\odot$ in G45.07+0.13. The UC HII regions are connected through a cold ( $T_d = 19\,\mathrm{K}$ ) bridge. The radial surface density distribution of the identified 518 YSOs exhibits dense clusters in the vicinity of both IRAS sources. The parameters of YSOs in the IRAS clusters (124 objects) and 394 non-cluster objects surrounding them show some differences. About 75% of the YSOs belonging to the IRAS clusters have an evolutionary age greater than $10^6$ yr. Their slope $\alpha$ of the KLF agrees well with a Salpeter-type initial mass function (IMF) ( $\gamma = 1.35$ ) for a high mass range (O–F stars, $\beta \sim 2$ ) at 1 Myr. The non-cluster objects are uniformly distributed in the molecular cloud, 80% of which are located to the right of the 0.1 Myr isochrone. The slope $\alpha$ of the KLF of non-cluster objects is $0.55\,\pm\,0.09$ , corresponding better to a Salpeter-type IMF for low-mass objects (G–M stars, $\beta \sim 1$ ). Our results show that two dense stellar clusters are embedded in these two physically connected UC HII regions. The clusters include several high- and intermediate-mass zero-age main sequence stellar objects. Based on the small age spread of the stellar objects, we suggest that the clusters originate from a single triggering shock. The extended emission observed in both UC HII regions is likely due to the stellar clusters.

We report the results of a sensitive search for water maser emission in the Local Group Galaxy NGC 6822 with the Karl G. Jansky Very Large Array. The observations provide tentative single-epoch detections of four candidates, associated with two infrared-bright star formation regions (Hubble I/III and Hubble IV). The candidate maser detections are all offset from the velocity range where strong emission from Hi neutral gas is observed towards NGC 6822, with the closest offset by $\sim\!40\, \mathrm{kms}^{-1}$ . Our observations include the location of NL1K, a previous tentative water maser detection in NGC 6822. We do not detect any emission from this location with a sensitivity limit approximately a factor of 5 better than the original Sardina Radio Telescope observations.

We describe a new polarised imaging pipeline implemented in the fhd software package. The pipeline is based on the optimal mapmaking imaging approach and performs horizon-to-horizon image reconstruction in all polarisation modes. We discuss the formalism behind the pipeline’s polarised analysis, describing equivalent representations of the polarised beam response, or Jones matrix. We show that, for arrays where antennas have uniform polarisation alignments, defining a non-orthogonal instrumental polarisation basis enables accurate and efficient image reconstruction. Finally, we present a new calibration approach that leverages widefield effects to perform fully polarised calibration. This analysis pipeline underlies the analysis of Murchison Widefield Array data in Byrne et al. (2022, MNRAS, 510, 2011).

The Murchison Widefield Array (MWA) is a low-frequency aperture array capable of high-time and frequency resolution astronomy applications such as pulsar studies. The large field-of-view of the MWA (hundreds of square degrees) can also be exploited to attain fast survey speeds for all-sky pulsar search applications, but to maximise sensitivity requires forming thousands of tied-array beams from each voltage-capture observation. The necessity of using calibration solutions that are separated from the target observation both temporally and spatially makes pulsar observations vulnerable to uncorrected, frequency-dependent positional offsets due to the ionosphere. These offsets may be large enough to move the source away from the centre of the tied-array beam, incurring sensitivity drops of ${\sim}30{-}50\%$ in Phase II extended array configuration. We analyse these offsets in pulsar observations and develop a method for mitigating them, improving both the source position accuracy and the sensitivity. This analysis prompted the development of a multi-pixel beamforming functionality that can generate dozens of tied-array beams simultaneously, which runs a factor of ten times faster compared to the original single-pixel version. This enhancement makes it feasible to observe multiple pulsars within the vast field of view of the MWA and supports the ongoing large-scale pulsar survey efforts with the MWA. We explore the extent to which ionospheric offset correction will be necessary for the MWA Phase III and the low-frequency square kilometre array (SKA-low).

We study the radio power of the core and its relation to the optical properties of the host galaxy in samples of high-excitation (HERG) and low-excitation (LERG) Fanaroff–Riley type II (FRII) radio galaxies. The radio galaxy sample is divided into two groups of core/non-core FRII, based on the existence of strong, weak or lack of single radio core component. We show that FRII LERGs with radio emission of the core have significantly higher [O III] line luminosities compared to the non-core LERG FRIIs. There is no significant difference between the hosts of the core and non-core FRIIs of LERG type in galaxy sizes, concentration indices, star formation rates, 4000-Å break strengths, colours, black hole masses, and black hole to stellar masses. We show that the results are not biased by the stellar masses, redshifts, and angular sizes of the radio galaxies. We argue that the detection of higher [O III] luminosities in the core FRIIs may indicate the presence of higher amounts of gas, very close to the active galactic nuclei (AGN) nucleus in the core FRIIs compared to the non-core FRIIs or may result from the interaction of the radio jets with this gas. The core and non-core FRIIs of the HERG type show no significant differences perhaps due to our small sample size. The effect of relativistic beaming on the radio luminosities and the contribution of restating AGN activity have also been considered.

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries.

Most massive stars (up to 100%) are thought to be in binary systems. The multiplicity of massive stars seems to be intrinsically linked to their formation and evolution, and so Massive Young Stellar Objects are key in observing this early stage of star formation. We have surveyed hundreds of MYSOs across the Galaxy from the RMS catalogue, using UKIDSS and VVV point source data. Preliminary results show binary fractions of 44±3% for the UKIDSS sample and 32±3% for the VVV sample. In addition we use the K-band magnitudes as a proxy for the companion mass, and find a significant fraction of the detected companions have estimated mass ratios greater than 0.5, which suggests a deviation from the capture formation scenario.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected spectra of about 100,000 stars. The survey provides not only the reduced spectra, but also the radial velocities, stellar parameters and surface abundances resulting from the analysis of the spectra. We present the work of the groups that analysed the spectra of the hottest stars in that Survey. The large temperature range that is covered (Teff = 7,000 to 50,000 K) requires the use of different analysis codes by the different groups. Eight groups each analysed part of the data, with significant overlap that allowed cross-checks. In total 17,693 spectra of 6,462 stars were analysed, most of them in 37 open star clusters. The homogenisation of all this information led to stellar parameters for 5,584 stars. Abundances for at least one of the elements He, C, N, O, Ne, Mg, Al, Si and Sc were determined for 292 stars. The GES hot star data, as well as the Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters.

Over the last two decades there have been considerable advances in modelling the spectra of massive stars and supernovae (SNe). Despite this progress, there are still numerous uncertainties that affect the accuracy of models. For massive stars, convection, instabilities, clumping, and our inability to model stellar winds self-consistently likely introduce systematic errors into our analyses. For SNe, and particularly for core-collapse SNe, departures from spherical symmetry strongly affect observed spectra and need to be taken into account. There are also issues with clumping, and mixing processes (both in the progenitor and the SN explosion) that need to be resolved. For both massive stars and SNe, the accuracy and availability of atomic data continues to be an ongoing issue influencing analyses.

In recent years, it has been discovered that massive stars commonly exhibit a non-coherent form of variability in their light curves referred to as stochastic low frequency (SLF) variability. Various physical mechanisms can produce SLF variability in such stars, including stochastic gravity waves excited at the interface of convective and radiative regions, dynamic turbulence generated in the near-surface layers, and clumpy winds. Gravity waves in particular are a promising candidate for explaining SLF variability as they can be ubiquitously generated in main sequence stars owing to the presence of a convective core, and because they provide the large-scale predominantly tangential velocity field required to explain macroturbulence in spectral line fitting. Here, I provide an overview of the methods and results of studying SLF variability in massive stars from time series photometry and spectroscopy.

In discussing open question in the field of massive stars, I consider their evolution from birth to death. After touching upon massive star formation, which may be bi-modal and not lead to a zero-age main sequence at the highest masses, I consider the consequences of massive stars being close to their Eddington limit. Then, when discussing the effects of a binary companion, I highlight the importance of massive Algols and contact binaries for understanding the consequences of mass transfer, and the role of binaries in forming Wolf-Rayet stars. Finally, a discussion on pair instability supernovae and of superluminous supernovae is provided.

The Alicante Survey of MAssive Stars in Hii Regions (A-SMASHeR) is aimed at finding the ratio of massive stars that are born in isolation. We present LIRIS/WHT images and EMIR/GTC spectra of the massive stellar content in A-SMASHeR regions. Our preliminary analysis yields ∼20% of regions hosting relatively (or truly) isolated massive stars.

We perform spectral fittings for O-type stars based on self-consistent wind solutions, providing Ṁ and ν(r directly derived from the initial stellar parameters. We introduce our two methods: m-CAK prescription and Lambert-procedure.

The Lambert-procedure allows the calculation of consistent v(r) that reduce the number of free parameters when a spectral fitting using CMFGEN is performed, even without recalculation of the Ṁ. Spectra calculated from our Lambert-solutions show significant differences compared to the initial β-law CMFGEN models. For m-CAK prescription, self-consistent solutions provide values for theoretical Ṁ on the order of the most recent predictions from other studies. Later, we find a global fit with the RT code FASTWIND. This is an important step towards the determination of stellar and wind parameters without using β-law. Our m-CAK prescription is valid for the O-type stars with Teff ≥ 30 kK and log g ≥ 3.2.

We expect that solutions introduced here to be extended to numerous studies about massive stars in future.

An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/optical spectroscopic libraries for population synthesis models and a large number of interstellar sight-lines towards the Magellanic Clouds.

Wolf-Rayet stars are regarded as candidates for progenitors of core-collapse supernovae, and they are expected to be progenitors of long gamma-ray bursts. These types of stars are considered to be fast rotators. Their high rotation speed breaks the sphericity of the star and leads to an axisymmetric wind density structure. In such a case, the electron scattering takes place in a nonspherical environment, and as a result, we might expect an intrinsic polarization. We present a 2.5D radiation hydrodynamic stellar wind model of these stars. The model simulations account for the deformation of the stellar surface due to rotation, gravity darkening, and nonradial forces. We computed the polarization from the density variable of the hydrodynamic model, derived the upper limit of rotational velocities, and found no conflict with the previous studies of Wolf-Rayet stars.

Gravitational-wave (GW) observations are revealing the population of compact objects from a new angle. Yet their stellar progenitors remain uncertain because few observational clues on their progenitors exist. Theoretical models typically assume that the progenitor evolution can be approximated with single-star models. We explore how binary evolution affects the pre-supernova (SN) structure of stars, and the resulting distribution of compact object remnants. We focus on the differences in the core properties of single stars and of donor stars that transfer their outer layers in binary systems and become binary-stripped. We show that the final structures of binary-stripped stars that lose their outer layers before the end of core helium burning are systematically different compared to single stars. As a result, we find that binary-stripped stars tend to explode more easily than single stars and preferentially produce neutron stars and fewer black holes, with consequences for GW progenitors.