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.

We study the behaviour of shocked wind-electrons leaving wind-driving stars after undergoing the outer wind termination shock. As an example, we describe the evolution of the keV-energetic electron distribution function downstream of the heliospheric termination shock. We start from a kinetic transport equation in the bulk frame of the heliosheath plasma flow taking into account shock-induced electron injection, convective changes, cooling processes, and whistler wave-induced energy diffusion. From this equation we proceed to an associated pressure moment of the electron distribution function arriving at a corresponding pressure transport equation which describes the evolution of the electron pressure in the bulk frame of the plasma along the plasma flow lines. We assume that the local distribution function, in view of the prevailing non-LTE conditions, is represented by a local kappa function with local kappa parameters that vary with the streamline coordinate $s$ downstream of the solar wind termination shock. We obtain the solution for the electron pressure as a function of the streamline coordinate $s$ from the pressure transport equation and demonstrate that, connected with this pressure, one obtains an expression for the entropy of the electron fluid which can also be derived as a streamline function. We show that the heliosheath electron fluid can essentially be characterized as an isobaric and isentropic flow. These results allow to generally conclude that astrotail plasma flows are characterized as such flows.

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 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 consider a fast magnetosonic multifluid shock as a representation of the solar-wind termination shock. We assume the action of the transition happens in a three-step process: In the first step, the upstream supersonic solar-wind plasma is subject to a strong electric field that flashes up on a small distance scale $\Delta z\simeq U_1/ \Omega _{\mathrm e}$ (first part of the transition layer), where $\Omega_{\mathrm e}$ is the electron gyro-frequency and $U_1$ is the upstream speed. This electric field both decelerates the supersonic ion flow and accelerates the electrons up to high velocities. In this part of the transition region, the electric forces connected with the deceleration of the ion flow strongly dominate over the Lorentz forces. We, therefore, call this part the demagnetization region. In the second phase, Lorentz forces due to convected magnetic fields compete with the electric field, and the highly anisotropic and energetic electron distribution function is converted into a shell distribution with energetic shell electrons storing about 3/4 of the upstream ion kinetic energy. In the third phase, the plasma particles thermalize due to the relaxation of free energy by plasma instabilities. The first part of the transition region opens up a new thermodynamic degree of freedom never before taken into account for the electrons, since the electrons are usually considered to be enslaved to follow the behavior of the protons in all velocity moments like density, bulk velocity, and temperature. We show that electrons may be the downstream plasma fluid that dominates the downstream plasma pressure.

In the majority of the literature on plasma shock waves, electrons play the role of "ghost particles," since their contribution to mass and momentum flows is negligible, and they have been treated as only taking care of the electric plasma neutrality. In some more recent papers, however, electrons play a new important role in the shock dynamics and thermodynamics, especially at the solar-wind termination shock. They react on the shock electric field in a very specific way, leading to suprathermal nonequilibrium distributions of the downstream electrons, which can be represented by a kappa distribution function. In this paper, we discuss why this anticipated hot electron population has not been seen by the plasma detectors of the Voyager spacecraft downstream of the solar-wind termination shock. We show that hot nonequilibrium electrons induce a strong negative electric charge-up of any spacecraft cruising through this downstream plasma environment. This charge reduces electron fluxes at the spacecraft detectors to nondetectable intensities. Furthermore, we show that the Debye length $\lambda _{\mathrm D}^{\kappa}$ grows to values of about $\lambda _{\mathrm D}^{\kappa}/\lambda _{\mathrm D}\simeq 10^{6}$ compared to the classical value $\lambda _{\mathrm D}$ in this hot-electron environment. This unusual condition allows for the propagation of a certain type of electrostatic plasma waves that, at very large wavelengths, allow us to determine the effective temperature of the suprathermal electrons directly by means of the phase velocity of these waves. At moderate wavelengths, the electron-acoustic dispersion relation leads to nonpropagating oscillations with the ion-plasma frequency $\omega _{\mathrm p}$ , instead of the traditional electron plasma frequency.

