Cooperative interaction of n-butylammonium ion with 1,3-alternate tetrapropoxycalix [4]arene: NMR and theoretical study

The interaction of 1,3-alternate tetrapropoxycalix[4]arene (1) with n-butylammonium ion (2) in CD(2)Cl(2) was examined using (1)H, (13)C and (14)N NMR spectroscopy and DFT (density functional theory) calculations. NMR shows that 1 forms with 2 an equimolecular hydrogen-bonded complex with the equilibrium constant 5.91 x 10(3) l/mol at 296 K. The structure of the complex can be shown to be asymmetric at 203 K, with 2 interacting by hydrogen bonds with the two ethereal oxygen atoms of one half of 1 and with the pi system of the other half, but is rapidly averaged to an apparent C(4h) symmetry by chemical exchange at higher temperatures. Using two related but independent techniques based on transverse and rotating-frame proton relaxation, it is shown that only an intermolecular exchange of 2 between the bound and free states takes place, in contrast to previously studied interaction of 1 with H(3)O(+). Its correlation time is 0.169 ms. It is shown by DFT calculations that such swift exchange is not possible without a cooperative interaction of both 2 and 1 with several molecules of water present. Similarities and contrasts between the exchange processes of 2 and H(3)O(+) bound to 1 are discussed, in particular with respect to the apparent quantum tunneling of the latter inside the molecule of the complex.     

Interaction of hydronium ion with dibenzo-18-crown-6: NMR, IR, and theoretical study

Interaction of dibenzo-18-crown-6 (DBC) with H 3O (+) (HP) in nitrobenzene- d 5 and dichloromethane- d 2 was studied by using (1)H and (13)C NMR spectra and relaxations, FTIR spectra, and quantum chemical DFT calculations. NMR shows that the DBC*HP complex is in a dynamic equilibrium with the reactants, the equilibrium constant K being 0.66 x 10 (3), 1.16 x 10 (4), and 1.03 x 10 (4) L x mol (-1) in CD 2Cl 2, nitrobenzene, and acetonitrile, respectively. The complex appears to have a C 2 v symmetry in NMR, but FTIR combined with DFT normal mode calculations suggest that such high symmetry is only apparent and due to exchange averaging of the structure. FTIR spectra as well as energy-optimized DFT calculations show that the most stable state of the complex in solution is that with three linear hydrogen bonds of HP with one CH 2-O-CH 2 and two Ar-O-Ar oxygen atoms. The structure is similar to that found in solid state but adopts a somewhat different conformation in solution. The dynamics of exchange between bound and free DBC was studied by NMR transverse relaxation. It was found to be too fast to give reproducible results when measured with the ordinary CPMG sequence or its variant DIFTRE removing residual static dipolar interaction, but it could be established by rotating-frame measurements with high intensity of the spin-lock field. The correlation time of exchange was found to be 5.6 x 10 (-6) and 3.8 x 10 (-6) s in dichloromethane and nitrobenzene, respectively. Such fast exchange can be explained by cooperative assistance of present water molecules.     

The entrance complex, transition state, and exit complex for the F + H2O → HF + OH reaction. Definitive predictions. Comparison with popular density functional methods

Following the H + H(2) and F + H(2) reactions, the fluorine atom - water system has the potential to become one of the best understood chemical reactions. Stationary points for the F + H(2)O potential energy surface have been located with the "Gold Standard" CCSD(T) method using the Dunning correlation consistent basis sets through quintuple zeta. The CCSD(T)/cc-pV5Z barrier height is prediced to be 2.5 kcal mol(-1), less than previous estimates of 4-7 kcal mol(-1). From higher level theoretical studies of the prototypical F + H(2) reaction, this barrier should be less than 0.5 kcal mol(-1) above the exact, nonrelativistic classical barrier height. 41 of the 49 DFT methods applied to F + H(2)O predict no barrier at all. The eight DFT methods that do predict a barrier show exothermicities that are somewhat too small. The CCSD(T)/cc-pV5Z entrance complex is bound by 3.4 kcal mol(-1) relative to separated F + H(2)O. The analogous exit complex is bound by 5.9 kcal mol(-1) relative to separated HF + OH.     

Ab initio/density functional theory and multichannel RRKM study for the ClO + CH2O reaction

