On the strength of a double Nobel prize winning experiment (Super)Kamiokande and an extremely successful long baseline neutrino programme, the third generation Water Cherenkov detector, Hyper-Kamiokande, is being developed by an international collaboration as a leading worldwide experiment based in Japan. The Hyper-Kamiokande detector will be hosted in the Tochibora mine, about 295 km away from the J-PARC proton accelerator research complex in Tokai, Japan. The currently existing accelerator will be steadily upgraded to reach a MW beam by the start of the experiment. A suite of near detectors will be vital to constrain the beam for neutrino oscillation measurements. A new cavern will be excavated at the Tochibora mine to host the detector. The experiment will be the largest underground water Cherenkov detector in the world and will be instrumented with new technology photosensors, faster and with higher quantum efficiency than the ones in Super-Kamiokande. The science that will be developed will be able to shape the future theoretical framework and generations of experiments. Hyper-Kamiokande will be able to measure with the highest precision the leptonic CP violation that could explain the baryon asymmetry in the Universe. The experiment also has a demonstrated excellent capability to search for proton decay, providing a significant improvement in discovery sensitivity over current searches for the proton lifetime. The atmospheric neutrinos will allow to determine the neutrino mass ordering and, together with the beam, able to precisely test the three-flavour neutrino oscillation paradigm and search for new phenomena. A strong astrophysical programme will be carried out at the experiment that will detect supernova neutrinos and will measure precisely solar neutrino oscillation.

Software that processes real-world data or that models a physical system must have some way of managing units. While simple approaches like the understood convention that all data are in a unit system (such as the MKS SI unit system) do work in practice, they are fraught with possible sources of error both by developers and users of the software. In this paper we present unyt, a Python library based on NumPy and SymPy for handling data that has units. It is designed both to aid quick interactive calculations and to be tightly integrated into a larger Python application or library. We compare unyt with two other Python libraries for handling units, Pint and astropy.units, and find that unyt is faster, has higher test coverage, and has fewer lines of code.

We present the latest results of our search for planets with HARPS-N at the 3.6 m Telescopio Nazionale Galileo under the Tracking Advanced Planetary Systems project: an in-depth study of the 15 most Li abundant giants from the PennState - Toru\'n Planet Search sample. Our goals are first, to obtain radial velocities of the most Li-rich giants we identified in our sample to search for possible low-mass substellar companions, and second, to perform an extended spectral analysis to define the evolutionary status of these stars. Methods. This work is based on high-resolution spectra obtained with the Hobby-Eberly Telescope and its High Resolution Spectro- graph, and with the HARPS-N spectrograph at the Telescopio Nazionale Galileo. Two stars, HD 181368 and HD 188214 , were also observed with UVES at the VLT to determine beryllium abundances. We report i) the discovery of two new planetary systems around the Li-rich giant stars: HD 238914 and TYC 3318-01333- 1 (a binary system); ii) reveal a binary Li-rich giant, HD 181368 ; iii) although our current phase coverage is not complete, we suggest the presence of planetary mass companions around TYC 3663-01966-1 and TYC 3105-00152-1 ; iv) we confirm the previous result for BD+48 740 and present updated orbital parameters, and v) we find a lack of a relation between the Li enhancement and the Be abundance for the stars HD 181368 and HD 188214 , for which we acquired blue spectra. We found seven stars with stellar or potential planetary companions among the 15 Li-rich giant stars. The binary star frequency of the Li-rich giants in our sample appears to be normal, but the planet frequency is twice that of the general sample, which suggests a possible connection between hosting a companion and enhanced Li abundance in giant stars. We also found most of the companions orbits to be highly eccentric.

