Glauco F. Bauerfeldt

Ozonolysis Reactions of Monoterpenes: A Variational Transition State Investigation.

The O3-initiated oxidation reactions of α-pinene ([1S,5S]-2,6,6-trimethylbicyclo[3.1.1]hept-2-ene), β-pinene ([1R,5R]-6,6-dimethyl-2-methylenebicyclo[3.1.1]heptane), camphene ([1R,4S]-2,2-dimethyl-3-methylenebicyclo[2.2.1]heptane) and sabinene ([1R,5R]-4-methylene-1-(1-methylethyl)bicycle[3.1.0]hexane), four monoterpenes typically emitted into the atmosphere, were studied at the B3LYP/6-31+G(2d,2p) level of theory. The rate coefficients were calculated on the basis of the variational transition state theory for two kinetic models, in order to investigate the reaction mechanism: first assuming a direct bimolecular reaction and the second, assuming the formation of a prebarrier-complex, which further reacts forming the corresponding molozonide. The barrier heights leading to the formation of exo-conformers of the molozonides of α-pinene and camphene are lower than the barrier heights for the formation of the endo-conformers of these molozonides, whereas the inverse trend is observed for β-pinene and sabinene. The canonical variational rate coefficients are found in reasonable agreement with the experimental data, especially when the prebarrier complexes are considered. Microcanonical variational rate coefficients are also calculated, as a final validation test, being found in an expected agreement with the canonical rate coefficients. The best predictions for the rate coefficients at 298 K, based on the microcanonical variational method, for α-pinene, β-pinene, camphene, and sabine are (in units cm(3) molecule(-1) s(-1)): 6.92 × 10(-17), 1.06 × 10(-17), 4.61 × 10(-19), and 4.81 × 10(-17), respectively. Our results suggest that the prebarrier complex is an important specie in the ozone addition mechanism and should be taken into account for the proper description of the overall kinetics.

Experimental and theorical study of the air quality in a suburban industrial-residential area in Rio de Janeiro, Brazil

Neste trabalho, dados de poluentes legislados e nao legislados foram obtidos em uma regiao industrial-residencial no Rio de Janeiro, no periodo de abril a novembro de 2002. As concentracoes de formaldeido e acetaldeido encontradas estavam na faixa de 25,25-231,73 ppb e 9,48-75,74 ppb respectivamente. Estas altas concentracoes podem ser atribuidas a diversas fontes: emissoes veiculares e industriais, a presenca de incineradores de lixo, efluentes e rejeitos despejados no rio Faria Timbo e emissoes de uma industria de metanol. As concentracoes medias de BTEX foram: 1,92, 9,02, 6,37, 21,65 e 5,99 µg m-3 para benzeno, tolueno, etilbenzeno, m,p-xileno e o-xileno, respectivamente. As concentracoes medias de PM10 mostram valores maximos para o mes de outubro, abril e junho. As concentracoes obtidas nas simulacoes foram para o mes de agosto usando um modelo de trajetorias. A concentracao maxima de ozonio calculada foi de 39,8 ppb, as 15:00 horas, obtendo uma boa concordância com os dados experimentais (38,5 ppb). A isopleta de ozonio mostra que o local esta em uma regiao alta de VOC/NOx, tipica de areas suburbanas e a favor do vento.

Thermochemical analysis and kinetics aspects for a chemical model for camphene ozonolysis.

