Jovian electrons serve as an important test-particle distribution in the inner heliosphere and have been used extensively in the past to study the (diffusive) transport of cosmic rays in the inner heliosphere. With new limits on the Jovian source function (i.e. the particle intensity just outside the Jovian magnetosphere), and a new set of in-situ observations at 1 AU for both cases of good and poor magnetic connection between the source and observer, we revisit some of these earlier simulations. We aim to find the optimal numerical set-up that can be used to simulate the propagation of 6 MeV Jovian electrons in the inner heliosphere. Using such a set-up, we further aim to study the residence (propagation) times of these particles for different levels of magnetic connection between Jupiter and an observer at Earth (1 AU). Using an advanced Jovian electron propagation model based on the stochastic differential equation (SDE) approach, we calculate the Jovian electron intensity for different model parameters. A comparison with observations leads to an optimal numerical set-up, which is then used to calculate the so-called residence (propagation) times of these particles. Comparing to in-situ observations, we are able to derive transport parameters that are appropriate to study the propagation of 6 MeV Jovian electrons in the inner heliosphere. Moreover, using these values, we show that the method of calculating the residence time applied in former literature is not suited to being interpreted as the propagation time of physical particles. This is due to an incorrect weighting of the probability distribution. We propose and apply a new method, where the results from each pseudo-particle are weighted by its resulting phase-space density (i.e. the number of physical particles that it represents). Thereby we obtain more reliable estimates for the propagation time.

Terzan 5 is the only Galactic globular cluster that has plausibly been detected at very-high energies by the High Energy Stereoscopic System. It has an unexpectedly asymmetric very-high-energy morphology that is offset from the cluster center, in addition to a large-scale, offset radio structure and compact diffuse X-ray emission associated with this cluster. We present new data from the Fermi Large Area Telescope on this source. We model the updated broadband spectral energy distribution, attributing this to cumulative pulsed emission from a population of embedded millisecond pulsars as well as unpulsed emission from the interaction of their leptonic winds with the ambient magnetic and soft-photon fields. In particular, our model invokes unpulsed synchrotron and inverse Compton components to model the radio and TeV data, cumulative pulsed curvature radiation to fit the Fermi data, and explains the hard Chandra X-ray spectrum via a "new" cumulative synchrotron component from electron-positron pairs within the pulsar magnetospheres that has not been implemented before. We find reasonable spectral fits for plausible model parameters. We also derive constraints on the millisecond pulsar luminosity function using the diffuse X-ray data and the Chandra sensitivity. Future higher-quality spectral and spatial data will help discriminate between competing scenarios (such as dark matter annihilation, white dwarf winds or hadronic interactions) proposed for the broadband emission as well as constrain degenerate model parameters.

Insights into the processes of Solar Energetic Particle (SEP) propagation are essential for understanding how solar eruptions affect the radiation environment of near-Earth space. SEP propagation is influenced by turbulent magnetic fields in the solar wind, resulting in stochastic transport of SEPs to Earth. Multi-spacecraft observations suggest that the cross-field propagation shapes the SEP fluxes at Earth strongly. However, modelling SEP cross-field transport as spatial diffusion has been shown to be insufficient without use of unrealistically large cross-field diffusion coefficients. Recent work has shown that the early-time propagation of energetic particles across the mean field direction in turbulent fields is not diffusive, as the particles propagating along meandering field lines. This early-time transport mode results in fast access of the particles across the mean field direction, in agreement with the SEP observations. In this work, we demonstrate the significance of turbulence strength on evolution of the SEP radiation environment near Earth. We calculate the transport parameters with a turbulence transport model, parametrised by the SEP parallel scattering mean free path at 1~AU, $\lambda_{\parallel}^{*}$, and show that the parallel and cross-field transport are connected, with conditions resulting in slow parallel transport corresponding to wider events. We find a scaling $\sigma_{\phi,\, \mathrm{max}}\propto (1/\lambda_{\parallel}^{*})^{1/4}$ for the Gaussian fitting of the longitudinal distribution of maximum intensities. The longitudes with highest intensities are shifted towards the west for strong scattering conditions. Our results emphasise the importance of understanding both the SEP transport and the interplanetary turbulence conditions for modelling and predicting the SEP radiation environment at Earth.

