APP下载

Monte Carlo study of the neutron ambient dose equivalent at the heavy ion medical machine in Wuwei

2022-11-21ShengCongHuangHuiZhangKaiBaiYingLuoHaiJunMaoZhongYingDai

Nuclear Science and Techniques 2022年9期

Sheng-Cong Huang• Hui Zhang • Kai Bai • Ying Luo • Hai-Jun Mao •Zhong-Ying Dai

Abstract The use of carbon-ion radiotherapy (CIRT) is gradually increasing. Owing to the generation of high-energy secondary neutrons during CIRT,its use presents new challenges in radiation protection.Thus,secondary neutron dose distributions must be explored and evaluated under clinical scenarios based on different treatment configurations. However, neutron dose and energy spectrum measurements are often difficult. This can be primarily attributed to the inherent limitations of most neutron detectors, such as their unsuitability for spectral measurements and inaccurate responses to neutrons with energies above 20 MeV. Numerical calculation methods based on probabilistic statistical theory are fast and convenient for neutron dose evaluation. In this study, external secondary neutron doses at the heavy ion medical machine in Wuwei,which is equipped with a passive beam delivery system,were calculated using the Monte Carlo method. The dependence of neutron doses on various treatment parameters (incident carbon-ion beam energy, spatial location,field size, and spread-out Bragg peak (SOBP) width) was investigated. Furthermore, the feasibility of applying an analytical model to predict the ambient dose equivalent was verified. For the combination involving an energy of 400 MeV/u and SOBP width of 6 cm, the ambient dose equivalent per therapeutic dose(H/D)at the isocenter was 79.87 mSv/Gy. The H/D value decreased rapidly with increasing spatial distance and slightly with increasing aperture size and SOBP width. The H/D values derived from the Monte Carlo simulations were in good agreement with the results reported in the literature. The analytical model could be used to quickly predict the H/D value along the incidence direction of the beam with an error of less than 20%. Thus, our study contributes to the understanding of the relationship between neutron radiation and treatment configuration parameters, which establishes a basis for predicting non-therapeutic radiation doses in CIRT.

Keywords Carbon-ion radiotherapy ∙Passive beam delivery ∙Secondary neutrons ∙Monte Carlo method

1 Introduction

Heavy-ion therapy plays an important role in the treatment of radiation-resistant tumors. This procedure heavily relies on the inverted depth dose distribution method and the superior radiobiological effect characteristics of heavyion beam radiation in reducing the dose exposure to normal tissues. Among the various types of heavy-ion beams,carbon-ion beams are currently the most suitable for clinical treatment [1]. The Institute of Modern Physics (IMP)used the Lanzhou heavy ion research facility (HIRFL) to conduct preclinical trials on heavy-ion cancer treatment in China. Subsequently, IMP, in cooperation with Lanzhou Kejintaiji Co., Ltd., independently developed the first medical heavy-ion accelerator demonstration device in China (product model: HIMM-GS1-WW1). The device is in service in Wuwei city of Gansu province,whereas it has been installed in Lanzhou and is in the commissioning,monitoring,and registration phases.In addition,carbon-ion particle therapy systems have been set up in Fujian, Zhejiang,and other regions[2].With the increased deployment of carbon-ion radiotherapy (CIRT), research on the associated radiation problems is required.

The carbon-ion beam profile is modulated by the beam delivery system before irradiation owing to the fact that the profile appears diminished after it is withdrawn from the accelerator. This passive beam delivery method uses configured hardware components to make the beam more conformal and thus demonstrates improved reliability,lower time consumption, insensitivity to target movement,and strong anti-interference ability[3].However,the use of hardware components has significantly increased the secondary radiation dose outside the primary field because secondary particles,such as neutrons,are generated when a beam interacts with a phantom or patient. Several studies have indicated that secondary neutrons in radiotherapy have a broad energy spectrum: for high-energy photon(>10 MeV) radiotherapy, the maximum energy produced by a device is 19–20 MeV, with a greater concentration at the 1 MeV level. Only a small fraction of neutrons have energy levels greater than 10 MeV, and for particle radiotherapy, the highest energy of secondary neutrons is up to twice the energy of the incident beam,which is up to several hundred mega-electron-volt [4, 5].Considering the relatively greater biological effects of neutrons, low-dose neutrons continue to have a greater impact on the probability of normal tissue complications and the risk of secondary tumors [6, 7].

