In the literature different so-called $\kappa$-distribution functions are discussed to fit and model the velocity (or energy) distributions of solar wind species, pickup ions or magnetospheric particles. Here we introduce a generalized (isotropic) $\kappa$-distribution as a "cookbook", which admits as special cases, or "recipes", all the other known versions of $\kappa$-models. A detailed analysis of the generalized distribution function is performed, providing general analytical expressions for the velocity moments, Debye length, and entropy, and pointing out a series of general requirements that plasma distribution functions should satisfy. From a contrasting analysis of the recipes found in the literature, we show that all of them lead to almost the same macroscopic parameters with a small standard deviation between them. However, one of these recipes called the regularized $\kappa$-distribution provides a functional alternative for macroscopic parameterization without any constraint for the power-law exponent $\kappa$.

Observations of the time-dependent cosmic-ray Sun shadow have been proven as a valuable diagnostic for the assessment of solar magnetic field models. In this paper, seven years of IceCube data are compared to solar activity and solar magnetic field models. A quantitative comparison of solar magnetic field models with IceCube data on the event rate level is performed for the first time. Additionally, a first energy-dependent analysis is presented and compared to recent predictions. We use seven years of IceCube data for the Moon and the Sun and compare them to simulations on data rate level. The simulations are performed for the geometrical shadow hypothesis for the Moon and the Sun and for a cosmic-ray propagation model governed by the solar magnetic field for the case of the Sun. We find that a linearly decreasing relationship between Sun shadow strength and solar activity is preferred over a constant relationship at the 6.4sigma level. We test two commonly used models of the coronal magnetic field, both combined with a Parker spiral, by modeling cosmic-ray propagation in the solar magnetic field. Both models predict a weakening of the shadow in times of high solar activity as it is also visible in the data. We find tensions with the data on the order of $3\sigma$ for both models, assuming only statistical uncertainties. The magnetic field model CSSS fits the data slightly better than the PFSS model. This is generally consistent with what is found previously by the Tibet AS-gamma Experiment, a deviation of the data from the two models is, however, not significant at this point. Regarding the energy dependence of the Sun shadow, we find indications that the shadowing effect increases with energy during times of high solar activity, in agreement with theoretical predictions.

With upcoming missions such as the James Webb Space Telescope (JWST), the European Extremely Large Telescope (ELT), and the Atmospheric Remote-sensing Infrared Exoplanet Large-survey (ARIEL), we soon will be on the verge of detecting and characterizing Earth-like exoplanetary atmospheres for the first time. These planets are most likely to be found around smaller and cooler K- and M-type stars. However, recent observations showed that their radiation environment might be much harsher than that of the Sun. Thus, the exoplanets are most likely exposed to an enhanced stellar radiation environment, which could affect their habitability, for example, in the form of a hazardous flux of energetic particles. Knowing the stellar radiation field, and being able to model the radiation exposure on the surface of a planet is crucial to assess its habitability. In this study, we present 3D magnetohydrodynamic (MHD)-based model efforts investigating M-stars, focusing on V374 Peg, Proxima Centauri, and LHS 1140, chosen because of their diverse astrospheric quantities. We show that V374 Peg has a much larger astrosphere (ASP) than our Sun, while Proxima Centauri and LHS 1140 most likely have ASPs comparable or even much smaller than the heliosphere, respectively. Based on a 1D transport model, for the first time, we provide numerical estimates of the modulation of Galactic cosmic rays (GCRs) within the three ASPs. We show that the impact of GCRs on the Earth-like exoplanets Proxima Centauri b and LHS 1140 b cannot be neglected in the context of exoplanetary habitability.

A new class of one-dimensional solar wind models is developed within the general polytropic, single-fluid hydrodynamic framework. The particular case of quasi-adiabatic radial expansion with a localized heating source is considered. We consider analytical solutions with continuous Mach number over the entire radial domain while allowing for jumps in the flow velocity, density, and temperature, provided that there exists an external source of energy in the vicinity of the critical point which supports such jumps in physical quantities. This is substantially distinct from both the standard Parker solar wind model and the original nozzle solutions, where such discontinuous solutions are not permissible. We obtain novel sample analytic solutions of the governing equations corresponding to both slow and fast wind.

The interpretation of recent observations of bow shocks around O-stars and the creation of corresponding models require a detailed understanding of the associated (magneto-)hydrodynamic structures. We base our study on three-dimensional numerical magneto-)hydrodynamical models, which are analyzed using the dynamically relevant parameters, in particular, the (magneto)sonic Mach numbers. The analytic Rankine-Hugoniot relation for HD and MHD are compared with those obtained by the numerical model. In that context we also show that the only distance which can be approximately determined is that of the termination shock, if it is a hydrodynamical shock. For MHD shocks the stagnation point does not, in general, lie on the inflow line, which is the line parallel to the inflow vector and passing through the star. Thus an estimate via the Bernoulli equation as in the HD case is, in general, not possible. We also show that in O-star astrospheres, distinct regions exist in which the fast, slow, Alfv\'enic, and sonic Mach numbers become lower than one, implying sub-slow magnetosonic as well as sub-fast and sub-sonic flows. Nevertheless, the analytic MHD Rankine Hugoniot relations can be used for further studies of turbulence and cosmic ray modulation.

