Simulation code for modeling of coherent effects of radiation generation in oriented crystals

Abstract

We present the CRYSTALRAD simulation code, combining all the features of the CRYSTAL simulation code for simulations of charged particles trajectories in a bent crystal and the RADCHARM++ code for calculation of the radiation spectrum. The CRYSTALRAD code is based on Monte Carlo simulations of trajectories in the planar and axial electric field either in a straight, bent, or periodically bent crystal taking into account multiple and single Coulomb scattering on nuclei and electrons, nuclear scattering and ionization energy losses. The trajectories simulated are used for calculation of radiation spectra by the Baier-Katkov method. We compare our simulations with experimental data taken at MAMI (MAinzer MIkrotron) as well as give an example for a possible future study with sub-GeV electrons interacting with Si bent crystals.

Improved experimental layout for dipole moment measurements at the LHC

“…Regarding the specifications for the two crystal lenses, the optimal crystal parameters for the plain-crystal scheme have been updated recently in [35] using state-of-the-art simulations for crystal channeling [41,42]. An optimal crystal with L C2 = 10 cm and θ C2 = 16 mrad was found to maximize the sensitivity to the electric and magnetic dipole moments.…”

MBN explorer atomistic simulations of 855 MeV electron propagation and radiation emission in oriented silicon bent crystal: theory versus experiment

“…[10].) This scheme has been implemented in a number of codes [25,[28][29][30][31][32][33][34]. Some of these codes allow to account for the effects, which are beyond the continuous potential model by including a stochastic force due to the random scattering of the projectile by lattice electrons and nuclei [30][31][32][33]35,36] and radiative damping force [32,33,37].…”

“…Within the framework of the quasi-classical method, the radiative recoil is accounted for making possible correct calculation of the spectral-angular distribution of the radiation in the whole photon energy range hω ε with exception for the extreme high-energy tail of the spectrum, (ε − hω)/ε 1. Quasi-classical computation of the emission spectrum is implemented in a number of aforementioned codes or their extensions [30,[32][33][34]40,44].…”

Layouts for fixed-target experiments and dipole moment measurements of short-lived baryons using bent crystals at the LHC

“…An agreement at the level of few % is also present between predictions reported in [55] and in [56], where the tools used in this paper and a yet different simulation approach are used, respectively. In particular, the latter is based on solving the equation of motion for a charged particle interacting with the electric field generated by crystalline lattice accounting also for multiple scattering [57][58][59].…”