A first-principles study on the adsorption behaviour of methanol and ethanol over C59B heterofullerene

In the present study, the adsorption behaviour of methanol (CH₃OH) and ethanol (C₂H₅OH) molecules over heterofullerene C₅₉B surface is studied by density functional theory calculations. This heterofullerene is obtained from C₆₀ by substituting a carbon atom with a boron atom and relaxing self-consistently the structure to the local minimum. The adsorption of CH₃OH and C₂H₅OH on the C₅₉B is exothermic and the relaxed geometries are stable. The CH₃OH and C₂H₅OH adsorption can also induce a change in the highest occupied molecular orbital and the lowest unoccupied molecular orbital energy gap of the nanocage. The dehydrogenation pathways of CH₃OH and C₂H₅OH via O–H and C–H bonds scission are also examined. The results indicate that O–H bond scission is the most favourable pathway on the C₅₉B surface.     

In situ XPS and MS study of methanol decomposition and oxidation on Pd(111) under millibar pressure range

The methanol decomposition and oxidation on a Pd(111) single crystal have been investigated in situ using ambient-pressure X-ray photoelectron spectroscopy (XPS) and mass-spectrometry (MS) in the temperature range of 300–600 K. It was found that even in the oxygen presence the methanol decomposition on palladium proceeds through two competitive routes: fast dehydrogenation to CO and H₂, and slow decomposition of methanol via the C–O bond scission. The rate of the second route is significant even in the millibar pressure range, which leads to a blocking of the palladium surface by carbon and to a prevention of the further methanol conversion. As a result, no gas phase products of methanol decomposition were detected by mass-spectrometry at 0.1 mbar CH₃OH in the whole temperature range. The methanol C–O bond scission produces CH_x species, which fast dehydrogenate to atomic carbon even at room temperature and further partially dissolve in the palladium bulk at 400 K with the formation of the PdC_x phase. According to in situ XPS data, the PdC_x phase forms even in the oxygen excess. The application of an in situ XPS–MS technique unambiguously shows a good correlation between a decrease in the surface concentration of all carbon-containing species and the rate of methanol conversion. Since these carbon species have a high reactivity towards oxygen, heating of Pd(111) above 450 K in a methanol–oxygen mixture yields CO, CO₂, and water. The product distribution indicates that the main route of methanol conversion is the dehydrogenation of methanol to CO and hydrogen. However, under the experimental conditions used, hydrogen is completely oxidized to water, while CO is partially oxidized to CO₂. No palladium oxide was detected by XPS in these conditions.     

Activation of methyl acetate on Pd(111)

The adsorption and activation of methyl acetate (CH₃COOCH₃), one of the simplest carboxylic esters, on Pd(111) have been studied using self-consistent periodic density functional theory calculations. Methyl acetate adsorbs weakly through the carbonyl oxygen. Its activation occurs via dehydrogenation, instead of direct C–O bond dissociation, on clean Pd(111): It is much more difficult to dissociate the C–O bonds (E_a ≈ 2.0 eV for the carbonyl and acetate–methyl bonds; E_a = 1.0 eV for the acetyl–methoxy bond) than to dissociate the C–H bonds to produce enolate (CH₂COOCH₃; E_a = 0.74 eV) or methylene acetate (CH₃COOCH₂; E_a = 0.82 eV). The barriers for C–H and C–O bond dissociation are directly calculated for enolate and methylene acetate, and estimated for further dehydrogenated derivatives (CH₃COOCH, CH₂COOCH₂, and CHCOOCH₃) based on the Brønsted–Evans–Polanyi linear energy relations formed by the calculated steps. The enolate pathway leads to successive dehydrogenation to CCOOCH₃, whereas methylene acetate readily dissociates to yield acetyl. The selectivity for dissociating the acyl–alkoxy C–O bond, which is desired for alcohol formation, is therefore fundamentally limited by the facility of dehydrogenation under vacuum/low-pressure conditions on Pd(111).     

Thermal evolution of acetylene and ethylene on Pd(111)

