The cosmic origin of fluorine is still under debate. Asymptotic giant branch (AGB) stars are among the few suggested candidates to efficiently synthesis F in our Galaxy, however their relative contribution is not clear. In this paper, we briefly review the theoretical studies from stellar yield models of the F synthesis and chemical equilibrium models of the F-containing molecules in the outflow around AGB stars. Previous detections of the F-bearing species towards AGB and post-AGB stars are also highlighted. We suggest that high-resolution ALMA observations of the AlF, one of the two main carriers of F in the outflow of AGB stars, can provide a reliable tracer of the F-budget in AGB stars. This will be helpful to quantify the role of AGB stars in the Galactic F budget.

Binary interaction with a stellar or planetary companion has been proposed to be the driving mechanism behind large-scale asymmetries, such as spirals and disks, observed within AGB outflows. We developed the first chemical kinetics model that takes the effect of a stellar companions’s UV radiation into account. The presence of a stellar companion can initiate a rich photochemistry in the inner wind. Its impact is determined by the intensity of the UV radiation and the extinction the radiation experiences. The outcome of the inner wind photochemistry depends on the balance between two-body reactions and photoreactions. If photoreactions dominate, the outflow can appear molecule-poor. If two-body reactions dominate, chemical complexity within the outflow can increase, yielding daughter species with a large inner wind abundance. A comprehensive view on the molecular content of the outflow, especially combined with abundance profiles, can point towards the presence of a stellar companion.

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.

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.

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 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.

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.

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.

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.

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.

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.

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.

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.

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 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.

The Sun moves with respect to the local interstellar medium (LISM) and modifies its properties to heliocentric distances as large as 1 pc. The solar wind (SW) is affected by penetration of the LISM neutral particles, especially H and He atoms. Charge exchange between the LISM atoms and SW ions creates pickup ions (PUIs) and secondary neutral atoms that can propagate deep into the LISM. Neutral atoms measured at 1 au can provide us with valuable information on the properties of pristine LISM. Voyager 1 and 2 spacecraft perform in situ measurements of the LISM perturbed by the presence of the heliosphere and relate them to the unperturbed region. We discuss observational data and numerical simulations that shed light onto the mutual influence of the SW and LISM. Physical phenomena accompanying the SW–LISM interaction are discussed, including the coupling of the heliospheric and interstellar magnetic field at the heliopause.

Transit observations have revealed the existence of atmospheric escape in several hot Jupiters. High energy photons from the host star heat the upper atmosphere and drive the hydrodynamic escape. The escaping atmosphere can interact with the stellar wind from the host star. We run radiation hydrodynamics simulations with non-equilibrium chemistry to investigate the wind effects on the escape and the transit signature. Our simulations follow the planetary outflow driven by the photoionization heating and the wind interaction in a dynamically coupled, self-consistent manner. We show that the planetary mass-loss rate is almost independent of the wind strength, which however affects the Ly-α transit depth considerably. But the Hα transit depth is almost independent of the wind strength because it is largely caused by the lower hot layer. We argue that observations of both lines can solve the degeneracy between the EUV flux from the host and the wind strength.