In this work, we will be testing four different general \textit{f(R)}-gravity models, two of which are the more realistic models (namely the Starobinsky and the Hu-Sawicki models), to determine if they are viable alternative models to pursue a more vigorous constraining test upon them. For the testing of these models, we use 359 low- and intermediate-redshift Supernovae Type 1A data obtained from thRede SDSS-II/SNLS2 Joint Light-curve Analysis (JLA). We develop a Markov Chain Monte Carlo (MCMC) simulation to find a best-fitting function within reasonable ranges for each \textit{f(R)}-gravity model, as well as for the Lambda Cold Dark Matter ($\Lambda$CDM) model. For simplicity, we assume a flat universe with a negligible radiation density distribution. Therefore, the only difference between the accepted $\Lambda$CDM model and the \textit{f(R)}-gravity models will be the dark energy term and the arbitrary free parameters. By doing a statistical analysis and using the $\Lambda$CDM model as our "true model", we can obtain an indication whether or not a certain \textit{f(R)}-gravity model shows promise and requires a more in-depth view in future studies. In our results, we found that the Starobinsky model obtained a larger likelihood function value than the $\Lambda$CDM model, while still obtaining the cosmological parameters to be $\Omega_{m} = 0.268^{+0.027}_{-0.024}$ for the matter density distribution and $\bar{h} = 0.690^{+0.005}_{-0.005}$ for the Hubble uncertainty parameter. We also found a reduced Starobinsky model that are able to explain the data, as well as being statistically significant.

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.

The interpretation of recent observations of bow shocks around O-stars and the creation of corresponding models require a detailed understanding of the associated (magneto-)hydrodynamic structures. We base our study on three-dimensional numerical magneto-)hydrodynamical models, which are analyzed using the dynamically relevant parameters, in particular, the (magneto)sonic Mach numbers. The analytic Rankine-Hugoniot relation for HD and MHD are compared with those obtained by the numerical model. In that context we also show that the only distance which can be approximately determined is that of the termination shock, if it is a hydrodynamical shock. For MHD shocks the stagnation point does not, in general, lie on the inflow line, which is the line parallel to the inflow vector and passing through the star. Thus an estimate via the Bernoulli equation as in the HD case is, in general, not possible. We also show that in O-star astrospheres, distinct regions exist in which the fast, slow, Alfv\'enic, and sonic Mach numbers become lower than one, implying sub-slow magnetosonic as well as sub-fast and sub-sonic flows. Nevertheless, the analytic MHD Rankine Hugoniot relations can be used for further studies of turbulence and cosmic ray modulation.

Huge astrospheres or stellar wind bubbles influence the propagation of cosmic rays at energies up to the TeV range and can act as small-scale sinks decreasing the cosmic ray flux. We model such a sink (in 2D) by a sphere of radius 10\,pc embedded within a sphere of a radius of 1\,kpc. The cosmic ray flux is calculated by means of backward stochastic differential equations from an observer, which is located at $r_{0}$, to the outer boundary. It turns out that such small-scale sinks can influence the cosmic ray flux at the observer's location by a few permille (i.e\ a few 0.1\%), which is in the range of the observations by IceCube, Milagro and other large area telescopes.

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.

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.

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a time-dependent modulation model. By introducing recent theoretical advances in the transport coefficients in the model, computed intensities are compared with Voyager 1, International Monitoring Platform (IMP) 8, and Ulysses proton observations in search of compatibility. The effect of different modulation parameters on computed intensities is also illustrated. It is shown that this approach produces, on a global scale, realistic cosmic-ray proton intensities along the Voyager 1 spacecraft trajectory and at Earth upto ~2004, whereafter the computed intensities recovers much slower towards solar minimum than observed in the inner heliosphere. A modified time dependence in the diffusion coefficients is proposed to improve compatibility with the observations at Earth after ~2004. This modified time dependence led to an improved compatibility between computed intensities and the observations along the Voyager 1 trajectory and at Earth even after ~2004. An interesting result is that the cosmic-ray modulation during the present polarity cycle is not determined only by changes in the drift coefficient and tilt angle of the wavy current sheet, but is also largely dependent on changes in the diffusion coefficients.

A two-dimensional, time-dependent numerical model is used to calculate the modulation of cosmic rays in the heliosphere. Computations are compared to spacecraft observations in the inner and outer heliosphere. It is shown that the model produces cosmic ray proton intensities compatible to different spacecraft observations on a global scale, at Earth and along both Voyager spacecraft trajectories. The study reveals that when the same modulation parameters, which resulted in compatible intensities along Voyager 1, were assumed along the Voyager 2 trajectory, the model failed to reproduce the observations. The study also found that any change in diffusion parameters alone could not reproduce the cosmic ray observations along Voyager 2 so that changes to the heliospheric geometry were necessary i.e the computed intensities along both Voyager trajectories suggest that the heliosphere is asymmetric. Furthermore, $E >$ 70 MeV and 133-242 MeV proton intensities along Voyager 1 and 2 trajectories are predicted from end of 2012 onwards. It is shown that the computed intensities along Voyager 1 increase with an almost constant rate up to the heliopause. However, the model shows that Voyager 2 is still under the influence of temporal solar activity changes because of its relatively large distance to the heliopause. Along the Voyager 2 trajectory, the intensities remained generally constant for some time and should soon start to increase steadily.

