En ce qui concerne les gaz polyatomiques, nous proposons deux hiérarchies distinctes formées d'équations de moments, qui permettent d'obtenir des lois de conservation de la densité de masse, de la quantité de mouvement et de l'énergie totale du gaz. Ces hiérarchies sont généralement coupées à un certain ordre. Une méthode qui fournit une solution appropriée au problème de fermeture est la méthode de la maximisation d'entropie. Nous formulons un problème variationnel et nous explorons en détail le cas physique de 14 moments. On étudie un mélange de gaz polyatomiques dans lequel la fonction de distribution de chaque espèce converge vers une Maxwellienne, chacune avec sa propre vitesse moyenne et température. Les lois pour la densité de masse, de quantité de mouvement et d'énergie peuvent être obtenues. En particulier, les coefficients phénoménologiques de la thermodynamique étendue peuvent être déterminés à partir des termes sources. On présente pour les mélanges de gaz monoatomiques l'asymptotique diffusive des équations de Boltzmann. Le développement de Hilbert de chaque fonction de distribution donne deux équations. La première équation permet d'affirmer que le mélange est proche de l'équilibre. La deuxième équation est une équation fonctionnelle linéaire en la variable de vitesse. Nous prouvons l'existence d'une solution de cette équation. D'une part, lorsque les masses moléculaires sont égales, les techniques introduites par Grad peuvent être utilisés. D'autre part, nous proposons une nouvelle approche qui est valable lorsque les masses moléculaires sont différentes. / Considering polyatomic gases, we first propose two independent hierarchies of the moment equations, which allow to obtain conservation laws for mass density, momentum and total energy of a gas. Such hierarchies are usually truncated at some order. A method which provides an appropriate solution to the closure problem is the maximization of entropy method. We formulate a variational problem and explore in detail the physical case of 14 moments. We study mixtures of polyatomic gases in which the distribution function of each species converges towards a Maxwellian distribution function, each with its own bulk velocity and temperature. Balance laws for mass density, momentum and energy can be obtained. In particular, the phenomenological coefficients of extended thermodynamics can be determined from the source terms. Regarding mixtures of monatomic gases, we discuss the diffusion asymptotics of the Boltzmann equations. The Hilbert expansion yields two equations. The first equation allows to state that the mixture is close to equilibrium. The second equation is a linear functional equation in the velocity variable. We prove the existence of a solution to this equation. On the one hand, when molecular masses are equal, the techniques introduced by Grad can be used. On the other hand, we propose a new approach, which only holds when molecular masses are different.

A high-order macroscopic model for the accurate description of rarefied polyatomic gas flows is introduced based on a simplified kinetic equation. The different energy exchange processes are accounted for with a two term collision model. The order of magnitude method is applied to the primary moment equations to acquire the optimized moment definitions and the final scaled set of Grad's 36 moment equations for polyatomic gases. The proposed kinetic model, which is an extension of the S-model, predicts correct relaxation of higher moments and delivers the accurate Prandtl (Pr) number. Also, the model has a proven H-theorem. At the first order, a modification of the Navier-Stokes-Fourier (NSF) equations is obtained, which shows considerable extended range of validity in comparison to the classical NSF equations in modeling sound waves. At third order of accuracy, a set of 19 regularized PDEs (R19) is obtained. Furthermore, the terms associated with the internal degrees of freedom yield various intermediate orders of accuracy, a total of 13 different orders. Attenuation and speed of linear waves are studied as the first application of the many sets of equations. For frequencies were the internal degrees of freedom are effectively frozen, the equations reproduce the behavior of monatomic gases. Thereafter, boundary conditions for the proposed macroscopic model are introduced. The unsteady heat conduction of a gas at rest and steady Couette flow are studied numerically and analytically as examples of boundary value problems. The results for different gases are given and effects of Knudsen numbers, degrees of freedom, accommodation coefficients and temperature dependent properties are investigated. For some cases, the higher order effects are very dominant and the widely used first order set of the Navier Stokes Fourier equations fails to accurately capture the gas behavior and should be replaced by a higher order set of equations. / Graduate / 0346, 0791, 0548, 0759 / behnamr@uvic.ca

