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.

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.

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.

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.

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.

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.

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.

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.

Calculations of the spectral shifts and line profiles of alkaline earth atoms in liquid helium environment

Liquid configurations generated by Metropolis Monte Carlo simulations are used in time-dependent density functional theory calculations of the spectral line shifts and line profiles of the lowest lying excitation of the alkaline earth atoms, Be, Mg, Ca, Sr and Ba embedded in liquid helium. The results are in very good agreement with the available experimental data. Special attention is given to the calculated spectroscopic shift and the associated line broadening. The analysis specifies the inhomogeneous broadening of the three separate contributions due to the splitting of the s -> p transition of the alkaline earth atom in the liquid environment. (C) 2012 Elsevier B. V. All rights reserved.

X-ray diffraction and theoretical investigation of the Gedunin crystal structure

The Gedunin compound (C28H34O6) is a natural product extracted from Trichilia pallida that has shown a wide activity. The crystallographic structure shows two conformers in the asymmetric unit, which differ in a rotation of the furan group. To understand this molecular arrangement, the density functional calculations. Molecular Electrostatic Potential (MEP) and thermodynamic function calculation have been performed at the B3LYP/6-311++g(d,p) level. Both conformers were optimized and the agreement with the experimental structure was very good, making possible further theoretical analysis of the structure. The inter-conversion between two conformers depends on the energy barrier. This process is studied in the vacuum and shows two transition states with a low energetic barrier for a potential energy curve scanning rigid around furan group: 4.37 kcal/mol and 16.52 kcal/mol. As the first transition state has a notably lower energetic barrier, the preferred inter-conversion pathway between the conformers involves the first rather than the second transition state. Understanding this transition state in detail led us to perform its optimization, showing an energetic barrier around 3.66 kcal/mol. The negative free energy and low enthalpy confirm that the process is spontaneous and exothermic. The results show that this requirement makes the existence of the two conformers in the asymmetric unit possible. The structure of molecules in the asymmetric unit is better understood when the MEP is used on the interaction between molecules. For Gedunin, both molecules have shown MEP with well-defined regions, and this behavior contributes to the observed link between molecules and for the negative regions complementing positive regions of another molecule. (C) 2011 Elsevier B.V. All rights reserved.

Spectroscopic and Catalytic Characterization of a Functional (FeFeII)-Fe-III Biomimetic for the Active Site of Uteroferrin and Protein Cleavage

A mixed-valence complex, [Fe(III)Fe(II)L1(mu-OAc)(2)]BF4 center dot H2O, where the ligand H(2)L1 = 2-{[[3-[((bis-(pyridin-2-ylmethyl)amino)methyl)-2-hydroxy-5-methylbenzyl](pyridin-2-ylmethyl)amino]methyl]phenol}, has been studied with a range of techniques, and, where possible, its properties have been compared to those of the corresponding enzyme system purple acid phosphatase. The (FeFeII)-Fe-III and Fe-2(III) oxidized species were studied spectroelectrochemically. The temperature-dependent population of the S = 3/2 spin states of the heterovalent system, observed using magnetic circular dichroism, confirmed that the dinuclear center is weakly antiferromagnetically coupled (H = -2JS(1).S-2, where J = -5.6 cm(-1)) in a frozen solution. The ligand-to-metal charge-transfer transitions are correlated with density functional theory calculations. The (FeFeII)-Fe-III complex is electron paramagnetic resonance (EPR)-silent, except at very low temperatures (<2 K), because of the broadening caused by the exchange coupling and zero-field-splitting parameters being of comparable magnitude and rapid spin-lattice relaxation. However, a phosphate-bound Fe-2(III) complex showed an EPR spectrum due to population of the S-tot = 3 state (J= -3.5 cm(-1)). The phosphatase activity of the (FeFeII)-Fe-III complex in hydrolysis of bis(2,4-dinitrophenyl)phosphate (k(cat.) = 1.88 x 10(-3) s(-1); K-m = 4.63 x 10(-3) mol L-1) is similar to that of other bimetallic heterovalent complexes with the same ligand. Analysis of the kinetic data supports a mechanism where the initiating nucleophile in the phosphatase reaction is a hydroxide, terminally bound to Fe-III. It is interesting to note that aqueous solutions of [Fe(III)Fe(II)L1(mu-OAc)(2)](+) are also capable of protein cleavage, at mild temperature and pH conditions, thus further expanding the scope of this complex's catalytic promiscuity.

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]

Polyazomethine with vinylene and phenantridine moieties in the main chain: Synthesis, characterization, opto(electrical) properties and theoretical calculations

