A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PurposeThe conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.MethodsFirst, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8–40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ResultsThe relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ConclusionsThe devised f(0)^CB_w value was calculated by averaging four “point doses” at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.     

Measurement-based model of a wide-bore CT scanner for Monte Carlo dosimetric calculations with GMCTdospp software

The aim of this work was to create a model of a wide-bore Siemens Somatom Sensation Open CT scanner for use with GMCTdospp, which is an EGSnrc-based software tool dedicated for Monte Carlo calculations of dose in CT examinations.The method was based on matching spectrum and filtration to half value layer and dose profile, and thus was similar to the method of Turner et al. (Med. Phys. 36, pp. 2154–2164). Input data on unfiltered beam spectra were taken from two sources: the TASMIP model and IPEM Report 78. Two sources of HVL data were also used, namely measurements and documentation. Dose profile along the fan-beam was measured with Gafchromic RTQA-1010 (QA+) film. Two-component model of filtration was assumed: bow-tie filter made of aluminum with 0.5 mm thickness on central axis, and flat filter made of one of four materials: aluminum, graphite, lead, or titanium.Good agreement between calculations and measurements was obtained for models based on the measured values of HVL. Doses calculated with GMCTdospp differed from the doses measured with pencil ion chamber placed in PMMA phantom by less than 5%, and root mean square difference for four tube potentials and three positions in the phantom did not exceed 2.5%. The differences for models based on HVL values from documentation exceeded 10%. Models based on TASMIP spectra and IPEM78 spectra performed equally well.     

Gold coated contact lens-type ocular in vivo dosimeter (CLOD) for monitoring of low dose in computed tomography: A Monte Carlo study

PurposeThis study reports a sensitivity enhancement of gold-coated contact lens-type ocular in vivo dosimeters (CLODs) for low-dose measurements in computed tomography (CT).MethodsMonte Carlo (MC) simulations were conducted to evaluate the dose enhancement from the gold (Au) layers on the CLODs. The human eye and CLODs were modeled, and the X-ray tube voltages were defined as 80, 120, and 140 kVp. The thickness of the Au layer attached to a CLOD ranged from 100 nm to 10 μm. The thickness of the active layer ranged from 20 to 140 μm. The dose ratio between the active layer of the Au-coated CLOD and a CLOD without a layer, i.e., the dose enhancement factor (DEF), was calculated.ResultsThe DEFs of the first 20-μm thick active layer of the 5-μm thick Au-coated CLOD were 18.4, 19.7, 20.2 at 80, 120, and 140 kVp, respectively. The DEFs decreased as the thickness of the active layer increased. The DEFs of 100-nm to 5-μm thick Au layers increased from 1.7 to 5.4 for 120-kVp X-ray tube voltage when the thickness of the active layer was 140 μm.ConclusionsThe MC results presented a higher sensitivity of Au-coated CLODs (∼20-times higher than that of CLODs without a gold layer). Au-coated CLODs can be applied to an evaluation of very low doses (a few cGy) delivered to patients during CT imaging.     

BioMOCA—a Boltzmann transport Monte Carlo model for ion channel simulation

With the recent availability of high-resolution structural information for several key ion channel proteins and large-scale computational resources, Molecular Dynamics has become an increasingly popular tool for ion channel simulation. However, the CPU requirements for simulating ion transport on time scales relevant to conduction still exceed the resources presently available. To address this problem, we have developed Biology Monte Carlo (BioMOCA), a three-dimensional (3D) coarse-grained particle ion channel simulator based on the Boltzmann Transport Monte Carlo (BTMC) methodology. Although this approach is widely employed in the engineering community to study charge transport in electron devices, its application to molecular biology and electrolytes in general is new and hence must be validated. The pair correlation function, which is a measure of the microscopic structure of matter, provides a suitable benchmark to compare the BTMC method against the well-established Equilibrium Monte Carlo (EMC) approach. For validation purposes BioMOCA is used to simulate several simple homogeneous equilibrium electrolytes at concentrations of physiological interest. The ion–ion pair correlation functions computed from these simulations compare very well with those obtained from EMC simulations. We also demonstrate several performance-improving techniques that result in a several-fold speed-up without compromising the pair correlation function. BioMOCA is then used to perform full 3D simulations of ion transport in the gramicidin A channel in situ in a membrane environment, as well as to study the link between the electrostatic and dielectric properties of the protein and the channel's selectivity.     

