Differential rotation in neutron stars allows for significantly larger masses than rigid rotation. Some of those hypermassive objects are, however, unstable and collapse to a black hole immediately after formation. Yet, the exact threshold of dynamical stability is still unknown.

In our work we explore the limits on masses of neutron stars with various degrees of differential rotation which could be stable against a prompt collapse to a black hole by using turning-point (j-constant) criterion. We considered both spheroidal and quasi-toroidal differentially rotating neutron stars described by the polytropic equation of state. We find that massive configurations could be temporarily stabilized by differential rotation. Such objects are important sources of gravitational waves. Our results are a starting point for more detailed studies of stability using hydrodynamical codes.

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.

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 decay of the magnetic field in the interior of a magnetar may trigger electron captures by nuclei in the stellar crust, thus providing an internal source of heating. In turn, the onset of electron captures and the heat released are altered by the magnetic field due to the Landau–Rabi quantization of electron motion. The loss of magnetic pressure might also lead to pycnonuclear fusions of the lightest elements. The maximum amount of heat that can be possibly released by each reaction and their location are calculated using nuclear data from both experiments and theoretical predictions of the Brussels-Montreal models based on self-consistent Hartree-Fock-Bogoliubov calculations. Results are found to be consistent with those inferred empirically by comparing neutron-star cooling simulations with observed thermal luminosity of soft gamma-ray repeaters and anomalous X-ray pulsars.

Identification of the electromagnetic-wave (EM) counterparts of gravitational-wave (GW) sources can significantly broaden the research scope of GW astronomy, by pinpointing the exact locations of GW events and their environments, and using GW sources as standard sirens for cosmology. Yet, only one GW event has been found to be associated with an EM counterpart so far. Here, we outline the challenges of identifying EM counterparts of GW events, and describe our global network of telescopes that has been used to uncover GW EM counterparts. We also introduce a new facility in construction, the 7-dimensional telescope (7DT). Our GECKO observations have demonstrated that we can cover 50 deg2 within one hour to find kilonovae at a few hundred Mpc away. Furthermore, 7DT will produce a low resolution spectral map of the GW localization area, facilitating the EM counterpart search.

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.