The Radio Neutrino Observatory Greenland (RNO-G) is currently being deployed and it is currently gathering data. As a precursor and complementary detector to the future radio array of IceCube-Gen2 in Antarctica, it will explore mainly the Northern sky via in-ice radio detection technique. The total array configuration includes 35 radio stations and will be fully completed within three years from now. The antennas will register the radio signals produced by the Askaryan effect in cascades generated in ice by neutrinos. RNO-G’s scientific purpose is to detect UHE neutrinos at energies above 10 PeV. Due to the attenuation length of radio waves in ice (order of 1 km) the radio detection allows to address neutrino energies above several PeV. The detector will reach unprecedented sensitivity in the scale from tens of PeV up to EeV. Models predict GRBs induced by binary neutron star mergers as likely transient sources of such highly energetic neutrinos. The current study of NS-NS mergers will therefore possibly be complemented by future RNO-G detections through multimessenger temporal and spatial coincidence, including an alert system. In this presentation, we will describe the instrument capabilities and explore the possibility of detection of such sources with RNO-G.

The Earth’s atmosphere is incessantly bombarded by energetic charged particles called cosmic rays (CR) which are having either solar or non-solar origin. Analysis based on information theoretic estimators can be effectively employed as a potential technique to analyze the dynamical changes in cosmic ray intensity during different solar cycles. In the present study, dynamical complexity based analysis using Jensen-Shannon divergence (JSD) has been employed which reveals the existence of some peculiar fluctuation properties in CRI flux at Jung neutron monitor station. JSD based dynamical complexity analyses confirm the existence of difference in dynamical properties of CR flux during solar cycles 20-21 and 22-23.

Multiphase outflows driven by active galactic nuclei (AGN) have a profound impact on the evolution of their host galaxies. The effects of AGN feedback are especially prominent in the brightest cluster galaxies (BCGs) of cool-core clusters, where there is a concentration of gas in all phases, ranging from cold molecular gas to hot, >107 K ionized plasma. In this proceeding I describe recent simulation efforts to understand the formation and evolution of the 10-kpc-scale Hα-emitting filaments driven by AGN activities. Combined with observed star formation regions co-spatial with the filaments, this feedback mechanism can directly contribute to the growth of the central galaxy, albeit delayed by the characteristic radiative cooling timescale, ∼10 Myr, of the outflowing plasma.

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.

We compare the properties of shocked gas in Sgr B2 with maps obtained from 3D simulations of a collision between two fractal clouds. In agreement with 13CO(1-0) observations, our simulations show that a cloud-cloud collision produces a region with a highly turbulent density substructure with an average . Similarly, our numerical multi-channel shock study shows that colliding clouds are efficient at producing internal shocks with velocities of 5 − 50 km s−1 and Mach numbers of ∼ 4 − 40, which are needed to explain the ∼ 10−9 SiO abundances inferred from our SiO(2-1) IRAM observations of Sgr B2. Overall, we find that both the density structure and the shocked gas morphology in Sgr B2 are consistent with a Myr-old cloud-cloud collision. High-velocity shocks are produced during the early stages of the collision and can ignite star formation, while moderate- and low-velocity shocks are important over longer time-scales and can explain the extended SiO emission in Sgr B2.

Stars form in clusters, while planets form in gaseous disks around young stars. Cluster dissolution occurs on longer time scales than disk dispersal. Planet formation thus typically takes place while the host star is still inside the cluster. We explore how the presence of other stars affects the evolution of circumstellar disks. Our numerical approach requires multi-scale and multi-physics simulations where the relevant components and their interactions are resolved. The simulations start with the collapse of a turbulent cloud, from which stars with disks form, which are able to influence each other. We focus on the effect of extinction due to residual cloud gas on the early evolution of circumstellar disks. We find that this extinction protects circumstellar disks against external photoevaporation, but these disks then become vulnerable to dynamic truncation by passing stars. We conclude that circumstellar disk evolution is heavily affected by the early evolution of the cluster.

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.

