Theoretical study of the O(3P) + HBr reaction: potential energy surfaces, knetics and dynamics

Neste trabalho, a reação O(3P) + HBr → OH + Br, que se insere no contexto da química atmosférica, mais especificamente nos ciclos de destruição catalítica do ozônio, foi estudada empregando superfícies de energia potencial ab initio de alta qualidade. Construímos superfícies para os estados eletrônicos 3A\" e 3A\', baseadas na expansão de muitos corpos, com potenciais de dois corpos ajustados por funções \"switched-MLJ\" e potenciais de três corpos interpolados utilizando o método RKHS. Os pontos ab initio utilizados foram obtidos no nível de teoria MRCISD+Q/CBS+SO. A superfície final para o estado eletrônico 3A\" foi corrigida por um fator multiplicativo de modo que reproduzisse cálculos benchmark para a altura de barreira, em que foram considerados efeitos relativísticos (escalares e spin-órbita), correlação caroço-valência e correlação de ordem alta (excitações triplas e quádruplas). As principais características da SEP 3A\" são dois mínimos de van der Walls no canal de entrada e saída e um estado de transição com uma barreira energética de 5,01 kcal mol-1. A SEP 3A\' tem um estado de transição linear em 6,45 kcal mol-1. O valor obtido para a entalpia de reação, a 0 K, foi de -15,7 kcal mol-1, muito próximo do valor experimental de -15,14 kcal mol-1. Calculamos constantes de velocidade, no intervalo de 200 a 1000 K, utilizando a teoria do estado de transição variacional com contribuições de tunelamento multidimensional (ICVT/µOMT) e com uma aproximação de espalhamento reativo quântico (QM/JS). Estes valores estão em ótima concordância com os dados experimentais da literatura, em todo intervalo de temperatura em que estão disponíveis: de 221 a 554 K para reação O + HBr e de 295 a 419 K para reação O + DBr. As constantes de velocidade, a 298 K e em cm3 molécula-1 s-1, obtidas para a reação O + HBr são: 3,62·10-14 (ICVT/µOMT) e 3,35·10-14 (QM/JS), enquanto que o valor experimental é 3,78·10-14. A qualidade destes resultados reforça nossa confiança nos procedimentos e aproximações utilizadas, abrindo caminho para a caracterização, em alto nível, de uma grande variedade de reações em fase gasosa.In this work, the O(3P) + HBr → OH + Br reaction, which is relevant to atmospheric chemistry, specially for the catalytic ozone depletion, was studied using high-level ab initio potential energy surfaces. We constructed surfaces for the 3A\" and 3A\' electronic states, based on the many-body expansion, with the two-body potentials adjusted by the switched-MLJ function and the three-body potentials interpolated using the RKHS method. The ab initio points were calculated at the MRCISD+Q/CBS+SO level of theory. The main features of the 3A\" are the presence of two van der Waals minima, at the entrance and exit channels, and a transition state with a barrier height of 5.01 kcal mol-1. The 3A\' PES has a linear transition state at 6.45 kcal mol-1. We obtained the enthalpy of reaction, at 0 K, of -15.7 kcal mol-1, in close agreement with the experimental value of -15.14 kcal mol-1. Rate constants, in the temperature range from 200 to 1000 K, were calculated using the variational transition state theory with contributions of multidimensional tunneling (ICVT/µOMT) and also a quantum reactive scattering approach (QM/JS). Their values are in fair agreement with the experimental data in the literature in the whole temperature range available: from 221 to 554 K for the O + HBr reaction, and from 295 to 419 K for the O + DBr reaction. At 298 K, the calculated rate constants, in cm3 molecule-1 s-1, for the O + HBr reaction are 3.62·10-14 (ICVT/µOMT) and 3.35·10-14 (QM/JS); and the experimental value is 3.78·10-14. The quality of these results reinforces our confidence in the procedures and approximations used, leading to the possibility of high-level characterization of a variety of gas phase reactions