Context. Current solar energetic particle (SEP) propagation models describe the effects of interplanetary plasma turbulence on SEPs as diffusion, using a Fokker-Planck (FP) equation. However, FP models cannot explain the observed fast access of SEPs across the average magnetic field to regions that are widely separated in longitude within the heliosphere without using unrealistically strong cross-field diffusion. Aims. We study whether the recently suggested early non-diffusive phase of SEP propagation can explain the wide SEP events with realistic particle transport parameters. Methods. We used a novel model that accounts for the SEP propagation along field lines that meander as a result of plasma turbulence. Such a non-diffusive propagation mode has been shown to dominate the SEP cross-field propagation early in the SEP event history. We compare the new model to the traditional approach, and to SEP observations. Results. Using the new model, we reproduce the observed longitudinal extent of SEP peak fluxes that are characterised by a Gaussian profile with $\sigma=30-50^\circ$, while current diffusion theory can only explain extents of 11$^\circ$ with realistic diffusion coefficients. Our model also reproduces the timing of SEP arrival at distant longitudes, which cannot be explained using the diffusion model. Conclusions. The early onset of SEPs over a wide range of longitudes can be understood as a result of the effects of magnetic field-line random walk in the interplanetary medium and requires an SEP transport model that properly describes the non-diffusive early phase of SEP cross-field propagation.

Pair cascades from millisecond pulsars (MSPs) may be a primary source of Galactic electrons and positrons that contribute to the increase in positron flux above 10 GeV as observed by PAMELA and AMS-02. The Fermi Large Area Telescope (LAT) has increased the number of detected gamma-ray MSPs tremendously. Light curve modelling furthermore favours abundant pair production in MSP magnetospheres, so that models of primary cosmic-ray positrons from pulsars should include the contribution from the larger numbers of MSPs and their potentially higher positron output per source. We model the contribution of Galactic MSPs to the terrestrial cosmic-ray electron / positron flux by using a population synthesis code to predict the source properties of present-day MSPs. We simulate pair spectra assuming an offset-dipole magnetic field which boosts pair creation rates. We also consider positrons and electrons that have additionally been accelerated to very high energies in the strong intrabinary shocks in black widow (BW) and redback (RB) binary systems. We transport these particles to Earth by calculating their diffusion and the radiative energy losses they suffer in the Galaxy using a model. Our model particle flux increases for non-zero offsets of the magnetic polar caps. We find that pair cascades from MSP magnetospheres contribute only modestly around a few tens of GeV to the measured fluxes. BW and RB fluxes may reach a few tens of percent of the observed flux up to a few TeV. Future observations should constrain the source properties in this case.

Observations by the Fermi Large Area Telescope of gamma-ray millisecond pulsar light curves imply copious pair production in their magnetospheres, and not exclusively in those of younger pulsars. Such pair cascades may be a primary source of Galactic electrons and positrons, contributing to the observed enhancement in positron flux above ~10 GeV. Fermi has also uncovered many new millisecond pulsars, impacting Galactic stellar population models. We investigate the contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray electrons and positrons. Our population synthesis code predicts the source properties of present-day millisecond pulsars. We simulate their pair spectra invoking an offset-dipole magnetic field. We also consider positrons and electrons that have been further accelerated to energies of several TeV by strong intrabinary shocks in black widow and redback systems. Since millisecond pulsars are not surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs freely escape and undergo losses only in the intergalactic medium. We compute the transported pair spectra at Earth, following their diffusion and energy loss through the Galaxy. The predicted particle flux increases for non-zero offsets of the magnetic polar caps. Pair cascades from the magnetospheres of millisecond pulsars are only modest contributors around a few tens of GeV to the lepton fluxes measured by AMS-02, PAMELA, and Fermi, after which this component cuts off. The contribution by black widows and redbacks may, however, reach levels of a few tens of percent at tens of TeV, depending on model parameters.

Millisecond pulsars occur abundantly in globular clusters. They are expected to be responsible for several spectral components in the radio through gamma-ray waveband (e.g., involving synchrotron and inverse Compton emission), as have been seen by Radio Telescope Effelsberg, Chandra X-ray Observatory, Fermi Large Area Telescope, and the High Energy Stereoscopic System (H.E.S.S.) in the case of Terzan 5 (with fewer spectral components seen for other globular clusters). H.E.S.S. has recently performed a stacking analysis involving 15 non-detected globular clusters and obtained quite constraining average flux upper limits above 230 GeV. We present a model that assumes millisecond pulsars as sources of relativistic particles and predicts multi-wavelength emission from globular clusters. We apply this model to the population of clusters mentioned above to predict the average spectrum and compare this to the H.E.S.S. upper limits. Such comparison allows us to test whether the model is viable, leading to possible constraints on various average cluster parameters within this framework.

