One of the big challenges for 21st century stellar astrophysics is the impact of binary interactions on stellar evolution. Such interactions are believed to play a key role in the death throes of 1-8 M⊙ stars, as they evolve from the AGB stars into Planetary Nebulae. X-ray surveys of UV-emitting AGB stars show that ∼40% of objects with FUV emission and GALEX FUV/NUV flux ratios ≳0.2 have variable X-ray emission characterized by very high temperatures (Tx∼35-160 MK) and luminosities (Lx∼0.002-0.2L⊙). We hypothesize that such AGB stars have accretion and (accretion-powered) outflows associated with a close binary companion. UV spectroscopy with HST/STIS of our brightest object (Y Gem) shows the presence of infalling and outflowing gas, providing direct kinematic confirmation of this hypothesis. However, the UV-emitting AGB star population is dominated by objects with little or no FUV emission, and we do not know whether the UV emission from these is intrinsic to the AGB star or extrinsic (i.e., due to binarity). Here we present the first results from a large grid of simple chromospheric models to help discriminate between the intrinsic and extrinsic mechanisms of UV emission for AGB stars.

Thanks to the long-term collaborations between nuclear and astrophysics, we have good understanding on the origin of elements in the universe, except for the elements around Ti and some neutron-capture elements. From the comparison between observations of nearby stars and Galactic chemical evolution models, a rapid neutron-capture process associated with core-collapse supernovae is required. The production of C, N, F and some minor isotopes depends on the rotation of massive stars, and the observations of distant galaxies with ALMA indicate rapid cosmic enrichment. It might be hard to find very metal-poor or Population III (and dust-free) galaxies at very high redshifts even with JWST.

We run numerical simulations of massive colliding wind binaries, and quantify the accretion onto the secondary under different conditions. We set 3D simulation of a LBV–WR system and vary the LBV mass loss rate to obtain different values of wind momentum ratio η. We show that the mean accretion rate for stationary systems fits a power law Macc∝ η–1.6 for a wide range of η, until for extremely small η saturation in the accretion is reached. We find that the stronger the primary wind, the smaller the opening angle of the colliding wind structure (CWS), and compare it with previous analytical estimates. We demonstrate the efficiency of clumpy wind in penetrating the CWS and inducing smaller scale clumps that can be accreted. We propose that simulations of colliding winds can reveal more relations as the ones we found, and can be used to constrain stellar parameters.

Maser properties can be measured with milli-arcsec precision over multiple epochs using ALMA, cm- and mm-wave VLBI and e-MERLIN. This allows: (i) Tracing SiO maser proper motions in the pulsation-dominated zone; (ii) Quantifying clumpiness, variability and asymmetry of the wind traced by masers; (iii) Contrasting behaviour from OH masers even at similar distances from the star; (iv) Measuring magnetic fields. Mass lost from the star, traced by SiO masers, is likely to take decades to reach ∼5 stellar radii. At 5–50 stellar radii, once dust is well formed, 22-GHz H2O masers show the wind accelerating through the escape velocity; its overall direction is away from the star but the velocity field is complex. In a few cases (so far), highly-directed, localised ejecta are seen. Magnetic fields appear to be stellar-centred and strong enough to influence wind kinematics. Recent ALMA and other observations have shown that otherwise inconspicuous companions shape a majority of evolved star winds, whilst advanced models demonstrate how, for some situations, this is compatible with masers showing negligible rotation proper motions. The long-term monitoring achievable at radio frequencies complements the multi-transition maser studies and analysis of thermal lines and dust at shorter wavelengths.

A recent analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way (MW) identified a kink in the initial-final mass relation (IFMR), located over a range of initial masses, 1.65 ≲ Mi/M⊙ ≲ 2.10, which unexpectedly interrupts the commonly assumed monotonic trend. The proposed interpretation links this observational fact to the formation of carbon stars and the modest outflows (with mass loss rate < 10−7 M⊙/yr) that are expected as long as the carbon excess remains too low to produce dust grains in sufficient amount. Under these conditions the mass of the carbon-oxygen core can grow more than is generally predicted by stellar models. We discuss these new findings also in light of a new systematic follow-up investigation, based on Gaia EDR3, of evolved giants (13 carbon stars, 3 S stars and 4 M stars) belonging to intermediate-age open clusters.

