PubMed Health⌕ Search

Biomedical subjects

Notker Rösch

Publications and source records attributed to Notker Rösch.

At least 19 recordsLinked to original sources

Density functional embedded cluster study of Cu(4), Ag(4) and Au(4) species interacting with oxygen vacancies on the MgO(001) surface.

Cu(4), Ag(4), and Au(4) species adsorbed on an MgO(001) surface that exhibits neutral (F(s)) and charged (F(s) (+)) oxygen vacancies have been studied using a density functional approach and advanced embedding models. The gas-phase rhombic-planar structure of the coinage metal tetramers is only moderately affected by adsorption. In the most stable surface configuration, the plane of the tetramers is oriented perpendicular to the MgO(001) surface; one metal atom is attached to an oxygen vacancy and another one is bound to a nearby surface oxygen anion. A very similar structural motif was recently found on defect-free MgO(001), where two O(2-) ions serve as adsorption sites. Following the trend of the interactions with the regular MgO(001) surface, Au(4) and Cu(4) bind substantially stronger to F(s) and F(s) (+) sites than Ag(4). This stronger adsorption interaction at oxygen vacancies, in particular at F(s), is partly due to a notable accumulation of electron density on the adsorbates. We also examined the propensity of small supported metal species to aggregate to adsorbed di-, tri- and tetramers. Furthermore, we demonstrated that core-level ionization potentials offer the possibility for detecting experimentally supported metal tetramers and characterizing them structurally with the help of calculated data.

Journal Article↗

Surface composition of materials used as catalysts for methanol steam reforming: a theoretical study.

PdZn (1:1) alloy is assumed to be the active component of a promising catalyst for methanol steam reforming. Using density functional calculations on periodic supercell slab models, followed by atomistic thermodynamics modeling, we study the chemical composition of the surfaces PdZn(111) and, as a reference, Cu(111) in contact with water and hydrogen at conditions relevant to methanol steam reforming. For the two surfaces, we determine similar maximum adsorption energies for the dissociative adsorption of H(2), O(2), and the molecular adsorption of H(2)O. These reactions are calculated to be exothermic by about -40, -320, and -20 kJ mol(-1), respectively. Using a thermodynamic analysis based on theoretically predicted adsorption energies and vibrational frequencies, we determine the most favorable surface compositions for given pressure windows. However, surface energy plots alone cannot provide quantitative information on individual coverages in a system of coupled adsorption reactions. To overcome this limitation, we employ a kinetic model, from which equilibrium surface coverages of H, O, OH, and H(2)O are derived. We also discuss the sensitivity of our results and the ensuing conclusions with regard to the model surfaces employed and the inaccuracies of our computational method. Our kinetic model predicts surfaces of both materials, PdZn and Cu, to be essentially adsorbate-free already from very low values of the partial pressure of H(2). The model surfaces PdZn(111) and Cu(111) are predicted to be free of water-related adsorbates for a partial H(2) pressure greater than 10(-8) and 10(-5) atm, respectively.

Journal Article↗

Comparative theoretical study of formaldehyde decomposition on PdZn, Cu, and Pd surfaces.

Methanol steam reforming, catalyzed by Pd/ZnO (PdZn alloy), is a potential source of hydrogen for on-board fuel cells. CO has been reported to be a minor side product of methanol decomposition that occurs in parallel to methanol steam reforming on PdZn catalysts. However, fuel cells currently used in vehicles are very sensitive to CO poisoning. To contribute to the understanding of pertinent reaction mechanisms, we employed density functional slab model calculations to study the decomposition of formaldehyde, a key intermediate in methanol decomposition and steam reforming reactions, on planar surfaces of Pd, Cu, and PdZn as well as on a stepped surface of PdZn. The calculated activation energies indicate that dehydrogenation of formaldehyde is favorable on Pd(111), but unfavorable on Cu(111) and PdZn(111). On the stepped PdZn(221) surface, the dehydrogenation process was calculated to be more competitive to formaldehyde desorption than on PdZn(111). Thus, we ascribe the experimentally observed small amount of CO, formed during steam reforming of methanol on the Pd/ZnO catalyst, to occur at metallic Pd species of the catalyst or at defect sites of PdZn alloy.