Three-dimensional models of astrospheres have recently become of interest. However, comparisons between these models and observations are non-trivial because of the two-dimensional nature of observations. By projecting selected physical values of three-dimensional models of astrospheres onto the surface of a sphere that is centred on a virtual all-sky observer, these models can be compared to observational data in different observables: the column density, bremsstrahlung flux, rotation measure, H$\alpha$ flux, and synchrotron or cyclotron flux. Projections were calculated by rotating and moving the astrosphere model to the desired position and orientation and by then computing the value of a given patch on the sphere by a modified line-of-sight integration. Contributions to the selected observable made by all model cells that are connected to the patch by the line of sight in question were taken into account. When the model produces a bow shock, a distinct parabolic structure produced by the outer astrosheath can be seen in every observable of the projection, the exact shape depending on the orientations of the line of sight and the stellar motion. Of all four examined astrosphere models, only that of $\lambda$ Cephei shows fluxes that are higher than current observational thresholds. This is due to the strong stellar wind and interstellar inflow of the $\lambda$ Cephei model.

Several concepts for heliospheric missions operating at heliocentric distances far beyond Earth orbit are currently investigated by the scientific community. The mission concept of the Interstellar Probe (McNutt et al. 2018), e.g., aims at reaching a distance of 1000 au away from the Sun within this century. This would allow the coming generation to obtain a global view of our heliosphere from an outside vantage point by measuring the Energetic Neutral Atoms (ENAs) originating from the various plasma regions. It would also allow for direct sampling of unperturbed interstellar medium, and for many observation opportunities beyond heliospheric science, such as visits to Kuiper Belt Objects, a comprehensive view on the interplanetary dust populations, and infrared astronomy free from the foreground emission of the Zodiacal cloud. In this study, we present a simple empirical model of ENAs from the heliosphere and derive basic requirements for ENA instrumentation onboard a spacecraft at great heliocentric distances. We consider the full energy range of heliospheric ENAs from 10 eV to 100 keV because each part of the energy spectrum has its own merits for heliospheric science. To cover the full ENA energy range, two or three different ENA instruments are needed. Thanks to parallax observations, some insights about the nature of the IBEX Ribbon and the dimensions of the heliosphere can already be gained by ENA imaging from a few au heliocentric distance. To directly reveal the global shape of the heliosphere, measurements from outside the heliosphere are, of course, the best option.

In the absence of efficient collisions, deviations from thermal equilibrium of plasma particle distributions are controlled by the self-generated instabilities. The whistler instability is a notorious example, usually responsible for the regulation of electron temperature anisotropy $A = T_{\perp}/T_\parallel>$ (with $\perp, \parallel$ respective to the magnetic field direction) observed in space plasmas, e.g., solar wind and planetary magnetospheres. Suprathermal electrons present in these environments change the plasma dispersion and stability properties, with expected consequences on the kinetic instabilities and the resulting fluctuations, which, in turn, scatter the electrons and reduce their anisotropy. In order to capture these mutual effects we use a quasilinear kinetic approach and PIC simulations, which provide a comprehensive characterization of the whistler instability under the influence of suprathermal electrons. Analysis is performed for a large variety of plasma conditions, ranging from low-beta plasmas encountered in outer corona or planetary magnetospheres to a high-beta solar wind characteristic to large heliospheric distances. Enhanced by the suprathermal electrons, whistler fluctuations stimulate the relaxation of temperature anisotropy, and this influence of suprathermals increases with plasma beta parameter.

The standard (non-relativistic) $\kappa$-distribution is widely used to fit data and to describe macroscopic thermodynamical behavior, e.g.\ the pressure (temperature) as the second moment of the distribution function. By contrast to a Maxwellian distribution, for small relevant values $\kappa < 2$ there exists a significant, but unphysical contribution to the pressure from unrealistic, superluminal particles with speeds exceeding the speed of light. Similar concerns exist for the entropy. We demonstrate here that by using the recently introduced regularized $\kappa$-distribution one can avoid such unphysical behaviour.

The cosmic-ray Sun shadow, which is caused by high-energy charged cosmic rays being blocked and deflected by the Sun and its magnetic field, has been observed by various experiments such as Argo-YBJ, HAWC, Tibet, and IceCube. Most notably, the shadow's size and depth was recently shown to correlate with the 11-year solar cycle. The interpretation of such measurements, which help to bridge the gap between solar physics and high-energy particle astrophysics, requires a solid theoretical understanding of cosmic-ray propagation in the coronal magnetic field. It is the aim of this paper to establish theoretical predictions for the cosmic-ray Sun shadow in order to identify observables that can be used to study this link in more detail. To determine the cosmic-ray Sun shadow, we numerically compute trajectories of charged cosmic rays in the energy range of 5 to 316 TeV for five different mass numbers. We present and analyse the resulting shadow images for protons and iron, as well as for typically measured cosmic-ray compositions. We confirm the observationally established correlation between the magnitude of the shadowing effect and both the mean sunspot number and the polarity of the magnetic field during the solar cycle. We also show that during low solar activity, the Sun's shadow behaves similarly to that of a dipole, for which we find a non-monotonous dependence on energy. In particular, the shadow can become significantly more pronounced than the geometrical disk expected for a totally unmagnetized Sun. For times of high solar activity, we instead predict the shadow to depend monotonously on energy, and to be generally weaker than the geometrical shadow for all tested energies. These effects should become visible in energy-resolved measurements of the Sun shadow, and may in the future become an independent measure for the level of disorder in the solar magnetic field.

