Search

Context. Star-forming regions are gaining considerable interest in the high-energy astrophysics community as possible Galactic particle accelerators. In general, the role of electrons has not been fully considered in this kind of cosmic-ray source. However, the intense radiation fields inside these regions might make electrons significant gamma-ray contributors. Aims. We study the young and compact star-forming region NGC 3603, a well known gamma-ray emitter. Our intention is to test whether its gamma-ray emission can be produced by cosmic-ray electrons. Methods. We build a novel model by creating realistic 3D distributions of the gas and the radiation field in the region. We introduce these models into PICARD to perform cosmic-ray transport simulations and produce gamma-ray emission maps. The results are compared with a dedicated Fermi Large Area Telescope data analysis at high energies. We also explore the radio and neutrino emissions of the system. Results. We improve the existing upper limits of the NGC 3603 gamma-ray source extension. Although the gamma-ray spectrum is well reproduced with the injection of CR protons, it requires nearly 30\% acceleration efficiency. In addition, the resulting extension of the simulated hadronic source is in mild tension with the extension data upper limit. The radio data disfavours the lepton-only scenario. Finally, combining both populations, the results are consistent with all observables, although the exact contributions are ambiguous., Comment: 12 pages, 10 figures

The spatial distribution of the gaseous components of the Milky Way is of great importance for a number of different fields, e.g. Galactic structure, star formation and cosmic rays. However, obtaining distance information to gaseous clouds in the interstellar medium from Doppler-shifted line emission is notoriously difficult given our unique vantage point in the Galaxy. It requires precise knowledge of gas velocities and generally suffers from distance ambiguities. Previous works often assumed the optically thin limit (no absorption), a fixed velocity field, and lack resolution overall. We aim to overcome these issues and improve previous reconstructions of the gaseous constituents of the interstellar medium of the Galaxy. We use 3D Gaussian processes to model correlations in the interstellar medium, including correlations between different lines of sight, and enforce a spatially coherent structure in the prior. For modelling the transport of radiation from the emitting gas to us as observers, we take absorption effects into account. A special numerical grid ensures high resolution nearby. We infer the spatial distributions of HI, CO, their emission line-widths, and the Galactic velocity field in a joint Bayesian inference. We further constrain these fields with complementary data from Galactic masers and young stellar object clusters. Our main result consists of a set of samples that implicitly contain statistical uncertainties. The resulting maps are spatially coherent and reproduce the data with high fidelity. We confirm previous findings regarding the warping and flaring of the Galactic disc. A comparison with 3D dust maps reveals a good agreement on scales larger than approximately 400 pc. While our results are not free of artefacts, they present a big step forward in obtaining high quality 3D maps of the interstellar medium., Comment: 21 pages, 22 figures

Cosmic rays (CRs) play a major role in the dynamics of the interstellar medium (ISM). Their interactions and transport ionize, heat, and push the ISM thereby coupling different regions of it. The spatial distribution of CRs depends on the distribution of their sources as well as the ISM constituents they interact with, such as gas, starlight, and magnetic fields. Particularly, gas interacts closely with CRs, influencing CR fluxes and gamma -ray emission. We illustrate the influence of 3D gas structures on CR transport and gamma -ray emission. We use the PICARD code and multiple samples of recent 3D reconstructions of the HI and H$_2$ Galactic gas constituents to investigate the impact on the transport of CRs and emission of gamma -rays. We find the necessary transport parameters to reproduce local measurements of CR fluxes, and see that they depend on the local distribution of gas density and structure. The distribution of CR fluxes exhibits energy-dependent structures that vary for all CR species due to their corresponding loss processes. Regions of enhanced secondary (primary) species are spatially correlated (anti-correlated) with the gas density. We observe a high sensitivity of the gamma -ray emission on the contrast of gas structures, as those determine the 3D spatial distributions of hadronic interactions and bremsstrahlung. We find that corresponding gas-induced structures in the distribution of CR electrons are also visible in Inverse Compton (IC) emission. Due to the aforementioned sensitivity, the analysis of CR data for CR sources and transport parameters requires the usage of accurate 3D gas maps.

