A theoretical investigation of the relative stability of hydrated glycine and methylcarbamic acid--from water clusters to interstellar ices

We have theoretically investigated how the low-energy conformers of the neutral and the zwitterionic forms of glycine as well as methylcarbamic acid are stabilized by the presence water. The MP2/6-311++G(d,p) method was utilized to conduct calculations on glycine and methylcarbamic acid in both isolated clusters and in clusters embedded in the conductor-like polarizable continuum model (C-PCM), where the clusters explicitly contain between one and ten water molecules. The neutral forms of glycine and methylcarbamic acid were found to have similar hydration energies, whereas the neutral methylcarbamic acid was determined to be approximately 32 kJ mol(-1) more stable than the neutral glycine in the isolated clusters and 30 kJ mol(-1) more stable in the C-PCM embedded clusters. Both the number and strength of the hydrogen bonding interactions between water and the zwitterions drive the stability. This lowers the relative energy of the glycine zwitterion from 50 kJ mol(-1) above neutral glycine, when there are two water molecules in the clusters to 11 kJ mol(-1) below for the clusters containing ten water molecules. For the methylcarbamic acid clusters with two water molecules, the zwitterion is 51 kJ mol(-1) higher in energy than the neutral form, but it remains 13 kJ mol(-1) above the neutral methylcarbamic acid in the clusters containing ten water molecules. When the bulk water environment is simulated by the C-PCM calculations, we find both the methylcarbamic acid and glycine zwitterionic forms have similar energies at 20 kJ mol(-1) above the neutral methylcarbamic acid energy and 10 kJ mol(-1) lower than the neutral glycine energy. Although neither methylcarbamic acid nor glycine have been detected in the interstellar medium yet, our findings indicate that methylcarbamic acid is the more stable product from methylamine and carbon dioxide reactions in a water ice. This suggests that methylcarbamic acid likely plays a role in the intermediate steps if glycine is formed in the interstellar medium.     

Strategic use of amino acid N-substituents to limit α-carbon-centered radical formation and consequent loss of stereochemical integrity

Ab initio calculations have been used to investigate the effect of N-substituents on the stability of α-carbon-centered amino acid radicals. Optimized structures of glycine derivatives and related substituted methanes, and the corresponding radicals, were determined with B3-LYP/6-31G(d). Single-point RMP2/6-31G(d) calculations on these structures were then used to obtain radical stabilization energies, which were compared with the relative rates of formation of the same or closely similar radicals in reactions with N-bromosuccinimde. These studies show that N-acylation and sulfonation decrease both the stability and the ease of formation of the α-carbon-centered radicals. Greater effects are seen with fluoroacyl, fluoroalkylsulfonyl and imido groups. The extent of the effect of the imido and fluoroalkylsulfonyl groups is such that N-phthaloyl- and trifluoromethanesulfonyl-protected amino acids react by hydrogen-atom abstraction from the side chain, thereby avoiding reaction at the chiral α-center and preserving its stereochemical integrity. The origins of these substituent effects are examined.     

Effect of side chains on competing pathways for beta-scission reactions of peptide-backbone alkoxyl radicals

High-level quantum chemistry calculations have been carried out to investigate beta-scission reactions of alkoxyl radicals located at the alpha-carbon of a peptide backbone. This type of alkoxyl radical may undergo three possible beta-scission reactions, namely C-C beta-scission of the backbone, C-N beta-scission of the backbone, and C-R beta-scission of the side chain. We find that the rates for the C-C beta-scission reactions are all very fast, with rate constants of the order 10(12) s(-1) that are essentially independent of the side chain. The C-N beta-scission reactions are all slow, with rate constants that range from 10(-0.7) to 10(-4.5) s(-1). The rates of the C-R beta-scission reactions depend on the side chain and range from moderately fast (10(7) s(-1)) to very fast (10(12) s(-1)). The rates of the C-R beta-scission reactions correlate well with the relative stabilities of the resultant side-chain product radicals (*R), as reflected in calculated radical stabilization energies (RSEs). The order of stabilities for the side-chain fragment radicals for the natural amino acids is found to be Ala < Glu < Gln approximately Leu approximately Met approximately Lys approximately Arg < Asp approximately Ile approximately Asn approximately Val < Ser approximately Thr approximately Cys < Phe approximately Tyr approximately His approximately Trp. We predict that for side-chain C-R beta-scission reactions to effectively compete with the backbone C-C beta-scission reactions, the side-chain fragment radicals would generally need an RSE greater than approximately 30 kJ mol(-1). Thus, the residues that may lead to competitive side-chain beta-scission reactions are Ser, Thr, Cys, Phe, Tyr, His, and Trp.     

