Potential energy surfaces are fundamental to materials simulation, providing acomprehensive landscape of the energy variations as a function of atomic positions within a system. These surfaces are crucial for understanding the behavior of materials at the atomic and molecular levels. By mapping out the potential energy surface, researchers can predict stable configurations, transition states, and thermodynamic and kinetic observables, essential for the design of novel materials. However, constructing accurate potential energy surfaces presents significant challenges due to the complexity constructing a potential with nearly the same accuracy of quantum mechanical calculations, the transferability to explore the high-dimensional space of atomic configurations, and the efficiency to be evaluated without the need for extensive computational resources. In this work, group theoretical irreducible derivativemethods are used to construct accurate Taylor series expansions of the potential energy surface. A new condition number optimized basis is developed which guarantees no amplification of error at the smallest computational cost allowed by group theory. Using this highly efficient method, the physics of the metal-insulator transition in LuNiO₃ is investigated within both a purely harmonic model and simplified anharmonic model, yielding reasonable and good agreement with the experimental phase transition temperature, respectively. Then the irreducible derivatives approach is used to benchmark the phonon anharmonicity of the several popular machine learning potentials: Gaussian approximation potentials, artificial neural network potentials, and graph neural network potentials. The benchmark indicated that the graph neural network potentials had the ability to accurate reproduce quantum mechanical derivatives up to 5th-order. Finally, insights from the benchmark are used to train an accurate graph neural network potential for strongly correlated UO₂ which yielded excellent agreement with quantum mechanical calculations and experiment. These methodological advancements underscore the potential of combining grouptheoretical approaches with cutting-edge machine learning techniques to enhance the precision and efficiency of materials simulations. By achieving high accuracy in modeling complex phenomena such as phase transitions and phonon anharmonicity, this work paves the way for future studies to not only explore the design of novel materials with unprecedented properties, but to design new machine learning techniques where group theory is built from the ground up.

In this thesis a reduced dimensionality quantum scattering model is applied to the study of polyatomic reactions of type X + CH4 <--> XH + CH3. Two dimensional quantum scattering of the symmetric hydrogen exchange reaction CH3+CH4 <--> CH4+CH3 is performed on an 18-parameter double-Morse analytical function derived from ab initio calculations at the CCSD(T)/cc-pVTZ//MP2/cc-pVTZ level of theory. Spectator mode motion is approximately treated via inclusion of curvilinear or rectilinear projected zero-point energies in the potential surface. The close-coupled equations are solved using R-matrix propagation. The state-to-state probabilities and integral and differential cross sections show the reaction to be primarily vibrationally adiabatic and backwards scattered. Quantum properties such as heavy-light-heavy oscillating reactivity and resonance features significantly influence the reaction dynamics. Deuterium substitution at the primary site is the dominant kinetic isotope effect. Thermal rate constants are in excellent agreement with experiment. The method is also applied to the study of electronically nonadiabatic transitions in the CH3 + HCl <--> CH4 + Cl(2PJ) reaction. Electrovibrational basis sets are used to construct the close-coupled equations, which are solved via Rmatrix propagation using a system of three potential energy surfaces coupled by spin-orbit interaction. Ground and excited electronic surfaces are developed using a 29-parameter double-Morse function with ab initio data at the CCSD(T)/ccpV( Q+d)Z-dk//MP2/cc-pV(T+d)Z-dk level of theory, and with basis set extrapolated data, both corrected via curvilinear projected spectator zero-point energies. Coupling surfaces are developed by fitting MCSCF/cc-pV(T+d)Z-dk ab initio spin orbit constants to 8-parameter functions. Scattering calculations are performed for the ground adiabatic and coupled surface models, and reaction probabilities, thermal rate constants and integral and differential cross sections are presented. Thermal rate constants on the basis set extrapolated surface are in excellent agreement with experiment. Characterisation of electronically nonadiabatic nonreactive and reactive transitions indicate the close correlation between vibrational excitation and nonadiabatic transition. A model for comparing the nonadiabatic cross section branching ratio to experiment is discussed.

Les processus élémentaires hétérogènes à l’interface gaz-solide présentent un intérêt fondamental dans de nombreux domaines tels que la catalyse hétérogène, la chimie atmosphérique et des milieux interstellaires, la rentrée atmosphérique de véhicules spatiaux ou encore la description des interactions plama-paroi. Cette thèse a pour objet l’étude de la dynamique des processus de collision non réactive de l’azote N2 sur une surface de tungstène W(100) et des processus de recombinaison moléculaire de l’hydrogène H2 sur des surfaces de tungstène W(100) et W(110). Leur dynamique quasi classique est simulée au moyen de surfaces d’énergie potentielle préalablement construites à partir de calculs de théorie de la fonctionnelle de la densité. Un potentiel multi-adsorbats est notamment développé pour tenir compte du taux de couverture de surface afin d’étudier la compétition entre la recombinaison directe, de type Eley-Rideal et la recombinaison par « atomes chauds » après diffusion hyperthermique d’un atome sur la surface. / Heterogeneous elementary processes at the gas-solid interface are ofgreat interest in many domains such as heterogeneous catalysis, atmospheric and interstellar media chemistry, spacecraft atmospheric re-entry and plasma-wall interactions description. This thesis focus on the dynamics of nitrogen, N2, non reactive scattering on a tungsten W(100) surface and hydrogen, H2, recombination processes on tungsten surfaces W(100) and W(110). The quasiclassical dynamics of these processes is simulated using potential energy surfaces based on density functional theory calculations. In particular, a multi-adsorbate potential is developed to include surface coverage in the dynamics simulation in order to scrutinize the interplay between both direct abstraction, the so-called Eley-Rideal recombination,and the Hot-Atom recombination process after hyperthermal diffusion on the surface

