Capture and immobilization of gaseous ruthenium tetroxide RuO4 in the UiO-66-NH2 metal-organic framework.

106Ru is a radioactive isotope usually generated by the nuclear industry within power plant reactors. During a nuclear accident, 106Ru reacts with oxygen, leading to the production of highly volatile ruthenium tetroxide RuO4. The combination of volatility and radioactivity makes 106RuO4, one of the most radiotoxic species and justifies the development of a specific setup for its capture and immobilization. In this study, we report for the first time the capture and immobilization of gaseous RuO4 within a porous metal-organic framework (UiO-66-NH2). We used specific installation for the production of gaseous RuO4 as well as for the quantification of this gas trapped within the filtering medium. We proved that UiO-66-NH2 has remarkable affinity for RuO4 capture, as this MOF exhibited the worldwide highest RuO4 decontamination factor (DF of 5745), hundreds of times higher than the DF values of sorbents daily used by the nuclear industry (zeolites or activated charcoal). The efficiency of UiO-66-NH2 can be explained by its pore diameters well adapted to the capture and immobilization of RuO4 as well as its conversion into stable RuO2 within the pores. This conversion corresponds to the reactivity of RuO4 with the MOF organic sub-network, leading to the oxidation of terephthalate ligands. As proved by powder X-ray diffraction and NMR techniques, these modifications did not decompose the MOF structure.

Capture of Gaseous Iodine in Isoreticular Zirconium-Based UiO-n Metal-Organic Frameworks: Influence of Amino Functionalization, DFT Calculations, Raman and EPR Spectroscopic Investigation.

A series of Zr-based UiO-n MOF materials (n=66, 67, 68) have been studied for iodine capture. Gaseous iodine adsorption was collected kinetically from a home-made set-up allowing the continuous measurement of iodine content trapped within UiO-n compounds, with organic functionalities (-H, -CH3 , -Cl, -Br, -(OH)2 , -NO2 , -NH2 , (-NH2 )2 , -CH2 NH2 ) by in-situ UV-Vis spectroscopy. This study emphasizes the role of the amino groups attached to the aromatic rings of the ligands connecting the {Zr6 O4 (OH)4 } brick. In particular, the preferential interaction of iodine with lone-pair groups, such as amino functions, has been experimentally observed and is also based on DFT calculations. Indeed, higher iodine contents were systematically measured for amino-functionalized UiO-66 or UiO-67, compared to the pristine material (up to 1211 mg/g for UiO-67-(NH2 )2 ). However, DFT calculations revealed the highest computed interaction energies for alkylamine groups (-CH2 NH2 ) in UiO-67 (-128.5 kJ/mol for the octahedral cavity), and pointed out the influence of this specific functionality compared with that of an aromatic amine. The encapsulation of iodine within the pore system of UiO-n materials and their amino-derivatives has been analyzed by UV-Vis and Raman spectroscopy. We showed that a systematic conversion of molecular iodine (I2 ) species into anionic I- ones, stabilized as I- ⋅⋅⋅I2 or I3 - complexes within the MOF cavities, occurs when I2 @UiO-n samples are left in ambient light.

Extrusion-Spheronization of UiO-66 and UiO-66_NH2 into Robust-Shaped Solids and Their Use for Gaseous Molecular Iodine, Xenon, and Krypton Adsorption.

The use of an extrusion-spheronization process was investigated to prepare robust and highly porous extrudates and granules starting from UiO-66 and UiO-66_NH2 metal-organic framework powders. As-produced materials were applied to the capture of gaseous iodine and the adsorption of xenon and krypton. In this study, biosourced chitosan and hydroxyethyl cellulose (HEC) are used as binders, added in low amounts (less than 5 wt % of the dried solids), as well as a colloidal silica as a co-binder when required. Characterizations of the final shaped materials reveal that most physicochemical properties are retained, except the textural properties, which are impacted by the process and the proportion of binders (BET surface area reduction from 5 to 33%). On the other hand, the mechanical resistance of the shaped materials toward compression is greatly improved by the presence of binders and their respective contents, from 0.5 N for binderless UiO-66 granules to 17 N for UiO-66@HEC granules. UiO-66_NH2-based granules demonstrated consequent iodine capture after 48 h, up to 527 mg/g, in line with the pristine UiO-66_NH2 powder (565 mg/g) and proportionally to the retaining BET surface area (-5% after shaping). Analogously, the shaped materials presented xenon and krypton sorption isotherms correlated to their BET surface area and high predicted xenon/krypton selectivity, from 7.1 to 9.0. Therefore, binder-aided extrusion-spheronization is an adapted method to produce shaped solids with adequate mechanical resistance and retained functional properties.

