Three-dimensional dose verification for clinical treatments of small intracranial tumours

Cancers of the brain and central nervous system account for 1.6% of new cancers and 1.8% of cancer deaths globally. The highest rates of all developed nations are observed in Australia and New Zealand. There are known complexities associated with dose measurement of very small radiation fields. Here, 3D dosimetric verification of treatments for small intracranial tumours using gel dosimetry was investigated.

Monte Carlo dosimetry for stereotactic radiosurgery

Stereotactic radiosurgery (SRS) treatments for brain cancers require small and precisely shaped photon beams. These beams can be generated by fitting a linear accelerator with a micro-multileaf collimator (mMLC) such as the BrainLAB m3, which offers greater flexibility for field shaping than standard SRS cone collimators

CTCombine : rotating and combining CT data with an accurate detector model to simulate radiotherapy portal imaging at non-zero beam angles

Established Monte Carlo user codes BEAMnrc and DOSXYZnrc permit the accurate and straightforward simulation of radiotherapy experiments and treatments delivered from multiple beam angles. However, when an electronic portal imaging detector (EPID) is included in these simulations, treatment delivery from non-zero beam angles becomes problematic. This study introduces CTCombine, a purpose-built code for rotating selected CT data volumes, converting CT numbers to mass densities, combining the results with model EPIDs and writing output in a form which can easily be read and used by the dose calculation code DOSXYZnrc...

Solid-state MAS NMR measurements of intact articular bovine cartilage

Articular cartilage (AC), an avascular connective tissue lining articulating surfaces of the long bones, comprises extracellular biopolymers. In functionally compromised states such as osteoarthritis, thinned or lost AC causes reduced mobility and increased health-care costs. Understanding of the characteristics responsible for the load bearing efficiency of AC and the factors leading to its degradation are incomplete. DTI shows the structural alignment of collagen in AC [1] and T2 relaxation measurements suggest that the average director of reorientational motion of water molecules depends on the degree of alignment of collagen in AC [2]. Information on the nature of the chemical interactions involved in functional AC is lacking. The need for AC structural integrity makes solid state NMR an ideal tool to study this tissue. We examined the contribution of water in different functional ‘compartments’ using 1H-MAS, 13C-MAS and 13C-CPMAS NMR of bovine patellar cartilage incubated in D2O. 1H-MAS spectra signal intensity was reduced due to H/D exchange without a measureable redistribution of relative signal intensity. Chemical shift anisotropy was estimated by lineshape analysis of multiple peaks in the 1H-MAS spinning sidebands. These asymmetrical sidebands suggested the presence of multiple water species in AC. Therefore, water was added in small aliquots to D2O saturated AC and the influence of H2O and D2O on organic components was studied with 13C-MAS-NMR and 13C-CPMAS-NMR. Signal intensity in 13C-MAS spectra showed no change in relative signal intensity throughout the spectrum. In 13C-CPMAS spectra, displacement of water by D2O resulted in a loss of signal in the aliphatic region due to a reduction in proton availability for cross-polarization. These results complement dehydration studies of cartilage using osmotic manipulation [3] and demonstrate components of cartilage that are in contact with mobile water.

Sensitivity of an Elekta iView GT a-Si EPID model to delivery errors for pre-treatment verification of IMRT fields

