F-atom concentration measurement via HCl titration

Titration of F atoms with HCl is shown to be described principally by the reaction, F + HCl → HF + Cl. This conclusion is arrived at by observing the effect of HCl titration on the intensity of an F-atom ESR signal, the signal intensity being proportional to F-atom concentration. The intensity decreases linearly with increasing HCl flow and (1.07 ± 0.12) HCl molecules are required for removal of each F atom. Since the titration volume contains F₂ as well as F, this latter result demonstrates the secondary importance of the competing reaction, Cl + F₂ → ClF + F. Accordingly, F-atom concentration may be determined in the presence of F₂ by addition of an excess amount of HCl, thereby converting the F present to HF, and measurement of the amount of HF present by laser optical absorption methods. This alternative (to ESR) technique is important with regard to evaluating HF/DF chemical laser systems where F and F₂ are present.     

Solvation of uranyl(II), europium(III) and europium(II) cations in "basic" room-temperature ionic liquids: a theoretical study

We report a molecular dynamics study of the solvation of UO2(2+), Eu3+ and Eu2+ ions in two "basic" (Lewis acidity) room-temperature ionic liquids (IL) composed of the 1-ethyl-3-methylimidazolium cation (EMI+) and a mixture of AlCl4- and Cl- anions, in which the Cl-/AlCl4- ratio is about 1 and 3, respectively. The study reveals the importance of the [UO2Cl4]2- species, which spontaneously form during most simulations, and that the first solvation shell of europium is filled with Cl- and AlCl4- ions embedded in a cationic EMI+ shell. The stability of the [UO2Cl4]2- and [Eu(III)Cl6]3- complexes is supported by quantum mechanical calculations, according to which the uranyl and europium cations intrinsically prefer Cl- to the AlCl4- ion. In the gas phase, however, [Eu(III)Cl6]3- and [Eu(II)Cl6]4- complexes are predicted to be metastable and to lose two to three Cl- ions. This contrasts with the results of simulations of complexes in ILs, in which the "solvation" of the europium complexes increases with the number of coordinated chlorides, leading to an equilibrium between different chloro species. The behavior of the hydrated [Eu(OH2)8]3+ complex is considered in the basic liquids; the complex exchanges H2O molecules with Cl- ions to form mixed [EuCl3(OH2)4] and [EuCl4(OH2)3]- complexes. The results of the simulations allow us to better understand the microscopic nature and solvation of lanthanide and actinide complexes in "basic" ionic liquids.     

Electrochemical kinetics of the hydrogen reaction on platinum in concentrated HCl(aq)

The hydrogen reaction in concentrated HCl(aq) solutions is a key reaction for the CuCl(aq)/HCl(aq) electrolytic cell. Here, electrochemical impedance spectroscopy (EIS) and linear sweep voltammetry (LSV) were used to obtain new data for the hydrogen reaction on platinum submerged in highly concentrated acidic solutions at 25 °C and 0.1 MPa. LSV and EIS data were collected for Pt in 0.5 mol/L H₂SO₄(aq), 1 mol/L HCl(aq) and 7.71 mol/L HCl(aq) solutions. It was found that exchange current density (j₀) values varied between 1 and 2 mA/cm². An equivalent circuit model was used to obtain comparable j₀ and limiting current density values from EIS data relative to values obtained with LSV data. It was found that as the concentration of acid increased, a noticeable decrease in the performance was observed.     

Calculating geochemical reaction pathways--exploration of the inner-sphere water exchange mechanism in Al(H2O)6(3+)(aq) + nH2O with ab Initio calculations and molecular dynamics

