Estudos do mecanismo de reações unimoleculares, das moléculas de acetato de metila e etila, induzidas por absorção multifotônica e de reações bimoleculares com os radicais OH e Cl, usando cálculos ab initio e de funcional de densidade foram realizados neste trabalho. A análise dos cálculos das barreiras de energia e das constantes de velocidade microcanônicas das reações unimoleculares com o formalismo da teoria RRKM, permitiu prever algumas reações cujos produtos não foram determinados experimentalmente. Além disso, semelhanças das geometrias e dos valores de energia de algumas estruturas entre diferentes superfícies de energia, analisadas com cálculos de IRC, permitiram a explicação da viabilidade de determinada via de reação em detrimento de outra. Em outro método, relatado por Forst(1), na determinação das constantes de velocidade variacionais, comparadas com a variação da energia de Gibbs, verificou-se constantes de velocidade microcanônicas menores e a localização dos estados de transição em comprimentos de ligação, também, menores. Nas reações com o radical OH e Cl, a via predominante de reação é o ataque aos hidrogênios do grupo ligado diretamente ao oxigênio da molécula, responsável pelo valor da constante final de mais de 93% nas reações com OH e 99% nas reações com Cl, a 298K, em concordância com os resultados experimentais, em detrimento do ataque ao grupo CH3 da porção carboxilato da molécula. Nas reações com OH, os valores de constante de velocidade calculados mais próximos dos valores experimentais foram obtidos a partir dos resultados com os funcionais mPW1B95-41, para o acetato de metila, e mPW1B95-44, para o acetato de etila, enquanto que nas reações com Cl foram obtidos com o os métodos CCSD(T)//B3LYP para o acetato de metila e CCSD(T)//MP2 para o acetato de etila. Os valores de constante de velocidade da reação com cloro são cerca de dez vezes maiores que aquelas para as reações com radicais OH. As reações com OH e Cl ocorrem em uma única etapa, sem estabilização do intermediário e sem efeito de tunelamento significante. Por sua vez, cálculos da constante de velocidade, pelo método da relação estrutura reatividade (SAR), sobreestimam a reatividade dos hidrogênios dos grupos CH3 na porção carboxilato, em ambos os acetatos, e na porção etóxido do acetato de etila. Este método falha ao descrever a participação de cada grupo na reação com Cl, pois prevê que a reação no grupo CH3 da porção alcóxido passa a ser predominante sobre a reação ao grupo CH2 nas reações do acetato de etila, oposto aos nossos cálculos onde a reação com o grupo CH2 é a predominante. Outra falha do método é na previsão de mesma reatividade para ambos os grupos CH3 no acetato de metila, pois por nossos cálculos a reação ocorre predominantemente no grupo CH3 na porção alcóxido da molécula. / Computational studies on methyl and ethyl acetates molecules using ab initio and density functional calculations exploring the unimolecular mechanism, induced by multiphoton absorption, and the bimolecular reactions with OH and Cl radicals, have been performed in this work. Analysis of the calculated energy barriers and rate constants of unimolecular reactions with the RRKM microcanonical theory, predicts the occurrence of some reactions whose products were not determined experimentally. Furthermore, similarities on geometrical and energetic of some structures between different energy surfaces, analyzed with IRC calculations, allowed the explanation of the viability of a particular reaction pathway over another. Forst\'s method, used to determine variational rate constants, when compared with the variation of Gibbs energy, generates microcanonical rate constants with smaller values and location of transition states in smaller bond lengths. With OH and Cl, the predominant reaction route is the attack on the hydrogens of the group bonded directly to oxygen in the molecule, responsible for more than 93% of the final constant value in reactions with OH and 99% in reactions with Cl , at 298 K, in agreement with the experimental results, while the attack to the CH3 group in the carboxylate portion of the molecule is the less likely to occur. In reactions with OH, the values of the rate constant calculated closer to the experimental values were obtained from the results with functional mPW1B95-41, for methyl acetate, and mPW1B95-44 for ethyl acetate, whereas in reactions with Cl, were obtained with the CCSD(T)//B3LYP method for the methyl acetate and the CCSD(T)//MP2 method for ethyl acetate. The rate constant for the reaction with chlorine are about ten times larger than those for reactions with OH radicals. Reactions with OH and Cl occur in a single step, without stabilization of the intermediary and without significant tunneling effect. The rate constant obtained by the structure-reactivity relationship (SAR) overestimates the reactivity of the hydrogens of the CH3 groups at the carboxylate portion, in both acetates, and at the ethoxide portion of ethyl acetate. This method fails to describe the participation of each group in the reaction with Cl, predicting that the reaction on the CH3 portion of the alkoxide becomes predominant over the reaction on the CH2 group of ethyl acetate, relative to our calculations where the reaction with the CH2 group is predominant. Another flaw by providing the same reactivity for both the CH3 in methyl acetate, is in disagreement with our calculations which indicate that the reaction occurs predominantly in the CH3 group at the alkoxide portion of the molecule.

