SUMRAY: R and Python codes for calculating cancer risk due to radiation exposure of a population

Sasaki, Michiya*; Furukawa, Kyoji*; Satoh, Daiki; Shimada, Kazumasa; Kudo, Shinichi*; Takagi, Shunji*; Takahara, Shogo; Kai, Michiaki*

Journal of Radiation Protection and Research, 48(2), p.90 - 99, 2023/06

https://doi.org/10.14407/jrpr.2022.00213

This paper reports on the calculation code that is the result of the activities of the "Task Group for Development of Cancer Risk Estimation Codes Associated with Radiation Exposure (FY2020-2021)" established by the Japan Health Physics Society. In order to promote research on the estimation of cancer risk associated with radiation exposure, the Task Group decided to disclose the source code, including the algorithm and parameters used in the calculations, and to release the code under a license that permits modification and redistribution of the code. The computational code was named SUMRAY and coded in two computer languages, that is R and Python. The code is capable of calculating the accumulated excess risk using Monte Carlo methods with a 95% confidence interval. The results of SUMRAY were compared with the results of the existing codes whose source code is not publicly available, under the same calculation conditions. From the results, it was found that they were in reasonable agreement within the confidence interval. It is expected that SUMRAY, an open-source software, will be used as a common basis for cancer risk estimation studies associated with radiation exposure.

Development of guidelines on radiation protection for the lens of the eye in Japan

Yokoyama, Sumi*; Iwai, Satoshi*; Tsujimura, Norio; Hashimoto, Makoto; Yoshitomi, Hiroshi; Kato, Masahiro*; Kurosawa, Tadahiro*; Tatsuzaki, Hideo*; Sekiguchi, Hiroshi*; Koguchi, Yasuhiro*; et al.

Proceedings of 15th International Congress of the International Radiation Protection Association (IRPA-15) (Internet), 8 Pages, 2022/00

https://irpa.net/irpa15/T3.1-O0124.pdf

Thermally altered subsurface material of asteroid (162173) Ryugu

Kitazato, Kohei*; Milliken, R. E.*; Iwata, Takahiro*; Abe, Masanao*; Otake, Makiko*; Matsuura, Shuji*; Takagi, Yasuhiko*; Nakamura, Tomoki*; Hiroi, Takahiro*; Matsuoka, Moe*; et al.

Nature Astronomy (Internet), 5(3), p.246 - 250, 2021/03

https://doi.org/10.1038/s41550-020-01271-2

Here we report observations of Ryugu's subsurface material by the Near-Infrared Spectrometer (NIRS3) on the Hayabusa2 spacecraft. Reflectance spectra of excavated material exhibit a hydroxyl (OH) absorption feature that is slightly stronger and peak-shifted compared with that observed for the surface, indicating that space weathering and/or radiative heating have caused subtle spectral changes in the uppermost surface. However, the strength and shape of the OH feature still suggests that the subsurface material experienced heating above 300 C, similar to the surface. In contrast, thermophysical modeling indicates that radiative heating does not increase the temperature above 200 C at the estimated excavation depth of 1 m, even if the semimajor axis is reduced to 0.344 au. This supports the hypothesis that primary thermal alteration occurred due to radiogenic and/or impact heating on Ryugu's parent body.

NIRS3; The Near Infrared Spectrometer on Hayabusa2

Iwata, Takahiro*; Kitazato, Kohei*; Abe, Masanao*; Otake, Makiko*; Arai, Takehiko*; Arai, Tomoko*; Hirata, Naru*; Hiroi, Takahiro*; Honda, Chikatoshi*; Imae, Naoya*; et al.

Space Science Reviews, 208(1-4), p.317 - 337, 2017/07

https://doi.org/10.1007/s11214-017-0341-0

