Modified Scarf-type interaction potential for vibrational energy eigenvalues and chemical equilibrium constants of selected diatomic molecules

Resumo Técnico: Potencial de Interação do tipo Scarf Modificado para Autovalores de Energia Vibracional e Constantes de Equilíbrio Químico

Definição do Problema
A determinação precisa das propriedades termodinâmicas de equilíbrio e das constantes de equilíbrio químico é crítica para a compreensão de processos químicos e físico-químicos. No entanto, as abordagens computacionais existentes, como métodos ab initio, da teoria do funcional da densidade e de cluster acoplado, frequentemente exigem cálculos extensos de estrutura eletrônica e frequências vibracionais altamente precisas. Esses métodos podem ser computacionalmente caros e podem sofrer de aproximações numéricas, particularmente para sistemas de interação fraca. Embora potenciais de interação molecular analiticamente solúveis (ex: Morse, Kratzer, Pöschl-Teller) ofereçam alternativas eficientes para derivar espectros vibracionais e quantidades termodinâmicas, investigações sistemáticas utilizando estruturas do tipo Scarf modificado para formulações termodinâmicas de equilíbrio e constantes de equilíbrio químico permanecem limitadas na literatura.

Metodologia
O estudo propõe um potencial de interação do tipo Scarf Modificado (MST) para preencher essas lacunas. A metodologia procede através dos seguintes passos analíticos e computacionais:

1. Definição do Potencial: Um potencial do tipo Scarf modificado, $U(\rho)$, é definido usando funções hiperbólicas envolvendo parâmetros de força ($U_0, U_1, U_2$) e um parâmetro de blindagem ($\alpha$). Os parâmetros do potencial são restringidos por condições de contorno moleculares, garantindo que o modelo satisfaça requisitos físicos para energia de dissociação ($D_e$), comprimento de ligação de equilíbrio ($\rho_e$), constante de rigidez ($k$) e acoplamento rotação-vibração ($\alpha_e$).
2. Solução Analítica: A equação de Schrödinger radial para o potencial MST é resolvida analiticamente utilizando o método Paramétrico Nikiforov-Uvarov (PNU). Esta abordagem transforma a equação diferencial em uma forma do tipo hipergeométrica para derivar uma expressão explícita para os autovalores de energia vibracional ($E_n$).
3. Formulação Termodinâmica: Os autovalores de energia derivados são utilizados para construir a função de partição vibracional ($Q_A$). Para melhorar a concordância com dados de referência, um parâmetro de ajuste ($\xi$) é introduzido ao nível da função de partição. A função de partição molecular total é obtida combinando os componentes vibracional, rotacional e translacional.
4. Derivação de Propriedades: Utilizando a função de partição total, o estudo deriva expressões para a entalpia molar ($H$), energia livre de Gibbs ($G$) e a constante de equilíbrio químico ($K_p$) para reações reversíveis que formam moléculas diatômicas heteronucleares.
5. Validação: O modelo é validado contra dados experimentais e de referência para dez moléculas diatômicas (Br$_2$, F$_2$, N$_2$, O$_2$, P$_2$, S$_2$, BrF, NO, PO, SO). As métricas de validação incluem:
   • Comparação das curvas de energia potencial e autovalores de energia vibracional contra dados RKR (Rydberg-Klein-Rees).
   • Comparação das constantes de equilíbrio químico contra dados termoquímicos NIST-JANAF.
   • Cálculo do Desvio Percentual Absoluto Médio (MPAD) para quantificar a precisão.

Principais Contribuições e Resultados
O artigo apresenta os seguintes achados específicos:

• Autovalores de Energia: A expressão analítica para os autovalores de energia vibracional reproduz com sucesso os dados experimentais RKR. Para Br$_2$, PO e SO, os valores de MPAD para as energias vibracionais são 2,669%, 0,547% e 1,86%, respectivamente.
• Superfícies de Energia Potencial: O potencial MST modela com precisão os perfis de energia potencial. Os valores de MPAD para as curvas de energia potencial são 2,009% (Br$_2$), 0,245% (PO) e 2,972% (SO), com o PO apresentando a maior concordância.
• Constantes de Equilíbrio Químico: O modelo demonstra forte capacidade em prever constantes de equilíbrio químico ($K_p$) na faixa de temperatura de 300–6000 K.
  • Para BrF e NO, os resultados teóricos aproximam-se estreitamente dos dados NIST-JANAF, com valores de MPAD de 1,33% e 1,93%, respectivamente.
  • Para PO e SO, observam-se desvios em baixas temperaturas, mas o modelo converge com os dados de referência em temperaturas moderadas e altas, mantendo uma concordância favorável.
• Consistência Termodinâmica: O estudo estabelece um arcabouço consistente ligando o potencial de interação MST a quantidades termodinâmicas (entalpia, energia livre de Gibbs) e constantes de equilíbrio, sem depender de cálculos de estrutura eletrônica de alto custo.

Significância e Alegações
Os autores alegam que o potencial de interação do tipo Scarf Modificado serve como um arcabouço útil e confiável para espectroscopia molecular, modelagem termodinâmica e estudos de equilíbrio químico de moléculas diatômicas. O trabalho demonstra que modelos analiticamente tratáveis podem fornecer descrições precisas de espectros vibracionais e propriedades de equilíbrio, evitando os custos computacionais associados a cálculos de estrutura eletrônica de larga escala. O estudo conclui que o potencial MST é particularmente eficaz para sistemas onde o comportamento de equilíbrio de moderada a alta temperatura é de interesse, oferecendo uma abordagem teórica transparente para investigar interações moleculares e propriedades termodinâmicas.