Estudos do mecanismo de reações unimoleculares, das moléculas de acetato de metila e etila, induzidas por absorção multifotônica e de reações bimoleculares com os radicais OH e Cl, usando cálculos ab initio e de funcional de densidade foram realizados neste trabalho. A análise dos cálculos das barreiras de energia e das constantes de velocidade microcanônicas das reações unimoleculares com o formalismo da teoria RRKM, permitiu prever algumas reações cujos produtos não foram determinados experimentalmente. Além disso, semelhanças das geometrias e dos valores de energia de algumas estruturas entre diferentes superfícies de energia, analisadas com cálculos de IRC, permitiram a explicação da viabilidade de determinada via de reação em detrimento de outra. Em outro método, relatado por Forst(1), na determinação das constantes de velocidade variacionais, comparadas com a variação da energia de Gibbs, verificou-se constantes de velocidade microcanônicas menores e a localização dos estados de transição em comprimentos de ligação, também, menores. Nas reações com o radical OH e Cl, a via predominante de reação é o ataque aos hidrogênios do grupo ligado diretamente ao oxigênio da molécula, responsável pelo valor da constante final de mais de 93% nas reações com OH e 99% nas reações com Cl, a 298K, em concordância com os resultados experimentais, em detrimento do ataque ao grupo CH3 da porção carboxilato da molécula. Nas reações com OH, os valores de constante de velocidade calculados mais próximos dos valores experimentais foram obtidos a partir dos resultados com os funcionais mPW1B95-41, para o acetato de metila, e mPW1B95-44, para o acetato de etila, enquanto que nas reações com Cl foram obtidos com o os métodos CCSD(T)//B3LYP para o acetato de metila e CCSD(T)//MP2 para o acetato de etila. Os valores de constante de velocidade da reação com cloro são cerca de dez vezes maiores que aquelas para as reações com radicais OH. As reações com OH e Cl ocorrem em uma única etapa, sem estabilização do intermediário e sem efeito de tunelamento significante. Por sua vez, cálculos da constante de velocidade, pelo método da relação estrutura reatividade (SAR), sobreestimam a reatividade dos hidrogênios dos grupos CH3 na porção carboxilato, em ambos os acetatos, e na porção etóxido do acetato de etila. Este método falha ao descrever a participação de cada grupo na reação com Cl, pois prevê que a reação no grupo CH3 da porção alcóxido passa a ser predominante sobre a reação ao grupo CH2 nas reações do acetato de etila, oposto aos nossos cálculos onde a reação com o grupo CH2 é a predominante. Outra falha do método é na previsão de mesma reatividade para ambos os grupos CH3 no acetato de metila, pois por nossos cálculos a reação ocorre predominantemente no grupo CH3 na porção alcóxido da molécula. / Computational studies on methyl and ethyl acetates molecules using ab initio and density functional calculations exploring the unimolecular mechanism, induced by multiphoton absorption, and the bimolecular reactions with OH and Cl radicals, have been performed in this work. Analysis of the calculated energy barriers and rate constants of unimolecular reactions with the RRKM microcanonical theory, predicts the occurrence of some reactions whose products were not determined experimentally. Furthermore, similarities on geometrical and energetic of some structures between different energy surfaces, analyzed with IRC calculations, allowed the explanation of the viability of a particular reaction pathway over another. Forst\'s method, used to determine variational rate constants, when compared with the variation of Gibbs energy, generates microcanonical rate constants with smaller values and location of transition states in smaller bond lengths. With OH and Cl, the predominant reaction route is the attack on the hydrogens of the group bonded directly to oxygen in the molecule, responsible for more than 93% of the final constant value in reactions with OH and 99% in reactions with Cl , at 298 K, in agreement with the experimental results, while the attack to the CH3 group in the carboxylate portion of the molecule is the less likely to occur. In reactions with OH, the values of the rate constant calculated closer to the experimental values were obtained from the results with functional mPW1B95-41, for methyl acetate, and mPW1B95-44 for ethyl acetate, whereas in reactions with Cl, were obtained with the CCSD(T)//B3LYP method for the methyl acetate and the CCSD(T)//MP2 method for ethyl acetate. The rate constant for the reaction with chlorine are about ten times larger than those for reactions with OH radicals. Reactions with OH and Cl occur in a single step, without stabilization of the intermediary and without significant tunneling effect. The rate constant obtained by the structure-reactivity relationship (SAR) overestimates the reactivity of the hydrogens of the CH3 groups at the carboxylate portion, in both acetates, and at the ethoxide portion of ethyl acetate. This method fails to describe the participation of each group in the reaction with Cl, predicting that the reaction on the CH3 portion of the alkoxide becomes predominant over the reaction on the CH2 group of ethyl acetate, relative to our calculations where the reaction with the CH2 group is predominant. Another flaw by providing the same reactivity for both the CH3 in methyl acetate, is in disagreement with our calculations which indicate that the reaction occurs predominantly in the CH3 group at the alkoxide portion of the molecule.