In collision-poor plasmas from space, e.g., solar wind, terrestrial magnetospheres, kinetic instabilities are expected to play a major role in constraining the temperature anisotropy of plasma particles, but a definitive answer can be given only after ascertaining their properties in these environments. Present study describes the full spectrum of electron firehose instabilities in the presence of suprathermal electron populations which are ubiquitous in space plasmas. Suprathermal electrons stimulate both the periodic and aperiodic branches, remarkable being the effects shown by the aperiodic mode propagating obliquely to the ambient magnetic field which markedly exceeds the growth rates of the parallel (periodic) branch reported recently in Lazar et al., (2017a, MNRAS 464, 564). Derived exclusively in terms of the plasma parameters, the anisotropy thresholds of this instability are also lowered in the presence of suprathermal electrons, predicting an enhanced effectiveness in the solar wind conditions. These results may also be relevant in various other astrophysical contexts where the firehose instabilities involve, e.g., solar flares, sites of magnetic field reconnection, accretion flows or plasma jets leading to shocks and co-rotating interactions in heliosphere, interstellar medium and galaxy clusters.

While several studies have investigated large-scale cluster winds resulting from an intra-cluster interaction of multiple stellar winds, as yet they have not provided details of the bordering flows inside a given cluster. The present work explores the principal structure of the combined flow resulting from the interaction of multiple stellar winds inside stellar clusters. The theory of complex potentials is applied to analytically investigate stagnation points, boundaries between individual outflows, and the hydrodynamic structure of the asymptotic large-scale cluster wind. In a second part, these planar considerations are extended to fully three-dimensional, asymmetric configurations of wind-driving stars. We find (i) that one can distinguish regions in the large-scale cluster wind that are determined by the individual stellar winds, (ii) that there are comparatively narrow outflow channels, and (iii) that the large-scale cluster wind asymptotically approaches spherical symmetry at large distances. The combined flow inside a stellar cluster resulting from the interaction of multiple stellar winds is highly structured.

The multitude of recent multi-point spacecraft observations of solar energetic particle (SEP) events have made it possible to study the longitudinal distribution of SEPs in great detail. SEPs, even those accelerated during impulsive events, show a much wider than expected longitudinal extent, bringing into question the processes responsible for their transport perpendicular to the local magnetic field. In this paper we examine some aspects of perpendicular transport by including perpendicular diffusion into a numerical SEP transport model that simulates the propagation of impulsively accelerated SEP electrons in the ecliptic plane. We find that: (i) The pitch-angle dependence of the perpendicular diffusion coefficient is an important, and currently mainly overlooked, transport parameter. (ii) SEP intensities are generally asymmetric in longitude, being enhanced towards the west of optimal magnetic connection to the acceleration region. (iii) The maximum SEP intensity may also be shifted (parameter dependently) away from the longitude of best magnetic connectivity at 1 AU. We also calculate the maximum intensity, the time of maximum intensity, the onset time and the maximum anisotropy as a function of longitude at Earth's orbit and compare the results, in a qualitative fashion, to recent spacecraft observations.

The Voyager plasma observations show that the physics of the heliosheath is rather complex, and especially that the temperature derived from observation differs from expectations. To explain this fact the temperature in the heliosheath should be based on $\kappa$ distributions instead of Maxwellians because the former allows for much higher temperature. Here we show an easy way to calculate the $\kappa$ temperatures when those estimated from the data are given as Maxwellian temperatures. We use the moments of the Maxwellian and $\kappa$ distributions to estimate the $\kappa$ temperature. Moreover, we show that the pressure (temperature) given by a truncated $\kappa$ distribution is similar to that given by a Maxwellian and only starts to increase for higher truncation velocities. We deduce a simple formula to convert the Maxwellian to $\kappa$ pressure or temperature. We apply this result to the Voyager-2 observations in the heliosheath.

The propagation of charged cosmic rays through the Galactic environment influences all aspects of the observation at Earth. Energy spectrum, composition and arrival directions are changed due to deflections in magnetic fields and interactions with the interstellar medium. Today the transport is simulated with different simulation methods either based on the solution of a transport equation (multi-particle picture) or a solution of an equation of motion (single-particle picture). We developed a new module for the publicly available propagation software CRPropa 3.1, where we implemented an algorithm to solve the transport equation using stochastic differential equations. This technique allows us to use a diffusion tensor which is anisotropic with respect to an arbitrary magnetic background field. The source code of CRPropa is written in C++ with python steering via SWIG which makes it easy to use and computationally fast. In this paper, we present the new low-energy propagation code together with validation procedures that are developed to proof the accuracy of the new implementation. Furthermore, we show first examples of the cosmic ray density evolution, which depends strongly on the ratio of the parallel $\kappa_\parallel$ and perpendicular $\kappa_\perp$ diffusion coefficients. This dependency is systematically examined as well the influence of the particle rigidity on the diffusion process.