Interaction process between gaseous CH3I and NaCl particles: implication for iodine dispersion in the atmosphere.

Gaseous iodomethane (CH3I) is naturally emitted into the atmosphere by biological activity in oceans and during severe accidents (SAs) in nuclear power plants. In this latter case, a part of radioactive iodine such as 131I may be released. Improving the knowledge of CH3I transport and reactivity in the atmosphere is important since they are strongly linked to first the cycle of ozone and second to the dispersion of radioactive CH3I with potential radiological consequences on both the environment and human health. Here, the interaction process of CH3I with NaCl as a surrogate of atmospheric aerosols was investigated under ambient air conditions by using Diffuse Reflectance Fourier Transform Spectroscopy (DRIFTS). The DRIFTS spectra of NaCl clearly evidenced CH3I adsorption on the NaCl particle surface. A part of CH3I ((1.68 ± 0.85) × 1014 molecule per mgNaCl) was found to be strongly bonded to NaCl since no desorption was observed. The CH3I adsorption on the NaCl surface presented a 1st order kinetics relative to its gas phase concentration. The uptake coefficient was determined to be in the order of 10-11. These results show a low probability of CH3I to be taken up by halide-containing aerosols. These data are crucial for completing the iodine atmospheric chemical scheme.

Stability and radioactive gaseous iodine-131 retention capacity of binderless UiO-66-NH2 granules under severe nuclear accidental conditions.

In the present work, we aim to investigate the ability of the zirconium-based MOF-type compound UiO-66-NH2, to immobilize molecular gaseous iodine under conditions analogous to those encountered in an operating Filtered Containment Venting System (FCVS) line. Typically, the UiO-66-NH2 particles were exposed to 131I (beta and gamma emitters) and submitted to air/steam at 120 °C, under gamma irradiation (1.9 kGy h-1). In parallel to this experiment under simulated accidental conditions, the stability of the binderless UiO-66-NH2 granules under steam and gamma irradiation was investigated. In order to fit with the specifications required by typical venting systems, and to compare the efficiency of the selected MOF to porous materials commonly used by the industry, scale-up syntheses and UiO-66-NH2 millimetric-size shaping were realized. For this task, we developed an original binderless method, in order to analyze solely the efficiency of the UiO-66-NH2 material. The shaped MOF particles were then submitted separately to gamma irradiation, steam and temperature, for confirming their viability in a venting process. Their structural, textural and mechanical behaviors were characterized by the means several techniques including gas sorption, powder X-ray diffraction, infrared spectroscopy and crushing tests. Promising results were obtained to trap gaseous molecular iodine in severe accidental conditions.

Influence of Light and Temperature on the Extractability of Cerium(IV) as a Surrogate of Plutonium(IV) and its Effect on the Simulation of an Accidental Fire in the PUREX Process.

