Binary neutron star (NS) mergers have been expected to synthesize r-process elements and cause electromagnetic radiation called kilonovae. Although r-process nucleosynthesis was confirmed by the observations of GW170817/AT2017gfo, individual elements have not been identified except for strontium. Toward identification of elements in kilonova spectra, we perform radiative transfer simulations in NS merger ejecta. We find that Sr II triplet lines appear in the spectrum, which is consistent with the absorption feature observed in GW170817/AT2017gfo. The synthetic spectrum also shows the strong Ca II triplet lines. Absence of the Ca II line features in GW170817/AT2017gfo implies that the Ca/Sr ratio is < 0.002 in mass fraction, which is consistent with nucleosynthesis for electron fraction ≥ 0.40 and entropy per nucleon (in units of Boltzmann constant) ≥ 25. Identification of absorption lines in near-infrared wavelengths which have not yet been decoded may lead to clarify the abundances synthesized in NS merger ejecta.

The formation of the first galaxies in the Universe is the new frontier of both galaxy formation and reionization studies. This creates a fierce new challenge, i.e. to simultaneously understand in a unique and coherent picture the processes of galaxy formation and reionization, and – crucially – their connection. To this end, we present the thesan suite of cosmological radiation-magneto-hydrodynamical simulations. They are unique since they: (i) cover a very broad range of spatial and temporal scales; (ii) include an unprecedentedly-broad range of physical processes for simulations of such scales and resolution; (iii) exploit knowledge accumulated at low redshift to minimize the number of free parameters in the physical model; (iv) use a variance-suppression technique in the production of initial conditions to increase their statistical fidelity. Finally, the thesan suite includes multiple runs of the same initial conditions, exploring current unknowns in the physics of dark matter and ionizing sources.

IXPE is a NASA/ASI Small Explorer Mission. It will probe the X-ray polarization properties of celestial sources. In particular, for Accretion-powered Millisecond Pulsars (AMPs), IXPE can provide us with unique information on their geometry. These information together with pulse shape modelling will strongly boost the achievable sensitivity on measuring the AMPs mass and the radius. As a case-study, we simulated an observation of SAXJ1808.4-3658 and studied the accuracy that can be achieved in the measured time-dependent Stokes profiles. From these data we estimated how well IXPE will be able to constraint NS geometrical parameters, such as the inclination and hot-spot co-latitude angles.

The purpose of the present work is a detailed investigation of the dynamical evolution of Collinder 135 and UBC 7 star clusters. We present a set of dynamical numerical simulations using realistic star cluster -body modeling technique with the forward integration of the star-by-star cluster models to the present day, based on best-available 3D coordinates and velocities obtained from the latest Gaia EDR3 data release. We have established that Collinder 135 and UBC 7 are probably a binary star cluster and have common origin. We carried out a full star-by-star N-body simulation of the stellar population of both clusters using the new algorithm of Single Stellar Evolution and performed a comparison of the results obtained in the observational data (like cumulative number counts), which showed a fairly good agreement.

The convection-enhanced neutrino-driven supernova engine’s success in explaining a myriad of supernova properties has set it as the standard engine behind supernova. However, due to the success of rotationally-powered engines in explaining astrophysical transients like gamma-ray bursts, these engines have been revived as possible drivers of normal supernovae, competing with this standard engine. In this paper, these competing engines, and the constraints placed by compact remnant observations on these engines, are reviewed. We find that, with these constraints, such rotationally-powered engines can explain less than 1% of the current supernova remnants. In addition, we find that the remnant mass distribution can be used to constrain properties of the convection-enhanced neutrino-driven engine, helping astronomers understand the nature of convection in this engine.

Five long gamma-ray bursts (GRBs) have been found to have very high energy (VHE, > 100GeV) counterparts. Interestingly, more than one emission mechanism has been invoked to explain the VHE counterpart from different events. As a result of this discovery, it has become apparent that we have been missing half of the energy produced in the afterglow of GRBs. We have been studying the radio afterglows in order to investigate whether these VHE GRBs have unusual jet properties. Studying these events in the radio waveband is advantageous as the emission at lower frequencies is brighter for longer enabling detailed, long term study of the jet evolution. The jet properties and environments of these GRBs vary hugely in a similar manner to that seen in the ‘regular’ long GRB population with evidence of bright reverse shock emission and multiple jet components. This work is presented on behalf of a much larger collaboration.

In the present work, we have developed a two-dimensional gravitational model of barred galaxies to analyse the fate of escaping stars from the central barred region. For that, the model has been analysed for two different bar profiles viz. strong and weak. Here the phenomena of stellar escape from the central barred region have been studied from the perspective of an open Hamiltonian dynamical system. We observed that the escape routes correspond to the escape basins of the two index-1 saddle points. Our results show that the formation of spiral arms is encouraged for the strong bars. Also, the formation of grand design spirals is more likely for strong bars if they host central super massive black holes (SMBHs). In the absence of central SMBHs, the formation of less-prominent spiral arms is more likely. Again, for weak bars, the formation of inner disc rings is more probable.

T Tauri Stars (TTSs) offer a unique chance to study the physics of non-relativistic accretion engines. In this invited talk, the current status of the field is presented with special emphasis on the predictive power of the numerical simulation of magnetospheric accretion and close binary systems and its impact on astronomical observations.

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.

We report the initial results of deep eROSITA monitoring of the magnificent seven isolated neutron stars (INSs). Thanks to a combination of high count statistics and good energy resolution, the eROSITA datasets unveil the increasingly complex energy distribution of these presumably simple thermal emitters. For three targets, we report the detection of multiple (in some cases, phase-dependent) spectral absorption features and deviations from the dominant thermal continuum. Unexpected long-term changes of spectral state and timing behaviour have additionally been observed for two INSs. The results pose challenging theoretical questions on the nature of the variations and absorption features and ultimately impact the modeling of the atmosphere and cooling of highly magnetised neutron stars.

