[en] Retained austenite(RA) usually presents in the quenched Nuclear Pressure-Vessel SA508 Gr.3 steel. In the present work, the characteristic of RA decomposition and its effect on the impact toughness were investigated by microstructure observation, dilatometric experiments and Charpy impact tests. The results show that the RA transformed into martensite and bainite during tempering at 230 °C and 400 °C respectively, while mixture of long rod carbides and ferrite formed at 650 °C. The long rod carbides formed from RA decomposition decrease the critical cleavage stress for initiation of micro-cracks, and deteriorate the impact toughness of the steel. Pre-tempering at a low temperature such as 230 °C or 400 °C leading to the decomposition of RA into martensite or baintie can eliminate the deterioration of the toughness caused by direct decomposition into long rod carbides. The absorbed energy indicate that pre-tempering at 400 °C can drive dramatically improvement in the toughness of the steel. - Highlights: • The products of RA decomposition were localization observed by SEM and TEM. • Decomposition characteristic of RA were revealed during tempering at different temperature. • Impact toughness was dramatically improved by pre-tempering treatment.

[en] The microstructural evolution of reactor pressure vessel (RPV) steel and its effect on the mechanical properties during tempering at 650 °C were studied to reveal the time-dependent toughness and temper embrittlement. The results show that the toughening of the material should be attributed to the decomposition of the martensite/austenite constituents and uniform distribution of carbides. When the tempering duration was 5 h, the strength of the investigated steel decreased to strike a balance with the material impact toughness that reached a plateau. As the tempering duration was further increased, the material strength was slightly reduced but the material impact toughness deteriorated drastically. This time-dependent temper embrittlement is different from traditional temper embrittlement, and it can be partly attributed to the softening of the matrix and the broadening of the ferrite laths. Moreover, the dimensions and distribution of the grain carbides are the most important factors of the impact toughness. - Highlights: • The fracture mechanism of reactor pressure vessel (RPV) steels under impact load was investigated. • The Charpy V-notch impact test and the hinge model were employed for the study. • Grain boundary carbides play a key role in the impact toughness and fracture toughness. • The dependence of the deterioration of impact toughness on tempering time was analyzed for the first time.

[en] Purpose: To propose a simple model to explain the origin of ghost markers in marker-based optical tracking systems (OTS) and to develop retrospective strategies to detect and eliminate ghost markers. Methods: In marker-based OTS, ghost markers are virtual markers created due to the cross-talk between the two camera sensors, which can lead to system execution failure or inaccuracy in patient tracking. As a result, the users have to limit the number of markers and avoid certain marker configurations to reduce the chances of ghost markers. In this work, the authors propose retrospective strategies to detect and eliminate ghost markers. The two camera sensors were treated as mathematical points in space. The authors identified the coplanar within limit (CWL) condition as the necessary condition for ghost marker occurrence. A simple ghost marker detection method was proposed based on the model. Ghost marker elimination was achieved through pattern matching: a ghost marker-free reference set was matched with the optical marker set observed by the OTS; unmatched optical markers were eliminated as either ghost markers or misplaced markers. The pattern matching problem was formulated as a constraint satisfaction problem (using pairwise distances as constraints) and solved with an iterative backtracking algorithm. Wildcard markers were introduced to address missing or misplaced markers. An experiment was designed to measure the sensor positions and the limit for the CWL condition. The ghost marker detection and elimination algorithms were verified with samples collected from a five-marker jig and a nine-marker anthropomorphic phantom, rotated with the treatment couch from −60° to +60°. The accuracy of the pattern matching algorithm was further validated with marker patterns from 40 patients who underwent stereotactic body radiotherapy (SBRT). For this purpose, a synthetic optical marker pattern was created for each patient by introducing ghost markers, marker position uncertainties, and marker displacement. Results: The sensor positions and the limit for the CWL condition were measured with excellent reproducibility (standard deviation ≤ 0.39 mm). The ghost marker detection algorithm had perfect detection accuracy for both the jig (1544 samples) and the anthropomorphic phantom (2045 samples). Pattern matching was successful for all samples from both phantoms as well as the 40 patient marker patterns. Conclusions: The authors proposed a simple model to explain the origin of ghost markers and identified the CWL condition as the necessary condition for ghost marker occurrence. The retrospective ghost marker detection and elimination algorithms guarantee complete ghost marker elimination while providing the users with maximum flexibility in selecting the number of markers and their configuration to meet their clinic needs

