The combination of strong magnetic fields and fast rotation is often invoked as a characteristic of the central engine for outstanding sources such as GRBs, hypernovae, and superluminous supernovae. However, the actual properties of the magnetic field during the collapse of the stellar progenitor are very uncertain, since they depend on the evolution of the star and can be affected by complex dynamo processes occurring in the central proto-neutron star. Using 3D relativistic MHD models we show that higher-order multipolar fields can lead to the onset of a supernova, although they tend to produce less energetic explosions and less collimated outflows. Quadrupolar fields efficiently extract angular momentum from the central core, but the rotational energy is partly stored in the equatorial regions, rather than powering up the polar outflows. Finally, our results show a strong magnetic quenching of the hydrodynamic non-axisymmetric instabilities that are associated to the emission of GWs.

We numerically investigate the gravitational collapse of rotating magnetic protostellar clouds. The simulations are performed using 2D MHD code ‘Enlil’. The code is based on TVD scheme of increased order of accuracy. We developed a model of the initially non-uniform cloud, which self-consistently treats gas density and large-scale magnetic field distribution. Simulation results for the typical parameters of a solar mass cloud are presented. In agreement with our previous results for the uniform cloud, the isothermal collapse of the non-uniform cloud results in formation of hierarchical structure of the cloud, consisting of flattened envelope and thin quasi-magnetostatic primary disk near its equatorial plane. The non-uniform cloud collapses longer than the uniform one, since the magnetic field is dynamically stronger at the periphery of the cloud in the former case.

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.

We present our findings on magnetar spectral line analysis in the context of upcoming high resolution, high effective area, high throughput X-ray telescopes for two cases: persistent magnetar emission and magnetar bursts. For magnetars in quiescence, we present our preliminary work on modelling for phase-resolved emission. Our results reveal the necessity of constraining line depth and width concurrently with line energy to conclusively determine hotspot emission and corresponding magnetic field geometry. We then present the results of our simulations using effective area and response of various current and upcoming X-ray telescopes for magnetar spectral line detection and expand on the exciting opportunities upcoming telescopes provide to probe quiescent and burst emission region geometry and propagation in the extreme magnetic field of a magnetar.

A numerical nonstationary two-dimensional MHD model of a protoplanetary disk near T Tauri star with a jet is developed. The model assumes consideration of pure gas ionization, optically thin cooling, anisotropic thermal conductivity and viscosity. The relaxation of gas-dynamic flows is analyzed. Based on the evolution of plasma flows, profiles of hydrogen spectral lines have been obtained.

PRESTALINE is a package allowing a user to simulate and analyse spectra of various astrophysical objects. The package is based on the numerical models PRESTA (Kochina & Wiebe (2017)) and RADEX (van der Tak et al. (2007)). PRESTALINE provides the direct comparison of theoretical models with observations and allows estimating physical conditions in a studied object, such as kinetic temperature and chemical composition. Here we present the results of applying PRESTALINE to the test object DR21(OH) and discuss possible applications and future extensions of the project.

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.

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.

The population of black widows, binary systems containing a millisecond pulsar and a very low-mass companion star exposed to the high-energy pulsar wind, has grown exponentially in the past few years. The number of black widow candidates is now over 30 systems, but only 14 have been confirmed so far. Their relevance in analysing the extremes of the neutron stars properties led to multiwavelength dedicated studies that revealed a rich phenomenology. In this work, we provide a glimpse into the black widow class through modelling of high-cadence multi-band light curves of 6 systems, accounting for almost half of the confirmed population. A better understanding of the black widow population, which hosts some of the most massive and fastest spinning neutron stars, will ultimately benefit future modelling of compact object mergers.

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.

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.

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.

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.

We investigate the broad-band behaviour of circular polarization in radio pulsar profiles and show the relationship between polarization fraction and what proportion of that polarization is circular, both across frequency, and for a large number of pulsars viewed collectively. The behaviour observed may be explained by pulsar polarization originating from the partially-coherent combination of two linearly-polarized orthogonal modes with different flux spectral indices. (See also the poster in the “supplementary information”.)

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.

The Sun moves with respect to the local interstellar medium (LISM) and modifies its properties to heliocentric distances as large as 1 pc. The solar wind (SW) is affected by penetration of the LISM neutral particles, especially H and He atoms. Charge exchange between the LISM atoms and SW ions creates pickup ions (PUIs) and secondary neutral atoms that can propagate deep into the LISM. Neutral atoms measured at 1 au can provide us with valuable information on the properties of pristine LISM. Voyager 1 and 2 spacecraft perform in situ measurements of the LISM perturbed by the presence of the heliosphere and relate them to the unperturbed region. We discuss observational data and numerical simulations that shed light onto the mutual influence of the SW and LISM. Physical phenomena accompanying the SW–LISM interaction are discussed, including the coupling of the heliospheric and interstellar magnetic field at the heliopause.

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.

In this poster we present the structure of an axisymmetric, force-free magnetosphere of a twisted, aligned rotating dipole within a corotating plasma-filled magnetosphere. We explore various profiles for the twist. We find that as the current increases more field lines cross the light cylinder leading to more efficient spin-down. Moreover, we notice that the twist cannot be increased indefinitely and after a finite twist of about π/2 the field becomes approximately radial. This could have implications for torque variations of magnetars related to outbursts.

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.

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.