Coordination Properties of Ionic Liquid-Mediated Chromium(II) and Copper(II) Chlorides and Their Complexes with Glucose

The structural and coordination properties of complexes formed upon the interaction of copper(II) and chromium(II) chlorides with diallrylimidazolium chloride (RMlm(+)Cl(-)) ionic liquids and glucose are studied by a combination of density functional theory (DFT) calculations and X-ray absorption spectroscopy (XAS). In the absence of the carbohydrate substrate, isolated mononuclear four-coordinated MeCl42- species (Me = Cu, Cr) dominate in the ionic liquid solution. The organic part of the ionic liquid does not directly interact with the metal centers. The interactions between the RMlm(+) cations and the anionic metal chloride complexes are limited to hydrogen bonding with the basic Cl- ligands and the overall electrostatic stabilization of the anionic metal complexes. Exchange of Cl ligands by a hydroxyl group of glucose is only favorable for CrCl42-. For Cu2+ complexes, the formation of hydrogen bonded complexes between CuCl42- and glucose is preferred. No preference for the coordination of metal chloride species to specific hydroxyl group of the carbohydrate is found. The formation of binuclear metal chloride complexes is also considered. The reactivity and selectivity patterns of the Lewis acid catalyzed reactions of glucose are discussed in the framework of the obtained results.

Origin of the Increase of Activity and Selectivity of Nickel Doped by Au, Ag, and Cu for Acetylene Hydrogenation

Activity and selectivity are both important issues in heterogeneous catalysis and recent experimental results have shown that Ni catalysts doped by gold exhibit high activity for the hydrogenation of acetylene with good selectivity of ethylene formation. To unravel the underlying mechanism for this observation, the general trend of activity and selectivity of Ni surfaces doped by Au, Ag, and Cu has been investigated using density functional theory calculations. Complete energy profiles from C2H2 to C2H4 on Ni(111), Au/Ni(111), Ag/Ni(111) and Cu/Ni(111) are obtained and their turnover frequencies (TOFs) are computed. The results show that acetylene adsorption on Ni catalyst is strong which leads to the low activity while the doping of Au, Ag, and Cu on the Ni catalyst weakens the acetylene adsorption, giving rise to the increase of activity. The selectivity of ethylene formation is also quantified by using the energy difference between the hydrogenation barriers and the absolute value of the adsorption energies of ethylene. It is found that the selectivity of ethylene formation increases by doping Au and Ag, while those of Cu/Ni and Ni are similar.

Electronic stopping power of H and He in Al and LiF from first principles

Non-linearities in the electronic stopping power of light projectiles in bulk Al and LiF are addressed from first principles using time-evolving time-dependent density functional theory. In the case of Al, the agreement of the calculations with experiments for H and He projectiles is fair, but a recently observed transition for He from one value of the electronic friction coefficient to a higher value at v ~ 0.3 a.u. is not reproduced by the calculations. For LiF, better accuracy is obtained as compared with previously published simulations, albeit the threshold remains overestimated. © 2013 Elsevier B.V.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 American Institute of Physics.

Magnetic properties of copper hexadecaphthalocyanine (F16CuPc) thin films and powders

The structural and magnetic properties of F16CuPc thin films and powder, including x-ray diffraction (XRD), superconducting quantum interference device (SQUID) magnetometry, and theoretical modelling of exchange interactions are reported. Analysis of XRD from films, with thickness ranging between 100 and 160 nm, deposited onto Kapton and a perylene-3,4,9,10-tetracarboxylic-3,4,9,10-dianhydride (PTCDA) interlayer shows that the stacking angle (defined in the text) of the film is independent of the thickness, but that the texture is modified by both film thickness and substrate chemistry. The SQUID measurements suggest that all samples are paramagnetic, a result that is confirmed by our theoretical modelling including density functional theory calculations of one-dimensional molecular chains and Green's function perturbation theory calculations for a molecular dimer. By investigating theoretically a range of different geometries, we predict that the maximum possible exchange interaction between F16CuPc molecules is twice as large as that in unfluorinated copper-phthalocyanine (CuPc). This difference arises from the smaller intermolecular spacing in F16CuPc. Our density functional theory calculation for isolated F16CuPc molecule also shows that the energy levels of Kohn-Sham orbitals are rigidly shifted similar to 1 eV lower in F16CuPc compared to CuPc without a significant modification of the intramolecular spin physics, and that therefore the two molecules provide a suitable platform for independently varying magnetism and charge transport. 

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. 

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.

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.