The post-main sequence evolutionary path of massive stars comprises various transition phases, in which the stars shed large amounts of material into their environments. Our studies focus on two of them: B[e] supergiants and yellow hypergiants, for which we investigate the structure and dynamics within their environments. We find that each B[e] supergiant is surrounded by a unique set of rings or arc-like structures. These structures are either stable over time or they display high variability, including expansion and dilution. In contrast, yellow hypergiants are embedded in multiple shells of gas and dust. These objects are famous for their outburst activity. Moreover, the dynamics in their extended atmospheres imply an enhanced pulsation activity prior to outburst. The physical mechanism(s) leading to episodic mass ejections in these two types of stars is still uncertain. We propose that strange-mode instabilities, excited in the inflated envelopes of these objects, play a significant role.

In the Galactic center, there are many massive stars blowing strong stellar winds, which will strongly influence the surrounding environment and even the Galactic feedback. The Galactic center is quiescent at present, so the unique continuous energy input source is the massive star, consequently giving rise to many special features, such as the radio bubbles, the X-ray chimneys, the non-thermal filaments and high-metallicity abundance. However, it is difficult to quantify their contributions due to the complex environment in this region, and the past supernovae and Sgr A* activity are also important factors shaping these features. In this work, we discuss some structures possibly related to the stellar winds and perform preliminary simulations to study their evolution. We conclude the stellar winds can obviously influence a large scale ∼ 100 pc, and can possibly influence a larger scale environment indirectly.

Fast Radio Bursts (FRBs) are millisecond-long bursts of radio emission of extragalactic origin. The nature or FRBs is still unknown. Whether all FRBs are representatives of the same source population, or whether multiple underlying populations exist, is also unknown. One class that stands out is that of the “repeaters”, i.e. FRBs from which multiple bursts have been detected. In these cases, appropriate models should be non-cataclysmic but yet being able to create powerful coherent radio emission. Magnetars are among those source types that are considered as possible explanation for (repeating) FRBs. This review will summarise the basic properties of FRBs and those of magnetars to provide a critical assessment of the possible physical connection between these classes of sources. We conclude that while magnetars may indeed be related to the FRB phenomenon, it is unlikely that they explain all FRBs, i.e. at least two classes of FRBs exist.

During the recombination of the universe, supersonic relative motion between baryons and dark matter (DM) generally existed. In the presence of such streaming motions, gas clumps can collapse outside of virial radii of their closest dark matter halos. Such baryon dominant objects are thought to be self-gravitating and are called supersonically induced gas objects; SIGOs. We perform three-dimensional hydrodynamical simulations by including H2 chemical reactions and stream velocity and follow SIGO’s formation from z = 200 to z = 25. SIGOs can be formed under the influence of stream velocity, and cooling is effective in contracting gas clouds. We follow its further evolution with higher resolution. We find that there are SIGOs which become Jeans unstable outside of the virial radius of the closest DM halos. Those SIGOs are gravitationally unstable and trigger star formation.

We use hydrodynamical simulations to study the evolution of baryonic gas in a cosmologically evolving dark matter halo. We model both the inner and outer regions of the halo using a density profile that transitions from an inner NFW profile to a flat profile far from the halo. Metallicity-dependent radiative cooling and AGN jet feedback are implemented, which lead to heating and cooling cycles in the core. We analyze the evolution of gas and the central supermassive black hole (SMBH) across cosmological time. We find that the properties of the gas and the SMBH are correlated across halo masses and feedback efficiencies.

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 argue that measurements of X-ray polarization using the recently launched Imaging X-ray Polarimetry Explorer will answer many open questions about magnetars in particular the physical state of their surfaces, whether vacuum birefringence exists, and the nature of the hard X-ray emission from these objects. We outline the capabilities of the instrument, specific models and the results of simulations for the magnetar 4U 0142+61.

Most stars with birth masses larger than that of our Sun belong to binary or higher order multiple systems. Similarly, most stars have stellar winds. Radiation pressure and multiplicity create outflows of material that remove mass from the primary star and inject it into the interstellar medium or transfer it to a companion. Both have strong impact on the subsequent evolution of the stars, yet they are often studied separately. In this short review, I will sketch part of the landscape of the interplay between stellar winds and binarity. I will present several examples where binarity shapes the stellar outflows, providing new opportunities to understand and measure mass loss properties. Stellar winds spectral signatures often help clearly identifying key stages of stellar evolution. The multiplicity properties of these stages then shed a new light onto evolutionary connections between the different categories of evolved stars.

