We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

We study the slowing down of particle beams passing through the dusty plasma with power-law kappa-distributions. Three plasma components, electrons, ions and dust particles, can have different kappa-parameter. We derive the deceleration factor (the velocity moment equation) and the slowing down time of a test particle, and numerically study the slowing down properties of an electron beam, a proton beam and a dust particle beam, respectively, in the kappa-distributed dusty plasma. We show that the slowing down properties of particle beams depend strongly on the kappa-parameters of the plasma components, and the dust component plays a dominant role in the slowing down. And the slowing down also depends on mass and charge of the dust particles in the dusty plasma. More detailed results are shown in 17 numerical graphs.

The equations governing motion of the satellite under the effect of oblateness of Earth and atmospheric drag have been simulated, for a fixed initial position and three different initial velocities, till satellite collapses on Earth. Simulation of motion of artificial Earth satellite subject to the combined effects of oblate Earth and atmospheric drag is presented. The atmospheric model considered here takes in to account of exponential variation of the density with initial distance of Satellite from Earth's surface, scale height and radial distance. The minimum and maximum values of orbital elements and their variation over a time for different initial velocities have been reported.

Although turbulence has been conjectured to be important for magnetic reconnection, still very little is known about its role in collisionless plasmas. Previous attempts to quantify the effect of turbulence on reconnection usually prescribed Alfv\'enic or other low-frequency fluctuations or investigated collisionless kinetic effects in just two-dimensional configurations and antiparallel magnetic fields. In view of this, we analyzed the kinetic turbulence self-generated by three-dimensional guide-field reconnection through force-free current sheets in frequency and wavenumber spaces, utilizing 3D particle-in-cell code numerical simulations. Our investigations reveal reconnection rates and kinetic turbulence with features similar to those obtained by current in-situ spacecraft observations of MMS as well as in the laboratory reconnection experiments MRX, VTF and \textsc{Vineta}-II. In particular we found that the kinetic turbulence developing in the course of 3D guide-field reconnection exhibits a broadband power-law spectrum extending beyond the lower-hybrid frequency and up to the electron frequencies. In the frequency space the spectral index of the turbulence appeared to be close to -2.8 at the reconnection X-line. In the wavenumber space it also becomes -2.8 as soon as the normalized reconnection rate reaches 0.1. The broadband kinetic turbulence is mainly due to current-streaming and electron-flow-shear instabilities excited in the sufficiently thin current sheets of kinetic reconnection. The growth of the kinetic turbulence corresponds to high reconnection rates which exceed those of fast laminar, non-turbulent reconnection.

Free-space optical communications have held the promise of revolutionizing space communications for a long time. The benefits of increasing the bitrate while reducing the volume, mass and energy of the space terminals have attracted the attention of many researchers for a long time. In the last few years, more and more technology demonstrations have been taking place with participants from both the public and the private sector. The National Institute of Information and Communications Technology (NICT) in Japan has a long experience in this field. SOTA (Small Optical TrAnsponder) was the last NICT space lasercom mission, designed to demonstrate the potential of this technology applied to microsatellites. Since the beginning of SOTA mission in 2014, NICT regularly established communication using the Optical Ground Stations (OGS) located in the Headquarters at Koganei (Tokyo) to receive the SOTA signals, with over one hundred successful links. All the goals of the SOTA mission were fulfilled, including up to 10-Mbit/s downlinks using two different wavelengths and apertures, coarse and fine tracking of the OGS beacon, space-to-ground transmission of the on-board-camera images, experiments with different error correcting codes, interoperability with other international OGS, and experiments on quantum communications. The SOTA mission ended on November 2016, more than doubling the designed lifetime of 1-year. In this paper, the SOTA characteristics and basic operation are explained, along with the most relevant technological demonstrations.

