Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Efficient Fokker-Planck Solver and its Application to Stochastic Particle Acceleration in Galaxy Clusters
428
0
0.0
(
0
)
Authors
Julius Donnert
Ask ChatGPT about the research
No Arabic abstract
Particle acceleration by turbulence plays a role in many astrophysical environments. The non- linear evolution of the underlying cosmic-ray spectrum is complex and can be described by a Fokker-Planck equation, which in general has to be solved numerically. We present here an implementation to compute the evolution of a cosmic-ray spectrum coupled to turbulence considering isotropic particle pitch-angle distributions and taking into account the relevant particle energy gains and losses. Our code can be used in run time and post-processing to very large astrophysical fluid simulations. We also propose a novel method to compress cosmic- ray spectra by a factor of ten, to ease the memory demand in very large simulations. We show a number of code tests, which firmly establish the correctness of the code. In this paper we focus on relativistic electrons, but our code and methods can be easily extended to the case of hadrons. We apply our pipeline to the relevant problem of particle acceleration in galaxy clusters. We define a sub-grid model for compressible MHD-turbulence in the intra- cluster-medium and calculate the corresponding reacceleration timescale from first principles. We then use a magneto-hydrodynamic simulation of an isolated cluster merger to follow the evolution of relativistic electron spectra and radio emission generated from the system over several Gyrs.
rate research
Read More
Magnetic Tunnel Junctions (MTJs) constitute the novel memory element in STT-MRAM, which is ramping to production at major foundries as an eFlash replacement. MTJ switching exhibits a stochastic behavior due to thermal fluctuations, which is modeled by s-LLGS and Fokker-Planck (FP) equations. This work implements and benchmarks Finite Volume Method (FVM) and analytical solvers for the FP equation. To deploy an MTJ model for circuit design, it must be calibrated against silicon data. To address this challenge, this work presents a regression scheme to fit MTJ parameters to a given set of measured current, switching time and error rate data points, yielding a silicon-calibrated model suitable for MRAM macro transient simulation.
Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected throughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{rm radio}proptosigma_v^{3.3pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.
We present MultiModeCode, a Fortran 95/2000 package for the numerical exploration of multifield inflation models. This program facilitates efficient Monte Carlo sampling of prior probabilities for inflationary model parameters and initial conditions and is the first publicly available code that can efficiently generate large sample-sets for inflation models with $mathcal O(100)$ fields. The code numerically solves the equations of motion for the background and first-order perturbations of multi-field inflation models with canonical kinetic terms and arbitrary potentials, providing the adiabatic, isocurvature, and tensor power spectra at the end of inflation. For models with sum-separable potentials MultiModeCode also computes the slow-roll prediction via the $delta N$ formalism for easy model exploration and validation. We pay particular attention to the isocurvature perturbations as the system approaches the adiabatic limit, showing how to avoid numerical instabilities that affect some other approaches to this problem. We demonstrate the use of MultiModeCode by exploring a few toy models. Finally, we give a concise review of multifield perturbation theory and a users manual for the program.
We use the Matryoshka run to study the time dependent statistics of structure-formation driven turbulence in the intracluster medium of a 10$^{15}M_odot$ galaxy cluster. We investigate the turbulent cascade in the inner Mpc for both compressional and incompressible velocity components. The flow maintains approximate conditions of fully developed turbulence, with departures thereof settling in about an eddy-turnover-time. Turbulent velocity dispersion remains above $700$ km s$^{-1}$ even at low mass accretion rate, with the fraction of compressional energy between 10% and 40%. Normalisation and slope of compressional turbulence is susceptible to large variations on short time scales, unlike the incompressible counterpart. A major merger occurs around redshift $zsimeq0$ and is accompanied by a long period of enhanced turbulence, ascribed to temporal clustering of mass accretion related to spatial clustering of matter. We test models of stochastic acceleration by compressional modes for the origin of diffuse radio emission in galaxy clusters. The turbulence simulation model constrains an important unknown of this complex problem and brings forth its dependence on the elusive micro-physics of the intracluster plasma. In particular, the specifics of the plasma collisionality and the dissipation physics of weak shocks affect the cascade of compressional modes with strong impact on the acceleration rates. In this context radio halos emerge as complex phenomena in which a hierarchy of processes acting on progressively smaller scales are at work. Stochastic acceleration by compressional modes implies statistical correlation of radio power and spectral index with merging cores distance, both testable in principle with radio surveys.
We derive the stochastic description of a massless, interacting scalar field in de Sitter space directly from the quantum theory. This is done by showing that the density matrix for the effective theory of the long wavelength fluctuations of the field obeys a quantum version of the Fokker-Planck equation. This equation has a simple connection with the standard Fokker-Planck equation of the classical stochastic theory, which can be generalised to any order in perturbation theory. We illustrate this formalism in detail for the theory of a massless scalar field with a quartic interaction.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
