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$.

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.

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.

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.

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 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.

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 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.

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.

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.

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.

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.

We present Hubble Space Telescope optical coronagraphic polarization imaging observations of the dusty debris disk HD 61005. The scattered light intensity image and polarization structure reveal a highly inclined disk with a clear asymmetric, swept back component, suggestive of significant interaction with the ambient interstellar medium. The combination of our new data with the published 1.1 micron discovery image shows that the grains are blue scattering with no strong color gradient as a function of radius, implying predominantly sub-micron sized grains. We investigate possible explanations that could account for the observed swept back, asymmetric morphology. Previous work has suggested that HD 61005 may be interacting with a cold, unusually dense interstellar cloud. However, limits on the intervening interstellar gas column density from an optical spectrum of HD 61005 in the Na I D lines render this possibility unlikely. Instead, HD 61005 may be embedded in a more typical warm, low-density cloud that introduces secular perturbations to dust grain orbits. This mechanism can significantly distort the ensemble disk structure within a typical cloud crossing time. For a counterintuitive relative flow direction--parallel to the disk midplane--we find that the structures generated by these distortions can very roughly approximate the HD 61005 morphology. Future observational studies constraining the direction of the relative interstellar medium flow will thus provide an important constraint for future modeling. Independent of the interpretation for HD 61005, we expect that interstellar gas drag likely plays a role in producing asymmetries observed in other debris disk systems, such as HD 15115 and Delta-Velorum.

Modeling of the global heliosphere seeks to investigate the interaction of the solar wind with the partially ionized local interstellar medium. Models that treat neutral hydrogen self-consistently and in great detail, together with the plasma, but that neglect magnetic fields, constitute a sub-category within global heliospheric models. There are several different modeling strategies used for this sub-category in the literature. Differences and commonalities in the modeling results from different strategies are pointed out. Plasma-only models and fully self-consistent models from four research groups, for which the neutral species is modeled with either one, three, or four fluids, or else kinetically, are run with the same boundary parameters and equations. They are compared to each other with respect to the locations of key heliospheric boundary locations and with respect to the neutral hydrogen content throughout the heliosphere. In many respects, the models' predictions are similar. In particular, the locations of the termination shock agree to within 7% in the nose direction and to within 14% in the downwind direction. The nose locations of the heliopause agree to within 5%. The filtration of neutral hydrogen from the interstellar medium into the inner heliosphere, however, is model dependent, as are other neutral results including the hydrogen wall. These differences are closely linked to the strength of the interstellar bow shock. The comparison also underlines that it is critical to include neutral hydrogen into global heliospheric models.