Barrierless Pathway Research Articles

Damage mechanism of hydroxyl radicals toward adenine—thymine base pair

The adenine—thymine base pair was studied in the presence of hydroxyl radicals in order to probe the hydrogen bond effect. The results show that the hydrogen bonds have little effect on the hydroxylation and dehydrogenation happened at the sites, which are not involved in a hydrogen bond, while at the sites involved in hydrogen bond formation in the base pair, the reaction becomes more difficult, both in view of the free energy barrier and the exothermicity. With a 6–311++G(d,p) level of description, both B3LYP and MP2 methods confirm that the C8 site of isolated adenine has the highest possibility to form covalent bond with the hydroxyl radicals, though with different energetics: B3LYP predicts a barrierless pathway, while MP2 finds a transition state with an energy of 106.1 kJ/mol. For the dehydrogenation reactions, B3LYP method predicts that the free energy barrier increases in the order of HN9 < HN61 < HN62 < H2 < H8.

Probing the Reactivity of Hydroxyl Radicals toward Isolated Thymine Using Theoretical Calculations

B3LYP and MP2 methods are used to study the hydroxylation and dehydrogenation of thymine at different sites in the presence of hydroxyl radicals (⋅OH). With a 6-311++G(d,p) level of description, the consistent results thermodynamically show that the relative free energy increases in the order of CM5 < N1 < C6 < N3 on the hydrogen abstraction. However, toward the hydroxylation, though the two methods give a similar trend that free energy barrier increases in the order of C6 < C5 < C4 < C2, they describe the reaction at C6 site in a different way: B3LYP calculation gives a barrierless pathway, whereas MP2 calculation gives a reactant complex with an energy of 19.5 kJ/mol and a pathway with a barrier energy of 64.9 kJ/mol. The calculated results are consistent with experimental data and the preference for the C5 site is due to the facility of ⋅OH in approaching the substrate. © 2013 Wiley Periodicals, Inc.

Formation of Acene-Based Polymers: Mechanistic Computational Study

Understanding the mechanism of linear acene decomposition and its reactivity is a prerequisite for controlling the stability of acenes and their future applications. Previously, we suggested that long acenes may undergo polymerization since the polymerization product is thermodynamically more stable than the dimerization product. However, due to kinetic considerations, the most thermodynamically stable product, the polymer, might not necessarily be formed. To elucidate the situation, we investigated the mechanisms of acene polymerization computationally, using pentacene, hexacene, and heptacene as representative examples. Similarly to dimerization, acene polymerization follows a stepwise biradical pathway. Structural and steric hindrance of the polymer backbone forces acene polymerization to proceed via the less reactive noncentral benzene rings. Consequently, dimerization is always kinetically more favorable than polymerization, irrespective of acene length. Although, for long acenes starting from hexacene, both polymerization and dimerization are barrierless pathways relative to the reactants, polymerization is thermodynamically preferred for hexacene and heptacene and even more so for longer acenes (since polymerization forms four new C-C bonds while dimerization forms only two). Indeed, reinvestigation of available experimental data suggests that acene-based polymers were probably obtained experimentally previously.

Excited-State Tautomerization of Gas-Phase Cytosine

In order to investigate experimentally observed phototautomerization of gas-phase cytosine, several excited-state tautomerization mechanisms were characterized at the EOM-CCSD and TDDFT levels. All pathways that took place exclusively on the S1 surface were found to have significant barriers that were much higher than the barriers involved in radiationless decay of cytosine tautomers through conical intersections back to the ground state; tautomerization in this fashion cannot compete with radiationless relaxation. However, an alternative possibility is that the conical intersections that facilitate radiationless decay could also facilitate tautomerization. Barrierless pathways indicate that it is energetically possible that bifurcation at the conical intersections can lead to a subset of the population reaching different tautomers. This could be an explanation for the observed tautomerization of keto cytosine after exposure to low-energy UV light.

Ab initio insight into ultrafast nonadiabatic decay of hypoxanthine: keto-N7H and keto-N9H tautomers

