We present a fully analytical, time-dependent leptonic one-zone model that describes a simplified radiation process of multiple interacting ultrarelativistic electron populations, accounting for the flaring of GeV blazars. In this model, several mono-energetic, ultrarelativistic electron populations are successively and instantaneously injected into the emission region, i.e., a magnetized plasmoid propagating along the blazar jet, and subjected to linear, time-independent synchrotron radiative losses, which are caused by a constant magnetic field, and nonlinear, time-dependent synchrotron self-Compton radiative losses in the Thomson limit. Considering a general multiple-injection scenario is, from a physical point of view, more realistic than the usual single-injection scenario invoked in common blazar models, as blazar jets may extend over tens of kiloparsecs and, thus, most likely pick up several particle populations from intermediate clouds. We analytically compute the electron number density by solving a kinetic equation using Laplace transformations and the method of matched asymptotic expansions. Moreover, we explicitly calculate the optically thin synchrotron intensity, the synchrotron self-Compton intensity in the Thomson limit, as well as the associated total fluences. In order to mimic injections of finite duration times and radiative transport, we model flares by sequences of these instantaneous injections, suitably distributed over the entire emission region. Finally, we present a parameter study for the total synchrotron and synchrotron self-Compton fluence spectral energy distributions for a generic three-injection scenario, varying the magnetic field strength, the Doppler factor, and the initial electron energy of the first injection in realistic parameter domains, demonstrating that our model can reproduce the typical broad-band behavior seen in observational data.

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.

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.

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.

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.

Context: Since a few years, signatures of supernova remnants associated with molecular clouds have been detected in gamma rays. Whether these gamma rays are generated by cosmic ray electrons or by cosmic ray protons is usually not known. The detection of hadronic ionization signatures in spatial coincidence with gamma ray signatures can help to unambiguously identify supernova remnants as sources of cosmic ray protons. Methods: In order to calculate hadronic signatures from cosmic ray-induced ionization for an examination of the origin of the observed gamma rays, the transport equation for cosmic ray protons propagating in a molecular cloud, including the relevant momentum loss processes, is solved analytically and the proton flux at any position in the cloud is determined. Results: Since the solution of the transport equation is obtained for arbitrary source functions, it can be used for a variety of supernova remnants. The corresponding theoretical ionization rate, as a function of the penetration depth, is derived and compared to photoinduced ionization profiles in a case study with four supernova remnants associated with molecular clouds. Three of the remnants show a clear dominance of the hadronically induced ionization rate, while for one remnant, X-ray emission seems to dominate by a factor of 10. Conclusions: This is the first derivation of position-dependent profiles for cosmic ray-induced ionization with an analytic solution for arbitrary cosmic ray source spectra. The cosmic ray-induced ionization has to be compared to X-ray ionization for strong X-ray sources. For sources dominated by cosmic ray-induced ionization (e.g., W49B), the ionization profiles can be used in the future to map the spatial structure of hadronic gamma rays and rotation-vibrational lines induced by cosmic ray protons, helping to identify sources of hadronic cosmic rays.

Astrospheres and wind bubbles of massive stars are believed to be sources of cosmic rays with energies $E\lesssim 1\,$TeV. These particles are not directly detectable, but their impact on surrounding matter, in particular ionisation of atomic and molecular hydrogen, can lead to observable signatures. A correlation study of both gamma ray emission, induced by proton-proton interactions of cosmic ray protons with kinetic energies $E_\mathrm{p}\ge 280\,$MeV with ambient hydrogen, and ionisation induced by cosmic ray protons of kinetic energies $E_\mathrm{p}< 280\,$MeV can be performed in order to study potential sources of (sub)TeV cosmic rays.

The proper quantum plasma treatment of the electron gas in degenerate stars such as white dwarfs provides an additional quantum contribution to the electron pressure. The additional pressure term modifies the equation for hydrostatic equilibrium, resulting in the quantum modified Lane-Emden equation for polytropic equation of states. The additional pressure term also modifies the expression for the limiting Chandrasekhar mass of white dwarfs. An approximate solution is derived of the quantum modified Lane-Emden equation for general polytropic indices, and it is demonstrated that the quantum corrections reduce the standard Chandrasekhar mass and enhance the white dwarf radius by negligibly small values only.