@article{ThomasPfrommerPakmor2021, author = {Thomas, Timon and Pfrommer, Christoph and Pakmor, R{\"u}diger}, title = {A finite volume method for two-moment cosmic ray hydrodynamics on a moving mesh}, series = {Monthly notices of the Royal Astronomical Society}, volume = {503}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab397}, pages = {2242 -- 2264}, year = {2021}, abstract = {We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magnetohydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.}, language = {en} } @article{WerhahnPfrommerGirichidis2021, author = {Werhahn, Maria and Pfrommer, Christoph and Girichidis, Philipp}, title = {Cosmic rays and non-thermal emission in simulated galaxies - III. Probing cosmic-ray calorimetry with radio spectra and the FIR-radio correlation}, series = {Monthly notices of the Royal Astronomical Society}, volume = {508}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab2535}, pages = {4072 -- 4095}, year = {2021}, abstract = {An extinction-free estimator of the star formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR), and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic-ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magnetohydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star light and cosmic microwave background (CMB) photons while less energy is converted into synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission flattens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts.}, language = {en} } @article{WerhahnPfrommerGirichidisetal.2021, author = {Werhahn, Maria and Pfrommer, Christoph and Girichidis, Philipp and Puchwein, Ewald and Pakmor, R{\"u}diger}, title = {Cosmic rays and non-thermal emission in simulated galaxies}, series = {Monthly notices of the Royal Astronomical Society}, volume = {505}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford University Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stab1324}, pages = {3273 -- 3294}, year = {2021}, abstract = {Current-day cosmic ray (CR) propagation studies use static Milky Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magnetohydrodynamic (MHD) simulations of isolated galaxies with the moving-mesh code arepo that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into account all relevant loss processes. We show that this steady-state assumption is well justified in the disc and generally for regions that emit non-thermal radio and gamma rays. Additionally, we model the spectra of primary electrons, accelerated by supernova remnants, and secondary electrons and positrons produced in hadronic CR proton interactions with the gas. We find that proton spectra above 10 GeV only weakly depend on galactic radius, while they acquire a radial dependence at lower energies due to Coulomb interactions. Radiative losses steepen the spectra of primary CR electrons in the central galactic regions, while diffusive losses dominate in the outskirts. Secondary electrons exhibit a steeper spectrum than primaries because they originate from the transported steeper CR proton spectra. Consistent with Voyager-1 and AMS-02 data, our models (i) show a turnover of proton spectra below GeV energies due to Coulomb interactions so that electrons start to dominate the total particle spectra and (ii) match the shape of the positron fraction up to 10 GeV. We conclude that our steady-state CR modelling in MHD CR galaxy simulations is sufficiently realistic to capture the dominant transport effects shaping their spectra, arguing for a full MHD treatment to accurately model CR transport in the future.}, language = {en} } @article{WilhelmTelezhinskyDwarkadasetal.2020, author = {Wilhelm, Alina and Telezhinsky, Igor and Dwarkadas, Vikram V. and Pohl, Martin}, title = {Stochastic re-acceleration and magnetic-field damping in Tycho's supernova remnant}, series = {Astronomy and astrophysics}, volume = {639}, journal = {Astronomy and astrophysics}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201936079}, pages = {14}, year = {2020}, abstract = {Context. Tycho's supernova remnant (SNR) is associated with the historical supernova (SN) event SN 1572 of Type Ia. The explosion occurred in a relatively clean environment, and was visually observed, providing an age estimate. This SNR therefore represents an ideal astrophysical test-bed for the study of cosmic-ray acceleration and related phenomena. A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvenic drift are needed to explain the rather soft radio spectrum and the narrow rims observed in X-rays. Aims. We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary process. The re-acceleration of particles in the turbulent region immediately downstream of the shock should be efficient enough to impact particle spectra over several decades in energy. The so-called Alfvenic drift and particle feedback on the shock structure are not required in this scenario. Additionally, we investigate whether synchrotron losses or magnetic-field damping play a more profound role in the formation of the non-thermal filaments. Methods. We solved the full particle transport equation in test-particle mode using hydrodynamic simulations of the SNR plasma flow. The background magnetic field was either computed from the induction equation or follows analytic profiles, depending on the model considered. Fast-mode waves in the downstream region provide the diffusion of particles in momentum space. Results. We show that the broadband spectrum of Tycho can be well explained if magnetic-field damping and stochastic re-acceleration of particles are taken into account. Although not as efficient as standard diffusive shock acceleration, stochastic acceleration leaves its imprint on the particle spectra, which is especially notable in the emission at radio wavelengths. We find a lower limit for the post-shock magnetic-field strength similar to 330 mu G, implying efficient amplification even for the magnetic-field damping scenario. Magnetic-field damping is necessary for the formation of the filaments in the radio range, while the X-ray filaments are shaped by both the synchrotron losses and magnetic-field damping.}, language = {en} } @phdthesis{Ehlert2023, author = {Ehlert, Kristian}, title = {Simulations of active galactic nuclei feedback with cosmic rays and magnetic fields}, doi = {10.25932/publishup-57816}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-578168}, school = {Universit{\"a}t Potsdam}, pages = {155}, year = {2023}, abstract = {The central gas in half of all galaxy clusters shows short cooling times. Assuming unimpeded cooling, this should lead to high star formation and mass cooling rates, which are not observed. Instead, it is believed that condensing gas is accreted by the central black hole that powers an active galactic nuclei jet, which heats the cluster. The detailed heating mechanism remains uncertain. A promising mechanism invokes cosmic ray protons that scatter on self-generated magnetic fluctuations, i.e. Alfv{\´e}n waves. Continuous damping of Alfv{\´e}n waves provides heat to the intracluster medium. Previous work has found steady state solutions for a large sample of clusters where cooling is balanced by Alfv{\´e}nic wave heating. To verify modeling assumptions, we set out to study cosmic ray injection in three-dimensional magnetohydrodynamical simulations of jet feedback in an idealized cluster with the moving-mesh code arepo. We analyze the interaction of jet-inflated bubbles with the turbulent magnetized intracluster medium. Furthermore, jet dynamics and heating are closely linked to the largely unconstrained jet composition. Interactions of electrons with photons of the cosmic microwave background result in observational signatures that depend on the bubble content. Those recent observations provided evidence for underdense bubbles with a relativistic filling while adopting simplifying modeling assumptions for the bubbles. By reproducing the observations with our simulations, we confirm the validity of their modeling assumptions and as such, confirm the important finding of low-(momentum) density jets. In addition, the velocity and magnetic field structure of the intracluster medium have profound consequences for bubble evolution and heating processes. As velocity and magnetic fields are physically coupled, we demonstrate that numerical simulations can help link and thereby constrain their respective observables. Finally, we implement the currently preferred accretion model, cold accretion, into the moving-mesh code arepo and study feedback by light jets in a radiatively cooling magnetized cluster. While self-regulation is attained independently of accretion model, jet density and feedback efficiencies, we find that in order to reproduce observed cold gas morphology light jets are preferred.}, language = {en} } @phdthesis{Werhahn2023, author = {Werhahn, Maria}, title = {Simulating galaxy evolution with cosmic rays: the multi-frequency view}, doi = {10.25932/publishup-57285}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-572851}, school = {Universit{\"a}t Potsdam}, pages = {5, 220}, year = {2023}, abstract = {Cosmic rays (CRs) constitute an important component of the interstellar medium (ISM) of galaxies and are thought to play an essential role in governing their evolution. In particular, they are able to impact the dynamics of a galaxy by driving galactic outflows or heating the ISM and thereby affecting the efficiency of star-formation. Hence, in order to understand galaxy formation and evolution, we need to accurately model this non-thermal constituent of the ISM. But except in our local environment within the Milky Way, we do not have the ability to measure CRs directly in other galaxies. However, there are many ways to indirectly observe CRs via the radiation they emit due to their interaction with magnetic and interstellar radiation fields as well as with the ISM. In this work, I develop a numerical framework to calculate the spectral distribution of CRs in simulations of isolated galaxies where a steady-state between injection and cooling is assumed. Furthermore, I calculate the non-thermal emission processes arising from the modelled CR proton and electron spectra ranging from radio wavelengths up to the very high-energy gamma-ray regime. I apply this code to a number of high-resolution magneto-hydrodynamical (MHD) simulations of isolated galaxies, where CRs are included. This allows me to study their CR spectra and compare them to observations of the CR proton and electron spectra by the Voyager-1 satellite and the AMS-02 instrument in order to reveal the origin of the measured spectral features. Furthermore, I provide detailed emission maps, luminosities and spectra of the non-thermal emission from our simulated galaxies that range from dwarfs to Milk-Way analogues to starburst galaxies at different evolutionary stages. I successfully reproduce the observed relations between the radio and gamma-ray luminosities with the far-infrared (FIR) emission of star-forming (SF) galaxies, respectively, where the latter is a good tracer of the star-formation rate. I find that highly SF galaxies are close to the limit where their CR population would lose all of their energy due to the emission of radiation, whereas CRs tend to escape low SF galaxies more quickly. On top of that, I investigate the properties of CR transport that are needed in order to match the observed gamma-ray spectra. Furthermore, I uncover the underlying processes that enable the FIR-radio correlation (FRC) to be maintained even in starburst galaxies and find that thermal free-free-emission naturally explains the observed radio spectra in SF galaxies like M82 and NGC 253 thus solving the riddle of flat radio spectra that have been proposed to contradict the observed tight FRC. Lastly, I scrutinise the steady-state modelling of the CR proton component by investigating for the first time the influence of spectrally resolved CR transport in MHD simulations on the hadronic gamma-ray emission of SF galaxies revealing new insights into the observational signatures of CR transport both spectrally and spatially.}, language = {en} } @phdthesis{Thomas2022, author = {Thomas, Timon}, title = {Cosmic-ray hydrodynamics: theory, numerics, applications}, doi = {10.25932/publishup-56384}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-563843}, school = {Universit{\"a}t Potsdam}, pages = {334}, year = {2022}, abstract = {Cosmic rays (CRs) are a ubiquitous and an important component of astrophysical environments such as the interstellar medium (ISM) and intracluster medium (ICM). Their plasma physical interactions with electromagnetic fields strongly influence their transport properties. Effective models which incorporate the microphysics of CR transport are needed to study the effects of CRs on their surrounding macrophysical media. Developing such models is challenging because of the conceptional, length-scale, and time-scale separation between the microscales of plasma physics and the macroscales of the environment. Hydrodynamical theories of CR transport achieve this by capturing the evolution of CR population in terms of statistical moments. In the well-established one-moment hydrodynamical model for CR transport, the dynamics of the entire CR population are described by a single statistical quantity such as the commonly used CR energy density. In this work, I develop a new hydrodynamical two-moment theory for CR transport that expands the well-established hydrodynamical model by including the CR energy flux as a second independent hydrodynamical quantity. I detail how this model accounts for the interaction between CRs and gyroresonant Alfv{\´e}n waves. The small-scale magnetic fields associated with these Alfv{\´e}n waves scatter CRs which fundamentally alters CR transport along large-scale magnetic field lines. This leads to the effects of CR streaming and diffusion which are both captured within the presented hydrodynamical theory. I use an Eddington-like approximation to close the hydrodynamical equations and investigate the accuracy of this closure-relation by comparing it to high-order approximations of CR transport. In addition, I develop a finite-volume scheme for the new hydrodynamical model and adapt it to the moving-mesh code Arepo. This scheme is applied using a simulation of a CR-driven galactic wind. I investigate how CRs launch the wind and perform a statistical analysis of CR transport properties inside the simulated circumgalactic medium (CGM). I show that the new hydrodynamical model can be used to explain the morphological appearance of a particular type of radio filamentary structures found inside the central molecular zone (CMZ). I argue that these harp-like features are synchrotron-radiating CRs which are injected into braided magnetic field lines by a point-like source such as a stellar wind of a massive star or a pulsar. Lastly, I present the finite-volume code Blinc that uses adaptive mesh refinement (AMR) techniques to perform simulations of radiation and magnetohydrodynamics (MHD). The mesh of Blinc is block-structured and represented in computer memory using a graph-based approach. I describe the implementation of the mesh graph and how a diffusion process is employed to achieve load balancing in parallel computing environments. Various test problems are used to verify the accuracy and robustness of the employed numerical algorithms.}, language = {en} } @article{WarrenEllisonBarkovetal.2017, author = {Warren, Donald C. and Ellison, Donald C. and Barkov, Maxim V. and Nagataki, Shigehiro}, title = {Nonlinear Particle Acceleration and Thermal Particles in GRB Afterglows}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {835}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa56c3}, pages = {15}, year = {2017}, abstract = {The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, N(E) proportional to E-p for p greater than or similar to 2. However, in standard Fermi shock acceleration, a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future gamma-ray observatories, such as CTA, very likely.}, language = {en} } @article{AceroAloisioAmansetal.2017, author = {Acero, F. and Aloisio, R. and Amans, J. and Amato, Elena and Antonelli, L. A. and Aramo, C. and Armstrong, T. and Arqueros, F. and Asano, Katsuaki and Ashley, M. and Backes, M. and Balazs, C. and Balzer, A. and Bamba, Aya and Barkov, Maxim and Barrio, J. A. and Benbow, Wystan and Bernloehr, K. and Beshley, V. and Bigongiari, C. and Biland, A. and Bilinsky, A. and Bissaldi, Elisabetta and Biteau, J. and Blanch, O. and Blasi, P. and Blazek, J. and Boisson, C. and Bonanno, G. and Bonardi, A. and Bonavolonta, C. and Bonnoli, G. and Braiding, C. and Brau-Nogue, S. and Bregeon, J. and Brown, A. M. and Bugaev, V. and Bulgarelli, A. and Bulik, T. and Burton, Michael and Burtovoi, A. and Busetto, G. and Bottcher, M. and Cameron, R. and Capalbi, M. and Caproni, Anderson and Caraveo, P. and Carosi, R. and Cascone, E. and Cerruti, M. and Chaty, Sylvain and Chen, A. and Chen, X. and Chernyakova, M. and Chikawa, M. and Chudoba, J. and Cohen-Tanugi, J. and Colafrancesco, S. and Conforti, V. and Contreras, J. L. and Costa, A. and Cotter, G. and Covino, Stefano and Covone, G. and Cumani, P. and Cusumano, G. and Daniel, M. and Dazzi, F. and De Angelis, A. and De Cesare, G. and De Franco, A. and De Frondat, F. and Dal Pino, E. M. de Gouveia and De Lisio, C. and Lopez, R. de los Reyes and De Lotto, B. and de Naurois, M. and De Palma, F. and Del Santo, M. and Delgado, C. and della Volpe, D. and Di Girolamo, T. and Di Giulio, C. and Di Pierro, F. and Di Venere, L. and Doro, M. and Dournaux, J. and Dumas, D. and Dwarkadas, Vikram V. and Diaz, C. and Ebr, J. and Egberts, Kathrin and Einecke, S. and Elsaesser, D. and Eschbach, S. and Falceta-Goncalves, D. and Fasola, G. and Fedorova, E. and Fernandez-Barral, A. and Ferrand, Gilles and Fesquet, M. and Fiandrini, E. and Fiasson, A. and Filipovic, Miroslav D. and Fioretti, V. and Font, L. and Fontaine, Gilles and Franco, F. J. and Freixas Coromina, L. and Fujita, Yutaka and Fukui, Y. and Funk, S. and Forster, A. and Gadola, A. and Lopez, R. Garcia and Garczarczyk, M. and Giglietto, N. and Giordano, F. and Giuliani, A. and Glicenstein, J. and Gnatyk, R. and Goldoni, P. and Grabarczyk, T. and Graciani, R. and Graham, J. and Grandi, P. and Granot, Jonathan and Green, A. J. and Griffiths, S. and Gunji, S. and Hakobyan, H. and Hara, S. and Hassan, T. and Hayashida, M. and Heller, M. and Helo, J. C. and Hinton, J. and Hnatyk, B. and Huet, J. and Huetten, M. and Humensky, T. B. and Hussein, M. and Horandel, J. and Ikeno, Y. and Inada, T. and Inome, Y. and Inoue, S. and Inoue, T. and Inoue, Y. and Ioka, K. and Iori, Maurizio and Jacquemier, J. and Janecek, P. and Jankowsky, D. and Jung, I. and Kaaret, P. and Katagiri, H. and Kimeswenger, S. and Kimura, Shigeo S. and Knodlseder, J. and Koch, B. and Kocot, J. and Kohri, K. and Komin, N. and Konno, Y. and Kosack, K. and Koyama, S. and Kraus, Michaela and Kubo, Hidetoshi and Mezek, G. Kukec and Kushida, J. and La Palombara, N. and Lalik, K. and Lamanna, G. and Landt, H. and Lapington, J. and Laporte, P. and Lee, S. and Lees, J. and Lefaucheur, J. and Lenain, J. -P. and Leto, Giuseppe and Lindfors, E. and Lohse, T. and Lombardi, S. and Longo, F. and Lopez, M. and Lucarelli, F. and Luque-Escamilla, Pedro Luis and Lopez-Coto, R. and Maccarone, M. C. and Maier, G. and Malaguti, G. and Mandat, D. and Maneva, G. and Mangano, S. and Marcowith, Alexandre and Marti, J. and Martinez, M. and Martinez, G. and Masuda, S. and Maurin, G. and Maxted, N. and Melioli, Claudio and Mineo, T. and Mirabal, N. and Mizuno, T. and Moderski, R. and Mohammed, M. and Montaruli, T. and Moralejo, A. and Mori, K. and Morlino, G. and Morselli, A. and Moulin, Emmanuel and Mukherjee, R. and Mundell, C. and Muraishi, H. and Murase, Kohta and Nagataki, Shigehiro and Nagayoshi, T. and Naito, T. and Nakajima, D. and Nakamori, T. and Nemmen, R. and Niemiec, Jacek and Nieto, D. and Nievas-Rosillo, M. and Nikolajuk, M. and Nishijima, K. and Noda, K. and Nogues, L. and Nosek, D. and Novosyadlyj, B. and Nozaki, S. and Ohira, Yutaka and Ohishi, M. and Ohm, S. and Okumura, A. and Ong, R. A. and Orito, R. and Orlati, A. and Ostrowski, M. and Oya, I. and Padovani, Marco and Palacio, J. and Palatka, M. and Paredes, Josep M. and Pavy, S. and Persic, M. and Petrucci, P. and Petruk, Oleh and Pisarski, A. and Pohl, Martin and Porcelli, A. and Prandini, E. and Prast, J. and Principe, G. and Prouza, M. and Pueschel, Elisa and Puelhofer, G. and Quirrenbach, A. and Rameez, M. and Reimer, O. and Renaud, M. and Ribo, M. and Rico, J. and Rizi, V. and Rodriguez, J. and Fernandez, G. Rodriguez and Rodriguez Vazquez, J. J. and Romano, Patrizia and Romeo, G. and Rosado, J. and Rousselle, J. and Rowell, G. and Rudak, B. and Sadeh, I. and Safi-Harb, S. and Saito, T. and Sakaki, N. and Sanchez, D. and Sangiorgi, P. and Sano, H. and Santander, M. and Sarkar, S. and Sawada, M. and Schioppa, E. J. and Schoorlemmer, H. and Schovanek, P. and Schussler, F. and Sergijenko, O. and Servillat, M. and Shalchi, A. and Shellard, R. C. and Siejkowski, H. and Sillanpaa, A. and Simone, D. and Sliusar, V. and Sol, H. and Stanic, S. and Starling, R. and Stawarz, L. and Stefanik, S. and Stephan, M. and Stolarczyk, T. and Szanecki, M. and Szepieniec, T. and Tagliaferri, G. and Tajima, H. and Takahashi, M. and Takeda, J. and Tanaka, M. and Tanaka, S. and Tejedor, L. A. and Telezhinsky, Igor O. and Temnikov, P. and Terada, Y. and Tescaro, D. and Teshima, M. and Testa, V. and Thoudam, S. and Tokanai, F. and Torres, D. F. and Torresi, E. and Tosti, G. and Townsley, C. and Travnicek, P. and Trichard, C. and Trifoglio, M. and Tsujimoto, S. and Vagelli, V. and Vallania, P. and Valore, L. and van Driel, W. and van Eldik, C. and Vandenbroucke, Justin and Vassiliev, V. and Vecchi, M. and Vercellone, Stefano and Vergani, S. and Vigorito, C. and Vorobiov, S. and Vrastil, M. and Vazquez Acosta, M. L. and Wagner, S. J. and Wagner, R. and Wakely, S. P. and Walter, R. and Ward, J. E. and Watson, J. J. and Weinstein, A. and White, M. and White, R. and Wierzcholska, A. and Wilcox, P. and Williams, D. A. and Wischnewski, R. and Wojcik, P. and Yamamoto, T. and Yamamoto, H. and Yamazaki, Ryo and Yanagita, S. and Yang, L. and Yoshida, T. and Yoshida, M. and Yoshiike, S. and Yoshikoshi, T. and Zacharias, M. and Zampieri, L. and Zanin, R. and Zavrtanik, M. and Zavrtanik, D. and Zdziarski, A. and Zech, Alraune and Zechlin, Hannes and Zhdanov, V. and Ziegler, A. and Zorn, J.}, title = {Prospects for Cherenkov Telescope Array Observations of the Young Supernova Remnant RX J1713.7-3946}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {840}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa6d67}, pages = {14}, year = {2017}, abstract = {We perform simulations for future Cherenkov Telescope Array (CTA) observations of RX J1713.7-3946, a young supernova remnant (SNR) and one of the brightest sources ever discovered in very high energy (VHE) gamma rays. Special attention is paid to exploring possible spatial (anti) correlations of gamma rays with emission at other wavelengths, in particular X-rays and CO/H I emission. We present a series of simulated images of RX J1713.7-3946 for CTA based on a set of observationally motivated models for the gamma-ray emission. In these models, VHE gamma rays produced by high-energy electrons are assumed to trace the nonthermal X-ray emission observed by XMM-Newton, whereas those originating from relativistic protons delineate the local gas distributions. The local atomic and molecular gas distributions are deduced by the NANTEN team from CO and H I observations. Our primary goal is to show how one can distinguish the emission mechanism(s) of the gamma rays (i.e., hadronic versus leptonic, or a mixture of the two) through information provided by their spatial distribution, spectra, and time variation. This work is the first attempt to quantitatively evaluate the capabilities of CTA to achieve various proposed scientific goals by observing this important cosmic particle accelerator.}, language = {en} } @article{delValleMuellerRomero2018, author = {del Valle, Maria Victoria and M{\"u}ller, A. L. and Romero, G. E.}, title = {High-energy radiation from collisions of high-velocity clouds and the Galactic disc}, series = {Monthly notices of the Royal Astronomical Society}, volume = {475}, journal = {Monthly notices of the Royal Astronomical Society}, number = {4}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx2984}, pages = {4298 -- 4308}, year = {2018}, abstract = {High-velocity clouds (HVCs) are interstellar clouds of atomic hydrogen that do not follow normal Galactic rotation and have velocities of a several hundred kilometres per second. A considerable number of these clouds are falling down towards the Galactic disc. HVCs form large and massive complexes, so if they collide with the disc a great amount of energy would be released into the interstellar medium. The cloud-disc interaction produces two shocks: one propagates through the cloud and the other through the disc. The properties of these shocks depend mainly on the cloud velocity and the disc-cloud density ratio. In this work, we study the conditions necessary for these shocks to accelerate particles by diffusive shock acceleration and we study the non-thermal radiation that is produced. We analyse particle acceleration in both the cloud and disc shocks. Solving a time-dependent two-dimensional transport equation for both relativistic electrons and protons, we obtain particle distributions and non-thermal spectral energy distributions. In a shocked cloud, significant synchrotron radio emission is produced along with soft gamma rays. In the case of acceleration in the shocked disc, the non-thermal radiation is stronger; the gamma rays, of leptonic origin, might be detectable with current instruments. A large number of protons are injected into the Galactic interstellar medium, and locally exceed the cosmic ray background. We conclude that under adequate conditions the contribution from HVC-disc collisions to the galactic population of relativistic particles and the associated extended non-thermal radiation might be important.}, language = {en} }