[en] Purpose: The purpose of this work is to investigate the impact of small rotational errors on the magnitudes and distributions of spatial dose variations for intracranial stereotactic radiotherapy (SRT) treatment setups, and to assess the feasibility of using the original dose map overlaid with rotated contours (ODMORC) method as a fast, online evaluation tool to estimate dose changes (using DVHs) to clinical target volumes (CTVs) and organs-at-risks (OARs) caused by small rotational setup errors. Methods: Fifteen intracranial SRT cases treated with either three-dimensional conformal radiation therapy (3DCRT) or intensity-modulated radiation therapy (IMRT) techniques were chosen for the study. Selected cases have a variety of anatomical dimensions and pathologies. Angles of ±3 deg. and ±5 deg. in all directions were selected to simulate the rotational errors. Dose variations in different regions of the brain, CTVs, and OARs were evaluated to illustrate the various spatial effects of dose differences before and after rotations. DVHs accounting for rotations that were recomputed by the treatment planning system (TPS) and those generated by the ODMORC method were compared. A framework of a fast algorithm for multicontour rotation implemented by ODMORC is introduced as well. Results: The average values of relative dose variations between original dose and recomputed dose accounting for rotations were greater than 4.0% and 10.0% in absolute mean and in standard deviation, respectively, at the skull and adjacent regions for all cases. They were less than 1.0% and 2.5% in absolute mean and in standard deviation, respectively, for dose points 3 mm away from the skull. The results indicated that spatial dose to any part of the brain organs or tumors separated from the skull or head surface would be relatively stable before and after rotations. Statistical data of CTVs and OARs indicate the lens and cochleas have the large dose variations before and after rotations, whereas the remaining ROIs have insignificant dose differences. DVH comparisons suggest that the ODMORC method is able to estimate the DVH of CTVs fairly accurately (within 1.5% of relative dose differences for evaluation volumes). The results also show that most of the OARs including the brain stem, spinal cord, chiasm, hippocampuses, optic nerves, and retinas, which were relatively distal from the skull and surface, had good agreement (within 2.0% of relative dose differences for 0.1 cc of the volumes ) between the ODMORC method and the recomputation, whereas OARs more proximate to the bone-tissue interface or surface, such as the lenses and cochlea, had larger dose variations (greater than 5.0%) for some cases due to the incapability of the ODMORC to account for scatter contribution variations proximate to interfaces and intrinsic dose calculation uncertainties for ROIs with small volumes. Conclusions: The ODMORC method can be implemented as an online evaluation system for rotation-induced dose changes of CTVs and most OARs and for other related dose consequence analyses.

[en] The ionization chamber volume averaging effect is a well-known issue without an elegant solution. The purpose of this study is to propose a novel convolution-based approach to address the volume averaging effect in model-based treatment planning systems (TPSs). Ionization chamber-measured beam profiles can be regarded as the convolution between the detector response function and the implicit real profiles. Existing approaches address the issue by trying to remove the volume averaging effect from the measurement. In contrast, our proposed method imports the measured profiles directly into the TPS and addresses the problem by reoptimizing pertinent parameters of the TPS beam model. In the iterative beam modeling process, the TPS-calculated beam profiles are convolved with the same detector response function. Beam model parameters responsible for the penumbra are optimized to drive the convolved profiles to match the measured profiles. Since the convolved and the measured profiles are subject to identical volume averaging effect, the calculated profiles match the real profiles when the optimization converges. The method was applied to reoptimize a CC13 beam model commissioned with profiles measured with a standard ionization chamber (Scanditronix Wellhofer, Bartlett, TN). The reoptimized beam model was validated by comparing the TPS-calculated profiles with diode-measured profiles. Its performance in intensity-modulated radiation therapy (IMRT) quality assurance (QA) for ten head-and-neck patients was compared with the CC13 beam model and a clinical beam model (manually optimized, clinically proven) using standard Gamma comparisons. The beam profiles calculated with the reoptimized beam model showed excellent agreement with diode measurement at all measured geometries. Performance of the reoptimized beam model was comparable with that of the clinical beam model in IMRT QA. The average passing rates using the reoptimized beam model increased substantially from 92.1% to 99.3% with 3%/3 mm and from 79.2% to 95.2% with 2%/2 mm when compared with the CC13 beam model. These results show the effectiveness of the proposed method. Less inter-user variability can be expected of the final beam model. It is also found that the method can be easily integrated into model-based TPS. (paper)

