Computational design of S‐nitrosothiol “click” reactions

To address a long‐standing problem of finding efficient reactions for chemical labeling of protein‐based S‐nitrosothiols (RSNOs), we computationally explored hitherto unknown (3+2) cycloaddition RSNO reactions with alkynes and alkenes. Nonactivated RSNO cycloaddition reactions have high activation enthalpy (>20 kcal/mol at the CBS‐QB3 level) and compete with alternative S—N bond insertion pathway. However, the (3+2) cycloaddition reaction barriers can be dramatically lowered by coordination of a Lewis acid to the N atom of the —SNO group. To exploit this effect, we propose to use reagents with Lewis acid and a strain‐activated carbon–carbon multiple bond linked by a rigid scaffold, which can react with RSNOs with small activation enthalpies (～5 kcal/mol) and high reaction exothermicities (～40 kcal/mol). The proposed efficient RSNO cycloaddition reactions can be used for future development of practical RSNO labeling reactions. © 2013 Wiley Periodicals, Inc.     

Theory of 1,3-dipolar cycloadditions: distortion/interaction and frontier molecular orbital models

Quantum chemical calculations of activation barriers and reaction energies for 1,3-dipolar cycloadditions by the high-accuracy CBS-QB3 method reveal previously unrecognized quantitative trends in activation barriers. The distortion/interaction model of reactivity explains why (1) there is a monotonic decrease of approximately 6 kcal/mol in the activation energy along the series oxides, imine, and ylide for the diazonium, nitrilium, and azomethine betaine classes of 1,3-dipoles; (2) nitrilium and azomethine betaines with the same trio of atoms have almost identical cycloaddition barrier heights; (3) barrier heights for the cycloadditions of a given 1,3-dipole with ethylene and acetylene have the same activation energies (mean absolute deviation of 0.6 kcal/mol) in spite of very different reaction thermodynamics (Delta DeltaH(rxn) range = 14-43 kcal/mol) and frontier molecular orbital (FMO) energy gaps. The energy to distort the 1,3-dipole and dipolarophile to the transition state geometry, rather than FMO interactions or reaction thermodynamics, controls reactivity for cycloadditions of 1,3-dipoles with alkenes or alkynes. A distortion/interaction energy analysis was also carried out on the transition states for the cycloadditions of diazonium dipoles with a set of substituted alkenes (CH2CHX, X = OMe, Me, CO 2Me, Cl, CN) and reveals that FMO interaction energies between the 1,3-dipole and the dipolarophile differentiate reactivity when transition state distortion energies are nearly constant.     

Activation Energies and Reaction Energetics for 1,3‐Dipolar Cycloadditions of Hydrazoic Acid with C?C and C?N Multiple Bonds from High‐Accuracy and Density Functional Quantum Mechanical Calculations

The reactions of hydrazoic acid (HN₃) with ethene, acetylene, formaldimine (H₂CNH), and HCN were explored with the high‐accuracy CBS‐QB3 method, as well as with the B3LYP and mPW1K density functionals. CBS‐QB3 predicts that the activation energies for the reactions of hydrazoic acid with ethylene, acetylene, formaldimine, and HCN have remarkably similar activation enthalpies of 19.0, 19.0, 21.6, and 25.2 kcal/mol, respectively. The reactions are calculated to have reaction enthalpies of −21.5 for triazoline formation from ethene, and −63.7 kcal/mol for formation of the aromatic triazole from acetylene. The reaction to form tetrazoline from formaldimine has a reaction enthalpy of −8 kcal/mol (ΔG_rxn=+5.6 kcal/mol), and the formation of tetrazole from HCN has a reaction enthalpy of −23.0 kcal/mol. The trends in the energetics of these processes are rationalized by differences in σ‐bond energies in the transition states and adducts, and the energy required to distort hydrazoic acid to its transition‐state geometry. The density functionals predict activation enthalpies that are in relatively good agreement with CBS‐QB3, the results differing from CBS‐QB3 results by ca. 1–2 kcal/mol. Significant errors are revealed for mPW1K in predicting the reaction enthalpies for all reactions.     

Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1,3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions

Activation barriers and reaction energetics for the three main classes of 1,3-dipolar cycloadditions, including nine different reactions, were evaluated with the MPW1K and B3LYP density functional methods, MP2, and the multicomponent CBS-QB3 method. The CBS-QB3 values were used as standards for 1,3-dipolar cycloaddition activation barriers and reaction energetics, and the density functional theory (DFT) and MP2 methods were benchmarked against these values. The MPW1K/6-31G* method and basis set performs best for activation barriers, with a mean absolute deviation (MAD) value of 1.1 kcal/mol. The B3LYP/6-31G* method and basis set performs best for reaction enthalpies, with a MAD value of 2.4 kcal/mol, while the MPW1K method shows large errors for reaction energetics. The MP2 method gives the expected systematic underestimation of barriers. Concerted and nearly synchronous transition structures are predicted by all DFT and MP2 methods. Also reported are revised estimated 0 K experimental activation enthalpies for a standard set of hydrocarbon pericyclic reactions and updated comparisons to experiment for DFT, ab initio, and multicomponent methods. B3LYP and MPW1K methods with MAD values of 1.5 and 2.1 kcal/mol, respectively, fortuitously outperform the multicomponent CBS-QB3 method, which has a MAD value of 2.3. The MAD value of the O3LYP functional improves to 2.4 kcal/mol from the previously reported 3.0 kcal/mol.     

Concerted vs stepwise mechanisms in dehydro-Diels-Alder reactions

The Diels-Alder reaction is not limited to 1,3-dienes. Many cycloadditions of enynes and a smaller number of examples with 1,3-diynes have been reported. These "dehydro"-Diels-Alder cycloadditions are one class of dehydropericyclic reactions which have long been used to generate strained cyclic allenes and other novel structures. CCSD(T)//M05-2X computational results are reported for the cycloadditions of vinylacetylene and butadiyne with ethylene and acetylene. Both concerted and stepwise diradical routes have been explored for each reaction, with location of relevant stationary points. Relative to 1,3-dienes, replacement of one double bond by a triple bond adds 6-6.5 kcal/mol to the activation barrier; a second triple bond adds 4.3-4.5 kcal/mol to the barrier. Product strain decreases the predicted exothermicity. In every case, a concerted reaction is favored energetically. The difference between concerted and stepwise reactions is 5.2-6.6 kcal/mol for enynes but diminishes to 0.5-2 kcal/mol for diynes. Experimental studies on intramolecular diyne + ene cycloadditions show two distinct reaction pathways, providing evidence for competing concerted and stepwise mechanisms. Diyne + yne cycloadditions connect with arynes and ethynyl-1,3-cyclobutadiene. This potential energy surface appears to be flat, with only a minute advantage for a concerted process; many diyne cycloadditions or aryne cycloreversions will proceed by a stepwise mechanism.     

Three characteristic reactions of alkynes with metal compounds in organic synthesis

Alkynes have two sets of mutually orthogonal π‐bonds that are different from the π‐bonds of alkenes. These π‐bonds are able to bond with transition metal compounds. Alkynes easily bond with the various kinds of compounds having a π‐bond such as carbon monoxide, alkenes, other alkynes and nitriles in the presence of the transition metal compounds. The most representative reaction of alkynes is called the Pauson–Khand reaction. The Pauson–Khand reactions include the cyclization of alkynes with alkenes and carbon monoxide in the presence of cobalt carbonyls. Similar Pauson–Khand reactions also proceed in the presence of other transition metal compounds. These reactions are the first type of characteristic reaction of alkynes. Other various kinds of cyclizations with alkynes also proceed in the presence of the transition metal compounds. These reactions are the second type of characteristic reaction of alkynes. These include cyclooligomerizations and cycloadditions. The cyclooligomerizations include mainly cyclotrimerizations and cyclotetramerizations, and the cycloadditions are [2 + 2], [2 + 2 + 1], [2 + 2 + 2], [3 + 2], [4 + 2], etc., type cycloadditions. Alkynes are fairly reactive because of the high s character of their σ‐bonds. Therefore, simple coupling reactions with alkynes also proceed besides the cyclizations. The coupling reactions are the third type of characteristic reactions of alkynes in the presence of, mainly, the transition metal compounds. These reactions include carbonylations, dioxycarbonylations, Sonogashira reactions, coupling reactions with aldehydes, ketones, alkynes, alkenes and allyl compounds. Copyright © 2008 John Wiley & Sons, Ltd.     