Experimental studies have indicated that the ambient dose equivalent per therapeutic dose (H/D) ranges from less than 0.1 to as high as 80 mSv/Gy [8–10]. The instruments usually used for measuring the ambient neutron dose equivalent during particle therapy are wide-energy neutron detection instruments (WENDI-II) or extended-range Bonner spheres (ERBS) [5, 11–14]. The former only provides dose equivalents without the neutron spectrum, and it is impossible to obtain detailed information regarding secondary neutrons. In the latter measurement, the moderator must be replaced according to the differences in the neutron energy region. Because analyzing the secondary neutron spectrum is tedious and unsuitable for large-scale measurements at multiple locations,researchers have begun to consider using the Monte Carlo method to obtain simulated data[15–18].Based on this,we attempted to establish an empirical formula for the secondary neutron dose equivalent [16, 19]. Furthermore,most studies on secondary neutrons produced in radiotherapeutic processes have focused on photon and proton radiotherapy, and only limited studies have focused on secondary neutrons produced in the process of CIRT.

In this study, we used the Monte Carlo method to estimate neutron ambient dose equivalents in a patient’s environment for passive CIRT using a heavy-ion medical machine.The variation law of H/D values in the isocenter and patient environment areas under different treatment parameter configuration conditions, such as beam energy,multileaf collimator (MLC) aperture size, and spread-out Bragg peak (SOBP) width, was studied. In addition, the study verified the feasibility of applying an analytical model to predict the neutron dose equivalent to the heavy ion medical machine (HIMM) in Wuwei [16]. Thus, this study aimed to study the characteristics of external neutrons introduced by the passive beam delivery system,evaluate the proposed analytical model, and analyze the feasibility of neutron exposure in the HIMM.

2 Materials and methods

2.1 Monte Carlo simulations

The HIMM has four treatment terminals with active and passive beam delivery methods. The treatment terminals were used for irradiation from different directions: horizontal, vertical, horizontal + vertical, and 45° inclination.The dose rate at the isocenter was 0.001–1 Gy/s [2, 20].The maximum treatment depth in human tissue is 27 cm,which can be used for whole-body irradiation. The main hardware components in the nozzle include a primary collimator, ridge filter unit, range shifter unit, MLC,compensator, and two monitors (dose and beam position detectors). The beam was scanned by the scanning magnets, expanded into a 15×15 cm2irradiation field in the lateral direction, and then collimated by the primary collimator after passing through the vacuum window. When passing through the ridge filter, the pristine Bragg peak spread in the longitudinal direction, forming an SOBP.After the beam passed through the range shifter to modulate the range, the MLC was used for the lateral conformation. The final three-dimensional distribution of the high-dose zone was consistent with the shape of the target region. The use of the compensator made the dose distribution on the trailing edge of the target area steeper,which helped to protect endangered organs. The shape and position of the beam were measured using a beam position detector,and the intensity of the beam was monitored using a dose detector.

The beamline setting of the HIMM Room 2 in the Monte Carlo simulation is shown in Fig. 1. Prior research has shown that H/D values are strongly influenced by the presence of the patient or phantom, phantom size,anatomical location of the treatment site,and orientation of the patient or phantom[16,21].To eliminate the effects of internal neutrons and focus on machine-only behaviors,the free-in-the-air (without phantom) option was used in most of the simulations. The range shifters and compensators were simplified owing to their patient specificity. In addition, the monitors were not included in the simulation because their effects on secondary neutrons are negligible owing to their thinness compared to other components[17].The materials and corresponding distances of all the hardware components to the isocenter included in the calculational geometry are listed in Table 1.

The neutron spectral fluence was tallied into a set of 12-cm-diameter spherical receptors. A total of 16 receptors were located at the isocenter at distances of 25, 50, 100,150, and 200 cm from the isocenter and angles of 0°, 45°,and 90° for the incident beam direction (see Fig. 1)[15,16].The neutron spectral fluence φ(E)was tallied into 480 logarithmically spaced bins from 1 to 1GeV. In most simulation settings,a closed MLC was performed to obtain a conservative H/D value. The SOBP width was set to 6 cm, which is the standard setting for particle therapy[22].