We have simulated exchange of inner-sphere and bulk water molecules for different sizes of Al3+(aq) clusters, Al(H2O)63+ + nH2O for n = 0, 1, 6, or 12, with ab initio and molecular dynamics simulations, in order to understand how robust the ab initio method is for identifying hydrolytic reaction pathways of particular importance to geochemistry. In contrast to many interfacial reactions, this particular elementary reaction is particularly simple and well-constrained by experiment. Nevertheless, we find that a rich array of parallel reaction pathways depend sensitively on the details of the solvation sphere and structure and that larger clusters are not necessarily better. Inner-sphere water exchange in Al3+(aq) may occur through two Langford-Gray dissociative pathways, one in which the incoming and outgoing waters are cis, the other in which they are trans to one another. A large majority of exchanges in the molecular dynamics simulations occurred via the trans mechanism, in contrast to the predictions of the ab initio method. In Al(H2O)63+ + H2O, the cis mechanism has a transition state of 84.3 kJ/mol, which is in good agreement with previous experimental and ab initio results, while the trans mechanism has only a saddle point with two negative frequencies, not a transition state, at 89.7 kJ/mol. In addition to the exchange mechanisms, dissociation pathways could be identified that were considerably lower in energy than experiment and varied considerably between 60 and 100 kJ/mol, depending on the particular geometry and cluster size, with no clear relation between the two. Ab initio calculations using large clusters with full second coordination spheres (n = 12) were unable to find dissociation or exchange transition states because the network of hydrogen bonds in the second coordination sphere was too rigid to accommodate the outgoing inner-sphere water. Our results indicate that caution should surround ab initio simulation of complicated dynamic processes such as hydrolysis, ion exchange, and interfacial reactions that involve several steps. Dynamic methods of simulation need to accompany static methods such as ab initio calculation, and it is best to consider simulated pathways as hypotheses to be tested experimentally rather than definitive properties of the reaction.     

Rates of oxygen-isotope exchange between sites in the [HxTa6O19](8-x)-(aq) Lindqvist ion and aqueous solutions: comparisons to [HxNb6O19](8-x)-(aq)

Rates of steady oxygen-isotope exchange differ in interesting ways for two sets of structural oxygens in the [HxTa6O19](8-x)-(aq) Lindqvist ion when compared to published data on the [HxNb6O19]8-x(aq) version. Because of the lanthanide contraction, the [HxTa6O19](8-x)-(aq) and [HxNb6O19](8-x)-(aq) ions are virtually isostructural and differ primarily in a full core (Kr vs Xe) and the 4f14 electrons in the [HxTa6O19](8-x)-(aq) ion. For both molecules, both pH-dependent and -independent pathways are evident in isotopic exchange of the 12 mu2-O(H) and 6 eta=O sites. Rate parameters for eta=O exchange at conditions where there is no pH dependence are, for the Ta(V) and Nb(V) versions respectively, K(298)(0) = 2.72 x 10(-5) s(-1) and 9.7 x 10(-6) s(-1), DeltaH = 83.6 +/- 3.2 and 89.4 kJ.mol(-1), and DeltaS = -51.0 +/- 10.6 and -42.9 J.mol(-1).K-1. For the mu2-O sites, K(298)(0) = 1.23 x 10(-6) s(-1), DeltaH = 70.3 +/- 9.7 and 88.0 kJ.mol(-1), and DeltaS = -116.1 +/- 32.7 and -29.4 J.mol(-1).K-1. Protonation of the 6 eta=O sites is energetically unfavored relative to the 12 mu2-O bridges in both molecules, although not equally so. Experimentally, protonation labilizes both the mu2-O(H) and eta=O sites to isotopic exchange in both molecules. Density-functional electronic-structure calculations indicate that proton affinities of structural oxygens in the two molecules differ with the [HxTa6O19](8-x)-(aq) anion having a smaller affinity to protonate than the [HxNb6O19]8-x(aq) ion. This difference in proton affinities is evident in the solution chemistry as pKa = 11.5 for the [HTa6O19]7-(aq) ion and pKa = 13.6 for the [HNb6O19]7-(aq) ion. Most striking is the observation that eta=O sites isotopically equilibrate faster than the mu2-O sites for the [HxTa6O19](8-x)-(aq) Lindqvist ion but slower for the [HxNb6O19](8-x)-(aq) ion, indicating that predictions about site reactivities in complicated structures, such as the interface of aqueous solutions and oxide solids, should be approached with great caution.     

(HCl)2 and (HF)2 in small helium clusters: quantum solvation of hydrogen-bonded dimers

