Accurate computations of the structures and binding energies of the imidazole ... benzene and pyrrole ... benzene complexes

Using explicitly-correlated coupled-cluster theory with single and double excitations, the intermolecular distances and interaction energies of the T-shaped imidazole⋯⋯benzene and pyrrole⋯⋯benzene complexes have been computed in a large augmented correlation-consistent quadruple-zeta basis set, adding also corrections for connected triple excitations and remaining basis-set-superposition errors. The results of these computations are used to assess other methods such as Møller–Plesset perturbation theory (MP2), spin-component-scaled MP2 theory, dispersion-weighted MP2 theory, interference-corrected explicitly-correlated MP2 theory, dispersion-corrected double-hybrid density-functional theory (DFT), DFT-based symmetry-adapted perturbation theory, the random-phase approximation, explicitly-correlated ring-coupled-cluster-doubles theory, and double-hybrid DFT with a correlation energy computed in the random-phase approximation.

Modeling the Histidine–Phenylalanine Interaction: The NH···π Hydrogen Bond of Imidazole·Benzene

NH···π hydrogen bonds occur frequently between the amino acid side groups in proteins and peptides. Data-mining studies of protein crystals find that ~80% of the T-shaped histidine···aromatic contacts are CH···π, and only ~20% are NH···π interactions. We investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet-cooled imidazole·benzene (Im·Bz) complex as a model for the NH···π interaction between histidine and phenylalanine. Ground- and excited-state dispersion-corrected density functional calculations and correlated methods (SCS-MP2 and SCS-CC2) predict that Im·Bz has a Cs-symmetric T-shaped minimum-energy structure with an NH···π hydrogen bond to the Bz ring; the NH bond is tilted 12° away from the Bz C₆ axis. IR depletion spectra support the T-shaped geometry: The NH stretch vibrational fundamental is red shifted by −73 cm⁻¹ relative to that of bare imidazole at 3518 cm⁻¹, indicating a moderately strong NH···π interaction. While the Sₒ(A1g) → S₁(B₂u) origin of benzene at 38 086 cm⁻¹ is forbidden in the gas phase, Im·Bz exhibits a moderately intense Sₒ → S₁ origin, which appears via the D₆h → Cs symmetry lowering of Bz by its interaction with imidazole. The NH···π ground-state hydrogen bond is strong, De=22.7 kJ/mol (1899 cm⁻¹). The combination of gas-phase UV and IR spectra confirms the theoretical predictions that the optimum Im·Bz geometry is T shaped and NH···π hydrogen bonded. We find no experimental evidence for a CH···π hydrogen-bonded ground-state isomer of Im·Bz. The optimum NH···π geometry of the Im·Bz complex is very different from the majority of the histidine·aromatic contact geometries found in protein database analyses, implying that the CH···π contacts observed in these searches do not arise from favorable binding interactions but merely from protein side-chain folding and crystal-packing constraints. The UV and IR spectra of the imidazole·(benzene)₂ cluster are observed via fragmentation into the Im·Bz+ mass channel. The spectra of Im·Bz and Im·Bz₂ are cleanly separable by IR hole burning. The UV spectrum of Im·Bz₂ exhibits two 000 bands corresponding to the Sₒ → S₁ excitations of the two inequivalent benzenes, which are symmetrically shifted by −86/+88 cm⁻¹ relative to the 000 band of benzene.

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Attempted hydrogen–deuterium exchange of the protio-trimethyloxonium dication (CH3)3OH2+, study of methylating ability of (CH3)3O+ in superacids and theoretical investigations

