Monte Carlo dose verification for intensity-modulated arc therapy

Intensity-modulated arc therapy (IMAT), a technique which combines beam rotation and dynamic multileaf collimation, has been implemented in our clinic. Dosimetric errors can be created by the inability of the planning system to accurately account for the effects of tissue inhomogeneities and physical characteristics of the multileaf collimator (MLC). The objective of this study is to explore the use of Monte Carlo (MC) simulation for IMAT dose verification. The BEAM/DOSXYZ Monte Carlo system was implemented to perform dose verification for the IMAT treatment. The implementation includes the simulation of the linac head/MLC (Elekta SL20), the conversion of patient CT images and beam arrangement for 3D dose calculation, the calculation of gantry rotation and leaf motion by a series of static beams and the development of software to automate the entire MC process. The MC calculations were verified by measurements for conventional beam settings. The agreement was within 2%. The IMAT dose distributions generated by a commercial forward planning system (RenderPlan. Elekta) were compared with those calculated by the MC package. For the cases studied, discrepancies of over 10% were found between the MC and the RenderPlan dose calculations. These discrepancies were due in part to the inaccurate dose calculation of the RenderPlan system. The computation time for the IMAT MC calculation was in the range of 20-80 min on 15 Pentium-Ill computers. The MC method was also useful in verifying the beam apertures used in the IMAT treatments.     

Analysing collimator structure effects in head-scatter calculations for IMRT class fields using scatter raytracing

The frequent blocking of the irradiated volume in intensity modulated radiation therapy (IMRT) makes the head-scatter fraction of the incident photon fluence more significant than that in conventional therapy with open fields. On the other hand. certain collimator configurations block scatter photons directed to a given observation point while allowing primary photons to be transmitted. The 'anomalous blocking' makes the primary field a poor indicator of the scatter fluence. Since large MU-to-cGy ratios in IMRT can magnify head-scatter uncertainties, it becomes necessary to accurately model both the effective scatter source and the collimator structure that limits the scatter reaching the irradiated volume. First we obtain a dual-source model, using a Taylor series expansion to derive the effective scatter source distribution from the data measured for the Elekta SL20 linac equipped with a multi-leaf collimator (MLC). Then, using a raytracing algorithm, we calculate the transmission of scatter rays from the effective scatter source plane to points in the patient plane. The method can account for the anomalous blocking of scatter by the MLC leaves and the backup diaphragms. For a variety of collimator settings tested, the calculations agree with measurements to an accuracy of 0.002psi10 x 10, where psi10 x 10 is the total (primary + scatter) photon fluence of an open 10 x 10 cm2 field for the same MU delivered. Although the significance of collimator structure in IMRT depends strongly on fields shapes employed for the delivery, potential cumulative errors on the order of a few per cent can be avoided in fluence calculations if the proposed method is used.     

Using a photon phase-space source for convolution/superposition dose calculations in radiation therapy

For a given linac design, the dosimetric characteristics of a photon beam are determined uniquely by the energy and radial distributions of the electron beam striking the x-ray target. However, in the usual commissioning of a beam from measured data, a large number of variables can be independently tuned, making it difficult to derive a unique and self-consistent beam model. For example, the measured dosimetric penumbra in water may be attributed in various proportions to the lateral secondary electron range, the focal spot size and the transmission through the tips of a non-divergent collimator; the head-scatter component in the tails of the transverse profiles may not be easy to resolve from phantom scatter and head leakage; and the head-scatter tails corresponding to a certain extra-focal source model may not agree self-consistently with in-air output factors measured on the central axis. To reduce the number of adjustable variables in beam modelling, we replace the focal and extra-focal sources with a single phase-space plane scored just above the highest adjustable collimator in a EGS/BEAM simulation of the linac. The phase-space plane is then used as photon source in a stochastic convolution/superposition dose engine. A photon sampled from the uncollimated phase-space plane is first propagated through an arbitrary collimator arrangement and then interacted in the simulation phantom. Energy deposition kernel rays are then randomly issued from the interaction points and dose is deposited along these rays. The electrons in the phase-space file are used to account for electron contamination. 6 MV and 18 MV photon beams from an Elekta SL linac are used as representative examples. Except for small corrections for monitor backscatter and collimator forward scatter for large field sizes (<0.5% with <20 x 20 cm2 field size), we found that the use of a single phase-space photon source provides accurate and self-consistent results for both relative and absolute dose calculations.     

国内六大行政区域六城市中老年人群膝关节骨性关节炎患病危险因素比较

目的：了解中国不同地区间中老年人群膝关节骨性关节炎患病危险因素。方法：调查时间为2005—07／08。①从中国六大行政区（西北，华北，华东。中南，东北，西南）选出六城市（西安，石家庄，上海。广州，哈尔滨市，成都），用分层多阶段整群抽样方法，抽取6218名40岁及以上具有正式户口常住男女人群进行膝关节骨性关节炎的流行病学问卷调查（包括一般情况、现病史、既往史、体格检查、X射线片检查情况和疾病诊断6个方面，共计94个问题141个变量指标），并对其中4808名有症状者进行X射线平片膝正侧位投照。②膝关节骨性关节炎诊断标准为临床症状阳性加X射线Kellgren ＆ Lawrence分级二级及以上者。③计算患病率，并采用Epilnf06．0和SPSS 10．0软件对其中83个变量进行多因素非条件Logistfc回归分析，表示疾病与暴露因素之间联系强度的指标用比值比（OR），若OR〉1，说明疾病发生危险性增加，与暴露因素呈正关联；若OR〈1，说明疾病发生危险性减少，与暴露因素呈负关联。 结果：①六城市膝关节骨性关节炎总患病率为15．6％，其中西安7．7％，石家庄11．2％，上海9．8％。广州30．5％，哈尔滨16．9％，成都17．5％，各城市患病率比较差异显著（P〈0．01）。②Logistic回归分析膝关节骨性关节炎在大部分城市有共同的危险因素如年龄大（OR=1．032—1．181），使用蹲坑排便年限长（OR=1．021-1．077），体质量高（OR=1．048—1．073），和开始饮酒年龄大（OR=1．008～1．028）；而从事专职体育运动（OR=1．651，西安），骨质疏松病史（OR=3．311，石家庄），吸烟（OR=2．654，石家庄），类风湿关节炎病史（OR=4．964，上海），文化程度高（OR=2．593，上海），女性（OR=2．510，广州），姐妹骨关节炎史（OR=13．251，哈尔滨），母亲骨关节炎史（OR=5．683，成都）等危险因素分别在不同地区出现． 结论：年龄大、使用蹲坑排便年限长、体质量高和开始饮酒年龄大是中国六地区膝关节骨性关节炎患病的共同危险因素，同时，不同地区主要危险因素又有一定差异。