Theoretical evaluation of adiabatic and vertical electron affinity of some radiosensitizers in solution using FEP, ab initio and DFT methods

The biological activity of radiosensitizers is associated to their electron affinity (EA), which can be divided in two main processes: vertical and adiabatic. In this work, we calculated the EAs of nitrofurans and nitroimidazoles (Fig. 2) using Hartree–Fock (HF) and density functional theory (DFT) methods and evaluated solvent effects (water and carbon tetrachloride) on EAs. For water, we combined the polarized continuum model (PCM) and free energy perturbation (FEP) (finite difference thermodynamic integration, FDTI) methods. For carbon tetrachloride, we used the FDTI method. The values of adiabatic EA obtained are in agreement with experimental data (deviations of 0.013 eV). The vertical EAs were calculated according to Cederbaum's outer valence Green function (OVGF) method. This methodology, which relies on theoretical aspects of free energy calculations on charged molecules in solution, was used to select potential selective radiosensitizers from recently reported compounds and could be helpful in the rational design of new and more selective bioreductive anticancer drugs.     

Identification of Potential Carbonic Anhydrase Inhibitors for Glaucoma Treatment Through an In-Silico Approach

Glaucoma is a neurodegenerative disease and second leading cause of blindness in western world. The disease is characterized by an elevated intraocular pressure. Carbonic anhydrase plays a major role by forming aqueous humor and its inhibition can reduce intraocular pressure by partially suppressing the secretion of aqueous humor. Thus in this study, we proposed to identify the potential novel compounds targeting the carbonic anhydrase. The diversity set-II molecules library consisting of 1880 compounds from National Cancer Institute were virtually screened (molecular docking) against human carbonic anhydrase protein. For the obtained best compounds, the nature of the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO), which determine nucleophilic and electrophilic activity, were calculated by using density functional theory (DFT). The in silico screening suggested 5 best compounds that are effective in comparison to the dorzolamide, a widely used carbonic anhydrase inhibitor for glaucoma treatment. Of the five compounds, 4-nitro-7-[(1-oxidopyridin-1-ium-2-yl) thio] benzofurazan (ZINC01757986) exhibited the better binding affinity (??9.2 cal/mol) in comparison to dorzolamide (??7.2 kcal/mol). The DFT studies on novel identified compound, ZINC01757986 exhibited less HOMO–LUMO energy gap, low hardness and more softness (0.2305 eV, 0.1152 eV and 8.6805 eV) when compared to dorzalamide (0.9536 eV, 0.4768 eV and 2.0973 eV). These studies emphasize that ZINC01757986 can be used as potential carbonic anhydrase inhibitor and lead compounds for the development of an effective anti-glaucoma drug. The results emphasize that these compounds could be potential lead molecules for further structure-based discovery of antiglaucoma drugs.

     

A study of hydration of sodium hyaluronate from compressibility and high precision densitometric measurements

Sodium hyaluronate samples of various molecular weights were prepared and characterized by size exclusion chromatography and dilution viscometry. Densitometry was used for determination of the densities of sodium hyaluronate solutions and the results expressed as partial specific volumes in the respective solvents. Solution adiabatic compressibilities were measured. Hydration parameters of sodium hyaluronate were consequently determined and the values obtained discussed in relation to existing hydration models for sodium hyaluronate in water.     

Quantification of thermal ring flexibilities of aromatic and heteroaromatic compounds

The consequences of thermal fluctuations occurring at room temperatures on the aromatic character of a broad group of compounds were analyzed in three distinct ways. First of all, the ring deformations were modeled along normal coordinates coming from quantum thermo-chemistry computations. The amplitudes of vibrations were estimated according to absorbed energies at room temperature. Alternatively, in-plane and out-of-plane ring deformations were modeled via scanning procedure with partial relaxation of the molecular geometry. The influence of ring deformations on π–electron delocalization was expressed in terms of HOMA values. Besides, the ring deformability was defined as the averaged change of bond angles or dihedral angles constituting the ring that was associated with 1.5 kcal mol^-1 increase of the system energy. The molecules structures adopted during vibrations at room temperature can lead to significant heterogeneity of structural index of aromaticity. The broad span of HOMA values was obtained for analyzed five- or six-membered aromatic and heteroaromatic rings. However, the averaged values obtained for such fluctuations almost perfectly match HOMA values of molecule in the ground state. It has been demonstrated that the ring deformability imposed by bond angle changes is much smaller than for dihedral angles with the same rise of system energy. Interestingly in the case of out-of-plane vibrations modeled by scanning procedure there is observed linear correlation between ring deformability and HOMA values. Proposed method for inclusion of thermal vibrations in the framework of π–electron delocalization provides natural shift of the way of thinking about aromaticity from a static quantity to a dynamic and heterogeneous one due to inclusion of a more realistic object of analysis – thermally deformed structures. From this perspective the thermal fluctuations are supposed to be non-negligible contributions to aromaticity phenomenon.     