Journal Article↗

The heat of formation of gaseous PuO(2)2+ from relativistic density functional calculations.

Using a set of model reactions, we estimated the heat of formation of gaseous PuO2(2+) from quantum-chemical reaction enthalpies and experimental heats of formation of reference species. To this end, we carried out relativistic density functional calculations on the molecules PuO(2)2+, PuO2, PuF6, and PuF4. We used a revised variant (PBEN) of the Perdew-Burke-Ernzerhof gradient-corrected exchange-correlation functional, and we accounted for spin-orbit interaction in a self-consistent fashion. As open-shell Pu species with two or more unpaired 5f electrons are involved, spin-orbit interaction significantly affects the energies of the model reactions. Our theoretical estimate for the heat of formation DeltafH degree 0(PuO2(2+),g), 418+/-15 kcal mol-1, evaluated using plutonium fluorides as references, is in good agreement with a recent experimental result, 413+/-16 kcal mol-1. The theoretical value connected to the experimental heat of formation of PuO2(g) has a notably higher uncertainty and therefore was not included in the final result.

Journal Article↗

Molecular-dynamics simulations of pyronine 6G and rhodamine 6G dimers in aqueous solution.

We have carried out molecular-dynamics (MD) simulations on dimers of the positively charged laser dyes pyronine 6G (P6G) and rhodamine 6G (R6G) in aqueous solution, generating trajectories of 2.5 ns for various computational protocols. We discuss how the choice of atomic partial charges and the length of the trajectories affect the predicted structures of the dimers and compare our results to those of earlier MD-simulations, which were restricted to only 0.7 ns. Our results confirm that monomers of P6G easily undergo relative rotations within the dimer, but we found new conformations of the R6G dimer at longer simulation times. In addition, we analyzed in detail the energy change during the formation of dimers. With suitable corrections, the electrostatic energy from an Ewald treatment agrees with the results from an approach relying on a residue-based cutoff. For P6G, we show that the strong solvent-mediated electrostatic attraction between the monomers is counteracted by an almost equally large solvent-induced entropy contribution to yield a small driving force to dimer formation, in very good agreement with the free-energy change from a thermodynamic-integration procedure. Thus, earlier rationalizations of the dimer formation, based only on energy arguments, yield a qualitatively wrong picture.

Computer Simulation↗

Two hydrogen ligands on tetrairidium clusters: a relativistic density functional study.

Structural and energetic properties of Ir(4)H(2) have been determined by applying a relativistic density functional method. As previously obtained for Ir(4)H, terminal coordination of H ligands is preferred, in contrast to some other transition metals. Square-planar Ir(4) isomers with an H binding energy of up to 318 kJ mol(-1) per atom were determined as the most stable structures of Ir(4)H(2). Isomers with a tetrahedral or butterfly structure of the metal framework exhibit average H atom binding energies of up to approximately 300 kJ mol(-1). For all three types of isomers, a surprisingly large number of stable minima was identified. Unexpectedly, structural as well as energetic properties of Ir(4)H(2) complexes are very similar to Ir(4)H. Thus binding of an H atom to Ir(4) is only slightly affected by the presence of a second H ligand. In all cases examined, the reaction H(2)+ Ir(4)--> H(2)Ir(4) was found to be exothermic with reaction energies of up to 170 kJ mol(-1).

Computer Simulation↗

C-O bond scission of methoxide on Pd nanoparticles: a density functional study.

C-O bond scission of methoxide species adsorbed at the surface of Pd nanoparticle was studied by DF calculations for the example of cuboctahedral Pd(79). To investigate different locations of adsorbed intermediates as well as the transition state of C-O bond scission, a substrate model was used, which allows one to consider adsorbates without any local geometry restrictions. In contrast to reaction sites on the flat Pd(111) surface and on extended facets, scission of the C-O bond of methoxide at cluster edges is exothermic by approximately 40 kJ mol(-1) and the decomposition product CH(3) is found to be stabilized there. However, the high calculated activation barrier, approximately 140 mol(-1), implies only a very slow reaction compared to dehydrogenation of CH(3)O.