In this work, a chemical model for the camphene ozonolysis, leading to carbonyl final products, is proposed and discussed on the basis of the thermochemical properties and kinetic data obtained at density functional theory levels of calculation. The mechanism is initiated by the electrophilic attack of ozone to the double bond in camphene leading to a 1,2,3-trioxolane intermediate, which decomposes to peroxy radicals and carbonyl compounds in a total of 10 elementary reactions. The thermodynamic properties (enthalpy and entropies differences) are calculated at 298 K. For the thermochemical evaluation, theoretical calculations are performed with the B3LYP, MPW1PW91, and mPW1K density functionals and the basis sets 6-31G(d), 6-31G(2d,2p), 6-31+G(d,p), and 6-31+G(2d,2p). Eventually, single point calculations adopting the 6-311++G(2d,2p) basis set are performed in order to improve the electronic energies. The enthalpy profiles suggest highly exothermic reactions for the individual steps, with a global enthalpy difference of -179.18 kcal mol(-1), determined at the B3LYP∕6-31+G(2d,2p) level. The Gibbs free energy differences for each step, at 298 K, calculated at the B3LYP∕6-311++G(2d,2p)∕∕B3LYP∕6-31+G(2d,2p) level, are used to estimate the composition of a final product mixture under equilibrium conditions as 58% of camphenilone and 42% of 6,6-dimethyl-ɛ-caprolactone-2,5-methylene. For the reaction kinetics, the bimolecular O(3) + camphene step is assumed to be rate determining in the global mechanism. A saddle point for the ozone addition to the double bond is located and rate constants are determined on the basis of the transition state theory. This saddle point is well represented by a loosely bound structure and corrections for the basis set superposition error (BSSE) are calculated, either by considering the effect over the geometry optimization procedure (here referred as CP1 procedure), or the effect of the BSSE over the electronic energy of a previously optimized geometry, included a posteriori (here referred as CP2). The rate constants, calculated at 298 K from the data obtained at the mPW1K∕6-31+G(d,p), CP1∕B3LYP∕∕6-31+G(2d,2p), and CP2∕B3LYP∕∕6-31+G(2d,2p) levels (3.62 × 10(-18), 1.12 × 10(-18), and 1.39 × 10(-18) cm(3) molecule(-1) s(-1)), are found in good agreement with the available experimental data at the same temperature, 0.9 × 10(-18) cm(3) molecule(-1) s(-1) [R. Atkinson, S. M. Aschmann, and J. Arey, Atmos. Environ. 24, 2647 (1990)]. The importance of the BSSE corrections for the final rate constants must be pointed out. Furthermore, this work will contribute to a better understanding of the chemistry of monoterpenes in the atmosphere, as well as the implications for the phenomena of pollution.

Study of the Stepwise Deprotonation Reactions of Glyphosate and the Corresponding pKa Values in Aqueous Solution

Glyphosate (N-(phosphonomethyl)glycine) (Gph) is a herbicide that is broadly used in several countries. Its application to eliminate weeds may have the undesired effect of diminishing the metallic cations found in the soil (e.g., Ni(2+) and Zn(2+)), due to a complexation reaction that depends on the soils pH. To better understand the molecular structures of glyphosate that are involved in such a complexation reaction, we have studied all possible glyphosate conformations in aqueous solution that may be involved in deprotonation reactions in the pH range from 2 to 11 using the polarizable continuum method (PCM). We have also compared direct (or absolute) methods to calculate pKa values, the cluster-continuum model and the proton-exchange scheme, using different thermodynamic cycles. The best result was achieved when using a proton-exchange scheme, which was able to properly reproduce three glyphosate experimental pKa values predicted for the glyphosate structures and conformations previously determined.

Theoretical calculations of the kinetics of the OH reaction with 2-methyl-2-propen-1-ol and its alkene analogue

In this work, the first and rate determining steps of the mechanism of the OH addition to 2-methyl-2-propen-1-ol (MPO221) and methylpropene (M2) have been studied at the DFT level, employing the BH and HLYP functional and the cc-pVDZ and aug-cc-pVDZ basis sets. The thermochemical properties of equilibrium (enthalpy, entropy and Gibbs free energies) have been determined within the conventional statistical thermodynamics relations and the rate coefficients have been determined on the basis of the variational transition state theory. The adoption of the microcanonical variational transition state theory was proved to be crucial for the description of the kinetics of OH addition to these unsaturated compounds. The rate coefficients obtained for the OH reactions with MPO221 and M2 at 298.15 K deviate, respectively, 27% and 13% from the experimental rate coefficient available in the literature. A non-Arrhenius profile is observed for the rate coefficients. Moreover, the values of the rate coefficients for the MPO221 + OH reaction are greater than those for the M2 + OH reaction, suggesting that the substitution of the hydrogen atom in an alkene by the –OH functional group increases the reactivity with respect to the hydroxyl radical.

Comparative kinetics of the 3-buten-1-ol and 1-butene reactions with OH radicals: a density functional theory/RRKM investigation.

The compared kinetics of the reactions of unsaturated alcohols and alkenes with OH radicals is a topic of great interest from both the theoretical chemistry and the atmospheric chemistry points of view. The enhanced reactivity of an unsaturated alcohol, with respect to its alkene analogue, toward OH radicals has been previously demonstrated, at 298 K, by experimental and theoretical research. In this work, a new comparative investigation of such reactions is performed for 3-buten-1-ol and 1-butene. The model assumes that the overall kinetics is governed by the first OH addition steps of the mechanism. Calculations have been performed at the DFT level, employing the BHandHLYP functional and the cc-pVDZ and aug-cc-pVDZ basis sets, and the rate coefficients have been determined on the basis of the microcanonical variational transition state theory. The rate coefficients obtained for the OH reactions with 3-buten-1-ol (kOH(31BO)) and 1-butene (kOH(1B)) at 298.15 K are lower than the experimental rate coefficient available in the literature, showing deviations of 18% and 25%, respectively. Negative temperature dependence is verified for these rate coefficients. The kOH(31BO)/kOH(1B) ratios have also been investigated as a function of the temperature, suggesting that at room temperature the unsaturated alcohol reacts with the OH radicals faster than 1-butene, by a factor of 1.2, but at higher temperatures (400-500 K), the alkene should react faster, and that the stabilization of prebarrier complexes and saddle points due to hydrogen bonds is no longer an important factor to govern the reactivity of the unsaturated alcohol toward OH radicals, with respect to the alkene analogue.