The opto(electrical) properties and theoretical calculations of polyazomethine with vinylene and phenantridine moieties in the main chain were investigated in the present study. 2,5-Bis(hexyloxy)-1,4-bis[(2,5-bis(hexyloxy)-4-formyl-phenylenevinylene]benzene was polymerized in solution with 3,8-diamino-6-phenylphenanthridine (PAZ-PV-Ph). The temperatures of 5% weight loss (T-5%) of the polyazomethine was observed at 356 degrees C in nitrogen. Electrochemical properties of thin film of the polymer were studied by differential pulse voltammetry. The HOMO level of the PAZ-PV-Ph was at -4.97 eV. The energy band gap (E-g) was detected of approximately similar to 1.9 eV. Energy band gap (E-gopt) was additionally calculated from absorption spectrum and absorption coefficient alpha. The absorption UV-vis spectra of polyazomethine recorded in solution showed a blue shift in comparison with the solid state. HOMO-LUMO levels and E-g were additionally calculated theoretically by density functional theory and molecular simulations of PAZ-PV-Ph are presented. Current density-voltage (J-U) measurements were performed on ITO/PAZ-PV-Ph/Al, ITO/TiO2/PAZ-PV-Ph/Al and ITO/PEDOT/PAZ-PV-Ph:TiO2/Al devices in the dark and during irradiation with light (under illumination of 1000 W m(-2)). The polymer was tested using AFM technique and roughness (R-a, R-ms) along with skew and kurtosis are presented.

Raman spectroscopy and DFT calculations of para-coumaric acid and its deprotonated species

An electronic and vibrational spectroscopic analysis of p-coumaric acid (HCou) and its deprotonated species was performed by UV-vis and Raman, respectively, and the results were supported by density functional theory (OFT) calculations. Electronic UV-vis spectral data of HCou solutions show that the deprotonation of the carboxyl group (Cou(-)) leads to a blue shift of the lowest energy electronic transition in comparison to the neutral species, whereas the subsequent deprotonation of the phenolic moiety (Cou(2-)) carries out to a more delocalized chromophore. The DFT geometric parameters calculations suggest that the variation in the electronic delocalization for the three organic species is due to different contribution of a quinoid structure that is significantly distorted in the case of Cou(2-). The Raman data of HCou and its sodium salts show that the main spectral features that allow to differentiate the three organic species are those involving the styrene nu(C=C)(sty) vibration at 1600cm(-1) region. Even though the Raman spectra of the sodium salts of Cou(-) and Cou(2-) anions show subtle differences, the appearing of a band at ca. 1598cm(-1) in the Na(2)Cou spectrum, assigned to a mode involving the carboxylate asymmetric stretching, nu(as)(COO), and the styrene stretching, nu(C=C)(sty), is quite characteristic, as confirmed by the theoretical Raman spectrum. Considering that p-coumaric acid is an archetypical phenolic compound with several biological activities that essentially depend upon the medium pH, Raman spectroscopy results reported in this work can provide a proper way to characterize such important phytochemical compound in different protonation states. In order to complement the characterization of the sodium salts, X-ray diffraction (XRD) and thermal analysis were performed. (C) 2011 Elsevier B.V. All rights reserved.

First principles calculations for vacancies and antisites in PbSe and PbTe: bulk and nanowire.

The energetic stability and the electronic properties of vacancies (VX) and antisites (XY) in PbSe and PbTe are investigated. PbSe and PbTe are narrow band gap semiconductors and have the potential to be used in infrared detectors, laser, and diodes. They are also of special interest for thermoelectric devices (TE). The calculations are based in the Density Functional Theory (DFT) and the General Gradient Approximation (GGA) for the exchange-correlation term, as implemented in the VASP code. The core and valence electrons are described by the Projected Augmented Wave (PAW) and the Plane Wave (PW) methods, respectively. The defects are studied in the bulk and nanowire (NW) system. Our results show that intrinsec defects (vacancies and antisites) in PbTe have lower formation energies in the NW as compared to the bulk and present a trend in migrate to the surface of the NW. For the PbSe we obtain similar results when compare the formation energy for the bulk and NW. However, the Pb vacancy and the antisites are more stable in the core of the NW. The intrinsec defects are shallow defects for the bulk system. For both PbSe and PbTe VPb is a shallow acceptor defect and VSe and VT e are shallow donor defects for the PbSe and PbTe, respectively. Similar electronic properties are observed for the antisites. For the Pb in the anion site we obtain an n-type semiconductor for both PbSe and PbTe, SeP b is a p-type for the PbSe, and T eP b is a n-type for PbTe. Due the quantum con¯nement effects present in the NW (the band gap open), these defects have different electronic properties for the NW as compared to the bulk. Now these defects give rise to electronic levels in the band gap of the PbTe NW and the VT e present a metallic character. For the PbSe NW a p-type and a n-type semiconductor is obtained for the VP b and P bSe, respectively. On the other hand, deep electronic levels are present in the band gap for the VSe and SePb. These results show that due an enhanced in the electronic density of states (DOS) near the Fermi energy, the defective PbSe and PbTe are candidates for efficient TE devices.