A Monte Carlo simulation study for cosmic-ray chemical composition measurement with Cherenkov Telescope Array

Our Galaxy is filled with cosmic-ray particles and more than 98% of them are atomic nuclei. In order to clarify their origin and acceleration mechanism, chemical composition measurements of these cosmic rays with wide energy coverage play an important role. Imaging Atmospheric Cherenkov Telescope (IACT) arrays are designed to detect cosmic gamma-rays in the very-high-energy regime ($\sim$TeV). Recently these systems proved to be capable of measuring cosmic-ray chemical composition in the sub-PeV region by capturing direct Cherenkov photons emitted by charged primary particles. Extensive air shower profiles measured by IACTs also contain information about the primary particle type since the cross section of inelastic scattering in the air depends on the primary mass number. The Cherenkov Telescope Array (CTA) is the next generation IACT system, which will consist of multiple types of telescopes and have a km$^2$-scale footprint and extended energy coverage (20 GeV to 300 TeV). In order to estimate CTA potential for cosmic ray composition measurement, a full Monte Carlo simulation including a description of extensive air shower and detector response is needed. We generated a number of cosmic-ray nuclei events (8 types selected from H to Fe) for a specific CTA layout candidate in the southern-hemisphere site. We applied Direct Cherenkov event selection and shower profile analysis to these data and preliminary results on charge number resolution and expected event count rate for these cosmic-ray nuclei are presented.Comment: All CTA contributions at arXiv:1709.03483 , Proceedings of the 35th International Cosmic Ray Conferenc

Approximate Solutions of Mathematical Models of Supercooling Solidification

In [1], [2], some one-dimensional mathematical models of supercooling solidification have been established, and some existence theorems have been proven by a difference method. In this paper, the models and the method are shown again in §§ 1, 2 and another analytic method of approximate solution is proposed in § 3. It is based on the assumption that a profile of temperature distribution at any time may be considered linear in every inner region. By some numerical examples, solutions by the approximate method are compared with some solutions given by the difference method. It is then realized that the approximate solutions come sufficiently close to the difference solutions

Electron acceleration with improved Stochastic Differential Equation method: cutoff shape of electron distribution in test-particle limit

We develop a method of stochastic differential equation to simulate electron acceleration at astrophysical shocks. Our method is based on It\^{o}'s stochastic differential equations coupled with a particle splitting, employing a skew Brownian motion where an asymmetric shock crossing probability is considered. Using this code, we perform simulations of electron acceleration at stationary plane parallel shock with various parameter sets, and studied how the cutoff shape, which is characterized by cutoff shape parameter $a$, changes with the momentum dependence of the diffusion coefficient $\beta$. In the age-limited cases, we reproduce previous results of other authors, $a\approx2\beta$. In the cooling-limited cases, the analytical expectation $a\approx\beta+1$ is roughly reproduced although we recognize deviations to some extent. In the case of escape-limited acceleration, numerical result fits analytical stationary solution well, but deviates from the previous asymptotic analytical formula $a\approx\beta$.Comment: corrected typos, 10 pages, 4 figures, 2 tables, JHEAp in pres