Nonadiabatic dynamics simulations at the SA-CASSCF level were performed for the two most stable keto-N7H and keto-N9H tautomers of hypoxanthine in order to obtain deep insight into the lifetime of the optically bright S1((1)ππ*) excited state and the relevant decay mechanisms. Supporting calculations on the ground-state (S0) equilibrium structures and minima on the crossing seams of both tautomers were carried out at the MR-CIS and CASSCF levels. These studies indicate that there are four slightly different kinds of conical intersections in each tautomer, exhibiting a chiral character, each of which dominates a barrierless reaction pathway. Moreover, both tautomers reveal the ultrafast S1→ S0 decay, in which the S1 state of keto-N9H in the gas phase has a lifetime of 85.5 fs, whereas that of keto-N7H has a longer lifetime of 137.7 fs. An excellent agreement is found between the present results and the experimental value of 130 ± 20 fs in aqueous solution (Chen and Kohler, Phys. Chem. Chem. Phys., 2012, 14, 10677-10689).

Nonadiabatic ab initio molecular dynamics of photoisomerization in bridged azobenzene

The photoisomerization mechanisms of bridged azobenzene are investigated by means of surface hopping dynamics simulations based on the Zhu-Nakamura theory. In the geometry optimizations and potential energy surface calculations, four minimum-energy conical intersections between the ground state and the lowest excited state are found to play important roles in the trans-cis and cis-trans isomerization processes. The trans-cis photoisomerization proceeds through two minimum-energy conical intersections. Ultrafast pedal motion of the N atoms and twisting of phenyl rings around their N-C bonds allows the molecule to move to a minimum-energy conical intersection, after which surface hopping from S(1) to S(0) occurs. In the S(0) state, further rotation occurs around the N=N bond and two N-C bonds until the azo moiety and phenyl rings complete their isomerization. Finally, the cis form is achieved by subsequent adjustment of the ethylene bridge. In the cis-trans photodynamics, there is one rotational pathway, in the middle of which two CIs are responsible for the surface hopping to the S(0) state. After the nonadiabatic transition, the molecule reaches the trans form through a barrierless pathway and the two phenyl rings and the additional bridge complete their reorientation almost at the same time.

Some new reaction pathways for the formation of cytosine in interstellar space – A quantum chemical study

The detection of nucleic acid bases in carbonaceous meteorites suggests that their formation and survival is possible outside of the Earth. Small N-heterocycles, including pyrimidine, purines and nucleobases, have been extensively sought in the interstellar medium. It has been suggested theoretically that reactions between some interstellar molecules may lead to the formation of cytosine, uracil and thymine though these processes involve significantly high potential barriers. We attempted therefore to use quantum chemical techniques to explore if cytosine can possibly form in the interstellar space by radical–radical and radical-molecule interaction schemes, both in the gas phase and in the grains, through barrier-less or low barrier pathways. Results of DFT calculations for the formation of cytosine starting from some of the simple molecules and radicals detected in the interstellar space are being reported. Global and local descriptors such as molecular hardness, softness and electrophilicity, and condensed Fukui functions and local philicity indices were used to understand the mechanistic aspects of chemical reaction. The presence and nature of weak bonds in the molecules and transition states formed during the reaction process have been ascertained using Bader’s quantum theory of atoms in molecules (QTAIMs). Two exothermic reaction pathways starting from propynylidyne (CCCH) and cyanoacetylene (HCCCN), respectively, have been identified. While the first reaction path is found to be totally exothermic, it involves a barrier of 12.5kcal/mol in the gas phase against the lowest value of about 32kcal/mol reported in the literature. The second path is both exothermic and barrier-less. The later has, therefore, a greater probability of occurrence in the cold interstellar clouds (10–50K).

Unraveling Ultrafast Dynamics in Photoexcited Aniline