Anisotropic diffusion is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and the interplay of cosmic rays with the Galactic magnetic field. This diffusion term contributes to the highly stiff nature of the cosmic ray transport equation. To conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators (namely, Crank-Nicolson (CN)) have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods to treat the linear anisotropc diffusion equation in the presence of advection and time-independent and time-dependent sources. These methods allow us to take even larger step sizes that can substantially speed-up the simulations whilst generating highly accurate solutions. In or subsequent work, we will use these exponential solvers in the Picard code to study anisotropic cosmic ray diffusion and we will consider additional physical processes such as continuous momentum losses and reacceleration., Comment: The 38th International Cosmic Ray Conference (ICRC2023)

Context. We present an analysis of our high-resolution relativistic-hydrodynamics model of the stellar- and pulsar-wind interaction in the LS-5039 system. Aims. With our high-resolution simulation covering three orbital periods, we analyse the impact of turbulence with a particular focus on short-term and orbit-to-orbit variations. Methods. Our model uses a relativistic hydrodynamics description of the wind interaction in the LS-5039 system assuming a pulsar-wind driven scenario. The corresponding system of equations is solved using the finite-volume code Cronos. We compute statistical quantities, also relevant for particle acceleration in this system, from results of multiple consecutive timesteps. Results. In our simulation we find the previously observed shock structures related to the wind-collision region (WCR), including the pulsar-wind termination, being dynamically influenced by orbital motion. In our high-resolution simulation we find high turbulence levels following from instabilities driven at the WCR. These instabilities lead to strong fluctuations of several dynamical quantities especially around and after apastron. These fluctuations are expected to impact the particle transport and also especially the related emission of non-thermal radiation. As an important example, the region from which gamma-ray emission has been found to be boosted due to relativistic beaming in previous studies shows strong variations in size both on short and on orbital timescales. Conclusions. Using a large computational domain together with high spatial resolution allowed a detailed study of fluctuations in the stellar- and pulsar-wind interaction. The results indicate a possible influence on the non-thermal emission from this system, which will be analysed with dedicated simulations in a forthcoming publication., Comment: 11 pages, 12 figures, Accepted by A&A

The anisotropic diffusion equation is imperative in understanding cosmic ray diffusion across the Galaxy, the heliosphere, and its interplay with the ambient magnetic field. This diffusion term contributes to the highly stiff nature of the CR transport equation. In order to conduct numerical simulations of time-dependent cosmic ray transport, implicit integrators have been traditionally favoured over the CFL-bound explicit integrators in order to be able to take large step sizes. We propose exponential methods that directly compute the exponential of the matrix to solve the linear anisotropic diffusion equation. These methods allow us to take even larger step sizes; in certain cases, we are able to choose a step size as large as the simulation time, i.e., only one time step. This can substantially speed-up the simulations whilst generating highly accurate solutions (l2 error $\leq 10^{-10}$). Additionally, we test an approach based on extracting a constant diffusion coefficient from the anisotropic diffusion equation, where the constant coefficient term is solved implicitly or exponentially and the remainder is treated using some explicit method. We find that this approach, for homogeneous linear problems, is unable to improve on the exponential-based methods that directly evaluate the matrix exponential., Comment: in submission

Cosmic rays (CRs) are an integral component of the interstellar medium, producing broadband emission while interacting with other Galactic matter components like the interstellar gas or magnetic fields. In addition to observations, numerical simulations of CR propagation through the Galaxy help to increase the level of understanding of Galactic CR transport and diffuse $\gamma$-ray emission as seen by different experiments. Up to now, the standard approach at modelling source distributions used as input for such transport simulations often rely on radial symmetry and analytical functions rather than individual, observation-based sources. We aim at a redefinition of existing CR source distributions by combining sources observed with the H.E.S.S. experiment and simulated random sources, which follow the matter density in the Milky Way. As a result, H.E.S.S.-inspired Galactic CR source distributions are inferred. We use the PICARD code to perform 3D-simulations of nuclei and electrons in CR propagation using our hybrid source distribution models. Furthermore, also gamma-ray maps and spectra, simulated with the redefined source models, are evaluated in different regions in the Galaxy and compared with each other to determine the statistical scatter of the underlying distributions. We find global consistency between our models and in comparison to previous simulations, with only some localised fluctuations, e.g. in the spiral arms. This implementation of a three-dimensional source model based on observations and simulations enables a new quality of propagation modelling. It offers possibilities for more realistic CR transport scenarios beyond radial symmetry and delivers meaningful results in both the arm and interarm regions of the Galaxy. This gives a more realistic picture of the Galactic $\gamma$-ray sky by including structures from the source model and not just the gas distributions., Comment: 32 pages, 14 figures, accepted for publication in Astroparticle Physics