This letter presents the results of an advanced parametrization of the solar wind electron temperature anisotropy and the instabilities resulting from the interplay of the (bi-)Maxwellian core and (bi-)Kappa halo populations in the slow solar wind. A large set of observational data (from the Ulysses, Helios and Cluster missions) is used to parametrize these components and establish their correlations. The instabilities are significantly stimulated in the presence of suprathermals, and the instability thresholds shape the limits of the temperature anisotropy for both the core and halo populations re-stating the incontestable role that the selfgenerated instabilities can play in constraining the electron anisotropy. These results confirm a particular implication of the suprathermal electrons which are less dense but hotter than thermal electrons.

A previously published analytical magnetohydrodynamic model for the local interstellar magnetic field in the vicinity of the heliopause (R\"oken et al. 2015) is extended from incompressible to compressible, yet predominantly subsonic flow, considering both isothermal and adiabatic equations of state. Exact expressions and suitable approximations for the density and the flow velocity are derived and discussed. In addition to the stationary induction equation, these expressions also satisfy the momentum balance equation along stream lines. The practical usefulness of the corresponding, still exact analytical magnetic field solution is assessed by comparing it quantitatively to results from a fully self-consistent magnetohydrodynamic simulation of the interstellar magnetic field draping around the heliopause.

Context. Recent studies on Kappa distribution functions invoked in space plasma applications have emphasized two alternative approaches which may assume the temperature parameter either dependent or independent of the power-index $\kappa$. Each of them can obtain justification in different scenarios involving Kappa-distributed plasmas, but direct evidences supporting any of these two alternatives with measurements from laboratory or natural plasmas are not available yet. Aims. This paper aims to provide more facts on this intriguing issue from direct fitting measurements of suprathermal electron populations present in the solar wind, as well as from their destabilizing effects predicted by these two alternating approaches. Methods. Two fitting models are contrasted, namely, the global Kappa and the dual Maxwellian-Kappa models, which are currently invoked in theory and observations. The destabilizing effects of suprathermal electrons are characterized on the basis of a kinetic approach which accounts for the microscopic details of the velocity distribution. Results. In order to be relevant, the model is chosen to accurately reproduce the observed distributions and this is achieved by a dual Maxwellian-Kappa distribution function. A statistical survey indicates a $\kappa$-dependent temperature of the suprathermal (halo) electrons for any heliocentric distance. Only for this approach the instabilities driven by the temperature anisotropy are found to be systematically stimulated by the abundance of suprathermal populations, i.e., lowering the values of $\kappa$-index.

This paper summarizes the results obtained by the team "Heliosheath Processes and the Structure of the Heliopause: Modeling Energetic Particles, Cosmic Rays, and Magnetic Fields" supported by the International Space Science Institute in Bern, Switzerland. We focus on the physical processes occurring in the outer heliosphere, especially at its boundary called the heliopause (HP), and in the LISM. The importance of magnetic field, charge exchange between atoms and ions, and solar cycle on the heliopause topology and observed heliocentric distances to different heliospheric discontinuities are discussed. It is shown that time-dependent boundary conditions are necessary to describe the heliospheric asymmetries detected by the Voyager spacecraft. We also discuss the structure of the HP, especially due to its instability and magnetic reconnection. It is demonstrated that the Rayleigh-Taylor instability of the nose of the HP creates consecutive layers of the interstellar and heliospheric plasma which are magnetically connected to different sources. This may be a possible explanation of abrupt changes in the galactic and anomalous cosmic ray fluxes observed by Voyager 1 when it was crossing the HP structure for a period of about one month in the summer of 2012. This paper also discusses the plausibility of fitting simulation results to a number of observational data sets obtained by in situ and remote measurements. The distribution of magnetic field in the vicinity of the HP is discussed in the context of Voyager measurements. We discuss the transport of energetic particles in the inner and outer heliosheath, concentrating on the anisotropic spatial diffusion diffusion tensor and the pitch-angle dependence of perpendicular diffusion and demonstrate that the latter can explain the observed pitch-angle anisotropies of both the anomalous and galactic cosmic rays in the outer heliosheath.

The study aims at a derivation of generalized \RH relations, especially that for the entropy, for the case of different upstream/downstream polytropic indices and their implications. We discuss the solar/stellar wind interaction with the interstellar medium for different polytropic indices and concentrate on the case when the polytropic index changes across hydrodynamical shocks. We use first a numerical mono-fluid approach with constant polytropic index in the entire integration region to show the influence of the polytropic index on the thickness of the helio-/astrosheath and on the compression ratio. Second, the Rankine-Hugonoit relations for a polytropic index changing across a shock are derived analytically, particularly including a new form of the entropy condition. In application to the/an helio-/astrosphere, we find that the size of the helio-/astrosheath as function of the polytropic index decreases in a mono-fluid model for indices less than $\gamma=5/3$ and increases for higher ones and vice versa for the compression ratio. Furthermore, we demonstrate that changing polytropic indices across a shock are physically allowed only for sufficiently high Mach numbers and that in the hypersonic limit the compression ratio depends only on the downstream polytropic index, while the ratios of the temperature and pressure as well as the entropy difference depend on both, the upstream and downstream polytropic indices.

We extend a two-component model for the evolution of fluctuations in the solar wind plasma so that it is fully three-dimensional (3D) and also coupled self-consistently to the large-scale magnetohydrodynamic (MHD) equations describing the background solar wind. The two classes of fluctuations considered are a high-frequency parallel-propagating wave-like piece and a low-frequency quasi-two-dimensional component. For both components, the nonlinear dynamics is dominanted by quasi-perpendicular spectral cascades of energy. Driving of the fluctuations, by, for example, velocity shear and pickup ions, is included. Numerical solutions to the new model are obtained using the Cronos framework, and validated against previous simpler models. Comparing results from the new model with spacecraft measurements, we find improved agreement relative to earlier models that employ prescribed background solar wind fields. Finally, the new results for the wave-like and quasi-two-dimensional fluctuations are used to calculate ab initio diffusion mean free paths and drift lengthscales for the transport of cosmic rays in the turbulent solar wind.

Huge astrospheres or stellar wind bubbles influence the propagation of cosmic rays at energies up to the TeV range and can act as small-scale sinks decreasing the cosmic ray flux. We model such a sink (in 2D) by a sphere of radius 10\,pc embedded within a sphere of a radius of 1\,kpc. The cosmic ray flux is calculated by means of backward stochastic differential equations from an observer, which is located at $r_{0}$, to the outer boundary. It turns out that such small-scale sinks can influence the cosmic ray flux at the observer's location by a few permille (i.e\ a few 0.1\%), which is in the range of the observations by IceCube, Milagro and other large area telescopes.

The generally accepted representation of $\kappa$-distributions in space plasma physics allows for two different alternatives, namely assuming either the temperature or the thermal velocity to be $\kappa$-independent. The present paper aims to clarify the issue concerning which of the two possible choices and the related physical interpretation is the correct one. A quantitative comparison of the consequences of the use of both distributions for specific physical systems leads to their correct interpretation. It is found that both alternatives can be realized, but are valid for principally different physical systems. The investigation demonstrates that, before employing one of the two alternatives, one should be conscious about the nature of the physical system one intends to describe, otherwise one would possibly obtain unphysical results.

Both physical arguments and simulations of the global heliosphere indicate that the tailward heliopause is flattened considerably in the direction perpendicular to both the incoming flow and the large-scale interstellar magnetic field. Despite this fact, all of the existing global analytical models of the outer heliosheath's magnetic field assume a circular cross section of the heliotail. To eliminate this inconsistency, we introduce a mathematical procedure by which any analytically or numerically given magnetic field can be deformed in such a way that the cross sections along the heliotail axis attain freely prescribed, spatially dependent values for their total area and aspect ratio. The distorting transformation of this method honors both the solenoidality condition and the stationary induction equation with respect to an accompanying flow field, provided that both constraints were already satisfied for the original magnetic and flow fields prior to the transformation. In order to obtain realistic values for the above parameters, we present the first quantitative analysis of the heliotail's overall distortion as seen in state-of-the-art three-dimensional hybrid MHD-kinetic simulations.

The interaction between a supersonic stellar wind and a (super-)sonic interstellar wind has recently been viewed with new interest. We here first give an overview of the modeling, which includes the heliosphere as an example of a special astrosphere. Then we concentrate on the shock structures of fluid models, especially of hydrodynamic (HD) models. More involved models taking into account radiation transfer and magnetic fields are briefly sketched. Even the relatively simple HD models show a rich shock structure, which might be observable in some objects. We employ a single fluid model to study these complex shock structures, and compare the results obtained including heating and cooling with results obtained without these effects. Furthermore, we show that in the hypersonic case valuable information of the shock structure can be obtained from the Rankine-Hugoniot equations. We solved the Euler equations for the single fluid case and also for a case including cooling and heating. We also discuss the analytical Rankine-Hugoniot relations and their relevance to observations. We show that the only obtainable length scale is the termination shock distance. Moreover, the so-called thin shell approximation is usually not valid. We present the shock structure in the model that includes heating and cooling, which differs remarkably from that of a single fluid scenario in the region of the shocked interstellar medium. We find that the heating and cooling is mainly important in this region and is negligible in the regions dominated by the stellar wind beyond an inner boundary.

The Energetic Neutral Atom (ENA) full-sky maps obtained with the Interstellar Boundary Explorer (IBEX) show an unexpected bright narrow band of increased intensity. This so-called ENA ribbon results from charge exchange of interstellar neutral atoms with protons in the outer heliosphere or beyond. Amongst other hypotheses it has been argued that this ribbon may be related to a neutral density enhancement, or H-wave, in the local interstellar medium. Here we quantitatively demonstrate, on the basis of an analytical model of the principal large-scale heliospheric structure, that this scenario for the ribbon formation leads to results that are fully consistent with the observed location of the ribbon in the full-sky maps at all energies detected with high-energy sensor IBEX-Hi.

In-situ spacecraft observations recently suggested that the transport of energetic particles accelerated at heliospheric shocks can be anomalous, i.e. the mean square displacement can grow non-linearly in time. In particular, a new analysis technique has permitted the study of particle transport properties from energetic particle time profiles upstream of interplanetary shocks. Indeed, the time/spatial power laws of the differential intensity upstream of several shocks are indicative of superdiffusion. A complete determination of the key parameters of superdiffusive transport comprises the power-law index, the superdiffusion coefficient, the related transition scale at which the energetic particle profiles turn to decay as power laws, and the energy spectral index of the shock accelerated particles. Assuming large-scale spatial homogeneity of the background plasma, the power-law behaviour can been derived from both a (microscopic) propagator formalism and a (macroscopic) fractional transport equation. We compare the two approaches and find a relation between the diffusion coefficients used in the two formalisms. Based on the assumption of superdiffusive transport, we quantitatively derive these parameters by studying energetic particle profiles observed by the Ulysses and Voyager 2 spacecraft upstream of shocks in the heliosphere, for which a superdiffusive particle transport has previously been observed. Further, we have jointly studied the electron energy spectra, comparing the values of the spectral indices observed with those predicted by the standard diffusive shock acceleration theory and by a model based on superdiffusive transport. For a number of interplanetary shocks and for the solar wind termination shock, for the first time we obtain the anomalous diffusion constants and the scale at which the probability of particle free paths changes to a power-law...

We present the implementation of turbulence transport equations in addition to the Reynolds-averaged MHD equations within the Cronos framework. The model is validated by comparisons with earlier findings before it is extended to be applicable to regions in the solar wind that are not highly super-Alfv\'enic. We find that the respective additional terms result in absolute normalized cross-helicity to decline more slowly, while a proper implementation of the mixing terms can even lead to increased cross-helicities in the inner heliosphere. The model extension allows to place the inner boundary of the simulations closer to the Sun, where we choose its location at 0.1 AU for future application to the Wang-Sheeley-Arge model. Here, we concentrate on effects on the turbulence evolution for transient events by injecting a coronal mass ejection (CME). We find that the steep gradients and shocks associated with these structures result in enhanced turbulence levels and reduced cross-helicity. Our results can now be used straightforwardly for studying the transport of charged energetic particles, where the elements of the diffusion tensor can now benefit from the self-consistently computed solar wind turbulence. Furthermore, we find that there is no strong back-reaction of the turbulence on the large-scale flow so that CME studies concentrating on the latter need not be extended to include turbulence transport effects.

An analytical representation of the interstellar magnetic field in the vicinity of the heliosphere is derived. The three-dimensional field structure close to the heliopause is calculated as a solution of the induction equation under the assumption that it is frozen into a prescribed plasma flow resembling the characteristic interaction of the solar wind with the local interstellar medium. The usefulness of this analytical solution as an approximation to self-consistent magnetic field configurations obtained numerically from the full MHD equations is illustrated by quantitative comparisons.

Cosmic rays passing through large astrospheres can be efficiently cooled inside these "cavities" in the interstellar medium. Moreover, the energy spectra of these energetic particles are already modulated in front of the astrospherical bow shocks. We study the cosmic ray flux in and around lambda Cephei as an example for an astrosphere. The large-scale plasma flow is modeled hydrodynamically with radiative cooling. We studied the cosmic ray flux in a stellar wind cavity using a transport model based on stochastic differential equations. The required parameters, most importantly, the elements of the diffusion tensor, are based on the heliospheric parameters. The magnetic field required for the diffusion coefficients is calculated kinematically. We discuss the transport in an astrospheric scenario with varying parameters for the transport coefficients. We show that large stellar wind cavities can act as sinks for the galactic cosmic ray flux and thus can give rise to small-scale anisotropies in the direction to the observer. Small-scale cosmic ray anisotropies can naturally be explained by the modulation of cosmic ray spectra in huge stellar wind cavities.

The transport of energetic particles such as Cosmic Rays is governed by the properties of the plasma being traversed. While these properties are rather poorly known for galactic and interstellar plasmas due to the lack of in situ measurements, the heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the 3D heliospheric transport of energetic particles are structures such as corotating interaction regions (CIRs), which, due to strongly enhanced magnetic field strengths, turbulence, and associated shocks, can act as diffusion barriers on the one hand, but also as accelerators of low energy CRs on the other hand as well. In a two-fold series of papers we investigate these effects by modeling inner-heliospheric solar wind conditions with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve as an input to a transport code employing a stochastic differential equation (SDE) approach (second paper). In this first paper we present results from 3D MHD simulations with our code CRONOS: for validation purposes we use analytic boundary conditions and compare with similar work by Pizzo. For a more realistic modeling of solar wind conditions, boundary conditions derived from synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where the potential field modeling is performed with a finite-difference approach (FDIPS) in contrast to the traditional spherical harmonics expansion often utilized in the WSA model. Our results are validated by comparing with multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses) regions.