Ambient plasma electrons undergo strong heating in regions associated with compressive traveling interplanetary solar-wind bulk-velocity jumps due to their specific interactions with the jump-inherent electric fields. After thermalization of this energy gain per shock passage through the operation of the Buneman instability, strong electron heating occurs that substantially influences the radial electron temperature profile. We describe the reduction of the jump amplitude due to energy expended by the traveling jump structure. We consider three effects; namely energy loss due to heating of electrons, energy loss due to work done against the pick-up-ion pressure gradient, and an energy gain due to nonlinear jump steepening. Taking these effects into account, we show that the decrease in jump amplitude with solar distance is more pronounced when the initial jump amplitude is higher in the inner solar system. Independent of the initial jump amplitude, it eventually decreases with increasing distance to a value of the order of $\Delta U/U\simeq 0.1$ at the position of the heliospheric termination shock, where $\Delta U$ is the jump amplitude, and $U$ is the average solar-wind bulk velocity.The electron temperature, on the other hand, is strongly correlated with the initial jump amplitude, leading to electron temperatures between 6000 K and 20 000 K at distances beyond 50 AU. We compare our results with in-situ measurements of the electron-core temperature from the Ulysses spacecraft in the plane of the ecliptic for $1.5\, \mathrm{AU}\leq r \leq 5\,\mathrm{AU}$, where $r$ is the distance from the Sun. We find a very good agreement between our results and these observations, which corroborates our extrapolated predictions beyond $r=5\,\mathrm{AU}$.

Understanding the outer heliospheric interface is a major challenge, especially in the light of recent observations by the IBEX and Voyager missions. We present further details on a new required theoretical correction that has been identified as substantial in a recent paper, the so-called source depletion effect. These results complement and expand earlier calculations of transit-time delays by presenting global skymaps of Energetic Neutral Atoms (ENAs) calculated with the depletion correction, comparing them with skymaps calculated without these corrections. Our results demonstrate that the depletion correction is critical for interpreting IBEX-Hi ENA fluxes generated in the inner heliosheath, and that any attempt to reconstruct the shape of the heliospheric boundary needs to include the depletion correction, unless arriving at considerably erroneous results.

