The impact of dose rate on responses of human lens epithelial cells to ionizing irradiation.

The knowledge on responses of human lens epithelial cells (HLECs) to ionizing radiation exposure is important to understand mechanisms of radiation cataracts that are of concern in the field of radiation protection and radiation therapy. However, biological effects in HLECs following protracted exposure have not yet fully been explored. Here, we investigated the temporal kinetics of γ-H2AX foci as a marker for DNA double-strand breaks (DSBs) and cell survival in HLECs after exposure to photon beams at various dose rates (i.e., 150 kVp X-rays at 1.82, 0.1, and 0.033 Gy/min, and 137Cs γ-rays at 0.00461 Gy/min (27.7 cGy/h) and 0.00081 Gy/min (4.9 cGy/h)), compared to those in human lung fibroblasts (WI-38). In parallel, we quantified the recovery for DSBs and cell survival using a biophysical model. The study revealed that HLECs have a lower DSB repair rate than WI-38 cells. There is no significant impact of dose rate on cell survival in both cell lines in the dose-rate range of 0.033-1.82 Gy/min. In contrast, the experimental residual γ-H2AX foci showed inverse dose rate effects (IDREs) compared to the model prediction, highlighting the importance of the IDREs in evaluating radiation effects on the ocular lens.

Monte carlo simulation study on the dose and dose-averaged linear energy transfer distributions in carbon ion radiotherapy.

Dose-averaged linear energy transfer (LETd) is conventionally evaluated from the relative biological effectiveness (RBE)-LETd fitted function used in the treatment planning system. In this study, we calculated the physical doses and their linear energy transfer (LET) distributions for patterns of typical CIRT beams using Monte Carlo (MC) simulation. The LETd was then deduced from the MC simulation and compared with that obtained from the conventional method. The two types of LETd agreed well with each other, except around the distal end of the spread-out Bragg peak. Furthermore, an MC simulation was conducted with the material composition of water and realistic materials. The profiles of physical dose and LETd were in good agreement for both techniques. These results indicate that the previous studies to analyze the minimum LETd in CIRT cases are valid for practical situations, and the material composition conversion to water little affects the dose distribution in the irradiation field.

Comparison of dose and risk estimates between ISS Partner Agencies for a 30-day lunar mission.

The International Partner Agencies of the International Space Station (ISS) present a comparison of the ionizing radiation absorbed dose and risk quantities used to characterize example missions in lunar space. This effort builds on previous collaborative work that characterizes radiation environments in space to support radiation protection for human spaceflight on ISS in low-Earth orbit (LEO) and exploration missions beyond (BLEO). A "shielded" ubiquitous galactic cosmic radiation (GCR) environment combined with--and separate from--the transient challenge of a solar particle event (SPE) was modelled for a simulated 30-day mission period. Simple geometries of relatively thin and uniform shields were chosen to represent the space vehicle and other available shielding, and male or female phantoms were used to represent the body's self-shielding. Absorbed dose in organs and tissues and the effective dose were calculated for males and females. Risk parameters for cancer and other outcomes are presented for selected organs. The results of this intracomparison between ISS Partner Agencies itself provide insights to the level of agreement with which space agencies can perform organ dosimetry and calculate effective dose. This work was performed in collaboration with the advisory and guidance efforts of the International Commission on Radiological Protection (ICRP) Task Group 115 and will be presented in an ICRP Report.

Subfoveal perfluorocarbon liquid droplet can be removed without puncturing the retina.

Subfoveal perfluorocarbon liquid (PFL) is a vision-threatening complication that requires removal as safely as possible. We experienced a surgical case of proliferative vitreoretinopathy, in which a subfoveal PFL droplet was removed through the fovea, without puncturing the retina. In this case, although the retina was completely attached after primary vitrectomy and cataract removal, a subfoveal PFL droplet was found. At the second surgery, after peeling the internal limiting membrane, the droplet was passively aspirated by placing a 25-gauge blunt needle on the surface of the retina. This useful technique avoids puncturing the retina and does not require the patient to maintain a specific position postoperatively.

