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.

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.

The processes responsible for the effective longitudinal transport of solar energetic particles (SEPs) are still not completely understood. We address this issue by simulating SEP electron propagation using a spatially 2D transport model that includes perpendicular diffusion. By implementing, as far as possible, the most reasonable estimates of the transport (diffusion) coefficients, we compare our results, in a qualitative manner, to recent observations {at energies of 55 -- 105 keV}, focusing on the longitudinal distribution of the peak intensity, the maximum anisotropy and the onset time. By using transport coefficients which are derived from first principles, we limit the number of free parameters in the model to: (i) the probability of SEPs following diffusing magnetic field lines, quantified by $a \in [0,1]$, and (ii) the broadness of the Gaussian injection function. It is found that the model solutions are extremely sensitive to the magnitude of the {perpendicular} diffusion coefficient and relatively insensitive to the form of the injection function as long as a reasonable value of $a=0.2$ is used. We illustrate the effects of perpendicular diffusion on the model solutions and discuss the viability of this process as a dominant mechanism by which SEPs are transported in longitude. Lastly, we try to quantity the effectiveness of perpendicular diffusion as an interplay between the magnitude of the relevant diffusion coefficient and the SEP intensity gradient driving the diffusion process. It follows that perpendicular diffusion is extremely effective early in a SEP event when large intensity gradients are present, while the effectiveness quickly decreases with time thereafter.

We study the transport of high-energy particles in pulsar wind nebulae (PWN) using three-dimensional MHD (see Porth et al. (2014) for details) and test-particle simulations, as well as a Fokker-Planck particle transport model. The latter includes radiative and adiabatic losses, diffusion, and advection on the background flow of the simulated MHD nebula. By combining the models, the spatial evolution of flux and photon index of the X-ray synchrotron emission is modelled for the three nebulae G21.5-0.9, the inner regions of Vela, and 3C 58, thereby allowing us to derive governing parameters: the magnetic field strength, average flow velocity and spatial diffusion coefficient. For comparison, the nebulae are also modelled with the semi-analytic Kennel & Coroniti (1984) model but the Porth et al. (2014) model generally yields better fits to the observational data. We find that high velocity fluctuations in the turbulent nebula (downstream of the termination shock) give rise to efficient diffusive transport of particles, with average P\'eclet number close to unity, indicating that both advection and diffusion play an important role in particle transport. We find that the diffusive transport coefficient of the order of $\sim2\times 10^{27} (L_{\rm s}/0.42\rm Ly) cm^{2}s^{-1}$ ($L_{\rm s}$ is the size of the termination shock) is independent of energy up to extreme particle Lorentz factors of $\gamma_{p}\sim10^{10}$.