The thermal evolution of acetylene and ethylene and their deuterated counterparts on a palladium (111) surface has been studied by high-resolution electron energy loss spectroscopy in the temperature range 150–500 K. Analysis of the vibrational spectra indicates that chemisorbed acetylene evolves at 300 K in the presence of surface hydrogen to mainly ethylidyne, CCH₃, and a small amount of residual acetylene. Spectra obtained with and without preadsorbed hydrogen provide evidence for a 〉C CH₂ intermediate in the reaction. Chemisorbed ethylene also evolves to ethylidyne after heating from 150 to 300 K but much of the ethylene desorbs. The high temperature (400–500 K) behavior of C₂H₂ and C₂H₄ involves formation of a CH species. Although a small amount of the CH species may be formed from the dehydrogenation of ethylidyne, it is found that carbon-carbon bond scission of acetylene near 400 K is the dominant mechanism in CH formation.     

Hydroxyl end groups influence in vibrational and transport properties in polymer/monomer solutions: the PEO/EG case

A study has been made of vibrational properties in ethylene glycol (EG; H(OCH₂CH₂)OH) and EG monomethyl ether (EGmE; CH₃(OCH₂CH₂)OH) in solution together with poly(ethylene oxide) (PEO; H(OCH₂CH₂)_n,OH) at different concentrations, performed by Fourier transform infrared absorbance (FT-IR) spectroscopy. The results ae compared with previous viscometry and photon correlation spectroscopy (PCS) studies, using EG dimethyl ether (EGdE; CH₃(OCH₂CH₂)OCH₃) as solvent as well. These homologous systems differ from each other in the number of OH end groups, in particular two for EG, one for EgmE and zero for EGdE. Combining analysis of the vibrational and transport properties of EG, EGmE and EGdE in solution with PEO over a wide range of concentration made it possible to check the quality (good theta or poor) of these three different solvents and the role played by the hydrogen bond in the various solute-solvent interaction mechanisms, resulting in the well known de Gennes scaling law.     

Comparison of O–H and C–H activation of methanol on Ni-based cluster: a DFT investigation

Using the density functional theory, the initial dehydrogenation of methanol on Ni_xM_y (M?=?Ni, Co, Fe, Mn, Cr, x?+?y?=?4, y?=?1, 2) clusters is investigated. Two adsorption and dehydrogenation mechanisms of methanol are studied: one proceeds along the C–H scission and another begins with the breaking of the O-H bond. The adsorption sites of methanol on the Ni or M sites of the Ni_xM_y clusters are considered. The adsorption of methanol on Ni₄ cluster is stronger than those on bimetallic clusters, while the initial dehydrogenation barriers on Ni_xM_y clusters are lower than that on Ni₄ cluster. The comparable energy barriers of two pathways (O–H or C–H dissociation) on Ni-based clusters indicate that these two paths are quite competitive. In addition, the Ni₂M₂ clusters show superior activation performance compared with the Ni₃M clusters, especially for Ni₂Mn₂ and Ni₂Cr₂ clusters. The effects of alloyed metal on the catalytic activity of Ni for methanol initial dehydrogenation, including the adsorption energy, O–H or C–H bond scission barrier and frontier molecular orbital levels, are discussed. It can be concluded that the addition of Co, Fe, Mn and Cr to Ni catalyst is able to enhance the activity of the methanol dehydrogenation reaction.     

Adsorption and decomposition of γ-butyrolactone on Pd(1 1 1) and Pt(1 1 1)

The adsorption and thermal chemistry of γ-butyrolactone (GBL) on the (1 1 1) surface of Pd and Pt has been investigated using a combination of high resolution electron energy loss spectroscopy (HREELS) and temperature programmed desorption (TPD). HREELS results indicate that GBL adsorbs at 160 K on both surfaces through its oxygenate functionality. On Pd(1 1 1), adsorbed GBL undergoes ring-opening and decarbonylation by 273 K to produce adsorbed CO and surface hydrocarbon species. On Pt(1 1 1), very little dissociation is observed using HREELS, with almost all of the GBL simply desorbing. TPD results are consistent with decarbonylation and subsequent dehydrogenation reactions on Pd(1 1 1), although small amounts of CO₂ are also detected. TPD results from Pt(1 1 1) indicate that a small proportion of adsorbed GBL (perhaps on defect sites) does undergo ring-opening to produce CO, CO₂, and H₂. These results suggest that the primary dissociation pathway for GBL on Pd(1 1 1) is through O-C scission at the carbonyl position. Through comparisons with previously published studies of cyclic oxygenates, these results also demonstrate how ring strain and functionality affect the ring-opening rate and mechanism.     