Ab initio and CVTST investigation of the gas phase nucleophilic substitution CH3Cl + OH- → CH3OH + Cl-

In this work the theoretical study of the gas-phase bimolecular nucleophilic substitution reaction, CH3Cl + OH-→ CH3OH + Cl-, is introduced aiming the description of the reaction path and the calculation of rate coefficients with the canonical variational transition state (CVTST) method. The calculations were performed at the MP2/6-31+G(d) level. The calculated enthalpy difference for the reaction at 298.15 K (-49.93 kcal mol-1) is in good agreement with the literature value: -50.4 kcal mol-1. The calculated rate coefficient, 1.94 × 10-9 cm3 molecule-1 s-1 at 298.15 K, also shows good agreement with the experimental data: 1.3-1.6 × 10-9 cm3 molecule-1 s-1. Moreover, the rate coefficients show non-Arrhenius behavior, decreasing as the temperature increases, which is consistent with the experimental expectation. In this way, the performance of the variational transition state theory for this reaction can be considered satisfactory.

Biodiesel Synthesis: Influence of Alkaline Catalysts in Methanol-Oil Dispersion

Biodiesel synthesis from soybean oil using methanol and alkaline catalysts occurs in the following two consecutive steps: dispersion of methanol in the oil and methanolysis. The effect of the alkaline catalysts NaOCH3, KOCH3, NaOH, and KOH in the dispersion step at 30-60 °C and under mechanical stirring at 400 rpm was evaluated. The dispersion step accounts for 44.6-73.3% of the total synthesis time and was poorly favored compared to methanolysis due to the increase in temperature. The catalysts decreased the dispersion time, although most of them increased the methanol-oil interfacial tension. K-containing catalysts were more active than their Na analogues due to higher adsorption of K in the methanol-oil interface and the higher production of methyl esters (which act as emulsifying agents), which promote a more favorable interfacial tension. The alkaline cation effect was more significant in the dispersion step than in the methanolysis step.

Determination of 4-Nonylphenol in Surface Waters of the Guandu River Basin by High Performance Liquid Chromatography with Ultraviolet Detection

4-Nonylphenol, a degradation product of ethoxylated alkylphenols, due to its harmful effects on the environment, has been banned in European Union countries, alongside their precursors. The guide on quality of drinking water from the United States Environmental Protection Agency (US EPA) recommends a maximum concentration of 28 μg L for fresh water. In Brazil, there is no clear legislation containing values of maximum concentration of 4-nonylphenol. Due to this lack of regulation, a continuous monitoring is necessary for this pollutant in environmental samples. The occurrence of 4-nonylphenol (4-NP) in the surface waters of Guandu River in the state of Rio de Janeiro, Brazil, was studied by using solid-phase extraction and reversed phase liquid chromatography separation with UV detection. The analytical method satisfies these requirements, being able to detect and quantify 4-NP in a desired concentration range. Of the 19 samples analyzed, 4-nonylphenol was detected in 12, quantified in 2, showing concentration levels of 1.73 and 2.32 μg L in Santa Cruz and Paracambi, respectively. This is the main hydrographic basin in the state of Rio de Janeiro, where water is collected for treatment and later distributed to most cities in the metropolitan region, including Rio de Janeiro City, and these results are therefore alarming.

Rate Coefficient for the Reaction of Cl Atoms with cis-3-Hexene at 296 ± 2 K

The rate coefficient of the cis-3-hexene + Cl atoms reaction at 296 ± 2 K and 750 ± 10 Torr was determined using the relative rate technique. The reaction was investigated using an 80 L Teflon reaction bag and a gas chromatograph coupled with flame-ionization detection. Chlorine atoms were produced by the photolysis of trichloroacetyl chloride. No previous experimental data was available in the literature, to the best of our knowledge. The mean second-order rate coefficient value found was (4.13 ± 0.51) × 10 cm molecule s. The experimental value agrees with the rate coefficient estimated by structure-reactivity analysis, 4.27 × 10 cm molecule s. Moreover, both addition and hydrogen abstraction channels contribute to the global kinetics, with branching ratios 70:30. Effective lifetime with respect to Cl atoms is predicted as 67.2 hours; however, the cis-3-hexene + Cl channel is suggested to be non-negligible at atmospheric conditions. Other atmospheric implications are discussed.