The Bisous model is a tool that uses stochastic methods to detect the network of galactic filaments. This model is explicitly developed to detect the structure from observational data, using only galaxy positions as input. This paper shows that the Bisous model gives reliable results and including photometric data improves the resulting filamentary network. We used MultiDark-Galaxies catalogue to create a mock with photometric redshifts and samples with different galaxy number densities. We found that the filaments detected with the Bisous model are reliable; 85% of the detected filaments are unchanged compared to results with more complete input data. Adding photometric data improves the fraction of galaxies in filaments. Using the confusion matrix technique, we found the false discovery rate to always be below 5% when using photometric data.

GW170817, the merger of two neutron stars witnessed through both its gravitational wave siren and its glow at all wavelengths of light, represents the first multi-messenger detection of a compact binary merger. The association of the GW in-spiral signal from GW170817 with a γ-ray burst, a kilonova, and a non-thermal afterglow spanning all bands of the electromagnetic spectrum, has provided rich constraints on the physics and astrophysics of neutron stars. Starting from the example of GW170817, I briefly summarize recent results on observations of electromagnetic afterglows from gravitational wave triggers. In the light of these results, I highlight some key questions that are yet to be answered after the GW170817 discovery. I conclude by commenting briefly on some opportunities that lie in front of us, as improvements in ground-based gravitational wave detectors’ sensitivities will transform a trickle of multi-messenger discoveries into a flood, bringing the field of gravitational wave astronomy from its infancy to its maturity.

The link between hot and cool stellar outflows is shown to be critical for correctly predicting the masses of the most massive black holes (BHs) below the so-called pair-instability supernova (PISN) mass gap. Gravitational Wave (GW) event 190521 allegedly hosted an “impossibly” heavy BH of 85 M ⊙. Here we show how our increased knowledge of both metallicity Z and temperature dependent mass loss is critical for our evolutionary scenario of a low-Z blue supergiant (BSG) progenitor of an initially approx 100 M ⊙ star to work. We show using MESA stellar evolution modelling experiments that as long as we can keep such stars above 8000 K such low-Z BSGs can avoid strong winds, and keep a very large envelope mass intact before core collapse. This naturally leads to the Cosmic Time dependent maximum BH function below the PISN gap.

GW170817/GRB 170817A has had a huge impact on our understanding of gamma-ray burst (GRB) afterglows, and has prompted a huge sustained effort at modeling the details of the geometry and emission from GRB jets. While no additional electro-magnetic counterparts have been detected to gravitational wave emission from neutron star mergers so far, it is certainly reasonable to expect further detections in the future. Whether these will be very similar in nature to GRB 170817A or instead will provide us with samplings of afterglow model parameters across a wide parameter space remains an open question. In this presentation I will survey some of the work done by the various active groups worldwide in theoretical modeling and understanding afterglows post 170817A.

Galactic winds probe how feedback regulates the mass and metallicity of galaxies. Galactic winds have cold gas, which is mainly observable with absorption and emission lines. Theoretically studying how absorption lines are produced requires numerical simulations and realistic starburst UV backgrounds. We use outputs from a suite of 3D PLUTO simulations of wind-cloud interactions to first estimate column densities and temperatures. Then, to create synthetic spectra, we developed a python interface to link our PLUTO simulations to TRIDENT via the YT-package infrastructure. We produce UV backgrounds accounting for the star formation rate of starbursts. For this purpose, we use fluxes generated by STARBURST99, which are then processed through CLOUDY to create customised ion tables. Such tables are subsequently read into TRIDENT to generate absorption spectra. We explain how the various packages and tools communicate with each other to create ion spectra consistent with spectral energy distributions of starburst systems.

Evolved massive stars are major cosmic engines, providing strong mechanical and radiative feedback on their host environment. They contribute to the enrichment of their environment through a strong stellar winds, still poorly understood. Wind physics across the life cycle of these stars is the key ingredient to accomplish a complete understanding of their evolution in the near and distant Universe. Nowadays, the development of the observational instruments is so advanced that the observations became very sensitive to the details of the stellar surface making possible to quantitatively study what happens on their surfaces and above where the stellar winds become dominant. Three-dimensional radiative hydrodynamics simulations of evolved stars are essential to a proper and quantitative analysis of these observations. This work presents how these simulations have been (and will be) crucial to prepare and interpret a multitude of observations and how they are important to achieve the knowledge of the mass-loss mechanism.

