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

Sintesi Tecnica: Potenziale di Interazione di tipo Scarf Modificato per gli Autovalori dell'Energia Vibrazionale e le Costanti di Equilibrio Chimico

Definizione del Problema
La determinazione accurata delle proprietà termodinamiche di equilibrio e delle costanti di equilibrio chimico è critica per la comprensione dei processi chimici e fisico-chimici. Tuttavia, gli approcci computazionali esistenti, come i metodi ab initio, la teoria del funzionale della densità e i metodi coupled-cluster, richiedono spesso estese calcoli della struttura elettronica e frequenze vibrazionali altamente accurate. Tali metodi possono essere computazionalmente costosi e possono soffrire di approssimazioni numeriche, in particolare per sistemi con interazioni deboli. Sebbene i potenziali di interazione molecolare analiticamente risolvibili (ad esempio, Morse, Kratzer, Pöschl-Teller) offrano alternative efficienti per derivare spettri vibrazionali e quantità termodinamiche, le investigazioni sistematiche che utilizzano quadri di tipo Scarf modificato per formulazioni termodinamiche di equilibrio e costanti di equilibrio chimico rimangono limitate nella letteratura.

Metodologia
Lo studio propone un potenziale di interazione di tipo Scarf modificato (MST) per colmare tali lacune. La metodologia procede attraverso i seguenti passaggi analitici e computazionali:

1. Definizione del Potenziale: Un potenziale di tipo Scarf modificato, $U(\rho)$, viene definito utilizzando funzioni iperboliche che coinvolgono parametri di forza ($U_0, U_1, U_2$) e un parametro di screening ($\alpha$). I parametri del potenziale sono vincolati dalle condizioni al contorno molecolari, garantendo che il modello soddisfi i requisiti fisici per l'energia di dissociazione ($D_e$), la lunghezza di legame all'equilibrio ($\rho_e$), la costante di rigidezza ($k$) e l'accoppiamento rotazione-vibrazione ($\alpha_e$).
2. Soluzione Analitica: L'equazione di Schrödinger radiale per il potenziale MST viene risolta analiticamente utilizzando il metodo Parametrico Nikiforov-Uvarov (PNU). Questo approccio trasforma l'equazione differenziale in una forma di tipo ipergeometrico per derivare un'espressione esplicita per gli autovalori dell'energia vibrazionale ($E_n$).
3. Formulazione Termodinamica: Gli autovalori dell'energia derivati vengono utilizzati per costruire la funzione di ripartizione vibrazionale ($Q_A$). Per migliorare l'accordo con i dati di riferimento, viene introdotto un parametro di fitting ($\xi$) a livello di funzione di ripartizione. La funzione di ripartizione molecolare totale si ottiene combinando le componenti vibrazionali, rotazionali e traslazionali.
4. Derivazione delle Proprietà: Utilizzando la funzione di ripartizione totale, lo studio deriva le espressioni per l'entalpia molare ($H$), l'energia libera di Gibbs ($G$) e la costante di equilibrio chimico ($K_p$) per reazioni reversibili che formano molecole diatomic etonucleari.
5. Validazione: Il modello è validato rispetto a dati sperimentali e di riferimento per dieci molecole diatomic (Br$_2$, F$_2$, N$_2$, O$_2$, P$_2$, S$_2$, BrF, NO, PO, SO). La validazione include:
   • Confronto tra le curve di energia potenziale e gli autovalori dell'energia vibrazionale rispetto ai dati Rydberg-Klein-Rees (RKR).
   • Confronto tra le costanti di equilibrio chimico e i dati termochimici NIST-JANAF.
   • Calcolo della Deviazione Percentuale Assoluta Media (MPAD) per quantificare l'accuratezza.

Contributi Chiave e Risultati
Il documento presenta i seguenti risultati specifici:

• Autovalori dell'Energia: L'espressione analitica per gli autovalori dell'energia vibrazionale riproduce con successo i dati sperimentali RKR. Per Br$_2$, PO e SO, i valori MPAD per le energie vibrazionali sono rispettivamente 2,669%, 0,547% e 1,86%.
• Superfici di Energia Potenziale: Il potenziale MST modella accuratamente i profili di energia potenziale. I valori MPAD per le curve di energia potenziale sono 2,009% (Br$_2$), 0,245% (PO) e 2,972% (SO), con PO che mostra l'accordo più elevato.
• Costanti di Equilibrio Chimico: Il modello dimostra una forte capacità di predire le costanti di equilibrio chimico ($K_p$) nell'intervallo di temperatura 300–6000 K.
  • Per BrF e NO, i risultati teorici corrispondono strettamente ai dati NIST-JANAF con valori MPAD rispettivamente di 1,33% e 1,93%.
  • Per PO e SO, si osservano deviazioni alle basse temperature, ma il modello converge con i dati di riferimento a temperature moderate e alte, mantenendo un accordo favorevole.
• Consistenza Termodinamica: Lo studio stabilisce un quadro coerente che collega il potenziale di interazione MST alle quantità termodinamiche (entalpia, energia libera di Gibbs) e alle costanti di equilibrio senza fare affidamento su calcoli ad alto costo della struttura elettronica.

Significatività e Rivendicazioni
Gli autori affermano che il potenziale di interazione di tipo Scarf modificato funge da quadro utile e affidabile per la spettroscopia molecolare, la modellazione termodinamica e gli studi sull'equilibrio chimico delle molecole diatomic. Il lavoro dimostra che i modelli analiticamente trattabili possono fornire descrizioni accurate degli spettri vibrazionali e delle proprietà di equilibrio, evitando i costi computazionali associati ai calcoli su larga scala della struttura elettronica. Lo studio conclude che il potenziale MST è particolarmente efficace per sistemi in cui è di interesse il comportamento all'equilibrio a temperatura moderata o alta, offrendo un approccio teorico trasparente per investigare le interazioni molecolari e le proprietà termodinamiche.