Fetal and Maternal Atomic Bomb Survivor Dosimetry Using the J45 Pregnant Female Phantom Series: Considerations of the Kneeling and Lying Posture with Comparisons to the DS02 System.

ABSTRACT: Organ dosimetry data of the atomic bomb survivors and the resulting cancer risk models derived from these data are currently assessed within the DS02 dosimetry system developed through the Joint US-Japan Dosimetry Working Group. In DS02, the anatomical survivor models are limited to three hermaphroditic stylized phantoms-an adult (55 kg), a child (19.8 kg), and an infant (9.7 kg)-that were originally designed for the preceding DS86 dosimetry system. As such, organ doses needed for assessment of in-utero cancer risks to the fetus have continued to rely upon the use of the uterine wall in the adult non-pregnant stylized phantom as the dose surrogate for all fetal organs regardless of gestational age. To address these limitations, the Radiation Effects Research Foundation (RERF) Working Group on Organ Dose (WGOD) has established the J45 (Japan 1945) series of high-resolution voxel phantoms, which were derived from the UF/NCI series of hybrid phantoms and scaled to match mid-1940s Japanese body morphometries. The series includes male and female phantoms-newborn to adult-and four pregnant female phantoms at gestational ages of 8, 15, 25, and 38 wk post-conception. In previous studies, we have reported organ dose differences between those reported by the DS02 system and those computed by the WGOD using 3D Monte Carlo radiation transport simulations of atomic bomb gamma-ray and neutron fields for the J45 phantoms series in their traditional "standing" posture, with some variations in their facing direction relative to the bomb hypocenter. In this present study, we present the J45 pregnant female phantoms in both a "kneeling" and "lying" posture and assess the dosimetric impact of these more anatomically realistic survivor models in comparison to current organ doses given by the DS02 system. For the kneeling phantoms facing the bomb hypocenter, organ doses from bomb source photon spectra were shown to be overestimated by the DS02 system by up to a factor of 1.45 for certain fetal organs and up to a factor of 1.17 for maternal organs. For lying phantoms with their feet in the direction of the hypocenter, fetal organ doses from bomb source photon spectra were underestimated by the DS02 system by factors as low as 0.77, while maternal organ doses were overestimated by up to a factor of 1.38. Organs doses from neutron contributions to the radiation fields exhibited an increasing overestimation by the DS02 stylized phantoms as gestational age increased. These discrepancies are most evident in fetal organs that are more posterior within the mother's womb, such as the fetal brain. Further analysis revealed that comparison of these postures to the original standing posture indicate significant dose differences for both maternal and fetal organ doses depending on the type of irradiation. Results from this study highlight the degree to which the existing DS02 system can differ from organ dosimetry based upon 3D radiation transport simulations using more anatomically realistic models of those survivors exposed during pregnancy.

Modelling oxygen effects on the in- and out-of-field radiosensitivity of cells exposed to intensity-modulated radiation fields.

Objective. The delivery of intensity-modulated radiation fields has improved the conformity of dose to tumour targets during radiotherapy (RT). Previously, it has been shown that intercellular communication between cells positioned in- and outside of the radiation field impacts cellular radiosensitivity under hypoxic and normoxic conditions. However, the mechanism of intercellular communication in hypoxia remains to be fully understood. In this study, the cell-killing effects of intercellular communication in hypoxia were modelled in an effort to better understand the underlying mechanisms of response.Approach. By irradiating a 50% area of the culture dish (half-field exposure), experimental dose-response curves for cell survival and residual DNA double-strand breaks (DSBs) were generated in prostate (DU145) and non-small cell lung cancer (H1299) cells. The oxygen enhancement ratio (OER) was determined from early DSB yields (corresponding to relative direct damage) and used to model the in- and out-of-field radiosensitivity.Main results. The developed integrated microdosimetric-kinetic (IMK) model successfully predicted the experimental dose responses for survival and lethal lesions, and provides a mechanistic interpretation that the probability of hits for releasing cell-killing signals is dependent on oxygen. This experimental and modelling study also suggests that residual DSBs correspond to logarithmic survival fraction (meaning lethal lesions) for in- and out-of-field cells. Our data suggest that the OER value determined using uniform-field exposure can be applied to predict the in- and out-of-field radiosensitivity of cells following exposure to intensity modulated beams.Significance. The developed IMK model facilitates a more precise understanding of intercellular signalling following exposure to intensity-modulated radiation fields.

