Accelerating Coupled-Cluster Calculations with GPUs: An Implementation of the Density-Fitted CCSD(T) Approach for Heterogeneous Computing Architectures Using OpenMP Directives

An algorithm is presented for the coupled-cluster singles, doubles, and perturbative triples correction [CCSD(T)] method based on the density fitting or the resolution-of-the-identity (RI) approximation for performing calculations on heterogeneous computing platforms composed of multicore CPUs and graphics processing units (GPUs). The directive-based approach to GPU offloading offered by the OpenMP application programming interface has been employed to adapt the most compute-intensive terms in the RI-CCSD amplitude equations with computational costs scaling as O ( N O 2 N V 4 ) , O ( N O 3 N V 3 ) , and O ( N O 4 N V 2 ) (where NO and NV denote the numbers of correlated occupied and virtual orbitals, respectively) and the perturbative triples correction to execute on GPU architectures. The pertinent tensor contractions are performed using an accelerated math library such as cuBLAS or hipBLAS. Optimal strategies are discussed for splitting large data arrays into tiles to fit them into the relatively small memory space of the GPUs, while also minimizing the low-bandwidth CPU-GPU data transfers. The performance of the hybrid CPU-GPU RI-CCSD(T) code is demonstrated on pre-exascale supercomputers composed of heterogeneous nodes equipped with NVIDIA Tesla V100 and A100 GPUs and on the world's first exascale supercomputer named "Frontier", the nodes of which consist of AMD MI250X GPUs. Speedups within the range 4-8× relative to the recently reported CPU-only algorithm are obtained for the GPU-offloaded terms in the RI-CCSD amplitude equations. Applications to polycyclic aromatic hydrocarbons containing 16-66 carbon atoms demonstrate that the acceleration of the hybrid CPU-GPU code for the perturbative triples correction relative to the CPU-only code increases with the molecule size, attaining a speedup of 5.7× for the largest circumovalene molecule (C66H20). The GPU-offloaded code enables the computation of the perturbative triples correction for the C60 molecule using the cc-pVDZ/aug-cc-pVTZ-RI basis sets in 7 min on Frontier when using 12,288 AMD GPUs with a parallel efficiency of 83.1%.

The General Atomic and Molecular Electronic Structure System (GAMESS): Novel Methods on Novel Architectures

The primary focus of GAMESS over the last 5 years has been the development of new high-performance codes that are able to take effective and efficient advantage of the most advanced computer architectures, both CPU and accelerators. These efforts include employing density fitting and fragmentation methods to reduce the high scaling of well-correlated (e.g., coupled-cluster) methods as well as developing novel codes that can take optimal advantage of graphical processing units and other modern accelerators. Because accurate wave functions can be very complex, an important new functionality in GAMESS is the quasi-atomic orbital analysis, an unbiased approach to the understanding of covalent bonds embedded in the wave function. Best practices for the maintenance and distribution of GAMESS are also discussed.

Dose Escalation Using Magnetic Resonance Guided High-Dose-Rate Endorectal Brachytherapy to Enhance Clinical Response after Neoadjuvant Radiotherapy in Rectal Adenocarcinoma

PURPOSE/OBJECTIVE(S)

