Examination of the Silver Colloid Binding Behavior of Disulfide-Tethered Bipyridine Ligands and Their fac-Tricarbonylrhenium(I) Complexes

The syntheses of 2,2'-bipyridin-5-ylmethyl-5-(1,2-dithiolan-3-yl)pentanoate (L1) and N-(2,2'-bipyridin-5-ylmethyl)-5-(1,2-dithiolan-3-yl)pentanamide (L2) and their neutral fac carbonylrhenium(I) complexes [Re(L1)(CO)(3)Br] and [Re(L2)(CO)(3)Br] are reported. The. electronic absorption and emission spectra of the complexes are similar to the spectrum of the reference compound [Re(bipy)(CO)(3)Br] and correlate well with the density functional theory calculations undertaken. The surface-enhanced Raman spectroscopy (SERS) spectra (excited at both 532 and 785 nm) of the ligands and complexes were examined and compared to the spectrum of ethyl 5-(1,2-dithiolan-3-yl)pentanoate (L3), revealing that there is very little contribution to the spectra of these species from the dithiolated alkyl chains. The spectra are dominated by the characteristic peaks of a metalated 2,2'-bipyridyl group,arising from the silver colloid/ion complexation, and the rhenium center. The rhenium complexes show weak SERS bands related to the CO stretches and a broad band at 510 cm(-1) assigned to Re-CO stretching. Concentration dependent studies, measured by the relative intensity of several assigned peaks, indicate that, as the surface coverage increases, the bipyridine moiety lifts off the surface In the case of L1 and L2, this gives rise to complexes with silver at low concentration, enhancing the signals observed, while for the tricarbonylbromorhenium complexes of these ligands, the presence of the disulfide tether allows an enhancement in the limits of detection of these surface-borne species of 20 times in the case of [ReL2(CO)(3)Br] over [Re(bipy)(CO)(3)Br].

Ab initio potentials of F + Li-2 accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li-2 + F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li-2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Effect of positron-atom interactions on the annihilation gamma spectra of molecules

Calculations of ?-spectra for positron annihilation on a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The annihilation ?-spectra characterise the momentum distribution of the electron-positron pair at the instant of annihilation. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum densities calculated using modern computational quantum chemistry density functional theory tools. The calculation, in its simplest form, effectively treats the low-energy (thermalised, room-temperature) positron as a plane wave and gives annihilation ?-spectra that are about 40% broader than experiment, although the main chemical trends are reproduced. We show that this effective 'narrowing' of the experimental spectra is due to the action of the molecular potential on the positron, chiefly, due to the positron repulsion from the nuclei. It leads to a suppression of the contribution of small positron-nuclear separations where the electron momentum is large. To investigate the effect of the nuclear repulsion, as well as that of short-range electron-positron and positron-molecule correlations, a linear combination of atomic orbital description of the molecular orbitals is employed. It facilitates the incorporation of correction factors which can be calculated from atomic many-body theory and account for the repulsion and correlations. Their inclusion in the calculation gives -spectrum linewidths that are in much better agreement with experiment. Furthermore, it is shown that the effective distortion of the electron momentum density, when it is observed through positron annihilation -spectra, can be approximated by a relatively simple scaling factor. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.

Calculation of gamma spectra for positron annihilation on molecules

Calculations of gamma spectra for positron annihilation for a selection of molecules, including methane and its fluoro-substitutes, ethane, propane, butane and benzene are presented. The contribution to the ?-spectra from individual molecular orbitals is obtained from electron momentum distributions calculated using the density functional theory (DFT) based B3LYP/TZVP model. For positrons thermalised to room temperature, the calculation, in its simplest form, effectively treats the positron as a plane wave and gives positron annihilation ?-spectra linewidths that are broader (30-40%) than experiment, although the main chemical trends are reproduced. The main physical reason for this is the neglect of positron repulsion from the nuclei. We show that this effect can be incorporated through momentum-dependent correction factors, determined from positron-atom calculations, e.g., many-body perturbation theory. Inclusion of these factors in the calculation gives linewidths that are in improved agreement with experiment.

On the Mechanism of Lewis Acid Catalyzed Glucose Transformations in Ionic Liquids

A complementary computational and experimental study of the reactivity of Lewis acidic CrCl2, CuCl2 and FeCl2 catalysts towards glucose activation in dialkylimidazolium chloride ionic liquids is performed. The selective dehydration of glucose to 5-hydroxymethylfurfural (HMF) proceeds through the intermediate formation of fructose. Although chromium(II) and copper(II) chlorides are able to dehydrate fructose with high HMF selectivity, reasonable HMF yields from glucose are only obtained with CrCl2 as the catalyst. Glucose conversion by CuCl2 is not selective, while FeCl2 catalyst does not activate sugar molecules. These differences in reactivity are rationalized on the basis of in situ X-ray absorption spectroscopy measurements and the results of density functional theory calculations. The reactivity in glucose dehydration and HMF selectivity are determined by the behavior of the ionic liquid-mediated Lewis acid catalysts towards the initial activation of the sugar molecules. The formation of a coordination complex between the Lewis acidic Cr2+ center and glucose directs glucose transformation into fructose. For Cu2+ the direct coordination of sugar to the copper(II) chloride complex is unfavorable. Glucose deprotonation by a mobile Cl- ligand in the CuCl42- complex initiates the nonselective conversion. In the course of the reaction the Cu2+ ions are reduced to Cu+. Both paths are prohibited for the FeCl2 catalyst.

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]