Attempted hydrogen–deuterium exchange of trimethyloxonium ion, (CH3)3O+ with excess of 1:1 2HF/SbF5 superacid at −30°C over a period of 30 days showed no exchange. Theoretical calculations at the MP2/6–31G** level are in accord with the lack of hydrogen–deuterium exchange in the methyl group of the (CH3)3O+ cation as protonation (protosolvation) prefers the oxygen lone pair of electrons, instead of a C—H bond. Methylation of aromatics with the (CH3)3O+CF3SO3− in CF3SO3H and 2CF3SO3H:B(O3SCF3)3 was also studied. Whereas in triflic acid no alkylation was observed, in triflatoboric acid, a powerful superacid, alkylation takes place, indicating protolytic activation of the trimethyloxonium ion.

Calculated 11B–13C NMR chemical shift relationship in hypercoordinate methonium and boronium ions

The boronium-carbonium continuum was extended to include hypercoordinated protonated methanes and their boron analogs. The 11B NMR chemical shifts of the hypercoordinated hydriodo boron compounds and the 13C NMR chemical shifts of the corresponding isoelectronic and isostructural carbocations were calculated by using the GIAO-MP2 method. The data show good linear correlation between 11B and 13C NMR chemical shifts, which indicates that the same factors that determine the chemical shifts of the boron nuclei also govern the chemical shifts of carbon nuclei of these hypercoordinated hydriodo compounds.

Desenvolvimento e avaliação de métodos de extração e separação cromatográfica em colunas monolíticas e superficialmente porosas para determinação de herbicidas triazínicos em solos e águas

O uso de pesticidas levou ao aumento da produtividade e qualidade dos produtos agrícolas, porém o seu uso acarreta na intoxicação dos seres vivos pela ingestão gradativa de seus resíduos que contaminam o solo, a água e os alimentos. Dessa forma, há a necessidade do monitoramento constante de suas concentrações nos compartimentos ambientais. Para isto, busca-se o desenvolvimento de métodos de extração e enriquecimento de forma rápida, com baixo custo, gerando um baixo volume de resíduos, contribuindo com a química verde. Dentre estes métodos destacam-se a extração por banho de ultrassom e a extração por ponto nuvem. Após o procedimento de extração, o extrato obtido pode ser analisado por técnicas de Cromatografia a Líquido de Alta Eficiência (HPLC) e a Cromatografia por Injeção Sequencial (SIC), empregando fases estacionárias modernas, tais como as monolíticas e as partículas superficialmente porosas. O emprego de SIC com coluna monolítica (C18, 50 x 4,6 mm) e empacotada com partículas superficialmente porosas (C18, 30 x 4,6 mm, tamanho de partícula 2,7 µm) foi estudado para separação de simazina (SIM) e atrazina (ATR), e seus metabólitos, desetilatrazina (DEA), desisopropilatrazina (DIA) e hidroxiatrazina (HAT). A separação foi obtida por eluição passo-a-passo, com fases móveis compostas de acetonitrila (ACN) e tampão Acetato de Amônio/Ácido acético (NH4Ac/HAc) 2,5 mM pH 4,2. A separação na coluna monolítica foi realizada com duas fases móveis: MP1= 15:85 (v v-1) ACN:NH4Ac/HAc e MP2= 35:65 (v v-1) ACN:NH4Ac/HAc a uma vazão de 35 µL s-1. A separação na coluna com partículas superficialmente porosas foi efetivada com as fases móveis MP1= 13:87 (v v-1) ACN: NH4Ac/HAc e MP2= 35:65 (v v-1) ACN:NH4Ac/HAc à vazão de 8 µL s-1. A extração por banho de ultrassom em solo fortificado com os herbicidas (100 e 1000 µg kg-1) resultou em recuperações entre 42 e 160%. A separação de DEA, DIA, HAT, SIM e ATR empregando HPLC foi obtida por um gradiente linear de 13 a 35% para a coluna monolítica e de 10 a 35% ACN na coluna com partículas superficialmente porosas, sendo a fase aquosa constituída por tampão NH4Ac/HAc 2,5 mM pH 4,2. Em ambas as colunas a vazão foi de 1,5 mL min-1 e o tempo de análise 15 min. A extração por banho de ultrassom das amostras de solo com presença de ATR, fortificadas com concentrações de 250 a 1000 µg kg-1, proporcionou recuperações entre 40 e 86%. A presença de ATR foi confirmada por espectrometria de massas. Foram realizados estudos de fortificação com ATR e SIM em amostras de água empregando a extração por ponto nuvem com o surfactante Triton-X114. A separação empregando HPLC foi obtida por um gradiente linear de 13 a 90% de ACN para a coluna monolítica e de 10 a 90% de ACN para a coluna empacotada, sempre em tampão NH4Ac/HAc 2,5 mM pH 4,2. Em ambas as colunas a vazão foi de 1,5 mL min-1 e o tempo de análise 16 min. Fortificações entre 1 e 50 µg L-1 resultaram em recuperações entre 65 e 132%.