Nas últimas décadas, o desenvolvimento das técnicas de ionização à pressão atmosférica impulsionou a espectrometria de massas, na caracterização e elucidação estrutural de compostos de grande massa molecular. O surgimento dessas técnicas foi o responsável pela amplitude nas aplicações e estudos de espectrometria de massas, sendo a ionização por electrospray a mais versátil dentre essas fontes de ionização. O caráter eletrolítico da fonte de ionização por electrospray permite-se obter íons provenientes de três processos químicos: i) ácido-base; ii) redox e iii) complexação. A extensão com que cada um desses processos ocorrerá dependerá de fatores relacionados à operação da fonte de ionização e grandezas termoquímicas do analito. O notável progresso em técnicas experimentais, processamento de dados e integração entre as mais diversas áreas de aplicação da química, tem estimulado e beneficiado a aplicação da química teórica em estudos de reações em fase gasosa. A aplicação da química computacional fornece uma compreensão quantitativa das variações estruturais e energéticas dos possíveis íons formados durante a ionização da amostra, permitindo também a compreensão das possíveis vias de dissociação. É neste sentido, que o sinergismo entre a aplicação de conceitos derivados da química quântica pode auxiliar nas análises de espectrometria de massas. O objetivo desta tese foi o de se aplicar os modelos fundamentados na mecânica quântica para obtenção de grandezas termoquímicas relacionadas aos fenômenos que ocorrem durante as análises de espectrometria de massas por electrospray. Inicialmente, a comparação entre métodos ab initio, modelos compostos e aqueles embasados na teoria do funcional de densidade foram empregados nos cálculos das grandezas termoquímicas para a -butirolactona e 2-pirrolidinona, com a finalidade de se obter parâmetros termoquímicos de alta qualidade. Os modelos compostos G2, G2(MP2), CBS-Q, CBS-QB3 e os métodos B3LYP, B3P86, B98, PW91PW91 e MP2 foram testados. Os valores obtidos para a entalpia de formação, afinidade protônica e basicidade em fase gasosa para essas duas moléculas foram comparados aos dados experimentais disponíveis na literatura. Os melhores resultados para os valores de entalpia de formação foram obtidos ao se empregar o modelo B3LYP/6-31+G(d,p). A afinidade protônica e basicidade em fase gasosa foram mais bem descritas por B3LYP e G2(MP2). Posteriormente, foram estudadas a 1,4-benzoquinona, 1,4-naftoquinona, bem como seus derivados (2-acetilamina-1,4-naftoquinona; 2-propionilamina-1,4-naftoquinona; 2-butirilamina-1,4-naftoquinona, 2-benzoilamina-1,4-naftoquinona, 2-succinilamina-1,4-naftoquinona e lapachol). A escolha do modelo teórico empregado foi realizada por se comparar as geometrias, afinidade protônica, basicidade em fase gasosa, energia de ionização e afinidade eletrônica obtidos para a 1,4-benzoquinona com aqueles disponíveis na literatura. O modelo B3LYP/6-31+G(d,p) foi o mais exato com relação as grandezas termoquímicas supracitadas, assim, este modelo foi empregado no cálculo dessas mesmas grandezas para a 1,4-naftoquinona e seus derivados. A influência do substituinte na estrutura eletrônica dessas moléculas protonadas, desprotonadas, oxidadas, reduzidas e cationizadas foi estudada utilizando-se análises energéticas, geométricas, eletrônicas e topológicas. O desenvolvimento desses estudos compreendeu as análises das grandezas termoquímicas e análises da densidade eletrônica pelos métodos NBO (Natural Bond Orbital), NSA (Natural Steric Analysis), NRT (Natural Resonance Theory) e AIM (Atoms in molecules). A ionização por electrospray e a dissociação induzida por colisão foram realizadas para os derivados da 1,4-naftoquinona, sendo analisados suas moléculas protonadas, desprotonadas, reduzidas, oxidadas e cationizadas com Na+. As vias mecanísticas de dissociação foram embasadas nas análises da superfície de energia obtidas pelo cálculo das energias de Gibbs e entalpias. / In recent decades, the development of atmospheric ionization techniques improved mass spectrometry, principally for characterization and structural elucidation of high-molecular weight compounds. The development of spray ionization was responsible for the spread of applications and studies of mass spectrometry, where the electrospray ionization is the most versatile among the ionization sources. The electrolytic character of electrospray source allows obtaining ions by three different chemical processes: i) acid-base; ii) redox and, iii) metal complexation. These processes will occur through several factors which can be related to the ionization source and thermochemical parameters of analyte. The notable progress of experimental analyses, computational data and the integration between several areas of chemical application have stimulated the use of theoretical chemistry at gas-phase studies. Computational chemistry can furnish a quantitative understanding of the structure and energy of possible ions during the ionization process. For this reason, the synergism between the concepts from quantum chemistry and gas-phase chemistry can help mass spectrometry analysis. The main purpose of this thesis was the application of several quantum mechanical models to obtain thermochemical parameters which can be related to mass spectrometry phenomena. Firstly, the comparison between ab initio, composite model and DFT methods were employed to obtain the thermochemical parameters to -butyrolactone e 2-pyrrolidinona, in order to obtain high performance of thermochemical parameters. The composite G2, G2(MP2), CBS-Q and, the B3LYP, B3P86, B98, PW91PW91 and MP2 methods were tested. The calculated values were compared to experimental values reported in the literature. The best results for enthalpies of formation were obtained when B3LYP/6-31+G(d,p) model was employed. The proton affinity and gas-phase basicity were better described by using of B3LYP and G2(MP2). Secondly, the studies with quinones were performed, where the 1,4-benzoquinone, 1,4-naphthoquinone and its derivatives (2-acylamino-1,4-naphthoquinone; 2-propyonilamino-1,4-naphthoquinone; 2-butyrilamino-1,4-naphthoquinone; 2-benzoylamino-1,4-naphthoquinone; 2-succynilamino-1,4-naphtoquinone and, lapachol) were studied. A search for a theoretical model was made to compare the geometries, proton affinity, gas-phase basicity, ionization energy and electron affinity to 1,4-benzoquinone with those reported in the literature. The most accurate results were obtained by using of B3LYP/6-31+G(d,p). Thus, this model was applied in all studies with 1,4-naphthoquinone derivatives. The influence of substituent groups on electronic structure of protonated, deprotonated, reduced, oxidized and cationized molecules were studied by energetic, geometrics, electronics and topological analyses. The development of these studies and the determination of the thermochemical parameters and wave function analysis was achieved by means of Natural Bond Orbitals, Natural Steric Analysis, Natural Resonance Theory and Atoms in Molecules The electrospray ionization and gas-phase collision-induced dissociation were made for the 1,4-naphtoquinone derivatives by analyzing their protonated, deprotonated species, the radicalar and sodiated ones. The main fragmentation pathways were elucidated on the basis of the energy surface by using Gibbs energies and enthalpies.