Intermolecular and intramolecular Diels-Alder cycloadditions of 3-ylidenepiperazine-2,5-diones and 5-acyloxy-2(1h)-pyrazinones

The 3-ylidenepiperazine-2,5-diones 16 and 39 and 5-acyloxy-2(1H)-pyrazinones 17 can serve as starting materials for the Diels-Alder reactions of alkenes and alkynes to the piperazine ring, under acidic conditions or in the presence of acetyl chloride, to afford tricyclic piperazine-2,5-diones 19, 20, 23-25, 27, 44, and 45. Intramolecular cycloadditions occur if 3-ylidenepiperazine-2,5-diones 30 and 32 are used as the starting materials. This procedure is a convenient path to bridged bicyclo[2.2.2]diazaoctane ring systems such as 31 and 33, the former being found in biologically active secondary mold metabolites, such as VM55599 (1) or brevianamide A (5), which have been isolated from various fungi. The synthesis of the indole compound 31 provided evidence for the proposed biochemical pathway with a Diels-Alder reaction as key step. Quantum chemical calculations have revealed that piperazinones with a cationic azadiene moiety are the most reactive species in Diels-Alder cycloadditions.     

Meeting the Challenge of Intermolecular Gold(I)‐Catalyzed Cycloadditions of Alkynes and Allenes

The development of gold(I)‐catalyzed intermolecular carbo‐ and hetero‐cycloadditions of alkynes and allenes has been more challenging than their intramolecular counterparts. Here we review, with a mechanistic perspective, the most fundamental intermolecular cycloadditions of alkynes and allenes with alkenes.     

Computational methods to calculate accurate activation and reaction energies of 1,3-dipolar cycloadditions of 24 1,3-dipoles

Theoretical calculations were performed on the 1,3-dipolar cycloaddition reactions of 24 1,3-dipoles with ethylene and acetylene. The 24 1,3-dipoles are of the formula X≡Y(+)-Z(-) (where X is HC or N, Y is N, and Z is CH(2), NH, or O) or X═Y(+)-Z(-) (where X and Z are CH(2), NH, or O and Y is NH, O, or S). The high-accuracy G3B3 method was employed as the reference. CBS-QB3, CCSD(T)//B3LYP, SCS-MP2//B3LYP, B3LYP, M06-2X, and B97-D methods were benchmarked to assess their accuracies and to determine an accurate method that is practical for large systems. Several basis sets were also evaluated. Compared to the G3B3 method, CBS-QB3 and CCSD(T)/maug-cc-pV(T+d)Z//B3LYP methods give similar results for both activation and reaction enthalpies (mean average deviation, MAD, < 1.5 kcal/mol). SCS-MP2//B3LYP and M06-2X give small errors for the activation enthalpies (MAD < 1.5 kcal/mol), while B3LYP has MAD = 2.3 kcal/mol. SCS-MP2//B3LYP and B3LYP give the reasonable reaction enthalpies (MAD < 5.0 kcal/mol). The B3LYP functional also gives good results for most 1,3-dipoles (MAD = 1.9 kcal/mol for 17 common 1,3-dipoles), but the activation and reaction enthalpies for ozone and sulfur dioxide are difficult to calculate by any of the density functional methods.     

Cyclic (Aryl)(Amido)Carbenes: NHCs with Triplet‐like Reactivity