A combination of ultrafast time-resolved velocity map imaging (TR-VMI) methods and complete active space self-consistent field (CASSCF) ab initio calculations are implemented to investigate the electronic excited-state dynamics in aniline (aminobenzene), with a perspective for modeling (1)πσ* mediated dynamics along the amino moiety in the purine derived DNA bases. This synergy between experiment and theory has enabled a comprehensive picture of the photochemical pathways/conical intersections (CIs), which govern the dynamics in aniline, to be established over a wide range of excitation wavelengths. TR-VMI studies following excitation to the lowest-lying (1)ππ* state (1(1)ππ*) with a broadband femtosecond laser pulse, centered at wavelengths longer than 250 nm (4.97 eV), do not generate any measurable signature for (1)πσ* driven N-H bond fission on the amino group. Between wavelengths of 250 and >240 nm (<5.17 eV), coupling from 1(1)ππ* onto the (1)πσ* state at a 1(1)ππ*/(1)πσ* CI facilitates ultrafast nonadiabatic N-H bond fission through a (1)πσ*/S(0) CI in <1 ps, a notion supported by CASSCF results. For excitation to the higher lying 2(1)ππ* state, calculations reveal a near barrierless pathway for CI coupling between the 2(1)ππ* and 1(1)ππ* states, enabling the excited-state population to evolve through a rapid sequential 2(1)ππ* → 1(1)ππ* → (1)πσ* → N-H fission mechanism, which we observe to take place in 155 ± 30 fs at 240 nm. We also postulate that an analogous cascade of CI couplings facilitates N-H bond scission along the (1)πσ* state in 170 ± 20 fs, following 200 nm (6.21 eV) excitation to the 3(1)ππ* surface. Particularly illuminating is the fact that a number of the CASSCF calculated CI geometries in aniline bear an exceptional resemblance with previously calculated CIs and potential energy profiles along the amino moiety in guanine, strongly suggesting that the results here may act as an excellent grounding for better understanding (1)πσ* driven dynamics in this ubiquitous genetic building block.

Structure, energetics, and reactions of alkali tetramers

Electronic structure calculations have been carried out for all possible alkali tetramers that can be formed from X(2) + X(2) → X(2)X(2), X(2) + Y(2) → X(2)Y(2), and XY + XY → X(2)Y(2) alkali dimer association reactions. Vibrationally stable rhombic (D(2h)) and planar (C(s)) structures are found for all possible tetramers formed from the alkali metals, Li to Cs. All tetramer formation reactions (from ground state singlet homonuclear or heteronuclear dimers) are found to be exothermic with binding energies ranging from 6282 cm(-1) for Li(2)Li(2) to 1985 cm(-1) for Cs(2)Cs(2). Extensive calculations, carried out at long-range for several reactant pairs, indicate that there are barrier-less pathways for the formation of tetramers from dimer association reactions. At low temperatures, direct formation of tetramers is unlikely, owing to the large exothermicity associated with these association reactions, but atom exchange reactions (X(2) + Y(2) ↔ XY + XY) are possible for some species.

Essential on the Photophysics and Photochemistry of the Indole Chromophore by Using a Totally Unconstrained Theoretical Approach.

Indole is a chromophore present in many different molecules of biological interest, such as the essential amino acid tryptophan and the neurotransmitter serotonin. On the basis of CASPT2//CASSCF quantum chemical calculations, the photophysical properties of the system after UV irradiation have been studied through the exploration of the potential energy hypersurfaces of the singlet and triplet low-lying valence excited states. In contrast to previous studies, the present work has been carried out without imposing any restriction to the geometry of the molecule (C1 symmetry) and by performing minimum energy path calculations, which is the only instrument able to provide the lowest-energy evolution of the system. Relevant findings to the photophysics of bare indole have been obtained, which compete with the currently accepted mechanism for the energy decay in the molecule. The results show the presence of a conical intersection (CI) between the initially populated (1)(La ππ*) and the (1)(Lb ππ*) state, easily accessible through a barrierless pathway from the Franck-Condon region. At this CI region, part of the population is switched from the bright (1)(La ππ*) state to the (1)(Lb ππ*) state, and the system evolves toward a minimum structure from which the expected fluorescence takes place. The reported low values of the fluorescence quantum yield are explained by means of a new nonradiative mechanism specific for the (1)(Lb ππ*) state, in which the presence of an ethene-like CI between the (Lb ππ*) and ground states is the main feature.

Tuning of switching properties and excited-state dynamics of fulgides by structural modifications

The ultrafast photo-induced dynamics of the E-isomers of four selected photochromic fulgides with distinct structural motifs have been elucidated by femtosecond broadband transient absorption spectroscopy in n-hexane as solvent. E→C and E→Z isomerisations, respectively, with time constants of ∼0.12 ± 0.02 ps and ∼0.34 ± 0.03 ps taking place in parallel were found for derivatives with a methyl substituent at the central hexatriene (HT) unit. In contrast, fulgides with increased steric constraints by an iso-propyl substituent or by intramolecular bridging displayed virtually zero E→Z isomerisation, but instead a desired accelerated and more efficient ring closure in a reaction time of only ∼50 ± 10 fs. Both photoisomerisations appear to follow excited-state pathways with distinctive conical intersections. For the ring closure, direct barrierless pathways with steep downhill gradients are likely. Furthermore, the results indicate conformer-specific reactions, with ring closure exclusively by the E(α) conformer and E→Z isomerisation predominantly by the E(β) conformer, because the E(α)→Z channel is unfavoured by the faster and kinetically more competitive E(α)→C reaction. DFT calculations of the equilibrium structures showed that the sterically demanding groups at the HT unit shift the conformer equilibria towards the E(α) conformers. At the same time, they appear to cause a favourable pre-orientation of the furyl unit that accelerates the conrotatory ring closure in the E(α)→C reaction. Benzo-annulation of the furyl unit has little effect on the observed dynamics. Overall, the results demonstrate how the excited-state dynamics and thereby the photoswitching properties of fulgides can be successfully tuned and improved by structural modifications at the chromophores.