Gas phase molecular clusters present an ideal medium for observing factors that drive chemical reactions without outside interferences from excessive solvent molecules. Introducing an ion into the cluster promotes ion-molecule interactions that may manifest in a variety of non-covalent or even covalent binding motifs and are of significant importance in many fields including atmospheric and astronomical sciences. For instance, in outer space, molecules are subject to ionizing radiation where ion-molecule reactions become increasingly competitive to molecule-molecule interactions. To elucidate individual ion-molecule interaction information, mass spectrometry was used in conjunction with appropriate theoretical calculations. Three main categories of experiment were conducted in this dissertation. The first of which were thermochemical equilibrium measurements where an ion was introduced to an ion mobility drift cell wherein thermalizing collisions occur with helium buffer gas facilitating a reversible reaction with a neutral molecule allowing the standard changes in enthalpy and entropy to be determined. The second type of experiment was an ion mobility experiment where an ionized homo- or hetero-cluster was injected into the drift cell at specific conditions allowing the reduced mobility and collisional cross-section to be evaluated. Thirdly, kinetics measurements were taken following injection of an ion into the drift cell were an irreversible reaction ensued with the neutral species hindering equilibrium, but prompting rate constant assessment. Previous research has laid the groundwork for this dissertation as the results and discussion contained herein will build upon existing data while maintaining originality. For example, past work has given support for ion-molecule reactions involving precursor species such as acetylene and hydrogen cyanide to form more complex organics, perhaps leading to biologically relevant species. The chemical systems studied for this research are either ionized substituted benzenes like phenylacetylene and benzonitrile or polycyclic aromatic nitrogen-containing hydrocarbons like quinoline and quinoxaline interacting with a variety of neutral species. Hydrogen bonding and its many sub-sections are of the utmost importance to the kinds of reactions studied here. Past work has shown the tendency of organic radical cations to form conventional and unconventional ionic hydrogen bonds with gas phase solvents. Other non-covalent modes of interaction have also been detected in addition to the formation of covalently bound species. Gas phase reactions studied here will explore, via mass-selected ion mobility, reversible and irreversible reactions leading to binding enthalpy and entropy and rate constant determination, respectively, in addition to collisional cross-section determination.

