Ethanol and Water Adsorption on Close-Packed 3d, 4d, and 5d Transition-Metal Surfaces: A Density Functional Theory Investigation with van der Waals Correction

Nowadays, there is a great interest in the economic success of direct ethanol fuel cells; however, our atomistic understanding of the designing of stable and low-cost catalysts for the steam reforming of ethanol is still far from satisfactory, in particular due to the large number of undesirable intermediates. In this study, we will report a first-principles investigation of the adsorption properties of ethanol and water at low coverage on close-packed transition-metal (TM) surfaces, namely, Fe(110), Co(0001), Ni(111), Cu(111), Ru(0001), Rh(111), Pd(111), Ag(111), Os(0001), Ir(111), Pt(111), and Au(111), employing density functional theory (DFT) calculations. We employed the generalized gradient approximation with the formulation proposed by Perdew, Burke, and Erzenholf (PBE) to the exchange correlation functional and the empirical correction proposed by S. Grimme (DFT+D3) for the van der Waals correction. We found that both adsorbates binds preferentially near or on the on top sites of the TM surfaces through the 0 atoms. The PBE adsorption energies of ethanol and water decreases almost linearly with the increased occupation of the 4d and 5d d-band, while there is a deviation for the 3d systems. The van der Waals correction affects the linear behavior and increases the adsorption energy for both adsorbates, which is expected as the van der Waals energy due to the correlation effects is strongly underestimated by DFT-PBE for weak interacting systems. The geometric parameters for water/TM are not affected by the van der Waals correction, i.e., both DFT and DFT+D3 yield an almost parallel orientation for water on the TM surfaces; however, DFT+D3 changes drastically the ethanol orientation. For example, DFT yields an almost perpendicular orientation of the C-C bond to the TM surface, while the C-C bond is almost parallel to the surface using DFT +D3 for all systems, except for ethanol/Fe(110). Thus, the van der Waals correction decreases the distance of the C atoms to the TM surfaces, which might contribute to break the C-C bond. The work function decreases upon the adsorption of ethanol and water, and both follow the same trends, however, with different magnitude (larger for ethanol/TM) due to the weak binding of water to the surface. The electron density increases mainly in the region between the topmost layer and the adsorbates, which explains the reduction of the substrate work function.

Platinum-Based Nanoalloys PtnTM55-n (TM=Co, Rh, Au): A Density Functional Theory Investigation

Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.

Interactions between 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and gamma-Al2O3 (100) surface calculated by density functional theory

The main aim of this work is to investigate the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4C1Im]+[BF4]-) ionic liquid (IL) adsorption on the gamma-Al2O3 (100) by density functional theory calculations to try to rationalize the adsorption as an electrostatic phenomenon. Optimized geometries and interaction energies of IL one-monolayer on the gamma-Al2O3 were obtained on high surface coverage (one cationanion pair per 94.96 nm2). A study of dispersion force was made to estimate its contribution to the adsorption. Overall, the process is ruled by electrostatic interaction between ions and surface. Adsorption of the anion [BF4]- and cation [C4C1Im]+ was also studied by Bader charge analysis and charge density difference for supported and unsupported situations. It is suggested that the IL ions have their charges maintained with significant anion cloud polarization inward to the acid aluminum sites. (c) 2012 Wiley Periodicals, Inc.

Encapsulation of small magnetic clusters in fullerene cages: A density functional theory investigation within van der Waals corrections

The encapsulation of magnetic transition-metal (TM) clusters inside carbon cages (fullerenes, nanotubes) has been of great interest due to the wide range of applications, which spread from medical sensors in magnetic resonance imaging to photonic crystals. Several theoretical studies have been reported; however, our atomistic understanding of the physical properties of encapsulated magnetic TM 3d clusters is far from satisfactory. In this work, we will report general trends, derived from density functional theory within the generalized gradient approximation proposed by Perdew, Burke, and Ernzerhof (PBE), for the encapsulation properties of the TMm@C-n (TM = Fe, Co, Ni; m = 2-6, n = 60,70,80,90) systems. Furthermore, to understand the role of the van der Waals corrections to the physical properties, we employed the empirical Grimme's correction (PBE + D2). We found that both PBE and PBE + D2 functionals yield almost the same geometric parameters, magnetic and electronic properties, however, PBE + D2 strongly enhances the encapsulation energy. We found that the center of mass of the TMm clusters is displaced towards the inside C-n surfaces, except for large TMm clusters (m = 5 and 6). For few cases, e. g., Co-4 and Fe-4, the encapsulation changes the putative lowest-energy structure compared to the isolated TMm clusters. We identified few physical parameters that play an important role in the sign and magnitude of the encapsulation energy, namely, cluster size, fullerene equatorial diameter, shape, curvature of the inside C-n surface, number of TM atoms that bind directly to the inside C-n surface, and the van der Waals correction. The total magnetic moment of encapsulated TMm clusters decreases compared with the isolated TMm clusters, which is expected due to the hybridization of the d-p states, and strongly depends on the size and shape of the fullerene cages.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

Adsorption of NO on the Rh-13, Pd-13, Ir-13, and Pt-13 Clusters: A Density Functional Theory Investigation