Electron impact spectra of some alkyl derivatives of water and related compounds

Electron impact energy loss spectra at impact energies of 100 eV and 50 eV and a scattering angle of 2 degrees are presented for alkyl derivatives of H₂O and related compounds. Spectra of H₂O, CH₃OH, CH₃OCH₃, ethylene oxide, CH₃CH₂OH, (CH₃)₂CHOH, (CH₃)₃C C₂H₅OC₂H₅ and tetrahydrofuran are tentatively assigned using derived term values and ionization potentials from photoelectron spectrosc Substituent effects on Rydberg orbital energies are discussed using Taft σ* values.     

Atmospheric oxidation chemistry of 1–methoxy 2–propyl acetate initiated by OH radicals: kinetics and mechanisms

Kinetics and mechanism of the gas-phase reaction of CH₃C(O)OCH(CH₃)CH₂OCH₃ (MPA) with OH radicals in the presence of O₂ and NO have been investigated theoretically by performing a high and reliable level of theory, viz., CCSD(T)/6-311?+?G(d,p)//BH&HLYP/6-311++G(d,p)?+?0.9335×ZPE. The calculations predict that the H-abstraction from the ?CH₂?O? position of MPA is the most facile channel, which leads to the formation of the corresponding alkoxy radicals CH₃C(O)OCH(CH₃)C(O ?)HOCH₃ under atmospheric conditions. This activated radicals CH₃C(O)OCH(CH₃)C(O ?)HOCH₃ will undergo further rearrangement, fragmentation and oxidative reactions and predominantly leads to the formation of various products (methyl formate HC(O)OCH₃ and acetic anhydride CH₃C(O)OC(O)CH₃). In the presence of water, acetic anhydride can convert into acetic acid CH₃C(O)OH via the hydrolysis reaction. The calculated total rate constants over the temperature range 263–372?K are used to derive a negative activation energy (E_a= ?5.88 kJ/mol) and an pre-exponential factor (A?=?1.78×10^?12 cm³ molecule^?1 s^?1). The obtained Arrhenius parameters presented here are in strong agreement with the experimental values. Moreover, the temperature dependence of the total rate constant over a temperature range of 263?1000?K can be described by k?=?5.60 × 10^?14×(T/298?K)^3.4×exp(1725.7?K/T) cm³ molecule^?1 s^?1.     

A comprehensive study of (CH3)2CHOC(S)SC(O)OCH3 using matrix isolation technique,X‐ray analysis,spectroscopic studies and theoretical calculations

Potassium isopropyl xanthate, (CH₃)₂CHOC(S)SK, reacts with methyl chloroformiate ClC(O)OCH₃ to yield (methoxycarbonyl) (2‐propoxythiocarbonyl) sulfide, (CH₃)₂CHOC(S)SC(O)OCH₃. This novel xanthogen formate was characterized by ¹H and ¹³C{¹H} NMR spectroscopy, mass spectrometry and IR and Raman spectroscopy. The structure of a single crystal of (CH₃)₂CHOC(S)SC(O)OCH₃ was determined by X‐ray diffraction analysis at 173 K. The conformational properties have been studied by liquid IR and Raman spectroscopy, matrix isolation spectroscopy together with photochemical studies and quantum chemical calculations (HF and B3LYP methods with the 6‐31+G* basis set). The analysis of the IR spectrum of liquid (CH₃)₂CHOC(S)SC(O)OCH₃ suggests the presence of two conformers in equilibrium at room temperature. However, in the photochemical matrix study, an equilibrium of three conformers was detected. These forms were further characterized by theoretical calculations. Different photolysis products, such as CH₃OC(O)SCH(CH₃)₂, OCS, CO, CO₂ and CS₂, were identified by matrix spectroscopy. The IR absorptions of CH₃OC(O)SCH(CH₃)₂, for which literature data are scarce, were analysed in the light of the results of appropriate theoretical calculations. Copyright © 2009 John Wiley & Sons, Ltd.     

Determination of the interaction using FTIR within the components of GPE based on the synthesized grafted polymer

