
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
fifthordermonotonicitypreserving 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
twodimensional problems showed the advantages of using the highorder scheme
by comparing with results from a standard secondorder TVD MUSCL scheme. The
present code enabled us to explore longterm evolution of a threedimensional
accretion disk around a black hole, in which compressible MHD turbulence caused
continuous mass accretion via nonlinear growth of the magnetorotational
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 powerlaw kappadistributions. Three plasma components, electrons, ions
and dust particles, can have different kappaparameter. 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
kappadistributed dusty plasma. We show that the slowing down properties of
particle beams depend strongly on the kappaparameters 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 lowfrequency fluctuations or
investigated collisionless kinetic effects in just twodimensional
configurations and antiparallel magnetic fields. In view of this, we analyzed
the kinetic turbulence selfgenerated by threedimensional guidefield
reconnection through forcefree current sheets in frequency and wavenumber
spaces, utilizing 3D particleincell code numerical simulations. Our
investigations reveal reconnection rates and kinetic turbulence with features
similar to those obtained by current insitu 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 guidefield reconnection exhibits a broadband
powerlaw spectrum extending beyond the lowerhybrid 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 Xline. 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
currentstreaming and electronflowshear 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, nonturbulent reconnection.

Freespace 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 10Mbit/s downlinks using two different
wavelengths and apertures, coarse and fine tracking of the OGS beacon,
spacetoground transmission of the onboardcamera 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 1year. In this
paper, the SOTA characteristics and basic operation are explained, along with
the most relevant technological demonstrations.

The largescale 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 geoeffectiveness.
The complex and merged multiple magnetic clouds appear as the insitu 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 CMECME 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 superelastic 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 largescale magnetic clouds collisions.
The present study has significant implications not only in CMEmagnetosphere
interactions but also in the interstellar medium where interactions of
largescale magnetic clouds are possible.

Using spatiotemporal spectra we show direct evidence of excitation of
magnetosonic and Alfv\'en waves in threedimensional 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
twodimensional turbulent eddies.

It is very important to understand stochastic diffusion of energetic charged
particles in nonuniform background magnetic field in plasmas of astrophysics
and fusion devices. Using different methods considering alongfield 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 pitchangle 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 twostage mechanism of electron acceleration near
Xlines of 3D collisionless guidefield magnetic reconnection in the
nonrelativistic regime typical, e.g., for stellar coronae. We found that after
electrons are first preaccelerated during the linear growth of reconnection,
they become additionally accelerated in the course of the nonlinear stage of 3D
guidefield 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 curvaturedriven electron acceleration in the
guidefield direction. As a result, part of the the accelerated electron
spectra becomes a power law with a spectral index of $\sim1.6$ near the
Xline. 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/LASCOC2 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 timefrequency (wavelet)
analysis. We find that these different processes, in addition to following the
$\approx$11year 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.13.2
months found in the global occurrence rates of CMEs and in the sunspot area
(SSA) and those of 5.96.1 months found in the northern hemisphere. Midrange
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 particlescale simulations of nonsuspended 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 Shieldslike
threshold diagram in which a Stokeslike 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 particlebed 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 particledependent 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 timeresolved
cosmicray data from recent space missions, we construct a simple predictive
model of solar modulation which depends on direct solarphysics inputs: the
number of solar sunspots and the tilt angle of the heliospheric current sheet.
Under this framework, we present calculations of cosmicray proton spectra,
positron/electron and antiproton/proton ratios and their time dependence in
connection with the evolving solar activity. We report evidence for a timelag
$\Delta{T}=8.1\pm\,1.2$ months, between solar activity data and cosmicray flux
measurements in space, which reflects the dynamics of the formation of the
modulation region. This result enables us to forecast the cosmicray 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 Xclass solar flares. The main elements of the
model are the relative ecliptic motion of the planets Jupiter and Saturn, and
its synergy with a quasiperiodic component of solar activity. Using as input
the temporal distribution of flares during cycle 21, the general evolution of
cycles 2224 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 highresolution 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 selfgenerated, outofplane (Hall) magnetic field
plays the dominant role. Electrons flow along the magnetosheath separatrices,
converge in the diffusion region, and jet past the Xpoint into the
magnetosphere. The resulting accumulation of negative charge generates intense
parallel electric fields that eject electrons along the magnetospheric
separatrices and produce fieldaligned 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 semimajor axis,
eccentricity and inclination and we identify the orbits that can be compliant
with the 25year rule as the target of a singleburn 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 areatomass 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 onboard 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 Trappist1 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 Trappist1 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 waveturbulencedriven
(WTD) models and reconnectionloopopening (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 wideangle 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
GradShafranov technique to reconstruct the flux rope structure from in situ
measurements in order to understand the geoeffectiveness 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 fastmoving 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.3class 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 inertialrange
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 selfgenerated 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 solarwind 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 crosshelicity), the magnetic and velocity spectral indices are
not equal, which suggests that our observations of negative residual energy and
intermittent structures are related to nonlinear turbulent interactions rather
than the presence of advected preexisting fluxtube 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 nearrelativistic electrons data measured by the ACE spacecraft, the
onset times of protons with different energies and nearrelativistic electrons
have been estimated and compared with the time of solar soft and hard Xray and
radio burst data. The results show that first arriving relativistic and
nonrelativistic 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 CMEdriven 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 longterm 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 threedimensional study of the plasma dynamics at the flank
magnetopause of the Earth's magnetosphere during mainly northward
interplanetary magnetic field (IMF) periods. Twofluid 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
KelvinHelmholtz (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 midlatitude regions, but only
in the hemisphere that turns out to be the less KH unstable. These results
suggest that KHinduced 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 threedimensional
particleincell 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 outofplane 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
outofplane 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 2D
reconnection simulations, and produces significant anomalous resistivity and
viscosity in the electron diffusion region.