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]

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.

Modeling study of woody biomass:Interactions of cellulose, hemicellulose, and lignin

Lignocellulosic biomass pretreatment and the subsequent thermal conversion processes to produce solid, liquid, and gas biofuels are attractive solutions for today's energy challenges. The structural study of the main components in biomass and their macromolecular complexes is an active and ongoing research topic worldwide. The interactions among the three main components, cellulose, hemicellulose, and lignin, are studied in this paper using electronic structure methods, and the study includes examining the hydrogen bond network of cellulose-hemicellulose systems and the covalent bond linkages of hemicellulose-lignin systems. Several methods (semiempirical, Hartree-Fock, and density functional theory) using different basis sets were evaluated. It was shown that theoretical calculations can be used to simulate small model structures representing wood components. By comparing calculation results with experimental data, it was concluded that B3LYP/6-31G is the most suitable basis set to describe the hydrogen bond system and B3LYP/6-31G(d,p) is the most suitable basis set to describe the covalent system of woody biomass. The choice of unit model has a much larger effect on hydrogen bonding within cellulose-hemicellulose system, whereas the model choice has a minimal effect on the covalent linkage in the hemicellulose-lignin system. © 2011 American Chemical Society.

Formation mechanism of levoglucosan and formaldehyde during cellulose pyrolysis

Biomass pyrolysis is an efficient way to transform raw biomass or organic waste materials into useable energy, including liquid, solid, and gaseous materials. Levoglucosan (1,6-anhydro-β-d-glucopyranose) and formaldehyde are two important products in biomass pyrolysis. The formation mechanism of these two products was investigated using the density functional theory (DFT) method based on quantum mechanics. It was found that active anhydroglucose can be obtained from a cellulose homolytic reaction during high-temperature steam gasification of the biomass process. Anhydroglucose undergoes a hydrogen-donor reaction and forms an intermediate, which can transform into the products via three pathways, one (path 1) for the formation of levoglucosan and two (paths 2 and 3) for formaldehyde. A total of six elementary reactions are involved. At a pressure of 1 atm, levoglucosan can be formed at all of the temperatures (450-750 K) considered in this simulation, whereas formaldehyde can be formed only when the temperature is higher than 475 K. Moreover, the energy barrier of levoglucosan formation is lower than that of formaldehyde, which is in agreement with the mechanism proposed in the experiments. © 2011 American Chemical Society.

Thermal decomposition mechanism of levoglucosan during cellulose pyrolysis

Levoglucosan (1,6-anhydro-β-d-glucopyranose) decomposition is an important step during cellulose pyrolysis and for secondary tar reactions. The mechanism of levoglucosan thermal decomposition was studied in this paper using density functional theory methods. The decomposition included direct CO bond breaking, direct CC bond breaking, and dehydration. In total, 9 different pathways, including 16 elementary reactions, were studied, in which levoglucosan serves as a reactant. The properties of the reactants, transition states, intermediates, and products for every elementary reaction were obtained. It was found that 1-pentene-3,4-dione, acetaldehyde, 2,3-dihydroxypropanal, and propanedialdehyde can be formed from the CO bond breaking decomposition reactions. 1,2-Dihydroxyethene and hydroxyacetic acid vinyl ester can be formed from the CC bond breaking decomposition reactions. It was concluded that CO bond breaking is easier than CC bond breaking due to a lower activation energy and a higher released energy. During the 6 levoglucosan dehydration pathways, one water molecule which composed of a hydrogen atom from C3 and a hydroxyl group from C2 is the preferred pathway due to a lower activation energy and higher product stability. © 2012 Elsevier B.V. All rights reserved.