Prospective comparison of accuracy of intraocular lens calculation formulas in phacovitrectomy: a pilot study in a real-world clinical practice.

PURPOSE: To compare the accuracy of intraocular lens (IOL) power calculations among IOL formulas after phacovitrectomy. METHODS: We prospectively enrolled 206 eyes of 206 patients who underwent 25-gauge phacovitrectomy, without gas tamponade, for macular pathology. Pre-operative optical biometry used the IOLMaster 700 to calculate the IOL power with the new formulas, i.e. the Barrett Universal II (BU II), Emmetropia Verifying Optical version 2.0, Hill-Radial Basis Function (RBF) version 3.0, Kane, and Ladas Super Formula, and conventional formulas, i.e. Haigis, Hoffer Q, Holladay 1, Holladay 2, and Sanders-Retzlaff-Kraff/T (SRK/T). A single-piece foldable IOL was implanted in all cases. Manifest refractions were measured before and 3 months after surgery. RESULTS: The BU II formula showed the lowest standard deviation and mean and median absolute errors and had the highest percentage of eyes with a refractive prediction error within ± 0.25 D. The absolute error was significantly lower with the four new formulas, except the Hill-RBF, than with the Hoffer Q (all p = ≤ 0.010) and Holladay 1 formulas (all p = < 0.010). The absolute error with the BU II formula was also lower than that with the Holladay 2 (p = 0.012) and SRK/T (p = 0.024) formulas. CONCLUSION: Overall, the new IOL formulas, except the Hill-RBF, were superior to some of the conventional formulas for calculating IOL power in phacovitrectomy.

Application of a simple DNA damage model developed for electrons to proton irradiation.

Proton beam therapy allows irradiating tumor volumes with reduced side effects on normal tissues with respect to conventional x-ray radiotherapy. Biological effects such as cell killing after proton beam irradiations depend on the proton kinetic energy, which is intrinsically related to early DNA damage induction. As such, DNA damage estimation based on Monte Carlo simulations is a research topic of worldwide interest. Such simulation is a mean of investigating the mechanisms of DNA strand break formations. However, past modellings considering chemical processes and DNA structures require long calculation times. Particle and heavy ion transport system (PHITS) is one of the general-purpose Monte Carlo codes that can simulate track structure of protons, meanwhile cannot handle radical dynamics simulation in liquid water. It also includes a simple model enabling the efficient estimation of DNA damage yields only from the spatial distribution of ionizations and excitations without DNA geometry, which was originally developed for electron track-structure simulations. In this study, we investigated the potential application of the model to protons without any modification. The yields of single-strand breaks, double-strand breaks (DSBs) and the complex DSBs were assessed as functions of the proton kinetic energy. The PHITS-based estimation showed that the DSB yields increased as the linear energy transfer (LET) increased, and reproduced the experimental and simulated yields of various DNA damage types induced by protons with LET up to about 30 keVµm-1. These results suggest that the current DNA damage model implemented in PHITS is sufficient for estimating DNA lesion yields induced after protons irradiation except at very low energies (below 1 MeV). This model contributes to evaluating early biological impacts in radiation therapy.

Comparison of the Therapeutic Effects of [211At]NaAt and [131I]NaI in an NIS-Expressing Thyroid Cancer Mouse Model.

