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]

ESTIMATING REACTION CONSTANTS BY AB INITIO MOLECULAR MODELING: A STUDY ON THE OXIDATION OF PHENOL TO CATECHOL AND HYDROQUINONE IN ADVANCED OXIDATION PROCESSES

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

Estimating reaction constants by ab initio molecular modeling: a study on the oxidation of phenol to catechol and hydroquinone in advanced oxidation processes

Molecular modeling is growing as a research tool in Chemical Engineering studies, as can be seen by a simple research on the latest publications in the field. Molecular investigations retrieve information on properties often accessible only by expensive and time-consuming experimental techniques, such as those involved in the study of radical-based chain reactions. In this work, different quantum chemical techniques were used to study phenol oxidation by hydroxyl radicals in Advanced Oxidation Processes used for wastewater treatment. The results obtained by applying a DFT-based model showed good agreement with experimental values available, as well as qualitative insights into the mechanism of the overall reaction chain. Solvation models were also tried, but were found to be limited for this reaction system within the considered theoretical level without further parameterization.

Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

New insights on reaction pathway selectivity promoted by crown ether phase-transfer catalysis: Model ab initio calculations of nucleophilic fluorination

Crown ethers have the ability of solubilizing inorganic salts in apolar solvents and to promote chemical reactions by phase-transfer catalysis. However, details on how crown ethers catalyze ionic S(N)2 reactions and control selectivity are not well understood. In this work, we have used high level theoretical calculations to shed light on the details of phase-transfer catalysis mechanism of KF reaction with alkyl halides promoted by 18-crown-6. A complete analysis of the of the model reaction between KF(18-crown-6) and ethyl bromide reveals that the calculations can accurately predict the product ratio and the overall kinetics. Our results point out the importance of the K* ion and of the crown ether ring in determining product selectivity. While the K* ion favors the S(N)2 over the E2 anti pathway, the crown ether ring favors the S(N)2 over E2 syn route. The combination effects lead to a predicted 94% for the S(N)2 pathway in excellent agreement with the experimental value of 92%. A detailed analysis of the overall mechanism of the reaction under phase-transfer conditions also reveals that the KBr product generated in the nucleophilic fluorination acts as an inhibitor of the 18-crown-6 catalyst and it is responsible for the observed slow reaction rate. (C) 2012 Elsevier B.V. All rights reserved.

Ab initio studies of pristine and oxidized Cu3Au(100) and (111) surfaces

We review the previous literature and our recent work on first-principles studies of Cu3Au(100) and (111) surfaces, with focus on the segregation of atomic species to the surface at pristine conditions and in the presence of oxygen. In particular, the combined use of experimental and theoretical tools to achieve chemical identification at an atomic level of the surface species is emphasized and discussed.

The Structural and Electronic Properties of Tin Oxide Nanowires: An Ab Initio Investigation

We performed an ab initio investigation on the properties of rutile tin oxide (SnOx) nanowires. We computed the wire properties determining the equilibrium geometries, binding energies, and electronic band structures for several wire dimensions and surface facet configurations. The results allowed us to establish scaling laws for the structural properties, in terms of the nanowire perimeters. The results also showed that the surface states control most of the electronic properties of the nanowires. Oxygen incorporation in the nanowire surfaces passivated the surface-related electronic states, and the resulting quantum properties and scaling laws were fully consistent with electrons confined inside the nanowire. Additionally, oxygen incorporation in the wire surfaces generated an unbalanced concentration of spin up and down electrons, leading to magnetic states for the nanowires.

The T-2 phase in the Nb-Si-B system studied by ab initio calculations and synchrotron X-ray diffraction

The solid solution based on Nb5Si3 (Cr5B3 structure type, D8(l), tl32, 14/mcm, No140, a=6.5767 angstrom, c=11.8967 angstrom) in the Nb-Si-B system was studied from the structural and thermodynamic point of view both experimentally and by ab initio calculations. Rietveld refinement of powder X-ray synchrotron data allowed to determine the boron to silicon substitution mechanism and the structural parameters. Ab initio calculations of different ordered compounds and selected disordered alloys allowed to obtain in addition to the enthalpy of formation of the solution, substitution mechanism and structural parameters which are in excellent agreement with the experimental data. The stability of the phase is discussed. (C) 2012 Elsevier Inc. All rights reserved.