The problem of a vibrating-rotating polyatomic molecule is treated, with emphasis given to the case of molecules with C_3v symmetry. It is shown that several of the gross features of the rotational spectra of polyatomic molecules in excited vibrational states can be predicted by group theoretical considerations. Expressions for the rotational transition frequencies of molecules of C_3v symmetry in the ground vibrational state, singly excited degenerate vibrational states, and doubly excited degenerate vibrational states are given. The derivation of these expressions by fourth order perturbation theory as given by Amat, Nielsen, and Tarrago is discussed. The ground and V_10=1 rotational spectra of propyne have been investigated in the 17 to 70 GHz, and 17 to 53 GHz regions, respectively, and compared with predictions based on higher frequency measurements. The V_9=1 and V_10=2 rotational spectra of propyne have been investigated and assigned for the first time. A perturbation of the V_9=1 rotational spectra for K=-l has been discovered and discussed.

Some collision cascades, induced by keV polyatomic projectiles, result in the emission of multiple secondary ions. Such co-emissions imply that the ejecta originate from molecules co-located within the nano-volume perturbed by a single projectile impact. The relevance for the chemical analysis of nano-domains depends on the effectiveness of the projectile to cause co-emission of two or more secondary ions. This research examines how projectile characteristics, i.e. the energy and number of constituent atoms in the projectile, influence multiple secondary ion emission, or "superefficient" events. In addition we examine the relevance of this technique for nanostructure investigation. Yields have been measured for multi-ion emission events as a function of projectile characteristics. The data show that some collision cascades are "superefficient". For example, in a four-ion emission event, the yield for the phenylalanine quasi-molecular ion is two orders of magnitude larger from Au4+ impacts than from equal velocity Au+ projectiles. Yields for the co-emission of two phenylalanine quasi-molecular ions from "super-efficient" events have been measured. This case is particularly productive in that the detection of two analytically significant ions is recorded from a single event. Large increases (one to two orders of magnitude) in co-emitted ion yields were observed with increasing projectile energy and complexity. Correlation coefficients were calculated for the co-emission of two Ph ions, their behavior suggests differences in emission pathways for bombardment by atomic and polyatomic projectiles. Finally, we use this methodology to investigate surface structural effects on the occurrence of "super-efficient" events. The results indicate that it is possible to distinguish between two phases of a chemical compound although the stoichiometry remains the same. These results confirm previous predictions concerning the chemical nature of these "super-efficient" events. Also shown is that they are sensitive to the surface nanoenvironment. This approach extends the technology of Secondary Ion Mass Spectrometry by providing a methodology for probing surface nano-domains at the sub100 nm level.

Some collision cascades, induced by keV polyatomic projectiles, result in the emission of multiple secondary ions. Such co-emissions imply that the ejecta originate from molecules co-located within the nano-volume perturbed by a single projectile impact. The relevance for the chemical analysis of nano-domains depends on the effectiveness of the projectile to cause co-emission of two or more secondary ions. This research examines how projectile characteristics, i.e. the energy and number of constituent atoms in the projectile, influence multiple secondary ion emission, or "superefficient" events. In addition we examine the relevance of this technique for nanostructure investigation. Yields have been measured for multi-ion emission events as a function of projectile characteristics. The data show that some collision cascades are "superefficient". For example, in a four-ion emission event, the yield for the phenylalanine quasi-molecular ion is two orders of magnitude larger from Au4+ impacts than from equal velocity Au+ projectiles. Yields for the co-emission of two phenylalanine quasi-molecular ions from "super-efficient" events have been measured. This case is particularly productive in that the detection of two analytically significant ions is recorded from a single event. Large increases (one to two orders of magnitude) in co-emitted ion yields were observed with increasing projectile energy and complexity. Correlation coefficients were calculated for the co-emission of two Ph ions, their behavior suggests differences in emission pathways for bombardment by atomic and polyatomic projectiles. Finally, we use this methodology to investigate surface structural effects on the occurrence of "super-efficient" events. The results indicate that it is possible to distinguish between two phases of a chemical compound although the stoichiometry remains the same. These results confirm previous predictions concerning the chemical nature of these "super-efficient" events. Also shown is that they are sensitive to the surface nanoenvironment. This approach extends the technology of Secondary Ion Mass Spectrometry by providing a methodology for probing surface nano-domains at the sub100 nm level.