Computational heliophysics has shed light on the fundamental physical processes inside the Sun, such as the differential rotation, meridional circulation, and dynamo-generation of magnetic fields. However, despite the substantial advances, the current results of 3D MHD simulations are still far from reproducing helioseismic inferences and surface observations. The reason is the multi-scale nature of the solar dynamics, covering a vast range of scales, which cannot be solved with the current computational resources. In such a situation, significant progress has been achieved by the mean-field approach, based on the separation of small-scale turbulence and large-scale dynamics. The mean-field simulations can reproduce solar observations, qualitatively and quantitatively, and uncover new phenomena. However, they do not reveal the complex physics of large-scale convection, solar magnetic cycles, and the magnetic self-organization that causes sunspots and solar eruptions. Thus, developing a synergy of these approaches seems to be a necessary but very challenging task.

The hosts of binary neutron star (BNS) mergers and hence short GRBs are not only galaxies with old stellar populations, infact most short GRB hosts are at least mildly star-forming galaxies. According to theoretical studies of merger populations, both short and long merging time-scales are expected. The immediate environments of BNS mergers are not as directly related to the property of the progenitor system as for long GRBs, since the system usually travelled a significant distance from their birth place. However, studying the stellar population properties across the host can still give us vital information on the contribution of formation channels and on merger timescales. Here we review the properties of NS merger hosts in emission using integrated-light and resolved observations. The afterglows of short GRBs furthermore serve to study the interstellar medium in their host galaxies in absorption. We present our best example to date, GRB 160410A at z = 1.7, one of the highest redshift short GRBs.

We join gravitational-wave and electromagnetic data to implement a combined simultaneous fit of the GW170817 event. The LIGO-Virgo analysis includes the estimation of the inclination, the angle of the binary with respect to the gravitationa-wave detector network line of sight. From the observations of the afterglow, instead, we can recover the viewing angle. The inclination and the viewing angle are supplementary angles, and can be treated as a single parameter. The value of the inclination that we recover from the fit is in agreement with the LIGO-Virgo previous works, with an uncertainty that is 10-fold smaller, thanks to contribution of the electromagnetic data. Moreover, with the inclusion of the gravitational-wave data, the degeneracy between the viewing angle and the jet opening angle is broken. This procedure is useful not only for analyzing GW170817, but any gravitational-wave event with an electromagnetic counterpart.

A four-year sky survey with the use of the eROSITA telescope on board the Spektr-RG observatory will provide the best coverage in the soft (0.5–2 keV) and standard (2–10 keV) X-ray ranges, both in terms of sensitivity and angular resolution. We have analysed the possibility of detecting various types of isolated neutron stars with eROSITA. Among already known objects, eROSITA will be able to detect more than 160 pulsars, 21 magnetars, 7 central compact objects, all seven sources of the Magnificent Seven, and two other X-ray isolated neutron stars during the four-year survey mission.

We present a new algorithm to solve the equations of radiation hydrodynamics (RHD) in a frequency-integrated, two-moment formulation. Novel features of the algorithm include i) the adoption of a non-local Variable Eddington Tensor (VET) closure for the radiation moment equations, computed with a ray-tracing method, ii) support for adaptive mesh refinement (AMR), iii) use of a time-implicit Godunov method for the hyperbolic transport of radiation, and iv) a fixed-point Picard iteration scheme to accurately handle the stiff nonlinear gas-radiation energy exchange. Tests demonstrate that our scheme works correctly, yields accurate rates of energy and momentum transfer between gas and radiation, and obtains the correct radiation field distribution even in situations where more commonly used – but less accurate – closure relations like the Flux-limited Diffusion and Moment-1 approximations fail. Our scheme presents an important step towards performing RHD simulations with increasing spatial and directional accuracy, effectively improving their predictive capabilities.