Thermal energies deposited by OB stellar clusters in starburst galaxies lead to the formation of galactic superwinds. Multi-wavelength observations of starburst-driven superwinds pointed at complex thermal and ionization structures which cannot adequately be explained by simple adiabatic assumptions. In this study, we perform hydrodynamic simulations of a fluid model coupled to radiative cooling functions, and generate time-dependent non-equilibrium photoionization models to predict physical conditions and ionization structures of superwinds using the maihem atomic and cooling package built on the program flash. Time-dependent ionization states and physical conditions produced by our simulations are used to calculate the emission lines of superwinds for various parameters, which allow us to explore implications of non-equilibrium ionization for starburst regions with potential radiative cooling.

The spectacular detection of the first electromagnetic counterpart of a gravitational wave event detected by the LIGO/Virgo interferometers and originated by the coalescence of a double neutron star (NS) system (GW 170817) marked the dawn of a new era for astronomy. The short GRB 170817A associated to the gravitational wave event provided the long-sought evidence that at least a fraction of short GRBs are originated by NS-NS merging and suggested the intriguing possibility that relativistic jets can be launched in the process of a NS-NS merger. The wealth of data collected provided the first compelling observational evidence for the existence of kilonovae, i.e. the emission due to radioactive decay of heavy nuclei produced through rapid neutron capture. Besides the remarkable event associated to GW 170817, kilonova signatures have been identified in a few short GRBs light curves, supporting a scenario where kilonovae are ubiquitous and can probe neutron star mergers well beyond the horizon of the gravitational wave detectors. In this paper I will review the situation and perspectives of our understanding of short GRBs progenitors and kilonovae in the multi-messenger era.

Coalescence of binary neutron stars gives rise to kilonova, thermal emission powered by radioactive decays of newly synthesized r-process nuclei. Observational properties of kilonova are largely affected by bound-bound opacities of r-process elements. It is, thus, important to understand atomic properties of heavy elements to link the observed signals with nucleosynthesis of neutron star mergers. In this paper, we introduce the latest status of kilonova modeling by focusing on the aspects of atomic physics. We perform systematic atomic structure calculations of r-process elements to understand element-to-element variation in the opacities. We demonstrate that the properties of the atomic structure of heavy elements are imprinted in the opacities of the neutron star merger ejecta and consequently in the kilonova light curves and spectra. Using this latest opacity dataset, we briefly discuss implications for GW170817, expected diversity of kilonova emission, and prospects for element identification in kilonova spectra.

The focus of this work is to comprehensively understand hydro-dynamical back-flows and their role in dynamics and non-thermal spectral signatures particularly during the initial phase of X-shaped radio galaxies. In this regard, we have performed axisymmetric (2D) and three dimensional (3D) simulations of relativistic magneto-hydrodynamic jet propagation from tri-axial galaxies. High-resolution dynamical modelling of axisymmetric jets has demonstrated the effect of magnetic field strengths on lobe and wing formation. Distinct X-shape formation due to back-flow and pressure gradient of ambient is also observed in our 3D dynamical run. Furthermore, the effect of radiative losses and diffusive shock acceleration on the particle spectral evolution is demonstrated, which particularly highlights how crucial their contributions are in the emission signature of these galaxies. This imparts a significant effect on the galaxy’s equipartition condition, indicating that one must be careful in extending its use in estimating other parameters, as the criterion evolves with time.

From November 2019 to April 2020, the prototypical red supergiant Betelgeuse experienced an unexpected and historic dimming. This event was observed worldwide by astrophysicists, and also by the general public with the naked eye. We present here the results of our observing campaign with ESO’s VLT and VLTI in the visible and infrared domains. The observations with VLT/SPHERE-ZIMPOL, VLT/SPHERE-IRDIS, VLTI/GRAVITY and VLTI/MATISSE provide spatially resolved diagnostics of this event. Using PHOENIX atmosphere models and RADMC3D dust radiative transfer simulations, we built a consistent model reproducing the images and the photometry.

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.