The synthesis and study of a library of cyclic (aryl)(amido)carbenes (CArAmCs), which represent a class of electrophilic NHCs that feature low calculated singlet‐triplet gaps (ΔE_ST=19.9 kcal mol^?1; B3LYP/def2‐TZVP) and exhibit reactivity profiles expected from triplet carbenes, are described. The electrophilic properties of the CArAmCs were quantified by analyzing their respective selenium adducts, which exhibited the largest downfield ⁷⁷Se NMR chemical shifts (up to 1645 ppm) measured for any NHC derivative known to date, as well as their Ir carbonyl complexes, from which large Tolman electronic parameter (TEP) values (up to 2064 cm^?1) were ascertained. The CArAmCs were found to engage in reactions that are typically observed with triplet carbenes, including C?H insertions, [2+1] cycloadditions with alkenes as well as alkynes, and spontaneous oxidation upon exposure to oxygen.     

Nonradical mechanisms for the uncatalyzed thermal functionalization of silicon surfaces by alkenes and alkynes: a density functional study

We propose a new concerted mechanism for the uncatalyzed hydrosilylation of terminal alkenes and alkynes, alternative to the conventional radical-based mechanism. Density functional calculations have been carried out on these and on previously proposed alternative mechanisms for the hydrosilylation of ethylene and acetylene by suitable finite size clusters as models of the thermal functionalization of -SiH3, =SiH2, and [triple bound] SiH groups in flat Si(100) and Si(111) and porous silicon surfaces by alkenes and alkynes. For each step involved in the considered hydrosilylation pathways, we optimized the geometries of reactants and products and located the corresponding transition states. The calculated activation energies for the concerted pathways of ethylene and acetylene are, respectively, 57.6 and 60.9 kcal mol(-1) on -SiH3 and in the ranges 62-63 and 58-61 kcal mol-1 on =SiH2 and 64-66 and 56-61 kcal mol(-1) on SiH. These values are much lower than the activation energies calculated for the corresponding homolytic dissociation of the Si-H bond, which is the preliminary step in the radical path, 85.6, 82-83, and 79-81 kcal mol(-1), respectively, for -SiH3, =SiH2, and [triple bound] SiH groups. Our results thus suggest that the thermal hydrosilylation of alkenes and alkynes on silicon surfaces, for which a radical-based mechanism is currently accepted, may occur through a concerted mechanism.     

Stopped-flow kinetics of tetrazine cycloadditions; experimental and computational studies toward sequential transition states

The Diels-Alder cycloadditions of tetrazines (1) with alkynes (2) are expected to give bicyclic adducts (3). Kinetic measurements of the cycloadditions of 1a and 1b with 2a give DeltaG(++) = 19.2 +/- 1.0 and 11.5 +/- 1.2 kcal/mol, respectively. Stopped-flow UV studies on the reaction of 1b with 2a show an isosbestic point at 428 nm; this places an upper limit of 11.6 +/- 2.6 kcal/mol on DeltaG(++) for loss of N(2) from the putative bicyclic intermediate 3b. Calculations (B3LYP/6-31G(d,p) + ZPVE) of transition structures for the reaction of tetrazinediacid 1d with propynylamine 2c are consistent with the experimental results for the reaction of 1b with 2a. This and several related model systems reveal two interesting features of the calculated energy surfaces. First, there may be no barrier for the loss of nitrogen from structures 3 and thus there may be two sequential transition states. This also extends Berson's correlation of activation energy with reaction energy in pericyclic reactions to significantly lower barriers. Second, for the cycloaddition of 4e and 2c, there is neither an intermediate nor a transition state between TS3e and the final product 6e. It appears that the energy surface "turns a corner" in the vicinity of a structure resembling 5e. This is not a mathematically well-defined point but has chemical consequences in that the overall exothermicity of the reaction from 4e to 6e is not felt in TS3e.     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Temperature dependence of polar and steric effects in CCl2 cycloadditions

In cycloadditions of CCl₂ activation enthalpies increase with bulky substituents at the alkenes (Table I, Figure III), whereas by polar groups only activation entropies are changed.     