Assessment of MIRD data for internal dosimetry using the GATE Monte Carlo code

GATE/GEANT is a Monte Carlo code dedicated to nuclear medicine that allows calculation of the dose to organs of voxel phantoms. On the other hand, MIRD is a well-developed system for estimation of the dose to human organs. In this study, results obtained from GATE/GEANT using Snyder phantom are compared to published MIRD data. For this, the mathematical Snyder phantom was discretized and converted to a digital phantom of 100 × 200 × 360 voxels. The activity was considered uniformly distributed within kidneys, liver, lungs, pancreas, spleen, and adrenals. The GATE/GEANT Monte Carlo code was used to calculate the dose to the organs of the phantom from mono-energetic photons of 10, 15, 20, 30, 50, 100, 200, 500, and 1000 keV. The dose was converted into specific absorbed fraction (SAF) and the results were compared to the corresponding published MIRD data. On average, there was a good correlation (r ²>0.99) between the two series of data. However, the GATE/GEANT data were on average −0.16 ± 6.22% lower than the corresponding MIRD data for self-absorption. Self-absorption in the lungs was considerably higher in the MIRD compared to the GATE/GEANT data, for photon energies of 10–20 keV. As for cross-irradiation to other organs, the GATE/GEANT data were on average +1.5 ± 8.1% higher than the MIRD data, for photon energies of 50–1000 keV. For photon energies of 10–30 keV, the relative difference was +7.5 ± 67%. It turned out that the agreement between the GATE/GEANT and the MIRD data depended upon absolute SAF values and photon energy. For 10–30 keV photons, where the absolute SAF values were small, the uncertainty was high and the effect of cross-section prominent, and there was no agreement between the GATE/GEANT results and the MIRD data. However, for photons of 50–1,000 keV, the bias was negligible and the agreement was acceptable.     

Modeling winter wheat phenology and carbon dioxide fluxes at the ecosystem scale based on digital photography and eddy covariance data

Recent studies have shown that the greenness index derived from digital camera imagery has high spatial and temporal resolution. These findings indicate that it can not only provide a reasonable characterization of canopy seasonal variation but also make it possible to optimize ecological models. To examine this possibility, we evaluated the application of digital camera imagery for monitoring winter wheat phenology and modeling gross primary production (GPP).By combining the data for the green cover fraction and for GPP, we first compared 2 different indices (the ratio greenness index (green-to-red ratio, G/R) and the relative greenness index (green to sum value, G%)) extracted from digital images obtained repeatedly over time and confirmed that G/R was best suited for tracking canopy status. Second, the key phenological stages were estimated using a time series of G/R values. The mean difference between the observed phenological dates and the dates determined from field data was 3.3 days in 2011 and 4 days in 2012, suggesting that digital camera imagery can provide high-quality ground phenological data.Furthermore, we attempted to use the data (greenness index and meteorological data in 2011) to optimize a light use efficiency (LUE) model and to use the optimal parameters to simulate the daily GPP in 2012. A high correlation (R² = 0.90) was found between the values of LUE-based GPP and eddy covariance (EC) tower-based GPP, showing that the greenness index and meteorological data can be used to predict the daily GPP. This finding provides a new method for interpolating GPP data and an approach to the estimation of the temporal and spatial distributions of photosynthetic productivity.In this study, we expanded the potential use of the greenness index derived from digital camera imagery by combining it with the LUE model in an analysis of well-managed cropland. The successful application of digital camera imagery will improve our knowledge of ecosystem processes at the temporal and spatial levels.     

A molecular orbital treatment on the reactivity of some polycyclic aromatic hydrocarbon diol-epoxides