Modeling of plutonium(IV) behavior during an accidental fire in a reprocessing plant was considered using various non-radioactive metallic surrogates. Among those elements, cerium(IV) was supposed to be a suitable candidate due to possible formation of a complex with TPB, but its extractability and stability have not been studied previously under representative plutonium uranium reduction extraction (PUREX) conditions. In this work, we investigated the chemical analogy between cerium(IV) and plutonium(IV) in this extractive process and combustion thereof. Distribution ratios are reported for acidities of 1-4 mol L-1 in equal volumes of nitric acid and a 30:70 mixture of tributylphosphate and hydrogenated tetrapropylene. The influences of light, temperature, and extraction time were studied by UV-vis spectroscopy. The results showed that cerium(IV) is extracted quantitatively but is reduced over time to cerium(III) in the organic mixture. Spectrophotometric investigations of this reaction kinetics revealed an apparent rate constant k of 0.021 ± 0.002 mol0.5 L0.5 min-1 at 298 K and an apparent fractional reaction order of 0.5. The activation energy of this reduction was found to be around 82 ± 2 kJ mol-1 by the Arrhenius plot method. The combustion of mono- and biphasic solutions prepared with a cerium(IV) concentration of 10 g L-1 revealed that the extracted complexes, Ce2O·6NO3·3TBP(org) or Ce4O4·8NO3·6TBP(org), are reduced during the combustion. Compositions of the resulting ashes and soot were analyzed and highlighted the presence of pyrophosphates and polycyclic aromatic hydrocarbons, with some traces of cerium. Ce(IV) is not suitable to represent Pu(IV) from a chemical point of view in HNO3/TBP-HTP solutions.

Microhydration of caesium metaborate: structural and thermochemical properties of CsBO2 + n H2O (n = 1-4) aggregates.

The structures and thermodynamic properties of microhydrates of caesium metaborate (CsBO2) of nuclear safety interest are reported in this work. CsBO2 + n H2O (n = 1-4) molecular complexes were identified on the potential energy surface. The structures were optimized using the ωB97XD DFT method and the aug-cc-pVTZ basis set. Single-point energies were calculated at the CCSD(T)-F12a/awCVTZ and the ωB97XD/aug-cc-pVQZ levels of theory. The standard reaction enthalpies and the standard Gibbs free reaction energies were reported for all molecular complexes. The temperature dependence of ΔrG°(T) was evaluated for all studied structures over the temperature range 300-2000 K. Total hydration reactions were investigated. The results showed that the mono-hydrated form of CsBO2 exists only at temperatures lower than 720 K under standard conditions. The influence on the thermodynamic properties of the number of water molecules in the clusters was described, with successive dehydration from 720 to 480 K. In nuclear severe accident conditions, gaseous CsBO2 will remain unhydrated in the reactor coolant system.

A DFT study of RuO4 interactions with porous materials: metal-organic frameworks (MOFs) and zeolites.

Radioactive gaseous ruthenium tetroxide (RuO4) can be released into the environment in the case of a severe nuclear accident. Using periodic dispersion corrected density functional theory calculations, we have investigated for the first time the adsorption behavior of RuO4 into prototypical porous materials, Metal-Organic Frameworks (MOFs) and zeolites, with the aim of mitigating ruthenium releases to the outside. For the MOFs, we have screened a set of six structures (MIL-53(Al), MIL-120(Al), HKUST-1(Cu), UiO-66(Zr), UiO-67(Zr) and UiO-68(Zr)), while for the zeolites two structures have been selected: mordenite (MOR) with Si/Al ratios of 11 and 5, and faujasite (FAU) with a Si/Al ratio of 2.4. The DFT calculations show that the nature of the porous materials does not have a significant effect on the adsorption energy of RuO4 compounds and that the main interaction is due to the formation of hydrogen bonds. For the tested materials, computational results show that the interaction energies of RuO4 reach their maximum with the hydrated form of HKUST-1(Cu) (-114 kJ mol-1) due to the presence of strong hydrogen bonds between the water molecules and the oxygen atoms of RuO4.

Dissociative iodomethane adsorption on Ag-MOR and the formation of AgI clusters: an ab initio molecular dynamics study.

