Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Blazar Hadronic Code Comparison Project
266
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Blazar hadronic models have been developed in the past decades as an alternative to leptonic ones. In hadronic models the gamma-ray emission is associated with synchrotron emission by protons, and/or secondary leptons produced in proton-photon interactions. Together with photons, hadronic emission models predict the emission of neutrinos that are therefore the smoking gun for acceleration of relativistic hadrons in blazar jets. The simulation of proton-photon interactions and all associated radiative processes is a complex numerical task, and different approaches to the problem have been adopted in the literature. So far, no systematic comparison between the different codes has been performed, preventing a clear understanding of the underlying uncertainties in the numerical simulations. To fill this gap, we have undertaken the first comprehensive comparison of blazar hadronic codes, and the results from this effort will be presented in this contribution.
rate research
Read More
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.
Recently reported coincidences between high-energy neutrino events and major blazar outbursts reinforce the relevance of lepto-hadronic emission models for blazars. We study the influence of physical parameters on the neutrino output modeling blazar spectral energy distributions self-consistently assuming a relativistically propagating acceleration zone surrounded by a larger cooling zone. We find that the gross features of the spectral energy distribution can readily be explained with the model. A rigorous test requires time-resolved measurements of blazar spectral energy distributions during an outburst and high-statistics neutrino measurements to discriminate the leptonic and hadronic emission components.
The measured fluxes of secondary particles produced by the interactions of cosmic rays with the astronomical environment represent a powerful tool to infer some properties of primary cosmic rays. In this work we investigate the production of secondary particles in inelastic hadronic interactions between several cosmic rays species of projectiles and different target nuclei of the interstellar medium. The yields of secondary particles have been calculated with the FLUKA simulation package, that provides with very good accuracy the energy distributions of secondary products in a large energy range. An application to the propagation and production of secondaries in the Galaxy is presented.
Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $sim 10$% to $sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.
We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