We describe the special relativistic extension of the CRONOS code, which has been used for studies of gamma-ray binaries in recent years. The code was designed to be easily adaptable, allowing the user to easily change existing functionalities or introduce new modules tailored to the problem at hand. Numerically, the equations are treated using a finite-volume Godunov scheme on rectangular grids, which currently support Cartesian, spherical, and cylindrical coordinates. The employed reconstruction technique, the approximate Riemann solver and the equation of state can be chosen dynamically by the user. Further, the code was designed with stability and robustness in mind, detecting and mitigating possible failures early on. We demonstrate the code's capabilities on an extensive set of validation problems., Comment: 12 pages, 8 figures, accepted for publication in A&A

Context. We have presented a numerical model for the non-thermal emission of gamma-ray binaries in a pulsar-wind driven scenario. Aims. We apply this model to one of the best-observed gamma-ray binaries, the LS 5039 system. Methods. The model involves a joint simulation of the pulsar- and stellar-wind interaction and the transport of electronic pairs from the pulsar wind accelerated at the emerging shock structure. We compute the synchrotron and inverse Compton emission in a post-processing step, while consistently accounting for relativistic beaming and $\gamma \gamma$-absorption in the stellar radiation field. Results. The stellar- and pulsar-wind interaction leads to the formation of an extended, asymmetric wind collision region developing strong shocks, turbulent mixing, and secondary shocks in the turbulent flow. Both the structure of the collision region and the resulting particle distributions show significant orbital variation. Next to the acceleration of particles at the bow-like pulsar wind and Coriolis shock the model naturally accounts for the reacceleration of particles at secondary shocks contributing to the emission at very-high-energy (VHE) gamma-rays. The model successfully reproduces the main spectral features of LS 5039. While the predicted lightcurves in the high-energy and VHE gamma-ray band are in good agreement with observations, our model still does not reproduce the X-ray to low-energy gamma-ray modulation, which we attribute to the employed magnetic field model. Conclusions. We successfully model the main spectral features of the observed multiband, non-thermal emission of LS 5039 and thus further substantiates a wind-driven interpretation of gamma-ray binaries. Open issues relate to the synchrotron modulation, which might be addressed through a magnetohydrodynamic extension of our model., Comment: 14 pages, 13 figures. Accepted by A&A

Context. Gamma-ray binaries are systems that radiate the dominant part of their non-thermal emission in the gamma-ray band. In a wind-driven scenario, these binaries are thought to consist of a pulsar orbiting a massive star, accelerating particles in the shock arising in the wind collision. Aims. We develop a comprehensive, numerical model for the non-thermal emission of shock accelerated particles including the dynamical effects of fluid instabilities and orbital motion. We demonstrate the model on a generic binary system. Methods. The model is built on a dedicated three-dimensional particle transport simulation for the accelerated particles dynamically coupled to a simultaneous relativistic hydrodynamic simulation of the wind interaction. In a post-processing step, a leptonic emission model involving synchrotron and inverse Compton emission is evaluated based on resulting particle distributions and fluid solutions, consistently accounting for relativistic boosting and $\gamma\gamma$-absorption in the stellar radiation field. The model is implemented as an extension to the Cronos code. Results. In the generic binary, the wind interaction leads to the formation of an extended, asymmetric wind-collision region distorted by the effects of orbital motion, mixing and turbulence giving rise to strong shocks terminating the pulsar wind and secondary shocks in the turbulent fluid flow. With the presented approach it is, for the first time, possible to consistently account for the dynamical shock structure in particle transport processes, yielding a complex distribution of accelerated particles. The predicted emission extends over a broad region of energy, with significant orbital modulation in all bands., Comment: 12 pages, 5 figures. Accepted by A&A

