Water versus acetonitrile coordination to uranyl. Effect of chloride ligands

Optimizations at the BLYP and B3LYP levels are reported for the mixed uranyl chloro/water/acetonitrile complexes [UO(2)Cl(n)(H(2)O)(x)(MeCN)(5-n-x)](2-n) (n = 1-3) and [UO(2)Cl(n)(H(2)O)(x)(MeCN)(4-n-x)](2-n) (n = 2-4), in both the gas phase and a polarizable continuum modeling acetonitrile. Car-Parrinello molecular dynamics (CPMD) simulations have been performed for [UO(2)Cl(2)(H(2)O)(MeCN)(2)] in the gas phase and in a periodic box of liquid acetonitrile. According to population analyses and dipole moments evaluated from maximally localized Wannier function centers, uranium is less Lewis acidic in the neutral UO(2)Cl(2) than in the UO(2)(2+) moiety. In the gas phase the latter binds acetonitrile ligands more strongly than water, whereas in acetonitrile solution, the trend is reversed due to cooperative polarization effects. In the polarizable continuum the chloro complexes have a slight energetic preference for water over acetonitrile ligands, but several mixed complexes are so close in free energy ΔG that they should exist in equilibrium, in accord with previous interpretations of EXAFS data in solution. The binding strengths of the fifth neutral ligands decrease with increasing chloride content, to the extent that the trichlorides should be formulated as four-coordinate [UO(2)Cl(3)L](-) (L = H(2)O, MeCN). Limitations to their accuracy notwithstanding, density functional calculations can offer insights into the speciation of a complex uranyl system in solution, a key feature in the context of nuclear waste partitioning by complexant molecules.     

Redox energetics and kinetics of uranyl coordination complexes in aqueous solution

Detailed voltammetric results for five uranyl coordination complexes are presented and analyzed using digital simulations of the voltammetric data to extract thermodynamic (E(1/2)) and heterogeneous electron-transfer kinetic (k(0) and alpha) parameters for the one-electron reduction of UO(2)(2+) to UO(2)(+). The complexes and their corresponding electrochemical parameters are the following: [UO(2)(OH(2))(5)](2+) (E(1/2) = -0.169 V vs Ag/AgCl, k(0) = 9.0 x 10(-3) cm/s, and alpha = 0.50); [UO(2)(OH)(5)](3-) (-0.927 V, 2.8 x 10(-3) cm/s, 0.46); [UO(2)(C(2)H(3)O(2))(3)](-) (-0.396 V, approximately 0.1 cm/s, approximately 0.5); [UO(2)(CO(3))(3)](4-) (-0.820 V, 2.6 x 10(-5) cm/s, 0.41); [UO(2)Cl(4)](2-) (-0.065 V, 9.2 x 10(-3) cm/s, 0.30). Differences in the E(1/2) values are attributable principally to differences in the basicity of the equatorial ligands. Differences in rate constants are considered within the context of Marcus theory of electron transfer, but no specific structural change(s) in the complexes between the two oxidation states can be uniquely identified with the underlying variability in the heterogeneous rate constants and electron-transfer coefficients.     

Binding of pertechnetate to uranyl(VI) in aqueous solution. A density functional theory molecular dynamics study

According to constrained Car-Parrinello molecular dynamics simulations and thermodynamic integration, the free binding energy between uranyl hydrate and pertechnetate in aqueous solution is significantly lower than that between uranyl and nitrate, namely, by 1.7 kcal mol(-1). This is the first study of the differential binding of these two ligands to uranyl, which can have implications for the separability of uranium and technetium during the reprocessing of nuclear waste.     

Density functional study of the lithiation of cyclic vinyl ethers in solution

Lithiation of three cyclic vinyl ethers—2,3-dihydrofuran, 3,4-dihydro-2H-pyran, and 2,3-dihydrooxepin in ethereal solution are investigated at the density functional theory level of B3LYP/6-31++G(d,p). Several solvation models were used, including the microsolvation model, the continuum models such as the Onsager model, the polarized continuum models, and the isodensity polarized continuum model, and the mixed discrete-continuum model. Both the microsolvation and the mixed discrete-continuum model gave results consistent with experiments. Theoretical calculations also indicate that lithiation of 2,3-dihydrooxepin undergoes allylic lithiation concomitant with ring opening.     