An ab initio and DFT study of some halogen atom transfer reactions from alkyl groups to acyl radical

Ab initio calculations using the 6-311G**, cc-pVDZ, and (valence) double-zeta pseudopotential (DZP) basis sets, with (MP2, QCISD, CCSD(T)) and without (HF) the inclusion of electron correlation, and density functional (BHandHLYP, B3LYP) calculations predict that the transition states for the reaction of acetyl radical with several alkyl halides adopt an almost collinear arrangement of attacking and leaving radicals at the halogen atom. Energy barriers (DeltaE(double dagger)) for these halogen transfer reactions of between 89.2 (chlorine transfer from methyl group) and 25.3 kJ mol(-1) (iodine transfer from tert-butyl group) are calculated at the BHandHLYP/DZP level of theory. While the difference in forward and reverse energy barriers for iodine transfer to acetyl radical is predicted to be 15.1 kJ mol(-1) for primary alkyl iodide, these values are calculated to be 6.7 and -4.2 kJ mol(-1) for secondary and tertiary alkyl iodide respectively. These data are in good agreement with available experimental data in that atom transfer radical carbonylation reactions are sluggish with primary alkyl iodides, but proceed smoothly with secondary and tertiary alkyl iodides. These calculations also predict that bromine transfer reactions involving acyl radical are also feasible at moderately high temperature.     

Bond dissociation energies and radical stabilization energies associated with model peptide-backbone radicals

Bond dissociation energies (BDEs) and radical stabilization energies (RSEs) have been calculated for a series of models that represent a glycine-containing peptide-backbone. High-level methods that have been used include W1, CBS-QB3, U-CBS-QB3, and G3X(MP2)-RAD. Simpler methods used include MP2, B3-LYP, BMK, and MPWB1K in association with the 6-311+G(3df,2p) basis set. We find that the high-level methods produce BDEs and RSEs that are in good agreement with one another. Of the simpler methods, RBMK and RMPWB1K achieve good accuracy for BDEs and RSEs for all the species that were examined. For monosubstituted carbon-centered radicals, we find that the stabilizing effect (as measured by RSEs) of carbonyl substituents (CX=O) ranges from 24.7 to 36.9 kJ mol(-1), with the largest stabilization occurring for the CH=O group. Amino groups (NHY) also stabilize a monosubstituted alpha-carbon radical, with the calculated RSEs ranging from 44.5 to 49.5 kJ mol(-1), the largest stabilization occurring for the NH2 group. In combination, NHY and CX=O substituents on a disubstituted carbon-centered radical produce a large stabilizing effect ranging from 82.0 to 125.8 kJ mol(-1). This translates to a captodative (synergistic) stabilization of 12.8 to 39.4 kJ mol(-1). For monosubstituted nitrogen-centered radicals, we find that the stabilizing effect of methyl and related (CH2Z) substituents ranges from 25.9 to 31.7 kJ mol(-1), the largest stabilization occurring for the CH3 group. Carbonyl substituents (CX=O) destabilize a nitrogen-centered radical relative to the corresponding closed-shell molecule, with the calculated RSEs ranging from -30.8 to -22.3 kJ mol(-1), the largest destabilization occurring for the CH=O group. In combination, CH2Z and CX=O substituents at a nitrogen radical center produce a destabilizing effect ranging from -19.0 to -0.2 kJ mol(-1). This translates to an additional destabilization associated with disubstitution of -18.6 to -7.8 kJ mol(-1).     