<p>We construct two independent hierarchies of&nbsp;moment equations and we apply the maximum&nbsp;entropy principle for polyatomic gases. We&nbsp;formulate multivelocity and multitemperature&nbsp;model of Eulerian polyatomic gases starting from&nbsp;kinetic theory, that is compared in the&nbsp;neighborhood of global equilibrium state to the&nbsp;models based on extended thermodynamics. We&nbsp;analyze diffusion asymptotics of the Boltzmann&nbsp;<br />equations for mixtures of monatomic gases.</p> / <p>Конструишу се две независне хијерархије<br />једначина момената и примењује се принцип<br />максимума ентропије за вишеатомске гасове.<br />Формира се вишебрзински и вишетемпературни<br />модел Ојлерових вишеатомских гасова полазећи<br />од кинетичке теорије и добијени модел се<br />пореди у околини стања глобалне равнотеже са<br />моделом проширене термодинамике. Анализира<br />се дифузиона асимптотика Болцманових<br />једначина за мешавине једноатомских гасова.</p> / <p>Konstruišu se dve nezavisne hijerarhije<br />jednačina momenata i primenjuje se princip<br />maksimuma entropije za višeatomske gasove.<br />Formira se višebrzinski i višetemperaturni<br />model Ojlerovih višeatomskih gasova polazeći<br />od kinetičke teorije i dobijeni model se<br />poredi u okolini stanja globalne ravnoteže sa<br />modelom proširene termodinamike. Analizira<br />se difuziona asimptotika Bolcmanovih<br />jednačina za mešavine jednoatomskih gasova.</p>

Thesis (MSc)--University of Stellenbosch, 2011. / ENGLISH ABSTRACT: The selective excitation of an arbitrary vibrational level of a polyatomic molecule, without passage through an intermediary electronic excited state is demonstrated. This was achieved by simulating the interaction of a shaped, femtosecond pulse with one vibrational mode of the molecule. The carrier frequency of the pulse is chosen near resonant to the ground-to- rst-excited vibrational transition of the mode, and the pulse shape is optimized via closed-loop feedback. The simulation concentrates on the rst few vibrationally excited states since the density of states is still low, thus ensuring that the inter-vibrational decoherence time is relatively long compared to the pulse length. While various molecules were investigated this study focuses onUF6 for which detailed spectroscopic data for the v3 vibrational mode is available in literature. A multilevel model was developed and can be adapted for any number of levels. The model reported here was limited to a vibrational quantum number of four. The spectroscopic data included anharmonic splitting as well as forbidden transitions. The effect of rotational levels was not included. A density matrix approach was followed because this will allow for the introduction of dephasing of the coherent excitation via thermalizing collisions with the reservoir, as well as inter-vibrational relaxation. The time evolution of the density matrix is given by the Von Neumann equations. / AFRIKAANSE OPSOMMING: Die selektiewe opwekking van 'n arbitrêre vibrasionele vlak van 'n poliatomies molekule sonder oorgang na 'n intermediëre elektroniese opgewekte toetstand word gedemonstreer. Dit was bereik deur die interaksie te simuleer van 'n gevormde, femtosekonde pulse met een vibrasionele mode van 'n molekule. Die draer frekwensie van die pulse is so gekies dat dit naby resonansie van die grond-tot-eerste-opgewekte vibrasionele oorgang van die mode is, die puls vorm word geoptimeer deur 'n geslote-lus terugvoer. Die simulasie konsentreer op die eerste paar vibrasionele opgewekte toestande, omdat die digtheid van toestande nog steeds laag is, dus verseker dit dat inter-vibrasionele de-koherensie tyd relatief lank is in vergelyking met die puls se lengte. Verskillende molekules was ondersoek vir die studie. Die fokus is op UF6 waarvoor gedetaileerde spektroskopiese data vir die v3 vibrasionele beskikbaar is in die literatuur. 'n Multivlak model was ontwikkel en kan aangepas word vir enige aantal van vlakke. Die model wat hier aangemeld is, is beperk tot die vibrasionele kwantum getal van vier. Die spektroskopiese data het anharmonies splitting so wel as nie toegelaatbare oorgange bevat. Die effek van rotasionele vlakke was nie in berekening geneem nie. 'n Digtheids matriks benadering was gevolg, omdat dit toelaat vir die dekoherensie. Die tyd evolusie van die digtheids matriks word gegee deur die Von Neumann vergelykings.

