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.

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.

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.

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.

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.

Spectroscopic, Structural, and Conformational Properties of (Z)-4,4, 4-Trifluoro-3-(2-hydroxyethylamino)-1-(2-hydroxypheny1)-2-buten-1-one, C12H12F3NO3: A Trifluoromethyl-Substituted beta-Aminoenone

The ( Z)-4,4,4-trifluoro-3-(2-hydroxyethylamino)-1-(2-hydroxyphenyl)-2-buten-1-one (C12H12F3NO3) compound was thoroughly studied by IR, Raman, UV-visible, and C-13 and F-19 NMR spectroscopies. The solid-state molecular structure was determined by X-ray diffraction methods. It crystallizes in the P2(1)/c space group with a = 12.1420(4) angstrom, b = 7.8210(3) angstrom, c := 13.8970(5) angstrom, beta = 116.162(2)degrees, and Z = 4 molecules per unit cell. The molecule shows a nearly planar molecular skeleton, favored by intramolecular OH center dot center dot center dot 0 and NH center dot center dot center dot 0 bonds, which are arranged in the lattice as an OH center dot center dot center dot 0 bonded polymer coiled around crystallographic 2-fold screw-axes. The three postulated tautomers were evaluated using quantum chemical calculations. The lowest energy tautomer (I) calculated with density functional theory methods agrees with the observed crystal structure. The structural and conformational properties are discussed considering the effect of the intra- and intermolecular hydrogen bond interactions.

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.

Theoretical models for the antitrypanosomal activity of thiosemicarbazone derivatives

Thiosemicarbazones are cruzain inhibitors which have been identified as potential antitrypanosomal agents. In this work, several molecular properties were calculated at the density functional theory (DFT)/B3LYP/6-311G* level for a set of 44 thiosemicarbazones. Unsupervised and supervised pattern recognition techniques (hierarchical cluster analysis, principal component analysis, kth-nearest neighbors, and soft independent modeling by class analogy) were used to obtain structureactivity relationship models, which are able to classify unknown compounds according to their activities. The chemometric analyses performed here revealed that 12 descriptors can be considered responsible for the discrimination between high and low activity compounds. Classification models were validated with an external test set, showing that predictive classifications were achieved with the selected variable set. The results obtained here are in good agreement with previous findings from the literature, suggesting that our models can be useful on further investigations on the molecular determinants for the antichagasic activity. (C) 2012 Wiley Periodicals, Inc.

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]

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

The 2:1 cycloadducts from [3+2] 1,3-dipolar cycloaddition of nitrile oxide and vinylacetic acid. Synthesis and liquid crystal behaviour

Four liquid crystals (LC) 3,7a-bis(4-alkyloxyphenyl)-7,7a-dihydro-6H-isoxazolo[2,3-d][1,2,4]oxadiazol-6-yl)acetic acid (7a-d) were synthesised and the mesomorphic behaviour reported. The LCs were characterised as 2: 1 bisadducts, which were obtained from a double [3+2] 1,3-dipolar cycloaddition. In the first step, the cycloaddition of 4-alkyloxyphenylnitrile oxide (4a-d) and vinylacetic acid (5) gave the initial unobserved 1:1 cycloadducts 2-[3-(4-alkyloxyphenyl)-4,5-dihydroisoxazol-5-yl]acetic acid (6a-d). In the second step, the addition of a second equivalent of 4 to 6 yielded the 2: 1 bisadducts 7a-d without any traces of 6. All compounds 7a-d were unstable during the transition from the mesophase to the isotropic state upon first heating as evidenced by the large peaks in the differential scanning calorimetry traces. Due to the chemical instability of the compounds upon heating, the transition temperature related to the smectic C to smectic A transitions was acquired by means of an image processing method. X-Ray diffraction experiments were also used to analyse the liquid-crystalline phases. A theoretical calculation was performed using density functional theory (DFT) methods at the PBE1PBE/6-311+G(2d,p) level (with solvent effect) in order to get information about the energetic profile of the 2: 1 cycloaddition. DFT studies revealed that the cycloaddition process is controlled by the HOMO(dipolarophile) - LUMO(1,3-dipole), and that the double [3+2] 1,3-dipolar cycloaddition reaction is quite possible.

Rare-earth impurities in gallium nitride: The role of the Hubbard potential

We performed a first principles investigation on the electronic properties of 4f-rare earth substitutional impurities in zincblende gallium nitride (GaN:REGa, with RE=Eu, Gd, Tb, Dy, Ho, Er and Tm). The calculations were performed within the all electron methodology and the density functional theory. We investigated how the introduction of the on-site Hubbard U potential (GGA + U) corrects the electronic properties of those impurities. We showed that a self-consistent procedure to compute the Hubbard potential provides a reliable description on the position of the 4f-related energy levels with respect of the GaN valence band top. The results were compared to available data coming from a recent phenomenological model. (C) 2012 Elsevier B.V. All rights reserved.