Outflows from starburst galaxies can be driven by thermal pressure, radiation and cosmic rays. We present an analytic phenomenological model that accounts for these contributions simultaneously to investigate their effects on the hydrodynamical properties of outflows. We assess the impact of energy injection, wind opacity, magnetic field strength and the mass of the host galaxy on flow velocity, temperature, density and pressure profiles. For an M82-like wind, a thermally-dominated driving mechanism is found to deliver the fastest and hottest wind. Radiation-driven winds in typical starburst-galaxy configurations are unable to attain the higher flow velocities and temperatures associated with thermal and cosmic ray-driven systems, leading to higher wind densities which would be more susceptible to cooling and fragmentation at lower altitudes. High opacity winds are more sensitive to radiative driving, but terminal flow velocities are still lower than those achieved by other driving mechanisms at realistic opacities. We demonstrate that variations in the outflow magnetic field can influence its coupling with cosmic rays, where stronger fields enable greater streaming but less driving near the base of the flow, instead with cosmic rays redirecting their driving impact to higher altitudes. The gravitational potential is less important in M82-like wind configurations, and substantial variations in the flow profiles only emerge at high altitude in massive haloes. This model offers a more generalised approach to examine the large scale hydrodynamical properties for a wide variety of starburst galaxies.

Large-scale outflows from starburst galaxies are multi-phase, multi-component fluids. Charge-exchange lines which originate from the interfacing surface between the neutral and ionised components are a useful diagnostic of the cold dense structures in the galactic outflow. From the charge-exchange lines observed in the nearby starburst galaxy M82, we conduct surface-to-volume analyses and deduce that the cold dense clumps in its galactic outflow have flattened shapes, resembling a hamburger or a pancake morphology rather than elongated shapes. The observed filamentary H$\alpha$ features are therefore not prime charge-exchange line emitters. They are stripped material torn from the slow moving dense clumps by the faster moving ionised fluid which are subsequently warmed and stretched into elongated shapes. Our findings are consistent with numerical simulations which have shown that cold dense clumps in galactic outflows can be compressed by ram pressure, and also progressively ablated and stripped before complete disintegration. We have shown that some clumps could survive their passage along a galactic outflow. These are advected into the circumgalactic environment, where their remnants would seed condensation of the circumgalactic medium to form new clumps. The infall of these new clumps back into the galaxy and their subsequent re-entrainment into the galactic outflow form a loop process of galactic material recycling.

Quenching of star-formation has been identified in many starburst and post-starburst galaxies, indicating burst-like star-formation histories (SFH) in the primordial Universe. We have investigated the role of high energy cosmic rays (CRs) in such environments, particularly how they could contribute to this burst-like SFH via quenching and feedback. These high energy particles interact with the baryon and radiation fields of their host via hadronic processes to produce secondary leptons. The secondary particles then also interact with ambient radiation fields to generate X-rays through inverse-Compton scattering. In addition, they can thermalise directly with the semi-ionised medium via Coulomb processes. Heating at a rate of $\sim 10^{-25} \; \text{erg}~\text{cm}^{-3}~\text{s}^{-1}$ can be attained by Coulomb processes in a star-forming galaxy with one core-collapse SN event per decade, and this is sufficient to cause quenching of star-formation. At high-redshift, a substantial amount of CR secondary electron energy can be diverted into inverse-Compton X-ray emission. This yields an X-ray luminosity of above $10^{41}~\text{erg}~\text{s}^{-1}$ by redshift $z=7$ which drives a further heating effect, operating over larger scales. This would be able to halt inflowing cold gas filaments, strangulating subsequent star-formation. We selected a sample of 16 starburst and post-starburst galaxies at $7\lesssim z \lesssim 9$ and determine the star-formation rates they could have sustained. We applied a model with CR injection, propagation and heating to calculate energy deposition rates in these 16 sources. Our calculations show that CR feedback cannot be neglected as it has the strength to suppress star-formation in these systems. We also show that their currently observed quiescence is consistent with the suffocation of cold inflows, probably by a combination of X-ray and CR heating.

We investigate the interactions of energetic hadronic particles with the media in outflows from star-forming protogalaxies. These particles undergo pion-producing interactions which can drive a heating effect in the outflow, while those advected by the outflow also transport energy beyond the galaxy, heating the circumgalactic medium. We investigate how this process evolves over the length of the outflow and calculate the corresponding heating rates in advection-dominated and diffusion-dominated cosmic ray transport regimes. In a purely diffusive transport scenario, we find the peak heating rate reaches $10^{-26}\;\! {\rm erg~cm}^{-3}\;\! {\rm s}^{-1}$ at the base of the outflow where the wind is driven by core-collapse supernovae at an event rate of 0.1 $\text{yr}^{-1}$, but does not extend beyond 2 kpc. In the advection limit, the peak heating rate is reduced to $10^{-28}\;\! {\rm erg~cm}^{-3}\;\! {\rm s}^{-1}$, but its extent can reach to tens of kpc. Around 10% of the cosmic rays injected into the system can escape by advection with the outflow wind, while the remaining cosmic rays deliver an important interstellar heating effect. We apply our cosmic ray heating model to the recent observation of the high-redshift galaxy MACS1149-JD1 and show that it could account for the quenching of a previous starburst inferred from spectroscopic observations. Re-ignition of later star-formation may be caused by the presence of filamentary circumgalactic inflows which are reinstated after cosmic ray heating has subsided.

