Combined studies of DFT atomistic modelling and in situ FTIR spectroscopy on surface oxidants and CO oxidation at Ru electrodes

We report the combined studies of density functional theory (DFT) calculations and electrochemical in situ FTIR spectroscopy on surface oxidants and mechanisms of CO oxidation at the Ru(0001) electrodes. It is shown that CO can co-adsorb with both O and OH species at lower potential region where a low coverage of the (2 x 2)-O/OH adlayer formed; the oxidation of CO adsorbates takes place at higher potentials where a high coverage of the (1 x 1)-O/OH adlayer formed. Surface O species are not the active oxidants under all coverages studied, due to the high reaction barriers between CO and O (>1 eV). However, surface OH species with higher coverage are identified as the active oxidants, and CO oxidation takes place via a two-steps' mechanism of CO + 3OH -> COOH + 2OH -> CO2 + H2O + OH, in which three nearby OH species are involved in the CO2 formation: CO reacts with OH, forming COOH; COOH then transfers the H to a nearby OH to form H2O and CO2, at the same time, another H in the H2O transfers to a nearby OH to form a weak adsorbed H2O and a new OH. The reaction barrier of these processes is reduced significantly to around 0.50 eV. These new results not only provide an insight into surface active oxidants on Ru, which is directly relevant to fuel cell catalysis, but also reveals the extra complexity of catalytic reactions taking place at solid/liquid electrochemical interface in comparison to the relatively simpler ones at solid/gas phase. 

Chiral thiouronium salts: synthesis, characterisation and application in NMR enantio-discrimination of chiral oxoanions

Chiral thioureas and functionalised chiral thiouronium salts were synthesised starting from the relatively cheap and easily available chiral amines: (S)-methylbenzylamine and rosin-derived (+)-dehydroabietylamine. The introduction of a delocalised positive charge to the thiourea functionality, by an alkylation reaction at the sulfur atom, enables dynamic rotameric processes: hindered rotations about the delocalised CN and CS bonds. Hence, four different rotamers/isomers may be recognised: syn-syn, syn-anti, anti-syn and anti-anti. Extensive H-1 and C-13 NMR studies have shown that in hydrogen-bond acceptor solvents, such as perdeuteriated dimethyl sulfoxide, the syn-syn conformation is preferable. On the other hand, when using non-polar solvents, such as CDCl3, the mixture of syn-syn and syn-anti isomers is detectable, with an excess of the latter. Apart from this, in the case of S-butyl-N,N'-bis(dehydroabietyl)thiouronium ethanoate in CDCl3, the H-1 NMR spectrum revealed that strong bifurcated hydrogen bonding between the anion and the cation causes global rigidity without signs of hindered rotamerism observable on the NMR time scale. This suggested that these new salts might be used as NMR discriminating agents for chiral oxoanions, and are indeed more effective than their archetypal guanidinium analogues or the neutral thioureas. The best results in recognition of a model substrate, mandelate, were obtained with S-butyl-N,N'-bis(dehydroabietyl) thiouronium bistriflamide. It was confirmed that the chiral recognition occurred not only for carboxylates but also for sulfonates and phosphonates. Further H-1 NMR studies confirmed a 1 : 1 recognition mode between the chiral agent (host) and the substrate (guest); binding constants were determined by H-1 NMR titrations in solutions of DMSO-d(6) in CDCl3. It was also found that the anion of the thiouronium salt had a significant influence on the recognition process: anions with poor hydrogen-bond acceptor abilities led to the best discrimination. The presence of host-guest hydrogen bonding was confirmed in the X-ray crystal structure of S-butyl-N,N'-bis(dehydroabietyl)thiouronium bromide and by computational studies (density functional theory).

Surface energy and surface proton order of the ice Ih basal and prism surfaces