Astatine (211At) is an alpha-emitter with a better treatment efficacy against differentiated thyroid cancer compared with iodine (131I), a conventional beta-emitter. However, its therapeutic comparison has not been fully evaluated. In this study, we compared the therapeutic effect between [211At]NaAt and [131I]NaI. In vitro analysis of a double-stranded DNA break (DSB) and colony formation assay were performed using K1-NIS cells. The therapeutic effect was compared using K1-NIS xenograft mice administered with [211At]NaAt (0.4 MBq (n = 7), 0.8 MBq (n = 9), and 1.2 MBq (n = 4)), and [131I]NaI (1 MBq (n = 4), 3 MBq (n = 4), and 8 MBq (n = 4)). The [211At]NaAt induced higher numbers of DSBs and had a more reduced colony formation than [131I]NaI. In K1-NIS mice, dose-dependent therapeutic effects were observed in both [211At]NaAt and [131I]NaI. In [211At]NaAt, a stronger tumour-growth suppression was observed, while tumour regrowth was not observed until 18, 25, and 46 days after injection of 0.4, 0.8, and 1.2 MBq of [211At]NaAt, respectively. While in [131I]NaI, this was observed within 12 days after injection (1, 3, and 8 MBq). The superior therapeutic effect of [211At]NaAt suggests the promising clinical applicability of targeted alpha therapy using [211At]NaAt in patients with differentiated thyroid cancer refractory to standard [131I]NaI treatment.

Development of the DICOM-based Monte Carlo dose reconstruction system for a retrospective study on the secondary cancer risk in carbon ion radiotherapy.

Objective.A retrospective study on secondary cancer risk on carbon ion radiotherapy (CIRT) is ongoing at the Heavy Ion Medical Accelerator in Chiba (HIMAC). The reconstruction of the whole-body patient dose distribution is the key issue in the study because dose distribution only around the planning target volume was evaluated in the treatment planning system.Approach.We therefore developed a new dose reconstruction system based on the Particle and Heavy Ion Transport code System (PHITS) coupled with the treatment plan DICOM data set by extending the functionalities of RadioTherapy package based on PHITS (RT-PHITS). In the system, the geometry of patient-specific beam devices such as the range shifter, range compensator, and collimators as well as the individual patient's body are automatically reconstructed. Various functions useful for retrospective analysis on the CIRT are implemented in the system, such as those for separately deducing dose contributions from different secondary particles and their origins.Main results.The accuracy of the developed system was validated by comparing the dose distribution to the experimental data measured in a water tank and using a treatment plan on an anthropomorphic phantom.Significance.The extended RT-PHITS will be used in epidemiological studies based on clinical data from HIMAC.

Longitudinal Change in Retinal Nerve Fiber Layer Thickness and Its Association With Central Retinal Sensitivity After Epiretinal Membrane Surgery.

PURPOSE: To investigate the longitudinal changes in the peripapillary retinal nerve fiber layer (pRNFL) thickness after epiretinal membrane (ERM) vitrectomy with internal limiting membrane (ILM) peeling, examine associations between pRNFL thickness and central retinal sensitivity, and identify predictors of postoperative pRNFL thickness. DESIGN: Prospective, observational, cohort study. METHODS: This study enrolled 82 eyes of 82 Japanese patients that underwent surgery for unilateral idiopathic ERM, with their fellow eyes as controls. pRNFL thickness was measured in 4 (superior, temporal, inferior, and nasal) quadrants preoperatively and at 1, 3, 6, and 12 months postoperatively. Microperimetry was performed at 12 months postoperatively to evaluate central retinal sensitivity. Regression tree analysis was performed to predict pRNFL thickness at 12 months postoperatively. RESULTS: The temporal quadrant showed continuous pRNFL thinning after surgery, reaching statistical significance at 3, 6, and 12 months postoperatively (all P < 0.001). The pRNFL thicknesses in the fellow eyes significantly increased at all postoperative time points (all P < 0.001). At 12 months postoperatively, the average central retinal sensitivity was significantly correlated with the temporal pRNFL thickness in the eyes with ERM (r = 0.372, P < 0.001); no significant correlation was found in the fellow eyes. Regression tree analysis showed that the preoperative pRNFL thickness in the temporal quadrant and patient age were the main determinants of the temporal pRNFL thickness at 12 months postoperatively. CONCLUSIONS: The risk of deterioration of central retinal sensitivity after ERM vitrectomy with internal limiting membrane peeling should be considered for patients with thin temporal pRNFLs and older adults.