Cette thèse porte sur l’analyse isotopique de poussières interplanétaires riches en carbone (des micrométéorites ultracarbonées), pour mieux comprendre les processus de fractionnement isotopique des éléments légers (hydrogène et azote) observés dans les phases organiques de la matière extraterrestre primitive. Il est possible, au sein d’une collection de micrométéorites provenant des régions centrales du continent Antarctique, d’identifier des micrométéorites ultracarbonées (UCAMMs – pour UltraCarbonaceous Antarctic MicroMeteorites) qui contiennent environ dix fois plus de carbone que les météorites primitives. Ces particules rares (environ 1% des particules collectées) ouvrent la possibilité d’étudier en laboratoire des particules cométaires de grandes tailles (100-200 microns). Des analyses minéralogiques, chimiques et structurales sur deux fragments de micrométéorites ultracarbonées ont été réalisées dans le cadre de collaborations à l’aide de techniques de microscopie électronique à balayage, microsonde électronique, microspectroscopie infrarouge (synchrotron SOLEIL). Le travail principal de cette thèse concerne l’analyse isotopique de deux UCAMMs par spectrométrie de masse à émission ionique secondaire (SIMS) à l’aide de la microsonde ionique NanoSIMS. Les développements instrumentaux effectués en collaboration entre les équipes du CSNSM et de l’Institut Curie permettent de résoudre les interférences moléculaires avec un pouvoir de séparation en masse supérieur à 20 000 en conservant une sensibilité compatible avec les mesures envisagées. Ainsi, il a été possible de mesurer pour la première fois la composition isotopique de l’hydrogène avec les ions poly-atomiques 12C2D- et 12C2H- et de comparer sa distribution spatiale avec celle de l’azote (12C15N-/12C14N-) obtenue sur la même surface avec le même champ magnétique. Un protocole de fabrication d’une série d’échantillons standards a été développé à partir de polymères isotopiquement marqués en deutérium afin de déterminer le fractionnement instrumental du NanoSIMS démontrant la possibilité d’effectuer des mesures précises de la composition isotopique de l’hydrogène dans la matière organique avec des ions polyatomiques (rapports CD-/CH- et C2D-/C2H-). Les cartographies isotopiques à haute résolution spatiale (200 nm) en C2D-/C2H- et C15N-/C14N- montrent une distribution très hétérogène du deutérium avec des valeurs extrêmes allant jusqu’à 20 fois la valeur des océans terrestres, alors que les rapports 15N/14N restent globalement proches (à 20% près) de la valeur de l’atmosphère terrestre. Les mesures des rapports N/C en imagerie ionique effectuées au NanoSIMS confirment les valeurs mesurées par microsonde électronique et montrent que la matière organique des UCAMMs est riche en azote. Les comparaisons des images isotopiques entre elles et avec les images obtenues par les autres techniques montrent la matière organique des micrométéorites ultracarbonées est constituée de différents composants. Le composant principal ne montre pas de corrélation entre les excès en D et ceux en 15N. Une fraction mineure de la surface analysée montre des excès corrélés en D et en 15N, sans que la composition élémentaire ni la structure de cette composante ne diffèrent significativement du reste de la particule. Enfin, une seconde composante mineure présentant des enrichissements modérés en D et des rapports 15N/14N inférieurs à la valeur de l’atmosphère terrestre. L’ensemble des données peut être expliqué en supposant que les micrométéorites ultracarbonées proviennent de la surface d’objets glacés transneptuniens. Les caractéristiques élémentaires et isotopiques observées dans la matière organique des UCAMMs pourraient résulter du mélange de différentes strates du corps parent ayant subi une irradiation par le rayonnement cosmique galactique à grandes distances héliocentriques. / This thesis covers the isotopic analysis of carbon-rich interplanetary dust (ultracarbonaceous micrometeorites), in order to better understand the isotopic fractionation process of light elements (hydrogen and nitrogen) observed in the organic phases of primitive extraterrestrial matter. It is possible, within a collection of micrometeorites coming from the central regions of the Antarctic continent, to identify ultra-carbonaceous micrometeorites (UCAMMs – for Ultra-Carbonaceous Antarctic MicroMeteorites) which contain about ten times more carbon than the primitive meteorites. These particles are extremely rare and open the possibility to study in laboratory cometary particles of large sizes (100-200 microns). This thesis focuses on the isotopic analyses of UCAMMs, in order to understand the isotopic fractionation processes of light elements (hydrogen and nitrogen) observed in the organic matter of the primitive solar system. Mineralogical, chemical and structural analyses on two fragments of ultra-carbonaceous micrometeorites were performed in the framework of collaborations using techniques of scanning electron microscopy, electron microprobe, infrared microspectroscopy (SOLEIL synchrotron). The main work of this thesis concerns the isotopic analysis of two UCAMMs by Secondary Ion Mass Spectrometry (SIMS) using the NanoSIMS ion microprobe. The instrumental developments performed in collaboration between the CSNSM and the Curie Institute teams allow to resolve molecular interferences with a mass resolving power higher than 20 000 retaining a sensibility compatible with the measurements considered here. So, it has been possible to measure for the first time the hydrogen isotopic composition with the polyatomic ions 12C2D- and 12C2H- and to compare its spatial distribution with that of the nitrogen (12C15N-/12C14N-) obtained on the same surface with the same magnetic field. A protocol for producing a series of standard samples was developed from polymers isotopically enriched in deuterium in order to determine the instrumental fractionation of the NanoSIMS showing the possibility to perform precise measurements of the hydrogen isotopic composition in the organic matter with polyatomic ions (CD-/CH- and C2D-/C2H- ratios). The isotopic maps at high spatial resolution (200 nm) in C2D/C2H and C15N/C14N show a very heterogeneous distribution in deuterium with extreme values going up to 20 times the terrestrial ocean value, whereas the 15N/14N ratios remain globally close (within 20%) to the terrestrial atmosphere value. Measurements of N/C ratios in ion imaging performed at the NanoSIMS confirm the values measured by electron microprobe and show that the organic matter of UCAMMs is rich in nitrogen. Comparisons between isotopic images with each other and with the images obtained by the other techniques show that the organic matter of ultracarbonaceous micrometeorites is constituted of different components. The main component doesn’t exhibit a correlation between the D and 15N excesses. A minor component of the analyzed surface exhibits correlated excesses in D and in 15N, without neither the elemental or structural composition of this component significantly differs from the rest of the grain. Finally, a minor component showing moderate enrichments in D and 15N/14N ratios lower than the terrestrial atmosphere value could have been identified.All the data can be explained assuming that the ultracarbonaceous micrometeorites come from the surface of transneptunian icy objects. The elemental and isotopic characteristics observed in the organic matter of UCAMMs might result from the mixing of different strata of the parent body having undergone irradiation by the galactic cosmic radiation at large heliocentric distances.

