Equilibrium configurations of the internal magnetic field of a pulsar play a key role in modelling astrophysical phenomena from glitches to gravitational wave emission. In this paper, we present a numerical scheme for solving the Grad–Shafranov equation and calculating equilibrium configurations of pulsars, accounting for superconductivity in the core of the neutron star, and for the Hall effect in the crust of the star. Our numerical code uses a finite difference method in which the source term appearing in the Grad–Shafranov equation, which is used to model the magnetic equilibrium is non-linear. We obtain solutions by linearising the source and applying an under-relaxation scheme at each step of computation to improve the solver’s convergence. We have developed our code in both C++ and Python, and our numerical algorithm can further be adapted to solve any non-linear PDEs appearing in other areas of computational astrophysics. We produce mixed toroidal–poloidal field configurations, and extend the portion of parameter space that can be investigated with respect to previous studies. We find that in even in the more extreme cases, the magnetic energy in the toroidal component does not exceed approximately 5% of the total. We also find that if the core of the star is superconducting, the toroidal component is entirely confined to the crust of the star, which has important implications for pulsar glitch models which rely on the presence of a strong toroidal field region in the core of the star, where superfluid vortices pin to superconducting fluxtubes.

In this study, a peridynamic material model for a polycrystalline ice is utilised to investigate its fracture behaviour under dynamic loading condition. First, the material model was validated by considering a single grain, double grains and polycrystalline structure under tension loading condition. Peridynamic results are compared against finite element analysis results without allowing failure. After validating the material model, dynamic analysis of a polycrystalline ice material with two pre-existing cracks under tension loading is performed by considering weak and strong grain boundaries with respect to grain interiors. Numerical results show that the effect of microstructure is significant for weak grain boundaries. On the other hand, for strong grain boundaries, the effect of microstructure is insignificant. The evaluated results have demonstrated that peridynamics can be a very good alternative numerical tool for fracture analysis of polycrystalline ice material.

The Laplace resonance is a configuration that involves the commensurability between the mean motions of three small bodies revolving around a massive central one. This resonance was first observed in the case of the three inner Galilean satellites, Io, Europa, and Ganymede. In this work the Laplace resonance is generalised by considering a system of three satellites orbiting a planet that are involved in mean motion resonances. These Laplace-like resonances are classified in three categories: first-order (2:1&2:1, 3:2&3:2, 2:1&3:2), second-order (3:1&3:1) and mixed-order resonances (2:1&3:1). In order to study the dynamics of the system we implement a model that includes the gravitational interaction with the central body, the mutual gravitational interactions of the satellites, the effects due to the oblateness of the central body and the secular interaction of a fourth satellite and a distant star. Along with these contributions we include the tidal interaction between the central body and the innermost satellite. We study the survival of the Laplace-like resonances and the evolution of the orbital elements of the satellites under the tidal effects. Moreover, we study the possibility of capture into resonance of the fourth satellite.

. We present a significantly improved version of our numerical code JASMINE, that can now solve the Jeans equations for axisymmetric models of stellar systems, composed of an arbitrary number of stellar populations, a Dark Matter halo, and a central Black Hole. The stellar components can have different structural (density profile, flattening, mass, scale length), dynamical (rotational support, velocity dispersion anisotropy), and population (age, metallicity, Initial Mass Function, mass-to-light ratio) properties. These models, when combined with observations, will allow to investigate important issues, such as quantifying the systematic effects of IMF variations, of mass-to-light ratio gradients, and of different stellar kinematic components (e.g. counter rotating disks, kinematically decoupled cores) on luminosity-weighted properties. The developed analytical and numerical framework aims at modeling Early-Type Galaxies, but it can also be applied to dwarf Spheroidal galaxies and Globular Clusters.

beans is a tool for distributed data analysis. It provides web and command line interface for data analysis and plotting for huge datasets. beans is written in a general form and can be used in any field of research to analyze the data. The main purpose of beans is to provide to the community a versatile tool to store, analyze and then visualize any amount of scientific data (e.g. numerical simulation, observations).