We trained a Neural Network that can obtain selected STARLIGHT parameters directly from S-PLUS photometry. The training set consisted of over 55 thousand galaxies with their stellar population parameters obtained from a STARLIGHT application by Cid Fernandes et al. (2005). These galaxies were crossmatched with the S-PLUS iDR 3 database, thus, recovering the photometry for the 12 band filters for 55803 objects. We also considered the spectroscopic redshift for each object which was obtained from the SDSS. Finally, we trained a fully connected Neural Network with the 12-band photometry + redshift as features, and targeted some of the STARLIGHT parameters, such as stellar mass and mean stellar age. The model performed very well for some parameters, for example, the stellar mass, with an error of 0.23 dex. In the future, we aim to apply the model to all S-PLUS galaxies, obtaining never-before-seen photometric synthesis for most objects in the catalogue.

We provide analysis of the baryon asymmetry generated in the Scalar Field Condensate (SFC) baryogenesis model obtained in new inflation, chaotic inflation, Starobinsky inflation, MSSM inflation, quintessential inflation, considering both cases of efficient thermalization after inflation and also delayed thermalization. We have found that baryon asymmetry generated in SFC baryogenesis model is considerably bigger than the observed one for the new inflation, new inflation model by Shafi and Vilenkin, MSSM inflation, chaotic inflation with high reheating temperature and the simplest Shafi-Vilenkin chaotic inflationary model. Therefore, strong diluting mechanisms are needed to reduce the baryon excess to its observational value today for these models. We have shown that for the SFC baryogenesis model a successful generation of the observed baryon asymmetry is possible in Modified Starobinsky inflation, chaotic inflation with low reheating temperature, chaotic inflation in SUGRA and quintessential inflationary model.

During the late stages of a neutron star binary inspiral finite-size effects come into play, with the tidal deformability of the supranuclear density matter leaving an imprint on the gravitational-wave signal. As demonstrated in the case of GW170817—the first direct detection of gravitational waves from a neutron star binary—this can lead to strong constraints on the neutron star equation of state. As detectors become more sensitive, effects which may have a smaller influence on the neutron star tidal deformability need to be taken into consideration. Dynamical effects, such as oscillation mode resonances triggered by the orbital motion, have been shown to contribute to the tidal deformability, especially close to the neutron star coalesence. We calculate the contribution of the various stellar oscillation modes to the tidal deformability and demonstrate the (anticipated) dominance of the fundamental mode, showing what the impact of the matter composition is on the tidal deformability.

We have studied the input of the exothermic photochemistry into the formation of the non-thermal escape flux in the transition H2 − H region of the extended upper atmosphere of the hot exoplanet - the sub-neptune π Men c. The formation rate and the energy spectrum of hydrogen atoms formed with an excess of kinetic energy due to the exothermic photochemistry forced by the stellar XUV radiation were calculated using a numerical kinetic Monte Carlo model of a hot planetary corona. The escape flux was estimated to be equal to 2.5×1012 cm−2s−1 for the mean level of stellar activity in the XUV radiation flux. This results in the mean estimate of the atmospheric loss rate due to the exothermic photochemistry equal to 6.7×108 g s−1. The calculated estimate is close to the observational estimates of the possible atmospheric loss rate for the exoplanet π Men c in the range less than 1.0×109 gs−1.

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?

The effect of a parallel velocity shear on the explosive phase of magnetic reconnection in a double tearing mode is investigated within the 2D resistive magneto-hydrodynamic framework. All the systems follow a three phase evolution pattern with the phases delayed in time for an increasing shear speed. We find that the theoretical dependence of the reconnection rate with shear remains true in more general scenarios such as that of a plasmoid dominated double current sheet system. We also find that the power-law distribution of plasmoid sizes become steeper with an increasing sub-Alfvénic shear. We further demonstrate the effect of a velocity shear on acceleration of test particles pertaining to the modification in the energy spectrum.

Astrophysical systems possess various sites of particle acceleration, which gives rise to the observed non-thermal spectra. Diffusive shock acceleration (DSA) and stochastic turbulent acceleration (STA) are the candidates for producing very high energy particles in weakly magnetized regions. While DSA is a systematic acceleration process, STA is a random energization process, usually modelled as a biased random walk in energy space with a Fokker-Planck equation. In astrophysical systems, different acceleration processes work in an integrated manner along with various energy losses.

Here we study the interplay of both STA and DSA in addition to various energy losses, in a simulated RMHD jet cocoon. Further, we consider a phenomenologically motivated STA timescale and discuss its effect on the emission profile of the RMHD jet. A parametric study on the turbulent acceleration timescale is also conducted to showcase the effect of turbulence damping on the emission structure of the simulated jet.

We investigate the combined evolution of the dipolar surface magnetic field (Bs) and the spin-period (Ps) of known Magnetars and high magnetic field ($${{\rm{B}}_s} \mathbin{\lower.3ex\hbox{$\buildrel>\over {\smash{\scriptstyle\sim}\vphantom{_x}}$}} {10^{13}}{\rm{G}}$$) radio pulsars. We study the long term behaviour of these objects assuming a simple Ohmic dissipation of the magnetic field. Identifying the regions (in the Ps-Bs plane) in which these neutron stars would likely move into, before crossing the death-line to enter the pulsar graveyard, we comment upon the possible connection between the Magnetars and other classes of neutron stars.