In the heliosphere, especially in the inner heliosheath, mass-, momentum-, and energy loading induced by the ionization of neutral interstellar species plays an important, but for some species, especially Helium, an underestimated role. We discuss the implementation of charge exchange and electron impact processes for interstellar neutral Hydrogen and Helium and their implications for further modeling. Especially, we emphasize the importance of electron impact and a more sophisticated numerical treatment of the charge exchange reactions. Moreover, we discuss the non-resonant charge exchange effects. The rate coefficients are discussed and the influence of the cross-sections in the (M)HD equations for different reactions are revised as well as their representation in the collision integrals. Electron impact is in some regions of the heliosphere, particularly in the heliotail, more effective than charge exchange, and the ionization of neutral interstellar Helium contributes about 40% to the mass- and momentum loading in the inner heliosheath. The charge exchange cross-sections need to be modeled with higher accuracy, especially in view of the latest developments in their description. The ionization of Helium and electron impact ionization of Hydrogen needs to be taken into account for the modeling of the heliosheath and, in general, astrosheaths. Moreover, the charge exchange cross-sections need to be handled in a more sophisticated way, either by developing better analytic approximations or by solving the collision integrals numerically.

Recently it has been discussed whether a bow shock ahead of the heliospheric stagnation region does exist or not. This discussion was triggered by measurements indicating that the Alfv\'en speed and that of fast magnetosonic waves are higher than the flow speed of the local interstellar medium (LISM) relative to the heliosphere and resulted in the conclusion that there might exist either a bow wave or a slow magnetosonic shock. We demonstrate here that including the He$^{+}$ component of the LISM yields both an Alfv\'en and fast magnetosonic wave speed lower than the LISM flow speed. Consequently, the scenario of a bow shock in front of the heliosphere as modelled in numerous simulations of the interaction of the solar wind with the LISM remains valid.