Radioactive iodine species belong to the most dangerous components of nuclear effluents and waste produced by nuclear facilities. In this work, we use computer simulations at the periodic DFT level to investigate dissociative adsorption of iodomethane on silver-exchanged mordenite, which is among the most effective sorbents of iodine species available today. The structure, energetics, and mobility of complexes Ag-(CH3I) and Ag-(CH3I)2 formed upon adsorption of iodomethane on Ag+ sites are investigated using the ab initio MD approach. The free-energy profiles for the reaction CH3I + Ag-MOR → AgI + CH3-MOR are determined using the blue moon ensemble technique. The AgI species formed as a product of dissociative adsorption are shown to combine spontaneously into small clusters (AgI)n with the dimensions restricted by the size and geometry of confining voids. The structure and energetics of the (AgI)n species are analysed in detail and compared with the available experimental and theoretical data. The internal energy of formation of clusters in mordenite is shown to contribute significantly to the shift of equilibrium from the undissociated to dissociated form of adsorbed CH3I.

Facing the challenge of predicting the standard formation enthalpies of n-butyl-phosphate species with ab initio methods.

Tributyl-phosphate (TBP), a ligand used in the PUREX liquid-liquid separation process of spent nuclear fuel, can form an explosive mixture in contact with nitric acid that might lead to a violent explosive thermal runaway. In the context of safety of a nuclear reprocessing plant facility, it is crucial to predict the stability of TBP at elevated temperatures. So far, only the enthalpies of formation of TBP are available in the literature with rather large uncertainties, while those of its degradation products, di-(HDBP) and mono-(H2MBP), are unknown. In this goal, we have used state-of-the art quantum chemical methods to compute the formation enthalpies and entropies of TBP and its degradation products di-(HDBP) and mono-(H2MBP) in gas and liquid phases. Comparisons of levels of quantum chemical theory revealed that there are significant effects of correlation on their electronic structures, pushing for the need of not only high level of electronic correlation treatment, namely, local coupled cluster with single and double excitation operators and perturbative treatment of triple excitations, but also extrapolations to the complete basis to produce reliable and accurate thermodynamics data. Solvation enthalpies were computed with the conductor-like screening model for real solvents [COSMO-RS], for which we observe errors not exceeding 22 kJ mol-1. We thus propose with final uncertainty of about 20 kJ mol-1 standard enthalpies of formation of TBP, HDBP, and H2MBP which amounts to -1281.7 ± 24.4, -1229.4 ± 19.6, and -1176.7 ± 14.8 kJ mol-1, respectively, in the gas phase. In the liquid phase, the predicted values are -1367.3 ± 24.4, -1348.7 ± 19.6, and -1323.8± 14.8 kJ mol-1, to which we may add about -22 kJ mol-1 error from the COSMO-RS solvent model. From these data, the complete hydrolysis of TBP is predicted as an exothermic phenomena but showing a slightly endergonic process.

Performance of CuII -, PbII -, and HgII -Exchanged Mordenite in the Adsorption of I2 , ICH3 , H2 O, CO, ClCH3 , and Cl2 : A Density Functional Study.

Periodic dispersion-corrected DFT is used to investigate the adsorption of I2 and ICH3 , which may be released during a severe nuclear accident, for three divalent cation (Cu2+ , Pb2+ and Hg2+ )-exchanged mordenites with an Si/Al ratio of 23. Gases such as H2 O, CO, ClCH3 , and Cl2 present in the containment atmosphere can inhibit the selective adsorption of iodine species. To identify the most promising adsorbents, a systematic study is performed in which all the possible cationic sites in the main channel of the mordenite structure are considered. For the energetically most stable sites, the divalent cation is located in the small rings (five- or six-membered) containing two Al atoms, while in the energetically less stable configurations, the two Al atoms are far apart (>7 Å) and the cation is close to only one Al atom. Upon adsorption of the various molecules, the coordination number of the cation decreases with increasing interaction energy, as the molecules can attract the divalent cations from the framework. Finally, the computed interaction energies show that Hg-mordenite (MOR) could be a suitable material for selective adsorption of volatile iodine species, contrary to Cu-MOR and Pb-MOR.

Impact of the Si/Al ratio on the selective capture of iodine compounds in silver-mordenite: a periodic DFT study.