Les travaux présentés dans cette thèse peuvent être divisés en deux parties. Dans une première partie, nous avons étudié le processus de photoisomérisation dans plusieurs systèmes. Une analyse de structure électronique accompagnée d’un calcul préliminaire de dynamique semi-classique ont été appliqué à un modèle minimal du rétinal afin d’extraire les degrés de libertés les plus importants lors de l’isomérisation. Cela dans le but de construire des surfaces d’énergie potentielle diabatiques pour effectuer une étude de dynamique quantique. Une approche de type dynamique semi-classique a été appliquée à un modèle de moteur moléculaire dans le but d’étudier l’origine de l’uni-directionalité de sa rotation. Finalement, une étude de structure électronique d’un complexe de Rhénium contenant un ligand de type rétinal a été effectué pour étudier l’influence du métal sur la spectroscopie du ligand rétinal. Dans une deuxième partie nous nous sommes intéressés à l’étude des croisements intersystème dans un complexe de Rhénium. Afin de pouvoir apporter une explication à un comportement contrintuitif de ce complexe, nous avons développé un Hamiltonien modèle capable de tenir compte des couplages vibroniques interétats et spin-orbit. Cet Hamiltonien a été testé sur ce-dit système, et nous a permis, grâce à une étude de structure électronique de proposer un mécanisme de relaxation différent de celui proposé expérimentalement. / This thesis can be divided in two parts.In the first one, we have studied the photoisomerization process in several systems. An electronic structure analysis mixed with a preliminary semi-classical dynamics investigation has been applied to a minimal model of the retinal chromophore in order to select the most important degrees of freedom involved in the process. The goal of this is to build diabatic potential energy surfaces in order to conduct quantum dynamics simulations. A semi-classical approach has also been applied to a molecular motor model to study the origin of the unidirectionality of its rotary motion. Finally, an electronic structure of a rhenium complex with a retinal-like ligand has been performed to study the effect of the coordination to a metallic atom on the spectroscopy of the retinal ligand. In the second part, we have investigated the intersystem crossings in a rhenium complex. In order to bring an explanation to an experimentally observed conterintuitive behavior of this complex, we have developed a model Hamiltonian that includes both interstate vibronic coupling and spin-orbit coupling. This Hamiltonian has been tested on the said complex and, in complement to an electronic structure study, allowed us to formulate a decay mechanism different from the one proposed based on experiments.

In this work, experimental and theoretical approaches are applied to the study of chemical reaction dynamics. In Chapter 2, two applications of transition state theory are presented: (1) Application of microcanonical transition state theory to determine the rate constant of dissociation of C2F3I after π∗ ← π excitation. It was found that this reaction has a very fast rate constant and thus is a promising system for testing the statistical assumption of molecular reaction dynamics. (2) A general rate constant expression for the reaction of atoms and molecules at surfaces was derived within the statistical framework of flexible transition state theory. In Chapter 4, a computationally efficient TDDFT approach was found to produce useful potential energy surface landscapes for application to non-adiabatic predissociative dynamics of the molecule CS2 after excitation from the ground state to the singlet C-state. In Chapter 5, ultrafast experimental results of excitation of CS2 to the predissociative neutral singlet C-state is presented. The bandwidth of the excitation laser was carefully tuned to span a two-component scattering resonance with each component differently evolving electronically with respect to excited state character during the quasi-bound oscillation. Scalar time-resolved photoelectron spectra (TRPES) and vector time-resolved photoelectron angular distribution (TRPAD) observables were recorded during the predissociation. The TRPES yield of photoelectrons was found to oscillate with a quantum beat pattern for the photoelectrons corresponding to ionization to the vibrationless cation ground state; this beat pattern was obscured for photoelectron energies corresponding to ionization from the vibrationally excited CS2 cation. The TRPAD data was recorded for two general molecular ensemble cases: with and without a pre-excitation alignment laser pulse. It was found that in the case of ensemble alignment (Chapter 6), the “molecular frame” TRPAD (i.e. TRMFPAD) was able to image the purely valence electronic dynamics of the evolving CS2 C-state. The unaligned ensemble TRPAD observable suffers from excessive orientational averaging and was unable to observe the quantum beat. Engineering efforts were also undertaken to eliminate scattered light background signal (Chapter 7, Appendix A) and improve laser stability as a function of ambient pressure (Appendix B) for TRMFPAD experiments. / Thesis (Ph.D, Chemistry) -- Queen's University, 2012-09-11 22:18:20.89