Marine palynomorph abundance estimates of sediment core CRP-2A (Fig. 2)

A diverse assemblage of marine palynomorphs was recovered from the Oligocene - Miocene section of CRP-2/2A. Most of the assemblage is composed of previously unrecognised species. Three distinct groups of marine palynomorph were recognised: (1) prasinophytes, mainly Cymatiosphaera, (2) acritarchs, mainly Leiosphaeridia and Sigmopollis although Leiofusa is an important component of the bottom half of the hole, and (3) dinoflagellate cysts. About 27 species of in situ dinoflagellate cysts were recorded, of which seven apparently undescribed species of Lejeunecysta form a prominent component. Reworked specimens of several species of the Paleogene Transantarctic Flora occur in CRP-2/2A sediments. Several abundance peaks of reworked taxa from the Transantarctic Flora are recorded. Three marine palynomorph zones were recognised (MP3, MP2, MP1), considered to be early Oligocene, late Oligocene, and late Oligocene/early Miocene in age respectively. Samples from the Quaternary and Pliocene part of CRP-2/2A were also examined. These proved either barren or yielded very sparse low diversity floras.

Iminopropadienones R-N=C=C=C=O and carbon suboxide, C3O2. Theoretical and experimental 13C NMR spectra

The C-13 NMR data of five iminopropadienones R-N=C=C=C=O as well as carbon suboxide, C3O2, have been examined theoretically and experimentally. The best theoretical results were obtained using the GIAO/B3LYP/6-31 +G**//MP2/6-31G* level of theory, which reproduces the chemical shifts of the iminopropadienone substituents extremely well while underestimating those of the cumulenic carbons by 5-10 ppm. The computationally faster GIAO/HF/6-31 + G**//B3LYP/6-31 G* level is also adequate. (C) 2004 Elsevier B.V. All rights reserved.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Formation of Polycyclic Aromatic Hydrocarbons and Nitrogen Containing Polycyclic Aromatic Compounds in Titan's Atmosphere, the Interstellar Medium and Combustion

Several different mechanisms leading to the formation of (substituted) naphthalene and azanaphthalenes were examined using theoretical quantum chemical calculations. As a result, a series of novel synthetic routes to Polycyclic Aromatic Hydrocarbons (PAHs) and Nitrogen Containing Polycyclic Aromatic Compounds (N-PACs) have been proposed. On Earth, these aromatic compounds originate from incomplete combustion and are released into our environment, where they are known to be major pollutants, often with carcinogenic properties. In the atmosphere of a Saturn's moon Titan, these PAH and N-PACs are believed to play a critical role in organic haze formation, as well as acting as chemical precursors to biologically relevant molecules. The theoretical calculations were performed by employing the ab initio G3(MP2,CC)/B3LYP/6-311G** method to effectively probe the Potential Energy Surfaces (PES) relevant to the PAH and N-PAC formation. Following the construction of the PES, Rice-Ramsperger-Kassel-Markus (RRKM) theory was used to evaluate all unimolecular rate constants as a function of collision energy under single-collision conditions. Branching ratios were then evaluated by solving phenomenological rate expressions for the various product concentrations. The most viable pathways to PAH and N-PAC formation were found to be those where the initial attack by the ethynyl (C2H) or cyano (CN) radical toward a unsaturated hydrocarbon molecule led to the formation of an intermediate which could not effectively lose a hydrogen atom. It is not until ring cyclization has occurred, that hydrogen elimination leads to a closed shell product. By quenching the possibility of the initial hydrogen atom elimination, one of the most competitive processes preventing the PAH or N-PAC formation was avoided, and the PAH or N-PAC formation was allowed to proceed. It is concluded that these considerations should be taken into account when attempting to explore any other potential routes towards aromatic compounds in cold environments, such as on Titan or in the interstellar medium.