The large-scale magnetic cloud such as coronal mass ejections (CMEs) is the fundamental driver of the space weather. The interaction of the multiple CMEs in interplanetary space affects their dynamic evolution and geo-effectiveness. The complex and merged multiple magnetic clouds appear as the in-situ signature of the interacting CMEs. The Alfven waves are speculated to be one of the major possible energy exchange/dissipation mechanism during the interaction. However, no such observational evidence has been found in the literature. The case studies of CME-CME collision events suggest that the magnetic and thermal energy of the CME is converted into the kinetic energy. Moreover, the magnetic reconnection process is justified to be responsible for the merging of multiple magnetic clouds. Here, we present unambiguous evidence of sunward torsional Alfven waves in the interacting region after the super-elastic collision of multiple CMEs. The Walen relation is used to confirm the presence of Alfven waves in the interacting region of multiple CMEs/magnetic clouds. We conclude that Alfven waves and magnetic reconnection are the possible energy exchange/dissipation mechanisms during large-scale magnetic clouds collisions. The present study has significant implications not only in CME-magnetosphere interactions but also in the interstellar medium where interactions of large-scale magnetic clouds are possible.

Using spatio-temporal spectra we show direct evidence of excitation of magnetosonic and Alfv\'en waves in three-dimensional compressible magnetohydrodynamic turbulence at small Mach numbers. For the plasma pressure dominated regime, or high $\beta$ regime (with $\beta$ the ratio between fluid and magnetic pressure), and for the magnetic pressure dominated regime, or low $\beta$ regime, we study magnetic field fluctuations parallel and perpendicular to a guide magnetic field $\textbf{B}_0$. In the low $\beta$ case we find excitation of compressible and incompressible fluctuations, with a transfer of energy towards Alfv\'enic modes and to a lesser extent towards magnetosonic modes. In particular, we find signatures of the presence of fast magnetosonic waves in a scenario compatible with that of weak turbulence. In the high $\beta$ case, fast and slow magnetosonic waves are present, with no clear trace of Alfv\'en waves, and a significant part of the energy is carried by two-dimensional turbulent eddies.

It is very important to understand stochastic diffusion of energetic charged particles in non-uniform background magnetic field in plasmas of astrophysics and fusion devices. Using different methods considering along-field adiabatic focusing effect, various authors derived parallel diffusion coefficient $\kappa_\parallel$ and its correction $T$ to $\kappa_{\parallel 0}$, where $\kappa_{\parallel 0}$ is the parallel diffusion coefficient without adiabatic focusing effect. In this paper, using the improved perturbation method developed by He \& Schlickeiser and iteration process, we obtain a new correction $T'$ to $\kappa_{\parallel 0}$. Furthermore, by employing the isotropic pitch-angle scattering model $D_{\mu\mu}=D(1-\mu^2)$, we find that $T'$ has the different sign as that of $T$. In this paper the spatial perpendicular diffusion coefficient $\kappa_\bot$ with the adiabatic focusing effect is also obtained.

We report a newly found two-stage mechanism of electron acceleration near X-lines of 3D collisionless guide-field magnetic reconnection in the non-relativistic regime typical, e.g., for stellar coronae. We found that after electrons are first pre-accelerated during the linear growth of reconnection, they become additionally accelerated in the course of the nonlinear stage of 3D guide-field magnetic reconnection. This additional acceleration is due to the filamentation of electric and magnetic fields caused by streaming instabilities. In addition to enhanced parallel electric fields, the filamentation leads to additional curvature-driven electron acceleration in the guide-field direction. As a result, part of the the accelerated electron spectra becomes a power law with a spectral index of $\sim-1.6$ near the X-line. This second stage of acceleration due to nonlinear reconnection is relevant for the production of energetic electrons in, e.g., thin current sheets of stellar coronae.