The potential energy surface for the CH(2)O + ClO reaction was calculated at the QCISD(T)/6-311G(2d,2p)//B3LYP/6-311G(d,p) level of theory. The rate constants for the lower barrier reaction channels producing HOCl + HCO, H atom, OCH(2)OCl, cis-HC(O)OCl and trans-HC(O)OCl have been calculated by TST and multichannel RRKM theory. Over the temperature range of 200-2000 K, the overall rate constants were k(200-2000K) = 1.19 x 10(-13)T(0.79) exp(-3000.00/T). At 250 K, the calculated overall rate constant was 5.80 x 10(-17) cm(3) molecule(-1) s(-1), which was in good agreement with the experimental upper limit data. The calculated results demonstrated that the formation of HOCl + HCO was the dominant reaction channel and was exothermic by 9.7 kcal/mol with a barrier of 5.0 kcal/mol. When it retrograded to the reactants CH(2)O + ClO, an energy barrier of 14.7 kcal/mol is required. Furthermore, when HOCl decomposed into H + ClO, the energy required was 93.3 kcal/mol. These results suggest that the decomposition in both the forward and backward directions for HOCl would be difficult in the ground electronic state.     

The impact of tripodal chelates on water exchange kinetics and mechanisms: a variable temperature and pressure 17O NMR study to clarify the structure-reactivity relationship in [Fe(II)(L)(H2O2](-) complexes

The effect of temperature and pressure on the water exchange reaction of [Fe(II)(NTA)(H2O)2](-) and [Fe(II)(BADA)(H2O)2](-) (NTA = nitrilotriacetate; BADA = beta-alanindiacetate) was studied by 17O NMR spectroscopy. The [Fe(II)(NTA)(H2O)2](-) complex showed a water exchange rate constant, k(ex), of (3.1 +/- 0.4) x 10(6) s(-1) at 298.2 K and ambient pressure. The activation parameters DeltaH( not equal), DeltaS( not equal) and DeltaV( not equal) for the observed reaction are 43.4 +/- 2.6 kJ mol(-1), + 25 +/- 9 J K(-1) mol(-1) and + 13.2 +/- 0.6 cm(3) mol(-1), respectively. For [Fe(II)(BADA)(H2O)2](-), the water exchange reaction is faster than for the [Fe(II)(NTA)(H2O)2](-) complex with k(ex) = (7.4 +/- 0.4) x 10(6) s(-1) at 298.2 K and ambient pressure. The activation parameters DeltaH( not equal), DeltaS( not equal) and DeltaV( not equal) for the water exchange reaction are 40.3 +/- 2.5 kJ mol(-1), + 22 +/- 9 J K(-1) mol(-1) and + 13.3 +/- 0.8 cm(3) mol(-1), respectively. The effect of pressure on the exchange rate constant is large and very similar for both systems, and the numerical values for DeltaV( not equal) suggest in both cases a limiting dissociative (D) mechanism for the water exchange process.     

Structure and reactivity of bis(silyl) dihydride complexes (PMe(3))(3)Ru(SiR(3))(2)(H)(2): model compounds and real intermediates in a dehydrogenative C-Si bond forming reaction

A series of stable complexes, (PMe(3))(3)Ru(SiR(3))(2)(H)(2) ((SiR(3))(2) = (SiH(2)Ph)(2), 3a; (SiHPh(2))(2), 3b; (SiMe(2)CH(2)CH(2)SiMe(2)), 3c), has been synthesized by the reaction of hydridosilanes with (PMe(3))(3)Ru(SiMe(3))H(3) or (PMe(3))(4)Ru(SiMe(3))H. Compounds 3a and 3c adopt overall pentagonal bipyramidal geometries in solution and the solid state, with phosphine and silyl ligands defining trigonal bipyramids and ruthenium hydrides arranged in the equatorial plane. Compound 3a exhibits meridional phosphines, with both silyl ligands equatorial, whereas the constraints of the chelate in 3c result in both axial and equatorial silyl environments and facial phosphines. Although there is no evidence for agostic Si-H interactions in 3a and 3b, the equatorial silyl group in 3c is in close contact with one hydride (1.81(4) A) and is moderately close to the other hydride (2.15(3) A) in the solid state and solution (nu(Ru.H.Si) = 1740 cm(-)(1) and nu(RuH) = 1940 cm(-)(1)). The analogous bis(silyl) dihydride, (PMe(3))(3)Ru(SiMe(3))(2)(H)(2) (3d), is not stable at room temperature, but can be generated in situ at low temperature from the 16e(-) complex (PMe(3))(3)Ru(SiMe(3))H (1) and HSiMe(3). Complexes 3b and 3d have been characterized by multinuclear, variable temperature NMR and appear to be isostructural with 3a. All four complexes exhibit dynamic NMR spectra, but the slow exchange limit could not be observed for 3c. Treatment of 1 with HSiMe(3) at room temperature leads to formation of (PMe(3))(3)Ru(SiMe(2)CH(2)SiMe(3))H(3) (4b) via a CH functionalization process critical to catalytic dehydrocoupling of HSiMe(3) at higher temperatures. Closer inspection of this reaction between -110 and -10 degrees C by NMR reveals a plethora of silyl hydride phosphine complexes formed by ligand redistribution prior to CH activation. Above ca. 0 degrees C this mixture converts cleanly via silane dehydrogenation to the very stable tris(phosphine) trihydride carbosilyl complex 4b. The structure of 4b was determined crystallographically and exhibits a tetrahedral P(3)Si environment around the metal with the three hydrides adjacent to silicon and capping the P(2)Si faces. Although strong Si.HRu interactions are not indicated in the structure or by IR, the HSi distances (2.00(4) - 2.09(4) A) and average coupling constant (J(SiH) = 25 Hz) suggest some degree of nonclassical SiH bonding in the RuH(3)Si moiety. The least hindered complex, 3a, reacts with carbon monoxide principally via an H(2) elimination pathway to yield mer-(PMe(3))(3)(CO)Ru(SiH(2)Ph)(2), with SiH elimination as a minor process. However, only SiH elimination and formation of (PMe(3))(3)(CO)Ru(SiR(3))H is observed for 3b-d. The most hindered bis(silyl) complex, 3d, is extremely labile and even in the absence of CO undergoes SiH reductive elimination to generate the 16e(-) species 1 (DeltaH(SiH)(-)(elim) = 11.0 +/- 0.6 kcal x mol(-)(1) and DeltaS(SiH)(-)(elim) = 40 +/- 2 cal x mol(-)(1) x K(-)(1); Delta = 9.2 +/- 0.8 kcal x mol(-)(1) and Delta = 9 +/- 3 cal x mol(-)(1).K(-)(1)). The minimum barrier for the H(2) reductive elimination can be estimated, and is higher than that for silane elimination at temperatures above ca. -50 degrees C. The thermodynamic preferences for oxidative additions to 1 are dominated by entropy contributions and steric effects. Addition of H(2) is by far most favorable, whereas the relative aptitudes for intramolecular silyl CH activation and intermolecular SiH addition are strongly dependent on temperature (DeltaH(SiH)(-)(add) = -11.0 +/- 0.6 kcal x mol(-)(1) and DeltaS(SiH)(-)(add) = -40 +/- 2 cal.mol(-)(1) x K(-)(1); DeltaH(beta)(-CH)(-)(add) = -2.7 +/- 0.3 kcal x mol(-)(1) and DeltaS(beta)(-CH)(-)(add) = -6 +/- 1 cal x mol(-)(1) x K(-)(1)). Kinetic preferences for oxidative additions to 1 - intermolecular SiH and intramolecular CH - have been also quantified: Delta = -1.8 +/- 0.8 kcal x mol(-)(1) and Delta = -31 +/- 3 cal x mol(-)(1).K(-)(1); Delta = 16.4 +/- 0.6 kcal x mol(-)(1) and Delta = -13 +/- 6 cal x mol(-)(1).K(-)(1). The relative enthalpies of activation (-)(1) x K(-)(1)). Kinetic preferences for oxidative additions to 1 - intermolecular SiH and intramolecular CH - have been also quantified: Delta (H)SiH(add) = 1.8 +/- 0.8 kcal x mol(-)(1) and Delta S((SiH-add) =31+/- 3 cal x mol(-)(1) x K(-)(1); Delta S (SiH -add) = 16.4 +/- 0.6 kcal x mol(-)(1) and =Delta S (SiH -CH -add) =13+/- 6 cal x mol(-)(1) x K(-)(1). The relative enthalpies of activation are interpreted in terms of strong SiH sigma-complex formation - and much weaker CH coordination - in the transition state for oxidative addition.     

Methane-to-methanol oxidation by the hydrated iron(IV) oxo species in aqueous solution: a combined DFT and car-parrinello molecular dynamics study

Previously, we have shown that the ferryl ion ([FeIVO]2+) is easily produced from Fenton's reagent (i.e., a mixture of Fe2+ ions and H2O2 in aqueous solution), using DFT and Car-Parrinello MD calculations. To verify that the ferryl ion can indeed act as the active species in oxidation reactions with Fenton's reagent, we study in the present paper the reactivity of the ferryl ion toward an organic substrate, in particular the oxidation of methane to methanol. In the first part of this paper, we perform static DFT calculations on the reaction of CH4 with the [(H2O)5FeIVO]2+ complex in vacuo that show a strong prevalence of the oxygen-rebound mechanism over the methane coordination mechanism. This is in agreement with the static DFT results for methane oxidation by biocatalysts MMO and P450, but not with those for methane oxidation by bare metal-oxo ions, where the methane coordination mechanism prevails. The highest energy barrier in the oxygen-rebound mechanism is only 3 kcal/mol in vacuo, whereas in the methane coordination mechanism the highest barrier is 23 kcal/mol. Overall the oxidation reaction energy is downhill by 47 kcal/mol. We conclude that the ferryl ion can indeed act as the oxidative intermediate in the Fenton oxidation of organic species. In the second part of this paper, we perform a preliminary assessment of solvent effects on the oxidation by the ferryl ion in aqueous solution using the method of constrained (first principles) molecular dynamics. The free energy barrier of the H-abstraction reaction from methane by the ferryl ion (i.e., the first step in the rebound mechanism) in aqueous solution is, with 22 kcal/mol in solution, significantly higher than in vacuo. Given the fact that methane has a relatively strong C-H bond (ca. 10 kcal/mol stronger than the C-H bonds in the more typical Fenton's reagent substrates), we infer that for many organic substrates oxidation with the ferryl ion as an active intermediate may be a perfectly viable route.     

NMR and theoretical study of the cooperative interaction of hydrated proton with dibenzo-24-crown-8

Interaction of hydrated protons (HPs) with dibenzo-24-crown-8 (DBC in nitrobenzene-d(5) was studied by (1)H and (13)C NMR under assistance of ab initio-density functional theory (DFT) quantum calculations. HPs were afforded by hydrogen bis(1,2-dicarbollyl) cobaltate (HDCC) with 3.5 M excess of H(2)O. The forming of a complex between HP and DBC leads to marked and additive relative shifts of both (1)H and (13)C signals. This was utilized for the estimation of the stabilization constant K of the complex. Its value is at least 10(6) l/mol, which agrees with the result of independent extraction method (log K = 6.3). Using absolute integral intensities of the HP signal in a water-saturated system, it was shown that the form of HP present in the complex must be H(5)O(2)(+), in accord with formerly published structure of the complex in crystalline form. The investigation of the dynamics of exchange between bound and free DBC by transverse relaxation using variably spaced pulses in the Carr-Purcell-Meiboom-Gill (CPMG) sequence or on-resonance rotating-frame relaxation with variable spin-lock field intensity was partly hampered by the fast relaxation of some signals in the complex because of relative immobilization of its internal motions. In order to remove these effects, a pulse sequence dipolar interaction-free transverse relaxation (DIFTRE) for static DIFTRE was devised and the MLEV17 pulse sequence with high intensity of electromagnetic field was used in a separate series of experiments. Using the results of these latter experiments, the correlation time of exchange was established to be about 0.8 ms, which complied with the shape of the spectra. The accompanying ab initio DFT calculations showed that the apparent symmetry of the molecules of both DBC and its complex with H(5)O(2)(+) was probably the result of averaging of energetically close conformations (five for DBC and four for the complex). Both NMR and the calculations show that the basic mode of binding of the ion in the complex is analogous to that found in crystal but the most pronounced conformation is slightly different.     

First investigation of non-classical dihydrogen bonding between an early transition-metal hydride and alcohols: IR, NMR, and DFT approach

The interaction of [NbCp(2)H(3)] with fluorinated alcohols to give dihydrogen-bonded complexes was studied by a combination of IR, NMR and DFT methods. IR spectra were examined in the range from 200-295 K, affording a clear picture of dihydrogen-bond formation when [NbCp(2)H(3)]/HOR(f) mixtures (HOR(f) = hexafluoroisopropanol (HFIP) or perfluoro-tert-butanol (PFTB)) were quickly cooled to 200 K. Through examination of the OH region, the dihydrogen-bond energetics were determined to be 4.5+/-0.3 kcal mol(-1) for TFE (TFE = trifluoroethanol) and 5.7+/-0.3 kcal mol(-1) for HFIP. (1)H NMR studies of solutions of [NbCp(2)H(2)(B)H(A)] and HFIP in [D(8)]toluene revealed high-field shifts of the hydrides H(A) and H(B), characteristic of dihydrogen-bond formation, upon addition of alcohol. The magnitude of signal shifts and T(1) relaxation time measurements show preferential coordination of the alcohol to the central hydride H(A), but are also consistent with a bifurcated character of the dihydrogen bonding. Estimations of hydride-proton distances based on T(1) data are in good accord with the results of DFT calculations. DFT calculations for the interaction of [NbCp(2)H(3)] with a series of non-fluorinated (MeOH, CH(3)COOH) and fluorinated (CF(3)OH, TFE, HFIP, PFTB and CF(3)COOH) proton donors of different strengths showed dihydrogen-bond formation, with binding energies ranging from -5.7 to -12.3 kcal mol(-1), depending on the proton donor strength. Coordination of proton donors occurs both to the central and to the lateral hydrides of [NbCp(2)H(3)], the former interaction being of bifurcated type and energetically slightly more favourable. In the case of the strong acid H(3)O(+), the proton transfer occurs without any barrier, and no dihydrogen-bonded intermediates are found. Proton transfer to [NbCp(2)H(3)] gives bis(dihydrogen) [NbCp(2)(eta(2)-H(2))(2)](+) and dihydride(dihydrogen) complexes [NbCp(2)(H)(2)(eta(2)-H(2))](+) (with lateral hydrides and central dihydrogen), the former product being slightly more stable. When two molecules of TFA were included in the calculations, in addition to the dihydrogen-bonded adduct, an ionic pair formed by the cationic bis(dihydrogen) complex [NbCp(2)(eta(2)-H(2))(2)](+) and the homoconjugated anion pair (CF(3)COO...H...OOCCF(3))(-) was found as a minimum. It is very likely that these ionic pairs may be intermediates in the H/D exchange between the hydride ligands and the OD group observed with the more acidic alcohols in the NMR studies.     

Conformational analysis: crystallographic,NMR, and molecular mechanics studies of flexible sulfonic esters

Two novel X-ray structures of the sulfonic ester derivatives 2-(6-iodo-1,3-benzodioxol-5-yl)ethyl 4-nitrobenzenesulfonate, 3, and 2-(6-iodo-1,3-benzodioxol-5-yl)ethyl 4-methylbenzenesulfonate, 4, have been obtained in a study aimed at analyzing the structures and conformations of sulfonic ester derivatives that are routinely used in alkaloid syntheses. The crystal structure of 4 is highly unusual, containing four independent molecules that belong to two distinct conformational types: (1) a hairpin conformation (stabilized mainly by intramolecular pi-stacking) and (2) a stepped conformation (stabilized mainly by intermolecular pi-stacking). Compound 3, on the other hand, crystallizes exclusively as the hairpin conformer. New MM+ force field parameters for sulfonic esters have been developed using the X-ray data, empirical rules, and DFT calculations to estimate the bond dipole parameters. Grid searches of conformational space for 3 and 4 using MM methods show that there are several gas-phase conformations within 5 kcal/mol of the global minimum and that the lowest energy conformations (by approximately 4.6 kcal/mol) are of the hairpin type. Analysis of the MM conformational energies suggests that the dominant intramolecular interaction stabilizing the hairpin conformations of 3 and 4 is van der Waals attraction. Moreover, the lattice energies for packing the hairpin conformations of 3 and 4 are approximately 4 kcal/mol more favorable than for the stepped conformations. Various intermolecular interactions contribute to the complexity of the observed crystal structures of 3 and 4, including electrostatic attraction between O and I atoms in neighboring molecules. Langevin dynamics (LD) simulations at several temperatures (6.0 ns, friction coefficient = 2.5 ps(-1)) indicate that the conformational exchange rates are approximately 10(10)-10(11) s(-1) over the temperature range 213-400 K, accounting for the temperature-independent (1)H NMR spectra of 3 and 4.     

Accurate computational determination of the binding energy of the SO3 x H2O complex

Reliable thermochemical data for the reaction SO3 + H2O<-->SO3 x H2O (1a) are of crucial importance for an adequate modeling of the homogeneous H2SO4 formation in the atmosphere. We report on high-level quantum chemical calculations to predict the binding energy of the SO3 x H2O complex. The electronic binding energy is accurately computed to De = 40.9+/-1.0 kJ/mol = 9.8+/-0.2 kcal/mol. By using harmonic frequencies from density functional theory calculations (B3LYP/cc-pVTZ and TPSS/def2-TZVP), zero-point and thermal energies were calculated. From these data, we estimate D0 = -Delta H(1a)0(0 K) = 7.7+/-0.5 kcal/mol and Delta H(1a)0(298 K) = -8.3+/-1.0 kcal/mol.     

New insight into the gas-phase bimolecular self-reaction of the HOO radical

The singlet and triplet potential energy surfaces (PESs) for the gas-phase bimolecular self-reaction of HOO*, a key reaction in atmospheric environments, have been investigated by means of quantum-mechanical electronic structure methods (CASSCF and CASPT2). All the reaction pathways on both PESs consist of a first step involving the barrierless formation of a prereactive doubly hydrogen-bonded complex, which is a diradical species lying about 8 kcal/mol below the energy of the reactants at 0 K. The lowest energy reaction pathway on both PESs is the degenerate double hydrogen exchange between the HOO* moieties of the prereactive complex via a double proton transfer mechanism involving an energy barrier of only 1.1 kcal/mol for the singlet and 3.3 kcal/mol for the triplet at 0 K. The single H-atom transfer between the two HOO* moieties of the prereactive complex (yielding HOOH + O2) through a pathway keeping a planar arrangement of the six atoms involves a conical intersection between either two singlet or two triplet states of A' and A" symmetries. Thus, the lowest energy reaction pathway occurs via a nonplanar cisoid transition structure with an energy barrier of 5.8 kcal/mol for the triplet and 17.5 kcal/mol for the singlet at 0 K. The simple addition between the terminal oxygen atoms of the two HOO* moieties of the prereactive complex, leading to the straight chain H2O4 intermediate on the singlet PES, involves an energy barrier of 7.3 kcal/mol at 0 K. Because the decomposition of such an intermediate into HOOH + O2 entails an energy barrier of 45.2 kcal/mol at 0 K, it is concluded that the single H-atom transfer on the triplet PES is the dominant pathway leading to HOOH + O2. Finally, the strong negative temperature dependence of the rate constant observed for this reaction is attributed to the reversible formation of the prereactive complex in the entrance channel rather than to a short-lived tetraoxide intermediate.     

Favorable pendant-amino metal chelation in VX nerve agent model systems

We have performed DFT computational studies [B3LYP, 6-31+G] to obtain metal ion coordination isomers of VX-Me [MeP(O)(OMe)(SCH2CH2NMe2)], a model of two of the most lethal nerve agents: VX [MeP(O)(OEt)(SCH2CH2N(iPr)2)] and Russian-VX [MeP(O)(OCH2CHMe2)(SCH2CH2N(Et)2)]. Our calculations involved geometry optimizations of the neutral VX-Me model as well as complexes with H+, Li+, Na+, K+, Be2+, Mg2+, and Ca2+ that yielded 2-8 different stable chelation modes for each ion that involved mainly mono- and bidentate binding. Importantly, our studies revealed that the [O(P),N] bidentate binding mode, long thought to be the active mode in differentiating the hydrolytic path of VX from other nerve agents, was the most stable for all ions studied here. Binding energy depended mainly on ionic size as well as charge, with binding energies ranging from 364 kcal mol(-1) for Be2+ to 33 kcal mol(-1) for K+. Furthermore, calculated NMR shifts for VX-Me correlate to experimental values of VX.     

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.     

The effects of water vapor on the CH3O2 self-reaction and reaction with HO2

The gas phase reactions of CH3O2 + CH3O2, HO2 + HO2, and CH3O2 + HO2 in the presence of water vapor have been studied at temperatures between 263 and 303 K using laser flash photolysis coupled with UV time-resolved absorption detection at 220 and 260 nm. Water vapor concentrations were quantified using tunable diode laser spectroscopy operating in the mid-IR. The HO2 self-reaction rate constant is significantly enhanced by water vapor, consistent with what others have reported, whereas the CH3O2 self-reaction and the cross-reaction (CH3O2 + HO2) rate constants are nearly unaffected. The enhancement in the HO2 self-reaction rate coefficient occurs because of the formation of a strongly bound (6.9 kcal mol(-1)) HO2 x H2O complex during the reaction mechanism where the H2O acts as an energy chaperone. The nominal impact of water vapor on the CH3O2 self-reaction rate coefficient is consistent with recent high level ab initio calculations that predict a weakly bound CH3O2 x H2O complex (2.3 kcal mol(-1)). The smaller binding energy of the CH3O2 x H2O complex does not favor its formation and consequent participation in the methyl peroxy self-reaction mechanism.     

A 3-fold "butterfly valve" in command of the encapsulation's kinetic stability. Molecular baskets at work

Molecular basket 1, composed of a semirigid tris-norbornadiene framework and three revolving pyridine-based gates at the rim, has been built to "dynamically" enclose space and as such regulate molecular encapsulation. The gates were shown to fold via intramolecular hydrogen bonding and thereby form a C3nu symmetrical receptor: the 1H NMR resonance for the amide N-H protons of the pyridine gates appeared downfield (delta= 10.98 ppm), and the N-H vibrational stretch (IR) was observed at 3176 cm(-1). Accordingly, density functional theory (DFT, B3LYP) investigations revealed for the closed conformers of 1 to be energetically the most stable and dominant. The gearing of the pyridine "gates", about their axis, led to the interconversion of two dynamic enantiomers 1A and 1B comprising the clockwise and counterclockwise seam of intramolecular hydrogen bonds. Dynamic 1H NMR spectroscopic measurements and line-shape simulations suggested that the energy barrier of 10.0 kcal/mol (DeltaG++(A/B), 298 K) is required for the 1A/B interconversion, when CCl4 occupies the cavity of 1. Likewise, the activation free energy for CCl4 departing the basket was found to be 13.1 kcal/mol (DeltaG++, 298 K), whereas the thermodynamic stability of 1:CCl4 complex was -2.7 kcal/mol (DeltaGdegrees, 298 K). In view of that, CCl4 (but also (CH3)3CBr) was proposed to escape from, and a molecule of solvent to enter, the basket when the gates rotate about their axis: the exit of CCl4 requires the activation energy of 12.7 kcal/mol (DeltaG++(A/B) + DeltaGdegrees), similar to the experimentally found 13.1 kcal/mol (DeltaG++).     

Triggering water exchange mechanisms via chelate architecture. Shielding of transition metal centers by aminopolycarboxylate spectator ligands

Paramagnetic effects on the relaxation rate and shift difference of the (17)O nucleus of bulk water enable the study of water exchange mechanisms on transition metal complexes by variable temperature and variable pressure NMR. The water exchange kinetics of [Mn(II)(edta)(H2O)](2-) (CN 7, hexacoordinated edta) was reinvestigated and complemented by variable pressure NMR data. The results revealed a rapid water exchange reaction for the [Mn(II)(edta)(H2O)](2-) complex with a rate constant of k(ex) = (4.1 +/- 0.4) x 10(8) s(-1) at 298.2 K and ambient pressure. The activation parameters DeltaH(double dagger), DeltaS(double dagger), and DeltaV(double dagger) are 36.6 +/- 0.8 kJ mol(-1), +43 +/- 3 J K(-1) mol(-1), and +3.4 +/- 0.2 cm(3) mol(-1), which are in line with a dissociatively activated interchange (I(d)) mechanism. To analyze the structural influence of the chelate, the investigation was complemented by studies on complexes of the edta-related tmdta (trimethylenediaminetetraacetate) chelate. The kinetic parameters for [Fe(II)(tmdta)(H2O)](2-) are k(ex) = (5.5 +/- 0.5) x 10(6) s(-1) at 298.2 K, DeltaH(double dagger) = 43 +/- 3 kJ mol(-1), DeltaS(double dagger) = +30 +/- 13 J K(-1) mol(-1), and DeltaV(double dagger) = +15.7 +/- 1.5 cm(3) mol(-1), and those for [Mn(II)(tmdta)(H2O)](2-) are k(ex) = (1.3 +/- 0.1) x 10(8) s(-1) at 298.2 K, DeltaH(double dagger) = 37.2 +/- 0.8 kJ mol(-1), DeltaS(double dagger) = +35 +/- 3 J K(-1) mol(-1), and DeltaV(double dagger) = +8.7 +/- 0.6 cm(3) mol(-1). The water containing species, [Fe(III)(tmdta)(H2O)](-) with a fraction of 0.2, is in equilibrium with the water-free hexa-coordinate form, [Fe(III)(tmdta)](-). The kinetic parameters for [Fe(III)(tmdta)(H2O)](-) are k(ex) = (1.9 +/- 0.8) x 10(7) s(-1) at 298.2 K, DeltaH(double dagger) = 42 +/- 3 kJ mol(-1), DeltaS(double dagger) = +36 +/- 10 J K(-1) mol(-1), and DeltaV(double dagger) = +7.2 +/- 2.7 cm(3) mol(-1). The data for the mentioned tmdta complexes indicate a dissociatively activated exchange mechanism in all cases with a clear relationship between the sterical hindrance that arises from the ligand architecture and mechanistic details of the exchange process for seven-coordinate complexes. The unexpected kinetic and mechanistic behavior of [Ni(II)(edta')(H2O)](2-) and [Ni(II)(tmdta')(H2O)](2-) is accounted for in terms of the different coordination number due to the strong preference for an octahedral coordination environment and thus a coordination equilibrium between the water-free, hexadentate [M(L)](n+) and the aqua-pentadentate forms [M(L')(H2O)](n+) of the Ni(II)-edta complex, which was studied in detail by variable temperature and pressure UV-vis experiments. For [Ni(II)(edta')(H2O)](2-) (CN 6, pentacoordinated edta) a water substitution rate constant of (2.6 +/- 0.2) x 10(5) s(-1) at 298.2 K and ambient pressure was measured, and the activation parameters DeltaH(double dagger), DeltaS(double dagger), and DeltaV(double dagger) were found to be 34 +/- 1 kJ mol(-1), -27 +/- 2 J K(-1) mol(-1), and +1.8 +/- 0.1 cm(3) mol(-1), respectively. For [Ni(II)(tmdta')(H2O)](2-), we found k = (6.4 +/- 1.4) x 10(5) s(-1) at 298.2 K, DeltaH(double dagger) = 22 +/- 4 kJ mol(-1), and DeltaS(double dagger) = -59 +/- 5 J K(-1) mol(-1). The process is referred to as a water substitution instead of a water exchange reaction, since these observations refer to the intramolecular displacement of coordinated water by the carboxylate moiety in a ring-closure reaction.     

Active Thermochemical Tables: accurate enthalpy of formation of hydroperoxyl radical, HO2

Through the use of the Active Thermochemical Tables approach, the best currently available enthalpy of formation of HO2 has been obtained as delta(f)H(o)298 (HO2) = 2.94 +/- 0.06 kcal mol(-1) (3.64 +/- 0.06 kcal mol(-1) at 0 K). The related enthalpy of formation of the positive ion, HO2+, within the stationary electron convention is delta(f)H(o)298 (HO2+) = 264.71 +/- 0.14 kcal mol(-1) (265.41 +/- 0.14 kcal mol(-1) at 0 K), while that for the negative ion, HO2- (within the same convention), is delta(f)H(o)298 (HO2-) = -21.86 +/- 0.11 kcal mol(-1) (-21.22 +/- 0.11 kcal mol(-1) at 0 K). The related proton affinity of molecular oxygen is PA298(O2) = 100.98 +/- 0.14 kcal mol(-1) (99.81 +/- 0.14 kcal mol(-1) at 0 K), while the gas-phase acidity of H2O2 is delta(acid)G(o)298 (H2O2) = 369.08 +/- 0.11 kcal mol(-1), with the corresponding enthalpy of deprotonation of H2O2 of delta(acid)H(o)298 (H2O2) = 376.27 +/- 0.11 kcal mol(-1) (375.02 +/- 0.11 kcal mol(-1) at 0 K). In addition, a further improved enthalpy of formation of OH is briefly outlined, delta(f)H(o)298 (OH) = 8.93 +/- 0.03 kcal mol(-1) (8.87 +/- 0.03 kcal mol(-1) at 0 K), together with new and more accurate enthalpies of formation of NO, delta(f)H(o)298 (NO) = 21.76 +/- 0.02 kcal mol(-1) (21.64 +/- 0.02 kcal mol(-1) at 0 K) and NO2, delta(f)H(o)298 (NO2) = 8.12 +/- 0.02 kcal mol(-1) (8.79 +/- 0.02 kcal mol(-1) at 0 K), as well as H(2)O(2) in the gas phase, delta(f)H(o)298 (H2O2) = -32.45 +/- 0.04 kcal mol(-1) (-31.01 +/- 0.04 kcal mol(-1) at 0 K). The new thermochemistry of HO2, together with other arguments given in the present work, suggests that the previous equilibrium constant for NO + HO2 --> OH + NO2 was underestimated by a factor of approximately 2, implicating that the OH + NO2 rate was overestimated by the same factor. This point is experimentally explored in the companion paper of Srinivasan et al. (next paper in this issue).     

Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

Computational evaluation of the evidence for tri-trans-[12]annulene

[reaction: see text] Automerization in tri-trans-[12]annulene (1) was investigated by DFT, MP2, and coupled-cluster methods. Using the highest level of theory employed here, CCSD(T)/cc-pVDZ//BHandHLYP/6-311+G(d,p), we located two low-energy pathways for degenerate conformational change from the lowest-energy conformer of 1 (1a): one with E(a) = 4.5 kcal/mol that interconverts the three inner trans hydrogens with the three outer trans hydrogens and one with E(a) = 2.7 kcal/mol that interconverts the three inner hydrogens with each other. These results are consistent with the experimental results of Oth and co-workers on [12]annulene 1a (Oth, J. F. M.; R?ttele, H.; Schr?der, G. Tetrahedron Lett. 1970, 61). The conformational exchange of the inner trans hydrogens with the outer ones is predicted to occur via a one-step process involving a C(2)-symmetric transition state and not via the D(3)-symmetric transition state (1b) that was postulated earlier. Conformer 1b was found to be a shallow minimum 6.7 kcal/mol above 1a with a barrier of 0.4 kcal/mol for conversion to 1a. Finally, GIAO-B3LYP/6-311+G(d,p) and BHandHLYP/6-311+G(d,p) computed (1)H NMR chemical shifts of 1a and three other low-lying isomers support Oth's original assignment of observed (1)H NMR peaks to 1a at both low and high temperature.