We present a rigorous theoretical study of the solvation of (HCl)(2) and (HF)(2) by small ((4)He)(n) clusters, with n=1-14 and 30. Pairwise-additive potential-energy surfaces of He(n)(HX)(2) (X=Cl and F) clusters are constructed from highly accurate four-dimensional (rigid monomer) HX-HX and two-dimensional (rigid monomer) He-HX potentials and a one-dimensional He-He potential. The minimum-energy geometries of these clusters, for n=1-6 in the case of (HCl)(2) and n=1-5 for (HF)(2), correspond to the He atoms in a ring perpendicular to and bisecting the HX-HX axis. The quantum-mechanical ground-state energies and vibrationally averaged structures of He(n)(HCl)(2) (n=1-14 and 30) and He(n)(HF)(2) (n=1-10) clusters are calculated exactly using the diffusion Monte Carlo (DMC) method. In addition, the interchange-tunneling splittings of He(n)(HCl)(2) clusters with n=1-14 are determined using the fixed-node DMC approach, which was employed by us previously to calculate the tunneling splittings for He(n)(HF)(2) clusters, n=1-10 [A. Sarsa et al., Phys. Rev. Lett. 88, 123401 (2002)]. The vibrationally averaged structures of He(n)(HX)(2) clusters with n=1-6 for (HCl)(2) and n=1-5 for (HF)(2) have the helium density localized in an effectively one-dimensional ring, or doughnut, perpendicular to and at the midpoint of the HX-HX axis. The rigidity of the solvent ring varies with n and reaches its maximum for the cluster size at which the ring is filled, n=6 and n=5 for (HCl)(2) and (HF)(2), respectively. Once the equatorial ring is full, the helium density spreads along the HX-HX axis, eventually solvating the entire HX dimer. The interchange-tunneling splitting of He(n)(HCl)(2) clusters hardly varies at all over the cluster size range considered, n=1-14, and is virtually identical to that of the free HCl dimer. This absence of the solvent effect is in sharp contrast with our earlier results for He(n)(HF)(2) clusters, which show a approximately 30% reduction of the tunneling splitting for n=4. A tentative explanation for this difference is proposed. The implications of our results for the interchange-tunneling dynamics of (HCl)(2) in helium nanodroplets are discussed.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.     

The low-temperature (5 to 310 K) heat capacity and thermodynamic functions of cesium fluoroxysulfate (CsSO4F) and the standard potential at 298.15 K for the half-reaction: SO4F−(aq) + 2H+(aq) + 2e− = HSO4−(aq) + HF(aq)

The low-temperature (5 to 310 K) heat capacity of cesium fluoroxysulfate, CsSO₄F, has been measured by adiabatic calorimetry. At T = 298.15 K, the heat capacity C_p^o(T) and standard entropy S^o(T) are (163.46±0.82) and (201.89±1.01) J · K^?1 · mol^?1, respectively. Based on an earlier measurement of the standard enthalpy of formation ΔH_f^o the Gibbs energy of formation ΔG_f^o(CsSO₄F, c, 298.15 K) is calculated to be ?(877.6±1.6) kJ · mol^?1. For the half-reaction: SO₄F^?(aq)+2H⁺(aq)+2e^? = HSO₄^?(aq)+HF(aq), the standard electrode potential E at 298.15 K, is (2.47±0.01) V.     

Unimolecular rate constants for HX or DX elimination (X = F, Cl) from chemically activated CF3CH2CH2Cl, C2H5CH2Cl, and C2D5CH2Cl: threshold energies for HF and HCl elimination

Vibrationally activated CF(3)CH(2)CH(2)Cl molecules were prepared with 94 kcal mol(-1) of vibrational energy by the combination of CF(3)CH(2) and CH(2)Cl radicals and with 101 kcal mol(-1) of energy by the combination of CF(3) and CH(2)CH(2)Cl radicals at room temperature. The unimolecular rate constants for elimination of HCl from CF(3)CH(2)CH(2)Cl were 1.2 x 10(7) and 0.24 x 10(7) s(-1) with 101 and 94 kcal mol(-1), respectively. The product branching ratio, k(HCl)/k(HF), was 80 +/- 25. Activated CH(3)CH(2)CH(2)Cl and CD(3)CD(2)CH(2)Cl molecules with 90 kcal mol(-1) of energy were prepared by recombination of C(2)H(5) (or C(2)D(5)) radicals with CH(2)Cl radicals. The unimolecular rate constant for HCl elimination was 8.7 x 10(7) s(-1), and the kinetic isotope effect was 4.0. Unified transition-state models obtained from density-functional theory calculations, with treatment of torsions as hindered internal rotors for the molecules and the transition states, were employed in the calculation of the RRKM rate constants for CF(3)CH(2)CH(2)Cl and CH(3)CH(2)CH(2)Cl. Fitting the calculated rate constants from RRKM theory to the experimental values provided threshold energies, E(0), of 58 and 71 kcal mol(-1) for the elimination of HCl or HF, respectively, from CF(3)CH(2)CH(2)Cl and 54 kcal mol(-1) for HCl elimination from CH(3)CH(2)CH(2)Cl. Using the hindered-rotor model, threshold energies for HF elimination also were reassigned from previously published chemical activation data for CF(3)CH(2)CH(3,) CF(3)CH(2)CF(3), CH(3)CH(2)CH(2)F, CH(3)CHFCH(3), and CH(3)CF(2)CH(3). In an appendix, the method used to assign threshold energies was tested and verified using the combined thermal and chemical activation data for C(2)H(5)Cl, C(2)H(5)F, and CH(3)CF(3).     