The secondary neutron fluence spectrum for the HIMM passive beam delivery system was calculated using the GATE v8.2/Geant4-10-05.Geant4(geometry and tracking)is a set of Monte Carlo toolkits developed by CERN,using object-oriented technology based on C++,and its purpose is to simulate the physical process of particle transportation in matter [23]. Owing to its flexibility, Geant4 is also widely used in high-energy physics,nuclear physics,space and astrophysics, and radiation protection and detection,and it can be applied to the field of medical physics with the addition of expansion packages. GATE was designed by the OpenGATE team, based on the Geant4 kernel.This simulation platform is mainly used for Monte Carlo simulations of computed tomography and can also simulate particle radiotherapy.

The electromagnetic and hadronic physical processes were implemented using the physics list QGSP_BERT_HP_EMV, according to the recommendation from the OpenGATE collaboration. The cut (production threshold) for particle transport was 1 mm. 1×107carbon-ion histories were tracked for each simulation.Statistical errors were less than 5% for most receptors. All calculations were conducted on a CPU server (Intel®Xeon® CPU E5-2690 V3 @2.60 GHz×48, OS: Scientific Linux release 7.6) and parallel [48] was used to improve computing efficiency.

2.2 Calculation of ambient dose equivalent and H/D values

To assess neutron exposure at various locations, the focus was first on the ambient dose equivalent obtained using the fluence spectrum and neutron fluence-to-ambient dose-equivalent conversion coefficients.

where i is the ith neutron energy bin in the neutron spectral fluence, φithe fluence in the ith energy bin, and hiis the conversion coefficient of the ith neutron energy. Figure 2 shows the response function of WENDI-II, a typical wideenergy neutron detection instrument, together with the ICRP74 fluence-to-ambient dose conversion coefficients and NCRP38 fluence-to-dose equivalent conversion coefficient [24–27]. In this study, we chose the coefficients recommended by the ICRP74 to convert the dose.However, its maximum energy is 201 MeV; hence, the calculated results for the neutron region between 201 and 1000 MeV are shown in Fig. 2. These were generated using the Monte Carlo code Fluka by Ferrari et al.

Fig. 1 Model of the passive beam delivery system

Table 1 Components included in the simulation with their material, density, and distance to the isocenter

Fig. 2 (Color online) Response function of WENDI-II with fluence-to-ambient dose conversion coefficients by ICRP74, Ferrari et al., and the NCRP38 fluence-to-dose equivalent conversion coefficient

The H/D value is defined as the ambient dose equivalent per therapeutic dose and is calculated using Eq. (2).It is an important quantitative index for evaluating secondary radiation levels around the patient.

where p is the number of carbon ions, and D is the therapeutic dose. The biologically effective dose was obtained by relative biological effectiveness (RBE), weighting the physically absorbed dose on the beam central axis at the depth corresponding to the center of the SOBP. Table 2 lists the RBE values corresponding to the depth at the center of the SOBP under different beam energies. The absorbed dose was determined from separate simulations in which the MLC aperture size was fixed at 10 × 10 cm2and a water phantom was introduced. The water phantom(size: 30 × 30 × 35 cm3) was located upstream at the isocenter, normal to the incident beam.

Table 2 RBE values corresponding to the depth at the center of the SOBP under different beam energies

2.3 H/D value analytical model from Newhauser and co-workers

The spatial distribution law of H/D values in the treatment room has been reported in previous studies [15],as follows:

where d is the distance between the position of the detector and the effective neutron point source, and α is the power exponent in which the H/D value decreases as the distance increases. The physical meaning of an effective neutron point source is a point in the MLC where most of the carbon ions have been blocked.

The effective neutron point source distance d is defined as:

where disois the distance from the effective neutron point source to the isocenter of the treatment room, and r and θ are as shown in Fig. 1.

Inspired by this conclusion, Zheng et al. proposed an ambient dose-equivalent analytical model for a treatment room [16]. Using Eq. (5), the H/D values at different spatial positions at distance d from the effective neutron point source can be estimated.

where arefand mrefare the aperture sizes of the MLC and SOBP, respectively. The reference values of arefand mrefwere 0 cm2and 6 cm,respectively.The coefficients Caand Cmwere introduced to indicate the effects of the field size and SOBP width on the H/D value.

3 Results and discussion

3.1 Characteristics of secondary neutron fluence spectrum and ambient dose equivalent spectrum