In a series of earlier papers, we have developed expressions for ion and electron velocity distribution functions and theirvelocity moments at the passage over the solar wind termination shock. As we have shown there, with introduction of appropriate particle invariants and the use of Liouville`s theorem one can get explicit solutions for the resulting total downstream pressure adding up from partial pressure contributions of solar wind protons, solar wind electrons and pick-up protons. These expressions deliver in a first step the main contributions to the total plasma pressure in the downstream plasma flow and consistently determine the shock compression ratio. Here now we start out from these individual fluid pressures downstream of the shock and thereafter evaluate for the first time the shock-induced entropy production of the different fluids, when they are passing over the shock to the downstream side. As is shown here, the resulting ion entropy production substantially deviates from earlier calculations using a pseudo-polytropic reaction of the ions to the shock compression, with polytropies selected to describe fluid-specific reactions at the shock passage similar to those seen by the VOYAGERs. From these latter models ion entropy jumps are derived that depend on the pick-up ion abundance, while our calculations, to the opposite, deliver an abundance-independent ion entropy production which only depends on the shock compression ratio and the tilt angle between the upstream magnetic field and the shock surface normal. We also do show here that only when including the strongly heated electrons into the entropy balance one then arrives at the total entropy production that just fulfills the thermodynamically permitted limit.

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.

A stellar wind passing through the reverse shock is deflected into the astrospheric tail and leaves the stellar system either as a sub-Alfvenic or as a super-Alfvenic tail flow. An example is our own heliosphere and its heliotail. We present an analytical method of calculating stationary, incompressible, and field-aligned plasma flows in the astrotail of a star. We present a recipe for constructing an astrosphere with the help of only a few parameters, like the inner Alfven Mach number and the outer Alfven Mach number, the magnetic field strength within and outside the stellar wind cavity, and the distribution of singular points of the magnetic field within these flows. Within the framework of a one-fluid approximation, it is possible to obtain solutions of the MHD equations for stationary flows from corresponding static MHD equilibria, by using noncanonical mappings of the canonical variables. The canonical variables are the Euler potentials of the magnetic field of magnetohydrostatic equilibria. Thus we start from static equilibria determined by the distribution of magnetic neutral points, and assume that the Alfven Mach number for the corresponding stationary equilibria is finite. The topological structure determines the geometrical structure of the interstellar gas - stellar wind interface. Additional boundary conditions like the outer magnetic field and the jump of the magnetic field across the astropause allow determination of the noncanonical transformations. This delivers the strength of the magnetic field at every point in the astrotail region beyond the reverse shock. The mathematical technique for describing such a scenario is applied to astrospheres in general, but is also relevant for the heliosphere. It shows the restrictions of the outer and the inner magnetic field strength in comparison with the corresponding Alfven Mach numbers in the case of subalfvenic flows.

We investigated the impact of three different modifications of Newtonian gravity on motions of Keplerian objects within the Solar System. These objects are located at distances of the order of the distance to the Oort cloud. With these three modifications we took into account a heliocentric Dark-Matter halo as was indicated by Diemand et al, Modified Newtonian Dynamics (MOND) and a vacuum-induced force due to a locally negative cosmological constant $\Lambda_-$ derived by Fahr & Siewert. In gravitationally bound systems it turns out that all three modifications deliver the same qualitative results: Initially circular orbits for the pure Newtonian case are forced to convert into ellipses with perihelion migrations. The quantitative consideration, however, of the orbital parameters showed strong differences between MOND on the one side, and Dark-Matter and $\Lambda_-$ effects on the other side.

Though pick-up ions (PUIs) are a well known phenomenon in the inner heliosphere their phase-space distribution nevertheless is a theoretically unsettled problem. Especially the question of how pick-up ions form their suprathermal tails, extending to far above their injection energies, still now is unsatistactorily answered. Though Fermi-2 velocity diffusion theories have revealed that such tails are populated, they nevertheless show that resulting population densities are much less than seen in observations showing power-laws with a velocity index of "-5". We first investigate here, whether or not observationally suggested power-laws can be the result of a quasi-equilibrium state between suprathermal ions and magnetohydrodynamic turbulences in energy exchange with eachother. We demonstrate that such an equilibrium cannot be established, since it would require too high pick-up ion pressures enforcing a shock-free deceleration of the solar wind. We furthermore show that Fermi-2 type energy diffusion in the outer heliosphere is too inefficient to determine the shape of the distribution function there. As we can show, however, power-laws beyond the injection threshold can be established, if the injection takes place at higher energies of the order of 100 keV. As we demonstrate here, such an injection is connected with modulated anomalous cosmic ray (ACR) particles at the lower end of their spectrum when they again start being convected outwards with the solar wind. Therefore, we refer to these particles as ACR-PUIs. In our quantitative calculation of the pick-up ion spectrum resulting under such conditions we in fact find again power-laws, however with a velocity power index of "-4" and fairly distance-independent spectral intensities.

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.

In this article we want to answer the cosmologically relevant question what, with some good semantic and physical reason, could be called the mass of an infinitely extended, homogeneously matter-filled and expanding universe. To answer this question we produce a space-like sum of instantaneous cosmic energy depositions surrounding equally each spacepoint in the homogeneous universe. We calculate the added-up instantaneous cosmic energy per volume around an arbitrary space point in the expanding universe. To carry out this sum we use as basic metrics an analogy to the inner Schwarzschild metric applied to stars, but this time applied to the spacepoint-related universe. It is then shown that this leads to the added-up proper energy within a sphere of a finite outer critical radius defining the point-related infinity. As a surprise this radius turns out to be reciprocal to the square root of the prevailing average cosmic energy density. The equivalent mass of the universe can then also be calculated and, by the expression which is obtained here, shows a scaling with this critical radius of this universe, a virtue of the universe which was already often called for in earlier works by E.Mach, H.Thirring and F.Hoyle and others. This radius on the other hand can be shown to be nearly equal to the Schwarzschild radius of the so-defined mass of the universe.

Cosmological consequences of a strictly valid total energy conservation for the whole universe are investigated in this paper. Interestingly enough as one consequence of ergodically behaving universes very specific scaling laws with the diameter R of the universe can be derived for relevant cosmic quantities. Especially the 1/R^2- scaling of mass - and vacuum energy - density then automatically leads to a vanishing cosmic curvature parameter k=0 and also reveals, that for such universes no horizon problem occurs. In addition the longstanding problem of observationally indicated very low cosmic vacuum energies in contrast to the very large quantumfield estimates is easily solved when the vacuum energy density decay with 1/R^2 is taken into account reconciling presently observed vacuum energy density values with those from theoretical expectations. We also suggest why the mass of the universe can permanently increase and can in fact grow up from a Planck mass as a pure vacuum fluctuation.

We present the basic equations for stationary, incompressible resistive MHD flows in two dimensions. This leads to a system of differential equations for two flux functions, one elliptic partial differential equation (Grad-Shafranov-like) for the magnetic flux function and one for the stream function of the flow. In these equations two potentials appear: one potential is a generalized pressure. The second potential couples the magnetic and the flow shear components of the system. With the restriction to flux or at least line conserving flows one has to solve a modified Ohm's law. For the two dimensional case these are two coupled differential equations, which represent the borderline between the resistive but flux conserving (or line conserving) case, and that of reconnective solutions. We discuss some simplified solutions of these equations.

In this MHD-model of the heliosphere, we assume a Parker-type flow, and a Parker-type spiral magnetic field, which is extrapolated further downstream from the termination shock to the heliopause. We raise the question whether the heliopause nose region may be leaky with respect to fields and plasmas due to nonideal plasma dynamics, implying a breakdown of the magnetic barrier. We analyse some simple scenarios to find reconnection rates and circumstances, under which the heliosphere can be an "open" or a "closed" magnetosphere. We do not pretend to offer a complete solution for the heliosphere, on the basis of nonideal MHD theory, but present a prescription to find such a solution on the basis of potential fields including the knowledge of neutral points. As an example we imitate the Parker spiral as a monopole with a superposition of homogeneous asymptotical boundary conditions. We use this toy model for x < -R where R = 100 AU is the distance of the termination shock to describe the situation in the nose region of the heliopause.

We make some remarks on reconnection in plasmas and want to present some calculations related to the problem of finding velocity fields which conserve magnetic flux or at least magnetic field lines. Hereby we start from views and definitions of ideal and non-ideal flows on one hand, and of reconnective and non-reconnective plasma dynamics on the other hand. Our considerations give additional insights into the discussion on violations of the frozen--in field concept which started recently with the papers by Baranov & Fahr (2003a; 2003b). We find a correlation between the nonidealness which is given by a generalized form of the Ohm's law and a general transporting velocity, which is field line conserving.

Various hydrodynamic models of the heliospheric interface have been presented meanwhile, numerically simulating the interaction of the solar wind plasma bubble with the counterstreaming partially ionized interstellar medium. In these model approaches the resulting interface flows are found by the use of hydrodynamic simulation codes trying to consistently describe the dynamic and thermodynamic coupling of the different interacting fluids of protons, H-atoms and pick-up ions. Within such approaches, the fluids are generally expected to be correctly described by the three lowest velocity moments, i.e., by shifted Maxwellians. We shall show that in these approaches the charge-exchange-induced momentum coupling is treated in an unsatisfactory representation valid only at supersonic differential flow speeds. Though this flaw can be removed by an improved coupling term, we shall further demonstrate that the assumption of shifted Maxwellians in some regions of the interface is insufficiently well fulfilled both for H-atoms and protons. Using a Boltzmann-kinetic description of the proton- and H-atom- distribution functions coupled by charge exchange processes we emphasize the fact that non-negligible deviations from shifted Maxwellians are generated in the interface. This has to be taken into account when interpreting inner heliospheric measurements in terms of interstellar parameters.

It is well known that the Solar System is presently moving through a partially ionized local interstellar medium. This gives rise to a counter-flow situation requiring a consistent description of behaviour of the two fluids -- ions and neutral atoms -- which are dynamically coupled by mutual charge exchange processes. Solutions to this problem have been offered in the literature, all relying on the assumption that the proton fluid, even under evidently nonequilibrium conditions, can be expected to stay in a highly-relaxated distribution function given by mono-Maxwellians shifted by the local proton bulk velocity. Here we check the validity of this assumption, calculating on the basis of a Boltzmann-kinetic approach the actually occurring deviations. As we show, especially for low degrees of ionization, $\xi \le 0.3$, both the H-atoms and protons involved do generate in the heliospheric interface clearly pronounced deviations from shifted Maxwellians with asymmetrically shaped distribution functions giving rise to non-convective transport processes and heat conduction flows. Also in the inner heliosheath region and in the heliotail deviations of the proton distribution from the hydrodynamic one must be expected. This sheds new light on the correctness of current calculations of H-atom distribution functions prevailing in the inner heliosphere and also of the Lyman-$\alpha $ absorption features in stellar spectra due to the presence of the hydrogen wall atoms. Deviations from LTE-functions would be even more pronounced in magnetic interfaces, which via CGL-effects cause temperature anisotropies to arise.