We present in this work the development of a solar data assimilation method based on an axisymmetric mean field dynamo model and magnetic surface data. Our mid-term goal is to predict the solar quasi cyclic activity. We focus on the ability of our variational data assimilation algorithm to constrain the deep meridional circulation of the Sun based on solar magnetic observations. Within a given assimilation window, the assimilation procedure minimizes the differences between data and the forecast from the model, by finding an optimal meridional circulation in the convection zone, and an optimal initial magnetic field, via a quasi-Newton algorithm. We demonstrate the capability of the technique to estimate the meridional flow by a closed-loop experiment involving 40 years of synthetic, solar-like data. We show that the method is robust in estimating a (stochastic) time-varying flow fluctuating 30% about the average, and that the horizon of predictability of the method is ~ 1 cycle length.

In this article we described CoLiTec software for full automated frames processing. CoLiTec software allows processing the Big Data of observation results as well as processing of data that is continuously formed during observation. The scope of solving tasks includes frames brightness equalization, moving objects detection, astrometry, photometry, etc. Along with the high efficiency of Big Data processing CoLiTec software also ensures high accuracy of data measurements. A comparative analysis of the functional characteristics and positional accuracy was performed between CoLiTec and Astrometrica software. The benefits of CoLiTec used with wide field and low quality frames were observed. The efficiency of the CoLiTec software was proved by about 700.000 observations and over 1.500 preliminary discoveries.

We describe the implementation and performance results of our massively parallel MPI†/OpenMP‡ hybrid TreePM code for large-scale cosmological N-body simulations. For domain decomposition, a recursive multi-section algorithm is used and the size of domains are automatically set so that the total calculation time is the same for all processes. We developed a highly-tuned gravity kernel for short-range forces, and a novel communication algorithm for long-range forces. For two trillion particles benchmark simulation, the average performance on the fullsystem of K computer (82,944 nodes, the total number of core is 663,552) is 5.8 Pflops, which corresponds to 55% of the peak speed.

An estimation of the mass of the main asteroid belt was made on the basis of the new version of EPM2014 ephemerides of the Institute of Applied Astronomy of Russian Academy of Sciences using about 800000 positional observations of planets and spacecraft. We obtained the individual estimations of masses of large asteroids from radar data, as well as estimates of the masses of asteroids by using known diameters and estimated average densities for the three taxonomic types (C, S, M), and used the known mass values of binary asteroids and asteroids to which spacecraft approached. A two-dimensional homogeneous annulus with dimensions corresponding observed width of the main asteroid belt (2.06 au and 3.27 au) was used instead of a previous massive one-dimensional ring for modeling total perturbations from small asteroids. The obtained value of the total mass of the main asteroid belt is (12.25 ± 0.19)10−10M⊙.

We use the Magneticum Pathfinder (www.magneticum.org) hydro-dynamical cosmological simulation set to investigate the buildup of the stellar component within cosmological structures. These simulations result in the self-consistent formation of ICM, AGNs, and both spheroidal and disk galaxy populations, which properly reproduce the observed properties.

The Magneticum Pathfinder (www.magneticum.org) cosmological, hydro-dynamical simulation (896h-1Mpc)3 follows in detail the thermal and chemical evolution of the ICM as well as the evolution of SMBHs and their associated feedback processes. We demonstrate that assuming cosmological parameters inferred from the CMB, the thermal SZ power spectrum as observed by PLANCK is well matched by the deep light-cones constructed from these cosmological simulations. The thermal SZ prediction from the full SZ maps are significantly exceeding previous templates at large l (e.g., l > 1000) and therefore predict a significantly larger contribution to the signal at l = 3000 compared to previous findings. The excess of positive values within the probability distribution of the thermal SZ signal within the simulated light-cone agrees with the one seen by PLANCK. This excess signal follows a power law shape with an index of roughly -3.2. The bulk of the thermal SZ signal originates from clusters and groups which form between z = 0 and z ≈ 2 where at high redshift (z > 1) significant part of the signal originates from proto-cluster regions, which are not yet virialized. The simulation predicts a mean fluctuating Compton Y value of 1.18 × 10-6, with a remaining contribution of almost 5 ×10-7 when removing contribution from halos above a virial mass of 1013 M⊙/h.