An Ab initio study of tautomerism between formhydroxamic acid and formhydroximic acid

Ab initio MO calculations including electron correlation with the 4-31G and 4-31G ** basis sets were performed in order to study the formhydroxamic acid-formhydroximic acid tautomerism. The geometries, relative energies, and activation energy of the tautomer and transition state were determined. Based on total-energy calculations at the MP 4/4-31G **//RHF /4-31G ** plus the scaled zero-point vibration energy level, the energy of formhydroxamic acid is determined to be lower than that of formhydroximic acid by 40.7 kJ/mol. The activation energy of the formhydroxamic acid-formhydroximic acid tautomerism via a 1,3-intramolecular hydrogen shift is 151.4 kJ/mol. © 1992 John Wiley & Sons, Inc.     

A DFT study of H-isomerisation in alkoxy-, alkylperoxy- and alkyl radicals: Some implications for radical chain reactions in polymer systems

Intramolecular 1-n H-shift (n = 2, 3… 7) reactions in alkoxy, alkyl and peroxy radicals were studied by density functional theory (DFT) at the B3LYP/6-311+G∗∗ level and compared with respective intermolecular H-transfers. It was found that starting from 1 to 3 H-shift the barrier heights stepwise decrease with increasing n reaching minimum for 1-5 and 1-6 H-shifts. This dependence can be ascribed to the decrease of the strain with increasing transition state (TS) ring size, which is minimal in six- and seven-member rings. The barrier heights of H-shifts in alkyl radicals are systematically larger than those in alkoxy radicals: the respective activation energies (E_a) of 1-5 and 1-6 H-shifts are about 59-67 kJ/mol for alkyl radical and 21-34 kJ/mol for alkoxy radicals. Further increase of the TS ring size in 1-7 H-shifts leads to the increase of the barrier to 44 kJ/mol in the hexyloxy radical and 84 kJ/mol for n-heptyl radical. We have also found that intermolecular H-transfer reactions in all three types of free radicals have smaller barriers than respective intramolecular 1-5 or 1-6 H-shifts by 4-25 kJ/mol. The mentioned difference can be explained in terms of enhanced nonbonding repulsion interaction in the cyclic TS structures compared to respective intermolecular TS. B3LYP/6-311+G∗∗ geometric parameters and imaginary frequencies for 1-n H-shifts TS are consistent with respective calculated barrier heights. Reactivity of some other radicals compared to alkoxy, peroxy and alkyl radicals as well as other factors influencing their reactivity (π-conjugation, steric effect and ring strain in cyclic TS, etc.) are also briefly discussed in relation to free radical reactions in polymer systems.     

Quantum mechanical investigation of the inner-sphere reorganization energy of cyclooctatetraene/cyclooctatetraene radical anion. Part I

The inner-sphere reorganization energy of the electron self-exchange of the couple cyclooctatetraene/cyclooctatetraene radical anion has been investigated by quantum mechanical calculations. The more stable Jahn Teller distorted B2g conformation of the radical anion has been used in this study. Two different theories have been applied in this first part. The harmonic approximation in the classical Marcus scheme has been modified by using projected force constants, which are obtained from the complete force constant matrix and the geometry changes of the molecule during the ET (introduced by Mikkelsen). A different approach (introduced by Nelsen) combines the different energies of the neutral and radical anion with and without relaxation corresponding to the vertical ionization potential and the vertical electron affinity. The electronic energies of the neutral molecule and the radical anion differ dramatically applying three different levels of quantum mechanical calculations (UAM1, UB3LYP, PMP2 with three different basis sets with and without diffuse functions). Nevertheless the Nelsen method gives almost consistent results for the inner-sphere reorganization energies: 120.1 kJ/mol for semiempirical UAM1 method, 159.3 kJ/mol, 156.4 kJ/mol and 158.3 kJ/mol for density functional UB3LYP/6-31G*, UB3LYP/6-31++G* and UB3LYP/AUG-cc-pVDZ calculations and 192.5 kJ/mol for ab-initio PMP2/6-31G* investigations, respectively. These values are in agreement with earlier experimental work supposing the total reorganization energy to be larger than 38 kcal/mol assuming an electron self-exchange rate of 10(4) M(-1) s(-1). The simple harmonic approximation of Marcus relation has not yet been applied for a molecule like cyclooctatetraene with large torsional geometry changes. Using the projected force constants after scaling, considerably different results for the inner-sphere reorganization energy have been calculated: 738.1 kJ/mol for the UB3LYP/6-31G*, 743.3 kJ/mol for UB3LYP/6-31++G* and 759.1 kJ/mol for UB3LYP/AUG-cc-pVDZ level of theory. Comparison with our concentration dependent EPR experiments are controversial to the earlier experimental results, but the latter supports the assumption that the electron self-exchange occurs in a time scale so that the molecules cannot complete their vibrational motions. Therefore the projected Marcus relation is not valid for cyclooctatetraene/cyclooctatetraene radical anion including a large torsional change during the electron transfer.     