We investigate the interactions of energetic hadronic particles (cosmic ray protons) with photons and baryons in protogalactic environments, where the target photons are supplied by the first generations of stars to form in the galaxy and the cosmological microwave background, while the target baryons are the interstellar and circumgalactic medium. We show that pair-production and photo-pion processes are the dominant interactions at particle energies above $10^{19}\;\! {\rm eV}$, while ${\rm pp}$-interaction pion-production dominates at the lower energies in line with expectations from, for example $\gamma$-ray observations of star-forming galaxies and dense regions of our own galaxy's interstellar medium. We calculate the path lengths for the interaction channels and determine the corresponding rates of energy deposition. We have found that protogalactic magnetic fields and their evolution can significantly affect the energy transport and energy deposition processes of cosmic rays. Within a Myr after the onset of star-formation the magnetic field in a protogalaxy could attain a strength sufficient to confine all but the highest energy particles within the galaxy. This enhances the cosmic ray driven self-heating of the protogalaxy to a rate of around $10^{-24}\;\! {\rm erg} \;\!{\rm cm}^{-3}\;\! {\rm s}^{-1}$ for a galaxy with strong star-forming activity that yields 1 core collapse SN event per year. This heating power exceeds even that due to radiative emission from the protogalaxy's stellar populations. However, in a short window before the protogalaxy is fully magnetised, energetic particles could stream across the galaxy freely, delivering energy into the circumgalactic and intergalactic medium.

Gammapy is a Python package for high-level gamma-ray data analysis built on Numpy, Scipy and Astropy. It enables us to analyze gamma-ray data and to create sky images, spectra and lightcurves, from event lists and instrument response information, and to determine the position, morphology and spectra of gamma-ray sources. So far Gammapy has mostly been used to analyze data from H.E.S.S. and Fermi-LAT, and is now being used for the simulation and analysis of observations from the Cherenkov Telescope Array (CTA). We have proposed Gammapy as a prototype for the CTA science tools. This contribution gives an overview of the Gammapy package and project and shows an analysis application example with simulated CTA data.

Galaxies at high redshifts with strong star formation are sources of high-energy cosmic rays. These cosmic rays interact with the baryon and radiation fields of the galactic environment via photo-pair, photo-pion and proton-proton processes to produce charged and neutral pions, neutrons and protons. The cosmic rays thereby deposit energy into the interstellar medium (ISM) as they propagate. We show how energy transport and deposition by ultra high-energy cosmic rays is regulated by the evolution of the galaxy, in particular by the development of the galactic magnetic field. We show how the particle-driven energy deposition can influence the thermal evolution of the host and its surroundings. Using a parametric protogalaxy model, we calculate the heating effect on the ISM as the cosmic rays are increasingly confined by the magnetic evolution of the galaxy.

The field of gamma-ray astronomy has seen important progress during the last decade, yet there exists so far no common software framework for the scientific analysis of gamma-ray telescope data. We propose to fill this gap by means of the GammaLib software, a generic library that we have developed to support the analysis of gamma-ray event data. GammaLib has been written in C++ and all functionality is available in Python through an extension module. On top of this framework we have developed the ctools software package, a suite of software tools that enables building of flexible workflows for the analysis of Imaging Air Cherenkov Telescope event data. The ctools are inspired by science analysis software available for existing high-energy astronomy instruments, and they follow the modular ftools model developed by the High Energy Astrophysics Science Archive Research Center. The ctools have been written in Python and C++, and can be either used from the command line, via shell scripts, or directly from Python. In this paper we present the GammaLib and ctools software versions 1.0 that have been released end of 2015. GammaLib and ctools are ready for the science analysis of Imaging Air Cherenkov Telescope event data, and also support the analysis of Fermi-LAT data and the exploitation of the COMPTEL legacy data archive. We propose to use ctools as the Science Tools software for the Cherenkov Telescope Array Observatory.

In the past decade imaging atmospheric Cherenkov telescope arrays such as H.E.S.S., MAGIC, VERITAS, as well as the Fermi-LAT space telescope have provided us with detailed images and spectra of the gamma-ray universe for the first time. Currently the gamma-ray community is preparing to build the next-generation Cherenkov Telecope Array (CTA), which will be operated as an open observatory. Gammapy (available at https://github.com/gammapy/gammapy under the open-source BSD license) is a new in-development Astropy affiliated package for high-level analysis and simulation of astronomical gamma-ray data. It is built on the scientific Python stack (Numpy, Scipy, matplotlib and scikit-image) and makes use of other open-source astronomy packages such as Astropy, Sherpa and Naima to provide a flexible set of tools for gamma-ray astronomers. We present an overview of the current Gammapy features and example analyses on real as well as simulated gamma-ray datasets. We would like Gammapy to become a community-developed project and a place of collaboration between scientists interested in gamma-ray astronomy with Python. Contributions welcome!

One of the most prominent features of the $\gamma$-ray sky is the emission from our own Galaxy. The Galactic plane has been observed by Fermi-LAT in GeV and H.E.S.S. in TeV light. Fermi has modeled the Galactic emission as the sum of a complex 'diffuse' emission model with the predominately point source catalogs of 1FHL and 2FGL, while H.E.S.S. has primarily detected extended TeV sources. At GeV energies, Galactic diffuse emission dominates the $\gamma$-ray Milky Way but, as sources have hard spectra, it is likely their emission dominates at TeV energies. Generally the spatial shape and fraction of source emission compared to diffuse emission in the Galactic plane is not well known and is dependent on the source detection method, threshold and diffuse emission modeling methods used. We present a simple image-analysis based method applied to Fermi-LAT data from 10 GeV to 500 GeV, covering a region of +/- 5 degrees in Galactic latitude and +/- 100 degrees in Galactic longitude, to separate source and diffuse emission. This method involves elongated filter smoothing, combined with significance clipping to exclude sources. We test the method against models based on the 1FHL catalog and very simple model Galaxies to evaluate the response for an input of known fraction and shape of diffuse and source emission.