Two-step validation of a Monte Carlo dosimetry framework for general radiology

The Monte Carlo technique is considered gold standard when it comes to patient-specific dosimetry. Any newly developed Monte Carlo simulation framework, however, has to be carefully calibrated and validated prior to its use. For many researchers this is a tedious work. We propose a two-step validation procedure for our newly built Monte Carlo framework and provide all input data to make it feasible for future related application by the wider community. The validation was at first performed by benchmarking against simulation data available in literature. The American Association of Physicists in Medicine (AAPM) report of task group 195 (case 2) was considered most appropriate for our application. Secondly, the framework was calibrated and validated against experimental measurements for trunk X-ray imaging protocols using a water phantom. The dose results obtained from all simulations and measurements were compared. Our Monte Carlo framework proved to agree with literature data, by showing a maximal difference below 4% to the AAPM report. The mean difference with the water phantom measurements was around 7%. The statistical uncertainty for clinical applications of the dosimetry model is expected to be within 10%. This makes it reliable for clinical dose calculations in general radiology. Input data and the described procedure allow for the validation of other Monte Carlo frameworks.     

Imaging dose from cone beam computed tomography in radiation therapy

Imaging dose in radiation therapy has traditionally been ignored due to its low magnitude and frequency in comparison to therapeutic dose used to treat patients. The advent of modern, volumetric, imaging modalities, often as an integral part of linear accelerators, has facilitated the implementation of image-guided radiation therapy (IGRT), which is often accomplished by daily imaging of patients. Daily imaging results in additional dose delivered to patient that warrants new attention be given to imaging dose. This review summarizes the imaging dose delivered to patients as the result of cone beam computed tomography (CBCT) imaging performed in radiation therapy using current methods and equipment. This review also summarizes methods to calculate the imaging dose, including the use of Monte Carlo (MC) and treatment planning systems (TPS). Peripheral dose from CBCT imaging, dose reduction methods, the use of effective dose in describing imaging dose, and the measurement of CT dose index (CTDI) in CBCT systems are also reviewed.     

生物组织中有限束宽光吸收的蒙特卡罗模拟

用蒙特卡罗法模拟了有限束宽均匀分布和高斯分布光在生物组织中的传播,分析了生物组织的光学参数及光源特性对光吸收分子的影响,结果表明：只要有效光吸收系数增加,最大光吸收率就会增加,光的侧向传输能力主要依赖于光的有效散射,光束宽度增加,辐照范围加宽,高斯光束的吸收分布的梯度更大。     

Distance-scaled Force Biased Monte Carlo Simulation for Solutions containing a Strongly Interacting Solute

The force-biased extension of the Metropolis Monte Carlo method [1] improves convergence by sampling moves preferentially along the directions of force (and torque) [2]. For solvated systems it is particularly effective [3] when coupled with the preferential sampling scheme [4] that attempts to move solvents near the solute more frequently. However, in recent force-biased simulations of aqueous ionic solutions [5] some of the water molecules in the vicinity of the solute remained essentially stationary. Only significant reduction in the stepsize produced some accepted moves.     

A Monte Carlo software bench for simulation of spectral k-edge CT imaging: Initial results

PurposeSpectral Computed Tomography (SCT) systems equipped with photon counting detectors (PCD) are clinically desired, since such systems provide not only additional diagnostic information but also radiation dose reductions by a factor of two or more. The current unavailability of clinical PCDs makes a simulation of such systems necessary.MethodsIn this paper, we present a Monte Carlo-based simulation of a SCT equipped with a PCD. The aim of this development is to facilitate research on potential clinical applications. Our MC simulator takes into account scattering interactions within the scanned object and has the ability to simulate scans with and without scatter and a wide variety of imaging parameters. To demonstrate the usefulness of such a MC simulator for development of SCT applications, a phantom with contrast targets covering a wide range of clinically significant iodine concentrations is simulated. With those simulations the impact of scatter and exposure on image quality and material decomposition results is investigated.ResultsOur results illustrate that scatter radiation plays a significant role in visual as well as quantitative results. Scatter radiation can reduce the accuracy of contrast agent concentration by up to 15%.ConclusionsWe present a reliable and robust software bench for simulation of SCTs equipped with PCDs.     

Construction of an intermediate-resolution lattice model and re-examination of the helix-coil transition: a dynamic Monte Carlo simulation