Intramolecular homolytic substitution of sulfinates and sulfinamides--a computational study

Ab initio and density functional theory (DFT) calculations predict that intramolecular homolytic substitution by alkyl radicals at the sulfur atom in sulfinates proceeds through a smooth transition state in which the attacking and leaving radicals adopt a near collinear arrangement. When forming a five-membered ring and the leaving radical is methyl, G3(MP2)-RAD//ROBHandHLYP/6-311++G(d,p) calculations predict that this reaction proceeds with an activation energy (ΔE(1)(?)) of 43.2 kJ mol(-1). ROBHandHLYP/6-311++G(d,p) calculations suggest that the formation of five-membered rings through intramolecular homolytic substitution by aryl radicals at the sulfur atom in sulfinates and sulfinamides, with expulsion of phenyl radicals, proceeds with the involvement of hypervalent intermediates. These intermediates further dissociate to the observed products, with overall energy barriers of 45-68 kJ mol(-1), depending on the system of interest. In each case, homolytic addition to the phenyl group competes with substitution, with calculated barriers of 51-78 kJ mol(-1). This computational study complements and provides insight into previous experimental observations.     

S-to-αC Radical Migration in the Radical Cations of Gly-Cys and Cys-Gly

The radical cations of Cys-Gly and Gly-Cys were studied using ion-molecule reactions (IMR), infrared multiple-photon dissociation (IRMPD) spectroscopy, and density functional theory (DFT) calculations. Homolytic cleavage of the S–NO bond of nitrosylated precursors generated radical cations with the radical site initially located on the sulfur atom. Time-resolved ion-molecule reactions showed that radical site migration via hydrogen atom transfer (HAT) occurred much more quickly in Gly-Cys^•+ than in Cys-Gly^•+. IRMPD and DFT calculations indicated that for Gly-Cys, the radical migrated from the sulfur atom to the α-carbon of glycine, which is lower in energy than the sulfur radical (–53.5 kJ/mol). This migration does not occur for Cys-Gly because the glycine α-carbon is higher in energy than the sulfur radical (10.3 kJ/mol). DFT calculations showed that the highest energy barriers for rearrangement are 68.2 kJ/mol for Gly-Cys and 133.8 kJ/mol for Cys-Gly, which is in agreement with both the IMR and IRMPD data and explains the HAT in Gly-Cys.     

Theoretical study on the isomerization behavior between alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals

[reaction: see text] Ab initio calculations using 6-311G**, cc-pVDZ, aug-cc-pVDZ, and a (valence) double-zeta pseudopotential (DZP) basis set, with (QCISD, CCSD(T)) and without (UHF) the inclusion of electron correlation, and density functional methods (BHandHLYP, B3LYP) predict that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals exist as isomers. At the CCSD(T)/cc-pVDZ//BHandHLY/cc-pVDZ level of theory, energy barriers of 15.1 and 17.7-21.7 kJ mol(-)(1) are calculated for the isomerization of s-trans-propenoyl and s-trans-crotonoyl radical to ketenylmethyl and 1-ketenylethyl radical, respectively. Similar results are obtained for the reactions of s-trans isomers involving silyl, germyl, and stannyl groups with energy barriers (DeltaE++) of 12.2-12.4, 13.1-13.9, and 12.9-18.2 kJ mol(-)(1) at the CCSD(T)/DZP//BHandHLYP/DZP calculation, respectively. These results suggest that alpha,beta-unsaturated acyl radicals and alpha-ketenyl radicals are not canonical forms but are isomeric species that can rapidly interconvert.     