The simulated neutron fluence spectrum and its corresponding ambient dose-equivalent spectrum for beam energies of 165,200,250,300,350,and 400 MeV/u at the isocenter are shown in Fig. 3. The simulation results are shown in a linear scale versus the logarithm of neutron energy, which ensures that the relative contribution from each bin is proportional to its area. As expected, the neutron fluence increased with the beam energy. The fluence spectrum contains two pronounced peaks: one low-energy evaporation peak with a position of approximately 1 MeV and one high-energy direct neutron peak that ranges from 15 MeV to twice the energy of the carbon ions. The distribution law of the dose-equivalent spectrum was analogous to that of the fluence spectrum. However, the contribution of low-energy neutrons decreases, and that of high-energy neutrons increases. This is explained by the fact that the ambient dose-equivalent conversion coefficient of neutrons with energy less than 100 keV is very low. Up to 23% of the ambient dose equivalent was from neutrons with energy ≥15 MeV.

These results are consistent with previously published results [5,28],except for the absence of a thermal neutron peak at approximately 10-7MeV. This difference might be related to the simulated parameter settings;the peak can be attributed to scattered neutrons (due to room return),whereas in this study, the shielding wall was simplified.Owing to the low fluence and small corresponding conversion factors of this peak,the effect of its absence on the ambient-dose equivalent calculations is negligible.

3.2 Effect of spatial distance on the H/D value

The H/D values for different energies for the 6-cm SOBP width,and closed MLC options were calculated and fitted according to the power-law formula. It can be found that the distribution of the H/D value at different spatial locations conforms to the power-law distribution, and the H/D value decreases rapidly with the increase in distance(Fig. 4a).Some data points deviated from the fitting curve;therefore, the data in different directions were normalized and fitted separately. The fitting parameters and goodnessof-fit values are listed in Table 3.

As shown in Fig. 4, the H/D values in all directions fit well with the power-law distribution. However, the H/D values of the locations in the 0° direction exhibited the smallest average error using the fitting formula (Fig. 4b).The H/D values in the 45°and 90°directions were slightly less consistent.This can be explained as follows:at 45°and 90°, the distance between the receptor and the isocenter increased; however, the increase in the distance between the receptor and MLC, which is the main producer of neutrons, was not synchronized with the former. This was particularly noticeable in the 90° direction. Therefore,when the subsequent analysis model predicted the H/D value at different spatial positions,the corresponding fitting parameters were used in different directions to improve accuracy.

3.3 Effect of SOBP width on H/D value

The 2-, 4-, 6-, 8-, 10-, and 12-cm SOBP widths and standardized simulation data based on the 6 cm result were selected.Figure 5 shows the variations in the relative H/D values for different SOBP widths. As the width of the SOBP gradually increased,the relative H/D value intensity at each location decreased almost linearly, and the drop range was not significantly different (approximately 30%).This phenomenon demonstrates that the ambient dose equivalent around the isocenter decreases with the thickening of the ridge filter when the MLC is fully closed.Although this finding is unexpected, it is understandable.The multileaf collimator is the main source of secondary neutrons downstream of the beamline. An increase in the thickness of the ridge filter reduced the energy of the beam,thereby reducing the neutrons produced by the interaction of the beam with the MLC.

Fig.3 (Color online)a Neutron fluence spectrum and b neutron ambient dose equivalent spectrum at the isocenter for different carbon ion beam energy. (SOBP width is 6 cm,MLC is closed)

Fig.4 (Color online)Monte Carlo simulation H/D values and fitting curves a at all locations,b in the 0°direction,c the 45°direction,d the 90°direction, and e the 90° direction with a linear scale (SOBP width: 6 cm, MLC: closed)

Table 3 Fitting parameters and goodness-of-fit in different directions

3.4 Effect of aperture size on H/D value

Based on simulations using aperture sizes α from 0 × 0 cm2(closed)to 15 × 15 cm2in increments of 3 cm per side, the effect of the open size of the aperture on the relative H/D value is shown in Fig. 6.

The secondary neutron dose equivalent is inversely proportional to the opening size of the MLC. The relative H/D values decreased with increasing open aperture, and the minimum was below 5%. Neutrons generated by the beam-limiting device are the main source of secondary neutrons in radiotherapy.The trend of the simulated data in Fig. 6 agrees with the expectations, which also indicates that the results obtained using a closed MLC are conservative. When the opening size was expanded from fully closed to 15 × 15 cm2, the neutron dose equivalent was reduced by approximately 90% in the 45° and 90° directions, respectively. At r = 200 cm and θ = 0°, the relative H/D value drop was only approximately 60%, which was significantly smaller than that at other locations. Based on the observations in Sect. 3.1,secondary neutrons consist of two parts: high-energy direct neutrons and low-energy evaporation neutrons. The adopted cutoff energy of 15 MeV subdivided the secondary neutron dose equivalent into a low-energy neutron dose equivalent,H1,and a highenergy neutron dose equivalent, H2.