Cu(II) in liquid ammonia: an approach by hybrid quantum-mechanical/molecular-mechanical molecular dynamics simulation.

To investigate the solvation structure of the Cu(II) ion in liquid ammonia, ab initio quantum-mechanical/molecular-mechanical (QM/MM) molecular dynamics (MD) simulations were carried out at Hartree Fock (HF) and hybrid density functional theory (B3 LYP) levels. A sixfold-coordinated species was found to be predominant in the HF case whereas five- and sixfold-coordinated complexes were obtained in a ratio 2:1 from the B3 LYP simulation. In contrast to hydrated Cu(II), which exhibits a typical Jahn-Teller distortion, the geometrical arrangement of ligand molecules in the case of ammonia can be described as a [2 + 4] ([2 + 3]) configuration with 4 (3) elongated copper-nitrogen bonds. First shell solvent exchange reactions at picosecond rate took place in both HF and B3 LYP simulations, again in contrast to the more stable sixfold-coordinated hydrate. NH3 ligands apparently lead to strongly accelerated dynamics of the Cu(II) solvate due to the "inverse" [2 + 4] structure with its larger number of elongated copper-ligand bonds. Several dynamical properties, such as mean ligand residence times or ion-ligand stretching frequencies, prove the high lability of the solvated complex.     

Microsolvation of HCl within cold NH3 clusters

The solid state solvation of HCl molecules with small ammonia clusters at an average temperature of 100 K was investigated by on-the-fly molecular dynamics methodology. Structures close to the proton jump from HCl molecule to the ammonia have been further checked with the MP2/aug-cc-pvDZ calculations. Ionization of HCl and/or sharing of the proton were found. Two Zundel-type ions were observedone with proton being shared between ammonium ion and Cl (-) anion (Cl (-)...H (+)...NH 3) in all complexes, and the second, between hydrogen chloride and Cl (-) anion in the HCl...Cl (-)...NH 4 (+)...(NH 3) 2 complex. However, in contrast to methanol clusters, ammonia clusters are not good for the proton wires since once the proton moves to ammonia, it is localized on the ammonium ion units.     

Osmotic coefficients of atomistic NaCl (aq) force fields

Solvated ions are becoming increasingly important for (bio)molecular simulations. But there are not much suitable data to validate the intermediate-range solution structure that ion-water force fields produce. We compare six selected combinations of four biomolecular Na-Cl force fields and four popular water models by means of effective ion-ion potentials. First we derive an effective potential at high dilution from simulations of two ions in explicit water. At higher ionic concentration multibody effects will become important. We propose to capture those by employing a concentration dependent dielectric permittivity. With the so obtained effective potentials we then perform implicit solvent simulations. We demonstrate that our effective potentials accurately reproduce ion-ion coordination numbers and the local structure. They allow us furthermore to calculate osmotic coefficients that can be directly compared with experimental data. We show that the osmotic coefficient is a sensitive and accurate measure for the effective ion-ion interactions and the intermediate-range structure of the solution. It is therefore a suitable and useful quantity for validating and parametrizing atomistic ion-water force fields.     

Classical and quantum studies of the photodissociation of a HX (X=Cl,F) molecule adsorbed on ice

The photodissociation dynamics of a HX (X = Cl,F) molecule adsorbed on a hexagonal ice surface at T = 0 K is studied using time-dependent quantum wave packets and quasiclassical trajectories. The relevant potential energy surfaces are calculated using high-level ab initio methods. We present here two dimensional calculations for the dynamics of the hydrogen photofragment for both HCl and HF molecules. The purpose of this paper is to compare the photodissociation dynamics of the two molecules which are adsorbed on the ice surface with different equilibrium geometries. The total photodissociation cross section and the angular distribution are calculated. The comparison with classical trajectory calculations provides evidence for typical quantum effects and reveals rainbow structures.     