An oxyethylene chain (–CH₂–CH₂–O–) grafted polymer (P(MMA–MAh)–PEGME) was synthesized by reacting poly(ethylene glycol) monomethyl ether (PEGME) with the copolymer of poly(methyl metacrylate–maleic anhydride) (P(MMA–MAh)) and endcapping the residual carboxylic acid with methanol. Rectorite modified with dodecyl benzyl dimethyl ammonium chloride (OREC) was used as a filler additive to modify gel polymer electrolytes (GPE) that consisted of P(MMA–MAh)–PEGME used as polymer matrix, propylene carbonate (PC) as plasticizer and LiClO₄ as lithium ion producer. Characterization of interaction of the polar group of C=O or C–O–C in PC and grafted polymer with Li⁺ and OH group on OREC surface has been thoroughly examined using Fourier transform infrared (FTIR), respectively. The quantitative analysis of FTIR shows that the absorptivity coefficients (a) of polymer/LiClO₄, PC/LiClO₄, PC/OREC and polymer/OREC are 0.705, 0.113, 0.430 and 0.500, respectively, which means that the Li⁺ or OH bonded polar group of C=O and C–O–C is more sensitive than free C=O and C–O–C in FTIR spectra. The limited values of bonded C=O and C–O–C equivalent fraction of polymer/LiClO₄, PC/LiClO₄, PC/OREC and polymer/OREC are 33%, 94%, 60% and 22%, respectively, which implies that the interaction within the components is reversible and the intensity of interaction is ordered as PC/LiClO₄, PC/OREC, polymer/OREC and polymer/LiClO₄.     

DFT study on complete ethylene decomposition on flat and stepped Pd

We applied density functional theory (DFT) calculations to study ethylidyne (CCH₃) adsorption and decomposition to C and H over flat and stepped Pd surfaces. Our calculations show that ethylidyne is the most stable molecule among all the possible dehydrogenation or decomposition residues of ethylene. We discuss various possible reaction pathways for ethylidyne decomposition and point out that the most probable one is via ethynyl (CCH) species suggested also by experimental observations. Our calculations indicate that the presence of steps modify the potential energy surface by increasing the binding of most of the species, and also lowering the activation barrier for several reactions. Furthermore we show that the energetics related to dehydrogenation of ethylene and its derivatives manifest a Brønsted–Evans–Polanyi type of behavior.     

Characterization of methylidyne on Pt(1 1 1) with infrared spectroscopy

Methylidyne (CH) was prepared on Pt(1 1 1) by three methods: thermal decomposition of diiodomethane (CH₂I₂), ethylene decomposition at temperatures above 450 K, and surface carbon hydrogenation. Methylidyne and its precursors are characterized by reflection absorption infrared spectroscopy (RAIRS). The C-I bond of diiodomethane breaks upon adsorption to produce methylene (CH₂), which decomposes to methylidyne at temperatures above 130 K. Above 200 K, methylidyne is the only hydrocarbon species observed with RAIRS, although reaction channels for the formation of methane (CH₄) and ethylene (C₂H₄) are indicated by temperature programmed desorption (TPD). As is well known from numerous previous studies, ethylene decomposes to ethylidyne (CCH₃) upon exposure to Pt(1 1 1) at 410 K. Upon annealing to 450 K, ethylidyne dissociates through two reaction pathways, dehydrogenation to ethynyl (CCH) and C-C bond scission to methylidyne. Ethylene dehydrogenation on the surface at 750 K and under low ethylene exposures produces surface carbon that can be hydrogenated to methylidyne with C-H and C-D stretch frequencies of 2956 and 2206 cm⁻¹, respectively. Hydrogen co-adsorption on the surface causes these frequencies to shift to higher values. Methylidyne is stable on Pt(1 1 1) to temperatures up to 500 K.     

Bond activation sequence observed in the chemisorption and surface reaction of ethanol on Ni(111)

The mechanism of ethanol decomposition on the Ni(111) surface has been investigated between 155 and 500 K. The sequence of bond scission steps which occur as ethanol undergoes dissociative reactions on this surface has been deduced using deuterium and ¹³C isotopic labels. Bond activation occurs in the order (1) OH, (2) CH₂ (methylene CH), (3) CC, (4) CH₃ (methyl CH). The products observed are CH₃CHO(g), CH₄(g), CO(g), H₂(g) and surface carbon, C(a). The latter species exhibits a carbidic AES lineshape in the temperature range 450 to 670 K, at which temperature it dissolves into the Ni bulk. Acetaldehyde, CH₃CHO, and methane, CH₄, desorb with the same threshold temperature (260–265 K), and the formation of both of these products is controlled by scission of the methylene CH bond (CH₂ group). The CH₃ group is cleaved from the intermediate surface CH₃CHO species to form CH₃(ads). H₂ exhibits a broad, doublet desorption peak from 300 to 450 K. The carbonoxygen bond in ethanol remains intact and CO ultimately desorbs in a single desorption limited process (T_p = 430 K). A small fraction of CO(a) species undergo exchange with the carbidic surface carbon in a minor process observed above 440 K.     

