Does Silica Surface Catalyse Peptide Bond Formation? New Insights from First‐Principles Calculations

The role that silica surface could have played in prebiotic chemistry as a catalyst for peptide bond formation has been addressed at the B3LYP/6-31+G(d,p) level for a model reaction involving glycine and ammonia on a silica cluster mimicking an isolated terminal silanol group present at the silica surface. Hydrogen-bond complexation between glycine and the silanol is followed by the formation of the mixed surface anhydride Si(surf)-O-C(=O)-R, which has been suggested in the literature to activate the C=O bond towards nucleophilic attack by a second glycine molecule, here simulated by the simpler NH3 molecule. However, B3LYP/6-31+G(d,p) calculations show that formation of the surface mixed anhydride Si(surf)-O-C(=O)-R is disfavoured (delta(r)G298 approximately 6 kcal mol(-1)), and that the surface bond only moderately lowers the free-energy barrier of the nucleophilic attack responsible for peptide bond formation (deltaG298(double dagger) approximately 48 kcal mol(-1)) in comparison with the uncatalysed reaction (deltaG298(double dagger) approximately 52 kcal mol(-1)). A further decrease of the free-energy barrier of peptide bond formation (deltaG298(double dagger) approximately 41 kcal mol(-1)) is achieved by a single water molecule close to the reaction centre acting as a proton-transfer helper in the activated complex. A possible role of strained silica surface defects on the formation of the surface mixed anhydride Si(surf)-O-C(=O)-R has also been addressed.     

Static and dynamic stereochemistry of the conformational atropisomers of tetra(o-tolyl)benzene

[graph: see text] Whereas only one atropisomer of 1,2,4,5-tetra(o-tolyl)benzene was observed by X-ray diffraction in the solid, five conformational atropisomers were detected by low-temperature NMR in solution. Their structures were assigned by a combination of NOE experiments, solvent effect, and ab initio calculations. Variable temperature dynamic NMR and bidimensional EXSY experiments allowed the barrier for the interconversion of these atropisomers to be determined (deltaG(double dagger) = 15.3 kcal mol(-1)).     

Dimers of and tautomerism between 2-pyrimidinethiol and 2(1H)-pyrimidinethione: a density functional theory (DFT) study

Density functional theory (BLYP, B3LYP, B3P86, B3PW91) with the 6-31+G(d,p), 6-311+G(d,p), and cc-pVTZ basis sets has been used to calculate structural parameters, relative energies, and vibrational spectra of 2-pyrimidinethiol (1) and 2(1H)-pyrimidinethione (2) and their hydrogen-bonded homodimers (C(2) 3, C(2h) [4](double dagger), C(2h) 5), monohydrates, and dihydrates and a heterodimer (6). Several transition state structures proposed for the tautomerization process have also been examined. At the B3PW91/6-311+G(d,p)//B3PW91/6-31+G(d,p) level of theory 2-pyrimidinethiol (1) is predicted to be 3.41 kcal/mol more stable (E(rel)) than 2(1H)-pyrimidinethione (2) in the gas phase and 2 is predicted to be 6.47 kcal/mol more stable than 1 in aqueous medium. An unfavorable planar intramolecular strained four center transition state (TS1) for the tautomerization of 1 and 2 in the gas-phase lies 29.07 kcal/mol higher in energy than 2-pyrimidinethiol (1). The C(2) 2-pyrimidinethiol dimer (3) is 6.84 kcal/mol lower in energy than the C(2) homodimer transition state structure ([11](double dagger)) that connects dimers 3 and 4. Transition state [11](double dagger) provides a facile pathway for tautomerization between 1 and 2 in the gas phase (monomer-dimer promoted tautomerization). The hydrogen bonded 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O structures are predicted to be 1.27 and 1.55 kcal/mol, respectively, higher in energy than 2(1H)-pyrimidinethione- - -H(2)O and 2(1H)-pyrimidinethione- - -2H(2)O. Water promoted tautomerization via cyclic transition states involving one water molecule (TS- - -H(2)O, [12](double dagger)) and two water molecules (TS- - -2H(2)O, [13](double dagger)) lie 11.42 and 11.44 kcal/mol, respectively, higher in energy than 2-pyrimidinethiol- - -H(2)O and 2-pyrimidinethiol- - -2H(2)O. Thus, the hydrated transition states [12](double dagger) and [13](double dagger) are involved in the tautomerism between 1 and 2 in aqueous medium.     

