Search

The solar wind (SW) and local interstellar medium (LISM) are turbulent media. Their interaction is governed by complex physical processes and creates heliospheric regions with significantly different properties in terms of particle populations, bulk flow and turbulence. Our knowledge of the solar wind turbulence \nature and dynamics mostly relies on near-Earth and near-Sun observations, and has been increasingly improving in recent years due to the availability of a wealth of space missions, including multi-spacecraft missions. In contrast, the properties of turbulence in the outer heliosphere are still not completely understood. In situ observations by Voyager and New Horizons, and remote neutral atom measurements by IBEX strongly suggest that turbulence is one of the critical processes acting at the heliospheric interface. It is intimately connected to charge exchange processes responsible for the production of suprathermal ions and energetic neutral atoms. This paper reviews the observational evidence of turbulence in the distant SW and in the LISM, advances in modeling efforts, and open challenges., Comment: 82 pages, 31 figures. Paper accepted by Space Science Reviews (collection: The Heliosphere in the Local Interstellar Medium: Into the Unknown), and presented at the ISSI Workshop on November 8-12, 2021

Rotationally symmetric bodies with longitudinal cross sections of parabolic shape are frequently used to model astrophysical objects, such as magnetospheres and other blunt objects, immersed in interplanetary or interstellar gas or plasma flows. We discuss a simple formula for the potential flow of an incompressible fluid around an elliptic paraboloid whose axis of symmetry coincides with the direction of incoming flow. Prescribing this flow, we derive an exact analytical solution to the induction equation of ideal magnetohydrodynamics for the case of an initially homogeneous magnetic field of arbitrary orientation being passively advected in this flow. Our solution procedure employs Euler potentials and Cauchy's integral formalism based on the flow's stream function and isochrones. Furthermore, we use a particular renormalization procedure that allows us to generate more general analytical expressions modeling the deformations experienced by arbitrary scalar or vector-valued fields embedded in the flow as they are advected first towards and then past the parabolic obstacle. Finally, both the velocity field and the magnetic field embedded therein are generalized from incompressible to mildly compressible flow, where the associated density distribution is found from Bernoulli's principle., Comment: 20 pages, 6 color figures; submitted to Physics of Fluids

We derive an exact, time-dependent analytical magnetic field solution for the inner heliosheath, which satisfies both the induction equation of ideal magnetohydrodynamics in the limit of infinite electric conductivity and the magnetic divergence constraint. To this end, we assume that the magnetic field is frozen into a plasma flow resembling the characteristic interaction of the solar wind with the local interstellar medium. Furthermore, we make use of the ideal Ohm's law for the magnetic vector potential and the electric scalar potential. By employing a suitable gauge condition that relates the potentials and working with a characteristic coordinate representation, we thus obtain an inhomogeneous first-order system of ordinary differential equations for the magnetic vector potential. Then, using the general solution of this system, we compute the magnetic field via the magnetic curl relation. Finally, we analyze the well-posedness of the corresponding Dirichlet-type initial-boundary value problem, specify compatibility conditions for the initial-boundary values, and outline the implementation of initial-boundary conditions., Comment: 13 pages, 1 figure, details added, minor corrections and improvements

Context. While the shapes of many observed bow shocks can be reproduced by simple astrosphere models, more elaborate approaches have recently been used to explain differing observable structures. Aims. By placing perturbations of an otherwise homogeneous interstellar medium in front of the astrospheric bow shock of the runaway blue supergiant $\lambda$ Cephei, the observable structure of the model astrosphere is significantly altered, providing insight into the origin of perturbed bow shock images. Methods. Three-dimensional single-fluid magnetohydrodynamic (MHD) models of stationary astrospheres were subjected to various types of perturbations and simulated until stationarity was reached again. As examples, simple perturbations of the available MHD parameters (number density, bulk velocity, temperature, and magnetic field) as well as a more complex perturbation were chosen. Synthetic observations were generated by line-of-sight integration of the model data, producing H$\alpha$, $70\,\mu$m dust emission, and bremsstrahlung maps of the perturbed astrosphere's evolution. Results. The resulting shock structures and observational images differ strongly depending on the type of the injected perturbation and the viewing angles, forming arc-like protrusions or bifurcations of the bow shock structure, as well as rings, arcs, and irregular structures detached from the bow shock., Comment: 13 pages, 12 figures. Accepted by and to be published in Astronomy & Astrophysics (A&A)