Context. Spacecraft observations have motivated the need for a refined description of the phase-space distribution function. Of particular importance is the pitch-angle diffusion coefficient that occurs in the Fokker-Planck transport equation. Aims. Simulations and analytical test-particle theories are compared to verify the diffusion description of particle transport, which does not allow for non-Markovian behavior. Methods. A Monte-Carlo simulation code was used to trace the trajectories of test particles moving in turbulent magnetic fields. From the ensemble average, the pitch-angle Fokker-Planck coefficient is obtained via the mean square displacement. Results. It is shown that, while excellent agreement with analytical theories can be obtained for slab turbulence, considerable deviations are found for isotropic turbulence. In addition, all Fokker-Planck coefficients tend to zero for high time values.

The transport environment for particles in the heliosphere, e.g. galactic cosmic rays (GCRs) and MeV electrons (including those originating from Jupiters magnetosphere), is defined by the solar wind flow and the structure of the embedded heliospheric magnetic field. Solar wind structures, such as co-rotating interaction regions (CIR), can result in periodically modulation of both particles species. A detailed analysis of this recurrent Jovian electron events and galactic cosmic ray decreases measured by SOHO EPHIN is presented here, showing clearly a change of phase between both phenomena during the cause of the years 2007 and 2008. This effect can be explained by the change of difference in heliolongitude between the Earth and Jupiter, which is of central importance for the propagation of Jovian electrons. Furthermore, the data can be ordered such that the 27-day Jovian electron variation vanishes in the sector which does not connect the Earth with Jupiter magnetically using observed solar wind speeds.