As shown in Fig. 7, the decreasing laws of the H/D values at different locations are distinct and are related to the neutron energy.At remote locations in the 0°direction,such as at 150 and 200 cm,the relative H1/D intensity first increased and then decreased. When the aperture size was at its maximum, the relative H1/D intensity was approximately 1 or greater than 1 (Fig. 7a, b). This phenomenon can be understood from the initial generation and eventual fate of the neutrons. As the aperture size increased, more carbon ions passed through the nozzle without creating neutrons in the MLC.In addition,a few secondary neutrons escaped from the upstream beamline. This difference can be attributed to two factors: a reduction in neutron production and an increase in neutron leakage. Because the leaking neutrons (mainly high-energy cascade neutrons)are forward-peaked,the high-energy neutron component at positions in the 0°direction increases.For locations further from the MLC, the change in relative value was more significant.

Because there is no interference from escaped high-energy neutrons in the 45° direction, the falling laws of the H1/D and H2/D at different locations are similar(Fig. 7c,d). In the 90° direction, both the relative intensities of H1/D and H2/D significantly decreased with an increase in the aperture size. However, the falling laws of H1/D and H2/D at the isocenter and 25 and 50 cm from the isocenterwere similar to those in the 45°direction.At 100,150,and 200 cm from the isocenter, the H1/D and H2/D gradually decreased with increasing aperture size (Fig. 7e, f). This phenomenon is considered to be related to the relative positions of the receptors and MLC. The locations of the isocenter and 25 and 50 cm from the isocenter in the 90°direction make a small angle with the MLC.Therefore,the relative intensities of the H1/D and H2/D falling laws are the same as those of the small-angle area. The locations 100, 150, and 200 cm from the isocenter relative to the MLC were large angles (above 45°), and the falling laws followed a new pattern. Most neutrons arriving at these locations are side-scattered from the MLC.The thicknesses of the shields for the side-scattering neutrons were different.The closer the secondary neutrons are generated to the side surface of the MLC, the thinner the shield is required to penetrate and the greater the dose contribution to the receptor.This also explains why the slope of the relative H/D intensity curve increases with aperture size.

Fig. 5 (Color online) H/D values as a function of SOBP width at isocenter θ = 0°, 45°,or 90° and r = 100 or 200 cm.The black line represented the binomial fit curve of the average H/D values for all simulation locations (the beam energy is 400 MeV/u, MLC is closed)

Fig. 6 (Color online) Ratio of neutron dose equivalent per therapeutic dose for an arbitrary aperture size (f)to that for a closed aperture,a/0,as a function of aperture size: a H/D value in 0°, b H/D value in 45°, and c H/D value in 90° (beam energy of 400 MeV/u,SOBP width of 6 cm)

Coefficient Cawas proposed to indicate the influence of the MLC aperture size.Considering that the falling laws of H/D in different directions were not completely consistent,we fitted the change in H/D values in different directions separately(Fig. 8).The fitting parameters and goodness-offit values are presented in Table 4.

Fig. 7 (Color online) Ratio of neutron dose equivalent per therapeutic dose for an arbitrary aperture size (f)to that for a closed aperture,() a/()0,as a function of aperture size:a H1/D value in 0°,b H2/D value in 0°,c H1/D value in 45°,d H2/D value in 45°,e H1/D value in 90°,and f H2/D value in 90°(beam energy of 400 MeV/u,SOBP width of 6 cm)

Fig.8 (Color online)a 0° (beam direction),b 45°,c and 90° (beam transverse direction) relative H/D value fitting curves.The linear formula was used to fit the change of the H/D value in each direction

Table 4 Fitting parameters and goodness-of-fit in different directions

3.5 Validation of the analytical model

SOBP width on the prediction accuracy was insignificant,whereas the error caused by the aperture size was relatively large. Figure 8c shows that the error of the analysis model was approximately 20% within the classical field size(10 × 10 cm2) in the 0° direction, which was within the typical neutron measurement uncertainty. The results illustrate the feasibility of the proposed analytical model.If the analytical model can be optimized or the predicted data corrected, the adaptability and accuracy of the evaluation model can be further improved.

3.6 Out-of-field neutron dose performance evaluation