We report on the parallel analysis of the periodic behaviour of coronal mass ejections (CMEs) based on 21 years [1996 -- 2016] of observations with the SOHO/LASCO--C2 coronagraph, solar flares, prominences, and several proxies of solar activity. We consider values of the rates globally and whenever possible, distinguish solar hemispheres and solar cycles 23 and 24. Periodicities are investigated using both frequency (periodogram) and time-frequency (wavelet) analysis. We find that these different processes, in addition to following the $\approx$11-year Solar Cycle, exhibit diverse statistically significant oscillations with properties common to all solar, coronal, and heliospheric processes: variable periodicity, intermittence, asymmetric development in the northern and southern solar hemispheres, and largest amplitudes during the maximum phase of solar cycles, being more pronounced during solar cycle 23 than the weaker cycle 24. However, our analysis reveals an extremely complex and diverse situation. For instance, there exists very limited commonality for periods of less than one year. The few exceptions are the periods of 3.1--3.2 months found in the global occurrence rates of CMEs and in the sunspot area (SSA) and those of 5.9--6.1 months found in the northern hemisphere. Mid-range periods of $\approx$1 and $\approx$2 years are more wide spread among the studied processes, but exhibit a very distinct behaviour with the first one being present only in the northern hemisphere and the second one only in the southern hemisphere. These periodic behaviours likely results from the complexity of the underlying physical processes, prominently the emergence of magnetic flux.