Binding Sites↗

Modeling adsorption of the uranyl dication on the hydroxylated alpha-Al2O3(0001) surface in an aqueous medium. Density functional study.

As a first step toward modeling the interaction of dissolved actinide contaminants with mineral surfaces, we studied low-coverage adsorption of aqueous uranyl, UO2(2+), on the hydroxylated alpha-Al2O3(0001) surface. We carried out density functional periodic slab model calculations and modeled solvation effects by explicit aqua ligands. We explored the formation of both inner- and outer-sphere complexes and estimated the corresponding adsorption energies. Effects of solvation were accounted for by explicit consideration of the first hydration shell of uranyl and by means of a posteriori corrections for long-range solvent effect. With energetics described at the GGA-PW91 level and under the assumption of a fully protonated ideal surface, we predict a weakly bound outer-sphere adsorption complex.

Journal Article↗

A density functional study of uranyl monocarboxylates.

We studied uranium(VI) monocarboxylate complexes by a relativistic density functional method using simple carboxylic acids as ligands, i.e. [UO2(OOCR)]+ (R = H, CH3, CH2CH3). These complexes exist in aqueous solution and, for R = CH3 and CH2CH3, may also be considered as models of uranyl complexated by humic substances. We investigated mono- and bidentate coordination modes. Short-range solvent effects were accounted for explicitly via aqua ligands of the first hydration shell and long-range electrostatic interactions were described via a polarizable continuum model. The calculated results for the uranyl U=O bond, the bond to aqua ligands, and the averaged uranium distances to equatorial oxygen atoms, U-Oeq, agreed quite well with EXAFS-derived interatomic distances. However, the uranyl-carboxylate bond was calculated to be notably shorter than the experimentally determined value. Experimental differences between mono- and bidentate coordination, obtained mainly from crystal structures, were qualitatively reproduced for the U-C distance but not for the average bond length, U-Oeq. We discuss these discrepancies between calculated and experimental results in some detail and suggest changes in the coordination number rather than variations of the coordination geometry as the main source of the experimentally observed variation of the U-Oeq distance.

Carboxylic Acids↗

Structure, stability, electronic and magnetic properties of Ni4 clusters containing impurity atoms.

Using a gradient-corrected density functional method, we studied computationally how single impurity atoms affect the structure and the properties of a Ni4 cluster. H and O atoms coordinate at a Ni-Ni bond, inducing small changes to the structure of bare Ni4 which is essentially a tetrahedron. For a C impurity, we found three stable structures at a Ni4 cluster. In the most stable geometry, the carbon atom cleaves a Ni-Ni bond of Ni4, binding to all Ni atoms. Inclusion of the impurity atom leads to a partial oxidation of the metal atoms and, in the most stable structures, reduces the spin polarization of the cluster compared to bare Ni4. An H impurity interacts mainly with the Ni 4s orbitals, whereas the Ni 3d orbitals participate strongly in the bonding with O and C impurity atoms. For these impurity atoms, Ni 3d contributions dominate the character of the HOMO of the ligated cluster, in contrast to the HOMO of bare Ni4 where Ni 4s orbitals prevail. We also discuss a simple model which relates the effect of a H impurity on the magnetic state of metal clusters to the spin character (minority or majority) of the LUMO or HOMO of the bare metal cluster.

Journal Article↗

Systematic DFT study of gas phase and solvated uranyl and neptunyl complexes [AnO2X4]n (An = U, Np; X = F, Cl, OH, n = -2; X = H2O, n = +2).