The smallest dark matter halos are formed first in the early universe. According to recent studies, the central density cusp is much steeper in these halos than in larger halos and scales as ρ ∝ r−(1.5−1.3). We present results of very large cosmological N-body simulations of the hierarchical formation and evolution of halos over a wide mass range, beginning from the formation of the smallest halos. We confirmed early studies that the inner density cusps are steeper in halos at the free streaming scale. The cusp slope gradually becomes shallower as the halo mass increases. The slope of halos 50 times more massive than the smallest halo is approximately -1.3. The concentration parameter is nearly independent of halo mass, and ruling out simple power law mass-concentration relations. The steeper inner cusps of halos near the free streaming scale enhance the annihilation luminosity of a Milky Way sized halo between 12 to 67%.

Using a set of simulations in the ΛCDM cosmology, and the mass of intermediate X-Ray Active Galactic Nuclei (AGN) dark matter halos (DMH), we study the Halo Occupation Distribution (HOD) of their host DMH and compare with recent observations. We assume that intermediate X-Ray AGNs have been triggered by a major merger of sub-halos inside larger host halos. We find that the binary major merger HOD slope (α ≈ 0.8) do not seem to reproduce the observed HOD slope (α < 0.6) for these type of AGNs. Other mechanisms may be igniting these AGNs.

The growing number of detected planets in binary star systems requires methods for a quick assessment of possible habitability of planets in such environments. We offer an analytic method to determine habitability of a terrestrial planet in binary star systems. In this context we give an answer to the most important question: Do the radiative and gravitational perturbations of a secondary influence the extent of the habitable zone (HZ)?

After we have defined the borders of the HZ, we will show the dynamical behaviour of a terrestrial planet in the HZ when adding a Jupiter to the system. In such a system the HZ shows signs of mean motion resonances and secular resonances, depending on the architecture of the planetary system.

Radial Velocity (RV) plots of emission lines, including Hα, from non-magnetic Cataclysmic Variables (CV), are usually fit with a sin curve. This sin fit sometimes does not prove to be the best fit for some non-magnetic CVs with accretion discs. An analytical model is created based on our 3D Smoothed Particle Hydrodynamic (SPH) numerical code in order to simulate the RV curves. The observational target is WW And, a long period non-magnetic CV. The model takes into account disc ellipticity and inclination angle that provide good non-sinusoidal fits to the observed RV data.

In this paper, an analytical approach is proposed to model chip formation in the case ofturning process. Numerical simulations of chip genesis are performed in order to fitefficiently the proposed analy- tical model. In particular, cutting edge radius influence,temperature and internal stresses distribution are studied using finite element modelling.Numerical model setting is made with experimental and literature data using forces andchip thickness.

Observations of about 60 binary star systems hosting exoplanets indicate the necessity of stability studies of planetary motion in such multi-stellar systems. For wide binary systems with separations between hundreds and thousands of AU, the results from single-star systems may be applicable but, in tight double stars systems, we have to take the stellar interactions into account which influences the planetary motion significantly.

This review discusses the different types of planetary motion in double stars and the stability of the planets for different binary configurations. An application to the most famous tight binary system (γ Cephei) is also shown. Finally, we analyze the habitability from the dynamical point of view in such systems, where we discuss the motion of terrestrial-like planets in the so-called habitable zone.

We propose a new approach for a solution of the planet transit problem that is based on numerical calculation of integrals. The paper presents our method for the case of linear limb-darkening law and orbital inclination i = 90°, and illustrates the work of the code PTS written on the basis of our approach.

Photospheric thermal radiation components from gamma-ray burst (GRB) jets are estimated based on relativistic hydrodynamic simulations of jet propagation. The light curves and spectra are derived, considering viewing angle effects. The light curves exhibit several seconds time variability and the luminosity is as large as that of GRB prompt emission. For observers at a viewing angle of several degrees the spectra below the peak energy are much softer than that of Planck distribution and close to typical GRB spectrum. Whereas the spectra for observers at small viewing angle are hard and close to Planck distribution. Numerical Amati and Yonetoku relations are reproduced.