A density functional study of uranyl monocarboxylates

We studied uranium(VI) monocarboxylate complexes by a relativistic density functional method using simple carboxylic acids as ligands, i.e. [UO2(OOCR)]+ (R = H, CH3, CH2CH3). These complexes exist in aqueous solution and, for R = CH3 and CH2CH3, may also be considered as models of uranyl complexated by humic substances. We investigated mono- and bidentate coordination modes. Short-range solvent effects were accounted for explicitly via aqua ligands of the first hydration shell and long-range electrostatic interactions were described via a polarizable continuum model. The calculated results for the uranyl U=O bond, the bond to aqua ligands, and the averaged uranium distances to equatorial oxygen atoms, U-Oeq, agreed quite well with EXAFS-derived interatomic distances. However, the uranyl-carboxylate bond was calculated to be notably shorter than the experimentally determined value. Experimental differences between mono- and bidentate coordination, obtained mainly from crystal structures, were qualitatively reproduced for the U-C distance but not for the average bond length, U-Oeq. We discuss these discrepancies between calculated and experimental results in some detail and suggest changes in the coordination number rather than variations of the coordination geometry as the main source of the experimentally observed variation of the U-Oeq distance.     

Raman study of vibrational relaxation of acetonitrile in solution

Line widths of isotropic Raman spectra of the ν₁ (a₁) CH or CD stretching bands of acetonitrile CH₃CN and CD₃CN were measured in a number of solvents. The vibrational dephasing theory, modified for use in binary mixtures, predicts quantitatively the solvent dependence of the Raman line widths.     

Speciation and coordination chemistry of uranyl(VI)-citrate complexes in aqueous solution

The pH dependence of uranyl(VI) complexation by citric acid was investigated using Raman and attenuated total reflection FTIR spectroscopies and electrospray ionization mass spectrometry. pH-dependent changes in the nu(s)(UO(2)) envelope indicate that three major UO(2)(2+)-citrate complexes with progressively increasing U=O bond lengths are present over a range of pH from 2.0 to 9.5. The first species, which is the predominant form of uranyl(VI) from pH 3.0 to 5.0, contains two UO(2)(2+) groups in spectroscopically equivalent coordination environments and corresponds to the [(UO(2))(2)Cit(2)](2)(-) complex known to exist in this pH range. At pH values >6.5, [(UO(2))(2)Cit(2)](2)(-) undergoes an interconversion to form [(UO(2))(3)Cit(3)](3)(-) and (UO(2))(3)Cit(2). ESI-MS studies on solutions of varying uranyl(VI)/citrate ratios, pH, and solution counteranion were successfully used to confirm complex stoichiometries. Uranyl and citrate concentrations investigated ranged from 0.50 to 50 mM.     

Density functional perturbational orbital theory of spin polarization in electronic systems. I. Formalism

A perturbational approach is presented for the general analysis of spin-polarization effect on electronic structures and energies within spin-density functional formalism. Explicit expressions for the changes in Kohn-Sham [Phys. Rev. 140, 1133 (1965)] orbital energies and coefficients as well as for the change in total electronic energy are derived upon using the local spin density and self-interaction-corrected exchange-correlation functionals. The application of the method for atoms provides analytical expressions for the exchange splitting energy and spin-polarization energy. The atomic exchange parameters are obtained from the expressions for the elements with Z=1-92 and they match well with Stoner exchange parameters for 3d metal elements.     

Rate-determining cooperative effects of bimolecular reactions in solution