Study of the oxygen vacancy influence on magnetic properties of Fe- and Co-doped SnO2 diluted alloys

Transition-metal (TM)-doped diluted magnetic oxides (DMOs) have attracted attention from both experimental and theoretical points of view due to their potential use in spintronics towards new nanostructured devices and new technologies. In the present work, we study the magnetic properties of Sn0.96TM0.04O2 and Sn0.96TM0.04O1.98(V (O))(0.02), where TM = Fe and Co, focusing in particular in the role played by the presence of O vacancies nearby the TM. The calculated total energy as a function of the total magnetic moment per cell shows a magnetic metastability, corresponding to a ground state, respectively, with 2 and 1 mu(B)/cell, for Fe and Co. Two metastable states, with 0 and 4 mu(B)/cell were found for Fe, and a single value, 3 mu(B)/cell, for Co. The spin-crossover energies (E (S)) were calculated. The values are E (S) (0/2) = 107 meV and E (S) (4/2) = 25 meV for Fe. For Co, E (S) (3/1) = 36 meV. By creating O vacancies close to the TM site, we show that the metastablity and E (S) change. For iron, a new state appears, and the state with zero magnetic moment disappears. The ground state is 4 mu(B)/cell instead of 2 mu(B)/cell, and the energy E (S) (2/4) is 30 meV. For cobalt, the ground state is then found with 3 mu(B)/cell and the metastable state with 1 mu(B)/cell. The spin-crossover energy E (S) (1/3) is 21 meV. Our results suggest that these materials may be used in devices for spintronic applications that require different magnetization states.

Effect of hydrogen bond formation on the elastic molecular scattering: a case study with methanol

Rayleigh optical activities of small hydrogen-bonded methanol clusters containing two to five molecules are reported. For the methanol trimer, tetramer, and pentamer both cyclic and linear structures are considered. After the geometry optimizations, the dipole moments and the dipole polarizabilities (mean, interaction, and anisotropic components) are calculated using HF, MP2 and DFT (B3LYP, B3P86 and BH&HLYP) with aug-cc-pVDZ extended basis set. The polarizabilities are used to analyse the depolarization ratios and the Rayleigh scattering activities. The variations in the activity and in the depolarization for Rayleigh scattered radiation with the increase in the cluster size for both cyclic and linear structures are analysed.

Spectroscopic Study on the Structural Differences of Thermally Induced Cross-Linking Segments in Emeraldine Salt and Base Forms of Polyaniline

This paper reports the spectroscopic study on the structural differences of thermally induced cross-linking segments in polyaniline in its emeraldine salt (PANI-ES) and base (PANI-EB) forms. Casting films of PANI-ES (ES-film) and PANI-EB (EB-film) were prepared and heated at 150 degrees C under atmospheric air for 30 min. Raman spectra excited at 632.8 nm of heated ES-film presented the characteristic bands of phenazine-like structures at 1638, 1392, and 575 cm(-1), whereas EB-film showed lower relative intensities for these bands. The lower content of phenazine-like segments in heated EB-film is related to residual polaronic segments from preparation procedures, as revealed by Raman. This statement was confirmed by a sequence of thermal and doping experiments in both films. Quantum-chemical calculations by density functional theory (DFT) and time-dependent density functional theory (TD-DFT) showed that the phenazine-like structure presents the intense Raman band at 1350 cm(-1) due to heterocycle breathing mode, and the non-phenazine-like structure (substituted hydrophenazine-type) presents higher energy for HOMO-LUMO transition, indicating the lack of conjugation in the heterocycle compared with the phenazine-like structure. According to experimental and theoretical data reported here, it is proposed that only thermally treated PANI-ES presents phenazine-like rings, whereas PANI-EB presents heterocyclic non-aromatic structures.

