The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.

The collision of magnetic reconnection jets is studied by means of a three dimensional numerical simulation at kinetic scale, in the presence of a strong guide field. We show that turbulence develops due to the jets collision producing several current sheets in reconnection outflows, aligned with the guide field direction. The turbulence is mainly two-dimensional, with stronger gradients in the plane perpendicular to the guide field and a low wave-like activity in the parallel direction. First, we provide a numerical method to isolate the central turbulent region. Second, we analyze spatial second-order structure function and prove that turbulence is confined in this region. Finally, we compute local magnetic and electric frequency spectra, finding a trend in the sub-ion range that differs from typical cases for which the Taylor hypothesis is valid, as well as wave activity in the range between ion and electron cyclotron frequencies. Our results are relevant to understand observations of reconnection jets collisions in space plasmas.

Reconnection outflows are regions of intense recent scrutiny, from in situ observations and from simulations. These regions are host to a variety of instabilities and intense energy exchanges, often even superior to the main reconnection site. We report here a number of results drawn from investigation of simulations. First, the outflows are observed to become unstable to drift instabilities. Second, these instabilities lead to the formation of secondary reconnection sites. Third, the secondary processes are responsible for large energy exchanges and particle energization. Finally, the particle distribution function are modified to become non-Maxwellian and include multiple interpenetrating populations.

We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic Particle-in-Cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind; and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3-9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and lunar magnetic anomaly simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross-sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.

We performed three-dimensional Particle-in-Cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy was about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy was transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in the cluster of null points evolved in three phases. During the first phase, ion beams were excited, that then gave part of their energy back to magnetic field in the second phase. In the third phase, magnetic reconnection occurs in many small patches around the current channels formed along the stripes of low magnetic field. Magnetic reconnection in null points presents essentially three-dimensional features, with no two dimensional symmetries or current sheets.

A widely accepted scenario of magnetic reconnection in collisionless space plasmas is the breakage of magnetic field lines in X-points. In laboratory, reconnection is commonly studied in pinches, current channels embedded into twisted magnetic fields. No model of magnetic reconnection in space plasmas considers both null-points and pinches as peers. We have performed a particle-in-cell simulation of magnetic reconnection in a three-dimensional configuration where null-points are present initially, and Z-pinches are formed during the simulation along the lines of spiral null-points. The non-spiral null-points are more stable than spiral ones, and no substantial energy dissipation is associated with them. On the contrary, turbulent magnetic reconnection in the pinches causes the magnetic energy to decay at a rate of ~1.5% per ion gyro period. Dissipation in similar structures is a likely scenario in space plasmas with large fraction of spiral null-points.

We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of `intermittent turbulence' within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5\---10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as Magnetospheric Multiscale Mission (MMS).

We adopted an unstructured hydrodynamical solver CharLES to the problem of global convection in the Sun. With the aim to investigate the properties of solar turbulent convection and reproduce differential rotation pattern. We performed simulations in two spherical shells, with 1.3 and 10 million cells. In the first, coarse mesh, the solution does not reproduce realistic convection, and is dominated by numerical effects. In the second mesh, thermal conduction leads to cooling of bottom layers, that could not be compensated by solar irradiance. More simulations in the 10M cells mesh should be performed to investigate the influence of transport coefficients and numerical effects. Our estimate of the code performance suggests, that realistic simulations in even finer grids could be performed for reasonable computational cost.

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

The aim of this paper is to analyse the validity of the Dopplergram and lambda-meter techniques for the Doppler diagnostics of solar photospheric velocities using the BaII 4554 A line. Both techniques are evaluated by means of NLTE radiative transfer calculations of the BaII 4554 A line in a three-dimensional hydrodynamical model of solar convection. We consider the cases of spatially unsmeared profiles and the profiles smeared to the resolution of ground-based observations. We find that: (i) Speckle-reconstructed Dopplergram velocities reproduce the ``true'' velocities well at heights around 300 km, except for intergranular lanes with strong downflows where the velocity can be overestimated. (ii) The lambda-meter velocities give a good representation of the ``true'' velocities through the whole photosphere, both under the original and reduced spatial resolutions. The velocities derived from the inner wing of smeared BaII 4554 A line profiles are more reliable than those for the outer wing. Only under high spatial resolution does the inner wing velocities calculated in intergranular regions give an underestimate (or even a sign reversal) compared with the model velocities. (iii) NLTE effects should be taken into account in modelling the BaII 4554 A line profiles. Such effects are more pronounced in intergranular regions. Our analysis supports the opinion that the Dopplergram technique applied to the BaII 4554 A line is a valuable tool for the Doppler diagnostics of the middle photosphere around 300 km. The \lambda-meter technique applied to this line gives us a good opportunity to ``trace'' the non-thermal motions along the whole photosphere up to the temperature minimum and lower chromosphere.

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

We analyse the skewness of the line-of-sight velocity distributions in model elliptical galaxies built through collisionless galaxy mergers. We build the models using large N-body simulations of mergers between either two spiral or two elliptical galaxies. Our aim is to investigate whether the observed ranges of skewness coefficient (h3) and the rotational support (V/sigma), as well as the anticorrelation between h3 and V, may be reproduced through collisionless mergers. Previous attempts using N-body simulations failed to reach V/sigma ~ 1-2 and corresponding high h3 values, which suggested that gas dynamics and ensuing star formation might be needed in order to explain the skewness properties of ellipticals through mergers. Here we show that high V/sigma and high h3 are reproduced in collisionless spiral-spiral mergers whenever a central bulge allows the discs to retain some of their original angular momentum during the merger. We also show that elliptical-elliptical mergers, unless merging from a high-angular momentum orbit, reproduce the strong skewness observed in non-rotating, giant, boxy ellipticals. The behaviour of the h3 coefficient therefore associates rapidly-rotating disky ellipticals to disc-disc mergers, and associates boxy, slowly-rotating giant ellipticals to elliptical-elliptical mergers, a framework generally consistent with the expectations of hierarchical galaxy formation.