The adsorption of NO on transition-metal (TM) surfaces has been widely studied by experimental and theoretical techniques; however, our atomistic understanding of the interaction of nitrogen monoxide (NO) with small TM clusters is far from satisfactory, which compromises a deep understanding of real catalyst devices. In this study, we report a density functional theory study of the adsorption properties of NO on the TM13 (TM = Rh, Pd, Ir, Pt) clusters employing the projected augmented wave method. We found that the interaction of NO with TM13 is much more complex than that for NO/TM(111). In particular, for low symmetry TM13 clusters, there is a strong rearrangement of the electronic charge density upon NO adsorption and, as a consequence, the adsorption energy shows a very complex dependence even for adsorption sites with the same local effective coordination. We found a strong enhancement of the binding energy of NO to the TM13 clusters compared with the TM(111) surfaces, as the antibonding NO states are not occupied for NO/TM13, and the general relationship based on the d-band model between adsorption energy and the center of gravity of the occupied d-states does not hold for the studied TM13 clusters, in particular, for clusters with low symmetry. In contrast with the adsorption energy trends, the geometric NO/TM13 parameters and the vibrational N-O frequencies for different coordination sites follow the same trend as for the respective TM(111) surfaces, while the changes in the frequencies between different surfaces and TM13 clusters reflect the strong NO-TM13 interaction.

Density Functional Theory (DFT) Study of Edaravone Derivatives as Antioxidants

Quantum chemical calculations at the B3LYP/6-31G* level of theory were employed for the structure-activity relationship and prediction of the antioxidant activity of edaravone and structurally related derivatives using energy (E), ionization potential (IP), bond dissociation energy (BDE), and stabilization energies(Delta E-iso). Spin density calculations were also performed for the proposed antioxidant activity mechanism. The electron abstraction is related to electron-donating groups (EDG) at position 3, decreasing the IP when compared to substitution at position 4. The hydrogen abstraction is related to electron-withdrawing groups (EDG) at position 4, decreasing the BDECH when compared to other substitutions, resulting in a better antioxidant activity. The unpaired electron formed by the hydrogen abstraction from the C-H group of the pyrazole ring is localized at 2, 4, and 6 positions. The highest scavenging activity prediction is related to the lowest contribution at the carbon atom. The likely mechanism is related to hydrogen transfer. It was found that antioxidant activity depends on the presence of EDG at the C-2 and C-4 positions and there is a correlation between IP and BDE. Our results identified three different classes of new derivatives more potent than edaravone.

Light emission and current rectification in a molecular device: Experiment and theory

We have investigated optical and transport properties of the molecular structure 2,3,4,5-tetraphenyl-1-phenylethynyl-cyclopenta-2,4-dienol experimentally and theoretically. The optical spectrum was calculated using Hartree-Fock-intermediate neglect of differential overlap-configuration interaction model. The experimental photoluminescence spectrum showed a peak around 470nm which was very well described by the modeling. Electronic transport measurements showed a diode-like effect with a strong current rectification. A phenomenological microscopic model based on non-equilibrium Green's function technique was proposed and a very good description electronic transport was obtained. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4767457]

Hybrid density functional study of small Rh-n (n=2-15) clusters

The physical properties of small rhodium clusters, Rh-n, have been in debate due to the shortcomings of density functional theory (DFT). To help in the solution of those problems, we obtained a set of putative lowest energy structures for small Rh-n (n = 2-15) clusters employing hybrid-DFT and the generalized gradient approximation (GGA). For n = 2-6, both hybrid and GGA functionals yield similar ground-state structures (compact), however, hybrid favors compact structures for n = 7-15, while GGA favors open structures based on simple cubic motifs. Thus, experimental results are crucial to indicate the correct ground-state structures, however, we found that a unique set of structures (compact or open) is unable to explain all available experimental data. For example, the GGA structures (open) yield total magnetic moments in excellent agreement with experimental data, while hybrid structures (compact) have larger magnetic moments compared with experiments due to the increased localization of the 4d states. Thus, we would conclude that GGA provides a better description of the Rh-n clusters, however, a recent experimental-theoretical study [ Harding et al., J. Chem. Phys. 133, 214304 (2010)] found that only compact structures are able to explain experimental vibrational data, while open structures cannot. Therefore, it indicates that the study of Rh-n clusters is a challenging problem and further experimental studies are required to help in the solution of this conundrum, as well as a better description of the exchange and correlation effects on the Rh n clusters using theoretical methods such as the quantum Monte Carlo method.

Ab-initio calculations for a realistic sensor: A study of CO sensors based on nitrogen-rich carbon nanotubes

The use of nanoscale low-dimensional systems could boost the sensitivity of gas sensors. In this work we simulate a nanoscopic sensor based on carbon nanotubes with a large number of binding sites using ab initio density functional electronic structure calculations coupled to the Non-Equilibrium Green's Function formalism. We present a recipe where the adsorption process is studied followed by conductance calculations of a single defect system and of more realistic disordered system considering different coverages of molecules as one would expect experimentally. We found that the sensitivity of the disordered system is enhanced by a factor of 5 when compared to the single defect one. Finally, our results from the atomistic electronic transport are used as input to a simple model that connects them to experimental parameters such as temperature and partial gas pressure, providing a procedure for simulating a realistic nanoscopic gas sensor. Using this methodology we show that nitrogen-rich carbon nanotubes could work at room temperature with extremely high sensitivity. Copyright 2012 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. [http://dx.doi.org/10.1063/1.4739280]

