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.

GW190425 is the second gravitational wave (GW) event caused by a binary neutron star (BNS) merger. We report the result of the follow-up observation of GW190425 by the Gravitational-wave Electromagnetic-wave (EM) Counterpart Korean Observatory (GECKO). Our observation demonstrates that GECKO can detect EM counterpart of a GW170817-like event.

Detection of Earth-mass planets with the radial velocity method requires a precision of about 10cm/s to identify a signal caused by such a planet. At the same time, noise originating in the photospheric and subphotospheric layers of the parent star is of the order of meters per second. Understanding the physical nature of the photospheric noise (so-called stellar jitter) and characterizing it are critical for developing techniques to filter out these unwanted signals. We take advantage of current computational and technological capabilities to create 3D realistic models of stellar subsurface convection and atmospheres to characterize the photospheric jitter. We present 3D radiative hydrodynamic models of several solar-type target stars of various masses and metallicities, discuss how the turbulent surface dynamics and spectral line characteristics depend on stellar properties, and provide stellar jitter estimates for these stars.

Photoelectric effect of dust grains by UV radiation is an important process for disk heating, but as a disk evolves, the amount of dust grains decreases. Photoeaporation is a disk dispersal process, which is caused by high-energy radiation. We perform a set of photoevaporation simulations solving hydrodynamics with radiative transfer and non-equilibrium chemistry in a self-consistent way. We run a series of simulations with varying the dust-to-gas mass ratio in a range . We show that H2 pumping and X-ray heating mainly contribute to the disk heating in case of and photoelectric effect mainly heats the gas in cases. The mass-loss profile changes significantly with respect to the main heating process. The outer disk is more efficiently dispersed when photoelectric effect is the main heating source.

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 introduce hydrodynamic simulations in which a protostar captures a cloudlet with a relatively small angular momentum. The cloudlet accretes onto the protostar and perturbs the gas disk rotating around the protostar. This cloudlet capture can reproduce some features observed in the molecular emission lines from TMC-1A. First, the cloudlet can reproduce the blue asymmetry observed in the CS emission. Second, the cloudlet can explain the slow infall observed in the C18O emission. Third, the impact of the cloudlet can explain the offset of the SO emission from the disk center. We also argue that a warm gas should confine the cloudlet through pressure. A protostar may obtain substantial mass by capturing cloudlets.

We present an overview of PION, an open-source software project for solving radiation-magnetohydrodynamics equations on a nested grid, aimed at modelling asymmetric nebulae around massive stars. A new implementation of hybrid OpenMP/MPI parallel algorithms is briefly introduced, and improved scaling is demonstrated compared with the current release version. Three-dimensional simulations of an expanding nebula around a Wolf-Rayet star are then presented and analysed, similar to previous 2D simulations in the literature. The evolution of the emission measure of the gas and the X-ray surface brightness are calculated as a function of time, and some qualitative comparison with observations is made.

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.

We analyzed 39 ks NuSTAR data of Cen X-3 through both orbital- and pulse-phase resolved spectroscopy. Orbital-phase resolved spectra show extrinsic fluctuations due to absorption by surrounding plasma, as the spectral fluctuation mainly emerges below 10 keV. Pulse-phase resolved spectra, on the other hand, show intrinsic fluctuations depending on effectiveness of Comptonization, since the spectrum becomes hard above 10 keV at the pulse peak.