Using particle-scale simulations of non-suspended sediment transport for a large range of Newtonian fluids driving transport, including air and water, we determine the bulk transport cessation threshold $\Theta^r_t$ by extrapolating the transport load as a function of the dimensionless fluid shear stress (`Shields number') $\Theta$ to the vanishing transport limit. In this limit, the simulated steady states of continuous transport can be described by simple analytical model equations relating the average transport layer properties to the law of the wall flow velocity profile. We use this model to calculate $\Theta^r_t$ for arbitrary environments and derive a general Shields-like threshold diagram in which a Stokes-like number replaces the particle Reynolds number. Despite the simplicity of our hydrodynamic description, the predicted cessation threshold, both from the simulations and analytical model, quantitatively agrees with measurements for transport in air and viscous and turbulent liquids despite not being fitted to these measurements. We interpret the analytical model as a description of a continuous rebound motion of transported particles and thus $\Theta^r_t$ as the minimal fluid shear stress needed to compensate the average energy loss of transported particles during an average rebound at the bed surface. This interpretation, supported by simulations near $\Theta^r_t$, implies that entrainment mechanisms are needed to sustain transport above $\Theta^r_t$. While entrainment by turbulent events sustains intermittent transport, entrainment by particle-bed impacts sustains continuous transport. Combining our interpretations with the critical energy criterion for incipient motion by Valyrakis and coworkers, we put forward a new conceptual picture of sediment transport intermittency.

The solar modulation effect of cosmic rays in the heliosphere is an energy-, time-, and particle-dependent phenomenon which arises from a combination of basic particle transport processes such as diffusion, convection, adiabatic cooling, and drift motion. Making use of a large collection of time-resolved cosmic-ray data from recent space missions, we construct a simple predictive model of solar modulation which depends on direct solar-physics inputs: the number of solar sunspots and the tilt angle of the heliospheric current sheet. Under this framework, we present calculations of cosmic-ray proton spectra, positron/electron and antiproton/proton ratios and their time dependence in connection with the evolving solar activity. We report evidence for a time-lag $\Delta{T}=8.1\pm\,1.2$ months, between solar activity data and cosmic-ray flux measurements in space, which reflects the dynamics of the formation of the modulation region. This result enables us to forecast the cosmic-ray flux near Earth well in advance by monitoring solar activity

A phenomenological model is presented for the quantitative description of individual solar cycles' features, such as onset, intensity, evolution, in terms of the number of M and X-class solar flares. The main elements of the model are the relative ecliptic motion of the planets Jupiter and Saturn, and its synergy with a quasi-periodic component of solar activity. Using as input the temporal distribution of flares during cycle 21, the general evolution of cycles 22-24 is reproduced in notable agreement with the observations, including the resurgence of activity in the last months of 2017, and further predictions are provided for cycle 25. This deterministic description could contribute to elucidating the responsible physical mechanisms and forecasting space weather.

We analyze a high-resolution simulation of magnetopause reconnection observed by the Magnetospheric Multiscale (MMS) mission and explain the occurrence of strongly localized dissipation with an amplitude more than an order of magnitude larger than expected. Unlike symmetric reconnection, wherein reconnection of the ambient reversed magnetic field drives the dissipation, we find the annihilation of the self-generated, out-of-plane (Hall) magnetic field plays the dominant role. Electrons flow along the magnetosheath separatrices, converge in the diffusion region, and jet past the X-point into the magnetosphere. The resulting accumulation of negative charge generates intense parallel electric fields that eject electrons along the magnetospheric separatrices and produce field-aligned beams. Many of these features match MMS observations.

As part of the dynamical analysis carried out within the Horizon 2020 ReDSHIFT project, this work analyzes the possible strategies to guide low altitude satellites towards an atmospheric reentry through an impulsive maneuver. We consider a fine grid of initial conditions in semi-major axis, eccentricity and inclination and we identify the orbits that can be compliant with the 25-year rule as the target of a single-burn strategy. Besides the atmospheric drag, we look for the aid provided by other dynamical perturbations - mainly solar radiation pressure - to facilitate a reentry. Indeed, in the case of typical area-to-mass ratios for objects in LEO, we observed that dynamical resonances can be considered only in combination with the atmospheric drag and for a very limited set of initial orbits. Instead, if an area augmentation device, as a solar sail, is available on-board the spacecraft, we verified that a wider range of disposal solutions become available. This information is exploited to design an improved mitigation scheme, that can be applied to any satellite in LEO.

The recent progress in laser propulsion research has advanced substantially the prospects to realize interstellar spaceflight within a few decades. Here we examine passive deceleration via momentum braking from ionized interstellar media. The very large area to mass relations needed as a consequence of the low interstellar densities, of the order of 0.1 particles per $\mathrm{cm}^{3}$, or lower, are potentially realizable with magnetic sails generated by superconducting coils. Integrating the equations of motion for interstellar protons hitting a Biot Savart loop we evaluate the effective reflection area $A(v)$ in terms of the velocity $v$ of the craft. We find that the numerical data is fitted over two orders of magnitude by the scaling relation $A(v)\ =\ 0.081A_R\log^3(I/(\beta I_c))$, where $A_R=\pi R^2$ is the bare sail area, $I$ the current and $\beta=v/c$. The critical current $I_c$ is $1.55\cdot10^6$ Ampere. The resulting universal deceleration profile can be evaluated analytically and mission parameters optimized for a minimal craft mass. For the case of a sample high speed transit to Proxima Centauri we find that magnetic momentum braking would involve daunting mass requirements of the order of $10^3$ tons. A low speed mission to the Trappist-1 system could be realized on the other side already with a 1.5 ton spacecraft, which would be furthermore compatible with the specifications of currently envisioned directed energy launch systems. The extended cruising times of the order of $10^4$ years imply however that a mission to the Trappist-1 system would be viable only for mission concepts for which time constrains are not relevant.

Two main models have been developed to explain the mechanisms of release, heating and acceleration of the nascent solar wind, the wave-turbulence-driven (WTD) models and reconnection-loop-opening (RLO) models, in which the plasma release processes are fundamentally different. Given that the statistical observational properties of helium ions produced in magnetically diverse solar regions could provide valuable information for the solar wind modelling, we examine the statistical properties of the helium abundance (A_He) and the speed difference between helium ions and protons (v_alpha,p) for coronal holes (CHs), active regions (ARs) and the quiet Sun (QS). We find bimodal distributions in the space of A_He and v_alpha,p/v_A (where v_A is the local Alfven speed)for the solar wind as a whole. The CH wind measurements are concentrated at higher A_He and v_alpha,p/v_A values with a smaller A_He distribution range, while the AR and QS wind is associated with lower A_He and v_alpha,p/v_A, and a larger A_He distribution range. The magnetic diversity of the source regions and the physical processes related to it are possibly responsible for the different properties of A_He and v_alpha,p/v_A. The statistical results suggest that the two solar wind generation mechanisms, WTD and RLO, work in parallel in all solar wind source regions. In CH regions WTD plays a major role, whereas the RLO mechanism is more important in AR and QS.

We investigate how a weak coronal mass ejection (CME) launched on 2016 October 8 without obvious signatures in the low corona produced a relatively intense geomagnetic storm. Remote sensing observations from SDO, STEREO and SOHO and in situ measurements from WIND are employed to track the CME from the Sun to the Earth. Using a graduated cylindrical shell (GCS) model, we estimate the propagation direction and the morphology of the CME near the Sun. CME kinematics are determined from the wide-angle imaging observations of STEREO A and are used to predict the CME arrival time and speed at the Earth. We compare ENLIL MHD simulation results with in situ measurements to illustrate the background solar wind where the CME was propagating. We also apply a Grad--Shafranov technique to reconstruct the flux rope structure from in situ measurements in order to understand the geo-effectiveness associated with the CME magnetic field structure. Key results are obtained concerning how a weak CME can generate a relatively intense geomagnetic storm: (1) there were coronal holes at low latitudes, which could produce high speed streams (HSSs) to interact with the CME in interplanetary space; (2) the CME was bracketed between a slow wind ahead and a HSS behind, which enhanced the southward magnetic field inside the CME and gave rise to the unexpected geomagnetic storm.

On 10 September 2017 Neutron Monitors (NMs) apparatus located at ground level and high latitudes detected an increase in the counting rate associated to solar energetic particles (SEP) emission from X8.2 class solar flare and its associated CME. This was the second highest flare of the current solar cycle. The origin was the active region AR 2673 when it was located at the edge of the west solar disc, without a direct magnetic connection with Earth. However, there was a peculiar condition, the solar protons accelerated by the CME shocks were injected within a heliospheric current sheet (HCS) region and when the Earth was crossing this region. We show that the HCS played the role of a magnetic connection between Sun and Earth. The fast-moving protons penetrated the magnetosphere near the north and south poles and their interaction in the atmosphere produced a ground level enhancement (GLE \#72), this was the second in the current solar cycle. In the NOAA Space Weather Scale, SEP is cataloged as radiation storms, from S1(minor) to S5 (extreme), and in the present case, it reached S3 (strong). In addition, four days before, on 6 September, AR 2673 erupted, was an X9.3-class flare (the strongest flare of the current cycle). It was associated with a halo CME toward Earth, triggering the second strongest geomagnetic storm of the current solar cycle on 7 September. We detail these observations.

It is widely reported that the power spectra of magnetic field and velocity fluctuations in the solar wind have power law scalings with inertial-range spectral indices of -5/3 and -3/2 respectively. Studies of solar wind turbulence have repeatedly demonstrated the impact of discontinuities and coherent structures on the measured spectral index. Whether or not such discontinuities are self-generated by the turbulence or simply observations of advected structures from the inner heliosphere has been a matter of considerable debate. This work presents a statistical study of magnetic field and velocity spectral indices over 10 years of solar-wind observations; we find that anomalously steep magnetic spectra occur in magnetically dominated intervals with negative residual energy. However, this increase in negative residual energy has no noticeable impact on the spectral index of the velocity fluctuations, suggesting that these intervals with negative residual energy correspond to intermittent magnetic structures. We show statistically that the difference between magnetic and velocity spectral indices is a monotonic function of residual energy, consistent with previous work which suggests that intermittency in fluctuations causes spectral steepening. Additionally, a statistical analysis of cross helicity demonstrates that when the turbulence is balanced (low cross-helicity), the magnetic and velocity spectral indices are not equal, which suggests that our observations of negative residual energy and intermittent structures are related to non-linear turbulent interactions rather than the presence of advected pre-existing flux-tube structures.

This paper provides an overview of the radiation aspects of manned space flight to Moon and Mars. The expected ionising radiation dose for an astronaut is assessed along the Apollo 11 flight path to the Moon. With the two dose values, the expected and the measured total dose, the radiation shielding and the activity of the Sun are estimated. To judge the risk or safety margin the radiation effects on humans are opposed. By only changing the flight duration an assessment for the total dose of a journey to Mars is provided.

According to the solar protons' data observed by Geostationary Operational Environmental Satellites (GOES) and neutron monitors on the ground of the Earth, and near-relativistic electrons data measured by the ACE spacecraft, the onset times of protons with different energies and near-relativistic electrons have been estimated and compared with the time of solar soft and hard X-ray and radio burst data. The results show that first arriving relativistic and non-relativistic protons and electrons may be accelerated by the concurrent flare. The results also suggest that release times of protons with different energies may be different, and the protons with lower energy may release earlier than those with higher energy. Some protons accelerated by concurrent flares may be further accelerated by CME-driven shock.

BOOST (BOOst Symmetry Test) is a proposed satellite mission to search for violations of Lorentz invariance by comparing two optical frequency references. One is based on a long-term stable optical resonator and the other on a hyperfine transition in molecular iodine. This mission will allow to determine several parameters of the standard model extension in the electron sector up to two orders of magnitude better than with the current best experiments. Here, we will give an overview of the mission, the science case and the payload.

We present a three-dimensional study of the plasma dynamics at the flank magnetopause of the Earth's magnetosphere during mainly northward interplanetary magnetic field (IMF) periods. Two-fluid simulations show that the initial magnetic shear at the magnetopause and the field line bending caused by the dynamics itself (in a configuration taken as representative of the properties of the flank magnetopause) influence both the location where the Kelvin-Helmholtz (KH) instability and the induced magnetic reconnection take place and their nonlinear development. The KH vortices develop asymmetrically with respect to the Earth's equatorial plane where the local KH linear growth rate is maximal. Vortex driven reconnection processes take place at different latitudes, ranging from the equatorial plane to mid-latitude regions, but only in the hemisphere that turns out to be the less KH unstable. These results suggest that KH-induced reconnection is not limited to specific regions around the vortices (inside, below or above), but may be triggered over a broad and continuous range of locations in the vicinity of the vortices.

We present detailed analysis of the turbulence observed in three-dimensional particle-in-cell simulations of magnetic reconnection at the magnetopause. The parameters are representative of an electron diffusion region encounter of the Magnetospheric Multiscale (MMS) mission. The turbulence is found to develop around both the magnetic X line and separatrices, is electromagnetic in nature, is characterized by a wave vector $k$ given by $k\rho_e\sim(m_eT_e/m_iT_i)^{0.25}$ with $\rho_e$ the electron Larmor radius, and appears to have the ion pressure gradient as its source of free energy. Taken together, these results suggest the instability is a variant of the lower hybrid drift instability. The turbulence produces electric field fluctuations in the out-of-plane direction (the direction of the reconnection electric field) with an amplitude of around $\pm 10$~mV/m, which is much greater than the reconnection electric field of around $0.1$~mV/m. Such large values of the out-of-plane electric field have been identified in the MMS data. The turbulence in the simulations controls the scale lengths of the density profile and current layers in asymmetric reconnection, driving them closer to $\sqrt{\rho_e\rho_i}$ than the $\rho_e$ or $d_e$ scalings seen in 2-D reconnection simulations, and produces significant anomalous resistivity and viscosity in the electron diffusion region.