A Monte Carlo model of an Elekta iViewGT amorphous silicon electronic portal imaging device (a-Si EPID) has been validated for pre-treatment verification of clinical IMRT treatment plans. The simulations involved the use of the BEAMnrc and DOSXYZnrc Monte Carlo codes to predict the response of the iViewGT a-Si EPID model. The predicted EPID images were compared to the measured images obtained from the experiment. The measured EPID images were obtained by delivering a photon beam from an Elekta Synergy linac to the Elekta iViewGT a-Si EPID. The a-Si EPID was used with no additional build-up material. Frame averaged EPID images were acquired and processed using in-house software. The agreement between the predicted and measured images was analyzed using the gamma analysis technique with acceptance criteria of 3% / 3 mm. The results show that the predicted EPID images for four clinical IMRT treatment plans have a good agreement with the measured EPID signal. Three prostate IMRT plans were found to have an average gamma pass rate of more than 95.0 % and a spinal IMRT plan has the average gamma pass rate of 94.3 %. During the period of performing this work a routine MLC calibration was performed and one of the IMRT treatments re-measured with the EPID. A change in the gamma pass rate for one field was observed. This was the motivation for a series of experiments to investigate the sensitivity of the method by introducing delivery errors, MLC position and dosimetric overshoot, into the simulated EPID images. The method was found to be sensitive to 1 mm leaf position errors and 10% overshoot errors.

Dosimetry of cone-defined stereotactic radiosurgery fields with a commercial synthetic diamond detector

Purpose Small field x-ray beam dosimetry is difficult due to a lack of lateral electronic equilibrium, source occlusion, high dose gradients and detector volume averaging. Currently there is no single definitive detector recommended for small field dosimetry. The objective of this work was to evaluate the performance of a new commercial synthetic diamond detector, namely the PTW 60019 microDiamond, for the dosimetry of small x-ray fields as used in stereotactic radiosurgery (SRS). Methods Small field sizes were defined by BrainLAB circular cones (4 – 30 mm diameter) on a Novalis Trilogy linear accelerator and using the 6 MV SRS x-ray beam mode for all measurements. Percentage depth doses were measured and compared to an IBA SFD and a PTW 60012 E diode. Cross profiles were measured and compared to an IBA SFD diode. Field factors, Ω_(Q_clin,Q_msr)^(f_clin,f_msr ), were calculated by Monte Carlo methods using BEAMnrc and correction factors, k_(Q_clin,Q_msr)^(f_clin,f_msr ), were derived for the PTW 60019 microDiamond detector. Results For the small fields of 4 to 30 mm diameter, there were dose differences in the PDDs of up to 1.5% when compared to an IBA SFD and PTW 60012 E diode detector. For the cross profile measurements the penumbra values varied, depending upon the orientation of the detector. The field factors, Ω_(Q_clin,Q_msr)^(f_clin,f_msr ), were calculated for these field diameters at a depth of 1.4 cm in water and they were within 2.7% of published values for a similar linear accelerator. The corrections factors, k_(Q_clin,Q_msr)^(f_clin,f_msr ), were derived for the PTW 60019 microDiamond detector. Conclusions We conclude that the new PTW 60019 microDiamond detector is generally suitable for relative dosimetry in small 6 MV SRS beams for a Novalis Trilogy linear equipped with circular cones.

Design and experimental testing of air slab caps which convert commercial electron diodes into dual purpose, correction-free diodes for small field dosimetry

