In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

We assess the ionising effect of low energy protostellar cosmic rays in protoplanetary disks around a young solar mass star for a wide range of disk parameters. We assume a source of low energy cosmic rays located close to the young star which travel diffusively through the protoplanetary disk. We use observationally inferred values from nearby star-forming regions for the total disk mass and the radial density profile. We investigate the influence of varying the disk mass within the observed scatter for a solar mass star. We find that for a large range of disk masses and density profiles that protoplanetary disks are "optically thin" to low energy ($\sim$3 GeV) cosmic rays. At $R\sim10$au, for all of the disks that we consider ($M_\mathrm{disk}=6.0\times10^{-4} - 2.4\times 10^{-2}M_\odot$), the ionisation rate due to low energy stellar cosmic rays is larger than that expected from unmodulated galactic cosmic rays. This is in contrast to our previous results which assumed a much denser disk which may be appropriate for a more embedded source. At $R\sim70$au, the ionisation rate due to stellar cosmic rays dominates in $\sim$50% of the disks. These are the less massive disks with less steep density profiles. At this radius there is at least an order of magnitude difference in the ionisation rate between the least and most massive disk that we consider. Our results indicate, for a wide range of disk masses, that low energy stellar cosmic rays provide an important source of ionisation at the disk midplane at large radii ($\sim$70au).

Emission from the radioactive trace element $^{26}$Al has been observed throughout the Milky Way with the COMPTEL and INTEGRAL satellites. In particular the Doppler shifts measured with INTEGRAL connect $^{26}$Al with superbubbles, which may guide $^{26}$Al flows off spiral arms in the direction of Galactic rotation. In order to test this paradigm, we have performed galaxy-scale simulations of superbubbles with $^{26}$Al injection in a Milky Way-type galaxy. We produce all-sky synthetic $\gamma-$ray emission maps of the simulated galaxies. We find that the 1809keV emission from the radioactive decay of $^{26}$Al is highly variable with time and the observer's position. This allows us to estimate an additional systematic variability of 0.2dex for a star formation rate derived from $^{26}$Al for different times and measurement locations in Milky Way-type galaxies. High-latitude morphological features indicate nearby emission with correspondingly high integrated $\gamma-$ray intensities. We demonstrate that the $^{26}$Al scale height from our simulated galaxies depends on the assumed halo gas density. We present the first synthetic 1809keV longitude-velocity diagrams from 3D hydrodynamic simulations. The line-of-sight velocities for $^{26}$Al can be significantly different from the line-of-sight velocities associated with the cold gas. Over time, $^{26}$Al velocities consistent with the INTEGRAL observations, within uncertainties, appear at any given longitude, broadly supporting previous suggestions that $^{26}$Al injected into expanding superbubbles by massive stars may be responsible for the high velocities found in the INTEGRAL observations. We discuss the effect of systematically varying the location of the superbubbles relative to the spiral arms.

We investigate the ionising effect of low energy cosmic rays (CRs) from a young star on its protoplanetary disk (PPD). We consider specifically the effect of $\sim3\,$GeV protons injected at the inner edge of the PPD. An increase in the ionisation fraction as a result of these CRs could allow the magnetorotational instability to operate in otherwise magnetically dead regions of the disk. For the typical values assumed we find an ionisation rate of $\zeta_\mathrm{CR} \sim 10^{-17}\mathrm{s^{-1}}$ at $1\,$au. The transport equation is solved by treating the propagation of the CRs as diffusive. We find for increasing diffusion coefficients the CRs penetrate further in the PPD, while varying the mass density profile of the disk is found to have little effect. We investigate the effect of an energy spectrum of CRs. The influence of a disk wind is examined by including an advective term. For advective wind speeds between $1-100\mathrm{km\,s^{-1}}$ diffusion dominates at all radii considered here (out to 10$\,$au) for reasonable diffusion coefficients. Overall, we find that low energy CRs can significantly ionise the midplane of PPDs out to $\sim\,1\,$au. By increasing the luminosity or energy of the CRs, within plausible limits, their radial influence could increase to $\sim2\,$au at the midplane but it remains challenging to significantly ionise the midplane further out.