An expression for the rate constant of condensed-phase bimolecular reactions is derived. The key feature of the proposed model is the formulation of the energy-dependent rate constant in terms of the diffusion rate and the ratio of the volume in phase space that leads to product over the total volume. The dependence of the bimolecular rate constant by the reduced barrier x = E thermod/kT is given in explicit form in terms of the incomplete and the complete gamma functions of Euler. The performance of the proposed model is tested against the experimental rate constants for the Menschutkin reaction by fitting the parameters of the expression for the rate constant to experimental data at various temperatures. The potential energy barrier obtained from the regression (16.75 kcal mol-1) is close to the independently computed value at the CPCM B3LYP/CRENBL.6-311(+)G(d) level of theory (16.84 kcal mol-1). The corresponding fitting to the transition state theory expression affords the lower value of 14.65 kcal mol-1.     

Aqueous coordination chemistry and photochemistry of uranyl(VI) oxalate revisited: a density functional theory study

Using density functional theory (DFT) calculations, we revisited a classical problem of uranyl(VI) oxalate photochemical decomposition. Photoreactivities of uranyl(VI) oxalate complexes are found to correlate largely with ligand-structural arrangements. Importantly, the intramolecular photochemical reaction is inhibited when oxalate is bound to uranium exclusively in chelate binding mode. Previously proposed mechanisms involving a UO(2)(C(2)O(4))(2)(2-) (1:2) complex as the main photoreactive species are thus unlikely to apply, because the two oxalic acids are bound to uranium in a chelating binding mode. Our DFT results suggest that the relevant photoreactive species are UO(2)(C(2)O(4))(3)(4-) (1:3) and (UO(2))(2)(C(2)O(4))(5)(6-) (2:5) complexes binding uranium in an unidentate fashion. These species go through decarboxylation upon excitation to the triplet state, which ensues the release of CO(2) and reduction of U(vi) to U(v). The calculations also suggest an alternative intermolecular pathway at low pH via an electron transfer between the excited state *UO(2)(2+) and hydrogen oxalate (HC(2)O(4)(-)) which eventually leads to the production of CO and OH(-) with no net reduction of U(VI). The calculated results are consistent with previous experimental findings that CO is only detected at low pH while U(IV) is detected only at high pH.     

Variable denticity in carboxylate binding to the uranyl coordination complexes

Tris-carboxylate complexes of uranyl [UO₂]²⁺ with acetate and benzoate were generated using electrospray ionization mass spectrometry, and then isolated in a Fourier transform ion cyclotron resonance mass spectrometer. Wavelength-selective infrared multiple photon dissociation (IRMPD) of the tris-acetato uranyl anion resulted in a redox elimination of an acetate radical, which was used to generate an IR spectrum that consisted of six prominent absorption bands. These were interpreted with the aid of density functional theory calculations in terms of symmetric and antisymmetric −CO₂ stretches of the monodentate and bidentate acetate, CH₃ bending and umbrella vibrations, and a uranyl O—U—O asymmetric stretch. The comparison of the calculated and measured IR spectra indicated that the predominant conformer of the tris-acetate complex contained two acetate ligands bound in a bidentate fashion, while the third acetate was monodentate. In similar fashion, the tris-benzoate uranyl anion was formed and photodissociated by loss of a benzoate radical, enabling measurement of the infrared spectrum that was in close agreement with that calculated for a structure containing one monodentate and two bidentate benzoate ligands.     

Density functional model studies of uranyl adsorption on (001) surfaces of kaolinite

The adsorption of uranyl on two types of neutral (001) surfaces of kaolinite, tetrahedral Si(t) and octahedral Al(o), was studied by means of density functional periodic slab model calculations. Various types of model surface complexes, adsorbed at different sites, were optimized and adsorption energies were estimated. As expected, the Si(t) surface was found to be less reactive than the Al(o) surface. At the neutral Al(o) surface, only adsorption at protonated sites is calculated to be exothermic for inner- as well as outer-sphere adsorption complexes, with monodentate coordination being preferred. Adsorption energies as well as structural features of the adsorption complexes are mainly determined by the number of deprotonated surface hydroxyl groups involved. Outer-sphere complexes on both surfaces exhibit a shorter U-O bond to the aqua ligand of uranyl that is in direct contact with the surface than to the other aqua ligands. This splitting of the shell of equatorial U-O bonds is at variance with common expectations for outer-sphere surface complexes of uranyl.     