Experimental and theoretical studies of sodium cation interactions with the acidic amino acids and their amide derivatives

The binding of Na+ to aspartic acid (Asp), glutamic acid (Glu), asparagine (Asn), and glutamine (Gln) is examined in detail by studying the collision-induced dissociation (CID) of the four sodiated amino acid complexes with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Analysis of the energy-dependent CID cross sections provides 0 K sodium cation affinities for the complexes after accounting for unimolecular decay rates, internal energy of the reactant ions, and multiple ion-molecule collisions. Quantum chemical calculations for a number of geometric conformations of each Na+(L) complex are determined at the B3LYP/6-311+G(d,p) level with single-point energies calculated at MP2(full), B3LYP, and B3P86 levels using a 6-311+G(2d,2p) basis set. This coordinated examination of both experimental work and quantum chemical calculations allows the energetic contributions of individual functionalities as well as steric influences of relative chain lengths to be thoroughly explored. Na+ binding affinities for the amide complexes are systematically stronger than those for the acid complexes by 14 +/- 1 kJ/mol, which is attributed to an inductive effect of the OH group in the carboxylic acid side chain. Additionally, the Na+ binding affinity for the longer-chain amino acids (Glx) is enhanced by 4 +/- 1 kJ/mol compared to the shorter-chain Asx because steric effects are reduced.     

Multi-component orbital interactions during oxyacyl radical addition reactions involving imines and electron-rich olefins

Ab initio and DFT calculations reveal that oxyacyl radicals add to imines and electron-rich olefins through simultaneous SOMO-pi*, SOMO-pi and pi*-HOMO interactions between the radical and the radicalophile. At the BHandHLYP/aug-cc-pVDZ level, energy barriers of 20.3 and 22.0 kJ mol(-1) are calculated for the attack of methoxycarbonyl radical at the carbon and nitrogen ends of methanimine, respectively. In comparison, barriers of 22.0 and 8.6 kJ mol(-1) are calculated at BHandHLYP/aug-cc-pVDZ for reaction of methoxycarbonyl radical at the 1- and 2-positions in aminoethylene, respectively. Natural bond orbital (NBO) analysis at the BHandHLYP/6-311G** level of theory reveals that SOMO-pi*, SOMO-pi and pi*-LP interactions are worth 111, 394 and 55 kJ mol(-1) respectively in the transition state (8) for reaction of oxyacyl radical at the nitrogen end of methanimine; similar interactions are observed for the chemistry involving aminoethylene. These multi-component interactions are responsible for the unusual motion vectors associated with the transition states involved in these reactions.     

An ab initio and DFT study of radical addition reactions of imidoyl and thioyl radicals to methanimine

Ab initio and DFT calculations reveal that both imidoyl and thioyl radicals add to the nitrogen end of methanimine through simultaneous SOMO-π*(imine), SOMO-π(imine), SOMO-LP(N) and π*(radical)-LP(N) interactions between the radical and the imine. At the CCSD(T)/cc-pVDZ//BHandHLYP/cc-pVTZ level of theory, barriers of 13.8 and 26.1 kJ mol(-1) are calculated for the attack of the methylimidoyl radical at the carbon- and nitrogen- end of methanimine, respectively, indicating that the imidoyl radial has a preference for addition to the nitrogen end of imine. On the other hand, barriers of 25.1 and 13.4 kJ mol(-1) are calculated at the same level of theory for the addition reaction of the methanethioyl radical at the carbon- and nitrogen- end of methanimine, respectively. Natural bond orbital (NBO) analysis at the BHandHLYP/6-311G** level of theory reveals that SOMO-π*(imine), SOMO-π(imine), SOMO-LP(N) and π*(radical)-LP(N) interactions are worth 111, 89, 115 and 17 kJ mol(-1), respectively, in the transition state (4) for the reaction of methylimidoyl radical at the nitrogen end of methanimine; similar interactions are observed for the chemistry involving all the radicals studied here. These multi-component interactions are responsible for the unusual motion vectors associated with the transition states involved in these reactions.     

