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.

Gaia data allows for search for extended stellar structures in phase (coordinates plus velocities) space. We describe a method of using DBSCAN clustering algorithm, which is used to group closely-packed-together data points, to a list of preliminary selected pairs of stars, with parameters expected to be found within stellar streams and comoving groups: loose structures in which stars are not gravitationally bound, but do share motion and evolutionary properties. To test our approach, we construct a model population of background stars, and use pair-constructing and clustering algorithms on it. Results show that transitioning to a list of pairs sharply reveals structures not presented in background model, which then become more apparent targets in coordinate-velocity phase space for DBSCAN algorithm thanks to now increased relative density of the extended stellar structure.

The INTEGRAL satellite, in orbit since October 2002, has significantly contributed to the study of magnetars and, thanks to its unique capabilities for the study of transient gamma-ray phenomena, it is now playing an important role in multimessenger astrophysics. The most recent results include the discovery of a peculiar burst from SGR J1935+2154, which gave the first observational evidence for the connection between magnetars and fast radio bursts, and extensive searches for bursting activity in peculiar sources, such as the repeating FRB 20200120E in M81 and the ultra-long period magnetar candidate GLEAMX J162759.5–523504.3.

In this paper I present a method to enhance the search sensitivity for long transient Gravitational Waves produced by Neutron Star instabilities. This method consists in a selective image filter, called Triangular Filter, to be applied to data spectrograms. It is shown that thanks to this implementation a ∼20% gain in sensitivity is achievable.

In this paper we have inferred the magnetic shielding characteristics and space weather hazards of selected potentially habitable extrasolar planets using a dynamical geophysical model from calculations of internal heat, phases of volcanism and planetary magnetic moments. The space weather hazards on the extrasolar planet Kepler-452b orbiting around a Sun-like star are found to be a minimum which enhances the habitability probability of this planet.

We present Ekster, a new method for simulating the formation and dynamics of individual stars in a relatively low-resolution gas background. Here, we use Ekster to simulate star cluster formation in two different regions from each of two galaxy models with different spiral potentials. We simulate these regions for 3 Myr to study where and how star clusters form. We find that massive GMC regions form more massive clusters than sections of spiral arms. Additionally we find that clusters form both by accreting gas and by merging with other proto-clusters, the latter happening more frequently in the denser GMC regions.

The Laser Interferometer Space Antenna (LISA) mission will observe from space gravitational waves emitted by neutron stars and white dwarfs within galactic binaries. These compact stars can have intense magnetic fields. Therefore, the impact of the magnetic fields on the orbital and the spins evolution of binary systems can potentially be detected by LISA through the GW’s strain. Within the magnetic dipole-dipole approximation, we found that magnetism generates a secular drift of the mean longitude which, in turn, shifts all the frequencies contained in the GW signal. For a quasi-circular orbit, the signal is mainly monochromatic and the magnetic shift is proportional to the product of the magnetic moments and is inversely proportional to the 7/2 power of the semi-major axis. Hence, for a highly magnetic binary system in compact orbit, a non-negligible amount of the frequency measured by LISA might have a magnetic origin.

One of the models explaining the high luminosity of pulsing ultra-luminous X-ray sources (pULXs) was suggested by Mushtukov et al. (2015). They showed that the accretion columns on the surfaces of highly magnetized neutron stars can be very luminous due to opacity reduction in the high magnetic field. However, a strong magnetic field leads also to amplification of the electron-positron pairs creation. Therefore, increasing of the electron and positron number densities compensates the cross-section reduction, and the electron scattering opacity does not decrease with the magnetic field magnification. As a result, the maximum possible luminosity of the accretion column does not increase with the magnetic field. It ranges between 1040 − 1041 erg s−1 depending only slightly on the magnetic field strength.