Determinação da concentração de metais no material particulado (PM10 e PM2,5) do ar atmosférico em uma área urbana da cidade do Natal – RN

This study aimed to assess ambient air quality in a urban area of Natal, capital of Rio Grande do Norte (latitude 5º49'29 '' S and longitude 35º13'34'' W), aiming to determine the metals concentration in particulate matter (PM10 and PM2,5) of atmospheric air in the urban area o the Natal city. The sampling period for the study consisted of data acquisition from January to December 2012. Samples were collected on glass fiber filters by means of two large volumes samplers, one for PM2,5 (AGV PM 2,5) and another for PM10 (PM10 AGV). Monthly averages ranged from 8.92 to 19.80 g.m-3 , where the annual average was 16,21 g.m-3 for PM10 and PM2,5 monthly averages ranged from 2,84 to 7,89 g.m -3 , with an annual average of 5,61 g.m-3 . The results of PM2,5 and PM10 concentrations were related meteorological variables and for information on the effects of these variables on the concentration of PM, an exploratory analysis of the data using Principal Component Analysis (PCA) was performed. The results of the PCA showed that with increasing barometric pressure, the direction of the winds, the rainfall and relative humidity decreases the concentration of PM and the variable weekday little influence compared the meteorological variables. Filters containing particulate matter were selected in six days and subjected to microwave digestion. After digestion samples were analyzed by with Inductively Coupled Plasma Mass Spectrometry (ICP-MS). The concentrations for heavy metals Vanadium, Chromium, Manganese, Nickel, Copper, Arsenic and lead were determined. The highest concentrations of metals were for Pb and Cu, whose average PM10 values were, respectively, 5,34 and 2,34 ng.m-3 and PM2,5 4,68 and 2,95 ng.m-3 . Concentrations for metals V, Cr, Mn, Ni, and Cd were respectively 0,13, 0,39, 0,48, 0,45 and 0,03 ng.m-3 for PM10 fraction and PM2,5 fraction, 0,05, 0,10, 0,10, 0,34 and 0,01 ng.m-3. The concentration for As was null for the two fractions

Estudo teórico de intermediários tetraédricos acidez / basicidade e estereosseletividade enzimáticos

The present work aimed first, the theoretical study of tetrahedral intermediate stability formed from carbonyl addition reactions using the second (MP2) and third (MP3) order Møller–Plesset perturbation theory. Linear correlations between electronic energy difference of reactions with Wiberg Indexes and C-O bond lengths were obtained, and was observed that the stability of adducts formed depends directly of electronic density involved between these atoms. The knowing of electronic parameters of these structures has an important hole due to the large use on reactions that in his course forms this tetrahedral intermediate. Employing the ONIOM (B3LYP:AMBER) methodology, was evaluated the stereoselectivity of a enzymatic reaction between CAL B enzyme and a long chain ester. In this study, were obtained the electronic energies of ground state and intermediate state of transesterification rate-determing step from two possible proquirals faces Re and Si. The objective was study the enantioselectivity of CAL B and rationalizes it using quantum theory of atoms in molecules (QTAIM). A theoretical study employing inorganic compounds was performed using ab initio CBS-QB3 method aiming to find a link between thermodynamic and equilibrium involving acids and bases. The results observed showed an excellent relationship between difference in Gibbs free energy, ΔG of acid dissociation reaction and ΔG of hydrolysis reaction of the corresponding conjugate base. It was also observed, a relationship between ΔG of hydrolysis reaction of conjugate acids and their corresponding atomic radius showing that stability plays an important role in hydrolysis reactions. The importance of solvation in acid/base behavior when compared to theoretical and experimental ΔG´s also was evaluated.