Density functional theory study of uranium(VI) aquo chloro complexes in aqueous solution

Mixed uranyl aquo chloro complexes of the type [UO2(H2O)xCly]2-y (y = 1, 2, 3, 4; x + y = 4, 5) have been optimized at the BLYP, BP86, and B3LYP levels of density functional theory in vacuo and in a polarizable continuum modeling bulk water (PCM) and have been studied at the BLYP level with Car-Parrinello molecular dynamics (MD) simulations in the gas phase and in explicit aqueous solution. Free binding energies were evaluated from static PCM data and from pointwise thermodynamic integration involving constrained MD simulations in water. The computations reveal significant solvent effects on geometric and energetic parameters. Based on the comparison of PCM-optimized or MD-averaged uranyl-ligand bond distances with EXAFS-derived values, the transition between five- and four-coordination about uranyl is indicated to occur at a Cl content of y = 2 or 3.     

Density functional theory study of difluorovinylidene isomerization to difluoroacetylene

The isomerization of difluorovinylidene: C = CF₂ into difluoroacetylene FCCF was systematically studied using density functional theory (DFT) methods with the 6-311 + G(2 p, 2 p) Gaussian-type basis set. The computed geometries for the reactant, product, and transition state as well as the enthalpy and reaction barrier were compared with the previously computed values obtained by the LSD and NLSD DFT methods and the ab initio CCSD method. The suitability of some of the DFT methods for the computation of this system was discussed and the best estimated computational data were selected. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 62: 515–520, 1997     

A comparative electrochemical study of diffusion in room temperature ionic liquid solvents versus acetonitrile.

Measurements on the diffusion coefficient of the neutral molecule N,N,N',N'-tetramethyl-para-phenylenediamine and the radical cation and dication generated by its one- and two-electron oxidation, respectively, are reported over the range 298-348 K in both acetonitrile and four room temperature ionic liquids (RTILs). Data were collected using single and double potential step chronoamperometry at a gold disk electrode of micrometer dimension, and analysed via fitting to the appropriate analytical expression or, where necessary, to simulation. The variation of diffusion coefficient with temperature was found to occur in an Arrhenius-type manner for all combinations of solute and solvent. For a given ionic liquid, the diffusional activation energies of each species were not only closely equivalent to each other, but also to the RTIL's activation energy of viscous flow. In acetonitrile supported with 0.1 M tetrabutylammonium perchlorate, the ratio in diffusion coefficients of the radical cation and dication to the neutral molecule were calculated as 0.89 +/- 0.05 and 0.51 +/- 0.03, respectively. In contrast, amongst the ionic liquids the same ratios were determined to be on average 0.53 +/- 0.04 and 0.33 +/- 0.03. The consequences of this dissimilarity are considered in terms of the modelling of voltammetric data gathered within ionic liquid solvents.     

Density functional perturbational orbital theory of spin polarization in electronic systems. II. Transition metal dimer complexes

We present a theoretical scheme for a semiquantitative analysis of electronic structures of magnetic transition metal dimer complexes within spin density functional theory (DFT). Based on the spin polarization perturbational orbital theory [D.-K. Seo, J. Chem. Phys. 125, 154105 (2006)], explicit spin-dependent expressions of the spin orbital energies and coefficients are derived, which allows to understand how spin orbitals form and change their energies and shapes when two magnetic sites are coupled either ferromagnetically or antiferromagnetically. Upon employment of the concept of magnetic orbitals in the active-electron approximation, a general mathematical formula is obtained for the magnetic coupling constant J from the analytical expression for the electronic energy difference between low-spin broken-symmetry and high-spin states. The origin of the potential exchange and kinetic exchange terms based on the one-electron picture is also elucidated. In addition, we provide a general account of the DFT analysis of the magnetic exchange interactions in compounds for which the active-electron approximation is not appropriate.