Understanding the nucleosynthetic origin of nitrogen and the evolution of the N/O ratio in the interstellar medium is crucial for a comprehensive picture of galaxy chemical evolution at high-redshift because most observational metallicity (O/H) estimates are implicitly dependent on the N/O ratio. The observed N/O at high-redshift shows an overall constancy with O/H, albeit with a large scatter. We show that these heretofore unexplained features can be explained by the pre-supernova wind yields from rotating massive stars (M≳10M⊙,ν/νcrit≳0.4). Our models naturally produce the observed N/O plateau, as well as the scatter at low O/H. We find the scatter to arise from varying star formation efficiency. However, the models that have supernovae dominated yields produce a poor fit to the observed N/O at low O/H. This peculiar abundance pattern at low O/H suggests that dwarf galaxies are most likely to be devoid of SNe yields and are primarily enriched by pre-supernova wind abundances.

In this poster we present the structure of an axisymmetric, force-free magnetosphere of a twisted, aligned rotating dipole within a corotating plasma-filled magnetosphere. We explore various profiles for the twist. We find that as the current increases more field lines cross the light cylinder leading to more efficient spin-down. Moreover, we notice that the twist cannot be increased indefinitely and after a finite twist of about π/2 the field becomes approximately radial. This could have implications for torque variations of magnetars related to outbursts.

We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-calledMulti-particle collision (MPC) method, alternative to Fokker-Planck and Monte Carlo approaches. MPC is related to particle-mesh approaches for the computation of self consistent long-range fields, ensuring that simulation time scales with N log N in the number of particles, as opposed to N2 for direct N-body. The collisional relaxation effects are modelled by computing particle interactions based on a collision operator approach that ensures rigorous conservation of energy and momenta and depends only on particles velocities and cell-based integrated quantities.

So far detached compact binaries containing neutron stars have been observed either at intermediate stages of the evolution by radio telescopes or at merger by ground-based gravitational wave detectors. Sensitive to gravitational waves from binaries millions to thousands years prior to the merger, the future Laser Interferometer Space Antenna (LISA) will be crucial for bridging the gap between the currently accessible regimes. Depending on the binary type, LISA could potentially discover from a few to several hundreds in the entirely new regime throughout the Milky Way. Here we provide a concise summary of the current expectation for the detection of Galactic binaries containing neutron stars with LISA, focusing on double neutron stars and neutron star - white dwarf binaries. We outline a few examples of science investigations that LISA data will enable for these binaries.

We investigate the broad-band behaviour of circular polarization in radio pulsar profiles and show the relationship between polarization fraction and what proportion of that polarization is circular, both across frequency, and for a large number of pulsars viewed collectively. The behaviour observed may be explained by pulsar polarization originating from the partially-coherent combination of two linearly-polarized orthogonal modes with different flux spectral indices. (See also the poster in the “supplementary information”.)

The detection of an electromagnetic counterpart to the gravitational-wave source GW 170817 marked year zero of the multi-messenger gravitational-wave era. This event was generated by the merger of two neutron stars and gave rise to an electromagnetic transient, dubbed a “kilonova”. In this proceeding article, I will show how radiative transfer simulations can illuminate neutron star mergers and provide a connection between numerical models and observational data. I will present viewing-angle dependent kilonova predictions made with the Monte Carlo radiative transfer code POSSIS and show how these can be used to interpret data, place constraints on models and guide future follow-up campaigns of gravitational-wave triggers.

Most stars form in clumpy and sub-structured clusters. These properties also emerge in hydro-dynamical simulations of star-forming clouds, which provide a way to generate realistic initial conditions for N-body runs of young stellar clusters. However, producing large sets of initial conditions by hydro-dynamical simulations is prohibitively expensive in terms of computational time. We introduce a novel technique for generating new initial conditions from a given sample of hydro-dynamical simulations, at a tiny computational cost. In particular, we apply a hierarchical clustering algorithm to learn a tree representation of the spatial and kinematic relations between stars, where the leaves represent the single stars and the nodes describe the structure of the cluster at larger and larger scales. This procedure can be used as a basis for the random generation of new sets of stars, by simply modifying the global structure of the stellar cluster, while leaving the small-scale properties unaltered.

To form stars in hydrodynamical simulations, we introduce the grouped star formation prescription to convert the grouped sink particles into stars that follow the IMF. We show that this method is robust in different physical scales. Such methods to form stars are likely to become more important as galactic or even cosmological scale simulations begin to probe sub-parsec scales.

Scientific synergies between Athena and some of the key multi-messenger facilities that should be operative concurrently with Athena are presented. These facilities include LIGO A+, Advanced Virgo+ and future detectors for ground-based observation of gravitational waves (GW), LISA for space-based observations of GW, IceCube and KM3NeT for neutrino observations, CTA for very high energy observations. Multimessenger synergy science themes discussed here include pressing issues in the field of Astrophysics, Cosmology and Fundamental physics such as: the central engine and jet physics in compact binary mergers, accretion processes and jet physics in SMBBHs and in compact stellar binaries, the equation of state in neutron stars, cosmic accelerators and the origin of cosmic rays, the origin of intermediate and high-Z elements in the Universe, the Cosmic distance scale and tests of General Relativity and Standard Model. Observational strategies for implementing the identified science topics are also discussed.