Using ionic liquids to trap unique coordination environments: polymorphic solvates of ErCl3(OH2)(4).2([C2mim]Cl)

Two polymorphs of ErCl(3)(OH(2))(4).2([C(2)mim]Cl) solvates were isolated from the same solution of 1-ethyl-3-methylimidazolium chloride when HCl(aq) was added, while [C(2)mim](3)[ErCl(6)] was isolated without HCl addition, illustrating how ionic liquids can be used to trap unusual coordination environments in the solid state.     

The nature of superacid electrophilic species in HF/SbF(5): a density functional theory study

A density functional theory study at the B3LYP/6-31++G** + RECP(Sb) level of the HF/SbF(5) superacid system was carried out. The geometries of possible electrophilic species, such as H(2)F(+).Sb(2)F(11)(-) and H(3)F(2)(+).Sb(2)F(11)(-), were calculated and correspond with available experimental results. Calculations of different equilibrium reactions involving HF and SbF(5) allowed the relative concentration of the most energetically favorable species present in 1:1 HF/SbF(5) solutions to be estimated. These species are H(+).Sb(2)F(11)(-), H(2)F(+).Sb(2)F(11)(-), H(3)F(2)(+).Sb(2)F(11)(-), and H(4)F(3)(+).Sb(2)F(11)(-), which correspond to 36.9, 16.8, 36.9, and 9.4%, respectively. Calculations of the acid strength of the electrophilic species were also performed and indicated that, for the same anion, the acid strength increases with the solvation degree. The entropic term also plays a significant role in proton-transfer reactions in superacid systems.     

Electron detachment and relaxation of OH-(aq)

Femtosecond transient absorption spectroscopy is used to study the primary reaction dynamics of photoinduced electron detachment of the hydroxide ion in water, OH- (aq). The electron is detached by excitation of OH- (aq) to the charge-transfer-to-solvent (CTTS) state at 200 nm. The subsequent relaxation processes are probed in the spectral range from 193 to 800 nm with femtosecond time resolution. We determine both the time-dependent quantum yields of OH- (aq), OH(aq), and e-(aq), and we observe a transient spectral signature which is assigned to relaxation of hot (OH-)* ions formed via solvent-assisted conversion of the excited CTTS state to OH-. The primary quantum yield of OH(aq) is 65 +/- 5%, while recombination with e-(aq) reduces the yield to 34% after 5 ps and 12% after 200 ps. The yield of hot (OH-)* ions is 35 +/- 5%. Rotational anisotropy measurements of OH- (aq) and OH(aq) indicate a reorientation time for OH- (aq) of 1.9 ps, while no rotational anisotropy is resolved for the OH(aq) radical within our time resolution of 0.3 ps. This is consistent with the notion that OH(aq) radicals formed after electron detachment are only weakly bound to the hydrogen bond network of water. The assignment of the experimental data is supported by a series of electronic structure calculations of simple complexes of OH- (H(2)O)(n).     

An ab initio study on HXC(double bond)O …︁ HY molecular complexes (X,Y = F,Cl)

Post-Hartree-Fock calculations were carried out to predict the stabilities and properties of four HClCO …︁ HCl, HClCO …︁ HF, HFCO …︁ HCl, and HFCO …︁ HF molecular complexes. Full geometry optimizations and vibrational frequency calculations were performed for all systems using standard 6-311G (d, p) and 6-311G (2d,2p) basis sets at the MP 2 level of theory. Single-point calculations of the interaction energies were carried out for all complexes at the MP 4(SDTQ ) level with the 6-311G (d, p) basis set. All systems were found to be stable and their predicted molecular parameters match well available (very scarce) experimental data. © 1996 John Wiley & Sons, Inc.     

Coupled reference interaction site model/simulation approach for thermochemistry of solvation: theory and prospects

We present a new methodology for computing solvation free energy, which is based upon the reference interaction site model (RISM)/hypernetted chain (HNC) solvation free energy expression, but which substitutes radial distribution functions taken from simulations for those calculated by simultaneous solution of the RISM and HNC equations. Consequently, solvation free energy can be obtained from a single molecular dynamics or Monte Carlo simulation. Here we describe in detail the coupled RISM/simulation approach, and offer some error analysis. Finally we give the results of its application to a set of small test molecules in aqueous solution. The success shown in some of our results demonstrates that the coupled RISM/simulation approach is worth considering further as a potentially useful tool in studies of solvated systems, such as aqueous molecular biosystems.