Pulsars are believed to be sources of relativistic electrons and positrons. The abundance of detections of gamma-ray millisecond pulsars by Fermi Large Area Telescope coupled with their light curve characteristics that imply copious pair production in their magnetospheres, motivated us to investigate this old pulsar population as a source of Galactic electrons and positrons and their contribution to the enhancement in cosmic-ray positron flux at GeV energies. We use a population synthesis code to predict the source properties (number, position, and power) of the present-day Galactic millisecond pulsars, taking into account the latest Fermi and radio observations to calibrate the model output. Next, we simulate pair cascade spectra from these pulsars using a model that invokes an offset-dipole magnetic field. We assume free escape of the pairs from the pulsar environment. We then compute the cumulative spectrum of transported electrons and positrons at Earth, following their diffusion and energy losses as they propagate through the Galaxy. Our results indicate that the predicted particle flux increases for non-zero offsets of the magnetic polar caps. Comparing our predicted spectrum and positron fraction to measurements by AMS-02, PAMELA, and Fermi, we find that millisecond pulsars are only modest contributors at a few tens of GeV, after which this leptonic spectral component cuts off. The positron fraction is therefore only slightly enhanced above 10 GeV relative to a background flux model. This implies that alternative sources such as young, nearby pulsars and supernova remnants should contribute additional primary positrons within the astrophysical scenario.

Recent observations by the Fermi Large Area Telescope (LAT) and the High Energy Stereoscopic System (H.E.S.S.) have revealed globular clusters (GC) to be sources of high-energy (HE) and very-high-energy (VHE) gamma-rays. It has been suggested that the presence of large numbers of millisecond pulsars (MSPs) within these clusters may be either directly responsible for these gamma-ray fluxes through emission of pulsed curvature radiation, or indirectly through the injection of relativistic leptons into the cluster. These relativistic particles are plausibly re-accelerated in shocks, created by the collision of stellar winds, before interacting with the soft-photon radiation field set up by the stellar population of the host cluster. Inverse Compton (IC) scattering then produces gamma-radiation in the TeV band. In order to calculate the IC spectrum, an accurate profile for the energy density of the soft-photon field is required. We construct such a profile by deriving a radially-dependent expression for the stellar energy density, and then solve it numerically. As a next step, the average energy density values for three different regions of the cluster (demarcated by its core, half-mass, and tidal radii) are determined, which we consequently import into an existing radiation code to predict the TeV gamma-ray spectrum. As an application, we consider the case of Terzan 5, boasting a population of 34 radio MSPs, and compare our predicted spectrum with that measured by H.E.S.S.

Globular clusters (GCs) are astronomical tapestries embroidered with an abundance of exotic stellar-type objects. Their high age promises a rich harvest of evolved stellar products, while the deep potential wells and high mass densities at their centres probably facilitate the formation of multiple-member stellar systems via increased stellar encounter rates. The ubiquity of low-mass X-ray binaries, thought to be the progenitors of millisecond pulsars (MSPs), explain the large number of observed GC radio pulsars and X-ray counterparts. The Fermi Large Area Telescope (LAT) recently unveiled the first gamma-ray GC pulsar (PSR J1823-3021A). The first observations of GCs in the GeV and TeV bands furthermore created much excitement, and in view of the above, it seems natural to explain these high-energy lanterns by investigating an MSP origin. An MSP population is expected to radiate several pulsed spectral components in the radio through gamma-ray wavebands, in addition to being sources of relativistic particles. The latter may interact with background photons in the clusters producing TeV excesses, while they may also radiate synchrotron photons as they traverse the cluster magnetic field. We present multiwavelength modelling results for Terzan 5. We also briefly discuss some alternative interpretations for the observed GC gamma-ray signals.

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.

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.