The Influence of Spin Effects on the Gas Phase Reactions of Carbanions with N and O Atoms

Molecular anions have been recently detected in the denser regions of the interstellar medium. However, the chemical reactions of molecular anions with atomic species that are abundant in the ISM remain largely unexplored. This work is an experimental and computational study of CH(2)CN(-), CH(3)CHCN(-), (CH(3))(2)CCN(-), and CH(2)CHO(-) reacting with N and O atoms. In all cases the reactions of anions with O atoms exhibit larger reaction rate constants compared to the corresponding reactions with N atoms. Our study indicates that spin-forbidden reactions are the probable pathways in the reactions with N atoms, whereas spin-allowed reactions are the dominant processes in the reactions with O atoms. The major factor influencing the reaction rate constants of anions with N and O atoms is whether a spin-allowed barrierless pathway exists. The rich chemistry observed in this work provides a greater understanding of the ion-atom reaction processes, as well as some new avenues for further spin chemistry research.

Computational studies of gas phase reactions of carbon chain anions with N and O atoms

Experimental studies of gas phase reactions of carbanions with N and O atoms have been reported previously to understand ion chemistry relevant to the interstellar medium. In all cases reactions of anions with O atoms exhibit larger reaction rate constants compared to the corresponding N atom reactions. In addition, the open-shell carbon chain anions exhibit higher reactivities than the corresponding closed-shell species in N atom reactions, whereas similar reactivities were observed for both open and closed-shell anions in O atom reactions. These trends are investigated by the current theoretical study of the reactions of HC(n)(-)(n = 2, 4, and 6) and C(n)(-) (n = 2, 4-7) with N and O atoms. Our results indicate that spin-forbidden processes are the probable pathways in reactions of closed-shell anions HC(n)(-) with N atoms, and spin conversion limits the reaction efficiency. In reactions of open-shell anions C(n)(-) with N atoms, about 50% of the collisions may proceed through spin-allowed barrierless pathways, which results in relatively higher reaction efficiencies than for the closed-shell reactions. For reactions of all anions with O atoms, the spin-allowed barrierless pathways are the only channels, such that all reactions occur with very high efficiencies. This work provides a greater understanding of the influence of spin effects on the reactivities of anion reactions involving N and O atoms that may be important in the interstellar medium.

The photophysics of alloxazine: a quantum chemical investigation in vacuum and solution

(Time-dependent) Kohn-Sham density functional theory and a combined density functional/multi-reference configuration interaction method (DFT/MRCI) were employed to explore the ground and low-lying electronically excited states of alloxazine, a flavin related molecule. Spin-orbit coupling was taken into account using an efficient, nonempirical mean-field Hamiltonian. Intersystem crossing (ISC) rate constants for S --> T transitions were computed, employing both direct and vibronic spin-orbit coupling. Solvent effects were mimicked by a conductor-like screening model and micro-hydration with up to six explicit water molecules. Multiple minima were found on the first excited singlet (S(1)) potential energy hypersurface (PEH) with electronic structures (1)(npi*) and (1)(pipi*), corresponding to the dark 1 (1)A'' (S(1)) state and the nearly degenerate, optically bright 2 (1)A' (S(2)) state in the vertical absorption spectrum, respectively. In the vacuum the minimum of the (1)(npi*) electronic structure is clearly found below that of the (1)(pipi*) electronic structure. Population transfer from (1)(pipi*) to (1)(npi*) may proceed along an almost barrierless pathway. Hence, in the vacuum, internal conversion (IC) between the 2 (1)A' and the 1 (1)A'' state is expected to be ultrafast and fluorescence should be quenched completely. The depletion of the (1)(npi*) state is anticipated to occur via competing IC and direct ISC processes. In aqueous solution this changes, due to the blue shift of the (1)(npi*) state and the red shift of the (1)(pipi*) state. However, the minimum of the (1)(npi*) state still is expected to be found on the S(1) PEH. For vibrationally relaxed alloxazines pronounced fluorescence and ISC by a vibronic spin-orbit coupling mechanism is expected. At elevated temperatures or excess energy of the excitation laser, the (1)(npi*) state is anticipated to participate in the deactivation process and to partially quench the fluorescence.