Theoretical insight into stereoselective reaction mechanisms of 2,4-pentanediol-tethered ketene-olefin [2 + 2] cycloaddition

We report ab initio molecular dynamics calculations based on density functional theory performed on an intramolecular [2 + 2] cycloaddition between ketene and olefin linked with a 2,4-pentanediol (PD) tether. We find that the encounter of the ketene and olefin moieties could be prearranged in the thermal equilibrated state before the cycloaddition. The reaction mechanism is found to be stepwise, similar to that of intermolecular ketene [2 + 2] cycloadditions with ordinary alkenes. A distinct feature of the reaction pathway for a major diastereoisomer is a differential activation free energy of about 1.5 kcal/mol, including 2.8 kcal/mol as the differential activation entropy, with a transition state consisting of a flexible nine-membered ring in the olefin-PD-ketene moiety. This theoretical study provides a reasonable explanation for the strict stereocontrollability of the PD-tethered ketene-olefin cycloaddition, irrespective of reaction types or conditions.     

An unusual exchange between alkylidyne alkyl and bis(alkylidene) tungsten complexes promoted by phosphine coordination: kinetic, thermodynamic, and theoretical studies

d0 Tungsten alkylidyne alkyl complex (Me3SiCH2)3W(CSiMe3)(PMe3) (4a) was found to undergo a rare, PMe3-promoted exchange with its bis(alkylidene) tautomer (Me3SiCH2)2W(=CHSiMe3)2(PMe3) (4b). Thermodynamic studies of the exchange showed that 4b is favored and gave Keq and the enthalpy and entropy of the equilibrium: DeltaH degrees = -1.8(0.5) kcal/mol and DeltaS degrees = -1.5(1.7) eu. Kinetic studies of the alpha-H migration between 4a and 4b by variable-temperature NMR gave rate constants k1 and k-1 for the reversible reactions and activation enthalpies and entropies: DeltaH1 = 16.2(1.2) kcal/mol and DeltaS1 = -22.3(4.0) eu for the forward reaction (4a --> 4b); DeltaH2 = 18.0(1.3) kcal/mol and DeltaS2 = -20.9(4.3) eu for the reverse reaction (4b --> 4a). Ab initio calculations at the B3LYP level revealed that PMe3 binds with the bis(alkylidene) tautomer relatively more strongly than with the alkylidyne tautomer and thus stabilizes the bis(alkylidene) tautomer.     

On the lack of conjugation stabilization in polyynes (polyacetylenes)

By analogy to conjugated polyenes, conjugative stabilization of polyynes with the -CC-CC- group might be expected to be substantial. On the contrary, consistent with our recent report of a surprising lack of conjugative stabilization in butadiyne, we find by G3(MP2) calculations and by comparisons with available experimental data from these and other laboratories that the ground-state stabilization of conjugated polyynes is in fact quite small, amounting to <1 kcal mol(-)(1). By similar calculations, the 2,4-pentadiyn-1-yl radical shows no enhanced stabilization relative to 2-propyn-1-yl radical, despite the potential stabilization of the odd electron by two conjugated triple bonds and unlike the behavior of 2,4-pentadien-1-yl radical. The thermochemistry of straight-chain alkynes and polyynes is very self-consistent. Enthalpies of hydrogenation, leading to enthalpies of formation, are predictable with a high degree of accuracy (absolute mean deviation = +/-0.39 kcal mol(-)(1) vs theoretical values and +/-0.52 vs experimental) from three molecular structure enthalpies and one conjugation stabilization parameter.     

Derivatization for electrospray ionization-mass spectrometry. 4. Alkenes and alkynes

A survey of derivatization strategies and prospective derivatization reactions for conversion of simple alkenes and alkynes to 'electrospray-active' species is presented. General synthetic strategies are discussed and illustrative examples of prospective derivatives prepared from model compounds are presented along with their electrospray ionization (ES) mass spectra. The identified derivatives of these neutral, nonpolar analytes are either ionic or are ionizable in solution through Bronsted acid/base chemistry, by Lewis acid/base chemistry, or by chemical or electrochemical electron-transfer chemistry. Once ionized, the derivatives are expected to be amenable to detection by electrospray ionization-mass spectrometry. Derivatives are identified for positive and negative ion analysis of both alkenes and alkynes. Copyright 2000 John Wiley & Sons, Ltd.     