The interaction between the strong winds in stellar colliding-wind binary (CWB) systems produces two shock fronts, delimiting the wind collision region (WCR). There, particles are expected to be accelerated mainly via diffusive shock acceleration (DSA). We investigate the injection and the acceleration of protons in typical CWB systems by means of Monte Carlo simulations, with both a test-particle approach and a non-linear method modelling a shock locally modified by the backreaction of the accelerated protons. We use magnetohydrodynamic simulations to determine the background plasma in the WCR and its vicinity. This allows us to consider particle acceleration at both shocks, on either side of the WCR, with a realistic large-scale magnetic field. We highlight the possible effects of particle acceleration on the local shock profiles at the WCR. We include the effect of magnetic field amplification due to resonant streaming instability (RSI), and compare results without and with the backreaction of the accelerated protons. In the latter case we find a lower flux of the non-thermal proton population, and a considerable magnetic field amplification. This would significantly increase the synchrotron losses of relativistic electrons accelerated in CWB systems, lowering the maximal energy they can reach and strongly reducing the inverse Compton fluxes. As a result, $\gamma$-rays from CWBs would be predominantly due to the decay of neutral pions produced in nucleon-nucleon collisions. This might provide a way to explain why, in the vast majority of cases, CWB systems have not been identified as $\gamma$-ray sources, while they emit synchrotron radiation.

We investigate the impact of the recently published ISRF model of the Milky Way by Popescu et al. on the CR electrons and positrons in the context of cosmic-ray transport modelling. We also study predictions for diffuse Galactic gamma radiation and underline the importance of the increased ISRF intensities.We use the PICARD code for solving the CR transport equation to obtain predictions based on Galactic cosmic-ray electron fluxes. We show that the new ISRF yields gamma-ray intensity increases, most particular at the Galactic Center and in the Galactic Plane. The impact is largest at energies around 220 GeV., Comment: 20 pages, 15 figures

The interaction between the strong winds of the stars in colliding-wind binary (CWB) systems produces two shock fronts, delimiting the wind collision region (WCR). There, particles are expected to be accelerated mainly via diffusive shock acceleration, and to produce $\gamma$-rays, in processes involving relativistic electrons and/or protons. We investigate the injection and the acceleration of protons in typical CWB systems by means of Monte Carlo simulations, with a test-particle approach. We use magnetohydrodynamic simulations to determine the background conditions in the wind collision region. This allows us to consider particle acceleration at both shocks, on either side of the WCR, with a self-consistently determined large-scale magnetic field, which has an impact on the shape of the WCR, and the topology of which plays an important role in particle acceleration at collisionless shocks. Such studies may contribute to improve $\gamma$-ray flux predictions for CWB systems.

An important contribution to the gamma-ray emission of our Galaxy is the Galactic Diffuse emission. We present specific developments within the PICARD code for modeling of Galactic cosmic ray propagation that are relevant for the computation of the gamma-ray emission to accommodate a new Interstellar Radiation Field (ISRF) model. We study the differences of the individual ISRF models on Galactic Diffuse gamma-ray emission predictions in the GeV to the TeV energy regime. Emphasis was laid on obtaining robust predictions for observable signatures in the very high energy gamma-ray regime with a special attention to the energy regime for H.E.S.S. Preliminary longitude-latitude profile and residuum maps of inverse Compton emissions are shown for axisymmetric and four-arm spiral Galaxy models.

We present results obtained from our newly developed Galactic cosmic-ray transport code PICARD, that solves the cosmic-ray transport equation. This code allows for the computation of cosmic-ray spectra and the resulting gamma-ray emission. Relying on contemporary numerical solvers allows for efficient computation of models with deca-parsec resolution. PICARD can handle locally anisotropic spatial diffusion acknowledging a full diffusion tensor. We used this framework to investigate the transition from axisymmetric to spiral-arm cosmic-ray source distributions. Wherever possible we compare model predictions with constraining observables in cosmic-ray astrophysics., Comment: 8 pages, 2 figures, presented at the ICRC2015 in Den Haag

We study the impact of possible spiral-arm distributions of Galactic cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our Galaxy. We investigate model cosmic-ray spectra at the nominal position of the sun and at different positions within the Galaxy. The modelling is performed using the recently introduced numerical cosmic ray propagation code \textsc{Picard}. Assuming non-axisymmetric cosmic ray source distributions yields new insights on the behaviour of primary versus secondary nuclei. We find that primary cosmic rays are more strongly confined to the vicinity of the sources, while the distribution of secondary cosmic rays is much more homogeneous compared to the primaries. This leads to stronger spatial variation in secondary to primary ratios when compared to axisymmetric source distribution models. A good fit to the cosmic-ray data at Earth can be accomplished in different spiral-arm models, although leading to decisively different spatial distributions of the cosmic-ray flux. This results in very different cosmic ray anisotropies, where even a good fit to the data becomes possible. Consequently, we advocate directions to seek best fit propagation parameters that take into account the higher complexity introduced by the spiral-arm structure on the cosmic-ray distribution. We specifically investigate whether the flux at Earth is representative for a large fraction of the Galaxy. The variance among possible spiral-arm models allows us to quantify the spatial variation of the cosmic-ray flux within the Galaxy in presence of non-axisymmetric source distributions., Comment: 38 pages, 16 figures, accepted for publication in Astroparticle Physics