Fully Reversible Transition from Wenzel to Cassie−Baxter States on Corrugated Superhydrophobic Surfaces

Liquid drops on textured surfaces show different dynamical behaviors depending on their wetting states. They are extremely mobile when they are supported by composite solid-liquid-air interfaces (Cassie-Baxter state) and immobile when they fully wet the textured surfaces (Wenzel state). By reversibly switching between these two states, it will be possible to achieve control over the fluid dynamics. Unfortunately, these wetting transitions are usually prevented by surface energy barriers. We demonstrate here a new, simple design paradigm consisting of parallel grooves with an appropriate aspect ratio that allows for the controlled, barrierless, reversible switching of the wetting states upon application of electrowetting. We report a direct observation of the barrierless dynamical pathway for the reversible transitions between the Wenzel (collapsed) and Cassie-Baxter (suspended) states and present a theory that accounts for these transitions, including detailed lattice Boltzmann simulations.

A Modified Four-State Model for the “Dual Fluorescence” of N6,N6-Dimethyladenine Derived from Femtosecond Fluorescence Spectroscopy

The radiationless deactivation of the excited electronic states of the dual fluorescence molecule N(6),N(6)-dimethyladenine (DMAde) was investigated using femtosecond time-resolved fluorescence up-conversion spectroscopy. The molecules were studied in solution in water and in dioxane. Fluorescence-time profiles were recorded in the wide wavelength range of 290 <or= lambda(fl) <or= 650 nm. The excitation wavelengths in the region of the first UV absorption band were tuned from close to the electronic origin (lambda(pump) = 294 nm) to excess energies of approximately 5400 cm(-1) above (lambda(pump) = 258 nm). Global fits to the measured curves turned out to reflect distinctive molecular relaxation processes on five well-defined time scales. Sub-100 fs and 0.52(3) ps lifetimes were found to predominate at the shortest UV and blue emission wavelengths in water, 1.5(1) and 3.0(2) ps components at intermediate wavelengths and a 62(1) ps value in the red region of the spectrum (2sigma error limits of the last digits in parentheses). In dioxane, these lifetimes changed to <or=0.27 and 0.63(4) ps in the UV, 1.5(1) and 10.9(10) ps in a wide range of intermediate, and 1.40(4) ns at the longest wavelengths. However, little dependence of the respective time constants on lambda(pump) was observed, indicating that the ensuing relaxation processes proceed via practically barrierless pathways through conical intersections. Building on the knowledge for the parent molecule adenine (Ade), the observations were rationalized with the help of a modified four-state model for the electronic dynamics in DMAde with the pipi*(L(a)), pipi*(L(b)), and npi* states similar to those in Ade and an intramolecular charge-transfer (ICT) state, which has no counterpart in Ade, responsible for the long-wavelength fluorescence.

Photochemical formation of HCO and CH3 on the ground S (A1′) state of CH3CHO

The dynamics of the photodissociation of CH(3)CHO into CH(3) + HCO products have been investigated at energies between 30,953 and 31,771 cm(-1), spanning the threshold for radical production on the triplet (T(1)) surface. A barrierless pathway to CH(3) + HCO radical products formed on the ground state (S(0)) surface was discovered and established to be an important reaction channel in acetaldehyde photodissociation throughout this wavelength range. HCO laser induced fluorescence (LIF) spectra recorded from CH(3)CHO dissociated above and below the T(1) barrier energy are quite different; HCO produced on S(0) yields a more congested LIF spectrum with sharp rotational transitions, while HCO formed on the T(1) surface displays fewer, more intense, Doppler-broadened lines. These differences have been further explored in the populations of the HCO K(a) = 1 doublets. Despite the upper and lower levels being almost isoenergetic, HCO formed on T(1) preferentially populates the upper K(c) state due to the geometry of the T(1) transition state structure. In contrast, HCO formed on S(0) produces equal population in each of the upper and lower K(a) = 1 components. Product state distributions (PSDs) showed that HCO formed on S(0) is born with an approximately statistical distribution of population in the available product states, modeled well by phase space theory. HCO formed on the T(1) surface, in contrast, has a PSD that can be characterized as arising from "impulsive" dynamics. Previous discrepancies in the height of the T(1) barrier are discussed following the observation that, once the T(1) channel is energetically accessible, there is competition between the S(0) and T(1) pathways, with the dominance of the triplet channel increasing with increasing photolysis energy.