Epoxide formation by ring closure of the cinnamyloxy radical

[formula: see text] The cinnamyloxy and oxiranyl benzyl radicals were generated by photolysis of alkyl 4-nitrobenzenesulfenates. The yet unprecedented epoxide ring formation from a primary alkoxy radical was observed. Experimental evidence supports the fact that the mode of ring opening of the oxiranyl carbinyl radical system is thermodynamically driven. B3LYP/6-31G* calculations indicate that the closed form of the radical is approximately 5 kcal/mol more stable than the open one.     

Thermochemistry, bond energies, and internal rotor potentials of dimethyl tetraoxide

Thermochemical properties of dimethyl tetraoxide (CH(3)OOOOCH(3)), the dimer of the methylperoxy radical, are studied using ab initio and density functional theory methods. Methylperoxy radicals are known to be important intermediates in the tropospheric ozone cycle, and the self-reaction of methylperoxy radicals, which is thought to proceed via dimethyl tetraoxide, leads to significant chain radical termination in this process. Dimethyl tetraoxide has five internal rotors, three of them unique; the potential energy profiles are calculated for these rotors, as well as for those in the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals. The dimethyl tetraoxide internal rotor profiles show barriers to rotation of 2-8 kcal mol(-1). Using B3LYP/6-31(d) geometries, frequencies, internal rotor potentials, and moments of inertia, we determine entropy and heat capacity values for dimethyl tetraoxide and its radicals. Isodesmic work reactions with the G3B3 and CBS-APNO methods are used; we calculate this enthalpy as -9.8 kcal mol(-1). Bond dissociation energies (BDEs) are calculated for all C-O and O-O bonds in dimethyl tetraoxide, again with the G3B3 and CBS-APNO theoretical methods, and we suggest the following BDEs: 46.0 kcal mol(-1) for CH(3)-OOOOCH(3), 20.0 kcal mol(-1) for CH(3)O-OOOCH(3), and 13.9 kcal mol(-1) for CH(3)OO-OOCH(3). From the BDE calculations and the isodesmic enthalpy of formation for dimethyl tetraoxide, we suggest enthalpies of 2.1, 5.8, and 1.4 kcal mol(-1) for the CH(3)OO, CH(3)OOO, and CH(3)OOOO radicals, respectively. We evaluate the suitability of 10 different density functional theory (DFT) methods for calculating thermochemical properties of dimethyl tetraoxide and its radicals with the 6-31G(d) and 6-311++G(3df,3pd) basis sets, using a variety of work reaction schemes. Overall, the best-performed DFT methods of those tested were TPSSh, BMK, and B1B95. Significant improvements in accuracy were made by moving from atomization to isodesmic work reactions, with most DFT methods yielding errors of less than 2 kcal mol(-1) with the 6-311++G(3df,3pd) basis set for isodesmic calculations on the dimethyl tetraoxide enthalpy. These isodesmic calculations were basis set consistent, with a considerable reduction in error found by using the 6-311++G(3df,3pd) basis set over the 6-31G(d) basis set. This was not the case, however, for atomization and bond dissociation work reactions, where the two basis sets returned similar results. Improved group additivity terms for the O-O-O moiety (O/O2 central atom group) are also determined.     

Unusual structural and energetic features of homolytic bond dissociation: from tetrakis(disyl)diphosphine to tetrakis(di-tert-butylsilyl)hydrazine