We present the implementation of turbulence transport equations in addition to the Reynolds-averaged MHD equations within the Cronos framework. The model is validated by comparisons with earlier findings before it is extended to be applicable to regions in the solar wind that are not highly super-Alfv\'enic. We find that the respective additional terms result in absolute normalized cross-helicity to decline more slowly, while a proper implementation of the mixing terms can even lead to increased cross-helicities in the inner heliosphere. The model extension allows to place the inner boundary of the simulations closer to the Sun, where we choose its location at 0.1 AU for future application to the Wang-Sheeley-Arge model. Here, we concentrate on effects on the turbulence evolution for transient events by injecting a coronal mass ejection (CME). We find that the steep gradients and shocks associated with these structures result in enhanced turbulence levels and reduced cross-helicity. Our results can now be used straightforwardly for studying the transport of charged energetic particles, where the elements of the diffusion tensor can now benefit from the self-consistently computed solar wind turbulence. Furthermore, we find that there is no strong back-reaction of the turbulence on the large-scale flow so that CME studies concentrating on the latter need not be extended to include turbulence transport effects.

The transport of energetic particles such as Cosmic Rays is governed by the properties of the plasma being traversed. While these properties are rather poorly known for galactic and interstellar plasmas due to the lack of in situ measurements, the heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the 3D heliospheric transport of energetic particles are structures such as corotating interaction regions (CIRs), which, due to strongly enhanced magnetic field strengths, turbulence, and associated shocks, can act as diffusion barriers on the one hand, but also as accelerators of low energy CRs on the other hand as well. In a two-fold series of papers we investigate these effects by modeling inner-heliospheric solar wind conditions with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve as an input to a transport code employing a stochastic differential equation (SDE) approach (second paper). In this first paper we present results from 3D MHD simulations with our code CRONOS: for validation purposes we use analytic boundary conditions and compare with similar work by Pizzo. For a more realistic modeling of solar wind conditions, boundary conditions derived from synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where the potential field modeling is performed with a finite-difference approach (FDIPS) in contrast to the traditional spherical harmonics expansion often utilized in the WSA model. Our results are validated by comparing with multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses) regions.

In this manuscript we present a new approach for the numerical solution of the Galactic Cosmic Ray propagation problem. We introduce a method using advanced contemporary numerical algorithms while retaining the general complexity of other established codes. In this paper we present the underlying numerical scheme in conjunction with tests showing the correctness of the scheme. Finally we show the solution of a first example propagation problem using the new code to show its applicability to Galactic Cosmic Ray propagation.

A standard planetary nebula stays more than 10 000 years in the state of a photoionized nebula. As long as the timescales of the most important ionizing processes are much smaller, the ionization state can be characterized by a static photoionization model and simulated with codes like CLOUDY (Ferland et al. 1998). When the star exhibits a late Helium flash, however, its ionizing flux stops within a very short period. The star then re-appears from itsopaque shell after a few years (or centuries) as a cold giant star without any hard ionizing photons. Describing the physics of such behavior requires a fully time-dependent radiative transfer model. Pollacco (1999), Kerber et al. (1999) and Lechner & Kimeswenger (2004) used data of the old nebulae around V605 Aql and V4334 Sgr to derive a model of the pre-outburst state of the CSPN in a static model. Their argument was the long recombination time scale for such thin media. With regard to these models Schoenberner (2008) critically raised the question whether a significant change in the ionization state (and thus the spectrum) has to be expected after a time of up to 80 years, and whether static models are applicable at all., Comment: (3 pages, 1 figure, to appear in proceedings of the IAU Symposium 283: "Planetary Nebulae: An Eye to the Future", Eds.: A. Manchado, L. Stanghellini and D. Schoenberner; presenting author: Stefan Kimeswenger)