The total energies of several aromatic hydrocarbon diol-epoxides and their carbonium cations were calculated by using the semi-empirical CNDO/2 SCF MO calculations. The reactivity of the diol-epoxides toward electrophilic reagents was evaluated from the energy difference between the diol-epoxide and its carbonium cation. Next, the total energies of the diol-epoxides and their carbonium cations were divided into several components, which represent the contributions of parts of molecules to the total energy. That is, the total energy was represented as the sum of the block energies and inter-block energies, which were defined in relation to the parts of the molecules. By using this concept of the block and the inter-block energies the reactivities of the diol-epoxides were interpreted in connection with the molecular structure.     

Nutrients and temperature additively increase stream microbial respiration

Rising temperatures and nutrient enrichment are co‐occurring global‐change drivers that stimulate microbial respiration of detrital carbon, but nutrient effects on the temperature dependence of respiration in aquatic ecosystems remain uncertain. We measured respiration rates associated with leaf litter, wood, and fine benthic organic matter (FBOM) across seasonal temperature gradients before (PRE) and after (ENR1, ENR2) experimental nutrient (nitrogen [N] and phosphorus [P]) additions to five forest streams. Nitrogen and phosphorus were added at different N:P ratios using increasing concentrations of N (~80–650 μg/L) and corresponding decreasing concentrations of P (~90–11 μg/L). We assessed the temperature dependence, and microbial (i.e., fungal) drivers of detrital mass‐specific respiration rates using the metabolic theory of ecology, before vs. after nutrient enrichment, and across N and P concentrations. Detrital mass‐specific respiration rates increased with temperature, exhibiting comparable activation energies (E, electronvolts [eV]) for all substrates (FBOM E = 0.43 [95% CI = 0.18–0.69] eV, leaf litter E = 0.30 [95% CI = 0.072–0.54] eV, wood E = 0.41 [95% CI = 0.18–0.64] eV) close to predicted MTE values. There was evidence that temperature‐driven increased respiration occurred via increased fungal biomass (wood) or increased fungal biomass‐specific respiration (leaf litter). Respiration rates increased under nutrient‐enriched conditions on leaves (1.32×) and wood (1.38×), but not FBOM. Respiration rates responded weakly to gradients in N or P concentrations, except for positive effects of P on wood respiration. The temperature dependence of respiration was comparable among years and across N or P concentration for all substrates. Responses of leaf litter and wood respiration to temperature and the combined effects of N and P were similar in magnitude. Our data suggest that the temperature dependence of stream microbial respiration is unchanged by nutrient enrichment, and that increased temperature and N + P availability have additive and comparable effects on microbial respiration rates.     

Recognition of cyclooxygenase-2 (COX-2) active site by NSAIDs: a computer modelling study

The energetics and models of COX-2 complexed with nonsteroidal anti-inflammatory drugs (NSAIDs) having different degrees of selectivity for two isoforms of COX (COX-2 and COX-1) have been studied using computer modelling approach. The models are obtained for complexes of NS398 (NS), a selective COX-2 inhibitor; indoprofen (Ind), a non-selective inhibitor; di-tert-butylbenzofurans (DHDMBFs) with substituents at the 5th position: CONH(CH2)2OMe (BF1), CONH-c-Pr (BF2), 3-methylene-gamma-butyrolactonyl (BF3) and oxicams namely, meloxicam (Mel), piroxicam (Pir) and tenoxicam (Ten). These were optimized using molecular mechanics (MM) and molecular dynamics (MD) techniques. The binding energies and structures were compared with pharmacological parameters and available results with COX-1. In case of NS a larger difference in the binding energies between COX-2 and COX-1 was noticed as compared to that of Ind. It also had stronger interaction with His90 and Tyr355 which is considered important for COX-2 selectivity. There was a difference in the compactness at the channel entrance between COX-2 selective and non-selective ligands. Models with DHDMBFs and oxicams showed a similar correlation. The results were used to design a peptide inhibitor, Tyr-Arg-Cys-Ala-delta Phe-Cys (Pept) which could fit better in the COX-2 cavity. As per our MD simulation results this peptide inhibitor showed both higher activity and COX-2 selectivity.     

A practical tool to evaluate dose distributions using radiochromic film in radiation oncology