Quantum chemical calculations of the structures and thermodynamics of homolytic dissociation of the central P-P and N-N bonds in tetrakis(disyl)diphosphine and tetrakis(di-tert-butylsilyl)hydrazine have been performed. The theory predicted negative standard enthalpies for homolytic bond dissociation in both cases, -71.0 and -108.4 kJ mol(-1) for the diphosphine and hydrazine, respectively, using the ONIOM (MP2/6-31+G*:B3LYP/3-21G*) level. The dissociation is accompanied by considerable structural changes in the radicals as compared to the corresponding fragments of the parent molecules, resulting in low dissociation enthalpies. The most pronounced changes in both radicals are the relaxation of bond angles in the substituents and a conformational change in the orientation of the substituent groups. In addition, the bis(di-tert-butylsilyl)aminyl radical displays a considerable increase in Si-N-Si angle and shortening of the Si-N bonds upon dissociation. These changes are not associated with any appreciable delocalisation of the lone electron, as the spin density is found from the B3LYP/3-21G* calculations to be largely concentrated on the nitrogen atom. It has been also shown that although the dissociation energies are low for both compounds, the intrinsic energies of the central bonds are still high, 140.6 kJ mol(-1) for the P-P bond in tetrakis(disyl)diphosphine and 490.6 kJ mol(-1) for the N-N bond in tetrakis(di-tert-butylsilyl)hydrazine, using the ONIOM method. The calculations predict that the dissociation of tetrakis(disyl)diphosphine would have negative free energy even without taking relaxation of the fragments into account, while the full potential of releasing about 306 kJ mol(-1) of energy stored in the ligands of tetrakis(di-tert-butylsilyl)hydrazine is only fully realised upon a considerable separation of the fragments.     

Modelling the hydrolysis of succinimide: formation of aspartate and reversible isomerization of aspartic acid via succinimide

In the present study, we have modelled the nucleophilic attack of water and a hydroxyl anion on the carbonyl carbon of a succinimide derivative leading to aspartate and aspartic acid. Calculations have been carried out at the B3LYP/6-31 +G* level in a vacuum. The IEF-PCM methodology has been used to carry out single point calculations in solution. In neutral medium, hydrolysis is facilitated by the presence of a polar continuum, whereas in basic medium the polar environment hinders the hydrolysis of succinimide. The deltaH degrees and deltaS degrees values for the cyclization reactions of aspartic acid yielding succinimide are 29.2 kJ mol(-1) and 133.5 kJ mol(-1) K(-1) respectively in accordance with the experimental results on the isomerization of the Ac-Asp-Gly-NHMe dipeptide unit. In a neutral medium, the isoaspartate: aspartate is found to be 2.2:1 in a vacuum and 3.4:1 in solution, in line with the experimental findings based on the hydrolysis of a tetrapeptide (Ac-Gly-Asn-Gly-Gly-NHMe) and a hexapeptide (Val-Tyr-Pro-Asn-Gly-Ala) where this ratio was found to be 3.1:1.     

Molecular orbital calculations of ring opening of the isoelectronic cyclopropylcarbinyl radical, cyclopropoxy radical, and cyclopropylaminium radical cation series of radical clocks