The H/D values (perpendicular to the beam direction)from the field were compared with those reported in the literature(Fig. 9).The 300 MeV/u beam energy and 6 cm SOBP widths were selected to be comparable to those found in published articles. Table 5 lists the beamline parameters used in previous studies. The H/D value is related to the specific model and design of the equipment and is affected by many factors, such as the beam energy,aperture size, air gap, phantom, location, and method of measurement or calculation. Therefore, these differences must be considered when comparing H/D values of various facilities (Fig. 10).

Fig.9 (Color online)Ratio of the predicted values,pred,to the Monte Carlo simulated values,mc.a different SOBP width in 0°,45°,and 90° directions, b different MLC aperture sizes in 0°, 45°, and 90° directions, c different MLC aperture sizes in the 0° direction

Table 5 Parameters in the beamline used in different studies

Fig. 10 (Color online)Relationship between the neutron dose equivalent per therapeutic dose and the distance from the edge of the field using different radiotherapy equipment (the upper and lower bounds describe the empirical range of doses associated with different intensity modulated radiation therapy techniques [30])

The neutron ambient dose equivalents in carbon and proton radiotherapy are comparable in the horizontal direction of the out-of-field area.A popular explanation for this phenomenon is that the cross section for nuclear interactions is higher for carbon-ion beams than for proton beams, but fewer carbon ions are required to deliver the same therapeutic dose:approximately 1010carbon ions per treatment, which is a factor of 100 fewer ions than proton therapy [31, 32].

Compared to the experimental measurement data of the carbon-ion radiotherapy device of the Hyogo Ion Beam Medical Center (HIBMC) by Yonai et al., the simulation results in this study were 20 times greater [11]. This discrepancy may originate from the difference between the beamline components and the calculation method. (1) The collimator system in this study is composed of a fixed primary collimator and a 6.48-cm thick tungsten copper alloy MLC. In reference 10, an adjustable four-leaf collimator was used for collimation in advance to reduce the beam profile to the required field size,and a 140-mm thick iron MLC was used, which could significantly reduce the neutron dose. (2) Reference 10 used the WENDI-II wideenergy neutron dose equivalent monitor for measurement,whereas this study used the Monte Carlo method for calculation. The Monte Carlo method results in a calculated value that is larger than the measured value [33].

4 Conclusion

In this study, the Monte Carlo method was used to calculate the external neutron ambient dose equivalent,which was generated in the passive beam delivery system for the HIMM facility in Wuwei, and an analytical model for the H/D values under different treatment parameter settings was established and compared with published data.Our simulation results showed that: (1) for CIRT with a passive beam delivery system,the ambient dose equivalent from external neutrons cannot be ignored; (2) various therapeutic parameter settings, such as locations, collimator aperture size, and SOBP width, bring diverse neutron exposures during CIRT; (3) the analytical model for predicting the external H/D values in the literature is feasible for the HIMM facility;and(4)some strategies,such as the installation of an adjustable pre-collimator and an increase in the MLC thickness, are applicable to optimize a beam line design to reduce secondary neutron radiation.

At present, most neutron detectors have limited responses at energies greater than 20 MeV.However,these instruments cannot meet the monitoring requirements of high-energy secondary neutrons during particle therapy.The present study implies that the Monte Carlo method can be employed to obtain the complete energy spectrum of secondary neutrons at various locations,thereby improving the efficiency and accuracy of acquiring the secondary neutron spectrum. Additionally, the use of the neutron fluence-dose conversion coefficient in the ICRP 74 report and reference 21 effectively avoids the errors caused by the energy response of the detectors.

Although the results of the present study are based on the HIMM facility in Wuwei, owing to the similarity of passive beam delivery systems, the simulation method is also applicable to other machines equipped with passive systems. The configurations in this work (free beam condition, closed MLC) are conservative, mainly focusing on the external neutrons produced by the beamline components.

There is growing evidence that the hazard of non-therapeutic neutron radiation generated during radiotherapy is nonnegligible and worthy of further investigation. This work contributes to our understanding of non-therapeutic radiation during CIRT and provides a basis for further reduction and accurate assessment of the neutron radiation received by patients, especially for pediatric and young patients with long life expectancy.

Author contributionsAll authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Sheng-Cong Huang, Hui Zhang, Kai Bai, Ying Luo, Hai-Jun Mao and Zhong-Ying Dai. The first draft of the manuscript was written by Sheng-Cong Huang, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.