NIRS3: The Near Infrared Spectrometer is installed on the Hayabusa2 spacecraft to observe the target C-type asteroid 162173 Ryugu at near infrared wavelengths of 1.8 to 3.2 micrometer. It aims to obtain reflectance spectra in order to detect absorption bands of hydrated and hydroxide minerals in the 3 micrometer-band. We adopted a linear-image sensor with indium arsenide (InAs) photo diodes and a cooling system with a passive radiator to achieve an optics temperature of 188 K, which enables to retaining sufficient sensitivity and noise level in the 3 micrometer wavelength region. We conducted ground performance tests for the NIRS3 flight model (FM) to confirm its baseline specifications. The results imply that the properties such as the signal-to-noise ratio (SNR) conform to scientific requirements to determine the degree of aqueous alteration, such as CM or CI chondrite, and the stage of thermal metamorphism on the asteroid surface.

Comparison between Monte Carlo simulation and measurement with a 3D polymer gel dosimeter for dose distributions in biological samples

Furuta, Takuya; Maeyama, Takuya*; Ishikawa, Kenichi*; Fukunishi, Nobuhisa*; Fukasaku, Kazuaki*; Takagi, Shu*; Noda, Shigeho*; Himeno, Ryutaro*; Hayashi, Shinichiro*

Physics in Medicine & Biology, 60(16), p.6531 - 6546, 2015/08

https://doi.org/10.1088/0031-9155/60/16/6531

Low reproducibility of dose distribution in inhomogeneous regions such as soft matter near bones is known with the simple dose analysis currently adopted in treatment planning of particle cancer therapy. Therefore a treatment planning system based on Monte Carlo simulation having better accuracy is highly desired. In order to assess the simulation accuracy of a Monte Carlo simulation code in situations closely related to medical application, we performed a comparison of dose distribution in a biological sample obtained by experiment and that by simulation. In particular, we irradiate a carbon beam on a biological sample composed of fresh chicken meat and bones, with a PAGAT gel dosimeter placed behind it, and compare the complex dose distribution in the gel dosimeter created by the beam passing through the inhomogeneous sample. Monte Carlo simulation using PHITS code was conducted by reconstructing the biological sample from its computed tomography images. The simulation accurately reproduced the experimental distal edge structure of the dose distribution with an accuracy under about 2 mm.

Parallel computing with Particle and Heavy Ion Transport code System (PHITS)

Furuta, Takuya; Sato, Tatsuhiko; Ogawa, Tatsuhiko; Niita, Koji*; Ishikawa, Kenichi*; Noda, Shigeho*; Takagi, Shu*; Maeyama, Takuya*; Fukunishi, Nobuhisa*; Fukasaku, Kazuaki*; et al.

Proceedings of Joint International Conference on Mathematics and Computation, Supercomputing in Nuclear Applications and the Monte Carlo Method (M&C + SNA + MC 2015) (CD-ROM), 9 Pages, 2015/04

In Particle and Heavy Ion Transport code System PHITS, two parallel computing functions are prepared to reduce the computational time. One is the distributed-memory parallelization using message passing interface (MPI) and the other is the shared-memory parallelization using OpenMP directives. Each function has advantages and disadvantages, and thus, by adopting both functions in PHITS, it is possible to conduct parallel computing suited for needs of users. It is also possible to conduct the hybrid parallelization by the intra-node OpenMP parallelization and the inter-node MPI parallelization in supercomputer systems. Each parallelization functions were explained together with some application results obtained using a workstation and a supercomputer system, K computer at RIKEN.

Beam test of a new radio frequency quadrupole linac for the Japan Proton Accelerator Research Complex

Kondo, Yasuhiro; Morishita, Takatoshi; Yamazaki, Saishun; Hori, Toshihiko; Sawabe, Yuki; Chishiro, Etsuji; Fukuda, Shimpei; Hasegawa, Kazuo; Hirano, Koichiro; Kikuzawa, Nobuhiro; et al.