PurposeTriple channel algorithm and specific procedures make more reliable radiochromic dosimetry for treatment planning verification and quality assurance in radiation therapy. A tool to obtain radiochromic dose distributions and compare them with the ones resulting from a treatment planning system was developed and applied.MethodsThe tool was developed as Microsoft Excel macro; it builds dose calibration curves against net optical density of Gafchromic EBT3 film, produces axial, coronal and sagittal dose maps and allows to evaluate them against dose distributions calculated by the Varian treatment planning system Eclipse using gamma index and gamma angle.ResultsThe net optical density standard errors of estimate of calibration curves at 6 MV Varian DBX600 linac energy were 0.2%, 0.4% and 0.2% for the red, green and blue channels. Tests of these curves by means of three independent eight dose points measurement series, at 15 MV and 6 MV Varian 2100C linac and at 6 MV DBX600 linac energies, showed less than 2% of dose errors for the red channel and less than 3% for the green channel in the range 100–450 cGy. The comparisons between dose distributions from Gafchromic EBT3 triple channel algorithm and the ones from Eclipse analytic anisotropic algorithm (AAA) showed values of gamma index 95th percentile between 0.6 and 1.0.ConclusionThe obtained results encourage the application of this tool in radiation therapy quality assurance.     

Global synthesis of the temperature sensitivity of leaf litter breakdown in streams and rivers

Streams and rivers are important conduits of terrestrially derived carbon (C) to atmospheric and marine reservoirs. Leaf litter breakdown rates are expected to increase as water temperatures rise in response to climate change. The magnitude of increase in breakdown rates is uncertain, given differences in litter quality and microbial and detritivore community responses to temperature, factors that can influence the apparent temperature sensitivity of breakdown and the relative proportion of C lost to the atmosphere vs. stored or transported downstream. Here, we synthesized 1025 records of litter breakdown in streams and rivers to quantify its temperature sensitivity, as measured by the activation energy (E_a, in eV). Temperature sensitivity of litter breakdown varied among twelve plant genera for which E_a could be calculated. Higher values of E_a were correlated with lower‐quality litter, but these correlations were influenced by a single, N‐fixing genus (Alnus). E_a values converged when genera were classified into three breakdown rate categories, potentially due to continual water availability in streams and rivers modulating the influence of leaf chemistry on breakdown. Across all data representing 85 plant genera, the E_a was 0.34 ± 0.04 eV, or approximately half the value (0.65 eV) predicted by metabolic theory. Our results indicate that average breakdown rates may increase by 5–21% with a 1–4 °C rise in water temperature, rather than a 10–45% increase expected, according to metabolic theory. Differential warming of tropical and temperate biomes could result in a similar proportional increase in breakdown rates, despite variation in E_a values for these regions (0.75 ± 0.13 eV and 0.27 ± 0.05 eV, respectively). The relative proportions of gaseous C loss and organic matter transport downstream should not change with rising temperature given that E_a values for breakdown mediated by microbes alone and microbes plus detritivores were similar at the global scale.     

DFT studies of hydrocarbon combustion on metal surfaces

Catalytic combustion of hydrocarbons is an important technology to produce energy. Compared to conventional flame combustion, the catalyst enables this process to operate at lower temperatures; hence, reducing the energy required for efficient combustion. The reaction and activation energies of direct combustion of hydrocarbons (CH?→?C?+?H) on a series of metal surfaces were investigated using density functional theory (DFT). The data obtained for the Ag, Au, Al, Cu, Rh, Pt, and Pd surfaces were used to investigate the validity of the Brønsted-Evans-Polanyi (BEP) and transition state scaling (TSS) relations for this reaction on these surfaces. These relations were found to be valid (R²?=?0.94 for the BEP correlation and R²?=?1.0 for the TSS correlation) and were therefore used to estimate the energetics of the combustion reaction on Ni, Co, and Fe surfaces. It was found that the estimated transition state and activation energies (E_TS?=??69.70 eV and E_a?=?1.20 eV for Ni, E_TS?=??87.93 eV and E_a?=?1.08 eV for Co and E_TS?=??92.45 eV and E_a?=?0.83 eV for Fe) are in agreement with those obtained by DFT calculations (E_TS?=??69.98 eV and E_a?=?1.23 eV for Ni, E_TS?=??87.88 eV and E_a?=?1.08 eV for Co and E_TS?=??92.57 eV and E_a?=?0.79 eV for Fe). Therefore, these relations can be used to predict energetics of this reaction on these surfaces without doing the time consuming transition state calculations. Also, the calculations show that the activation barrier for CH dissociation decreases in the order Ag ? Au ? Al ? Cu ? Pt ? Pd ? Ni?>?Co?>?Rh?>?Fe.     