Deactivation via ring opening: A quantum chemical study of the excited states of furan and comparison to thiophene

Minimum nuclear arrangements of the ground and low-lying excited electronic states of furan were obtained by means of (time dependent) Kohn–Sham density functional theory. A combined density functional/multi-reference configuration interaction method (DFT/MRCI) was employed to compute the spectral properties at these points. Multiple minima were found on the first excited singlet (S 1) potential energy hypersurface with electronic structures S1, S2, S3 corresponding to the 1 1A 2 (π → 3s-Ryd), 1 1B 2 (π → π ∗), and 2 1A 1 (π → π ∗) states in the vertical absorption spectrum, respectively. In analogy to recently published studies in thiophene [S. Salzmann, M. Kleinschmidt, J. Tatchen, R. Weinkauf, C.M. Marian, Phys. Chem. Chem. Phys. 10 (2008) 380] a deactivation mechanism for electronically excited furan was detected that involves the opening of the pentacyclic ring. We found a nearly barrierless relaxation pathway from the Franck–Condon region along a C–O bond-breaking coordinate. Hereby the initially excited 1B 2 (π → π ∗) state undergoes a conical intersection with a 1B 1 (π → σ ∗) state. The system can return to the electronic ground state through a second conical intersection of the 1(π → σ ∗) state before the minimum of that B 1 state is reached.

New Insights into the Reaction Mechanisms of Phenylium Ions with Benzene

The chemistry of phenylium (benzen-1-ylium) cations with benzene is investigated by using a guided ion beam tandem mass spectrometer. The main ionic products from the reaction of C6H5+ with C6H6 are observed at m/z 155 (covalent adduct C12H11+), 154 (C12H10+), 153 (C12H9+), 129 (C10H9+), and 115 (C9H7+). We propose a mechanism according to which channels at m/z 154-115 are formed by elimination of stable neutral molecules (such as H2, C2H2, C3H4) from the collision complex C12H11+, for which the most plausible structure is protonated biphenyl. The proposed mechanism is demonstrated by using partial isotopical labeling of reagents to look for possible H/D atom scrambling. Almost the same ions are produced when benzene is chemically ionized at atmospheric pressure in an APCI source from which oxygen is excluded. Because an ion trap analyzer is coupled to this source, tandem MS experiments can be performed, allowing structural details to be established. Moreover, the use of partially deuterated reagents has allowed the detection of minor reactive channels resulting from charge exchange and H-/D- hydride-transfer processes. Theoretical calculations show that the most stable structure for ions at m/z 129 C10H9+ is that of protonated naphthalene, resulting from the loss of an acetylene molecule by the condensation product, with a reaction exothermicity of 1.27 eV. We have found a possible barrierless pathway for such a channel that might be viable for the synthesis of naphthalene, the smallest PAH, even at low collision energies and therefore would be of particular astrochemical relevance.

Novel mechanism for oxidation of CO by Fe 2O 3 clusters

The oxidation of CO by Fe 2O 3 clusters has been studied using a first principles gradient corrected density functional approach. It is shown that the direct oxidation of CO via surface oxygen atoms is feasible with a low barrier of 0.4 eV. The present studies also suggest a new indirect mechanism for CO oxidation. Here, the first CO molecule adsorbed on a surface Fe site breaks one of the FeO bonds and a subsequent CO molecule is oxidized by the less coordinated O atom in an almost barrierless reaction pathway. Subsequent attachment of CO and O 2 to Fe sites results in the formation of a CO 3 complex where the heat of formation is enough to allow CO 2 to leave the surface reverting the system back to Fe 2O 3. The role of charge transfer and the nature of the transition states are highlighted. It is shown that small clusters are ferromagnetic indicating that super-exchange interactions of the antiferromagnetic bulk are modified at the reduced sizes.