Kinetics, mechanism, and computational studies of rhenium-catalyzed desulfurization reactions of thiiranes (thioepoxides)

The oxorhenium(V) dimer {MeReO(edt)}2 (1; where edt = 1,2-ethanedithiolate) catalyzes S atom transfer from thiiranes to triarylphosphines and triarylarsines. Despite the fact that phosphines are more nucleophilic than arsines, phosphines are less effective because they rapidly convert the dimer catalyst to the much less reactive catalyst [MeReO(edt)(PAr3)] (2). With AsAr3, which does not yield the monomer, the rate law is given by v = k[thiirane][1], independent of the arsine concentration. The values of k at 25.0 degrees C in CDCl3 are 5.58 +/- 0.08 L mol(-1) s(-1) for cyclohexene sulfide and ca. 2 L mol(-1) s(-1) for propylene sulfide. The activation parameters for cyclohexene sulfide are deltaH(double dagger) = 10.0 +/- 0.9 kcal mol(-1) and deltaS(double dagger) = -21 +/- 3 cal K(-1) mol(-1). Arsine enters the catalytic cycle after the rate-controlling release of alkene, undergoing a reaction with the Re(VII)(O)(S) intermediate that is so rapid in comparison that it cannot be studied directly. The use of a kinetic competition method provided relative rate constants and a Hammett reaction constant, rho = -1.0. Computations showed that there is little thermodynamic selectivity for arsine attack at O or S of the intermediate. There is, however, a large kinetic selectivity in favor of Ar3AsS formation: the calculated values of deltaH(double dagger) for attack of AsAr3 at Re=O vs Re=S in Re(VII)(O)(S) are 23.2 and 1.1 kcal mol(-1), respectively.     

Retro Diels-Alder reactions of 5,6-disubstituted-7-oxabicyclo

Several 5,6-disubstituted-7-oxabicyclo[2.2.1]hept-2-enes (1-4) were synthesized on > or = 0.1 mol scale. The heat-induced retro Diels-Alder (rDA) decomposition of these derivatives was studied by thermal analysis, and the kinetics of the rDA were measured for 4. First-order rate constants (k = 1.91-14.2 x 10(-5) s(-1)), measured at four temperatures between 124 and 150 degrees C, were used to calculate Arrhenius activation parameters Ea (34.5 +/- 0.5 kcal/mol) and ln A (1.77 +/- 0.03 x 10(4)). The observed activation energy was significantly larger (by 9.5 kcal/mol) than that previously measured for the maleic anhydride adduct 1, and this was attributed to the difference in LUMO energies for the two dienophiles. Modeling of the activation parameters found for 4 with density functional theory (DFT) calculations for similar compounds 5 and 6 gave close quantitative correlations for deltaH double dagger, deltaG double dagger, deltaS double dagger. The rDA reactions studied were found to be entropy-driven.     

An ab-initio study of some homolytic substitution reactions of acyl radicals at silicon, germanium and tint