Diffuse X-ray emission has recently been detected from the globular cluster (GC) Terzan 5, extending out to ~2.5' from the cluster centre. This emission may arise from synchrotron radiation (SR) by energetic leptons being injected into the cluster by the resident millisecond pulsar (MSP) population that interact with the cluster field. These leptons may also be reaccelerated in shocks created by collisions of pulsar winds, and may interact with bright starlight and cosmic microwave background (CMB) photons, yielding gamma rays at very high energies (VHE) through the inverse Compton (IC) process. In the GeV range, Fermi Large Area Telescope (LAT) has detected a population of GCs, very plausibly including Terzan 5, their spectral properties and energetics being consistent with cumulative magnetospheric emission from a population of MSPs. H.E.S.S. has furthermore detected a VHE excess in the direction of Terzan 5. One may derive constraints on the number of MSPs, N_tot, and the radial profiles of the GC B-field, stellar energy density, as well as the diffusion coefficient using the spatially-resolved X-ray, high-energy (HE), and VHE fluxes. If the Fermi LAT flux is due to magnetospheric processes, it will scale with the number of visible gamma-ray MSPs, N_vis. The HE spectrum therefore provides an independent way of constraining the number of MSPs (since N_tot >= N_vis). Consequently, the synthesis of available multiwavelength data presents a unique opportunity to constrain several parameters of the GC Terzan 5.

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.

The Fermi Large Area Telescope (LAT) has recently detected a population of globular clusters (GCs) in high-energy (HE) gamma-rays. Their spectral properties and energetics are consistent with cumulative emission from a population of millisecond pulsars (MSPs) hosted by these clusters. For example, the HE spectra exhibit fairly hard power-law indices and cutoffs around a few GeV, typical of pulsed spectra measured for the gamma-ray pulsar population. The energetics may be used to constrain the number of visible MSPs in the cluster (N_vis), assuming canonical values for the average gamma-ray efficiency and spin-down power. This interpretation is indeed strengthened by the fact that the first gamma-ray MSP has now been identified in the GC NGC 6624, and this MSP is responsible for almost all of the HE emission from this cluster. On the other hand, it has been argued that the MSPs are also sources of relativistic leptons which may be reaccelerated in shocks originating in collisions of stellar winds in the cluster core, and may upscatter bright starlight and cosmic microwave background photons to very high energies. Therefore, this unpulsed component may give an independent constraint on the total number of MSPs (N_tot) hosted in the GC, for a given cluster magnetic field B and diffusion coefficient k_0. Lastly, the transport properties of the energetic leptons may be further constrained using multiwavelength data, e.g., to infer the radial dependence of k_0 and B. We present results on our modeling of the pulsed and unpulsed gamma-ray fluxes from the GC Terzan 5.

In this work linear stability analysis of a four-fluid optically thin plasma consisting of electrons, ions, neutral atoms, and charged dust particles is performed with respect to the radiation-condensation (RC) instability. The energy budget of the plasma involves the input from heating through photo-electron emission by dust particles under external ultraviolet radiation as well as radiative losses in inelastic electron-neutral, electron-ion, neutral-neutral collisions. It is shown that negatively charged particles stimulate the RC instability in the sense that the conditions for the instability to hold are wider than similar conditions in a single-fluid description.

We have studied the impact of cosmic-ray acceleration in SNR on the spectra of cosmic-ray nuclei in the Galaxy using a series expansion of the propagation equation, which allows us to use analytical solutions for part of the problem and an efficient numerical treatment of the remaining equations and thus accurately describes the cosmic-ray propagation on small scales around their sources in three spatial dimensions and time. We found strong variations of the cosmic-ray nuclei flux by typically 20% with occasional spikes of much higher amplitude, but only minor changes in the spectral distribution. The locally measured spectra of primary cosmic rays fit well into the obtained range of possible spectra. We further showed that the spectra of the secondary element Boron show almost no variations, so that the above findings also imply significant fluctuations of the Boron-to-Carbon ratio. Therefore the commonly used method of determining CR propagation parameters by fitting secondary-to-primary ratios appears flawed on account of the variations that these ratios would show throughout the Galaxy.

The remanent magnetisation of meteorite material in the solar system indicates that magnetic fields of several Gauss are present in the protoplanetary disk. It is shown that such relativelystrong magnetic fields can be generated in dusty protoplanetary disks by relative shear motions of the charged dust and the neutral gas components. Self-consistent multi-fluid simulations show that for typical plasma parameters shear flows with collisional momentum transfer between the different components of the dusty plasma result in a very efficient generation of magnetic fields with strengths of about 0.1-1 Gauss on spatial scales of two astronomical units in about one year. Based on the simulations one may predict that future measurements will reveal the strong magnetisation of circumstellar disks around young stellar objects.