L’étude de la dissociation de molécules poly-atomiques en surface est d’une importance à la fois fondamentale et industrielle. La compréhension du mécanisme et la dynamique réactionnelle sous-jacents représente un défi. Comme un système modèle, la dissociation de méthane sur la surface de nickel a fait l’objet de nombreuses études pour élucider les chemins de réaction et le transfert d’énergie parmi les différents degrés de liberté durant la réaction. La mode-spécifique ou liaison-spécifique réactivité pour la dissociation de CH4 sur Ni(111) et Ni(100) ont été mise en évidence récemment par des expériences de pointe du jet moléculaire. Jusqu’à présent, les études théoriques de la dynamique réactionnelle ont été effectuées avec un modèle simplifié dans lequel CH4 est décrit comme une molécule pseudo-diatomique. Le concept d’un groupe méthyle spectateur introduit dans un tel modèle impose des contraintes drastiques. Par exemple, l’indiscernabilité des quatre liaisons C-H de méthane est violée par le fait que la liaison C-H capable de se dissocier se singularise par rapport aux trois autres liaisons inertes. En réalité, n’importe quelle des quatre liaisons est susceptible de se dissocier. Par ailleurs, l’unique mode vibrationnel du modèle pseudo-diatomique ne ressemble à aucun des quatre modes vibrationnels principaux du méthane, qui décrivent tous des mouvements collectifs de plusieurs atomes. Lorsque tous les degrés de liberté sont pris en compte, la dimensionnalité de la surface de l’énergie potentielle pour CH4/Ni(111) est très élevée (15 degrés de liberté pour CH4 et certains degrés de liberté du substrat). Construire une surface de l’énergie potentielle fiable à une telle grande dimension est, en soi, un grand défi. A notre connaissance, ce défi n’a jamais été tenté auparavant pour quelconque réaction d’une molécule poly-atomique sur une surface métallique. En utilisant un champ de force réactif, nous avons développé, dans le présent travail, une surface de l’énergie potentielle qui prend en compte tous les 15 degrés de liberté de CH4 ainsi que ceux des 3 premières couches de NI(111). Des simulations de dynamique moléculaire ont été effectuées pour étudier la dynamique réaction de CH4 sur Ni(111) aussi bien dans son état fondamental vibrationnel que dans un état excité. Ces simulations ont permis de révéler des comportements dynamiques inattendus et très intéressants. / In the present work, we undertook a challenging task, i.e., construction a full-dimension potential energy surface (PES) for a benchmark poly-atomic molecular surface reaction, CH4/Ni(111), by using a reactive force field. Careful appraisal of the PES was made in order to establish the validity of the PES. The differences between the results for the transition state (dissociation barriers and structures) given by our PES and those by DFT calculations do not exceed 15%. The molecular dynamics simulation results obtained by using our PES are compared to experimental results for CH4 dissociation probability on Ni(111). For the vibrationally excited state, v3 (v=1, J=0), the agreement between our simulation results and the experimental ones is excellent. For the ground state, the sticking coefficient is somehow over-estimated because of the under-estimation of the dissociation barrier by about 150 meV with our reactive force field. Nevertheless, the overall agreement between simulation and experiment is pretty good. Within the help of the full-dimensional PES, we have extensively studied some important aspects of reaction dynamics, e.g., the effects of surface impact position, surface temperature, vibrationally excited state, rotationally excited states etc. For CH4 in ground state (v=0, J=0), the investigation of the effect of CH4 impact position shows that the top site is the most reactive one. The surface temperature strongly affects the reactivity of methane, especially in the region of the low incident energy near to the dissociation threshold, while in the high incident energy region, the effect is less important. For CH4 v3 (v=1, J=0), an important coupling between rotation and vibration is found. The rotation of CH4 can enhance its reactivity in the following way. In its ground state (v=0, J=0), CH4 does not rotate during its flight to the surface. In this case, only one of the two lowest C-H bonds pointing initially toward the surface can be cleaved while the two other bonds never break. In v3 (v=1, J=0) vibrational state, due to the rotation induced by vibration-rotation coupling, any of the four H atoms can be dissociated even if it forms a C-H bond which has an unfavorable initial orientation (i.e., pointing away from the substrate). The rotation of CH4 induced by vibration-rotation coupling near the substrate allows for bringing an unfavorable initial orientation of C-H bond to the right one required by a transition state (TS) during the adsorbate’s approaching to the substrate. As the enhanced reactivity of vibrationally excited molecules is concerned, the intuitively limpid and overwhelmingly accepted explanation is that the vibration-induced bond stretching helps bond breaking. Our simulation results show clearly that the vibration-induced CH4 rotation contributes an important part to the enhanced reactivity of a v3 (v=1, J=0) vibrationally excited CH4. A series of simulations to determine the sticking curves for CH4 in the vibrational ground state (ν=0) but excited to higher rotational levels (J=0-12) have also been performed. Due to its small level spacing, the lowest rotational excited states (J=1-3) of CH4 do not affect its reactivity on Ni(111) as observed experimentally. We found that rotation enhances significantly CH4 reactivity on Ni(111) with a deposited rotational energy amounting only to 12% of the dissociation barrier. Moreover, in a hypothetic simulation, we found also very striking evidences that rotation can even promote better dissociation of CH4 on Ni(111) than vibration. In a vibrationally excited CH4, its C-H bonds undergo alternate stretching and compressing and the latter hinders dissociation. In this case, the reactivity is inevitably modulated by vibration phase. However, the centrifugal force due to rotation tends always to stretch the C-H bonds for CH4 in rotationally excited states. / 多原子气相分子的分解，不仅在物理，化学及相关学科有着基本的重要性，而且可以促进工业进程，如工业制氢气。对其涉及的反应，即化学键的断裂与形成，在理解其反应机制和动力学上更是一项挑战。作为多原子气相-固相化学反应中最为典型的反应，甲烷分子在金属镍表面的分解，已经被广泛地研究从而理解其在动力学过程中的能量转化和反应路径。最近，选态分子束实验报道了有关甲烷在镍表面分解反应的重要特征，即模式选择性和化学键选择性。从理论角度来看，以前大多数理论研究都是基于一个简化模型，即将甲烷分子看成是一个赝双原子分子(CH4=RH,其中R=CH3)。在该简化模型中，将甲基团当做一个“spectator”会导致严重的限制性，如四个碳氢键的不可分辨性就被破坏。因为在简化模型中，只有一个可分解的碳氢键而其他三个碳氢键则被保护起来；而在实际的分解反应过程中，甲烷分子的任何一个碳氢化学键应该都有概率被分解掉。此外，在该赝双原子分子模型中，单键伸长振动模式不能类比于甲烷的四个基本振动模式，因为其每种基本振动模式都涉及多个原子的复合运动。如果不将甲烷处理成赝双原子分子，那么该体系（CH4/Ni(111)）的势能面的维度会很高，即甲烷的15 个自由度加上部分基地原子的自由度。欲建立一个如此高维度而且又可靠的势能面，本身就是一个值得挑战的研究任务。据我们所知，目前对多原子分子在金属表面反应的高维度势能面的报道几乎没有。在本论文中，我们运用键序反应力场（REBO），为体系CH4/Ni(111)，首次建立起一个全维度的势能面。该势能面的维度包含甲烷的15 个自由度和3 层基地原子的自由度。在经典分子动力学（和准经典分子动力学）模拟下，我们研究了甲烷处于基态和激发态时在金属表面的分解活性，并发现了一些非常有趣的结果。本论文包含以下六章：第一章：简单介绍了甲烷在过渡金属表面分解的最新进展。在选态分子束试验报道中，我们介绍了一些有关该反应的重要特征，如模式选择性，化学键选择性，表面温度效应，空间效应，旋转激发效应等。在理论工作方面，主要介绍了两个理论研究小组近期在简化模型下的一些量子动力学结果。第二章：对本文所运用的理论方法和近似做了基本的介绍。这些方法主要归纳于两类：i)电子结构计算；ii) 分子动力学模拟。我们重点介绍了这些方法和近似的特征。第三章：我们运用二阶矩近似力场(SMA)和键序反应力场(REBO)模拟了氢分子在金属钯表面的分解反应, 从而验证反应力场在模拟表面化学反应的适用性。该章讨论了在参数化反应力场时的一些影响因素，如有效数据库大小，不同排斥势以及长程作用项等，为对复杂体系的研究提供了有效的帮助。第四章：基于键序反应力场(REBO)，我们首次为CH4/Ni(111) 体系建立起一个全维度势能面(PES)。同时我们对该势能面(REBO(PES))做了全面评估，如比较势能面(REBO(PES))与DFT计算得到的过渡状态结构和与之对应的分解势垒，比较两者对于不同形式相互作用给出的势能变化等。此外，我们还直接模拟了甲烷在基态时的活性，其模拟结果与实验有着很好的符合度，从而进一步地说明了该势能面(REBO(PES))的可靠性。第五章：在全维度势能面下，我们深入地研究了甲烷处于不同状态时在镍表面分解的反应活性，即基态(v=0,J=0)，反对称振动态v3 (v=1,J=0)和旋转激发态(v=0,J=1-12)。对于基态的甲烷，我们定性并定量地分析了表面碰撞位置，表面温度对其分解概率的影响。对于反对称振动态的甲烷，我们观察到振动激发态的甲烷分子反应活性比基态甲烷的反应活性要大大地增强。究其根源在于，平动能量不易转换至旋转自由度，而振动能量则非常容易转入到旋转自由度。我们利用三个定量参数详细地阐述了这种振动耦合转动的重要性。此外，对于甲烷处于旋转激发态时，我们发现其激发状态非常有利于甲烷的分解，尽管其旋转能量只有分解势垒的12%。更为惊奇的是，对于甲烷分子而言，其旋转激发态比振动激发态更有利于其分解。其相应的物理解释是，对于振动激发的甲烷，它的碳氢键处于伸长与收缩的交替中，而后者却阻止其分解。对于旋转激发中的甲烷，其离心力一直促使碳氢键的伸长。第六章：总结和展望。我们总结了本文的主要结论以及给出一些将来需要进行的工作，如同位素效应等。