Microdosimetric Modeling of Relative Biological Effectiveness for Skin Reactions: Possible Linkage Between In Vitro and In Vivo Data.

PURPOSE: Precise evaluation of the relative biological effectiveness (RBE) for skin reactions is indispensable for the treatment planning of particle therapy and boron-neutron capture therapy. The evaluation is also needed for the future radiation protection system for mixed radiation fields. Such RBE is often evaluated based on in vitro cell survival data, but its validity remains incompletely understood. This study aimed to develop a model for estimating RBE for skin reactions and dermal cell survival in the same framework and quantitatively discuss a possible link between them. METHODS AND MATERIALS: The microdosimetric kinetic model, which was originally developed for estimating cell surviving fractions for various radiations, was improved to be capable of estimating the mean and uncertainty of RBE for skin reactions. The parameter used in the model was independently determined from in vitro measurements of dermal cell survival and in vivo measurements of skin reactions taken from 8 and 23 articles, respectively. In the parameter determination, the characteristics of the radiation fields employed in each measurement were reproduced in detail by the Particle and Heavy Ion Transport Code System. RESULTS: Our model quantitatively revealed that RBE for skin reactions tend to be higher than those for dermal cell survival. RBE of various monoenergetic radiations calculated from this model confirmed that the past evaluations made by the International Commission on Radiological Protection and the National Council on Radiation Protection and Measurements a few decades ago are still supported by recent experimental data. CONCLUSIONS: Our model can play important roles not only in medical physics for avoiding unnecessary skin reactions in particle therapy and boron-neutron capture therapy but also in radiation protection for future decision-making of the recommended RBE values.

Inflammatory Signaling and DNA Damage Responses after Local Exposure to an Insoluble Radioactive Microparticle.

Cesium-bearing microparticles (Cs-BMPs) can reach the human respiratory system after inhalation, resulting in chronic local internal exposure. We previously investigated the spatial distribution of DNA damage induced in areas around a Cs-BMP; however, the biological impacts have not been fully clarified due to the limited amount of data. Here, we investigated the inflammatory signaling and DNA damage responses after local exposure to a Cs-BMP in vitro. We used two normal human lung cell lines, i.e., lung fibroblast cells (WI-38) and bronchial epithelial cells (HBEC3-KT). After 24 h exposure to a Cs-BMP, inflammation was evaluated by immunofluorescent staining for nuclear factor κB (NF-κB) p65 and cyclooxygenase 2 (COX-2). The number of DNA double-strand breaks (DSBs) was also detected by means of phospholylated histone H2AX (γ-H2AX) focus formation assay. Cs-BMP exposure significantly increased NF-κB p65 and COX-2 expressions, which were related to the number of γ-H2AX foci in the cell nuclei. Compared to the uniform (external) exposure to 137Cs γ-rays, NF-κB tended to be more activated in the cells proximal to the Cs-BMP, while both NF-κB p65 and COX-2 were significantly activated in the distal cells. Experiments with chemical inhibitors for NF-κB p65 and COX-2 suggested the involvement of such inflammatory responses both in the reduced radiosensitivity of the cells proximal to Cs-BMP and the enhanced radiosensitivity of the cells distal from Cs-BMP. The data show that local exposure to Cs-BMP leads to biological effects modified by the NF-κB pathway, suggesting that the radiation risk for Cs-BMP exposure can differ from that estimated based on conventional uniform exposure to normal tissues.