Maps of the radial magnetic field at a heliocentric distance of ten solar radii are used as boundary conditions in the MHD code CRONOS to simulate a 3D inner-heliospheric solar wind emanating from the rotating Sun out to 1 AU. The input data for the magnetic field are the result of solar surface flux transport modelling using observational data of sunspot groups coupled with a current sheet source surface model. Amongst several advancements, this allows for higher angular resolution than that of comparable observational data from synoptic magnetograms. The required initial conditions for the other MHD quantities are obtained following an empirical approach using an inverse relation between flux tube expansion and radial solar wind speed. The computations are performed for representative solar minimum and maximum conditions, and the corresponding state of the solar wind up to the Earths orbit is obtained. After a successful comparison of the latter with observational data, they can be used to drive outer-heliospheric models.

Galactic transport models for cosmic rays involve the diffusive motion of these particles in the interstellar medium. Due to the large-scale structured galactic magnetic field this diffusion is anisotropic with respect to the local field direction. We included this transport effect along with continuous loss processes in a quantitative model of galactic propagation for cosmic ray protons which is based on stochastic differential equations. We calculated energy spectra at different positions along the Sun's galactic orbit and compared them to the isotropic diffusion case. The results show that a larger amplitude of variation as well as different spectral shapes are obtained in the introduced anisotropic diffusion scenario and emphasize the need for accurate galactic magnetic field models.