Measurements made with the Voyager 1 spacecraft indicate that significant levels of compressive fluctuations exist in the inner heliosheath. Some studies have already been performed with respect to the mirror-mode instability in the downstream region close to the solar wind termination shock, and here we extend the investigation to the whole inner heliosheath. We employ quasilinear theory and results from a global magnetohydrodynamic model of the heliosphere to compute the time evolution of both the temperature anisotropy and the energy density of the corresponding magnetic fluctuations, and we demonstrate their likely presence in the inner heliosheath. Furthermore, we compute the associated, locally generated density fluctuations. The results can serve as inputs for future models of the transport of compressible turbulence in the inner heliosheath., Comment: 7 pages, 5 figures

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., Comment: accepted at ApJ Letters

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-316 TeV for five different mass numbers. We present and analyze 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., Comment: 18 pages, 88 figures

The popular JF12 analytic model by Jansson & Farrar provides a quantitative description of the Galaxy's large-scale magnetic field, which is widely used in various astrophysical applications. However, both the poloidal X-type component and the spiral disk component of JF12 exhibit regions in which the magnetic divergence constraint is violated. We first propose a cure for this problem, resulting in a truly solenoidal large-scale spiral field. Second, the otherwise straight field lines of the X-type component exhibit kinks in the Galactic plane that, in addition to implying the presence of a singular current sheet, may pose difficulties for e.g., numerical tracing of cosmic-ray particles. We propose and discuss two possible strategies to mitigate this problem. Although all corrections are kept as minimal as possible, the extended set of model parameters will have to be carefully readjusted in order to fully restore the agreement to observational data that the unmodified JF12 field is based on. Furthermore, the performance of our improved version of the field model is quantitatively assessed by test simulations using the CRPropa Galactic cosmic-ray propagation code., Comment: 12 pages, 6 figures

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., Comment: 8 pages, 8 Figures

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., Comment: 11 pages, 4 Figures

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., Comment: ApJ, accepted

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.

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., Comment: 15 pages, 8 figures. Published in The Astrophysical Journal

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.

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.

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., Comment: 13 pages, 4 figures, accepted for publication at ApJ. The MPEG animation referenced in Figure 3 may be accessed via [http://www.tp4.rub.de/heliophysics/movies/pflow_drape+btime_annot.mp4]

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.

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., Comment: 16 pages, 10 figures. Comments welcome

Magnetic flux tubes reaching from the solar convective zone into the chromosphere have to pass through the relatively cool, and therefore non-ideal (i.e. resistive) photospheric region enclosed between the highly ideal sub-photospheric and chromospheric plasma. It is shown that stationary MHD equilibria of magnetic flux tubes which pass through this region require an inflow of photospheric material into the flux tube and a deviation from iso-rotation along the tube axis. This means that there is a difference in angular velocity of the plasma flow inside the tube below and above the non-ideal region. Both effects increase with decreasing cross section of the tube. Although for characteristic parameters of thick flux tubes the effect is negligible, a scaling law indicates its importance for small-scale structures. The relevance of this "inflow effect" for the expansion of flux tubes above the photosphere is discussed., Comment: 17 pages, 9 figures. Comments welcome

Coronal mass ejections (CMEs), which dynamically connect the solar surface to the far reaches of interplanetary space, represent a major anifestation of solar activity. They are not only of principal interest but also play a pivotal role in the context of space weather predictions. The steady improvement of both numerical methods and computational resources during recent years has allowed for the creation of increasingly realistic models of interplanetary CMEs (ICMEs), which can now be compared to high-quality observational data from various space-bound missions. This review discusses existing models of CMEs, characterizing them by scientific aim and scope, CME initiation method, and physical effects included, thereby stressing the importance of fully 3-D ('4pi') spatial coverage., Comment: 14 pages plus references. Comments welcome. Accepted for publication in Solar Physics (SUN-360 topical issue)

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., Comment: 8 pages, 6 figures, accepted in ApJ