Relative reactivity of peracids versus dioxiranes (DMDO and TFDO) in the epoxidation of alkenes. A combined experimental and theoretical analysis

Comparative analysis of the calculated gas-phase activation barriers (DeltaE++) for the epoxidation of ethylene with dimethyldioxirane (DMDO) and peroxyformic acid (PFA) [15.2 and 16.4 kcal/mol at QCISD(T)// QCISD/6-31+G(d,p)] and E-2-butene [14.3 and 13.2 kcal/mol at QCISD(T)/6-31G(d)//B3LYP/6-311+G(3df,2p)] suggests similar oxygen atom donor capacities for both oxidants. Competition experiments in CH(2)Cl(2) solvent reveal that DMDO reacts with cyclohexene much faster than peracetic acid/acetic acid under scrupulously dried conditions. The rate of DMDO epoxidation is catalyzed by acetic acid with a reduction in the classical activation barrier of 8 kcal/mol. In many cases, the observed increase in the rate for DMDO epoxidation in solution may be attributed to well-established solvent and hydrogen-bonding effects. This predicted epoxidative reactivity for DMDO is not consistent with what has generally been presumed for a highly strained cyclic peroxide. The strain energy (SE) of DMDO has been reassessed and its moderated value (about 11 kcal/mol) is now more consistent with its inherent gas-phase reactivity toward alkenes in the epoxidation reaction. The unusual thermodynamic stability of DMDO is largely a consequence of the combined geminal dimethyl- and dioxa-substitution effects and unusually strong C-H and C-CH(3) bonds. Methyl(trifluoromethyl)dioxirane (TFDO) exhibits much lower calculated activation barriers than DMDO in the epoxidation reaction (the average DeltaDeltaE++ values are about 7.5 kcal/mol). The rate increase relative to DMDO of approximately 10(5), while consistent with the higher strain energy for TFDO (SE approximately 19 kcal/mol) is attributed largely to the inductive effect of the CF(3) group. We have also examined the effect of alkene strain on the rate of epoxidation with PFA. The epoxidation barriers are only slightly higher for the strained alkenes cyclopropene (DeltaE++ = 14.5 kcal/mol) and cyclobutene (DeltaE++ = 13.7 kcal/mol) than for cyclopentene (DeltaE++ = 12.1 kcal/mol), reflecting the fact there is little relief of strain in the transition state. Alkenes strained by twist or pi-bond torsion do exhibit much lower activation barriers.     

On the nature of the unusually long OO bond in HO(3) and HO(4) radicals

The HO(3) and HO(4) polyoxide radicals have attracted some attention due to their potential role in ozone chemistry. Experimentally, the geometrical structure of HO(3) is known whereas that of HO(4) is not. Moreover, the existence of the latter radical has been questioned. Theoretical calculations on the two species have been reported before, showing important structural differences depending on the computational level. Both radicals present an unusually long OO bond (around 1.7-1.8 A) that can be associated with an intricate interaction between HO, or HO(2), with O(2). The nature of such interaction is investigated in detail using large scale ab initio methods (CASSCF, CASPT2, MRCI, QCISD, CCSD(T)) and density functional techniques (B3LYP) in connection with extended basis sets. Stabilization enthalpies at 298 K with respect to HO (or HO(2)) and O(2) have been computed amounting to -3.21 kcal mol(-1) for HO(3) (trans conformation) and 11.33 kcal mol(-1) for HO(4) (cis conformation). The corresponding formation enthalpies are 6.12 and 11.83 kcal mol(-1). The trans conformation of HO(4) is less stable than the cis one by 6.17 kcal mol(-1). Transition states for HO(4) dissociation and for cis/trans conversion are also described.