Purpose Two diodes which do not require correction factors for small field relative output measurements are designed and validated using experimental methodology. This was achieved by adding an air layer above the active volume of the diode detectors, which canceled out the increase in response of the diodes in small fields relative to standard field sizes. Methods Due to the increased density of silicon and other components within a diode, additional electrons are created. In very small fields, a very small air gap acts as an effective filter of electrons with a high angle of incidence. The aim was to design a diode that balanced these perturbations to give a response similar to a water-only geometry. Three thicknesses of air were placed at the proximal end of a PTW 60017 electron diode (PTWe) using an adjustable “air cap”. A set of output ratios (ORfclin Det ) for square field sizes of side length down to 5 mm was measured using each air thickness and compared to ORfclin Det measured using an IBA stereotactic field diode (SFD). k fclin, f msr Qclin,Qmsr was transferred from the SFD to the PTWe diode and plotted as a function of air gap thickness for each field size. This enabled the optimal air gap thickness to be obtained by observing which thickness of air was required such that k fclin, f msr Qclin,Qmsr was equal to 1.00 at all field sizes. A similar procedure was used to find the optimal air thickness required to make a modified Sun Nuclear EDGE detector (EDGEe) which s “correction-free” in small field relative dosimetry. In addition, the feasibility of experimentally transferring k fclin, f msr Qclin,Qmsr values from the SFD to unknown diodes was tested by comparing the experimentally transferred k fclin, f msr Qclin,Qmsr values for unmodified PTWe and EDGEe diodes to Monte Carlo simulated values. Results 1.0 mm of air was required to make the PTWe diode correction-free. This modified diode (PTWeair) produced output factors equivalent to those in water at all field sizes (5–50 mm). The optimal air thickness required for the EDGEe diode was found to be 0.6 mm. The modified diode (EDGEeair) produced output factors equivalent to those in water, except at field sizes of 8 and 10 mm where it measured approximately 2% greater than the relative dose to water. The experimentally calculated k fclin, f msr Qclin,Qmsr for both the PTWe and the EDGEe diodes (without air) matched Monte Carlo simulated results, thus proving that it is feasible to transfer k fclin, f msr Qclin,Qmsr from one commercially available detector to another using experimental methods and the recommended experimental setup. Conclusions It is possible to create a diode which does not require corrections for small field output factor measurements. This has been performed and verified experimentally. The ability of a detector to be “correction-free” depends strongly on its design and composition. A nonwater-equivalent detector can only be “correction-free” if competing perturbations of the beam cancel out at all field sizes. This should not be confused with true water equivalency of a detector.

Force-deformation properties of the human heel pad during barefoot walking

Introduction: The plantar heel pad is a specialized fibroadipose tissue that attenuates and, in part, dissipates the impact energy associated with heel strike. Although near maximal deformation of the heel pad has been shown during running, in vivo measurement of the deformation and structural properties of the heel pad during walking remains largely unexplored. This study employed a fluoroscope, synchronized with a pressure platform, to obtain force–deformation data for the heel pad during walking. Methods: Dynamic lateral foot radiographs were acquired from 6 male and 10 female adults (age, 45 ± 10 yrs; height, 1.66 ± 0.10 m; and weight, 80.7 ± 10.8 kg), while walking barefoot at preferred speeds. The inferior aspect of the calcaneus was digitized and the sagittal thickness and deformation of the heel pad relative to the support surface calculated. Simultaneous measurement of the peak force beneath the heel was used to estimate the principal structural properties of the heel pad. Results: Transient loading profiles associated with walking induced rapidly changing deformation rates in the heel pad and resulted in irregular load–deformation curves. The initial stiffness (32 ± 11 N.mm-1) of the heel pad was an order of magnitude lower than its final stiffness (212 ± 125 N.mm-1) and on average, only 1.0 J of energy was dissipated by the heel pad with each step during walking. Peak deformation (10.3 mm) approached that predicted for the limit of pain tolerance (10.7 mm). Conclusion: These findings suggest the heel pad operates close to its pain threshold even at speeds encountered during barefoot walking and provides insight as to why barefoot runners may adopt ‘forefoot’ strike patterns that minimize heel loading.

Awards before and after the Nobel Prize : a Matthew effect and/or a ticket to one's own funeral?

The primary aim of this descriptive exploration of scientists’ life cycle award patterns is to evaluate whether awards breed further awards and identify researcher experiences after reception of the Nobel Prize. To achieve this goal, we collected data on the number of awards received each year for 50 years before and after Nobel Prize reception by all 1901–2000 Nobel laureates in physics, chemistry, and medicine or physiology. Our results indicate an increasing rate of awards before Nobel reception, reaching the summit precisely in the year of the Nobel Prize. After this pinnacle year, awards drop sharply. This result is confirmed by separate analyses of three different disciplines and by a random-effects negative binomial regression model. Such an effect, however, does not emerge for more recent Nobel laureates (1971–2000). In addition, Nobelists in medicine or physiology generate more awards shortly before and after prize reception, whereas laureates in chemistry attract more awards as time progresses.

Machine learning for activity recognition : hip versus wrist data