Transcatheter coil embolization for large pulmonary arteriovenous fistulae through an artificial tricuspid ball valve.

Pulmonary arteriovenous fistulae (PAFs) occur congenitally or are acquired. A PAF can cause hypoxemia, sudden death from rupture, abscess formation, and embolism. Treatment for PAF is transcatheter embolization or surgery. Transcatheter embolization is the first choice of treatment; however, this treatment is impossible to perform if a patient has had tricuspid or pulmonary valve replacement. In this paper, we describe a case of PAFs complicated with tricuspid valve replacement with a ball valve (which had been performed 40 years earlier) that was treated with transcatheter embolization. .

Track-structure modes in particle and heavy ion transport code system (PHITS): application to radiobiological research.

PURPOSE: In radiation physics, Monte Carlo radiation transport simulations are powerful tools to evaluate the cellular responses after irradiation. When investigating such radiation-induced biological effects, it is essential to perform track structure simulations by explicitly considering each atomic interaction in liquid water at the sub-cellular and DNA scales. The Particle and Heavy-Ion Transport code System (PHITS) is a Monte Carlo code which enables to calculate track structure at DNA scale by employing the track-structure modes for electrons, protons and carbon ions. In this paper, we review the recent development status and future prospects of the track-structure modes in the PHITS code. CONCLUSIONS: To date, the physical features of these modes have been verified using the available experimental data and Monte Carlo simulation results reported in literature. These track-structure modes can be used for calculating microdosimetric distributions to estimate cell survival and for estimating initial DNA damage yields. The use of PHITS track-structure mode is expected not only to clarify the underlying mechanisms of radiation effects but also to predict curative effects in radiation therapy. The results of PHITS simulations coupled with biophysical models will contribute to the radiobiological studies by precisely predicting radiation-induced biological effects based on the Monte Carlo approach.

Successful repeated uterine artery embolization in postpartum hemorrhage with disseminated intravascular coagulation: a case report and literature review.

BACKGROUND: Postpartum hemorrhage (PPH) is a potentially fatal condition requiring urgent and appropriate intervention. Uterine artery embolization (UAE) has a high hemostatic capacity for PPH, but it may fail. Disseminated intravascular coagulation (DIC) has been reported as a risk factor associated with the failure of UAE. CASE PRESENTATION: A 37-year-old primigravida with dichorionic diamniotic twins and placenta previa underwent cesarean section. The blood loss during surgery was 4950 mL. Hemostasis was achieved using an intrauterine balloon tamponade device. However, she lost a further 2400 mL of blood 5 h after surgery. We embolized both uterine arteries using gelatin sponges and confirmed hemostasis. She was suffering from DIC and received ample blood transfusions. However, a further 1300 mL of blood was lost 18 h after surgery and we performed repeated UAE, with complete recanalization of the uterine arteries on both sides and re-embolization with gelatin sponges. Her DIC was treated successfully by blood transfusions at this time, and she showed no further bleeding after the repeated UAE. CONCLUSIONS: DIC is a risk factor for the failure of UAE. Repeated UAE may be effective after sufficient improvement of the hematological status in patients with PPH and DIC.

Probabilistic risk assessment of solar particle events considering the cost of countermeasures to reduce the aviation radiation dose.