Theoretical study of the XP3 (X = Al, B, Ga) clusters

The lowest singlet and triplet states of AlP3, GaP3 and BP3 molecules with C-s, C-2v and C-3v symmetries were characterized using the B3LYP functional and the aug-cc-pVTZ and aug-cc-pVQZ correlated consistent basis sets. Geometrical parameters and vibrational frequencies were calculated and compared to existent experimental and theoretical data. Relative energies were obtained with single point CCSD(T) calculations using the aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets, and then extrapolating to the complete basis set (CBS) limit. (C) 2011 Elsevier B.V. All rights reserved.

High Viscosity of Imidazolium Ionic Liquids with the Hydrogen Sulfate Anion: A Raman Spectroscopy Study

Ionic liquids based on 1-alkyl-3-methylimidazolium cations and the hydrogen sulfate (or bisulfate) anion, HSO4-, are much more viscous than ionic liquids with alkyl sulfates, RSO4-. The structural origin of the high viscosity of HSO4- ionic liquids is unraveled from detailed comparison of the anion Raman bands in 1-ethyl-3-methylimidazolium hydrogen sulfate and 1-butyl-3-methylimidazolium hydrogen sulfate with available data for simple HSO(4)(-) salts in crystalline phase, molten phase, and aqueous solution. Two Raman bands at 1046 and 1010 cm(-1) have been assigned as symmetric stretching modes nu(s)(S = O) of HSO4-, the latter being characteristic of chains of hydrogen-bonded anions. The intensity of this component increases in the supercooled liquid phase. For comparison purposes, Raman spectra of 1-ethyl-3-methylimidazolium ethyl sulfate and 1-butyl-3-methylimidazolium methyl sulfate have been also obtained. There is no indication of difference in the strength of hydrogen bond interactions of imidazolium cations with HSO4- or RSO4- anions. Raman spectra at high pressures, up to 2.6 GPa, are also discussed. Raman spectroscopy provides evidence that hydrogen-bonded anions resulting in anion-anion interaction is the reason for the high viscosity of imidazolium ionic liquids with HSO4-. If the ionic liquid is exposed to moisture, these structures are disrupted upon absorption of water from the atmosphere.

A theoretical study of SnF2+, SnCl2+, and SnO2+ and their experimental search

We present a detailed theoretical study of the stability of the gas-phase diatomic dications SnF2+, SnCl2+, and SnO2+ using ab initio computer calculations. The ground states of SnF2+, SnCl2+, and SnO2+ are thermodynamically stable, respectively, with dissociation energies of 0.45, 0.30, and 0.42 eV. Whereas SnF2+ dissociates into Sn2+ + F, the long range behaviour of the potential energy curves of SnCl2+ and SnO2+ is repulsive and wide barrier heights due to avoided crossing act as a kind of effective dissociation energy. Their equilibrium internuclear distances are 4.855, 5.201, and 4.852 a(0), respectively. The double ionisation energies (T-e) to form SnF2+, SnCl2+, and SnO2+ from their respective neutral parents are 25.87, 23.71, and 25.97 eV. We combine our theoretical work with the experimental results of a search for these doubly positively charged diatomic molecules in the gas phase. SnO2+ and SnF2+ have been observed for prolonged oxygen (O-16(-)) ion beam sputtering of a tin metal foil and of tin (II) fluoride (SnF2) powder, respectively, for ion flight times of about 10(-5) s through a magnetic-sector mass spectrometer. In addition, SnCl2+ has been detected for O-16(-) ion surface bombardment of stannous (tin (II)) chloride (SnCl2) powder. To our knowledge, SnF2+ is a novel gas-phase molecule, whereas SnCl2+ had been detected previously by electron-impact ionization mass spectrometry, and SnO2+ had been observed before by spark source mass spectrometry as well as by atom probe mass spectrometry. We are not aware of any previous theoretical studies of these molecular systems. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4758475]

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).