Density-functional theory (DFT) is used to examine the basal and prism surfaces of ice Ih. Similar surface energies are obtained for the two surfaces; however, in each case a strong dependence of the surface energy on surface proton order is identified. This dependence, which can be as much as 50% of the absolute surface energy, is significantly larger than the bulk dependence (< 1%) on proton order, suggesting that the thermodynamic ground state of the ice surface will remain proton ordered well above the bulk order-disorder temperature of about 72 K. On the basal surface this suggestion is supported by Monte Carlo simulations with an empirical potential and solution of a 2D Ising model with nearest neighbor interactions taken from DFT. Order parameters that define the surface energy of each surface in terms of nearest neighbor interactions between dangling OH bonds (those which point out of the surface into vacuum) have been identified and are discussed. Overall, these results suggest that proton order-disorder effects have a profound impact on the stability of ice surfaces and will most likely have an effect on ice surface reactivity as well as ice crystal growth and morphology. S Supplementary data are available from stacks.iop.org/JPhysCM/22/074209/mmedia

The phase diagram of water at negative pressures: Virtual ices

The phase diagram of water at negative pressures as obtained from computer simulations for two models of water, TIP4P/2005 and TIP5P is presented. Several solid structures with lower densities than ice Ih, so-called virtual ices, were considered as possible candidates to occupy the negative pressure region of the phase diagram of water. In particular the empty hydrate structures sI, sII, and sH and another, recently proposed, low-density ice structure. The relative stabilities of these structures at 0 K was determined using empirical water potentials and density functional theory calculations. By performing free energy calculations and Gibbs-Duhem integration the phase diagram of TIP4P/2005 was determined at negative pressures. The empty hydrates sII and sH appear to be the stable solid phases of water at negative pressures. The phase boundary between ice Ih and sII clathrate occurs at moderate negative pressures, while at large negative pressures sH becomes the most stable phase. This behavior is in reasonable agreement with what is observed in density functional theory calculations.

Palladium-catalyzed liquid-phase hydrogenation/hydrogenolysis of disulfides

For the first time, the hydrogenation/hydrogenolysis of a range of disulfides has been achieved over a supported palladium catalyst using hydrogen under relatively benign conditions. These unexpected results demonstrate that it is possible to avoid the poisoning of the catalyst by either the nitrogen-containing groups or the sulfur species, allowing both efficient reaction and recycling of the catalyst under the proper conditions (e.g., at low temperatures). A slight loss in activity was found on recycling; however, the catalyst activity can be recovered using hydrogen pretreatment. The reaction mechanism for the hydrogenolysis and hydrogenation of ortho-, meta-, and para-dinitrodiphenyldisulfide to the corresponding aminothiophenol has been elucidated. Density functional theory calculations were used to investigate the adsorption mode of the dinitrodiphenyldisulfides; a clear dependence on adsorption geometry was found regarding whether the molecule is cleaved at the S-S bond before the reduction of the nitro group or vice versa. This study demonstrates the versatility of these catalysts for the hydrogenation/hydrogenolysis of sulfur-containing molecules, which normally are considered poisons, and will extend their use to a new family of substrates. (C) 2007 Elsevier Inc. All rights reserved.

The Role of Arg228 in the phosphorylation of galactokinase: The mechanism of GHMP kinases by QM/MM studies

GHMP kinases are a group of structurally-related small molecule kinases. They have been found in all kingdoms of life and are mostly responsible for catalysing the ATP-dependent phosphorylation of intermediary metabolites. Although the GHMP kinases are of clinical, pharmaceutical and biotechnological importance, the mechanism of GHMP-kinases is controversial. A catalytic base mechanism was suggested for mevalonate kinase that has a structural feature of the ?-phosphate of ATP close to an aspartate residue; however, for one GHMP member, homoserine kinase, where the residue acting as general base is absent, a direct phosphorylation mechanism was suggested. Furthermore, it has been proposed by some authors that all the GHMP kinases function via the direct phosphorylation mechanism. This controversy in mechanism has limited our ability to exploit these enzymes as drug targets and in biotechnology. Here the phosphorylation reaction mechanism of the human galactokinase, a member of GHMP kinase was investigated using molecular dynamics simulations and density functional theory-based QM/MM calculations (B3LYP-D/AMBER99). The reaction coordinates were localized by potential energy scan using adiabatic mapping method. Our results indicate that a highly conserved Glu174 captures Arg105 to the proximity of the a-phosphate of ATP forming a H-bond network, therefore the mobility of ATP in the large oxyanion hole is restricted. Arg228 functions to stabilize the negative charge developed at the ß,?-bridging oxygen of the ATP during bond cleavage. The reaction occurs via direct phosphorylation mechanism and the Asp186 in proximity of ATP does not directly participate in the reaction pathway. Since Arg228 is not conserved among GHMP kinases, reagents which form interactions with Arg228, and therefore can interrupt its function in phosphorylation may be developed into potential selective inhibitors for galactokinase.