[en] Purpose: The use of sophisticated dose calculation procedure in modern radiation therapy treatment planning is inevitable in order to account for complex treatment fields created by multileaf collimators (MLCs). As a consequence, independent volumetric dose verification is time consuming, which affects the efficiency of clinical workflow. In this study, the authors present an efficient adaptive beamlet-based finite-size pencil beam (AB-FSPB) dose calculation algorithm that minimizes the computational procedure while preserving the accuracy. Methods: The computational time of finite-size pencil beam (FSPB) algorithm is proportional to the number of infinitesimal and identical beamlets that constitute an arbitrary field shape. In AB-FSPB, dose distribution from each beamlet is mathematically modeled such that the sizes of beamlets to represent an arbitrary field shape no longer need to be infinitesimal nor identical. As a result, it is possible to represent an arbitrary field shape with combinations of different sized and minimal number of beamlets. In addition, the authors included the model parameters to consider MLC for its rounded edge and transmission. Results: Root mean square error (RMSE) between treatment planning system and conventional FSPB on a 10 × 10 cm² square field using 10 × 10, 2.5 × 2.5, and 0.5 × 0.5 cm² beamlet sizes were 4.90%, 3.19%, and 2.87%, respectively, compared with RMSE of 1.10%, 1.11%, and 1.14% for AB-FSPB. This finding holds true for a larger square field size of 25 × 25 cm², where RMSE for 25 × 25, 2.5 × 2.5, and 0.5 × 0.5 cm² beamlet sizes were 5.41%, 4.76%, and 3.54% in FSPB, respectively, compared with RMSE of 0.86%, 0.83%, and 0.88% for AB-FSPB. It was found that AB-FSPB could successfully account for the MLC transmissions without major discrepancy. The algorithm was also graphical processing unit (GPU) compatible to maximize its computational speed. For an intensity modulated radiation therapy (∼12 segments) and a volumetric modulated arc therapy fields (∼90 control points) with a 3D grid size of 2.0 × 2.0 × 2.0 mm³, dose was computed within 3–5 and 10–15 s timeframe, respectively. Conclusions: The authors have developed an efficient adaptive beamlet-based pencil beam dose calculation algorithm. The fast computation nature along with GPU compatibility has shown better performance than conventional FSPB. This enables the implementation of AB-FSPB in the clinical environment for independent volumetric dose verification

[en] The microstructure evolution of a Nuclear Pressure Vessel Steel (SA508 Gr.3) with macrosegregation was investigated during quenching and tempering. Macrosegregation is quantitatively analyzed by EPMA, and the results show that concentration of Mn, Mo and Ni in positive segregation zone may be 0.3%-0.5% higher than in negative segregation zone. A mixture microstructure of martensite, proeutectoid ferrite, carbide-free baintie, upper bainite, M-A island has been obtained during quenching from austenitization at a rate of 45 ℃/min. Long strip shape carbides between carbide-free bainite ferrite trend to spheroidized during tempering, but they are still existing in the material. Suggestions are put forward for how to suppress the formation of carbide-free bainite and long strip shape carbides by heat treatmens. (authors)