Ab initio calculations using the 6-311G**, cc-pVDZ, and (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional (BHandHLYP, B3LYP) calculations predict that homolytic substitution reactions of acetyl radicals at the silicon atoms in disilane can proceed via both backside and frontside attack mechanisms. At the highest level of theory (CCSD(T)/cc-pVDZ//MP2/cc-pVDZ), energy barriers (deltaE double dagger) of 77.2 and 81.9 kJ mol(-1) are calculated for the backside and frontside reactions respectively. Similar results are obtained for reactions involving germanium and tin with energy barriers (deltaE double dagger) of 53.7-84.2, and 55.0-89.7 kJ mol(-1) for the backside and frontside mechanisms, respectively. These data suggest that both homolytic substitution mechanisms are feasible for homolytic substitution reactions of acetyl radicals at silicon, germanium, and tin. BHandHLYP calculations provide geometries and energy barriers for backside and frontside transition states in good agreement with those obtained by traditional ab initio techniques.     

The mechanism of water oxidation catalysis promoted by [tpyRu(IV)=O]2L3+: a computational study

The resting state of the recently reported water oxidation catalyst [tpyRu(II)-OH(2)](2)L(3+) (tpy = terpyridine; L = bipyridylpyrazolylic anion) ([2,2](3+)) must be activated by a series of proton-coupled oxidations in which four protons and four electrons are removed overall to afford the catalytically competent species [tpyRu(IV)O](2)L(3+) ([4,4](3+)). We have examined all of the plausible redox intermediates utilizing density functional theory coupled to a continuum solvation model. Our calculations reproduce well the first three redox potentials under pH = 1 conditions, and a reasonable correlation between theory and experiment is found for the fourth irreversible redox process that accompanies O(2) generation. The computed oxidation potentials to access [5,4](4+) and [5,5](5+), 1.875 and 2.032 V vs NHE, respectively, exclude the otherwise plausible possibilities of the catalytically active species having a higher oxidation state. [4,4](3+) has an antiferromagnetically coupled ground state in which one ruthenium has two unpaired electrons antiparallel to those of the other ruthenium. As we found in our previous work, two radicaloid terminal oxygen moieties with different spin orientations that are induced by spin polarization from the electron-deficient Ru(IV) centers are found. Two mechanistic scenarios are relevant and interesting for the key O-O bond formation event: intramolecular oxo-oxo coupling and coupling between one terminal oxo and the oxygen atom of the incoming water substrate. The intramolecular oxo-oxo coupling is facile, with a low barrier of 13.9 kcal mol(-1), yielding a peroxo intermediate. The necessary subsequent addition of water in an associative substitution mechanism to cleave one of the Ru-peroxo bonds, however, is found to be impractical at room temperature, with a barrier of DeltaG(double dagger) = 30.9 kcal mol(-1). Thus, while plausible, the intramolecular oxo-oxo coupling is unproductive for generating molecular dioxygen. The intermolecular O-O coupling is associated with a high barrier (DeltaG(double dagger) = 40.2 kcal mol(-1)) and requires the assistance of an additional proton, which lowers the barrier dramatically to 24.5 kcal mol(-1).     

Cooperative and competitive effects of substituents at C1 and C4 on the barriers to ring inversion of 5,5-difluorobicyclo[2.1.0]pentanes

To identify the reasons for the very low barrier that has been measured for ring inversion of 1,4,5,5-tetrafluorobicyclo[2.1.0]pentane (deltaG(double dagger) = 6.8 +/- 0.2 kcal/mol), CASSCF and CASPT2 calculations have been performed on ring inversion in this and other bicyclo[2.1.0]pentanes. The results of the calculations show that a cooperative interaction between the geminal fluorines at C2 and the fluorines at C1 and C3 in the singlet cyclopentane-1,3-diyl transition structure (TS) contributes 3.7 kcal/mol to lowering the barrier to ring inversion in the tetrafluoro compound. In contrast, a competitive substituent effect in the TS for ring inversion of 1,4-dicyano-5,5-difluorobicyclo[2.1.0]pentane is predicted to raise the barrier height by 6.1 kcal/mol. The origin of these cooperative and competitive substituent effects is discussed.     

Alkaline hydrolysis of ethylene phosphate: an ab initio study by supermolecule model and polarizable continuum approach

The alkaline hydrolysis reaction of ethylene phosphate (EP) has been investigated using a supermolecule model, in which several explicit water molecules are included. The structures and single-point energies for all of the stationary points are calculated in the gas phase and in solution at the B3LYP/6-31++G(df,p) and MP2/6-311++G(df,2p) levels. The effect of water bulk solvent is introduced by the polarizable continuum model (PCM). Water attack and hydroxide attack pathways are taken into account for the alkaline hydrolysis of EP. An associative mechanism is observed for both of the two pathways with a kinetically insignificant intermediate. The water attack pathway involves a water molecule attacking and a proton transfer from the attacking water to the hydroxide in the first step, followed by an endocyclic bond cleavage to the leaving group. While in the first step of the hydroxide attack pathway the nucleophile is the hydroxide anion. The calculated barriers in aqueous solution for the water attack and hydroxide attack pathways are all about 22 kcal/mol. The excellent agreement between the calculated and observed values demonstrates that both of the two pathways are possible for the alkaline hydrolysis of EP.     

Neutral hydrolyses of carbon disulfide: An ab initio study of water catalysis

The water-catalyzed hydrolysis reaction of carbon disulfide (CS(2)) has been investigated at the levels of HF and MP2 with the basis set of 6-311++G(d,p) using the combined supramolecular/continuum models, in which up to six water molecules are involved in the hydrolysis and the effect of water bulk solvent is taken into account according to the polarizable continuum model (PCM). The activation Gibbs free energies in water solution, DeltaG(sol) (not equal) (298 K), for the rate-determining steps of one up to six water hydrolyses are 247.9, 184.2, 152.3, 141.8, 134.4, and 118.9 kJ/mol, respectively. The most favorable hydrolysis path of CS(2) involves a sort of eight-membered ring transition structure formed by six water molecules, among which three water molecules are not involved in the proton transfer, two near to the nonreactive sulfur atom, and one below the parent carbon disulfide. This suggests that the hydrolysis of CS(2) can be mediated with the water molecule(s) and be significantly facilitated by the cooperative effects of the water molecule(s) in the nonreactive region. The catalytic effects of water molecule(s) due to the alleviation of ring strain in the proton transfer process may result from the synergistic effects of rehybridization and charge reorganization from the prereaction complex to the rate-determining transition state structure induced by water molecule(s). PCM solvation models could significantly lower the rate-determining activation Gibbs free energies by 20-38 kJ/mol when two up to six explicit water molecules involved in the neutral hydrolysis of CS(2).     

Solvent effects and activation parameters in the competitive cleavage of C-CN and C-H bonds in 2-methyl-3-butenenitrile using [(dippe)NiH]2

The reaction of [(dippe)NiH]2 with 2-methyl-3-butenenitrile (2M3BN) in solvents spanning a wide range of polarities shows significant differences in the ratio of C-H and C-CN activated products. C-H cleavage is favored in polar solvents, whereas C-C cleavage is favored in nonpolar solvents. This variation is attributed to the differential solvation of the transition states, which was further supported through the use of sterically bulky solvents and weakly coordinating solvents. Variation of the temperature of reaction of [(dippe)NiH]2 with 2M3BN in decane and N,N-dimethylformamide (DMF) allowed for the calculation of Eyring activation parameters for the C-CN activation and C-H activation mechanisms. The activation parameters for the C-H activation pathway were DeltaH(double dagger) = 11.4 +/- 5.3 kcal/mol and DeltaS(double dagger) = -45 +/- 15 e.u., compared with DeltaH(double dagger) = 17.3 +/- 2.6 kcal/mol and DeltaS(double dagger) = -29 +/- 7 e.u. for the C-CN activation pathway. These parameters indicate that C-H activation is favored enthalpically, but not entropically, over C-C activation, implying a more ordered transition state for the former.     

Halogenation of cubane under phase-transfer conditions: single and double C-H-bond substitution with conservation of the cage structure.

The first highly selective C-H chlorination, bromination, and iodination of cubane (1) utilizing polyhalomethanes as halogen sources under phase-transfer (PT) conditions is described. Isomeric dihalocubanes with all possible combinations of chlorine, bromine, and iodine in ortho, meta, and para positions were also prepared by this method; m-dihalo products form preferentially. Ab initio and density functional theory (DFT) computations were used to rationalize the pronounced differences in the reactions of 1 with halogen (Hal(*)) vs carbon-centered trihalomethyl (Hal(3)C(*)) radicals (Hal = Cl, Br). For Hal(3)C radicals the C-H abstraction pathway is less unfavorable (DeltaG(double dagger)(298) = 21.6 kcal/mol for Cl(3)C(*) and 19.4 kcal/mol for Br(3)C(*) at B3LYP/6-311+G//B3LYP/6-31G) than the fragmentation of the cubane skeleton via S(H)2-attack on one of the carbon atoms of 1 (DeltaG(double dagger)(298) = 33.8 and 35.1 kcal/mol, respectively). In stark contrast, the reaction of 1 with halogen atoms preferentially follows the fragmentation pathway (DeltaG(double dagger)(298) = 2.1 and 7.5 kcal/mol) and C-H abstraction is more unfavorable (DeltaG(double dagger)(298) = 4.6 and 12.0 kcal/mol). Our computational results nicely agree with the behavior of 1 under PT halogenation conditions (where Hal(3)C(*) is involved in the activation step) and under free-radical photohalogenation with Hal(2) (Della, E. W., et al. J. Am. Chem. Soc. 1992, 114, 10730). The incorporation of a second halogen atom preferentially in the meta position of halocubanes demonstrates the control of the regioselectivity by molecular orbital symmetry.     

Transition states for the dimerization of 1,3-cyclohexadiene: a DFT, CASPT2, and CBS-QB3 quantum mechanical investigation

Quantum mechanical calculations using restricted and unrestricted B3LYP density functional theory, CASPT2, and CBS-QB3 methods for the dimerization of 1,3-cyclohexadiene (1) reveal several highly competitive concerted and stepwise reaction pathways leading to [4 + 2] and [2 + 2] cycloadducts, as well as a novel [6 + 4] ene product. The transition state for endo-[4 + 2] cycloaddition (endo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 28.7 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 19.0 kcal/mol) is not bis-pericyclic, leading to nondegenerate primary and secondary orbital interactions. However, the C(s) symmetric second-order saddle point on the B3LYP energy surface is only 0.3 kcal/mol above endo-2TS. The activation enthalpy for the concerted exo-[4 + 2] cycloaddition (exo-2TS, DeltaH(double dagger)(B3LYP(0K)) = 30.1 kcal/mol and DeltaH(double dagger)(CBS-QB3(0K)) = 21.1 kcal/mol) is 1.4 kcal/mol higher than that of the endo transition state. Stepwise pathways involving diallyl radicals are formed via two different C-C forming transition states (rac-5TS and meso-5TS) and are predicted to be competitive with the concerted cycloaddition. Transition states were located for cyclization from intermediate rac-5 leading to the endo-[4 + 2] (endo-2) and exo-[2 + 2] (anti-3) cycloadducts. Only the endo-[2 + 2] (syn-3) transition state was located for cyclization of intermediate meso-5. The novel [6 + 4] "concerted" ene transition state (threo-4TS, DeltaH(double dagger)(UB3LYP(0K)) = 28.3 kcal/mol) is found to be unstable with respect to an unrestricted calculation. This diradicaloid transition state closely resembles the cyclohexadiallyl radical rather than the linked cyclohexadienyl radical. Several [3,3] sigmatropic rearrangement transition states were also located and have activation enthalpies between 27 and 31 kcal/mol.     

Disproportionation of Gold(II) complexes. A density functional study of ligand and solvent effects

A computational study of gold(II) disproportionation is presented for the atomic ion as well as complexes with chloride and neutral ligands. The Au2+ atomic ion is stable to disproportionation, but the barrier is more than halved to 119 kcal/mol in an aqueous environment vs 283 kcal/mol in the gas phase. For dissociative disproportionation of chloride complexes, the loss of chlorine, either as an atom (Delta G(aq) = +20 kcal/mol) or as an anion (Delta G(aq) = +15 kcal/mol) represents the largest calculated barrier. The calculated transition state for associative disproportionation is only 9 kcal/mol above separated Au(II)Cl3(-) anions. For the disproportionation of Au(II)L3 complexes with neutral ligands, disproportionation is highly endergonic in the gas phase. Calculations imply that for synthesis of a monometallic Au(II) complex, a nonpolar solvent is preferred. With the exception of [Au(CO)3]2+, disproportionation of Au(II)L3 complexes to Au(I)L and Au(III)L3 is exergonic in solution phase for the ligands investigated. The driving force is provided by the very favorable solvation free energy of the trivalent gold complex. The solvation free energy contribution to the reaction (Delta G(solv)) is very large for small and polar ligands such as ammonia and water. Furthermore, calculations imply that choosing ligands that would yield neutral species upon disproportionation may provide an effective route to thwart this decomposition pathway for Au(II) complexes. Likewise, bulkier ligands that yield larger, more weakly solvated complex ions would appear to be desirable.     

Experimental and theoretical aspects of the haptotropic rearrangement of diiron and diruthenium carbonyl complexes bound to 4,6,8-trimethylazulene

The haptotropic rearrangement of dinuclear metal carbonyl species on the conjugate pi-ligand of (micro2,eta3:eta5-4,6,8-trimethylazulene)M2(CO)5 [M = Fe (3) and Ru (4)] was investigated in detail both experimentally and theoretically. The complexes, 3 and 4, were synthesized and characterized by spectroscopy and crystallography. The spin saturation transfer technique of 1H NMR was used to measure the rate constant k of the haptotropic isomerization between the two enantiomers of 3 and 4, from which thermodynamic parameters were determined: (3; deltaS(double dagger) = -7 +/- 1 cal K(-1) mol(-1), deltaH(double dagger) = 22 +/- 1 cal mol(-1), deltaG(double dagger)373 = 25 +/- 1 cal mol(-1)), (4; deltaS(double dagger) = 7 +/- 1 cal K(-1) mol(-1), deltaH(double dagger) = 25 +/- 1 cal mol(-1), deltaG(double dagger)373 = 23 +/- 1 cal mol(-1)). DFT calculations (the B3LYP, B1B95 and PBE1PBE methods) were also carried out using the CEP-31G and cc-pVDZ as the basis set of the transition metal and other elements, respectively, by which both ground state and transition state structures were optimized for the haptotropic rearrangement of 3 and 4. The potential energy surface for these reactions suggests that the reaction involves the conversion of the coordination mode from micro2eta3,eta5- (ground state) to micro2,eta1,eta5- (transition state). Mechanistic consideration, in particular that of differences in transition states between the diiron and diruthenium complexes, is also described.     

P-S Bond scission by bis(cyclopentadienyl)molybdenum(IV) dichloride, Cp(2)MoCl(2)(aq): first documented example of an organometallic complex hydrolyzing thiophosphinates

Thiophosphinate hydrolysis involving P-S bond scission is desirable for the degradation of organophosphate neurotoxins, and we report the first case for such a hydrolytic process by an organometallic compound. The metallocene, bis(cyclopentadienyl)molybdenum(IV) dichloride, Cp(2)MoCl(2) (Cp = eta(5)-C(5)H(5)), hydrolyzes a variety of thioaryl diphenylphosphinates in an aqueous THF solution. P-S scission of p-methoxythiophenyl diphenylphosphinate has a 500-fold rate of acceleration in the presence of Cp(2)MoCl(2)(aq) with activation parameters of 20(3) kcal mol(-)(1) and -15(3) cal mol(-)(1) K(-)(1) for DeltaH(double dagger) and DeltaS(double dagger), respectively. These activation parameters and the rate acceleration are consistent with an intermolecular hydrolytic process in which the Cp(2)Mo serves as a Lewis acid to activate the phosphinate for nucleophilic attack. Furthermore, rho = 2.3 (25 degrees C) which indicates a single nonconcerted mechanism in which the rate determining step is the nucleophilic attack on the activated phosphinate.     

Proton solvation in protic and aprotic solvents

Protonation pattern strongly affects the properties of molecular systems. To determine protonation equilibria, proton solvation free energy, which is a central quantity in solution chemistry, needs to be known. In this study, proton affinities (PAs), electrostatic energies of solvation, and pK_A values were computed in protic and aprotic solvents. The proton solvation energy in acetonitrile (MeCN), methanol (MeOH), water, and dimethyl sulfoxide (DMSO) was determined from computed and measured pK_A values for a specially selected set of organic compounds. pK_A values were computed with high accuracy using a combination of quantum chemical and electrostatic approaches. Quantum chemical density functional theory computations were performed evaluating PA in the gas‐phase. The electrostatic contributions of solvation were computed solving the Poisson equation. The computations yield proton solvation free energies with high accuracy, which are in MeCN, MeOH, water, and DMSO ?255.1, ?265.9, ?266.3, and ?266.4 kcal/mol, respectively, where the value for water is close to the consensus value of ?265.9 kcal/mol. The pK_A values of MeCN, MeOH, and DMSO in water correlates well with the corresponding proton solvation energies in these liquids, indicating that the solvated proton was attached to a single solvent molecule. © 2016 Wiley Periodicals, Inc.     

Predicting the energy of the water exchange reaction and free energy of solvation for the uranyl ion in aqueous solution

The structures and vibrational frequencies of UO2(H2O)4(2+) and UO2(H2O)5(2+) have been calculated using density functional theory and are in reasonable agreement with experiment. The energies of various reactions were calculated at the density functional theory (DFT) and MP2 levels; the latter provides the best results. Self-consistent reaction field calculations in the PCM and SCIPCM approximations predicted the free energy of the water exchange reaction, UO2(H2O)4(2+) + H2O <--> UO2(H2O)5(2+). The calculated free energies of reaction are very sensitive to the choice of radii (O and H) and isodensity values in the PCM and SCIPCM models, respectively. Results consistent with the experimental HEXS value of -1.19 +/- 0.42 kcal/mol (within 1-3 kcal/mol) are obtained with small cavities. The structures and vibrational frequencies of the clusters with second solvation shell waters: UO2(H2O)4(H2O)8(2+), UO2(H2O)4(H2O)10(2+), UO2(H2O)4(H2O)11(2+), UO2(H2O)5(H2O)7(2+), and UO2(H2O)5(H2O)10(2+), were calculated and are in better agreement with experiment as compared to reactions involving only UO2(H2O)4(2+) and UO2(H2O)5(2+). The MP2 reaction energies for water exchange gave gas-phase results that agreed with experiment in the range -5.5 to +3.3 kcal/mol. The results were improved by inclusion of a standard PCM model with differences of -1.2 to +2.7 kcal/mol. Rearrangement reactions based on an intramolecular isomerization leading to a redistribution of water in the two shells provide good values in comparison to experiment with values of Delta G(exchange) from -2.2 to -0.5 kcal/mol so the inclusion of a second hydration sphere accounts for most solvation effects. Calculation of the free energy of solvation of the uranyl cation yielded an upper bound to the solvation energy of -410 +/- 5 kcal/mol, consistent with the best experimental value of -421 +/- 15 kcal/mol.     

Conformational equilibrium of 1,2-dichloroethane in water: comparison of PCM and RISM-SCF methods

The RISM-SCF and polarizable continuum model (PCM) approaches have been applied to study the conformational equilibrium of 1,2-dichloroethane (DCE) in water. Both the electron correlation effect and basis sets play an important role in the relative energies of the gauche and trans conformers in gas and solution phases. Both PCM and RISM-MP2 methods resulted in a consistent trend with the previous experimental and theoretical studies that the population of the gauche conformer increases in going from the gas phase to the aqueous solution. However, the PCM treatment could not describe the solvent effect completely in that the sign of the relative free energy of the gauche and trans forms is opposite to the most recent experimental and theoretical data, while the RISM-MP2 gives the right sign in the free energy difference. We found that the larger excess chemical potential gain (by ca. -4.1 kcal/mol) for the gauche conformer is large enough to result in the gauche preference of DCE in water, though it has to compensate for more solute reorganization energy (approximately 1.6 kcal/mol) and overcome the energy difference (approximately 1.6 kcal/mol) in the gas phase. The radial distribution functions between DCE and the nearest water shows that the electrostatic repulsion between chlorine and oxygen atoms is higher in the trans conformer than in the gauche one, while the attractive interaction between chlorine and hydrogen of water is higher in the gauche conformer.     

Synthesis,structure, and dynamic behavior of cyclopentadienyl-lithium, -sodium,and -potassium annelated with bicyclo[2.2.2]octene units: a systematic study on site exchange of alkali metals on a cyclopentadienyl ring in tetrahydrofuran

Novel cyclopentadienyl (Cp)-alkali metal complexes 1-M and 2-M (M = Li, Na, K), in which the Cp ring is annelated with two bicyclo[2.2.2]octene units and substituted with a phenyl group for 1 and a tert-butyl group for 2, were synthesized, and their structures and dynamic behaviors were investigated by means of X-ray crystallography, dynamic (13)C NMR, and DFT calculations. The X-ray crystallography results indicated that 1-Li, 1-Na, and 2-Na form monomeric contact ion pairs (CIP) with three THF molecules coordinated to the metal atom. Also, in THF-d(8), all of the 1-M and 2-M form monomeric CIP in the ground state. However, variable-temperature (13)C NMR measurements of 1-M and 2-M in THF-d(8) demonstrated dynamic behavior in which the metal ion exchanges positions between the upper and lower faces of the Cp ring. From a study of the concentration dependence of the dynamic behavior, the exchange was found to proceed principally as an intramolecular process at concentration ranges lower than 0.2 M. The experimentally observed deltaG values for the intramolecular exchange process for all the 1-M and 2-M (except for 2-Li, whose intramolecular process was too slow to observe) were found to be quite similar in THF-d(8) solution and to fall within the range of 12-14 kcal mol(-)(1). Within this range, a tendency was observed for the deltaG values to increase as the size of the metal decreased. Theoretical calculations (B3LYP/6-31G(d)) afforded considerably large values as the gas-phase dissociation energy for 1-M (162.7 kcal mol(-)(1) for M = Li; 131.6 kcal mol(-)(1) for M = Na; 110.9 kcal mol(-)(1) for M = K) and for 2-M (170.0 kcal mol(-)(1) for M = Li; 137.5 kcal mol(-)(1) for M = Na; 115.4 kcal mol(-)(1) for M = K). These values should be compensated for by a decrease in the solvation energies for the metal ions with increasing size, as exemplified by the calculated solvation energy for M(+)(Me(2)O)(4), which serves as a model for metal ions solvated with four molecules of THF (-122.9 kcal mol(-)(1) for M = Li; -94.7 kcal mol(-)(1) for M = Na; -67.7 kcal mol(-)(1) for M = K). This compensation results in a small difference in the overall energy for dissociation of 1-M or 2-M in ethereal solutions, thus supporting the similar deltaG values observed for the intramolecular metal exchange.