A spatially independent model that calculates the time evolution of the electron spectrum in a spherically expanding pulsar wind nebula (PWN) is presented, allowing one to make broadband predictions for the PWN's non-thermal radiation. The source spectrum of electrons injected at the termination shock of the PWN is chosen to be a broken power law. In contrast to previous PWN models of a similar nature, the source spectrum has a discontinuity in intensity at the transition between the low and high-energy components. To test the model, it is applied to the young PWN G21.5--0.9, where it is found that a discontinuous source spectrum can model the emission at all wavelengths better than a continuous one. The model is also applied to the unidentified sources HESS J1427--608 and HESS J1507--622. Parameters are derived for these two candidate nebulae that are consistent with the values predicted for other PWNe. For HESS J1427--608 a present-day magnetic field of $B\s{age}=0.4\,\mu\text{G}$ is derived. As a result of the small present-day magnetic field, this source has a low synchrotron luminosity, while remaining bright at GeV/TeV energies. It is therefore possible to interpret HESS J1427--608 within the ancient PWN scenario. For the second candidate PWN HESS J1507--622, a present-day magnetic field of $B\s{age}=1.7\,\mu\text{G}$ is derived. Furthermore, for this candidate PWN a scenario is favoured in the present paper in which HESS J1507--622 has been compressed by the reverse shock of the supernova remnant.

Composite supernova remnants consist of a pulsar wind nebula located inside a shell-type remnant. The presence of a shell has implications on the evolution of the nebula, although the converse is generally not true. The purpose of this paper is two-fold. The first aim is to determine the effect of the pulsar's initial luminosity and spin-down rate, the supernova ejecta mass, and density of the interstellar medium on the evolution of a spherically-symmetric, composite supernova remnant expanding into a homogeneous medium. The second aim is to investigate the evolution of the magnetic field in the pulsar wind nebula when the the composite remnant expands into a non-uniform interstellar medium. The Euler conservation equations for inviscid flow, together with the magnetohydrodynamic induction law in the kinematic limit, are solved numerically for a number of scenarios where the ratio of magnetic to particle energy is $\sigma < 0.01$. The simulations in the first part of the paper is solved in a one-dimensional configuration. In the second part of the paper, the effect of an inhomogeneous medium on the evolution is studied using a two-dimensional, axis-symmetric configuration.

A spatially-independent model that can be used to calculate the temporal evolution of the electron/ positron spectrum in a spherically expanding pulsar wind nebula is presented. The model is applied to the young nebula G21.5-0.9, as well as to the unidentified TeV sources HESS J1427-608 and HESS J1507-622. The parameters derived from the model strengthens the idea that the unidentified sources can be identified as evolved pulsar wind nebulae.

The evolution of Pulsar Wind Nebulae (PWNe) plays a crucial role in interpreting the very high energy (VHE; > 10^11 eV) gamma-ray unidentified sources; and moreover it represents the only viable option to explain the discovery of several "dark sources" in the TeV gamma-ray (i.e. VHE gamma-ray sources without lower energies, radio or X-ray counterparts). The newest time-dependent modeling of PWNe presented in [1] and [2] has to be tested on a broader sample of young well-known PWNe and applied to the full-sample of "dark sources". The consequences of this interpretation go far beyond the interpretation of "dark sources": e.g. there could be strong implication in the origin of cosmic rays and (when considering a leptonic origin of the gamma-ray signal) they can be important for reinterpreting the detection of starburst galaxies in the TeV gamma-ray band.

In a Pulsar Wind Nebula (PWN), the lifetime of inverse Compton (IC) emitting electrons exceeds the lifetime of its progenitor pulsar (as well as its shell-type remnant), but it also exceeds the age of those that emit via synchrotron radiation. Therefore, during its evolution, the PWN can remain bright in IC so that its GeV-TeV gamma-ray flux remains high for timescales much larger (for 10^5 - 10^6 yrs) than the pulsar lifetime and the X-ray PWN lifetime. In this scenario, the magnetic field in the cavity induced by the wind of the progenitor star plays a crucial role. This scenario is in line with the discovery of several unidentified or "dark" sources in the TeV gamma-ray band without X-ray counterparts; and it is also finding confirmation in the recent discoveries at GeV gamma rays. Moreover, these consequences could be also important for reinterpreting the detection of starburst galaxies in the TeV gamma-ray band when considering a leptonic origin of the gamma-ray signal. Both theoretical aspects and their observational proofs will be discussed, as well as the first results of our new modeling code.

In this paper we explore the evolution of a PWN while the pulsar is spinning down. An MHD approach is used to simulate the evolution of a composite remnant. Particular attention is given to the adiabatic loss rate and evolution of the nebular field strength with time. By normalising a two component particle injection spectrum (which can reproduce the radio and X-ray components) at the pulsar wind termination shock to the time dependent spindown power, and keeping track with losses since pulsar/PWN/SNR birth, we show that the average field strength decreases with time as $t^{-1.3}$, so that the synchrotron flux decreases, whereas the IC gamma-ray flux increases, until most of the spindown power has been dumped into the PWN. Eventually adiabatic and IC losses will also terminate the TeV visibility and then eventually the GeV visibility.