[en] A sliding-window (SW) methodology for VMAT dose calculation was developed. For any two adjacent VMAT control points (CPs) n and n+1, the dose distribution was approximated by a 2-CP SW IMRT beam with the starting MLC positions at CP n and ending MLC positions at CP n+1, with the gantry angle fixed in the middle of the two VMAT CPs. Therefore, for any VMAT beam with N CPs, the dose is calculated with N-1 SW beams. VMAT plans were generated for ten patients in Pinnacle using 4° gantry spacing. For each patient, the VMAT plan was converted to a SW IMRT plan and dose was re-calculated. Another VMAT plan, with 1° gantry spacing, was created by interpolating the original VMAT beam. The original plans were delivered on an Elekta Versa HD and measured with Mapcheck2 using an in-house developed subarc method. For both the isodose distribution and DVH, there were significant differences between the original VMAT plan and either the SW or the interpolated plan. However, they were indistinguishable between the SW and interpolated plans. The average passing rate between the original VMAT plan and measurements was 84%. For both the interpolated and SW plans, the average passing rate was 96%. We conclude that the proposed SW approach improves VMAT dose calculation accuracy without increase in dose calculation time. (paper)

[en] The purpose of this study was to choose an appropriate head scatter source model for the fast and accurate independent planar dose calculation for intensity-modulated radiation therapy (IMRT) with MLC. The performance of three different head scatter source models regarding their ability to model head scatter and facilitate planar dose calculation was evaluated. A three-source model, a two-source model and a single-source model were compared in this study. In the planar dose calculation algorithm, in-air fluence distribution was derived from each of the head scatter source models while considering the combination of Jaw and MLC opening. Fluence perturbations due to tongue-and-groove effect, rounded leaf end and leaf transmission were taken into account explicitly. The dose distribution was calculated by convolving the in-air fluence distribution with an experimentally determined pencil-beam kernel. The results were compared with measurements using a diode array and passing rates with 2%/2 mm and 3%/3 mm criteria were reported. It was found that the two-source model achieved the best agreement on head scatter factor calculation. The three-source model and single-source model underestimated head scatter factors for certain symmetric rectangular fields and asymmetric fields, but similar good agreement could be achieved when monitor back scatter effect was incorporated explicitly. All the three source models resulted in comparable average passing rates (>97%) when the 3%/3 mm criterion was selected. The calculation with the single-source model and two-source model was slightly faster than the three-source model due to their simplicity

[en] The aim of this work is to investigate the clinical impact of detector size effect on patient specific intensity-modulated radiation therapy (IMRT) quality assurance (QA). Two photon beam models, BM6 and BM4, were commissioned using photon beam profiles measured with a 6 mm diameter and a 4 mm diameter ion chambers, respectively. A method was developed to extract the ''true'' cross beam profiles, free of volume averaging effect, using analytic fitting/deconvolution. The method was validated using beam profiles measured with a small (0.8 mm) diode detector for small (≤10x10 cm²) field sizes. These profiles were used to commission a third beam model (BM08). Planar dose distributions for eight IMRT plans (total of 53 fields) were calculated using the three beam models and measured with a two-dimensional detector array. Analysis using percent dose difference and distance-to-agreement criteria between the calculation and measurement was done to benchmark the performance of each beam model. The average passing rates between calculation and measurement were 93.8%, 98.9%, and 99.4% for BM6, BM4, and BM08, respectively, when 3%/3 mm criteria were used. A gradual increase in passing rates was noticed with the decrease in the size of the detector used to collect commissioning data. When 2%/2 mm criteria were used, the average passing rates increased significantly from 81.6% (BM6) to 92.6% (BM4) and 96.8% (BM08). These results quantify the enhancement of IMRT dose calculation accuracy with the reduction in detector size used for photon beam profiles measurement. Our study indicates that volume averaging effect can significantly affect the results of IMRT patient specific QA. By removing the detector size effect in beam commissioning, excellent passing rates can be achieved with more stringent criteria such as 2%/2 mm. The use of more stringent criteria for IMRT patient specific QA would likely result in higher chances of detecting any dosimetric errors arising from the treatment planning or delivery system