In protein modeling, spatial resolution and computational efficiency are always incompatible. As a compromise, an intermediate-resolution lattice model has been constructed in the present work. Each residue is decomposed into four basic units, i.e. the α-carbon group, the carboxyl group, the imino group, and the side-chain group, and each basic coarse-grained unit is represented by a minimum cubic box with eight lattice sites. The spacing of the lattice is about 0.56?Å, holding the highest spatial resolution for the present lattice protein models. As the first report of this new model, the helix-coil transition of a polyalanine chain was examined via dynamic Monte Carlo simulation. The period of formed α-helix was about 3.68 residues, close to that of a natural α-helix. The resultant backbone motion was found to be in the realistic regions of the conformational space in the Ramachandran plot. Helix propagation constant and nucleation constant were further determined through the dynamic hydrogen bonding process and torsional angle variation, and the results were used to make comparison between classical Zimm-Bragg theory and Lifson-Roig theory based on the Qian-Schellman relationship. The simulation results confirmed that our lattice model can reproduce the helix-coil transition of polypeptide and construct a moderately fine α-helix conformation without significantly weakening the priority in efficiency for a lattice model.     

Alternative implementations of Monte Carlo EM algorithms for likelihood inferences

Two methods of computing Monte Carlo estimators of variance components using restricted maximum likelihood via the expectation-maximisation algorithm are reviewed. A third approach is suggested and the performance of the methods is compared using simulated data.     

Monte Carlo simulation of the dose distribution of ICRP adult reference computational phantoms for acquisitions with a 320 detector-row cone-beam CT scanner

PurposeThe purpose of this study was to develop and validate a Monte Carlo (MC) simulation tool for patient dose assessment for a 320 detector-row CT scanner, based on the recommendations of International Commission on Radiological Protection (ICRP). Additionally, the simulation was applied on four clinical acquisition protocols, with and without automatic tube current modulation (TCM).MethodsThe MC simulation was based on EGS4 code and was developed specifically for a 320 detector-row cone-beam CT scanner. The ICRP adult reference phantoms were used as patient models. Dose measurements were performed free-in-air and also in four CTDI phantoms: 150 mm and 350 mm long CT head and CT body phantoms. The MC program was validated by comparing simulations results with these actual measurements acquired under the same conditions. The measurements agreed with the simulations across all conditions within 5%. Patient dose assessment was performed for four clinical axial acquisitions using the ICRP adult reference phantoms, one of them using TCM.ResultsThe results were nearly always lower than those obtained from other dose calculator tools or published in other studies, which were obtained using mathematical phantoms in different CT systems. For the protocol with TCM organ doses were reduced by between 28 and 36%, compared to the results obtained using a fixed mA value.ConclusionsThe developed simulation program provides a useful tool for assessing doses in a 320 detector-row cone-beam CT scanner using ICRP adult reference computational phantoms and is ready to be applied to more complex protocols.     

Optical parameters of embedded abnormalities in tissues as determined by Monte Carlo simulation

Photon propagation through tissue phantoms, made of heart, adipose, and spleen tissues was simulated by Monte Carlo procedure. To detect the presence of deep-seated abnormalities, phantoms of heart with adipose and spleen tissues embedded into it were created and simulations were performed to scan the tissue surface with a source and four detector model. Profiles drawn showed variation in parameters such as backscattered intensity in regions where adipose and spleen tissues were embedded. This study shows that depending on the type of embedded tissue the backscattered fraction as measured at 2 mm from the input fiber is altered. This is enhanced for adipose and decreased for spleen tissue. This is not only shown in scanned profiles on the surface but also in constructed images.     

A Monte Carlo model for the internal dosimetry of choroid plexuses in nuclear medicine procedures