Hydrogen bond and the resonance effect on the formamide-water complexes

The interaction of formamide and the two transition states of its amide group rotation with one, two, or three water molecules was studied in vacuum. Great differences between the electronic structure of formamide in its most stable form and the electronic structure of the transition states were noticed. Intermolecular interactions were intense, especially in the cases where the solvent interacted with the amide and the carbonyl groups simultaneously. In the transition states, the interaction between the lone pair of nitrogen and the water molecule becomes important. With the aid of the natural bond orbitals, natural resonance theory, and electron localization function (ELF) analyses an increase in the resonance of planar formamide with the addition of successive water molecules was observed. Such observation suggests that the hydrogen bonds in the formamidewater complexes may have some covalent character. These results are also supported by the quantitative ELF analyses. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

Isatin-Schiff base copper(II) complexes-A DFT study of the metal-ligand bonding situation

Herein, we report results of calculations based on density functional theory (BP86/TZVP) of a set of isatin-Schiff base copper(II) and related complexes, 1-12, that have shown significant pro-apoptotic activity toward diverse tumor cells. The interaction of the copper(II) cation with different ligands has been investigated at the same level of theory. The strength and character of the Cu(II)-L bonding was characterized by metal-ligand bond lengths, vibrational frequencies, binding energies, ligand deformation energies, and natural population analysis. The metal-ligand bonding situation was also characterized by using two complementary topological approaches, the quantum theory of atoms-in-molecules (QTAIM) and the electron localization function (ELF). The calculated electronic g-tensor and hyperfine coupling constants present significant agreement with the EPR experimental data. The calculated parameters pointed to complex 10 as the most stable among the isatin-Schiff base copper(II) species, in good agreement with experimental data that indicate this complex as the most reactive in the series. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012

The Role of Molecular Conformation and Polarizable Embedding for One- and Two-Photon Absorption of Disperse Orange 3 in Solution

Solvent effects on the one- and two-photon absorption (IPA and 2PA) of disperse orange 3 (DO3) in dimethyl sulfoxide (DMSO) are studied using a discrete polarizable embedding (PE) response theory. The scheme comprises a quantum region containing the chromophore and an atomically granulated classical region for the solvent accounting for full interactions within and between the two regions. Either classical molecular dynamics (MD) or hybrid Car-Parrinello (CP) quantum/classical (QM/MM) molecular dynamics simulations are employed to describe the solvation of DO3 in DMSO, allowing for an analysis of the effect of the intermolecular short-range repulsion, long-range attraction, and electrostatic interactions on the conformational changes of the chromophore and also the effect of the solute-solvent polarization. PE linear response calculations are performed to verify the character, solvatochromic shift, and overlap of the two lowest energy transitions responsible for the linear absorption spectrum of DO3 in DMSO in the visible spectral region. Results of the PE linear and quadratic response calculations, performed using uncorrelated solute-solvent configurations sampled from either the classical or hybrid CP QM/MM MD simulations, are used to estimate the width of the line shape function of the two electronic lowest energy excited states, which allow a prediction of the 2PA cross-sections without the use of empirical parameters. Appropriate exchange-correlation functionals have been employed in order to describe the charge-transfer process following the electronic transitions of the chromophore in solution.

Ligand- and Structure-Based Drug Design Strategies and PPAR delta/alpha Selectivity

Peroxisome-proliferator-activated receptors are a class of nuclear receptors with three subtypes: a, ? and d. Their main function is regulating gene transcription related to lipid and carbohydrate metabolism. Currently, there are no peroxisome-proliferator-activated receptors d drugs being marketed. In this work, we studied a data set of 70 compounds with a and d activity. Three partial least square models were created, and molecular docking studies were performed to understand the main reasons for peroxisome-proliferator-activated receptors d selectivity. The obtained results showed that some molecular descriptors (log P, hydration energy, steric and polar properties) are related to the main interactions that can direct ligands to a particular peroxisome-proliferator-activated receptors subtype.

Two-photon absorption in oxazole derivatives: An experimental and quantum chemical study

Experimental and theoretical studies on the two-photon absorption properties of two oxazole derivatives: 2,5-diphenyloxazole (PPO) and 2-(4-biphenylyI)-5-phenyl-1,3,4-oxadiazole (PBD) are presented. The two-photon absorption cross-section spectra were determined by means of the Z-scan technique, from 460 up to 650 nm, and reached peak values of 84 GM for PBD and 27 GM for PPO. Density Functional Theory and response function formalism are used to determine the molecular structures and the one- and two-photon absorption properties and to assist in the interpretation of the experimental results. The Polarizable Continuum Model in one-photon absorption calculations is used to estimate solvent effects. (C) 2011 Elsevier B.V. All rights reserved.