TinO2n-1 Magneli phases studied using density functional theory

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Capacitance spectroscopy and density functional theory

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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.

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.

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.

Comparison of density functional theory predictions of gas-phase deprotonation data

The SVWN, BVWN, BP86, BLYP, BPW91, B3P86, B3LYP, B3PW91, B1LYP, mPW1PW, and PBE1PBE density functionals, as implemented in Gaussian 98 and Gaussian 03, were used to calculate ΔG0 and ΔH0 values for 17 deprotonation reactions where the experimental values are accurately known. The PBE1PBE and B3P86 functionals are shown to compute results with accuracy comparable to more computationally intensive compound model chemistries. A rationale for the relative performance of various functionals is explored.

Comparison of Model Chemistry and Density Functional Theory Thermochemical Predictions with Experiment for Formation of Ionic Clusters of the Ammonium Cation Complexed with Water and Ammonia; Atmospheric Implications

The G2, G3, CBS-QB3, and CBS-APNO model chemistry methods and the B3LYP, B3P86, mPW1PW, and PBE1PBE density functional theory (DFT) methods have been used to calculate ΔH° and ΔG° values for ionic clusters of the ammonium ion complexed with water and ammonia. Results for the clusters NH4+(NH3)n and NH4+(H2O)n, where n = 1−4, are reported in this paper and compared against experimental values. Agreement with the experimental values for ΔH° and ΔG° for formation of NH4+(NH3)n clusters is excellent. Comparison between experiment and theory for formation of the NH4+(H2O)n clusters is quite good considering the uncertainty in the experimental values. The four DFT methods yield excellent agreement with experiment and the model chemistry methods when the aug-cc-pVTZ basis set is used for energetic calculations and the 6-31G* basis set is used for geometries and frequencies. On the basis of these results, we predict that all ions in the lower troposphere will be saturated with at least one complete first hydration shell of water molecules.

DENSITY OF AMORPHOUS CARBON BY USING DENSITY FUNCTIONAL THEORY

Amorphous carbon has been investigated for a long time. Since it has the random orientation of carbon atoms, its density depends on the position of each carbon atom. It is important to know the density of amorphous carbon to use it for modeling advance carbon materials in the future. Two methods were used to create the initial structures of amorphous carbon. One is the random placement method by randomly locating 100 carbon atoms in a cubic lattice. Another method is the liquid-quench method by using reactive force field (ReaxFF) to rapidly decrease the system of 100 carbon atoms from the melting temperature. Density functional theory (DFT) was used to refine the position of each carbon atom and the dimensions of the boundaries to minimize the ground energy of the structure. The average densities of amorphous carbon structures created by the random placement method and the liquid-quench method are 2.59 and 2.44 g/cm3, respectively. Both densities have a good agreement with previous works. In addition, the final structure of amorphous carbon generated by the liquid-quench method has lower energy.