To assess the proportion of patients with rectal adenocarcinoma achieving clinical complete response after neoadjuvant chemoradiation (NACTRT) and MR guided endorectal brachytherapy boost (MR-ERBT) MATERIALS/METHODS: Patients with rectal cancers (T2-T4/N0-N+) treated with concurrent chemoradiotherapy (50Gy/25# with Capecitabine) between June-2017 to April-2022. Post RT, patients having residual non-circumferential lesions <8cm in length were administered escalated-dose MR-ERBT with Ir192 HDR source. A median dose of 12Gy (8-15Gy) in 3 (2-3) fractions at 3-5 day intervals was delivered using MR-ERBT after external radiation. Data on near complete/complete clinical response (nCR/cCR) rates, local regrowth rates and clinical outcome were collected for analysis.

RESULTS

Of the 145 patients who received MR-ERBT, majority were staged as T3(78.7%) and N1-2 (77.5%) rectal cancers. Median tumor length was 4cm and 123 (85%) of the tumors were located in the lower rectum (0-5cm from anal verge). Seventy-six (52%) patients achieved cCR or nCR (37 cCR, 39 nCR) and were advised observation or watch and wait (WW) management. The 69 patients having partial response were advised surgery. The patients having nCR 16 (29%) underwent resection, of these 10 (62%) had pathological complete response (pCR). The patients with partial or poor response, 57 underwent resection and of these 11 (19%) had pCR, 12 patients refused surgery due to fear of permanent stoma and continued to be on follow up. Of the 79 patients undergoing resection, 36 (45.5%) had sphincter preserving surgeries. At the median follow up of 30 months, local regrowth was seen in 8 (10.5%) of patients on WW and 6 were surgically salvaged while other 2 had synchronous metastatic relapse. Thus, 56 (41%) achieved organ preservation and continued to be on WW management. Twelve (8,2%) patients developed distant metastasis in the entire cohort, 3 in the WW group and 9 in the resected group. There were no pelvic recurrences seen in the resected patients. The disease-free survival at 3 years were (96.1% vs 89% Observation vs. resected (p_0.05) respectively. The overall survival at 3 years were (93% vs 98% Observation vs. resected (p_0.44) respectively. Late rectal toxicity was observed in 16(11%) patients on observation CONCLUSION: Dose Escalated MR-EBRT is an effective and safe method to enhance complete clinical response, thus improving the rate of organ preservation for distal rectal cancers.

Porting fragmentation methods to GPUs using an OpenMP API: Offloading the resolution-of-the-identity second-order Møller-Plesset perturbation method

Using an OpenMP Application Programming Interface, the resolution-of-the-identity second-order Møller-Plesset perturbation (RI-MP2) method has been off-loaded onto graphical processing units (GPUs), both as a standalone method in the GAMESS electronic structure program and as an electron correlation energy component in the effective fragment molecular orbital (EFMO) framework. First, a new scheme has been proposed to maximize data digestion on GPUs that subsequently linearizes data transfer from central processing units (CPUs) to GPUs. Second, the GAMESS Fortran code has been interfaced with GPU numerical libraries (e.g., NVIDIA cuBLAS and cuSOLVER) for efficient matrix operations (e.g., matrix multiplication, matrix decomposition, and matrix inversion). The standalone GPU RI-MP2 code shows an increasing speedup of up to 7.5× using one NVIDIA V100 GPU with one IBM 42-core P9 CPU for calculations on fullerenes of increasing size from 40 to 260 carbon atoms using the 6-31G(d)/cc-pVDZ-RI basis sets. A single Summit node with six V100s can compute the RI-MP2 correlation energy of a cluster of 175 water molecules using the correlation consistent basis sets cc-pVDZ/cc-pVDZ-RI containing 4375 atomic orbitals and 14 700 auxiliary basis functions in ∼0.85 h. In the EFMO framework, the GPU RI-MP2 component shows near linear scaling for a large number of V100s when computing the energy of an 1800-atom mesoporous silica nanoparticle in a bath of 4000 water molecules. The parallel efficiencies of the GPU RI-MP2 component with 2304 and 4608 V100s are 98.0% and 96.1%, respectively.

Quantum Chemical Modeling of Propellant Degradation

An ab initio quantum chemical approach for the modeling of propellant degradation is presented. Using state-of-the-art bonding analysis techniques and composite methods, a series of potential degradation reactions are devised for a sample hydroxyl-terminated-polybutadiene (HTPB) type solid fuel. By applying thermochemical procedures and isodesmic reactions, accurate thermochemical quantities are obtained using a modified G3 composite method based on the resolution of the identity. The calculated heats of formation for the different structures produced presents an ∼2 kcal/mol average error when compared against experimental values.

PDG: A Composite Method Based on the Resolution of the Identity

The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller-Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles-doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the "standard" G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. The efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4-5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

A parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail O(NO2NV4) and O(NO3NV3)-scaling computational costs, where N_O and N_V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N_V⁴ and N_O²N_V²-scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24-51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014, 16, 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

Computation of epistemic uncertainty due to limited data samples in small field dosimetry using Fuzzy Set Theory

OBJECTIVE

To estimate the epistemic (or fuzzy) uncertainty, arising due to limited data samples in the measurement of the output factors (OFs) of the small fields using Fuzzy Set Theory (FST).

METHODS

EBT3 film samples of size 50 × 50 mm² were used for the measurement of the OF of stereotactic radiosurgery (SRS) cones of size 4, 6, 7.5, 10, 12.5 and 15 mm diameter, normalized with respect to the output of 100 × 100 mm² open field size. Three measurements were done per cone/field size. Red color channel was chosen for the dosimetry purpose, net optical density (NOD) was converted to the dose using non-linear relation. To estimate the epistemic uncertainty associated with the measured OFs due to limited number of data samples, a triangular fuzzy number (TFN) was assumed as the fuzziness in the dose delivered by the individual SRS cone/field. Uncertainty in the OF was estimated by applying the Fuzzy Vertex Method (FVM). The membership functions of the OF were constructed for each cone size and the nature of the uncertainty in the OF of the cones was expressed in the terms of its fuzziness. For the sake of completeness of the study, the statistical uncertainty involved in the procedure has also been calculated.

RESULTS

The statistical and fuzzy uncertainties in the measurement of OF of cones range from 3.28 to 6.25% and 2.58 to 5.44% respectively. The smallest cone of 4 mm has the largest values of statistical and fuzzy uncertainties. The membership functions of the OF for the studied cones were triangular in nature.

CONCLUSION

The epistemic uncertainty arising due to limited number of data samples holds a significant fraction of the prescribed dose, and therefore, should not be ignored in the total uncertainty estimation.

ADVANCES IN KNOWLEDGE

This study highlights the significance of epistemic component of measurement uncertainty arising out due to the insufficient/limited number of measurements of a quantity.

57Fe Mössbauer parameters from domain based local pair-natural orbital coupled-cluster theory

We report on applications of the domain based local pair-natural orbital (PNO) coupled-cluster method within the singles and doubles approximation (DLPNO-CCSD) to the calculation of ⁵⁷Fe isomer shifts and quadrupole splittings in a small training set of iron complexes consisting of large molecular ligands and iron atoms in varying charge, spin, and oxidation states. The electron densities and electric field gradients needed for these calculations were obtained within the recently implemented analytic derivative scheme. A method for the direct treatment of scalar relativistic effects in the calculation of effective electron densities is described by using the first-order Douglas-Kroll-Hess Hamiltonian and a Gaussian charge distribution model for the nucleus. The performance of DLPNO-CCSD is compared with four modern-day density functionals, namely, RPBE, TPSS, B3LYP, and B2PLYP, as well as with the second-order Møller-Plesset perturbation theory. An excellent correlation between the calculated electron densities and the experimental isomer shifts is attained with the DLPNO-CCSD method. The correlation constant a obtained from the slope of the linear correlation plot is found to be ≈-0.31 a.u.³ mm s^-1, which agrees very well with the experimental calibration constant α = -0.31 ± 0.04 a.u.³ mm s^-1. This value of a is obtained consistently using both nonrelativistic and scalar relativistic DLPNO-CCSD electron densities. While the B3LYP and B2PLYP functionals achieve equally good correlation between theory and experiment, the correlation constant a is found to deviate from the experimental value. Similar trends are observed also for quadrupole splittings. The value of the nuclear quadrupole moment for ⁵⁷Fe is estimated to be 0.15 b at the DLPNO-CCSD level. This is consistent with previous results and is here supported by a higher level of theory. The DLPNO-CCSD results are found to be insensitive to the intrinsic approximations in the method, in particular the PNO occupation number truncation error, while the results obtained with density functional theory (DFT) are found to depend on the choice of the functional. In a statistical sense, i.e., on the basis of the linear regression analysis, however, the accuracies of the DFT and DLPNO-CCSD results can be considered comparable.

Numerical simulation of 222Rn profiling in an experimental chamber using CFD technique

Measurement of indoor ²²²Rn concentration and interpretation of distribution patterns are important for inhalation dosimetry in occupational and residential areas. Experimental determination of ²²²Rn concentration distribution and estimation of inhalation doses depend on the underlying aspects such as calibration of the detectors, accuracy of the techniques etc. Therefore, ²²²Rn concentration distribution needs to be very well understood in a closed domain for the controlled studies. In the recent times, Computational fluid dynamics (CFD) technique has gained a lot of attention for the prediction and visualization of indoor ²²²Rn concentration profiles and their mixing ability in the domain. The present study aims to simulate the effect of forced mixing on the ²²²Rn concentration profile in a 22 m³ experimental chamber. This chamber is designed for carrying out the controlled experiments, calibration and inter-comparison studies of various types of ²²²Rn detectors. Effect of different parameters such as time, flow rates, fan-off and fan-on conditions have been studied on the transient response, extent of the air mixing patterns and subsequently on ²²²Rn concentration profile in the chamber. Further, Non uniformity index (NUI) is introduced as a measure of the uniformity of the distribution in the closed domain. NUI is estimated for different cases in order to efficiently interpret the effect of above mentioned parameters on ²²²Rn profile in the chamber. This study will be useful to represent the turbulent conditions in real indoor domains and occupational facilities as U-mines during calibration and inter-comparison exercises of different ²²²Rn detectors.

ESTIMATION OF UNCERTAINTY IN MEASUREMENT OF DOSE EQUIVALENT AT LABORATORY LEVEL USING CASO4 :Dy-BASED TLD BADGE SYSTEM IN INDIA

The objective of this paper is to estimate the combined uncertainty in the measurement of dose equivalent at laboratory level using CaSO4:Dy-based thermoluminescent dosemeter badge system by including variations in the components of the system. The variability of performance of the system is analysed using random effects one way analysis of variance model. The model enables estimation of the overall variance of the performance of the sampled population. The population in the study comprises all possible indicated dose equivalents on irradiation of dosemeters to a specific dose equivalent and radiation quality. Coefficient of variation and combined uncertainty at 95% level of confidence in the measurement of Hp(10) due to S-Cs radiation quality are found to be 6.6 and 14.3%, respectively, at the dose level of 5.31 mSv. The above parameters in the measurement of in-use quantity, i.e. whole body dose or photon dose equivalent are found to be 7.4 and 16.4%, respectively. The performance of the monitoring system on relative response has been observed to be satisfactory. Various factors affecting the variability of performance of the system are identified for further improvement in coefficient of variation.

Recent developments in the general atomic and molecular electronic structure system

A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.

Investigation of 220Rn emanation and exhalation from soil samples of Larsemann Hills region, Antarctica

In the present study, thoron exhalation flux density were measured in the soil samples collected around the Indian station namely Bharati (69° 24.41' S, 76° 11.72' E) and its nearby islands in the Larsemann hills region of Antarctica. Further, dependency of thoron mass emanation rate and emanation coefficient on the soil grain size was studied by segregating the soil samples into four different grain size groups: 50-100 μm, 100-200 μm, 200-500 μm and 500-1000 μm which showed that both of them follow a decreasing trend with increase in grain size. A comparison of measured mass emanation rate between different soil samples showed that it had a larger variation for the smaller grain size which eventually decreased as grain size increased while emanation coefficient was observed to be nearly constant for all the grain size groups. The variation in emanation coefficient with respect to mean grain size has been investigated and an empirical exponential model has been proposed for predicting emanation coefficient for different grain sizes.

MONTE CARLO-BASED INVESTIGATION OF MICRODOSIMETRIC DISTRIBUTION OF HIGH ENERGY BRACHYTHERAPY SOURCES

FLUKA-based Monte Carlo calculations were carried out to study microdosimetric distributions in air and in water for encapsulated high energy brachytherapy sources (60Co, 137Cs, 192Ir and 169Yb) by simulating a Tissue Equivalent Proportional Counter (Model LET1/2) having sensitive diameter of 1. 27 cm for a site size of 1 μm. The study also included microdosimetric distributions of bare sources. When the sources are in air, for a given source, the source geometry does not affect the y¯F and y¯D values significantly. When the encapsulated 192Ir, 137Cs and 60Co sources are in water, y¯F and y¯D values increase with distance in water which is due to degradation in the energy of photons. Using the calculated values of y¯D, relative biological effectiveness (RBE) was obtained for the investigated sources. When 60Co, 137Cs and 192Ir sources are in water, RBE increases from 1.03 ± 0.01 to 1.17 ± 0.01, 1.24 ± 0.01 to 1.46 ± 0.02 and 1.50 ± 0.01 to 1.75 ± 0.03, respectively, when the distance was increased from 3-15 cm, whereas for 169Yb, RBE is about 2, independent of distance in water.

STUDY OF EFFECT OF CONSECUTIVE HEATING ON THERMOLUMINESCENCE GLOW CURVES OF MULTI-ELEMENT TL DOSEMETER IN HOT GAS-BASED READER SYSTEM

The objective of this paper is to study the effect of consecutive heating of TL elements of a thermoluminescence dosemeter (TLD) card in hot N2 gas-based TLD badge reader. The effect is studied by theoretical simulations of clamped heating profiles of the discs and resulting TL glow curves. The simulated temperature profile accounts for heat transfer to disc from hot gas as well as radiative and convective heat exchanges between the disc and the surrounding. The glow curves are simulated using 10 component glow peak model for CaSO4:Dy using the simulated temperature profile. The shape of the simulated glow curves and trend in total TL signal of the three discs were observed to match closely with the experimental observations when elevated surrounding temperature was considered for simulation. It is concluded that the readout (heating) of adjacent TLD disc affects the surrounding temperature leading to the changes in temperature profile of the next disc.

MEASUREMENT OF OPERATIONAL QUANTITIES Hp(0.07) AND Hp(3) FOR INDIGENOUSLY DEVELOPED 106Ru/106Rh SOURCE USING AN EXTRAPOLATION CHAMBER

In the present study, a prototype 106Ru/106Rh source was fabricated using high level liquid waste from reactor fuel, fixed in a stainless steel housing with a window and backing made of silver. The study involves measurement of the operational quantities Hp(0.07), Hp(3) and the percentage depth dose (PDD) using an extrapolation chamber. It also involves determination of necessary correction factors to arrive at Hp(0.07) and Hp(3) following International Organisation for Standardisation (ISO) and methods suggested in literature. The study facilitates incorporation of the 106Ru/106Rh source as a beta reference source for quality assurance programme in TLD personnel monitoring as per the guidelines of ISO.

INVESTIGATION OF APPLICABILITY OF PURE PROPANE GAS FOR MICRODOSIMETRY AT NEUTRON FIELDS: A MONTE CARLO STUDY

Applicability of pure propane gas for microdosimetric measurements in neutron fields was investigated using the FLUKA Monte Carlo code. Monoenergetic neutrons in the energy range 1 keV-20 MeV and the ISO-neutron sources such as 241Am-Be, 241Am-B, 252Cf and 252Cf + D2O were considered in the present study. The tissue-equivalent proportional counter (TEPC) simulated in the study was LET-1/2 (by Far West Technology) with site sizes 1, 2 and 8 μm. The study demonstrates that for a given site size, the TEPC filled with tissue-equivalent propane and pure propane gases produce similar microdosimetric distributions when the density of pure propane gas is lowered appropriately. For the ISO-neutron sources, the density of propane gas requires scaling by a factor 0.85. For the monoenergetic neutrons, depending upon the neutron energy, the values of scaling factors are in the range of 0.58-0.93.

Phantom-Based Feasibility Studies on Phase-Contrast Mammography at Indian Synchrotron Facility Indus-2

Introduction:

Use of synchrotron radiation (SR) X-ray source in medical imaging has shown great potential for improving soft-tissue image contrast such as the breast. The present study demonstrates quantitative X-ray phase-contrast imaging (XPCI) technique derived from propagation-dependent phase change observed in the breast tissue-equivalent test materials.

Materials and Methods:

Indian synchrotron facility (Indus-2, Raja Ramanna Centre of Advanced Technology [RRCAT]) was used to carry out phantom feasibility study on phase-contrast mammography. Different phantoms and samples, including locally fabricated breast tissue-equivalent phantoms were used to perform absorption and phase mode imaging using 12 and 16 keV SR X-ray beam. Edge-enhancement index (EEI) and edge enhancement to noise ratio (EE/N) were measured for all the images. Absorbed dose to air values were calculated for 12 and 16 keV SR X-ray beam using the measured SR X-ray photon flux at the object plane and by applying the standard radiation dosimetry formalism.

Results and Conclusion:

It was observed in case of all the phantoms and test samples that EEI and EE/N values are relatively higher for images taken in the phase mode. The absorbed dose to air at imaging plane was found to be 75.59 mGy and 28.9 mGy for 12 and 16 keV SR energies, respectively. However, these dose values can be optimized by reducing the image acquisition time without compromising the image quality when clinical samples are imaged. This work demonstrates the feasibility of XPCI in mammography using 12 and 16 keV SR X-ray beams.

Accurate Prediction of Hyperfine Coupling Tensors for Main Group Elements Using a Unitary Group Based Rigorously Spin-Adapted Coupled-Cluster Theory

We present the development of a perturbative triples correction scheme for the previously reported unitary group based spin-adapted combinatoric open-shell coupled-cluster (CC) singles and doubles (COS-CCSD) approach and report on the applications of the newly developed method, termed "COS-CCSD(T)", to the calculation of hyperfine coupling (HFC) tensors for radicals consisting of hydrogen, second- and third-row elements. The COS-CCSD(T) method involves a single noniterative step with [Formula: see text] scaling of the computational cost for the calculation of triples corrections to the energy. The key feature of this development is the use of spatial semicanonical orbitals generated from standard restricted open-shell Hartree-Fock (ROHF) orbitals, which allows the unperturbed Hamiltonian operator to be defined in terms of a diagonal spin-free Fock operator. The HFC tensors are computed as a first-order property via implementation of an analytic derivative scheme. The required one-particle spin density matrix is computed by using one- and two-particle spin-free density matrices that are obtained from the analytic derivative implementation, in this way avoiding the use of any spin-dependent operator and maintaining spin adaptation of the CC wavefunction. Benchmark calculations of HFC tensors for a set of 21 radicals indicate reasonably good agreement of the COS-CCSD(T) results with experiment and a consistent improvement over the COS-CCSD method. We demonstrate that the accuracies of the isotropic hyperfine coupling constants obtained in unrestricted HF (UHF) reference based spin-orbital CCSD(T) calculations deteriorate when spin contamination in the UHF wavefunction is large, and the results may even become qualitatively incorrect when spin polarization is the driving mechanism. Within a similar noniterative perturbative treatment of triple excitations, the spin-adapted COS-CCSD(T) approach produces accurate results, thus ensuring cost-effectiveness together with reliability.

Dosimetry of indigenously developed 177Lu patch source for surface brachytherapy-Experimental and Monte Carlo methods

This paper describes the evaluation of dosimetry characteristics of an in-house developed ¹⁷⁷Lu skin patch source for treatment of non-melanoma skin cancer. A ¹⁷⁷Lu skin patch source based on Nafion-115 membrane backbone containing 3.46 ± 0.01 mCi of activity was used. Activity measurement of the patch source was based on gamma ray spectrometry using a HPGe detector. The efficiencies of the HPGe detector were fitted using an orthogonal polynomial function. The absorbed dose rate to water at 5 μm depth in water was determined using an extrapolation chamber, EBT3 Gafchromic film and compared with Monte Carlo methods. The correction factors such as Bragg-Gray stopping power ratio of water-to-air and chamber wall material being different from water, needed to be applied on measurements for establishing the dose rate at 5 μm depth, were calculated using the Monte Carlo method. Absorbed dose rate at 5 μm depth in water (surface dose rate) measured using an extrapolation chamber and EBT3 Gafchromic film were 9.9 ± 0.7 and 8.2 ± 0.1 Gy h^-1 mCi^-1 respectively for the source activity of 3.46 ± 0.01 mCi. The surface dose rate calculated using the Monte Carlo method was 8.7 ± 0.2 Gy h^-1 mCi^-1, which agrees reasonably well with measurement. The measured dose rate per mCi offers scope for ascertaining treatment time required to deliver the dose for propitious therapeutic outcome. Additionally, on-axis depth dose and lateral dose profiles at 5 μm and 1 mm depth in water phantom were also calculated using the Monte Carlo method.

Pretreatment Dose Verification in Volumetric Modulated Arc Therapy Using Liquid Ionization Chamber

Purpose:

The purpose of the present study was to evaluate the practicability of liquid ionization chamber (LIC) for pretreatment dose verification of the advanced radiotherapy techniques such as volumetric modulated arc therapy (VMAT).

Materials and Methods:

The dosimetric characteristics of LIC such as repeatability, sensitivity, monitor unit linearity, dose rate dependence, angular dependence, voltage-current response, and output factors were investigated in 6 MV therapeutic X-ray beams. The LIC was cross-calibrated against 0.125-cc air-filled thimble ionization chamber. A dedicated dosimetry insert made up of Perspex to incorporate the LIC at proper location in the intensity-modulated radiation therapy thorax phantom was locally fabricated. The collection efficiency and ion recombination correction factor was determined using the two-dose rate method. Pretreatment dose verification measurement of VMAT treatment plans were carried out using the liquid ionization chamber as well as small volume (0.125 cc) air-filled thimble ionization chamber. The measured dose values by the two dosimeters and TPS calculated dose at a given point were compared.

Results:

The relative percentage differences between the TPS calculated and measured doses were within ± 1.57% for LIC and ± 2.21% for 0.125 cc ionization chamber, respectively.

Conclusions:

The measured dose values by the two dosimeters and TPS calculated dose at a given point were found comparable suggesting that the LIC could be a good choice of dosimeter for pretreatment dose verification in VMAT.

Imaging and Dosimetric Study on Direct Flat-Panel Detector-Based Digital Mammography System

INTRODUCTION

Image quality of digital mammography system is generally defined by three primary physical parameters, namely, contrast, resolution, and noise. Quantification of these metrics can be done by measuring objective image quality parameters defined as contrast-to-noise ratio (CNR), modulation transfer function (MTF), and noise power spectra (NPS).

MATERIALS AND METHODS

In the present study, various imaging metrics such as CNR, contrast detail resolution, MTF, and NPS were evaluated for a direct flat-panel detector-based digital mammography system following the European Guidelines. Furthermore, system performance relating to both image quality and doses were evaluated using figure of merit (FOM) in terms of CNR2/mean glandular dose (MGD) under automatic exposure control (AEC) and clinically used OPDOSE operating mode.

RESULTS AND CONCLUSION

Under AEC mode, FOM values for the 4.5 cm thick BARC polymethyl methacrylate (PMMA) phantom were found to be 15.02, 15.88, and 19.82 at Mo/Mo, Mo/Rh, and W/Rh target/filter (T/F), respectively. Under OPDOSE mode, FOM values were found to 65.32, 11.80, and 1.14 for the BARC PMMA phantom thickness of 2, 4.5, and 8 cm, respectively. Under OPDOSE mode, the calculated MGD values for three Computerized Imaging Reference Systems slab phantoms having total thickness of 7 cm were observed to be 3.03, 2.32, and 1.75 mGy with glandular/adipose tissue compositions of 70/30, 50/50, and 30/70, respectively, whereas for the 2-8-cm thick BARC PMMA phantom, the calculated MGDs were found to be in the range of 0.57-3.32 mGy. All the calculated MGDs values were found to be lower than the acceptable level of dose limits provided in European Guidelines.

Effect of quantum confinement on lifetime of anharmonic decay of optical phonons in semiconductor nanostructures

The anharmonic decay of phonons underlies many important effects in semiconductors, e.g. hotspot formation, phonon bottleneck and thermal resistivity. In this article, we evaluate the effect of quantum confinement on the anharmonic decay transition probability in a cubic isotropic material. This article focuses on the anharmonic decay of longitudinal optical phonons, generated from hot electrons, are directly related to formation of hotspots in the active region of semiconductor devices. The confinement effect has been realized in double interface heterostructure quantum well (DHSQW) (e.g. AlAs/GaAs/AlAs) and free-standing quantum well (FSQW) (e.g. GaAs) structures as the confined phonon modes have different properties inside the structures. The longitudinal-optical phonon decay rate is reduced for the case of a DHSQW compared to bulk material and for a FSQW the decay rate has a strong dependence on wavevector value of the three phonons involved.

DEVELOPMENT OF AN ALGORITHM TO ESTIMATE EYE LENS DOSE IN TERMS OF OPERATIONAL QUANTITY Hp(3) USING HEAD TLD BADGE

In view of the recommendations of International Commission on Radiological Protection for reduction of the occupational annual dose limit for eye lens from 150 mSv to 20 mSv/y, questions have been raised on the adequacy of monitoring for the quantities Hp(10) and Hp(0.07). As an immediate requirement, in the present situation, where there is no exclusive eye lens dosemeter in India, the existing chest TLD badge was modified to be used as head badge (head dosemeter) by including a strap to enable wearing on the forehead. In order to estimate the eye lens dose in terms of the operational quantity Hp(3), the prevalent algorithm of chest badge was also modified. The modified algorithm was applied to estimate Hp(3) for dosemeters irradiated to various beta and photon radiations including mixtures. The Q values (estimated/delivered dose equivalent) were found to be within ±20% for most of the photon beams.

FLUKA-BASED MONTE CARLO INVESTIGATION OF MICRODOSIMETRIC DISTRIBUTIONS OF TELECOBALT BEAM

FLUKA-based Monte Carlo calculations of microdosimetric distributions in water phantom involving a walled spherical Tissue-Equivalent Proportional Counter filled with tissue-equivalent propane gas have been studied for an indigenously developed telecobalt machine. The simulated site size considered in the study was 2 μm. In the Monte Carlo calculations, field size was varied from 10 cm × 10 cm to 35 cm × 35 cm and the depth was varied as 5-20 cm. The study also includes calculation of microdosimetric distributions with a 30° wedge filter. The efficiency of the calculations was improved up to a factor of 26 by choosing appropriate cut off values for production and transport of electrons. The calculated microdosimetric distributions of telecobalt machine is distinctly different from that of a bare 60Co source which is attributed to the influence of scattered photons from the machine head and the water phantom.

Study of Residual TL of CaSO4:Dy-based Thermoluminescence Dosemeter

Thermoluminescence (TL) signal obtained during second readout of a TL dosemeter disc previously exposed to ionising radiation is termed as residual TL. The origin of residual TL has not been discussed in detail so far in the literature. In this work, experimentally obtained residual TL signal and its origin is studied for CaSO4:Dy-based TL dosemeter through numerical simulation and role of deeper traps has been discussed.

FEASIBILITY STUDY FOR IDENTIFICATION OF STATIC AND DYNAMIC EXPOSURE USING CCD IMAGING TECHNIQUE FOR Caso4:Dy TL DOSEMETERS

The occupational exposure incurred by the radiation workers due to the external radiation is estimated using personal dosemeter placed on the human body during the monitoring period. In certain situations, it is required to determine whether the dosemeter alone was exposed accidentally/intentionally in radiation field (static exposure) or was exposed while being worn by a worker moving in his workplace (dynamic exposure). The present thermoluminscent (TL) based personnel monitoring systems are not capable of distinguishing between the above stated (static and dynamic) exposure conditions. The feasibility of a new methodology developed using the charge coupled device based imaging technique for identification of the static/dynamic exposure of CaSO4:Dy based TL detectors for low energy photons has been investigated. The techniques for the qualitative and the quantitative assessments of the exposure conditions are presented in this paper.

DEEP, SHALLOW AND EYE LENS DOSES FROM 106Ru/106Rh-A COMPARSION

106Ru/106Rh is unique amongst other commonly used beta sources such as 147Pm, 85Kr, 204Tl, 32P, natU and 90Sr/90Y in the sense that it is capable of simultaneously delivering shallow/skin, eye lens and deep/whole body doses (WBDs) and they differ from each other substantially. In view of this, the investigation of various quantities defined for individual monitoring is possible and this makes 106Ru/106Rh beta source, a classical example in radiation protection and dosimetry. This led us to estimate skin, eye lens and WBDs for 106Ru/106Rh beta source. Optically stimulated luminescence based ultra-thin α-Al2O3:C disc dosimeters were used in the present study. Typical values (relative) of the eye lens and whole body/deep doses with respective to the skin dose (100%) were experimentally measured as ~66 ± 4.6% and 17 ± 3.9%, respectively. The study shows that 106Ru/106Rh beta source is capable of delivering even WBD which is not the case with other beta sources.

Optimization of beam dump shielding for K-130 cyclotron at VECC

A compact and efficient beam dump shield has been designed using Monte Carlo simulation code FLUKA to facilitate low background measurement of neutron and gamma rays using K130 cyclotron at Variable Energy Cyclotron Centre, Kolkata (VECC). Iron, lead and high density Polyethylene (HDPE) were considered in the design of the beam dump shield. Representative FLUKA simulation results have been validated using in-beam experiment performed on the same beam dump constituents. Experimental neutron and gamma-rays energy spectra have been found to be in fair agreement with the simulation results. Activation of various beam dump shield components were also carried out.

Quick, efficient and effective patient-specific intensity-modulated radiation therapy quality assurance using log file and electronic portal imaging device

AIM

The aim of work is to explore a quick, efficient, and effective patient-specific intensity-modulated radiation therapy (IMRT) quality assurance (QA).

MATERIALS AND METHODS

Software tools were developed to extract and analyze the multi-leaf collimator (MLC) leaf positions (LPs) from electronic portal imaging device (EPID) images for Varian C-series machine and TrueBeam, to extract useful data from MLC log file of C-series linear accelerator (LINAC), to extract useful information from the trajectory log binary file of TrueBeam LINAC, to compare LPs derived from EPID images with log file/trajectory log data, and to analyze IMRT treatment files using the MATLAB programming language. The difference in LP determined from the trajectory log and EPID images was proposed for patient-specific QA.

RESULTS

It was found that the differences in LP for regular radiation fields generated using stationary leaves are <0.5 mm for all the field sizes while for regular radiation fields generated using the moving leaves are more but <2 mm. The differences in LPs for IMRT field were also determined and found to be <2 mm.

CONCLUSIONS

The methodology demonstrated can be used for establishing the accuracy of trajectory log data and for independent routine IMRT QA by generating single number like gamma index to indicate pass or fail of an IMRT treatment plan. The QA indices such as numbers of occurrences of ≥2 mm error in LPS are found more than 5% of total number of occurrences; the dosimetric review of planned treatment is advisable.

BETA DOSE EVALUATION ALGORITHM FOR CaSO4:Dy BYSED TLD BADGE USED IN COUNTRYWIDE PERSONNEL MONITORING PROGRAMME IN INDIA

The prevalent algorithm for the estimation of beta dose, which in turn is used for estimation of skin dose for exposures involving beta radiations was observed to significantly overestimate the dose in individual monitoring based on CaSO4:Dy TLD badge in India. A new algorithm has been developed by estimating the correction factor from the response of dosemeter to different beta sources at various angles of incidence. The correction factor was observed to vary linearly with the ratio of the responses of dosemeter element without filter (D3) and dosemeter element under Perspex filter (D2). The correction factor determined using the ratio of D3 and D2 was applied to the response (D3) of dosemeter element without filter for estimation of beta dose. Protocol for identification of beta in the mixed gamma beta fields was defined such that it resulted in nearly same correction factor for given beta source in both gamma beta mixed fields and pure beta fields. The beta dose evaluation algorithm has provided the beta dose estimation within the required accuracy for >90% cases obtained from national quality assurance test data from different laboratories.

Estimation of Eye Lens Dose During Brain Scans Using Gafchromic Xr-QA2 Film in Various Multidetector CT Scanners

The purpose of this study was to estimate eye lens dose during brain scans in 16-, 64-, 128- and 256-slice multidetector computed tomography (CT) scanners in helical acquisition mode and to test the feasibility of using radiochromic film as eye lens dosemeter during CT scanning. Eye lens dose measurements were performed using Gafchromic XR-QA2 film on a polystyrene head phantom designed with outer dimensions equivalent to the head size of a reference Indian man. The response accuracy of XR-QA2 film was validated by using thermoluminescence dosemeters. The eye lens dose measured using XR-QA2 film on head phantom for plain brain scanning in helical mode ranged from 43.8 to 45.8 mGy. The XR-QA2 film measured dose values were in agreement with TLD measured dose values within a maximum variation of 8.9%. The good correlation between the two data sets confirms the viability of using XR-QA2 film for eye lens dosimetry.

Overview of the Quality Assurance Programme Implemented for TLD Based Individual Monitoring in India

The effective implementation of a comprehensive quality assurance (QA) programme in any individual monitoring service plays a key role in attaining and sustaining the level of performance at par with international standards. In India, individual monitoring of more than 120 000 radiation workers is provided through 16 laboratories using a CaSO4:Dy based thermoluminescence dosimetry (TLD) badge system. In such a wide-spread programme, the harmonization of procedures and regular QA check on the dosimetry system are utmost important to ensure the uniform standard of accuracy and reliability of the service. This paper discusses some aspects of the QA programme implemented at different stages of the TLD monitoring system and provides the results of the performance test of monitoring laboratories.