The six-valent uranyl and neptunyl complexes [An(VI)O2X4]n (An = U, Np; X = F, Cl, OH, n = -2; X = H2O, n = +2) have been studied within the framework of density functional theory. The relative stabilities of the cis and trans isomers, structural properties, charge distribution, and ligand binding energies have been determined using the modified Perdew-Burke-Ernzerhof functional at the all-electron scalar relativistic level. Uranyl and neptunyl complexes with different ligands have been compared in a systematic fashion, demonstrating close similarity of these actinides in oxidation state VI. In addition, the effect of an aqueous solution has been taken into account with the polarizable continuum model COSMO. Computed averaged ligand binding energies permit one to rationalize the observed different stabilities of the title species in aqueous media.

Journal Article↗

Density-functional model cluster studies of EPR g tensors of Fs+ centers on the surface of MgO.

We report g tensors of surface color centers, so-called F(s) (+) centers, of MgO calculated with two density-functional approaches using accurately embedded cluster models. In line with recent UHV measurements on single-crystalline MgO film, we determined only small g-tensor anisotropies and negative shifts Deltag identical with g-g(e) for all F(s) (+) sites considered, namely, (001)-terrace, step, edge, and corner sites. The g values are very sensitive to the local structure of the defect: relaxation reverses the sign of Deltag. However, accounting for the spin-orbit interaction either self-consistently or perturbatively yields very similar results. In addition to the values of the tensor components, their direction with respect to the surface was determined. In contrast to edges, significant deviations from ideal C(2v) symmetry were found for F(s) (+) centers at steps. Recent data on single-crystalline thin films are reevaluated in the light of these results.

Journal Article↗

Comparison of all sites for Ti substitution in zeolite TS-1 by an accurate embedded-cluster method.

We studied the preferential location of Ti centers in the framework of the Ti-containing MFI zeolite TS-1 using a hybrid DFT/MM embedding method developed recently. This "covalent elastic polarizable environment" (covEPE) cluster embedding allows a complete and self-consistent treatment of solid covalent systems such as zeolites. For the present study, we used a gradient-corrected density functional approach. The resulting structural features of both Si- and Ti-substituted forms of the zeolite framework fit well with available experimental information. The calculated substitution energy of Ti at the 12 crystallographically different tetrahedral sites of the MFI structure vary within 19 kJ/mol with T12 and T2 as most and least preferred sites, respectively. On the basis of these computational results and the preferential sites for Ti substitution reported from different experimental investigations, we concluded that the Ti distribution in the TS-1 framework is not governed by the thermodynamic stability of the pure material.

Journal Article↗

The heat of formation of the uranyl dication: theoretical evaluation based on relativistic density functional calculations.

By using a set of model reactions, we estimated the heat of formation of gaseous UO2(2+) from quantum-chemical reaction enthalpies and experimental heats of formation of reference species. For this purpose, we performed relativistic density functional calculations for the molecules UO2(2+), UO2, UF6, and UF5. We used two gradient-corrected exchange-correlation functionals (revised Perdew-Burke-Ernzerhof (PBEN) and Becke-Perdew (BP)) and we accounted for spin-orbit interaction in a self-consistent fashion. Indeed, spin-orbit interaction notably affects the energies of the model reactions, especially if compounds of U(IV) are involved. Our resulting theoretical estimates for delta fH(o)0 (UO2(2+)), 365+/-10 kcal mol(-1) (PBEN) and 370+/-12 kcal mol(-1) (BP), are in quantitative agreement with a recent experimental result, 364+/-15 kcal mol(-1). Agreement between the results of the two different exchange-correlation functionals PBEN and BP supports the reliability of our approach. The procedure applied offers a general means to derive unknown enthalpies of formation of actinide species based on the available well-established data for other compounds of the element in question.

Journal Article↗

Gold atoms and dimers on amorphous SiO(2): calculation of optical properties and cavity ringdown spectroscopy measurements.