Multidimensional mathematical analysis, like Machine Learning techniques, determines the different features of objects, which is difficult for the human mind. We create a machine learning model to predict galaxies’ detailed morphology (∼ 300000 SDSS-galaxies with z < 0.1) and train it on a labeled dataset defined within the Galaxy Zoo 2 (GZ2). We use convolutional neural networks (CNNs) to classify the galaxies into five visual types (completely rounded, rounded in-between, smooth cigar-shaped, edge-on, and spiral) and 34 morphological classes attaining >94% of accuracy for five-class morphology prediction except for the cigar-shaped (∼ 87%) galaxies.

The velocity-space distribution of the solar neighborhood stars shows complex substructures (moving groups) including the well-known Hercules stream. Recently, the Gaia observation revealed their detailed structures, but their origins are still in debate. We analyzed a high-resolution N-body simulation of a Milky Way (MW)-like galaxy. To find velocity-space distributions similar to that of the solar neighborhood stars, we used Kullback-Leibler divergence (KLD), which is a metric to measure similarities between probability distributions. The KLD analysis shows the time evolution and the spatial variation of the velocity-space distribution. Velocity-space distributions with small KLDs (i.e. high similarities) are frequently but not always detected around in the simulated MW. In the velocity-map with smallest KLD, the velocity-space substructures are made from bar resonances.

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.

Based on our current high resolution direct N-body modelling of the Milky Way typical Star Cluster systems dynamical evolution we try to numerically estimate the influence of individual spin values and orientations on gravitational wave (GW) waveforms and observed time-frequency maps during multiple cycles for binary black hole (BBH) mergers. In our up to date N-body dynamical simulations we use the high order relativistic post-Newtonian corrections for the BH binary particles (3.5 post-Newtonian (PN) terms including spin-spin and spin-orbit terms). In the current work, we present the GW waveforms catalogue which covers the large parameter space in mass ratios 0.05 - 0.82 and extreme possible individual spin cases.

The magnetar SGR J1830–0645 was discovered in outburst in October 2020. We studied its X-ray properties during the first month of the outburst using XMM–Newton, NuSTAR and Swift observations. The shape and amplitude of the pulse profile varied significantly with energy. The broadband spectrum was well described using two absorbed blackbody components plus a faint power law component at high energies. Phase-resolved spectral analysis of the data suggests that the emission could be attributed to thermal photons from a single heated region with a complex shape on the star surface undergoing resonant Compton scattering on charged particles located in the magnetosphere. Modelling the evolutionary path of the magnetar with our magneto-thermal evolutionary codes indicates that SGR J1830 was born ≈23 kyr ago with a dipolar magnetic field of ∼1015 G, slightly larger than the current value.

We continue studying convection as a possible factor of episodic accretion in protoplanetary disks. Within the model of a viscous disk, the accretion history is analyzed at different rates and regions of matter inflow from the envelope onto the disk. It is shown that the burst-like regime occurs in a wide range of parameters. The long-term evolution of the disk is modeled, including the decreasing-with-time matter inflow from the envelope. It is demonstrated that the disk becomes convectively unstable and maintains burst-like accretion onto the star for several million years. The general conclusion of the study is that convection can serve as one of the mechanisms of episodic accretion in protostellar disks, but this conclusion needs to be verified using more consistent hydrodynamic models.

Due to the rich phenomenology and extreme magnetic conditions, magnetars will be targets of great interest for the upcoming polarimetry space missions. In particular, the Imaging X-ray Polarimetry Explorer (IXPE), recently launched in December 2021, will operate in the 2–8 keV range. This will open a new window to study the polarized, persistent X-ray emission from magnetars. In this talk, I will present simulations of IXPE observations of magnetars using the IXPEObsSim package. I will discuss future prospect to discriminate between different magnetar’s emission mechanisms, as well as a potential detection of the signal of vacuum birefringence using IXPE.

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.