Investigations of the kinetic energy parameters of irradiated (La)‐doped phosphate glass

The current study characterizes and analyzes glow curves obtained from phosphate glass doped with different concentrations of lanthanum. Kinetic parameters of the glow curves obtained from beta‐irradiated phosphate glass samples doped with lanthanum were determined using a newly designed deconvoluted software. The obtained results from the analyses indicated that the glow curves of the phosphate glass samples were composed of five overlapping peaks. The activation energies of the five electron traps were located between 0.622 and 1.133 eV. The obtained kinetic parameters were evaluated using the designed software and another two methods and all revealed good agreement. The first three traps displayed non‐first‐order behaviour, while the two deep traps obeyed nearly first‐order kinetics.     

On the electronic properties of two-dimensional honeycomb GaInN and GaAlN alloys: a molecular analysis

We have performed first principles total energy calculations to investigate the structural and the electronic properties of two-dimensional honeycomb GaAlN and GaInN alloys. Calculations were done using a coronene-like (C₂₄H₁₂) cluster and for different numbers of Ga, Al, and In atoms. The exchange and correlation potential energies were treated within the generalized gradient approximation (GGA). The bond length, dipole moment, binding energy, and gap between the HOMO and the LUMO are reported as a function of x. The stability of the structures depends on the site of the substituted atom; for example, when three Ga atoms are substituted, the GaInN alloy becomes unstable. The gap in the GaAlN increases from 3.76 eV (GaN) to 4.51 eV (AlN), and in the GaInN decreases to 2.11 eV. The biggest polarity occurs when eight and four Ga atoms are substituted, for GaAlN and GaInN, respectively.     

A comparison of the CHARMM, AMBER and ECEPP potentials for peptides. II. Phi-psi maps for N-acetyl alanine N'-methyl amide: comparisons, contrasts and simple experimental tests

phi-psi maps of N-acetyl alanine N'-methyl amide have been computed using the CHARMM potential, the all-atom AMBER potential, and the ECEPP/2 potential, before and after adiabatic relaxation. Maps using the CHARMM and AMBER potentials were determined with values of 1.0 and 4.0 for the dielectric constant epsilon, and with a distance dependent dielectric constant. Adiabatic relaxation was carried out using flexible geometry for the CHARMM and AMBER potentials, and using rigid geometry for the AMBER and ECEPP potentials. In all cases, the lowest energy was found in the C7eq region (phi approximately -70 degrees, psi approximately 70 degrees). The maps with CHARMM and AMBER with epsilon = 4.0 and with ECEPP, without adiabatic relaxation, were broadly similar but differed in the relative energies allotted to high-energy regions of the map. After adiabatic relaxation with rigid geometry, the map with ECEPP, and the map with AMBER using a distance-dependent dielectric constant, agreed fairly well apart from differences in the relative energies of the alpha R, alpha L, and C7ax regions. After adiabatic relaxation with flexible geometry, the maps with CHARMM and AMBER became very similar; the lowest energies were observed in the C7eq region, the C5 region (phi approximately -150 degrees, psi approximately 150 degrees) and the C7ax region (phi approximately 70 degrees, psi approximately -70 degrees). Breakdown of the energies, after adiabatic relaxation, into electrostatic, nonbonded, and geometric (including torsional) contributions, showed that (1) with fixed geometry, the nonbonded and torsional contribution to the ECEPP and AMBER potentials were very similar, but the electrostatic contributions were markedly different; (2) with flexible geometry, the nonbonded contribution to the CHARMM and AMBER potentials did not vary greatly over the whole map. The phi-psi maps were subjected to three simple comparisons with experiment. (1) The maps were used to predict the characteristic ratio for poly-L-alanine, and the results were compared with experimental findings (D.A. Brant and P.J. Flory, J. Amer. Chem. Soc. 87, 2788-2791, 1965). The agreement with experiment was acceptable for ECEPP, and for CHARMM after adiabatic relaxation, marginal for AMBER after adiabatic relaxation, and unsatisfactory for CHARMM or AMBER without adiabatic relaxation. (2) Deviations of bond angles from their equilibrium values, in energy-minimized conformations, were compared with values deduced from crystals of terminally-blocked amino acids. With both the CHARMM and AMBER potentials using flexible geometry, one or more excessive deviations was observed in the C7ax local minimum.(ABSTRACT TRUNCATED AT 400 WORDS)     