Theoretical study on the reactions of CH3NHNH2 with ground state O(3P) atom and excited state O(1D) atom

The reaction mechanisms of methylhydrazine (CH₃NHNH₂) with O(³P) and O(¹D) atoms have been explored theoretically at the MPW1K/6-311+G(d,p), MP2/6-311+G(d,p), MCG3-MPWPW91 (single-point), and CCSD(T)/cc-pVTZ (single-point) levels. The triplet potential energy surface for the reaction of CH₃NHNH₂ with O(³P) includes seven stable isomers and eight transition states. When the O(³P) atom approaches CH₃NHNH₂, the heavy atoms, namely N and C atoms, are the favourable combining points. O(³P) atom attacking the middle-N atom in CH₃NHNH₂ results in the formation of an energy-rich isomer (CH₃NHONH₂) followed by migration of O(³P) atom from middle-N atom to middle-H atom leading to the product P6 (CH₃NNH₂+OH), which is one of the most favourable routes. The estimated major product CH₃NNH₂ is consistent with the experimental measurements. Reaction of O(¹D) + CH₃NHNH₂ presents different features as compared with O(³P) + CH₃NHNH₂. O(¹D) atom will first insert into C–H2, N1–H4, and N2–H5 bonds barrierlessly to form the three adducts, respectively. There are two most favourable paths for O(¹D) + CH₃NHNH₂. One is that the C–N bond cleavage accompanied by a concerted H shift from O atom to N atom (mid-N) leads to the product P_I (CH₂O + NH₂NH₂), and the other is that the N–N bond rupture along with a concerted H shift from O to N (end-N) forms P_IV (CH₃NH₂ + HNO). The similarities and discrepancies between two reactions are discussed.     

Infrared and Raman spectroscopic investigation of Eu3+-doped mono and Di-urethanesil hybrid siliceous materials

The general features of two series of sol-gel derived materials, designatedurethanesils (Ut), have been investigated by infrared and Raman spectroscopies with the goal of elucidating the chemical environment of the Eu³⁺ cations. The host frameworks of the two families of ormolytes studied have been represented by m-Ut(350) and d-Ut(300), where m stands for mono, d stands for di, 350 and 300 are the average molecular weights of the organic precursors (poly(ethylene glycol) methyl ether, PEGME, and poly(ethylene glycol), PEG, respectively). The hybrid matrix of the mono-xerogels is composed by a siliceous backbone bonded by means of urethane linkages (-NHC (=O)O-) to pendant methyl end capped oligopolymer chains with approximately 7 oxyethylene units, whereas that of the di-xerogels is based on a siliceous network grafted through urethane groups to both ends of poly(oxyethylene) segments containing about 6 (OCH₂CH₂) repeat units. Both classes of materials have been doped with europium triflate (Eu(CF₃SO₃)₃). The doped samples have been identified by m-Ut(350)_nEu(CF₃SO₃)₃ and d-Ut(300)_nEu(CF₃SO₃), where n is the molar ratio of (OCH₂CH₂) repeat units per Eu³⁺ ion. Materials with n ranging from ∞ to 5 have been analyzed. The spectral data obtained provide evidence that the cations begin to coordinate to the ether oxygen atoms of the oligopolymer chains at n=40 in the mono-urethanesils and at n=10 in the di-urethanesils. In mono-urethanesils samples with n>40 and in di-urethanesils materials with n>10, the Eu³⁺ coordinate exclusively to the carbonyl oxygen atoms of the urethane linkages. Paper presented at the 8th EuroConference on Ionics, Carvoeiro, Algarve, Portugal, Sept. 16–22, 2001.     

C2H4在Ru(1010)表面吸附与分解的研究