Generalized atoms in molecules approach (GAIM approach)

This thesis involves two parts. The first is a new-proposed theoretical approach called generalized atoms in molecules (GAIM). The second is a computational study on the deamination reaction of adenine with OH⁻/nH₂O (n=0, 1, 2, 3) and 3H₂O. The GAIM approach aims to solve the energy of each atom variationally in the first step and then to build the energy of a molecule from each atom. Thus the energy of a diatomic molecule (A-B) is formulated as a sum of its atomic energies, EA and EB. Each of these atomic energies is expressed as, EA = Hᴬ + Vₑₑᴬᴬ + 1/2Vₑₑᴬ<>ᴮ EB = Hᴮ + Vₑₑᴮᴮ + 1/2Vₑₑᴬ<>ᴮ where; Hᴬ and Hᴮ are the kinetic and nuclear attraction energy of electrons of atoms A and B, respectively; Vₑₑᴬᴬ and Vₑₑᴮᴮ are the interaction energy between the electrons on atoms A and B, respectively; and Vₑₑᴬ<>ᴮ is the interaction energy between the electrons of atom A with the electrons of atom B. The energy of the molecule is then minimized subject to the following constraint, |ρA(r)dr + |ρB(r)dr = N where ρA(r) and ρB(r) are the electron densities of atoms A and B, respectively, and N is the number of electrons. The initial testing of the performance of GAIM was done through calculating dissociation curves for H₂, LiH, Li₂, BH, HF, HCl, N₂, F₂, and Cl₂. The numerical results show that GAIM performs very well with H₂, LiH, Li₂, BH, HF, and HCl. GAIM shows convergence problems with N₂, F₂, and Cl₂ due to difficulties in reordering the degenerate atomic orbitals Pₓ, Py, and Pz in N, F, and Cl atoms. Further work for the development of GAIM is required. Deamination of adenine results in one of several forms of premutagenic lesions occurring in DNA. In this thesis, mechanisms for the deamination reaction of adenine with OH⁻/nH₂O, (n = 0, 1, 2, 3) and 3H₂O were investigated. HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), and B3LYP/6-31+G(d) levels of theory were employed to optimize all the geometries. Energies were calculated at the G3MP2B3 and CBS-QB3 levels of theory. The effect of solvent (water) was computed using the polarizable continuum model (PCM). Intrinsic reaction coordinate (IRC) calculations were performed for all transition states. Five pathways were investigated for the deamination reaction of adenine with OH⁻/nH₂O and 3H₂O. The first four pathways (A-D) begin with by deprotonation at the amino group of adenine by OH⁻, while pathway E is initiated by tautomerization of adenine. For all pathways, the next two steps involve the formation of a tetrahedral intermediate followed by dissociation to yield products via a 1,3-hydrogen shift. Deamination with a single OH⁻ has a high activation barrier (190 kJ mol⁻¹ using G3MP2B3 level) for the rate-determining step. Addition of one water molecule reduces this barrier by 68 kJ mol⁻¹ calculated at G3MP2B3 level. Adding more water molecules decreases the overall activation energy of the reaction, but the effect becomes smaller with each additional water molecule. The most plausible mechanism is pathway E, the deamination reaction of adenine with 3H₂O, which has an overall G3MP2B3 activation energy of 139 and 137 kJ mol⁻¹ in the gas phase and PCM, respectively. This barrier is lower than that for the deamination with OH⁻/3H₂O by 6 and 2 kJ mol⁻¹ in the gas phase and PCM, respectively.