Detailed molecular orbital calculations were directed to the cyclopropylcarbinyl radical (1), the cyclopropoxy radical (2), and the cyclopropylaminium radical cation (3) as well as their ring-opened products. Since a considerable amount of data are published about cyclopropylcarbinyl radicals, calculations were made for this species and related ring-opened products as a reference for 2 and 3 and their reactions. Radicals 1-3 have practical utility as "radical clocks" that can be used to time other radical reactions. Radical 3 is of further interest in photoelectron-transfer processes where the back-electron-transfer process may be suppressed by rapid ring opening. Calculations have been carried out at the UHF/6-31G*, MP4//MP2/6-31G*, DFT B3LYP/6-31G*, and CCSD(T)/cc-pVTZ//QCISD/cc-pVDZ levels. Energies are corrected to 298 K, and the barriers between species are reported in terms of Arrhenius E(a) and log A values along with differences in enthalpies, free energies, and entropies. The CCSD(T)-calculated energy barrier for ring opening of 1 is E(a) = 9.70, DeltaG* = 8.49 kcal/mol, which compares favorably to the previously calculated value of E(a) = 9.53 kcal/mol by the G2 method, but is higher than an experimental value of 7.05 kcal/mol. Our CCSD(T)-calculated E(a) value is also higher by 1.8 kcal/mol than a previously reported CBS-RAD//B3LYP/6-31G* calculation. The cyclopropoxy radical has a very small barrier to ring opening (CCSD(T), E(a) = 0.64 kcal/mol) and should be a very sensitive time clock. Of the three series studied, the cyclopropylaminium radical cation is most complex. In agreement with experimental data, bisected cyclopropylaminium radical cation is not found, but instead a ring-opened species is found. A perpendicular cyclopropylaminium radical cation (4) was found as a transition-state structure. Rotation of the 2p orbital in 4 to the bisected array results in ring opening. The minimum onset energy of photoionization of cyclopropylamine was calculated to be 201.5 kcal/mol (CCSD(T)) compared to experimental values of between about 201 and 204 kcal/mol. Calculations were made on the closely related cyclopropylcarbinyl and bicyclobutonium cations. Stabilization of the bisected cyclopropylcarbinyl conformer relative to the perpendicular species is much greater for the cations (29.1 kcal/ mol, QCISD) compared to the radicals (3.10 kcal/mol, QCISD). A search was made for analogues to the bicyclobutonium cation in the radical series 1 and 2 and the radical cation series 3. No comparable species were found. A rationale was made for some conflicting calculations involving the cyclopropylcarbinyl and bicyclobutonium cations. The order of stability of the cyclopropyl-X radicals was calculated to be X = CH2 > X = O > X = NH2+, where the latter species has no barrier for ring opening. The relative rate of ring opening for cyclopropyl-X radicals X = CH2 to X = O was calculated to be 3.1 x 10(6) s(-1) at 298 K (QCISD).     

The azulene-to-naphthalene rearrangement revisited: a DFT study of intramolecular and radical-promoted mechanisms

Intramolecular and radical-promoted mechanisms for the rearrangement of azulene to naphthalene are assessed with the aid of density functional calculations. All intramolecular mechanisms have very high activation energies (>/=350 kJ mol(-1) from azulene) and so can only be competitive at temperatures above 1000 degrees C. Two radical-promoted mechanisms, the methylene walk and spiran pathways, dominate the reaction below this temperature. The activation energy for an orbital symmetry-allowed mechanism via a bicyclobutane intermediate is 382 kJ mol(-1). The norcaradiene-vinylidene mechanism that has been proposed in order to explain the formation of small amounts of 1-phenyl-1-buten-3-ynes from flash thermolysis of azulene has an activation energy of 360 kJ mol(-1); subtle features of the B3LYP/6-31G(d) energy surface for this mechanism are discussed. All intermediates and transition states on the spiran and methylene walk radical-promoted pathways have been located at the B3LYP/6-31G(d) level. Interconversion of all n-H-azulyl radicals via hydrogen shifts was also examined, and hydrogen shifts around the five-membered ring are competitive with the mechanisms leading to rearrangement to naphthalene, but those around the seven-membered ring are not. Conversion of a tricyclic radical to the 9-H-naphthyl radical is the rate-limiting transition state on the spiran pathway, and lies 164.0 kJ mol(-1) above that of the 1-H-azulyl radical. The transition state for the degenerate hydrogen shift between the 9-H-azulyl and 10-H-azulyl radicals is 7.4 kJ mol(-1) lower. Partial equilibration of the intermediates in the spiran pathway via this shift may therefore occur, and this can account for the surprising formation of 1-methylnaphthalene from 2-methylazulene. The rate-limiting transition state for the methylene walk pathway involves the concerted transfer of a methylene group from one ring to the other and lies 182.3 kJ mol(-1) above that of the 1-H-azulyl radical. It is shown that rearrangement via a combination of 31% methylene walk and 69% spiran pathways can account semiquantitatively for all the products from 1-(13)C-azulene, 9-(13)C-azulene, and 4,7-(13)C(2)-azulene, in addition to accounting for the products from methylazulenes, and the formation of naphthalene-d(0) and -d(2) from azulene-4-d. It is also pointed out that a small extension to the spiran pathway could provide an alternative explanation for the formation of 1-phenyl-1-buten-3-ynes.