Physical Review Special Topics; Accelerators and Beams, 17(12), p.120101_1 - 120101_8, 2014/12

https://doi.org/10.1103/PhysRevSTAB.17.120101

We performed a beam test of a new radio frequency quadrupole linac (RFQ III) for the beam current upgrade of the Japan Proton Accelerator Research Complex. First, the conditioning of RFQ III was conducted, and after 20 h of conditioning, RFQ III became very stable with a nominal peak power and duty factor of 400 kW and 1.5%, respectively. An off-line beam test was subsequently conducted before installation in the accelerator tunnel. The transmission, transverse emittance, and energy spread of the 50-mA negative hydrogen beam from RFQ III were measured and compared with simulation results. The experiment and simulation results showed good agreement; therefore, we conclude that the performance of RFQ III conforms to its design.

Present status of J-PARC linac

Oguri, Hidetomo; Hasegawa, Kazuo; Ito, Takashi; Chishiro, Etsuji; Hirano, Koichiro; Morishita, Takatoshi; Shinozaki, Shinichi; Ao, Hiroyuki; Okoshi, Kiyonori; Kondo, Yasuhiro; et al.

Proceedings of 11th Annual Meeting of Particle Accelerator Society of Japan (Internet), p.389 - 393, 2014/10

http://www.pasj.jp/web_publish/pasj2014/proceedings/PDF/FSP0/FSP021.pdf

no abstracts in English

Development of a Cs-ceeded Rf-driven H ion source for the J-PARC

Koizumi, Isao; Ueno, Akira; Okoshi, Kiyonori; Ikegami, Kiyoshi*; Takagi, Akira*; Yamazaki, Saishun; Oguri, Hidetomo

Proceedings of 10th Annual Meeting of Particle Accelerator Society of Japan (Internet), p.920 - 923, 2014/06

http://www.pasj.jp/web_publish/pasj10/proceedings/PDF/SUP0/SUP032.pdf

In order to satisfy the J-PARC (Japan Proton Accelerator Research Complex) 2nd stage requirements the J-PARC rf-driven H ion source is under development. The source plasma is produced by an about 50 W CW 30 MHz-rf and an about 40 kW pulsed 2 MHz-rf. The Cs is installed into the plasma chamber by opening the Cs-valve between the plasma chamber and the Cs oven, which is heated to 200240 C. The H ion beam current is maximized by optimizing the Cs quantity depending upon the Cs-valve opening time, the PE temperature depending upon the flow rate of the air through the PE cooling plate, the AMFC by the AMFC coil current, and so on. The experimental results using the plasma chamber made of copper, which satisfied the above requirements, were reported last year. The experimental results using the plasma chamber made of stainless steel are presented in this paper. The observed several interesting differences between the two plasma chambers are also presented.

Verification and validation procedures of calculation codes for determining corrosive conditions in the BWR primary cooling system based on water radiolysis and mixed potential models

Uchida, Shunsuke*; Wada, Yoichi*; Yamamoto, Seiji*; Takagi, Junichi*; Hisamune, Kenji*

Journal of Nuclear Science and Technology, 51(1), p.24 - 36, 2014/01

https://doi.org/10.1080/00223131.2013.842945

ECP in the BWR primary cooling system can be measured only in the restricted region. In order to determine ECP at any location, ECP should be evaluated by computer simulation codes consisting of water radiolysis models to determine the concentrations of corrosive radiolytic species and mixed potential models to determine ECP based on corrosive species. Mitigation of SCC crack growth rate due to decreasing ECP can be authorized by the JSME Standards, while mitigation of ECP due to hydrogen addition has not been authorized yet. In the paper, standard procedures to authorize the computer simulation codes based on the verification and validation method are proposed. The numerical justification of every code applied as a standard code should be verified and its accuracy and applicability for plant analysis should be validated. Benchmark problems for verification processes are proposed and comparison of the calculated results with the measured ones for the plant of evaluation is required.