Silver modified zeolites with a mordenite structure can capture volatile iodine compounds (I2 and ICH3) which can be released during a severe nuclear accident. However under these particular conditions, molecules such as CO and H2O present in the containment atmosphere are expected to inhibit the adsorption of iodine compounds. In the present work, periodic density functional theory calculations have been carried out to investigate the interaction of I2, ICH3, H2O and CO molecules in silver-exchanged mordenite with various Si/Al ratios with the aim of finding values that favor a selective adsorption of I2 and ICH3. Computational results show that the interaction energies of CO and H2O remain of the same order of magnitude (from -120 to -140 kJ mol-1 for CO and from -90 to -120 kJ mol-1 for H2O) for all the investigated Si/Al ratios. In contrast, ICH3 is increasingly strongly adsorbed as the Si/Al ratio decreases, from around -145 kJ mol-1 when Si/Al = 47 to -190 kJ mol-1 for Si/Al = 5. The same trend is observed for I2 with a larger amplitude: from -135 kJ mol-1 for Si/Al = 47 to -300 kJ mol-1 for Si/Al = 5. Therefore, the use of silver-exchanged mordenite with Si/Al ratios of 5 or 11 would drastically limit the inhibiting effect of contaminants on the adsorption of volatile iodine species. Also for the same ratios, a spontaneous dissociation of I2 during its adsorption is observed, leading to the formation of AgI complexes which are prerequisite for the immobilization of iodine in the long term.

A Density Functional Theory and ab Initio Investigation of the Oxidation Reaction of CO by IO Radicals.

To get an insight into the possible reactivity between iodine oxides and CO, a first step was to study the thermochemical properties and kinetic parameters of the reaction between IO and CO using theoretical chemistry tools. All stationary points involved were optimized using the Becke's three-parameter hybrid exchange functional coupled with the Lee-Yang-Parr nonlocal correlation functional (B3LYP) and the Møller-Plesset second-order perturbation theory (MP2). Single-point energy calculations were performed using the coupled cluster theory with the iterative inclusion of singles and doubles and the perturbative estimation for triple excitations (CCSD(T)) and the aug-cc-pVnZ (n = T, Q, and 5) basis sets on geometries previously optimized at the aug-cc-pVTZ level. The energetics was then recalculated using the one-component DK-CCSD(T) approach with the relativistic ANO basis sets. The spin-orbit coupling for the iodine containing species was calculated a posteriori using the restricted active space state interaction method in conjunction with the multiconfigurational perturbation theory (CASPT2/RASSI) employing the complete active space (CASSCF) wave function as the reference. The CCSD(T) energies were also corrected for BSSE for molecular complexes and refined with the extrapolation to CBS limit while the DK-CCSD(T) values were refined with the extrapolation to FCI. The exploration of the potential energy surface revealed a two-steps mechanism with a trans and a cis pathway. The rate constants for the direct and complex mechanism were computed as a function of temperature (250-2500 K) using the canonical transition state theory. The three-parameter Arrhenius expressions obtained for the direct and indirect mechanism at the DK-CCSD(T)-cf level of theory is 1.49 × 10(-17) × T(1.77) exp(-47.4 (kJ mol(-1))/RT).

Thermochemistry of Ruthenium Oxyhydroxide Species and Their Impact on Volatile Speciations in Severe Nuclear Accident Conditions.

Literature thermodynamic data of ruthenium oxyhydroxides reveal large uncertainties in some of the standard enthalpies of formation, motivating the use of high-level relativistic correlated quantum chemical methods to reduce the level of discrepancies. Reaction energies leading to the formation of all possible oxyhydroxide species RuOx(OH)y(H2O)z have been calculated for a series of reactions combining DFT (TPSSh-5%HF) geometries and partition functions, CCSD(T) energies extrapolated to the complete basis set limits. The highly accurate ab initio thermodynamic data were used as input data of thermodynamic equilibrium computations to derive the speciation of gaseous ruthenium species in the temperature, pressure and concentration conditions of severe nuclear accidents occurring in pressurized water reactors. At temperatures lower than 1000 K, gaseous ruthenium tetraoxide is the dominating species, between 1000 and 2000 K ruthenium trioxide becomes preponderant, whereas at higher temperatures gaseous ruthenium oxide, dioxide and even Ru in gaseous phase are formed. Although earlier studies predicted the formation of oxyhydroxides in significant quantities, the use of highly accurate ab initio thermodynamic data for ruthenium gaseous species leads to a more reliable inventory of gaseous ruthenium species in which gaseous oxyhydroxide ruthenium molecules are formed only in negligible amounts.