We report on the optical absorption spectra of gold atoms and dimers deposited on amorphous silica in size-selected fashion. Experimental spectra were obtained by cavity ringdown spectroscopy. Issues on soft-landing, fragmentation, and thermal diffusion are discussed on the basis of the experimental results. In parallel, cluster and periodic supercell density functional theory (DFT) calculations were performed to model atoms and dimers trapped on various defect sites of amorphous silica. Optically allowed electronic transitions were calculated, and comparisons with the experimental spectra show that silicon dangling bonds [[triple bond]Si(.-)], nonbridging oxygen [[triple bond]Si-O(.-)], and the silanolate group [[triple bond]Si-O(-)] act as trapping centers for the gold particles. The results are not only important for understanding the chemical bonding of atoms and clusters on oxide surfaces, but they will also be of fundamental interest for photochemical studies of size-selected clusters on surfaces.

Journal Article↗

Efficient treatment of the Hartree interaction in the relativistic Kohn-Sham problem.

We elaborate the two-component Douglas-Kroll reduction of the Dirac-Kohn-Sham problem of relativistic density-functional theory as introduced by Matveev and Rosch [J. Chem. Phys. 118, 3997 (2003)]. That method retains corrections to the Coulomb self-interaction (or Hartree) term of the energy functional that are due to the picture change. Using analytic expressions for the matrix elements, one is able to abandon the resolution of the identity approach for a crucial step of the relativistic transformation. Thus, a major source of uncertainties of the method is eliminated because basis sets no longer have to be extended by functions of higher angular momentum, previously required to ensure kinetic balance. This approach also relies on the electron charge-density fitting scheme via an auxiliary basis set. An efficient approximate implementation results if one restricts the relativistic transformation to the spherically symmetric atom-centered auxiliary functions. It provides accurate results while simplifying greatly the expressions for the matrix elements of the relativistically transformed operators and significantly reducing the computational effort. We demonstrate the performance of the method for the fine structure of one-electron levels of the Hg atom, the g-tensor shifts of NO2, and the properties of the diatomic molecules Bi2, Pb2, PbO, and TlH.

Journal Article↗

Single d-metal atoms on F(s) and F(s+) defects of MgO(001): a theoretical study across the periodic table.

Single d-metal atoms on oxygen defects F(s) and F(s+) of the MgO(001) surface were studied theoretically. We employed an accurate density functional method combined with cluster models, embedded in an elastic polarizable environment, and we applied two gradient-corrected exchange-correlation functionals. In this way, we quantified how 17 metal atoms from groups 6-11 of the periodic table (Cu, Ag, Au; Ni, Pd, Pt; Co, Rh, Ir; Fe, Ru, Os; Mn, Re; and Cr, Mo, W) interact with terrace sites of MgO. We found bonding with F(s) and F(s+) defects to be in general stronger than that with O2- sites, except for Mn-, Re-, and Fe/F(s) complexes. In M/F(s) systems, electron density is accumulated on the metal center in a notable fashion. The binding energy on both kinds of O defects increases from 3d- to 4d- to 5d-atoms of a given group, at variance with the binding energy trend established earlier for the M/O2- complexes, 4d < 3d < 5d. Regarding the evolution of the binding energy along a period, group 7 atoms are slightly destabilized compared to their group 6 congeners in both the F(s) and F(s+) complexes; for later transition elements, the binding energy increases gradually up to group 10 and finally decreases again in group 11, most strongly on the F(s) site. This trend is governed by the negative charge on the adsorbed atoms. We discuss implications for an experimental detection of metal atoms on oxide supports based on computed core-level energies.

Elements↗

Optical absorption spectrum of gold atoms deposited on SiO(2) from cavity ringdown spectroscopy.

The optical properties of gold atoms supported on amorphous silica (alpha-SiO2) were studied experimentally and theoretically in the visible range. Samples were prepared in situ by depositing Au atoms at low coverages (5 x 10(12) cm(-2)) in UHV, and the optical absorption spectra were recorded by cavity ringdown spectroscopy. The atomic absorption bands can be attributed to gold atoms trapped at [triple bond] Si-O(.-) and [triple bond]Si-O(-) defect sites. The absence of optical transitions typical for Au(2) shows that the atoms are efficiently anchored at these defect sites, preventing diffusion and aggregation. Furthermore, these experimental results reveal that it is now possible to study optical properties of well-defined nanostructures at surface coverages as low as 5 x 10(11) cm(-2).

Journal Article↗