[en] ENGLISH SUMMARY: In this thesis we have studied the chemical-physical processes involved in the interfaces gas-surface since first principles. Then, we have carried out a study of the electronic structure and a kinetic and dynamic study of all the processes. In particular, we are interested in the theoretical study of some reactions that take place in the surfaces of the spatial vehicles during their Earth's atmospheric (78% N2, 21% O2 and 1% Ar) re-entry. These vehicles reach hypersonic flight velocities, with temperatures between 500 and 1600 K in different parts of the vehicle. So, in this work we have studied the interaction of atomic O and N (the molecular species dissociate because of the high flight velocities) with SiO2-based materials. This information is necessary for the study and design of new thermal protection systems (TPS) to be used in atmospheric terrestrial re-entry vehicles. TPS coatings for space vehicles must have very low catalytic efficiency to reduce the surface heat flux due to atomic recombination effects. Currently, silica-based materials (e.g., reaction-cured glass (RCG) with 94% of SiO2, 4% of B2O3, and 2% of SiB4) are primarily used in TPS surfaces, where catalytic atomic recombination of oxygen and nitrogen (via Eley-Rideal (ER) or Langmuir-Hinshelwood (LH) mechanisms) plays an important role of re-entry heating at thermal nonequilibrium conditions. However, other elementary processes need to be also included in the global kinetic modelling in these conditions: atomic adsorption -desorption, molecular adsorption -desorption, E-R or L-H reactions, etc. Therefore, we are interested in studying the main steps involved in the heterogeneous catalysis over a silica surface, in particular, beta-cristobalite, which is one of the polimorph of silica most stable at high temperatures and low pressures. Published kinetic models for air recombination over SiO2-based materials incorporate rough estimates of several important physicochemical quantities (sometimes used as adjustable parameters), such as, for example, atomic adsorption energies, atomic sticking coefficients, the number of adsorption sites by unit area, and others. In a first step of the work we have to study the O and N adsorption processes. Nevertheless, scarce experimental data are available on both processes, which are usually treated in a similar way by means of an average adsorption energy of 3.5 eV regardless of the kind of silica surface (e.g., Pyrex, quartz, RCG, etc.), although several used values are within the range of 2.6 - 5.5 eV. These data are indirectly derived from kinetic modelling of atomic recombination on SiO2 surfaces to fit some experimental data. Recent experimental measurements of the oxygen recombination coefficient (gamma-O) over silica have shown that this coefficient is 3 times higher for beta- cristobalite than for quartz, therefore this coefficient depends strongly on the silica surface. First of all, we summarized the theoretical concepts employed in this work, starting from the modelling of periodical systems, density functional theory, the corrugation reducing procedure method used to construct the analytical potential energy surface and the classical trajectory method employed to study the dynamics. Chapter 3 summarizes the theoretical study of some of the silica's polymorphs: alpha and beta-quartz, beta-tridimite and beta-cristobalite. A part of calculating the structural parameters we want to check if DFT with wave planes is quite accurate theory in gas surface interactions. Chapter 4 consist on the second stage of the heterogeneous catalysis: the atomic adsorption of O and N over beta-cristobalite. First of all, we show the adsorption results at DFT level. Secondly we show the construction of two multidimensional potential energy surfaces, which have been analyzed and finally we have carried out the dynamics study through classical trajectory calculations with the aim to get information about the sticking probability. In Chapter 5 we have studied the kinetics and dynamics of the Eley-Rideal reaction. We follow the same structure than in Chapter 4. So, first we have carried out some DFT calculations to get some adsorption energies and to construct the potential energy surface through the corrugation reducing procedure method. Then all the PESs are deeply analyzed and the dynamics is carried out. After the classical trajectory calculations we have obtained information about the atomic catalytic recombination coefficient (gamma-O(T)) and the atomic accommodation coefficient (beta-O(T)). Moreover, we have calculated the vibrorotational levels in which the molecule O2 in the gas phase is formed. Finally, we have started a kinetic studied of the E-R reaction employing some statistical methods (Transition State Theory, Variational Transition State Theory and some semiclassical tunnel effect corrections.

[es] El objetivo principal de esta tesis es el estudio de los procesos fisico-químicos que tienen lugar en las interfases gas-superficie sólida, partiendo desde primeros principios y haciendo un estudio de estructura electrónica, dinámico y cinético de los mismos. En particular, se pretende hacer un estudio teórico de algunas de las reacciones producidas en las superficies de los vehículos espaciales durante su reentrada en la atmósfera de la Tierra (78% N2, 21% O2, y 1% Ar). Estos vehículos llegan a velocidades de entre 22 y 30 veces la velocidad del sonido, con temperaturas que oscilan entre 500 y 1600 K en diferentes partes del vehículo, siendo la parte delantera y las alas, las zonas que alcanzan las temperaturas más elevadas. Así, este trabajo se centrará principalmente en el estudio de las reacciones con las especies atómicas de O y N (ya que las especies moleculares se disocian al colisionar con los vehículos espaciales) sobre superficies sólidas del tipo óxido (SiO2). Esta información tiene una gran trascendencia para el diseño de sistemas de protección térmica (TPS) de los vehículos espaciales, ya que estos TPSs deben tener un coeficiente catalítico muy bajo para reducir el flujo de calor que se transfiere a la superficie debido en parte al mecanismo de recombinación atómica sobre la misma. Con la finalidad de evitar estos calentamientos en la superficie se utilizan distintos tipos de materiales en forma de recubrimiento en los TPSs, sobre todo en la parte delantera del vehículo. El trasbordador espacial presenta un sistema de TPS basado en el uso de materiales de revestimiento con una gran capacidad térmica en combinación con un gran aislamiento térmico para evitar la conducción de calor al interior del vehículo. El mejoramiento de estos materiales ha sido el tema de varias investigaciones para obtener más durabilidad, mayor resistencia al choque térmico, menor conductividad térmica, . etc. La importancia de estos sistemas de protección térmica se ha hecho más relevante, desde las misiones Mercury, Gemini y Apollo que usaban unos sistemas de protección ablativos, que dejaban la cápsula en condiciones inutilizables para futuras expediciones. Su importancia se acentuó con el grave accidente que sufrió el Columbia en el año 2003. En el trasbordador americano (Shuttle) se usa una mezcla de materiales reaction-cured glass (RCG) con composición: 92% SiO2, 5% B2O3 y 3% SiB4 en algunas partes del fuselaje. Además de los mecanismos de recombinación atómica también pueden ser importantes los otros procesos elementales que intervienen en la catálisis heterogénea: la difusión, adsorción, las reacciones químicas, la desorción atómica y molecular y la difusión de los productos. Por tanto, nuestro objetivo se centra en el estudio teórico de las distintas etapas de la catálisis siendo nuestra superficie uno de los polimorfos de la sílice, en particular, la beta-cristobalita, ya que es el polimorfo de la sílice estable a temperaturas elevadas y a presiones bajas. Se deben mejorar los modelos cinéticos existentes, ya que tratan la superficie como si fuese totalmente catalítica. Los mecanismos y las constantes de los distintos procesos gas-superficie en condiciones de no-equilibrio no son bien conocidos. De hecho, muchos de estos procesos no se han caracterizado a bajas temperaturas. Por eso, se precisan datos experimentales y teóricos para simular (p. ej., mediante códigos de dinámica de fluidos) los vuelos de los trasbordadores con varios sistemas de protección térmica. Muchos de los modelos cinéticos experimentales para las reacciones en aire sobre varios tipos de sílices (p. ej., Pyrex, cuarzo, RCG,.) asumen un valor promedio de la energía de adsorción del O y N de Ead=3.5 eV, aunque también en otros se usan valores comprendidos entre 2.5 - 5.5 eV. Estas energías de adsorción derivan de modelos cinéticos que ajustan datos experimentales de los procesos Eley-Rideal/Langmuir-Hinshelwood. También existen medidas experimentales sobre el coeficiente de recombinación de oxígeno atómico (gamma-O) sobre sílice, indicando que este coeficiente es un factor de tres veces más bajo en beta-cuarzo que en beta-cristobalita, por tanto, muy dependiente de la estructura de sílice. En primer lugar, se presentan los conceptos fundamentales teóricos empleados en esta tesis, desde la modelización de sistemas periódicos, la teoría del funcional de la densidad, pasando por la explicación del método empleado en la construcción de las superficies de energía potencial interpoladas, terminando con una explicación del método de trayectorias clásicas. Las dos últimas secciones son las más elaboradas en esta tesis, y las que más esfuerzo de programación han requerido. El Capítulo 3 se centra en el estudio teórico de algunos polimorfos de la sílice, en particular alfa y beta-cuarzo, beta-tridimita y beta-cristobalita. El objetivo de este estudio es testear la metodología utilizada, la teoría del funcional de la densidad con ondas planas, y con ella determinar los parámetros estructurales y calcular las energías de cohesión de dichos sólidos, antes del estudio de los procesos reactivos sobre beta-cristobalita. El Capítulo 4 se refiere a la segunda etapa de la catálisis heterogénea, la adsorción atómica de O y N sobre beta-cristobalita; en primer lugar se presentan los cálculos de adsorción realizados a nivel DFT; en segundo lugar se construyen dos superficies de energía potencial multidimensionales (SEP1 y SEP2), las cuales son analizadas con mucho detalle y finalmente, se realiza el estudio dinámico sobre ambas superficies mediante trayectorias clásicas con el objetivo de obtener información detallada del proceso de sticking. El Capítulo 5 tiene como objetivo principal el estudio de la dinámica y la cinética de la reacción Eley-Rideal. Se sigue la misma estructura que en el Capítulo 4; primero se calcula toda la información a nivel DFT necesaria para construir las superficies de energía potencial. En particular, se construyen cuatro modelos de SEPs, teniendo en cuenta si la superficie analítica total es rígida (con dimensión 6) o no rígida (multidimensional). Así pues, se han construido 8 modelos de SEPs. En todas ellas se analiza la topología de las mismas y se comparan con la información DFT. La dinámica mediante trayectorias clásicas se lleva a cabo con todas las SEPs, tanto en condiciones de equilibrio térmico como en condiciones estado específicas. En condiciones de equilibrio térmico se calculan también el coeficiente de recombinación atómico (gamma-O(T)) y el coeficiente de acomodación energético (beta-O(T)), que da cuenta de la fracción de energía que se transfiere al sólido de la energía total accesible de la reacción. Aparte de los coeficientes térmicos, se realiza un estudio detallado de los niveles vibrorrotacionales con los que se forma la molécula de O2 gas, en todas las condiciones iniciales estudiadas. Finalmente, se presenta un estudio cinético preliminar de la reacción E-R, primero a partir de las probabilidades calculadas en la dinámica, y después empleando métodos estadísticos (teoría del estado de transición, teoría variacional y standard, con inclusión de un efecto túnel semiclásico.