用X射线电子能谱(XPS)、热脱附谱(TDS)和紫外光电子能谱(UPS)方法研究了乙烯(C_{2H4)在Ru(1010)表面的吸附,在低温下(200K以下)乙稀(C2Ｈ4})可以在Ru(1010)表面上以分子状态稳定吸附,在200K以上乙烯(C₂H ₄)则发生了脱氢分解反应.TDS结果表明乙烯(C₂H₄)分 解后的主要产物为乙炔(C< 关键词： 乙烯 钌(1010)表面 吸附与分解     

ZnO micro- and nano-structures: microwave-assisted solvothermal synthesis,morphology control and photocatalytic properties

ZnO with various morphologies have been prepared by a simple microwave-assisted solvothermal method in a short period of time. In this synthetic method, Zn(CH₃COO)₂⋅2H₂O is used as the reactant, a mixture of water and ethylene glycol is used as solvents, and no other additives are used. The morphology of the prepared samples can be controlled by changing the mixing ratio of ethylene glycol to water. The photocatalytic activities of the prepared samples are also investigated.     

Insight into the reaction mechanism of ethanol steam reforming catalysed by Co–Mo6S8

The production of hydrogen via steam reforming of ethanol (SRE) is favourable for the use of hydrogen as an alternative fuel. Co–Mo₆S₈ possesses high activity and stability for SRE to sustainably produce hydrogen. The competition among reaction pathways related to C–H, O–H, C–C, C–O cleavage and H₂ formation was studied. The adsorption and reaction of related intermediates in the ESR reaction pathway are described. The results indicated that the most feasible route for the decomposition of ethanol catalysed by Co–Mo₆S₈ is CH₃CH₂OH*→CH₃CH₂O*→CH₃CHO*→CH₂CHO*→CHCHO*→CHCO*→CH*+CO*. The CH* can be decomposed into C*+H*, and CO* can be oxidised via the redox mechanism of the water gas shift (WGS) reaction. Thus the final products are CO₂ and H₂. The present result may help people to design an SRE catalyst, which has the ability to break C–C to form CO and H₂, then CO react with H₂O in the WGS reaction generating CO₂ and H₂.     

Thermochemistry,bond energies and internal rotor barriers of methyl sulfinic acid,methyl sulfinic acid ester and their radicals

Sulfur–Oxygen containing hydrocarbons are formed in oxidation of sulfides and thiols in the atmosphere, on aerosols and in combustion processes. Understanding their thermochemical properties is important to evaluate their formation and transformation paths. Structures, thermochemical properties, bond energies, and internal rotor potentials of methyl sulfinic acid CH₃S(?O)OH, its methyl ester CH₃S(?O)OCH₃ and radicals corresponding to loss of a hydrogen atom have been studied. Gas phase standard enthalpies of formation and bond energies were calculated using B3LYP/6‐311G (2d, p) individual and CBS‐QB3 composite methods employing work reactions to further improve accuracy of the ${\Delta} _{{\bf f}} H_{{\bf 298}}^{{\bf o}} $ . Molecular structures, vibration frequencies, and internal rotor potentials were calculated. Enthalpies of the parent molecules CH₃S(?O)OH and CH₃S(?O)OCH₃ are evaluated as ?77.4 and ?72.7 kcal mol^?1 at the CBS? QB3 level; Enthalpies of radicals C^?H₂? S(?O)? OH, CH₃? S^?(?O)₂, C^?H₂? S(?O)? OCH₃ and CH₃? S(?O)? OC^?H₂ (CBS‐QB3) are ?25.7, ?52.3, ?22.8, and ?26.8 kcal mol^?1, respectively. The CH₃C(?O)O—H bond dissociation energy is of 77.1 kcal mol^?1. Two of the intermediate radicals are unstable and rapidly dissociate. The CH₃S(?O)? O. radical obtained from the parent CH₃? S(?O)? OH dissociates into methyl radical (${\bf CH}_{{\bf 3}}^{{\bf .}} $ ) plus SO₂ with endothermicity (ΔH_rxn) of only 16.2 kcal mol^?1. The CH₃? S(?O)? OC^?H₂ radical dissociates into CH₃? S^?=O and CH₂=O with little or no barrier and an exothermicity of ?19.9 kcal mol^?1. DFT and the Complete Basis Set‐QB3 enthalpy values are in close agreement; this accord is attributed to use of isodesmic work reactions for the analysis and suggests this combination of B3LYP/work reaction approach is acceptable for larger molecules. Copyright © 2010 John Wiley & Sons, Ltd.