A multiconfigurational time-dependent Hartree-Fock method for excited electronic states. I. General formalism and application to open-shell states

The solution of the time-dependent Schrodinger equation for systems of interacting electrons is generally a prohibitive task, for which approximate methods are necessary. Popular approaches, such as the time-dependent Hartree-Fock (TDHF) approximation and time-dependent density functional theory (TDDFT), are essentially single-configurational schemes. TDHF is by construction incapable of fully accounting for the excited character of the electronic states involved in many physical processes of interest; TDDFT, although exact in principle, is limited by the currently available exchange-correlation functionals. On the other hand, multiconfigurational methods, such as the multiconfigurational time-dependent Hartree-Fock (MCTDHF) approach, provide an accurate description of the excited states and can be systematically improved. However, the computational cost becomes prohibitive as the number of degrees of freedom increases, and thus, at present, the MCTDHF method is only practical for few-electron systems. In this work, we propose an alternative approach which effectively establishes a compromise between efficiency and accuracy, by retaining the smallest possible number of configurations that catches the essential features of the electronic wavefunction. Based on a time-dependent variational principle, we derive the MCTDHF working equation for a multiconfigurational expansion with fixed coefficients and specialise to the case of general open-shell states, which are relevant for many physical processes of interest. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3600397]

Performance of Nonlocal Optics When Applied to Plasmonic Nanostructures

Semiclassical nonlocal optics based on the hydrodynamic description of conduction electrons might be an adequate tool to study complex phenomena in the emerging field of nanoplasmonics. With the aim of confirming this idea, we obtain the local and nonlocal optical absorption spectra in a model nanoplasmonic device in which there are spatial gaps between the components at nanometric and subnanometric scales. After a comparison against time-dependent density functional calculations, we conclude that hydrodynamic nonlocal optics provides absorption spectra exhibiting qualitative agreement but not quantitative accuracy. This lack of accuracy, which is manifest even in the limit where induced electric currents are not established between the constituents of the device, is mainly due to the poor description of induced electron densities.

Electrooxidation of methanol in an alkaline fuel cell: determination of the nature of the initial adsorbate

It is essential to correctly determine the nature of the initial adsorbate in order to calculate the pathway for any given reaction. Recent literature provides conflicting information on the first step in the methanol decomposition pathway. This work sets out to establish what role the solution and the surface have to play in the initial adsorption-deprotonation process. Density functional theory (DFT) calculations, in combination with a cluster-continuum model approach are used to resolve the nature of the adsorbing species. We show that methanol is the dominant species in solution over methoxide, and also has a smaller barrier to adsorption. The nature of the surface species is revealed to be a methanol-OH complex.

The origin of high activity but low CO2 selectivity on binary PtSn in the direct ethanol fuel cell

The most active binary PtSn catalyst for direct ethanol fuel cell applications has been studied at 20 ^oC and 60 ^oC, using variable temperature electrochemical in-situ FTIR. In comparison with Pt, binary PtSn inhibits ethanol dissociation to CO(a), but promotes partial oxidation to acetaldehyde and acetic acid. Increasing the temperature from 20 ^oC to 60 ^oC facilitates both ethanol dissociation to CO(a) and their further oxidation to CO₂, leading to an increased selectivity towards CO₂; however, acetaldehyde and acetic acid are still the main products. Potential-dependent phase diagrams for surface oxidants of OH(a) formation on Pt(111), Pt(211) and Sn modified Pt(111) and Pt(211) surfaces have been determined using density functional theory (DFT) calculations. It is shown that Sn promotes the formation of OH(a) with a lower onset potential on the Pt(111) surface, whereas an increase in the onset potential is found on modification of the (211) surface. In addition, Sn inhibits the Pt(211) step edge with respect to ethanol C-C bond breaking compared with that found on the pure Pt, which reduces the formation of CO(a). Sn was also found to facilitate ethanol dehydrogenation and partial oxidation to acetaldehyde and acetic acid which, combined with the more facile OH(a) formation on the Pt(111) surface, gives us a clear understanding of the experimentally determined results. This combined electrochemical in-situ FTIR and DFT study, provides, for the first time, an insight into the long-term puzzling features of the high activity but low CO₂ production found on binary PtSn ethanol fuel cell catalysts.

