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.

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.

We present interferometric continuum and molecular line emission maps obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) of OH231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star that is key to understand the remarkable changes in nebular morphology and kinematics during the short transition from the AGB to the Planetary Nebula (PN) phase. The excellent angular resolution of our maps (∼20 mas ∼30 AU) allows us to scrutinize the central nebular regions of OH231.8+4.2, which hold the clues to understanding how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind (i.e. pre-PNe with massive bipolar outflows), the discovery of a rotating circumbinary disk of radius ∼30 AU traced by NaCl, KCl, and H2O emission lines. The disk lies at the base of a young bipolar wind with signs of rotation as well. A compact spatially resolved dust disk is found perpendicular to the bipolar outflow. We also identify a point-like continuum source, which likely represents the central Mira star enshrouded by a ∼3 R* shell or disk of hot (∼1400 K) freshly formed dust. The point source is slightly off-centre from the disk centroid, enabling us for the first time to place constraints to the orbital separation of the central binary system.

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.

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.

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

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.

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.

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.

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.

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.

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.

The detection of Earth-size exoplanets is a technological and data analysis challenge. Future progress in Earth-mass exoplanet detection is expected from the development of extreme precision radial velocity measurements. Increasing radial velocity precision requires developing a new physics-based data analysis methodology to discriminate planetary signals from host-star-related effects, taking stellar variability and instrumental uncertainties into account. In this work, we investigate and quantify stellar disturbances of the planet-hosting solar-type star HD121504 (G2V spectral type) from 3D radiative modeling obtained with the StellarBox code. The model has been used for determining statistical properties of the turbulent plasma and obtaining synthetic spectroscopic observations for several Fe I lines at different locations on the stellar disk to mimic high-resolution spectroscopic observations.