Problem addressed Wrist-worn accelerometers are associated with greater compliance. However, validated algorithms for predicting activity type from wrist-worn accelerometer data are lacking. This study compared the activity recognition rates of an activity classifier trained on acceleration signal collected on the wrist and hip. Methodology 52 children and adolescents (mean age 13.7 +/- 3.1 year) completed 12 activity trials that were categorized into 7 activity classes: lying down, sitting, standing, walking, running, basketball, and dancing. During each trial, participants wore an ActiGraph GT3X+ tri-axial accelerometer on the right hip and the non-dominant wrist. Features were extracted from 10-s windows and inputted into a regularized logistic regression model using R (Glmnet + L1). Results Classification accuracy for the hip and wrist was 91.0% +/- 3.1% and 88.4% +/- 3.0%, respectively. The hip model exhibited excellent classification accuracy for sitting (91.3%), standing (95.8%), walking (95.8%), and running (96.8%); acceptable classification accuracy for lying down (88.3%) and basketball (81.9%); and modest accuracy for dance (64.1%). The wrist model exhibited excellent classification accuracy for sitting (93.0%), standing (91.7%), and walking (95.8%); acceptable classification accuracy for basketball (86.0%); and modest accuracy for running (78.8%), lying down (74.6%) and dance (69.4%). Potential Impact Both the hip and wrist algorithms achieved acceptable classification accuracy, allowing researchers to use either placement for activity recognition.

The effect of repeated loading and freeze-thaw cycling on immature bovine thoracic motion segment stiffness

There is growing interest in the biomechanics of ‘fusionless’ implant constructs used for deformity correction in the thoracic spine, however, there are questions over the comparability of in vitro biomechanical studies from different research groups due to the various methods used for specimen preparation, testing and data collection. The aim of this study was to identify the effect of two key factors on the stiffness of immature bovine thoracic spine motion segments: (i) repeated cyclic loading and (ii) multiple freeze-thaw cycles, to aid in the planning and interpretation of in vitro studies. Two groups of thoracic spine motion segments from 6-8 week old calves were tested in flexion/extension, right/left lateral bending, and right/left axial rotation under moment control. Group (A) were tested with continuous repeated cyclic loading for 500 cycles with data recorded at cycles 3, 5, 10, 25, 50, 100, 200, 300, 400 and 500. Group (B) were tested after each of five freeze-thaw sequences, with data collected from the 10th load cycle in each sequence. Group A: Flexion/extension stiffness reduced significantly over the 500 load cycles (-22%; P=0.001), but there was no significant change between the 5th and 200th load cycles. Lateral bending stiffness decreased significantly (-18%; P=0.009) over the 500 load cycles, but there was no significant change in axial rotation stiffness (P=0.137). Group B: There was no significant difference between mean stiffness over the five freeze-thaw sequences in flexion/extension (P=0.813) and a near significant reduction in mean stiffness in axial rotation (-6%; P=0.07). However, there was a statistically significant increase in stiffness in lateral bending (+30%; P=0.007). Comparison of in vitro testing results for immature thoracic bovine spine segments between studies can be performed with up to 200 load cycles without significant changes in stiffness. However, when testing protocols require greater than 200 cycles, or when repeated freeze-thaw cycles are involved, it is important to account for the effect of cumulative load and freeze-thaw cycles on spine segment stiffness.

What patient information allows us to make accurate predictions of outcome?

Only some of the information contained in a medical record will be useful to the prediction of patient outcome. We describe a novel method for selecting those outcome predictors which allow us to reliably discriminate between adverse and benign end results. Using the area under the receiver operating characteristic as a nonparametric measure of discrimination, we show how to calculate the maximum discrimination attainable with a given set of discrete valued features. This upper limit forms the basis of our feature selection algorithm. We use the algorithm to select features (from maternity records) relevant to the prediction of failure to progress in labour. The results of this analysis motivate investigation of those predictors of failure to progress relevant to parous and nulliparous sub-populations.