Gas-Phase Reactivity of Cesium-Containing Species by Quantum Chemistry.

Thermodynamics and kinetics of cesium species reactions have been studied by using high-level quantum chemical tools. A systematic theoretical study has been done to find suitable methodology for calculation of reliable thermodynamic properties, allowing us to determine bimolecular rate constants with appropriate kinetic theories of gas-phase reactions. Four different reactions have been studied in this work: CsO + H2 = CsOH + H (R1), Cs + HI = CsI + H (R2), CsI + H2O = CsOH + HI (R3), and CsI + OH = CsOH + I (R4). All reactions involve steam, hydrogen, and iodine in addition of cesium. Most of the reactions are fast and (R3) and (R4) proceed even without energetic barrier. In terms of chemical reactivity in the reactor coolant system (RCS) in the case of severe accident, it can be expected that there will be no kinetic limitations for main cesium species (CsOH and CsI) transported along the RCS. Cs chemical speciation inside the RCS should be governed by the thermodynamics.

Thermodynamic properties of gaseous ruthenium species.

The review of thermodynamic data of ruthenium oxides reveals large uncertainties in some of the standard enthalpies of formation, motivating the use of high-level relativistic correlated quantum chemical methods to reduce the level of discrepancies. The reaction energies leading to the formation of ruthenium oxides RuO, RuO2, RuO3, and RuO4 have been calculated for a series of reactions. The combination of different quantum chemical methods has been investigated [DFT, CASSCF, MRCI, CASPT2, CCSD(T)] in order to predict the geometrical parameters, the energetics including electronic correlation and spin-orbit coupling. The most suitable method for ruthenium compounds is the use of TPSSh-5%HF for geometry optimization, followed by CCSD(T) with complete basis set (CBS) extrapolations for the calculation of the total electronic energies. SO-CASSCF seems to be accurate enough to estimate spin-orbit coupling contributions to the ground-state electronic energies. This methodology yields very accurate standard enthalpies of formations of all species, which are either in excellent agreement with the most reliable experimental data or provide an improved estimate for the others. These new data will be implemented in the thermodynamical databases that are used by the ASTEC code (accident source term evaluation code) to build models of ruthenium chemistry behavior in severe nuclear accident conditions. The paper also discusses the nature of the chemical bonds both from molecular orbital and topological view points.

Theoretical study of plutonium(IV) complexes formed within the PUREX process: a proposal of a plutonium surrogate in fire conditions.

We present a relativistic quantum chemical study to determine the best surrogate for plutonium(IV) to be used in experimental investigations of the behavior of plutonium-nitrate-TBP in fire conditions that might occur in the nuclear fuel refining process known as PUREX. In this study geometries and stabilities of Pu(NO3)6(2-) and Pu(NO3)4(TBP)2 complexes were compared to that of equivalent complexes of selected elements from the lanthanide and actinide series (Ce, Th, U) chosen on the basis of similar ionic radii and stability as tetravalent species. PBE and PBE0 DFT functionals have proven to be sufficient and affordable for qualitative studies, performing as good as the wave function based correlated method MP2. On the basis of our results, cerium(IV) appears to be a good surrogate for plutonium(IV).

Microhydration of caesium compounds: Cs, CsOH, CsI and Cs2I2 complexes with one to three H2O molecules of nuclear safety interest.