This thesis is concerned with a range of chemical problems which are amenable to theoretical investigation via the application of current methods of computational quantum chemistry. These problems include the calculation of accurate thermochemical data, the prediction of reaction kinetics, the study of molecular potential energy surfaces, and the investigation of molecular structures and binding. The heats of formation (from both atomisation energies and isodesmic schemes) of a set of approximately 120 C1 and C2 fluorocarbons and oxidised fluorocarbons (along with C3F6 and CF3CHFCF2) were calculated with the Gaussian-3 (G3) method (along with several approximations thereto). These molecules are of importance in the flame chemistry of 2-H-heptafluoropropane, which has been proposed as a potential fire retardant with which to replace chloro- and bromofluorocarbons (CFC�s and BFC�s). The calculation of the data reported here was carried out in parallel with the modelling studies of Hynes et al.1-3 of shock tube experiments on CF3CHFCF3 and on C3F6 with either hydrogen or oxygen atoms. G3 calculations were also employed in conjunction with the experimental work of Owens et al.4 to describe the pyrolysis of CFClBr2 (giving CFCl) at a radiation wavelength of 265 nm. The theoretical prediction of the dissociation energy of the two C-Br bonds was found to be consistent with the energy at which carbene production was first observed, thus supporting the hypothesis that the pyrolysis releases two bromine radicals (rather than a Br2 molecule). On the basis of this interpretation of the experiments, the heat of formation of CFClBr2 is predicted to be 184 � 5 kJ mol-1, in good agreement with the G3 value of 188 � 5 kJ mol-1. Accurate thermochemical data was computed for 18 small phosphorus containing molecules (P2, P4, PH, PH2, PH3, P2H2, P2H4, PO, PO2, PO3, P2O, P2O2, HPO, HPOH, H2POH, H3PO, HOPO and HOPO2), most of which are important in the reaction model introduced by Twarowski5 for the combustion of H2 and O2 in the presence of phosphine. Twarowski reported that the H + OH recombination reaction is catalysed by the combustion products of PH3 and proposed two catalytic cycles, involving PO2, HOPO and HOPO2, to explain this observation. Using our thermochemical data we computed the rate coefficients of the most important reactions in these cycles (using transition state and RRKM theories) and confirmed that at 2000K both cycles have comparable rates and are significantly faster than the uncatalysed H + OH recombination. The heats of formation used in this work on phosphorus compounds were calculated using the G2, G3, G3X and G3X2 methods along with the far more extensive CCSD(T)/CBS type scheme. The latter is based on the evaluation of coupled cluster energies using the correlation consistent triple-, quadruple- and pentuple-zeta basis sets and extrapolation to the complete basis set (CBS) limit along with core-valence correlation corrections (with counterpoise corrections for phosphorus atoms), scalar relativistic corrections and spin-orbit coupling effects. The CCSD(T)/CBS results are consistent with the available experimental data and therefore constitute a convenient set of benchmark values with which to compare the more approximate Gaussian-n results. The G2 and G3 methods were found to be of comparable accuracy, however both schemes consistently underestimated the benchmark atomisation energies. The performance of G3X is significantly better, having a mean absolute deviation (MAD) from the CBS results of 1.8 kcal mol-1, although the predicted atomisation energies are still consistently too low. G3X2 (including counterpoise corrections to the core-valence correlation energy for phosphorus) was found to give a slight improvement over G3X, resulting in a MAD of 1.5 kcal mol-1. Several molecules were also identified for which the approximations underlying the Gaussian-n methodologies appear to be unreliable; these include molecules with multiple or strained P-P bonds. The potential energy surface of the NNH + O system was investigated using density functional theory (B3LYP/6-31G(2df,p)) with the aim of determining the importance of this route in the production of NO in combustion reactions. In addition to the standard reaction channels, namely decomposition into NO + NH, N2 + OH and H + N2O via the ONNH intermediate, several new reaction pathways were also investigated. These include the direct abstraction of H by O and three product channels via the intermediate ONHN, giving N2 + OH, H + N2O and HNO + N. For each of the species corresponding to stationary points on the B3LYP surface, valence correlated CCSD(T) calculations were performed with the aug-cc-pVxZ (x = Q, 5) basis sets and the results extrapolated to the complete basis set limit. Core-valence correlation corrections, scalar relativistic corrections and spin orbit effects were also included in the resulting energetics and the subsequent calculation of thermochemical data. Heats of formation were also calculated using the G3X method. Variational transition state theory was used to determine the critical points for the barrierless reactions and the resulting B3LYP energetics were scaled to be compatible with the G3X and CCSD(T)/CBS values. As the results of modelling studies are critically dependent on the heat of formation of NNH, more extensive CCSD(T)/CBS calculations were performed for this molecule, predicting the heat of formation to be 60.6 � 0.5 kcal mol-1. Rate coefficients for the overall reaction processes were obtained by the application of multi-well RRKM theory. The thermochemical and kinetic results thus obtained were subsequently used in conjunction with the GRIMech 3.0 reaction data set in modelling studies of combustion systems, including methane / air and CO / H2 / air mixtures in completely stirred reactors. This study revealed that, contrary to common belief, the NNH + O channel is a relatively minor route for the production of NO. The structure of the inhibitor Nd-(N'-Sulfodiaminophosphinyl)-L-ornithine, PSOrn, and the nature of its binding to the OTCase enzyme was investigated using density functional (B3LYP) theory. The B3LYP/6-31G(d) calculations on the model compound, PSO, revealed that, while this molecule could be expected to exist in an amino form in the gas phase, on complexation in the active site of the enzyme it would be expected to lose two protons to form a dinegative imino tautomer. This species is shown to bind strongly to two H3CNHC(NH2)2+ moieties (model compounds for arginine residues), indicating that the strong binding observed between inhibitor and enzyme is partially due to electrostatic interactions as well as extensive hydrogen bonding (both model Arg+ residues form hydrogen bonds to two different sites on PSO). Population analysis and hydrogen bonding studies have revealed that the intramolecular bonding in this species consists of either single or semipolar bonds (that is, S and P are not hypervalent) and that terminal oxygens (which, being involved in semipolar bonds, carry negative charges) can be expected to form up to 4 hydrogen bonds with residues in the active site. In the course of this work several new G3 type methods were proposed, including G3MP4(SDQ) and G3[MP2(Full)], which are less expensive approximations to G3, and G3X2, which is an extension of G3X designed to incorporate additional electron correlation. As noted earlier, G3X2 shows a small improvement on G3X; G3MP4(SDQ) and G3[MP2(Full)], in turn, show good agreement with G3 results, with MAD�s of ~ 0.4 and ~ 0.5 kcal mol-1 respectively. 1. R. G. Hynes, J. C. Mackie and A. R. Masri, J. Phys. Chem. A, 1999, 103, 5967. 2. R. G. Hynes, J. C. Mackie and A. R. Masri, J. Phys. Chem. A, 1999, 103, 54. 3. R. G. Hynes, J. C. Mackie and A. R. Masri, Proc. Combust. Inst., 2000, 28, 1557. 4. N. L. Owens, Honours Thesis, School of Chemistry, University of Sydney, 2001. 5. A. Twarowski, Combustion and Flame, 1995, 102, 41.

In the Mars atmosphere, reactivity of trace components is as relevant as that of the major compounds if formation of complex molecules is considered. These are of great importance concerning the existence of a past or future life on Mars. In this context, the gas-phase chemistry of alkylcarbonate and alkyloxalate ions and radicals, ROCOO–/• and ROCOCOO–/•, is investigated for different alkyl substituents R (H, CH3, C2H5, i-C3H7, and t-C4H9). This study describes the structures, stability, and unimolecular dissociation behavior of these species and is achieved by means of mass spectrometric methods and ab initio calculations. Standard heats of formation of the ions and radicals are determined via computational methods, using atomization energies and bond-separation reactions. Vertical charge-transfer experiments are performed to provide evidence for the existence of the radicals under study and the NIDD (ion and neutral decomposition difference) method is employed to determine their reactivity. Typical processes observed involve direct bond cleavages leading to elimination of carbon dioxide. Concerning anionic compounds, classical metastable ion (MI) and collisional activation (CA) experiments enable the determination of their gas-phase behavior. This, in contrast to radicals, is not only constituted by direct bond cleavages, but also by hydride-transfer reaction or carbon monoxide expulsion involving formation of ion-neutral complexes as intermediates. Translational energy loss spectra are also employed to gain more insights concerning the dissociation energetics of CH3OCOO• and CH3OCOO+ formed by vertical charge-transfer of methylcarbonate. This rather unusual method for such study implies a careful evaluation of the error caused by the instrument which may otherwise generate dramatic deviations of the results compared to theory.

Radical reactions play central roles in regulating regional air quality and global climate. Some potentially important gas phase radical reactions are being investigated in this research project, including Cl reactions with acetone, butanone, 3-pentanone, pyridine, and dimethyl selenide (DMSe), HO2 complex formation and dissociation with formic and acetic acids, and reactive and non-reactive quenching of O(1D) by the potent greenhouse gases SO2F2, NF3, and SF5CF3. The involved radicals are generated by laser flash photolysis (LFP). Temporal profiles of either the radical reactant or a product are monitored in "real time" using atomic resonance fluorescence spectroscopy (RF), time-resolved UV-visible absorption spectroscopy (TRUVVAS), or tunable diode laser absorption spectroscopy (TDLAS), allowing kinetic and mechanistic information of these reactions to be obtained. These studies provide new knowledge of the investigated radical reactions and facilitate a better understanding of their significance in atmospheric chemistry.

Boron carbide is produced in a heat resistance furnace using boric oxide and petroleum coke as the raw materials. In this process, a large current is passed through the graphite rod located at the center of the cylindrical furnace, which is surrounded by the coke and boron oxide mixture. Heat generated due to resistance heating is responsible for the reaction of boron oxide with coke which results in the formation of boron carbide. The whole process is highly energy intensive and inefficient in terms of the production of boron carbide. Only 15% charge gets converted into boron carbide. The aim of the present work is to develop a mathematical model for this batch process and validate the model with experiments and to optimize the operating parameters to increase the productivity. To mathematically model the process and understand the influence of various operating parameters on the productivity, existing simple one-dimensional (1-D) mathematical model in radial direction is modified first. Two-dimensional (2-D) model can represent the process better; therefore in second stage of the project a 2-D mathematical model is also developed. For both, 1-D and 2-D models, coupled heat and mass balance equations are solved using finite volume technique. Both the models have been tested for time step and grid size independency. The kinetics of the reaction is represented using nucleation growth mechanism. Conduction, convection and radiation terms are considered in the formulation of heat transfer equation. Fraction of boron carbide formed and temperature profiles in the radial direction are obtained computationally. Experiments were conducted on a previously developed experimental setup consisting of heat resistance furnace, a power supply unit and electrode cooling device. The heating furnace is made of stainless steel body with high temperature ceramic wool insulation. In order to validate the mathematical model, experiments are performed in various conditions. Temperatures are measured at various locations in the furnace and samples are collected from the various locations (both in radial and angular directions) in the furnace for chemical analysis. Also, many experimental data are used from the previous work to validate the computed results. For temperatures measurement, pyrometer, C, B and K type thermocouple were used. It is observed that results obtained from both the models (1-D and 2-D) are in reasonable agreement with the experimental results. Once the models are validated with the experiments, sensitivity analysis of various parameters such as power supply, initial percentage of B4C in the charge, composition of the charge, and various modes of power supply, on the process is performed. It is found that linear power supply mode, presence of B4C in the initial mixture and increase in power supply give better productivity (fraction reacted). In order to have more confidence in the developed models, the parameters of one the computed results in the sensitivity analysis parameters are chosen (in present case, linear power supply is chosen) to perform the experiment. Results obtained from the experiment performed under the same simulated conditions as computed results are found in excellent match with each other.

This thesis is concerned with a range of chemical problems which are amenable to theoretical investigation via the application of current methods of computational quantum chemistry. These problems include the calculation of accurate thermochemical data, the prediction of reaction kinetics, the study of molecular potential energy surfaces, and the investigation of molecular structures and binding. The heats of formation (from both atomisation energies and isodesmic schemes) of a set of approximately 120 C1 and C2 fluorocarbons and oxidised fluorocarbons (along with C3F6 and CF3CHFCF2) were calculated with the Gaussian-3 (G3) method (along with several approximations thereto). These molecules are of importance in the flame chemistry of 2-H-heptafluoropropane, which has been proposed as a potential fire retardant with which to replace chloro- and bromofluorocarbons (CFC�s and BFC�s). The calculation of the data reported here was carried out in parallel with the modelling studies of Hynes et al.1-3 of shock tube experiments on CF3CHFCF3 and on C3F6 with either hydrogen or oxygen atoms. G3 calculations were also employed in conjunction with the experimental work of Owens et al.4 to describe the pyrolysis of CFClBr2 (giving CFCl) at a radiation wavelength of 265 nm. The theoretical prediction of the dissociation energy of the two C-Br bonds was found to be consistent with the energy at which carbene production was first observed, thus supporting the hypothesis that the pyrolysis releases two bromine radicals (rather than a Br2 molecule). On the basis of this interpretation of the experiments, the heat of formation of CFClBr2 is predicted to be 184 � 5 kJ mol-1, in good agreement with the G3 value of 188 � 5 kJ mol-1. Accurate thermochemical data was computed for 18 small phosphorus containing molecules (P2, P4, PH, PH2, PH3, P2H2, P2H4, PO, PO2, PO3, P2O, P2O2, HPO, HPOH, H2POH, H3PO, HOPO and HOPO2), most of which are important in the reaction model introduced by Twarowski5 for the combustion of H2 and O2 in the presence of phosphine. Twarowski reported that the H + OH recombination reaction is catalysed by the combustion products of PH3 and proposed two catalytic cycles, involving PO2, HOPO and HOPO2, to explain this observation. Using our thermochemical data we computed the rate coefficients of the most important reactions in these cycles (using transition state and RRKM theories) and confirmed that at 2000K both cycles have comparable rates and are significantly faster than the uncatalysed H + OH recombination. The heats of formation used in this work on phosphorus compounds were calculated using the G2, G3, G3X and G3X2 methods along with the far more extensive CCSD(T)/CBS type scheme. The latter is based on the evaluation of coupled cluster energies using the correlation consistent triple-, quadruple- and pentuple-zeta basis sets and extrapolation to the complete basis set (CBS) limit along with core-valence correlation corrections (with counterpoise corrections for phosphorus atoms), scalar relativistic corrections and spin-orbit coupling effects. The CCSD(T)/CBS results are consistent with the available experimental data and therefore constitute a convenient set of benchmark values with which to compare the more approximate Gaussian-n results. The G2 and G3 methods were found to be of comparable accuracy, however both schemes consistently underestimated the benchmark atomisation energies. The performance of G3X is significantly better, having a mean absolute deviation (MAD) from the CBS results of 1.8 kcal mol-1, although the predicted atomisation energies are still consistently too low. G3X2 (including counterpoise corrections to the core-valence correlation energy for phosphorus) was found to give a slight improvement over G3X, resulting in a MAD of 1.5 kcal mol-1. Several molecules were also identified for which the approximations underlying the Gaussian-n methodologies appear to be unreliable; these include molecules with multiple or strained P-P bonds. The potential energy surface of the NNH + O system was investigated using density functional theory (B3LYP/6-31G(2df,p)) with the aim of determining the importance of this route in the production of NO in combustion reactions. In addition to the standard reaction channels, namely decomposition into NO + NH, N2 + OH and H + N2O via the ONNH intermediate, several new reaction pathways were also investigated. These include the direct abstraction of H by O and three product channels via the intermediate ONHN, giving N2 + OH, H + N2O and HNO + N. For each of the species corresponding to stationary points on the B3LYP surface, valence correlated CCSD(T) calculations were performed with the aug-cc-pVxZ (x = Q, 5) basis sets and the results extrapolated to the complete basis set limit. Core-valence correlation corrections, scalar relativistic corrections and spin orbit effects were also included in the resulting energetics and the subsequent calculation of thermochemical data. Heats of formation were also calculated using the G3X method. Variational transition state theory was used to determine the critical points for the barrierless reactions and the resulting B3LYP energetics were scaled to be compatible with the G3X and CCSD(T)/CBS values. As the results of modelling studies are critically dependent on the heat of formation of NNH, more extensive CCSD(T)/CBS calculations were performed for this molecule, predicting the heat of formation to be 60.6 � 0.5 kcal mol-1. Rate coefficients for the overall reaction processes were obtained by the application of multi-well RRKM theory. The thermochemical and kinetic results thus obtained were subsequently used in conjunction with the GRIMech 3.0 reaction data set in modelling studies of combustion systems, including methane / air and CO / H2 / air mixtures in completely stirred reactors. This study revealed that, contrary to common belief, the NNH + O channel is a relatively minor route for the production of NO. The structure of the inhibitor Nd-(N'-Sulfodiaminophosphinyl)-L-ornithine, PSOrn, and the nature of its binding to the OTCase enzyme was investigated using density functional (B3LYP) theory. The B3LYP/6-31G(d) calculations on the model compound, PSO, revealed that, while this molecule could be expected to exist in an amino form in the gas phase, on complexation in the active site of the enzyme it would be expected to lose two protons to form a dinegative imino tautomer. This species is shown to bind strongly to two H3CNHC(NH2)2+ moieties (model compounds for arginine residues), indicating that the strong binding observed between inhibitor and enzyme is partially due to electrostatic interactions as well as extensive hydrogen bonding (both model Arg+ residues form hydrogen bonds to two different sites on PSO). Population analysis and hydrogen bonding studies have revealed that the intramolecular bonding in this species consists of either single or semipolar bonds (that is, S and P are not hypervalent) and that terminal oxygens (which, being involved in semipolar bonds, carry negative charges) can be expected to form up to 4 hydrogen bonds with residues in the active site. In the course of this work several new G3 type methods were proposed, including G3MP4(SDQ) and G3[MP2(Full)], which are less expensive approximations to G3, and G3X2, which is an extension of G3X designed to incorporate additional electron correlation. As noted earlier, G3X2 shows a small improvement on G3X; G3MP4(SDQ) and G3[MP2(Full)], in turn, show good agreement with G3 results, with MAD�s of ~ 0.4 and ~ 0.5 kcal mol-1 respectively. 1. R. G. Hynes, J. C. Mackie and A. R. Masri, J. Phys. Chem. A, 1999, 103, 5967. 2. R. G. Hynes, J. C. Mackie and A. R. Masri, J. Phys. Chem. A, 1999, 103, 54. 3. R. G. Hynes, J. C. Mackie and A. R. Masri, Proc. Combust. Inst., 2000, 28, 1557. 4. N. L. Owens, Honours Thesis, School of Chemistry, University of Sydney, 2001. 5. A. Twarowski, Combustion and Flame, 1995, 102, 41.

Deux axes de recherche relativement différents ont été explorés au cours de ce travail de thèse à la frontière entre ingénierie, géosciences et science des matériaux, par la valorisation de déchets ferreux et la synthèse de phosphates de fer pour le stockage de l’énergie. Le fil conducteur de ce travail a été l’utilisation de notre connaissance des processus RedOx et hydrothermaux naturels comme source d’inspiration pour le développement de nouveaux procédés à visée industrielle.Dans une première partie qui constitue l’axe principal de mes travaux de recherche, nous proposons un procédé de valorisation de déchets et coproduits industriels ferreux par l’extraction d’oxydes de fer et la production d’hydrogène. Ce procédé s’inspire des processus de production naturelle d’hydrogène observés au niveau de dorsales océaniques où le fer ferreux inclus dans des minéraux magmatiques est oxydé par l’eau à l’état liquide entrainant, en particulier, production d’hydrogène (H2) et formation d’un oxyde de fer, la magnétite Fe3O4. La transposition de ces processus à des déchets anthropiques contenant du fer réduit (Fe métal ou fer ferreux, FeO) ouvre la voie à un procédé de valorisation de ces matériaux. Une première étude théorique s’est concentrée sur l’identification des mécanismes d’oxydation de FeO en conditions hydrothermales et l’intérêt de l’utilisation d’acides organiques dilués. Le procédé a ensuite été appliqué aux laitiers d’aciérie, coproduits de l’industrie sidérurgique composés majoritairement de calcium mais contenant également 15 à 25 % de FeO. Le résultat majeur de cette étude réside dans la production de magnétite de tailles nanométriques, à très forte valeur ajoutée. Des propositions ont également été faites pour compléter la valorisation du laitier avec en particulier le stockage de CO2 sous forme minérale (à basse et haute température).La deuxième partie de cette thèse concerne la synthèse de phosphates de fer pour le stockage de l’énergie. Pour les applications mobiles, ce stockage est aujourd’hui majoritairement réalisé par des batteries électrochimiques au lithium pour lesquelles les phosphates de fer au lithium présentent un intérêt particulier. Le stockage et transfert de l’énergie se base sur la capacité de ces matériaux à insérer et extraire des éléments porteurs de charge (Li+, Na+) dans sa structure de type olivine, faisant varier sa composition chimique de LiFe2+PO4 (LFP) et Fe3+PO4 (FP). L’objectif principal de cette seconde partie est de proposer une nouvelle voie de synthèse de FP divisée en deux phases avec 1) la synthèse d’un précurseur de la famille du sarcopside, Fe2+3(PO4)2 (minéral rare à structure pseudo-olivine) et 2) son oxydation isostructurale sous air conduisant à la formation du matériau FP recherché. Cette réaction se caractérise par la migration à l’état solide (exsolution) d’un tiers du fer vers la surface et sa précipitation sous forme d’hématite, Fe2O3. Dans la recherche de nouveaux matériaux pour électrode positive de batterie, ce procédé a également été exporté à l’oxydation du minéral maricite NaFe2+PO4 vers la phase Na3Fe3+2(PO4)3. / Two areas of research relatively different were investigated during this Ph.D thesis at the boundary between process engineering, geosciences and material sciences through the valorization of ferrous byproducts and the synthesis of iron phosphates for energy storage applications. The common thread of this work was the use of our knowledge in redox and natural hydrothermal processes as a geo-inspiration source for the development of novel industrial processes.In the first and main part of this manuscript, we propose a novel valorization path for ferrous wastes and byproducts through the recovery of iron oxides and the production of hydrogen. This process is inspired from natural hydrogen production observed in mid-Atlantic ridges where ferrous iron content of magmatic minerals is oxidized by liquid water leading, among others, to hydrogen production and magnetite (Fe3O4) formation. Applying the same oxidation process to ferrous byproducts (metal Fe or ferrous iron FeO) enables the development of a novel valorization path. A first study is conducted on the identification of hydrothermal oxidation mechanism of reagent grade FeO and the influence of mild acetic acid on oxidation kinetic. This process is then applied to steel slag, a steel-making byproduct mainly composed of calcium but also of 15 to 25 %w FeO. The major result of this study was found in the characterization of magnetite as nanoparticles, a highly valuable product. In addition, we propose to complete the valorization process of steel slag by performing mineral CO2 sequestration (at room and high temperature).The second part of this Ph.D manuscript is dedicated to the synthesis of iron phosphates for energy storage applications. For mobile applications, this storage is nowadays mainly performed by lithium batteries. For these devices, a particular interest is given to lithium iron phosphates as positive electrode material for their ability to insert and disinsert lithium in its olivine related structure, modifying its chemical composition from LiFe2+PO4 (LFP) to Fe3+PO4 (FP). The main goal of this second part is to propose a novel synthesis path for FP through a two-step process with 1) the synthesis of a sarcopside related material Fe2+3(PO4)2 (a pseudo-olivine structured rare mineral) and 2) its isostructural oxidation in air to form the targeted FP material. This oxidation step leads to the partial solid migration of iron (exsolution mechanism) from the core to the surface and its precipitation as hematite, Fe2O3. For the sake of new materials as positive electrode, this process is then applied to the oxidation of maricite NaFe2+PO4 into Na3Fe3+2(PO4)3.