We present the "Galaxy Cluster Merger Catalog." This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers and clusters from cosmological simulations. These mock observations consist of projections of a number of important observable quantities in several different wavebands as well as along different lines of sight through each simulation domain. The web interface to the catalog consists of easily browseable images over epoch and projection direction, as well as download links for the raw data and a JS9 interface for interactive data exploration. The data is presented within a consistent format so that comparison between simulations is straightforward. All of the data products are provided in the standard FITS file format. Data is being stored on the yt Hub (http://hub.yt), which allows for remote access and analysis using a Jupyter notebook server. Future versions of the catalog will include simulations from a number of research groups and a variety of research topics related to the study of interactions of galaxy clusters with each other and with their member galaxies. The catalog is located at http://gcmc.hub.yt.

We present two evolved stars from the TAPAS (Tracking Advanced PlAnetary Systems) with HARPS-N project devoted to RV precision measurements of identified candidates within the PennState - Torun Centre for Astronomy Planet Search. Evolved stars with planets are crucial to understand the dependency of the planet formation mechanism on the mass and metallicity of the parent star and to study star-planet interactions. The paper is based on precise radial velocity (RV) measurements, for HD 103485 we collected 57 epochs over 3317 days with the Hobby-Eberly Telescope and its High Resolution Spectrograph and 18 ultra-precise HARPS-N data over 919 days. For BD+03 2562 we collected 46 epochs of HET data over 3380 days and 19 epochs of HARPS-N data over 919 days. We present the analysis of the data and the search for correlations between the RV signal and stellar activity, stellar rotation and photometric variability. Based on the available data, we interpret the RV variations measured in both stars as Keplerian motion. Both stars have masses close to Solar (1.11 and 1.14), very low metallicities ([Fe/H]=-0.50 and -0.71), and, both have Jupiter planetary mass companions (m sin i=7 and 6.4 Mj), in close to terrestrial orbits (1.4 and 1.3~au), with moderate eccentricities (e=0.34 and 0.2). However, we cannot totally exclude that the signal in the case of HD~103485 is due to rotational modulation of active regions. Based on the current data, we conclude that BD+03 2562 has a bona fide planetary companion while for HD 103485 we cannot totally exclude that the best explanation for the RV signal modulations is not the existence of a planet but stellar activity. If, the interpretation remains that both stars have planetary companions they represent systems orbiting very evolved stars with very low metallicities, a challenge to the conditions required for the formation of massive giant gas planets.

We present measurements of $\rho^0$, $\omega$ and K$^{*0}$ spectra in $\pi^{-} + $C production interactions at 158 GeV/c and $\rho^0$ spectra at 350 GeV/c using the NA61/SHINE spectrometer at the CERN SPS. Spectra are presented as a function of the Feynman's variable $x_\text{F}$ in the range $0 < x_\text{F} < 1$ and $0 < x_\text{F} < 0.5$ for 158 GeV/c and 350 GeV/c respectively. Furthermore, we show comparisons with previous measurements and predictions of several hadronic interaction models. These measurements are essential for a better understanding of hadronic shower development and for improving the modeling of cosmic ray air showers.

We present the latest result of the TAPAS project that is devoted to intense monitoring of planetary candidates that are identified within the PennState-Toru\'n planet search. We aim to detect planetary systems around evolved stars to be able to build sound statistics on the frequency and intrinsic nature of these systems, and to deliver in-depth studies of selected planetary systems with evidence of star-planet interaction processes. The paper is based on precise radial velocity measurements: 13 epochs collected over 1920 days with the Hobby-Eberly Telescope and its High-Resolution Spectrograph, and 22 epochs of ultra-precise HARPS-N data collected over 961 days. We present a warm-Jupiter ($T_{eq}=1350 K$, $m_{2} sin i=5.4\pm$0.4$M_{J}$) companion with an orbital period of 26.468 days in a circular ($e=0.036$) orbit around a giant evolved ($\log g=3.11\pm0.09$, $R=6.26\pm0.86R_{\odot}$) star with $M_{\star}=1.87\pm0.17M_{\odot}$. This is the most massive and oldest star found to be hosting a close-in giant planet. Its proximity to its host ($a=0.21au$) means that the planet has a $13.9\pm2.0\%$ probability of transits; this calls for photometric follow-up study. This massive warm Jupiter with a near circular orbit around an evolved massive star can help set constraints on general migration mechanisms for warm Jupiters and, given its high equilibrium temperature, can help test energy deposition models in hot Jupiters.

Evolved stars are crucial pieces to understand the dependency of the planet formation mechanism on the stellar mass and to explore deeper the mechanism involved in star-planet interactions. Over the past ten years, we have monitored about 1000 evolved stars for radial velocity variations in search for low-mass companions under the Penn State - Torun Centre for Astronomy Planet Search program with the Hobby-Eberly Telescope. Selected prospective candidates that required higher RV precision measurements have been followed with HARPS-N at the 3.6 m Telescopio Nazionale Galileo under the TAPAS project. We aim to detect planetary systems around evolved stars to be able to build sound statistics on the frequency and intrinsic nature of these systems, and to deliver in-depth studies of selected planetary systems with evidence of star-planet interaction processes. For HD 5583 we obtained 14 epochs of precise RV measurements collected over 2313 days with the Hobby-Eberly Telescope (HET), and 22 epochs of ultra-precise HARPS-N data collected over 976 days. For BD+15 2375 we collected 24 epochs of HET data over 3286 days and 25 epochs of HARPS-S data over 902 days. We report the discovery of two planetary mass objects orbiting two evolved Red Giant stars: HD~5583 has a m sin i = 5.78 M$_{J}$ companion at 0.529~AU in a nearly circular orbit (e=0.076), the closest companion to a giant star detected with the RV technique, and BD+15~2735 that with a m sin i= 1.06 M$_{J}$ holds the record of the lightest planet found so far orbiting an evolved star (in a circular e=0.001, 0.576~AU orbit). These are the third and fourth planets found within the TAPAS project, a HARPS-N monitoring of evolved planetary systems identified with the Hobby-Eberly Telescope.

Lithium rich giant stars are rare objects. For some of them, Li enrichment exceeds abundance of this element found in solar system meteorites, suggesting that these stars have gone through a Li enhancement process. We identified a Li rich giant HD 107028 with A(Li) > 3.3 in a sample of evolved stars observed within the PennState Torun Planet Search. In this work we study different enhancement scenarios and we try to identify the one responsible for Li enrichment for HD 107028. We collected high resolution spectra with three different instruments, covering different spectral ranges. We determine stellar parameters and abundances of selected elements with both equivalent width measurements and analysis, and spectral synthesis. We also collected multi epoch high precision radial velocities in an attempt to detect a companion. Collected data show that HD 107028 is a star at the base of Red Giant Branch. Except for high Li abundance, we have not identified any other anomalies in its chemical composition, and there is no indication of a low mass or stellar companion. We exclude Li production at the Luminosity Function Bump on RGB, as the effective temperature and luminosity suggest that the evolutionary state is much earlier than RGB Bump. We also cannot confirm the Li enhancement by contamination, as we do not observe any anomalies that are associated with this scenario. After evaluating various scenarios of Li enhancement we conclude that the Li-overabundance of HD 107028 originates from Main Sequence evolution, and may be caused by diffusion process.

We present three giant stars from the ongoing Penn State-Toru\'n Planet Search with the Hobby-Eberly Telescope, which exhibit radial velocity variations that point to a presence of planetary --mass companions around them. BD+49 828 is a $M=1.52 \pm 0.22$ $M_{\odot}$ K0 giant with a $m sini$=$1.6^{+0.4}_{-0.2}$ $M_{J}$ minimum mass companion in $a=4.2^{+0.32}_{-0.2}$ AU ($2590^{+300}_{-180}$d), $e=0.35^{+0.24}_{-0.10}$ orbit. HD 95127, a log$L$/$L_{\odot}$=$2.28 \pm 0.38$, $R = 20\pm 9$ $R_{\odot}$, $M=1.20 \pm 0.22$ $M_{\odot}$ K0 giant has a $m sini$=$5.01^{+0.61}_{-0.44}$ $M_{J}$ minimum mass companion in $a=1.28^{+0.01}_{-0.01}$ AU ($482^{+5}_{-5}$d), $e=0.11^{+0.15}_{-0.06}$ orbit. Finally, HD 216536, is a $M=1.36 \pm 0.38$ $M_{\odot}$ K0 giant with a $m sin i=1.47^{+0.20}_{-0.12}$ $M_{J}$ minimum mass companion in $a=0.609^{+0.002}_{-0.002}$ AU ($148.6^{+0.7}_{-0.7}$d), $e=0.38^{+0.12}_{-0.10}$ orbit. Both, HD 95127 b and HD 216536 b in their compact orbits, are very close to the engulfment zone and hence prone to ingestion in the near future. BD+49 828 b is among the longest period planets detected with the radial velocity technique until now and it will remain unaffected by stellar evolution up to a very late stage of its host. We discuss general properties of planetary systems around evolved stars and planet survivability using existing data on exoplanets in more detail.

Context. Stars that have evolved-off the Main Sequence are crucial in expanding the frontiers of knowledge on exoplanets toward higher stellar masses, and to constrain star-planet interaction mechanisms. These stars, however suffer from intrinsic activity that complicates the interpretation of precise radial velocity measurement and are often avoided in planet searches. We have, over the last 10 years, monitored about 1000 evolved stars for radial velocity variations in search for low-mass companions under the Penn State - Toru\'n Centre for Astronomy Planet Search with the Hobby-Eberly Telescope. Selected prospective candidates that required higher RV precision meassurements have been followed with HARPS-N at the 3.6 m Telescopio Nazionale Galileo. Aims. To detect planetary systems around evolved stars, to be able to build sound statistics on the frequency and intrinsic nature of these systems, and to deliver in-depth studies of selected planetary systems with evidences of star-planet interaction processes. Methods. We have obtained for TYC 1422-614-1 69 epochs of precise radial velocity measurements collected over 3651 days with the Hobby-Eberly Telescope, and 17 epochs of ultra precise HARPS-N data collected over 408 days. We have complemented these RV data with photometric time-series from the All Sky Automatic Survey archive. Results. We report the discovery of a multiple planetary system around the evolved K2 giant star TYC 1422-614-1. The system orbiting the 1.15 M$_\odot$ star is composed of a planet with mass m$sin i$=2.5 M$_J$ in a 0.69 AU orbit, and a planet/brown dwarf with m$sin i$=10 M$_J$ in a 1.37 AU orbit. The multiple planetary system orbiting TYC 1422-614-1 is the first finding of the TAPAS project, a HARPS-N monitoring of evolved planetary systems identified with the Hobby-Eberly Telescope.

We present simulations of the magnetized interstellar medium (ISM) in models of massive star forming (40 Msun / yr) disk galaxies with high gas surface densities (~100 Msun / pc^2) similar to observed star forming high-redshift disks. We assume that type II supernovae deposit 10 per cent of their energy into the ISM as cosmic rays and neglect the additional deposition of thermal energy or momentum. With a typical Galactic diffusion coefficient for CRs (3e28 cm^2 / s) we demonstrate that this process alone can trigger the local formation of a strong low density galactic wind maintaining vertically open field lines. Driven by the additional pressure gradient of the relativistic fluid the wind speed can exceed 1000 km/s, much higher than the escape velocity of the galaxy. The global mass loading, i.e. the ratio of the gas mass leaving the galactic disk in a wind to the star formation rate becomes of order unity once the system has settled into an equilibrium. We conclude that relativistic particles accelerated in supernova remnants alone provide a natural and efficient mechanism to trigger winds similar to observed mass-loaded galactic winds in high-redshift galaxies. These winds also help explaining the low efficiencies for the conversion of gas into stars in galaxies as well as the early enrichment of the intergalactic medium with metals. This mechanism can be at least of similar importance than the traditionally considered momentum feedback from massive stars and thermal and kinetic feedback from supernova explosions.

We investigate streaming instability using two-fluid approximation (neutral gas and dust) in a quasi-global, unstratified protoplanetary disc, with the help of PIERNIK code. We compare amplification rate of the eigen-mode in numerical simulations, with the corresponding growth resulting from the linear stability analysis of full system of Euler's equation including aerodynamic drag. Following Youdin & Goodman (2005) we show that (1) rapid dust clumping occurs due to the difference in azimuthal velocities of gas and dust, coupled by the drag force, (2) initial density perturbations are amplified by several orders of magnitude. We demonstrate that the multi-fluid extension of the simple and efficient Relaxing TVD scheme, implemented in PIERNIK, leads to results, which are compatible with those obtained with other methods.

We present three-dimensional global numerical simulations of the cosmic-ray (CR) driven dynamo in barred galaxies. We study the evolution of the interstellar medium of the barred galaxy in the presence of non-axisymmetric component of the potential, i.e., the bar. The magnetohydrodynamical dynamo is driven by CRs, which are continuously supplied to the disk by supernova (SN) remnants. No magnetic field is present at the beginning of simulations but one-tenth of SN explosions is a source of a small-scale randomly oriented dipolar magnetic field. In all models we assume that 10% of 10^51 erg SN kinetic energy output is converted into CR energy. To compare our results directly with the observed properties of galaxies, we construct realistic maps of polarized radio emission. The main result is that the CR-driven dynamo can amplify weak magnetic fields up to a few {\mu}G within a few Gyr in barred galaxies. The obtained e-folding time is equal to 300 Myr and the magnetic field reaches equipartition at time t ~ 4.0 Gyr. Initially, the completely random magnetic field evolves into large-scale structures. An even (quadrupole-type) configuration of the magnetic field with respect to the galactic plane can be observed. Additionally, the modeled magnetic field configuration resembles maps of the polarized intensity observed in barred galaxies. Polarization vectors are distributed along the bar and between spiral arms. Moreover, the drift of magnetic arms with respect to the spiral pattern in the gas density distribution is observed during the entire simulation time.

We present recent developments of global galactic-scale numerical models of the Cosmic Ray (CR) driven dynamo, which was originally proposed by Parker (1992). We conduct a series of direct CR+MHD numerical simulations of the dynamics of the interstellar medium (ISM), composed of gas, magnetic fields and CR components. We take into account CRs accelerated in randomly distributed supernova (SN) remnants, and assume that SNe deposit small-scale, randomly oriented, dipolar magnetic fields into the ISM. The amplification timescale of the large-scale magnetic field resulting from the CR-driven dynamo is comparable to the galactic rotation period. The process efficiently converts small-scale magnetic fields of SN-remnants into galactic-scale magnetic fields. The resulting magnetic field structure resembles the X-shaped magnetic fields observed in edge-on galaxies.

We study the global evolution of the magnetic field and interstellar medium (ISM) of the barred and ringed galaxies in the presence of non-axisymmetric components of the potential, i.e. the bar and/or the oval perturbations. The magnetohydrodynamical dynamo is driven by cosmic rays (CR), which are continuously supplied to the disk by supernova (SN) remnants. Additionally, weak, dipolar and randomly oriented magnetic field is injected to the galactic disk during SN explosions. To compare our results directly with the observed properties of galaxies we construct realistic maps of high-frequency polarized radio emission. The main result is that CR driven dynamo can amplify weak magnetic fields up to few $\mu$G within few Gyr in barred and ringed galaxies. What is more, the modelled magnetic field configuration resembles maps of the polarized intensity observed in barred and ringed galaxies.

We present a new particle module of the magnetohydrodynamic (MHD) Piernik code. The original multi-fluid grid code based on the Relaxing Total Variation Diminishing (RTVD) scheme has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. In this poster we introduce the scheme used to solve equations of motion for the particles and present the first results coming from the module. The results of test problems are also compared with the results coming from fluid simulations made with Piernik-MHD code. The comparison shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.

Piernik is a multi-fluid grid magnetohydrodynamic (MHD) code based on the Relaxing Total Variation Diminishing (RTVD) conservative scheme. The original code has been extended by addition of dust described within the particle approximation. The dust is now described as a system of interacting particles. The particles can interact with gas, which is described as a fluid. The comparison between the test problem results and the results coming from fluid simulations made with Piernik code shows the most important differences between fluid and particle approximations used to describe dynamical evolution of dust under astrophysical conditions.

We report first results of first global galactic-scale CR-MHD simulations of cosmic-ray-driven dynamo. We investigate the dynamics of magnetized interstellar medium (ISM), which is dynamically coupled with the cosmic-ray (CR) gas. We assume that exploding stars deposit small-scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. We conduct numerical simulations with the aid of a new parallel MHD code PIERNIK. We find that the initial magnetization of galactic disks by exploding magnetized stars forms a favourable conditions for the cosmic-ray-driven dynamo. We demonstrate that dipolar magnetic fields supplied on small SN-remnant scales, can be amplified exponentially, by the CR-driven dynamo, to the present equipartition values, and transformed simultaneously to large galactic-scales. The resulting magnetic field structure in an evolved galaxy appears spiral in the face-on view and reveals the so called X-shaped structure in the edge-on view.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present an extension of PIERNIK, which is designed for simulations of diffusive propagation of the Cosmic-Ray (CR) component in the magnetized ISM.

We present new developments on the Cosmic--Ray driven, galactic dynamo, modeled by means of direct, resistive CR--MHD simulations, performed with ZEUS and PIERNIK codes. The dynamo action, leading to the amplification of large--scale galactic magnetic fields on galactic rotation timescales, appears as a result of galactic differential rotation, buoyancy of the cosmic ray component and resistive dissipation of small--scale turbulent magnetic fields. Our new results include demonstration of the global--galactic dynamo action driven by Cosmic Rays supplied in supernova remnants. An essential outcome of the new series of global galactic dynamo models is the equipartition of the gas turbulent energy with magnetic field energy and cosmic ray energy, in saturated states of the dynamo on large galactic scales.

We conduct global galactic--scale magnetohydrodynamical (MHD) simulations of the cosmic--ray driven dynamo. We assume that exploding stars deposit small--scale, randomly oriented, dipolar magnetic fields into the differentially rotating ISM, together with a portion of cosmic rays, accelerated in supernova shocks. Our simulations are performed with the aid of a new parallel MHD code PIERNIK. We demonstrate that dipolar magnetic fields supplied on small SN--remnant scales, can be amplified exponentially by the CR--driven dynamo to the present equipartition values, and transformed simultaneously to large galactic--scales by an inverse cascade promoted by resistive processes.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin & Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we present Ohmic resistivity extension of the original Relaxing TVD MHD scheme, and show examples of magnetic reconnection in cases of uniform and current-dependent resistivity prescriptions.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme (Jin and Xin, 1995). The original scheme (see Trac & Pen (2003) and Pen et al. (2003)) has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we introduce the multifluid extension of Relaxing TVD scheme and present a test case of dust migration in a two-fluid disk composed of gas and dust. We demonstrate that due to the difference in azimuthal velocities of gas and dust and the drag force acting on both components dust drifts towards maxima of gas pressure distribution.

We present a new multi-fluid, grid MHD code PIERNIK, which is based on the Relaxing TVD scheme. The original scheme has been extended by an addition of dynamically independent, but interacting fluids: dust and a diffusive cosmic ray gas, described within the fluid approximation, with an option to add other fluids in an easy way. The code has been equipped with shearing-box boundary conditions, and a selfgravity module, Ohmic resistivity module, as well as other facilities which are useful in astrophysical fluid-dynamical simulations. The code is parallelized by means of the MPI library. In this paper we shortly introduce basic elements of the Relaxing TVD MHD algorithm, following Trac & Pen (2003) and Pen et al. (2003), and then focus on the conservative implementation of the shearing box model, constructed with the aid of the Masset's (2000) method. We present results of a test example of a formation of a gravitationally bounded object (planet) in a self-gravitating and differentially rotating fluid.

We conduct a series of magnetohydrodynamical (MHD) simulations of magnetized interstellar medium (ISM) disturbed by exploding stars. Each star deposits a randomly oriented, dipolar magnetic field into ISM. The simulations are performed in a Cartesian box, in a reference frame that is corotating with the galactic disk. The medium is stratified by vertical galactic gravity. The resulting turbulent state of ISM magnetized by the stellar explosions is processed with the aid of Fourier analysis. The results leads to the conclusion that the input of magnetic energy from exploding stars is additionally multiplied by differential rotation. The resulting magnetic field appears to grow up in small-scale component, while the total magnetic flux remains limited. Our results indicate that magnetic field originating from exploding stars can be a source of initial magnetic fields for a subsequent dynamo process.