Structure and thermodynamic properties (standard enthalpies of formation and Gibbs free energies) of hydrated caesium species of nuclear safety interest, Cs, CsOH, CsI and its dimer Cs2I2, with one up to three water molecules, are calculated to assess their possible existence in severe accident occurring to a pressurized water reactor. The calculations were performed using the coupled cluster theory including single, double and non-iterative triple substitutions (CCSD(T)) in conjunction with the basis sets (ANO-RCC) developed for scalar relativistic calculations. The second-order spin-free Douglas-Kroll-Hess Hamiltonian was used to account for the scalar relativistic effects. Thermodynamic properties obtained by these correlated ab initio calculations (entropies and thermal capacities at constant pressure as a function of temperature) are used in nuclear accident simulations using ASTEC/SOPHAEROS software. Interaction energies, standard enthalpies and Gibbs free energies of successive water molecules addition determine the ordering of the complexes. CsOH forms the most hydrated stable complexes followed by CsI, Cs2I2, and Cs. CsOH still exists in steam atmosphere even at quite high temperature, up to around 1100 K.

A theoretical study of the H-abstraction reactions from HOI by moist air radiolytic products (H, OH, and O (3P)) and iodine atoms (2P(3/2)).

The rate constants of the reactions of HOI molecules with H, OH, O ((3)P), and I ((2)P(3/2)) atoms have been estimated over the temperature range 300-2500 K using four different levels of theory. Geometry optimizations and vibrational frequency calculations are performed using MP2 methods combined with two basis sets (cc-pVTZ and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVTZ, cc-pVQZ, 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Reaction enthalpies at 0 K were calculated at the CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory and compared to the experimental values taken from the literature. Canonical transition-state theory with an Eckart tunneling correction is used to predict the rate constants as a function of temperature. The computational procedure has been used to predict rate constants for H-abstraction elementary reactions because there are actually no literature data to which the calculated rate constants can be directly compared. The final objective is to implement kinetics of gaseous reactions in the ASTEC (accident source term evaluation code) program to improve speciation of fission products, which can be transported along the reactor coolant system (RCS) of a pressurized water reactor (PWR) in the case of a severe accident.

Theoretical study of the gas-phase reactions of iodine atoms ((2)P(3/2)) with H(2), H(2)O, HI, and OH.

The rate constants of the reactions of iodine atoms with H(2), H(2)O, HI, and OH have been estimated using 39, 21, 13, and 39 different levels of theory, respectively, and have been compared to the available literature values over the temperature range of 250-2500 K. The aim of this methodological work is to demonstrate that standard theoretical methods are adequate to obtain quantitative rate constants for the reactions involving iodine-containing species. Geometry optimizations and vibrational frequency calculations are performed using three methods (MP2, MPW1K, and BHandHLYP) combined with three basis sets (cc-pVTZ, cc-pVQZ, and 6-311G(d,p)). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the cc-pVnZ (n = T, Q, and 5), aug-cc-pVnZ (n = T, Q, and 5), 6-311G(d,p), 6-311+G(3df,2p), and 6-311++G(3df,3pd) basis sets. Canonical transition state theory with a simple Wigner tunneling correction is used to predict the rate constants as a function of temperature. CCSD(T)/cc-pVnZ//MP2/cc-pVTZ (n = T and Q), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), and CCSD(T)/6-311++G(3df,3pd)//MP2/6-311G(d,p) levels of theory provide accurate kinetic rate constants when compared to available literature data. The use of the CCSD(T)/cc-pVQZ//MP2/cc-pVTZ and CCSD(T)/6-311++G(3df,3pd) levels of theory allows one to obtain a better agreement with the literature data for all reactions with the exception of the I + H(2) reaction R(1) . This computational procedure has been also used to predict rate constants for some reactions where no available experimental data exist. The use of quantum chemistry tools could be therefore extended to other elements and next applied to develop kinetic networks involving various fission products, steam, and hydrogen in the absence of literature data. The final objective is to implement the kinetics of gaseous reactions in the ASTEC (Accident Source Term Evaluation Code) code to improve speciation of fission transport, which can be transported along the Reactor Coolant System (RCS) of a Pressurized Water Reactor (PWR) in case of a severe accident.