Choroid plexuses are vascular structures located in the brain ventricles, showing specific uptake of some diagnostic and therapeutic radiopharmaceuticals currently under clinical investigation, such as integrin-binding arginine-glycine-aspartic acid (RGD) peptides. No specific geometry for choroid plexuses has been implemented in commercially available software for internal dosimetry.The aims of the present study were to assess the dependence of absorbed dose to the choroid plexuses on the organ geometry implemented in Monte Carlo simulations, and to propose an analytical model for the internal dosimetry of these structures for ¹⁸F, ⁶⁴Cu, ⁶⁷Cu, ⁶⁸Ga, ⁹⁰Y, ¹³¹I and ¹⁷⁷Lu nuclides. A GAMOS Monte Carlo simulation based on direct organ segmentation was taken as the gold standard to validate a second simulation based on a simplified geometrical model of the choroid plexuses. Both simulations were compared with the OLINDA/EXM sphere model.The gold standard and the simplified geometrical model gave similar dosimetry results (dose difference < 3.5%), indicating that the latter can be considered as a satisfactory approximation of the real geometry. In contrast, the sphere model systematically overestimated the absorbed dose compared to both Monte Carlo models (range: 4–50% dose difference), depending on the isotope energy and organ mass. Therefore, the simplified geometric model was adopted to introduce an analytical approach for choroid plexuses dosimetry in the mass range 2–16 g. The proposed model enables the estimation of the choroid plexuses dose by a simple bi-parametric function, once the organ mass and the residence time of the radiopharmaceutical under investigation are provided.     

A Monte Carlo method for the simulation of kinetie models

The purpose of this note is to illustrate the feasibility of simulating kinetic systems, such as commonly encountered in photosynthesis research, using the Monte Carlo (MC) method. In this approach, chemical events are considered at the molecular level where they occur randomly and the macroscopic kinetic evolution results from averaging a large number of such events. Their repeated simulation is easily accomplished using digital computing. It is shown that the MC approach is well suited to the capabilities and resources of modern microcomputers. A software package is briefly described and discussed, allowing a simple programming of any kinetic model system and its resolution. The execution is reasonably fast and accurate; it is not subject to such instabilities as found with the conventional analytical approach.Abbreviations MC Monte Carlo - RN random number - PSU photosynthetic unit Dedicated to Prof. L.N.M. Duysens on the occasion of his retirement.     

An automated Monte Carlo QC system for volumetric modulated arc therapy: Possibilities and challenges

PurposeTo develop and implement an automated Monte Carlo (MC) system for patient specific VMAT quality control in a patient geometry that generates treatment planning system (TPS) compliant DICOM objects and includes a module for 3D analysis of dose deviations. Also, the aims were to recommend diagnose specific tolerance criteria and an evaluation procedure.MethodsThe EGSnrc code package formed the basis for development of the MC system. The workflow consists of a number of modules connected to a TPS by means of manual DICOM exports and imports which were executed sequentially without user interaction. DVH comparison was performed in the TPS. In addition, MC- and TPS dose distributions were analysed by applying the normalized dose difference (NDD) formalism. NDD failure maps and a pass rate for a certain threshold were obtained. 170 clinical plans (prostate, thorax, head-and-neck and gynecological) were selected for analysis.ResultsAgreement within 1.5% was found between clinical- and MC data for the mean dose to the target volumes and within 3% for parameters more sensitive to the shape of the DVH e.g. D_98% PTV. Regarding the NDD analysis, tolerance criteria 2%/3 mm were established for prostate plans and 3%/3 mm for the rest of the cases.ConclusionsAn automated MC system was developed and implemented. Evaluation procedure is recommended with NDD-analysis as a first step. For pass rate < 95%, the evaluation continues with comparison of DVH parameters. For deviations larger than 2%, a visual inspection of the clinical- and MC dose distributions is performed.     

Monte Carlo commissioning of radiotherapy LINAC—Introducing an improved methodology

PurposeMonte Carlo (MC) commissioning of medical linear accelerator (LINAC) is a time-consuming process involving a comparison between measured and simulated cross beam/lateral profiles and percentage depth doses (PDDs) for various field sizes. An agreement between these two data sets is sought by trial and error method while varying the incident electron beam parameters, such as electron beam energy or width, etc. This study aims to improve the efficiency of MC commissioning of a LINAC by assessing the feasibility of using a limited number of simulated PDDs.Materials and methodsUsing EGSnrc codes, a Varian Clinac 2100 unit has been commissioned for 6 MV photon beam, and a methodology has been proposed to identify the incident electron beam parameters in a speedier fashion. Impact of voxel size in 3-dimensions and cost functions used for comparison of the measured and simulated data have been investigated along with the role of interpolation.ResultsA voxel size of 1 × 1×0.5 cm³ has been identified as suitable for accurate and fast commissioning of the LIANC. The optimum number of simulated PDDs (required for further interpolation) has been found to be five.ConclusionThe present study suggests that PDDs alone at times can be insufficient for an unambiguous commissioning process and should be supported by including the lateral beam profiles in the process.