Significance of β-dehydrogenation in ethanol electro-oxidation on platinum doped with Ru, Rh, Pd, Os and Ir

In the exploration of highly efficient direct ethanol fuel cells (DEFCs), how to promote the CO₂ selectivity is a key issue which remains to be solved. Some advances have been made, for example, using bimetallic electrocatalysts, Rh has been found to be an efficient additive to platinum to obtain high CO₂ selectivity experimentally. In this work, the mechanism of ethanol electrooxidation is investigated using first principles method. It is found that CH₃CHOH* is the key intermediate during ethanol electrooxidation and the activity of β-dehydrogenation is the rate determining factor that affects the completeness of ethanol oxidation. In addition, a series of transition metals (Ru, Rh, Pd, Os and Ir) are alloyed on the top layer of Pt(111) in order to analyze their effects. The elementary steps, α-, β-C-H bond and C-C bond dissociations are calculated on these bimetallic M/Pt(111) surfaces and the formation potential of OH* from water dissociation is also calculated. We find that the active metals increase the activity of β-dehydrogenation but lower the OH* formation potential resulting in the active site being blocked. By considering both β-dehydrogenation and OH* formation, Ru, Os and Ir are identified to be unsuitable for the promotion of CO₂ selectivity and only Rh is able to increase the selectivity of CO₂ in DEFCs.

The effects of stepped sites and ruthenium adatom decoration on methanol dehydrogenation over platinum-based catalyst surfaces

A density functional theory study of methanol dehydrogenation over stepped Pt(2 1 1) surfaces without and with Ru modification was carried out to understand fuel catalytic reactions on Pt-based catalysts. Two main pathways of the CH3OH dehydrogenation were examined: the O–H pathway which was initiated by O–H bond scission to form the methoxy (CH3O) intermediate followed by sequential cleavage of C–H bonds to CO, and the C–H pathway which was initiated by C–H bond scission to form the hydroxymethyl (CH2OH) followed by two C–H bond cleavages to COH and then CO. Possible crossover reactions between the O–H and C–H pathways were also computed. Compared to flat Pt(1 1 1), stepped Pt(2 1 1) increases the adsorption energies of intermediates, making no significant contribution to decreasing the reaction barriers of most elementary steps involved, except in the first hydrogen scission. However, on the Ru-modified surface, a significant reduction was found in reaction barriers for the first step of the C–H bond scission and a number of further dehydrogenation steps crossing over to the O–H pathway, with the most facile paths identified. Our data reveals the complexity of methanol catalytic reaction processes at the atomic level and contributes to a fundamental understanding of fuel reactions on Pt-based catalysts.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Resonantly Enhanced Multi-Photon Ionization Spectrum of the Neutral Green Fluorescent Protein Chromophore

The photophysics of the green fluorescent protein is governed by the electronic structure of the chromophore at the heart of its β-barrel protein structure. We present the first two-color, resonance-enhanced, multiphoton ionization spectrum of the isolated neutral chromophore in vacuo with supporting electronic structure calculations. We find the absorption maximum to be 3.65 ± 0.05 eV (340 ± 5 nm), which is blue-shifted by 0.5 eV (55 nm) from the absorption maximum of the protein in its neutral form. Our results show that interactions between the chromophore and the protein have a significant influence on the electronic structure of the neutral chromophore during photoabsorption and provide a benchmark for the rational design of novel chromophores as fluorescent markers or photomanipulators.