Cosmic-ray exposure to flight crews and passengers, which is called aviation radiation exposure, is an important topic in radiological protection, particularly for solar energetic particles (SEP). We therefore assessed the risks associated with the countermeasure costs to reduce SEP doses and dose rates for eight flight routes during five ground level enhancements (GLE). A four-dimensional dose-rate database developed by the Warning System for Aviation Exposure to Solar Energetic Particles, WASAVIES, was employed in the SEP dose evaluation. As for the cost estimation, we considered two countermeasures; one is the cancellation of the flight, and the other is the reduction of flight altitudes. Then, we estimated the annual occurrence frequency of significant GLE events that would bring the maximum flight route dose and dose rate over 1.0 mSv and 80 µSv/h, respectively, based on past records of GLE as well as historically large events observed by the cosmogenic nuclide concentrations in tree rings and ice cores. Our calculations suggest that GLE events of a magnitude sufficient to exceed the above dose and dose rate thresholds, requiring a change in flight conditions, occur once every 47 and 17 years, respectively, and their conservatively-estimated annual risks associated with the countermeasure costs are up to around 1.5 thousand USD in the cases of daily-operated long-distance flights.

Implementation of simplified stochastic microdosimetric kinetic models into PHITS for application to radiation treatment planning.

PURPOSE: The stochastic microdosimetric kinetic (SMK) model is one of the most sophisticated and precise models used in the estimation of the relative biological effectiveness of carbon-ion radiotherapy (CRT) and boron neutron capture therapy (BNCT). However, because of its complicated and time-consuming calculation procedures, it is nearly impractical to directly incorporate this model into a radiation treatment-planning system. MATERIALS AND METHODS: Through the introduction of Taylor expansion (TE) or fast Fourier transform (FFT), we developed two simplified SMK models and implemented them into the Particle and Heavy Ion Transport code System (PHITS). To verify the implementation, we calculated the photon isoeffective doses in a cylindrical phantom placed in the radiation fields of passive CRT and accelerator-based BNCT. RESULTS AND DISCUSSION: Our calculation suggested that both TE-based and FFT-based SMK models can reproduce the data obtained from the original SMK model very well for absorbed doses approximately below 5 Gy, whereas the TE-based SMK model overestimates the original data at higher doses. In terms of computational efficiency, the TE-based SMK model is much faster than the FFT-based SMK model. CONCLUSION: This study enables the instantaneous calculation of the photo isoeffective dose for CRT and BNCT, considering their cellular-scale dose heterogeneities. Treatment-planning systems that use the improved PHITS as a dose-calculation engine are under development.

Medical application of particle and heavy ion transport code system PHITS.

The Particle and Heavy Ion Transport code System (PHITS) is a general-purpose Monte Carlo simulation code that has been applied in various areas of medical physics. These include application in different types of radiotherapy, shielding calculations, application to radiation biology, and research and development of medical tools. In this article, the useful features of PHITS are explained by referring to actual examples of various medical applications.

Technical Note: validation of a material assignment method for a retrospective study of carbon-ion radiotherapy using Monte Carlo simulation.

We propose a two-step method to converse human tissue materials from patient computed tomography (CT) images, which is required in dose reconstructions for a retrospective study of carbon-ion radiotherapy (CIRT) using Monte Carlo (MC) simulation. The first step was to assign the standard tissues of the International Commission on Radiological Protection reference phantoms according to the CT-number. The second step was to determine the mass density of each material based on the relationship between CT-number and stopping power ratio (Hounsfield unit [HU]-SPR) registered in treatment planning system (TPS). Direct implementation of the well-calibrated HU-SPR curve allows the reproduction of previous clinical treatments recorded in TPS without uncertainty due to a mismatch of the CT scanner or scanning conditions, whereas MC simulation with realistic human tissue materials can fulfill the out-of-field dose, which was missing in the record. To validate our proposed method, depth-dose distributions in the homogenous and heterogeneous phantoms irradiated by a 400 MeV/u carbon beam with an 8 cm spread-out Bragg peak (SOBP) were computed by the MC simulation in combination with the proposed methods and compared with those of TPS. Good agreement of the depth-dose distributions between the TPS and MC simulation (within a 1% discrepancy in range) was obtained for different materials. In contrast, fluence distributions of secondary particles revealed the necessity of MC simulation using realistic human tissue. The proposed material assignment method will be used for a retrospective study using previous clinical data of CIRT at the National Institute of Radiological Sciences (NIRS).