A recent third-order, essentially non-oscillatory central scheme to advance the equations of single-fluid magnetohydrodynamics (MHD) in time has been implemented into a new numerical code. This code operates on a 3-D Cartesian, non-staggered grid, and is able to handle shock-like gradients without producing spurious oscillations. To demonstrate the suitability of our code for the simulation of coronal mass ejections (CMEs) and similar heliospheric transients, we present selected results from test cases and perform studies of the solar wind expansion during phases of minimum solar activity. We can demonstrate convergence of the system into a stable Parker-like steady state for both hydrodynamic and MHD winds. The model is subsequently applied to expansion studies of CME-like plasma bubbles, and their evolution is monitored until a stationary state similar to the initial one is achieved. In spite of the model's (current) simplicity, we can confirm the CME's nearly self-similar evolution close to the Sun, thus highlighting the importance of detailed modelling especially at small heliospheric radii. Additionally, alternative methods to implement boundary conditions at the coronal base, as well as strategies to ensure a solenoidal magnetic field, are discussed and evaluated.

New aspects of the slow solar wind turbulent heating and acceleration are investigated. A physical meaning of the lower boundary of the Alfv\'en wave turbulent spectra in the solar atmosphere and the solar wind is studied and the significance of this natural parameter is demonstrated. Via an analytical and quantitative treatment of the problem we show that a truncation of the wave spectra from the lower frequency side, which is a consequence of the solar magnetic field structure and its cyclic changes, results in a significant reduction of the heat production and acceleration rates. An appropriate analysis is presented regarding the link of the considered problem with existing observational data and slow solar wind initiation scenarios.

The spatial diffusion of cosmic rays in turbulent magnetic fields can, in the most general case, be fully anisotropic, i.e. one has to distinguish three diffusion axes in a local, field-aligned frame. We reexamine the transformation for the diffusion tensor from this local to a global frame, in which the Parker transport equation for energetic particles is usually formulated and solved. Particularly, we generalize the transformation formulas to allow for an explicit choice of two principal local perpendicular diffusion axes. This generalization includes the 'traditional' diffusion tensor in the special case of isotropic perpendicular diffusion. For the local frame, we motivate the choice of the Frenet-Serret trihedron which is related to the intrinsic magnetic field geometry. We directly compare the old and the new tensor elements for two heliospheric magnetic field configurations, namely the hybrid Fisk and the Parker field. Subsequently, we examine the significance of the different formulations for the diffusion tensor in a standard 3D model for the modulation of galactic protons. For this we utilize a numerical code to evaluate a system of stochastic differential equations equivalent to the Parker transport equation and present the resulting modulated spectra. The computed differential fluxes based on the new tensor formulation deviate from those obtained with the 'traditional' one (only valid for isotropic perpendicular diffusion) by up to 60% for energies below a few hundred MeV depending on heliocentric distance.

Magnetic clouds (MCs) are "magnetized plasma clouds" moving in the solar wind. MCs transport magnetic flux and helicity away from the Sun. These structures are not stationary but feature temporal evolution as they propagate in the solar wind. Simplified analytical models are frequently used for the description of MCs, and fit certain observational data well. The goal of the present study is to investigate numerically the validity of an analytical model which is widely used for the description of MCs, and to determine under which conditions this model's implied assumptions cease to be valid. A numerical approach is applied. Analytical solutions that have been derived in previous studies are implemented in a \textbf{3-D magnetohydrodynamic} simulation code as initial conditions. Initially, the analytical model represents the main observational features of the MCs. However, these characteristics prevail only if the structure moves with a velocity close to the velocity of the background flow. In this case an MC's evolution can quite accurately be described using an analytic, self-similar approach. The dynamics of the magnetic structures which move with a velocity significantly above or below that of the velocity of the solar wind is investigated in detail. Besides the standard case in which MCs only expand and propagate in the solar wind, the case of an MC rotating around its axis of symmetry is also considered, and the resulting influence on the MC's dynamics is studied.

In this paper we show results of numerical simulations for the turbulence in the interstellar medium. These results were obtained using a Riemann solver-free numerical scheme for high-Mach number hyperbolic equations. Here we especially concentrate on the physical properties of the ISM. That is, we do not present turbulence simulations trimmed to be applicable to the interstellar medium. The simulations are rather based on physical estimates for the relevant parameters of the interstellar gas. Applying our code to simulate the turbulent plasma motion within a typical interstellar molecular cloud, we investigate the influence of different equations of state (isothermal and adiabatic) on the statistical properties of the resulting turbulent structures. We find slightly different density power spectra and dispersion maps, while both cases yield qualitatively similar dissipative structures, and exhibit a departure from the classical Kolmogorov case towards a scaling described by the She-Leveque model. Solving the full energy equation with realistic heating/cooling terms appropriate for the diffuse interstellar gas, we are able to reproduce a realistic two-phase distribution of cold and warm plasma. When extracting maps of polarised intensity from our simulation data, we find encouraging similarity to actual observations. Finally, we compare the actual magnetic field strength of our simulations to its value inferred from the rotation measure. We find these to be systematically different by a factor of about 1.5, thus highlighting the often underestimated influence of varying line-of-sight particle densities on the magnetic field strength derived from observed rotation measures.