The redistribution of angular momentum is a long standing problem in our understanding of protoplanetary disk (PPD) evolution. The magnetorotational instability (MRI) is considered a likely mechanism. We present the results of a study involving multifluid global simulations including Ohmic dissipation, ambipolar diffusion and the Hall effect in a dynamic, self-consistent way. We focus on the turbulence resulting from the non-linear development of the MRI in radially stratified PPDs and compare with ideal MHD simulations. In the multifluid simulations the disk is initially set up to transition from a weak Hall dominated regime, where the Hall effect is the dominant non-ideal effect but approximately the same as or weaker than the inductive term, to a strong Hall dominated regime, where the Hall effect dominates the inductive term. As the simulations progress a substantial portion of the disk develops into a weak Hall dominated disk. We find a transition from turbulent to laminar flow in the inner regions of the disk, but without any corresponding overall density feature. We introduce a dimensionless parameter, $\alpha_\mathrm{RM}$, to characterise accretion with $\alpha_\mathrm{RM} \gtrsim 0.1$ corresponding to turbulent transport. We calculate the eddy turnover time, $t_\mathrm{eddy}$, and compared this with an effective recombination timescale, $t_\mathrm{rcb}$, to determine whether the presence of turbulence necessitates non-equilibrium ionisation calculations. We find that $t_\mathrm{rcb}$ is typically around three orders of magnitude smaller than $t_\mathrm{eddy}$. Also, the ionisation fraction does not vary appreciably. These two results suggest that these multifluid simulations should be comparable to single fluid non-ideal simulations.

Spectroscopic follow-up is a pre-requisite for studies of the formation and early evolution of brown dwarfs. Here we present IRTF/SpeX near-infrared spectroscopy of 30 candidate members of the young Upper Scorpius association, selected from our previous survey work. All 24 high confidence members are confirmed as young very low mass objects with spectral types from M5 to L1, 15-20 of them are likely brown dwarfs. This high yield confirms that brown dwarfs in Upper Scorpius can be identified from photometry and proper motions alone, with negligible contamination from field objects (<4%). Out of the 6 candidates with lower confidence, 5 might still be young very low mass members of Upper Scorpius, according to our spectroscopy. We demonstrate that some very low mass class II objects exhibit radically different near infrared (0.6 - 2.5micron) spectra from class III objects, with strong excess emission increasing towards longer wavelengths and partially filled in features at wavelengths shorter than 1.25micron. These characteristics can obscure the contribution of the photosphere within such spectra. Therefore, we caution that near infrared derived spectral types for objects with discs may be unreliable. Furthermore, we show that the same characteristics can be seen to some extent in all class II and even a significant fraction of class III objects (~40%), indicating that some of them are still surrounded by traces of dust and gas. Based on our spectra, we select a sample of objects with spectral types of M5 to L1, whose near-infrared emission represents the photosphere only. We recommend the use of these objects as spectroscopic templates for young brown dwarfs in the future.

We present the results of a multi-wavelength study of circumstellar discs around 44 young stellar objects in the 3 Myr old nearby Chamaeleon I star-forming region. In particular, we explore the far-infrared/submm regime using Herschel fluxes. We show that Herschel fluxes at 160-500$\,\mu$m can be used to derive robust estimates of the disc mass. The median disc mass is 0.005$M_{\odot}$ for a sample of 28 Class IIs and 0.006$M_{\odot}$ for 6 transition disks (TDs). The fraction of objects in Chamaeleon-I with at least the `minimum mass solar nebula' is 2-7%. This is consistent with previously published results for Taurus, IC348, $\rho$ Oph. Diagrams of spectral slopes show the effect of specific evolutionary processes in circumstellar discs. Class II objects show a wide scatter that can be explained by dust settling. We identify a continuous trend from Class II to TDs. Including Herschel fluxes in this type of analysis highlights the diversity of TDs. We find that TDs are not significantly different to Class II discs in terms of far-infrared luminosity, disc mass or degree of dust settling. This indicates that inner dust clearing occurs independently from other evolutionary processes in the discs.