We present our findings on magnetar spectral line analysis in the context of upcoming high resolution, high effective area, high throughput X-ray telescopes for two cases: persistent magnetar emission and magnetar bursts. For magnetars in quiescence, we present our preliminary work on modelling for phase-resolved emission. Our results reveal the necessity of constraining line depth and width concurrently with line energy to conclusively determine hotspot emission and corresponding magnetic field geometry. We then present the results of our simulations using effective area and response of various current and upcoming X-ray telescopes for magnetar spectral line detection and expand on the exciting opportunities upcoming telescopes provide to probe quiescent and burst emission region geometry and propagation in the extreme magnetic field of a magnetar.

The evolution of the magnetic field in neutron star crusts because of the Hall effect has received significant attention over the last two decades, which is strongly justified because of the dominance of this effect in highly magnetised neutron stars. However, the applicability of the Hall effect is based on the assumption that the crust does not fail and sustains its rigidity. This assumption can be violated for substantially strong magnetic fields. If this is the case, the evolution of the magnetic field is described by a different set of equations, which include the effects of a non-rigid crust. In this talk, after a brief review of the main characteristic of the Hall evolution, I will discuss the impact a plastic flow of the crust has on the magnetic field, studying axisymmetric models. Moreover, the way the crust fails impacts the overall evolution, with major differences appearing if the failure is local, intermediate or global. Quite remarkably, crustal failure and plasticity do not annul the Hall effect, and under certain circumstances they may even lead to a more dramatic evolution. I will discuss the impact of these effects in the context of neutron star timing behaviour, with special focus on timing noise, outbursts and glitches.

The INTEGRAL satellite, in orbit since October 2002, has significantly contributed to the study of magnetars and, thanks to its unique capabilities for the study of transient gamma-ray phenomena, it is now playing an important role in multimessenger astrophysics. The most recent results include the discovery of a peculiar burst from SGR J1935+2154, which gave the first observational evidence for the connection between magnetars and fast radio bursts, and extensive searches for bursting activity in peculiar sources, such as the repeating FRB 20200120E in M81 and the ultra-long period magnetar candidate GLEAMX J162759.5–523504.3.

The observational properties of isolated NSs are shaped by their magnetic field and surface temperature. They evolve in a strongly coupled fashion, and modelling them is key in understanding the emission properties of NSs. Much effort was put in tackling this problem in the past but only recently a suitable 3D numerical framework was developed. We present a set of 3D simulations addressing both the long-term evolution (≈ 104–106 yrs) and short-lived outbursts (≲ 1 yr). Not only a 3D approach allows one to test complex field geometries, but it is absolutely key to model magnetar outbursts, which observations associate to the appearance of small, inherently asymmetric hot regions. Even though the mechanism that triggers these phenomena is not completely understood, following the evolution of a localised heat injection in the crust serves as a model to study the unfolding of the event.

Motivated by their role as the direct or indirect source of many of the elements in the Universe, numerical modeling of core collapse supernovae began more than five decades ago. Progress toward ascertaining the explosion mechanism(s) has been realized through increasingly sophisticated models, as physics and dimensionality have been added, as physics and numerical modeling have improved, and as the leading computational resources available to modelers have become far more capable. The past five to ten years have witnessed the emergence of a consensus across the core collapse supernova modeling community that had not existed in the four decades prior. For the majority of progenitors – i.e., slowly rotating progenitors – the efficacy of the delayed shock mechanism, where the stalled supernova shock wave is revived by neutrino heating by neutrinos emanating from the proto-neutron star, has been demonstrated by all core collapse supernova modeling groups, across progenitor mass and metallicity. With this momentum, and now with a far deeper understanding of the dynamics of these events, the path forward is clear. While much progress has been made, much work remains to be done, but at this time we have every reason to be optimistic we are on track to answer one of the most important outstanding questions in astrophysics: How do massive stars end their lives?

Broadband X-ray data of five magnetars show that their hard X-ray pulses suffer periodic phase modulations, at a period ∼ 104 times their pulse period. The phenomenon is interpreted as a result of free precession of neutron stars that are prolately deformed to an asphericity of ∼ 10−4, by the magnetic stress of toroidal fields reaching ∼ 1016 G. The behavior is absent in their soft X-ray pulses, probably due to a higher emission symmetry. The ultra-high toroidal fields, considered common to magnetars, may persist longer than their dipole fields.