A computational study of the interaction between dopamine and DNA/RNA nucleosides

The interaction between protonated dopamine and neutral RNA and DNA nucleosides was studied by means of density functional theory calculations in vacuum and in implicit water. On the most stable complexes formed with each of the nucleosides, the vertical absorption excitation energies were evaluated and compared with the values of separated dopamine and corresponding nucleoside. The most stable complex was formed with guanosine and the spectral changes in this complex resulted in a significant reduction of the oscillator strength of the first dopamine’s transition. In the first guanosine’s transition, a redshift of 0.2 eV was found combined with a reduction of the oscillator strength.     

Improved reproducibility in the determination of the molecular weight of chitosan by analytical size exclusion chromatography

The reproducibility of the determination of the molecular weight of chitosans in the 90–210 kDa range (M_n) by analytical size exclusion chromatography with multi-angle laser light scattering (SEC-MALLS) was improved by reducing the salt concentration in the mobile phase from (0.3 M acetic acid, 0.2 M sodium acetate, and 0.8 mM sodium azide) to (0.15 M acetic acid, 0.1 M sodium acetate, and 0.4 mM sodium azide) using Tosoh TSKgel G6000PW_XL and G5000PW_XL columns in series. The variability of measured molecular weight was significantly reduced by lowering the acetate concentration in the mobile phase, while the average molecular weight did not change significantly. The coefficient of variation of the number-average molecular weight, CV(M_n), decreased from 7–12% to 3–6% upon mobile phase dilution. This reduced variability in molecular weight of chitosans obtained from SEC is a significant improvement when precise values of chitosan molecular weight are required, for example in stability studies where viscosity changes in concentrated chitosan solutions are assessed, and in gene delivery applications.     

Prokaryotic gene finding based on physicochemical characteristics of codons calculated from molecular dynamics simulations

An ab initio model for gene prediction in prokaryotic genomes is proposed based on physicochemical characteristics of codons calculated from molecular dynamics (MD) simulations. The model requires a specification of three calculated quantities for each codon: the double-helical trinucleotide base pairing energy, the base pair stacking energy, and an index of the propensity of a codon for protein-nucleic acid interactions. The base pairing and stacking energies for each codon are obtained from recently reported MD simulations on all unique tetranucleotide steps, and the third parameter is assigned based on the conjugate rule previously proposed to account for the wobble hypothesis with respect to degeneracies in the genetic code. The third interaction propensity parameter values correlate well with ab initio MD calculated solvation energies and flexibility of codon sequences as well as codon usage in genes and amino acid composition frequencies in ∼175,000 protein sequences in the Swissprot database. Assignment of these three parameters for each codon enables the calculation of the magnitude and orientation of a cumulative three-dimensional vector for a DNA sequence of any length in each of the six genomic reading frames. Analysis of 372 genomes comprising ∼350,000 genes shows that the orientations of the gene and nongene vectors are well differentiated and make a clear distinction feasible between genic and nongenic sequences at a level equivalent to or better than currently available knowledge-based models trained on the basis of empirical data, presenting a strong support for the possibility of a unique and useful physicochemical characterization of DNA sequences from codons to genomes.     

Thermoluminescence studies on photosynthetic energy conversion. II. Activation energies for three energy storage states associated with Photoreaction II of higher plants

Thermoluminescent glow curves of isolated chloroplasts were analyzed to determine activation energy, frequency factor and lifetime of the three glow peaks associated with Photoreaction II. Peak 1 had an activation energy of 0.8 eV, a frequency factor of 1.5 · 10¹⁴/s and a lifetime of 0.02 s. This glow peak appeared to be rather different from peaks 2 and 3, which were quite similar to each other. Peaks 2 and 3 gave activation energies respectively of 0.48 eV and 0.57 eV, frequency factors of 7 · 10⁷/s and 1 · 10⁹/s, and lifetimes of about 1.5 s. A second, less reliable method of analysis gave activation energies of 0.72 eV for peak 1 and 0.44 eV for peak 2. Our values are compared with those obtained by other workers, and the possible metastable states reflected by the three glow peaks are discussed.