İçeriğe atla

Gibbs Serbest Enerjisinin Değişimi, Basınca ve Sıcaklığa Bağlılık

Teradigma sitesinden

Gibbs Serbest Enerjisinin Değişimi: Basınca ve Sıcaklığa Bağlılık

Gibbs serbest enerjisi (G), sabit basınç ve sıcaklık koşullarındaki kimyasal sistemlerin evrimini yöneten termodinamik potansiyeldir. 1870’lerde Josiah Willard Gibbs tarafından formüle edilen bu işlev, entalpi (H), sıcaklık (T) ve entropi (S) arasındaki hassas ilişkiyi tek bir büyüklükte bir araya getirir. Sistem durumunu basınç (P) ve sıcaklığa (T) bağlı olarak tanımlayan bu potansiyelin, söz konusu değişkenlere olan bağımlılığı, termodinamiğin en temel ve verimli matematiksel yapılarından birini oluşturmaktadır.


1. Temel Kavramsal ve Matematiksel Zemin

1.1 Gibbs Serbest Enerjisinin Tanımı ve Doğal Değişkenleri

Gibbs serbest enerjisi şu bağıntıyla tanımlanır:

G=HTS

Bu ifade, iç enerji (U), hacim (V) ve basınç (P) üzerinden entalpinin H=U+PV biçiminde yazılmasıyla genişletildiğinde:

G=U+PVTS

elde edilir. Termodinamiğin temel bağıntısından (dU=TdSPdV) yola çıkılarak G’nin tam diferansiyeli türetilir[1]:

dG=SdT+VdP

Bu denklem, G’nin doğal değişkenlerinin T ve P olduğunu gösterir. Sabit bileşim koşullarında:

(GT)P=S

(GP)T=V

Maxwell bağıntısı olarak da bilinen bu iki kısmi türev ifadesi, Gibbs fonksiyonunun sıcaklık ve basınçla nasıl değiştiğini kesin biçimde ortaya koyar. Entropi her zaman pozitif bir büyüklük olduğundan, G’nin sıcaklıkla değişimi daima negatif eğimlidir; yani sıcaklık arttıkça Gibbs serbest enerjisi azalır. Öte yandan, G’nin basınçla olan türevi molar hacme eşittir; dolayısıyla basınç artışı Gibbs enerjisini arttırır.

1.2 Sıcaklığa Bağlılık: Gibbs-Helmholtz Denklemi

Sıcaklıkla değişimi daha ayrıntılı incelemek için G/T oranının sıcaklığa göre türevi alınır. Bu işlem uygulandığında Gibbs-Helmholtz denklemi elde edilir[2]:

((G/T)T)P=HT2

Kimyasal süreçlere uygulandığında bu denklem:

((ΔG/T)T)P=ΔHT2

biçimini alır. Bu ifade büyük bir pratik değer taşır: iki farklı sıcaklıktaki ΔG değerleri arasındaki ilişki, entalpi değişimine (ΔH) bağlanarak ifade edilebilir. Sabit ΔH varsayımı altında bu denklem entegre edildiğinde:

ΔG(T2)T2=ΔG(T1)T1ΔH(1T21T1)

bağıntısına ulaşılır. Bu denklem, yalnızca bir referans sıcaklığındaki ΔG ve ΔH verileri bilindiğinde, farklı sıcaklıklardaki serbest enerji değerlerinin hesaplanmasını mümkün kılar.

Sıcaklığa Bağlı Kendiliğindenlik Koşulları

ΔG=ΔHTΔS ifadesinden hareketle, bir sürecin kendiliğindenliği sıcaklığa bağlı olarak dört temel kategoride değerlendirilir[3]:

ΔH ΔS ΔG Yorumlama
Negatif Pozitif Her T’de negatif Her koşulda kendiliğinden
Pozitif Negatif Her T’de pozitif Hiçbir koşulda kendiliğinden değil
Negatif Negatif Düşük T’de negatif Düşük sıcaklıklarda kendiliğinden
Pozitif Pozitif Yüksek T’de negatif Yüksek sıcaklıklarda kendiliğinden

Kendiliğindenliğin sıcaklığa göre değiştiği eşik noktası ΔG=0 koşulundan:

Teq=ΔHΔS

olarak belirlenir. Bu sıcaklıkta sistem denge halindedir.

1.3 Basınca Bağlılık: İdeal Gazlar

Sabit sıcaklıkta dG=VdP ifadesinin entegrasyonu, başlangıç basıncı P’dan son basınç P’ye kadar yapılır. İdeal gaz için V=nRT/P olduğundan:

G(T,P)=G(T)+nRTlnPP

Bu bağıntı, ideal bir gazın molar Gibbs enerjisinin basınçla logaritmik olarak arttığını gösterir[4]. Logaritmik bağımlılığın kökeninde tamamen entropik etkiler yatmaktadır: basınç artıkça gazın hacmi azalır, erişilebilir mikro-durum sayısı düşer ve entropi değeri küçülür. Azalan entropi, G=HTS bağıntısı aracılığıyla Gibbs enerjisini artırır.

Tek bir mol gaz için molar Gibbs enerjisi:

G¯(T,P)=G¯(T)+RTlnPP

biçiminde yazılır; burada G¯(T) standart koşullardaki (P=1 bar) molar Gibbs enerjisidir.

1.4 Basınca Bağlılık: Gerçek Gazlar ve Kaçınıklık (Fugasity) Kavramı

Gerçek gazlar için moleküller arası çekim ve itme kuvvetleri, ideal gaz varsayımından sapmalara yol açar. Bu sapmalar, Gilbert N. Lewis tarafından 1901’de önerilen kaçınıklık (f) kavramıyla ele alınır[5]. Kaçınıklık, gerçek gazın etkin basıncını temsil eden bir büyüklüktür ve şu bağıntıyla tanımlanır:

G¯(T,P)=G¯(T)+RTlnff

Kaçınıklık katsayısı ϕ=f/P ise gerçek gazın ideal davranıştan sapma ölçüsünü verir. Düşük basınçlarda ϕ1 olup kaçınıklık basınca eşit hale gelir.

Gerçek bir gaz için kaçınıklık katsayısı, durum denklemleri (örneğin viral açılım) kullanılarak hesaplanır:

lnϕ=0P(Z1P)dP

burada Z=PV/nRT sıkışabilirlik faktörüdür. Azot gazı (N2) için 0°C ve 100 atm koşulunda kaçınıklık katsayısı ϕ=0,9703 olarak ölçülmüştür; bu durum, gerçek azot gazının Gibbs enerjisinin ideal değerden RTln(0,9703) kadar farklılaştığını ortaya koyar[6].

1.5 Katı ve Sıvılar İçin Basınca Bağlılık

Yoğun fazlar (sıvı ve katılar) için molar hacim yaklaşık olarak basınçtan bağımsızdır. Bu nedenle:

G¯(T,P)=G¯(T)+V¯(PP)

Sıkıştırılamaz madde yaklaşımı olarak da bilinen bu ifade, katı ve sıvıların molar hacminin (V¯) küçüklüğü nedeniyle basıncın Gibbs enerjisi üzerindeki etkisinin gazlara kıyasla çok daha sınırlı olduğunu gösterir. Ancak çok yüksek basınç koşullarında (jeofiziksel derinlik gibi) bu etki ihmal edilemez hale gelir.


2. Faz Dengesi ve Clapeyron Denklemi

2.1 Faz Geçişlerinde Gibbs Enerjisi

İki faz dengedeyken her fazın molar Gibbs enerjisi birbirine eşittir:

G1(T,P)=G2(T,P)

Bu koşul, faz diyagramlarındaki denge eğrilerinin matematiksel temelini oluşturur[7]. Denge eğrisi boyunca sıcaklık artışına karşılık basınçtaki değişim, Clapeyron denklemiyle verilir:

dPdT=ΔSΔV=ΔHTΔV

Bu denklem tamamen geneldir; katı-sıvı, sıvı-gaz ve katı-gaz dengelerine uygulanabilir[8]. Burada ΔS ve ΔV, sırasıyla faz geçişindeki entropi ve hacim değişmelerini gösterir; ΔH ise faz geçişinin molar entalpisidir (gizli ısı).

2.2 Clausius-Clapeyron Denklemi

Sıvı-gaz dengesi için Ayrıştırılamadı (sözdizim hatası): {\textstyle V_{gaz} \gg V_{sıvı}} ve gaz fazının ideal davranış gösterdiği varsayıldığında Clapeyron denkleminden Clausius-Clapeyron denklemi elde edilir[9]:

dlnPdT=ΔHvapRT2

İki sıcaklık arasında entegre edildiğinde:

lnP2P1=ΔHvapR(1T21T1)

Bu bağıntı, buhar basıncının sıcaklıkla üstel olarak arttığını, büyüme hızının ise buharlaşma entalpisi tarafından belirlendiğini ortaya koyar.

Somut örnek — Su: ΔHvap=40,7 kJ/mol ve T1=373 K’de P1=1,00 atm bilindiğinde, T2=383 K’deki buhar basıncı yaklaşık 1,41 atm, T2=363 K’deki ise yaklaşık 0,70 atm olarak hesaplanır. Bu hesap, 10°C’lik bir sıcaklık artışının buhar basıncını neredeyse iki katına çıkardığını gösterir.

2.3 Suyun Katı-Sıvı Dengesi: Negatif Eğim Durumu

Çoğu maddenin aksine, su katılaşırken hacmi artmaktadır (Ayrıştırılamadı (sözdizim hatası): {\textstyle \Delta V_{katılaşma} < 0} ). Bu nedenle:

dPdT=ΔHerimeTΔV<0

Buzun erimesi için tipik değerler kullanıldığında (ΔSerime22 J mol1 K1, ΔVerime1,7×106 m3/mol) eğim yaklaşık 1300 bar/K olarak hesaplanır[10]. Bu son derece dik ve negatif eğim, suyun erime noktasının basınçla çok yavaş düştüğünü; yalnızca 1°C’lik bir erime noktası düşüşü için yaklaşık 130 bar’lık bir basınç artışı gerektiğini gösterir. Aynı katılar için bu değer yaklaşık 55 bar/K düzeyindedir.


3. Güncel Araştırmalardan Bulgular

3.1 Gibbs-Helmholtz Denkleminin Endüstriyel Uygulamaları

Gibbs-Helmholtz denkleminin CO2 tutma sistemlerindeki uygulamaları, son yıllarda yoğun biçimde araştırılmıştır. Amin bazlı CO2 absorpsiyon sistemlerinde buhar-sıvı denge verileriyle kalorimetre ölçümleri arasındaki termodinamik tutarlılığın sınanmasında Gibbs-Helmholtz bağıntısı kilit bir araç olarak kullanılmaktadır[11]. Bu çerçevede gerçekleştirilen araştırmalarda, Gibbs-Helmholtz testinin denge verilerindeki sistematik hataları, basit Gibbs-Duhem tutarlılık testlerinden çok daha etkili biçimde tespit edebildiği ortaya konulmuştur[12].

3.2 Makine Öğrenmesiyle Gibbs Enerjisi Tahmini

2023 yılında Digital Discovery dergisinde yayımlanan bir çalışmada, Gibbs-Helmholtz Graph Neural Network (GH-GNN) adı verilen yeni bir makine öğrenimi modeli geliştirilmiştir[13]. Bu model, Gibbs-Helmholtz denkleminden türetilmiş ifadelerle grafik sinir ağlarını birleştirerek farklı sıcaklıklardaki sonsuz seyreltme aktivite katsayılarını (γ) yüksek doğrulukla tahmin etmektedir. Model; dispersiyon etkileşimleri ve hidrojen bağı oluşumu için açık moleküler betimleyiciler içermekte olup ayrışım süreçlerinin tasarımında kritik rol oynayan bu termodinamik büyüklüklerin geniş bir kimyasal uzayda tahmin edilmesini mümkün kılmıştır.

3.3 Metabolik Ağlarda Sıcaklığa Bağlı Gibbs Enerjisi

Grup-katkı yöntemleri kullanılarak metabolik tepkimelerin denge sabitlerinin sıcaklıkla değişiminin modellenmesi, biyokimyasal termodinamik alanında önemli bir gelişme olarak değerlendirilmektedir[14]. Sözü edilen çalışmada 600’den fazla enzimatik tepkime için denge sabitleri, entalpi ve entropi değişimleri derlenerek sıcaklığa bağlı ΔG hesaplamalarını mümkün kılan kapsamlı bir veritabanı oluşturulmuştur. Termodinamik açıdan uygulanabilir metabolik yolların belirlenmesinde Gibbs-Helmholtz bağıntısı merkezi bir rol üstlenmiştir.

3.4 Yüksek Basınçta Faz Dengesi Araştırmaları

Yüksek basınç koşullarında gerçekleştirilen son deneysel çalışmalar, Clapeyron eğrisi eğiminin özellikle faz sınırına yakın bölgelerde ve kritik nokta civarında karmaşık bir davranış sergilediğini göstermiştir. Jeofiziksel sistemlerde mantonun derin katmanlarına ait yüksek basınç-sıcaklık koşullarında mineral dönüşümlerinin belirlenmesi için Clapeyron denklemi uygulanmakta; 55 bar/K mertebesindeki dik eğim değerleri, önemli sıcaklık değişimlerine yol açmak için gereken muazzam basıncın altını çizmektedir[15].

3.5 Gibbs-Helmholtz Denkleminin Matematiksel Basitleştirilmesi

Journal of Chemical Education’da yayımlanan bir çalışma, Gibbs enerjisinin ve denge sabitlerinin sıcaklığa bağımlılığını hesaplamak için geleneksel Gibbs-Helmholtz çözümüne alternatif, basit bir cebirsel yöntem sunmuştur[16]. Polinom ısı kapasitesi ifadelerini kapsayan daha karmaşık durumlarda bile bu yöntemin yüksek doğrulukla sonuç verdiği gösterilmiş; yöntemin literatür değerleriyle karşılaştırılması olumlu sonuçlar vermiştir.


4. Kavramsal Analiz

4.1 Nizam, Gaye ve Sanat Analizi

Gibbs serbest enerjisinin basınç ve sıcaklıkla olan bağımlılığı, termodinamiğin en titiz matematiksel yapılarından birini barındırmaktadır. dG=SdT+VdP denkleminin tek bir özgün büyüklükte entalpi, entropi ve hacim bilgisini bir araya getirmesi; dahası, bu ilişkinin doğal değişkenler aracılığıyla tam diferansiyel bir yapıya oturması dikkat çekicidir.

G’nin sıcaklığa bağlılığı entropi aracılığıyla, basınca bağlılığı ise hacim aracılığıyla tanımlandığında ortaya çıkan simetrik yapı, termodinamik Maxwell bağıntıları ile de iç içe geçer. Bu bağıntıların birbirleriyle uyumlu olması, farklı termodinamik süreçler arasındaki kesin ve öngörülebilir ilişkileri mümkün kılar. Clapeyron denklemi bu uyumun bir ürünüdür: faz geçişinin eğimini (dP/dT), entalpi ve hacim değişimlerine bağlayan bu ifade, özünde iki koşulun eşzamanlı sağlandığı bir denge geometrisini tanımlar.

Suyun katı-sıvı faz sınırındaki negatif eğim durumu özellikle çarpıcıdır: artan basıncın erimeyi kolaylaştırdığı bu istisnai durum, ΔV<0 olan tek yaygın maddeyi temsil eder. Bu termodinamik özgünlüğün biyolojik süreçler ve iklim döngüleri için taşıdığı işlevsel önem—suyun 0°C’nin biraz üzerinde maksimum yoğunluğa ulaşması, buzun katman halinde yüzmesi—son derece dikkat çekicidir. Böylesine özelleşmiş termodinamik parametreler kombinasyonunun yaşam için bu denli kritik koşulları sağlaması düşündürücüdür.

4.2 İndirgemeci ve Materyalist Safsataların Eleştirisi

Gibbs serbest enerjisine ilişkin popüler anlatılar çoğu zaman bir kavramsal yanılgıyı barındırır: “Sistemler Gibbs enerjisini minimize eder” ya da “Moleküller en düşük enerji durumunu seçer” gibi ifadeler, soyut bir matematiksel büyüklüğe aktif bir fail rolü yükler. Oysa ΔG<0 koşulunun kendisi fiziksel bir nedenselllik taşımaz; bu koşul yalnızca gözlemlenen davranışın betimsel bir özetidir.

Benzer bir kavramsal hata, Gibbs-Helmholtz denkleminin açıklamalarında da karşılaşılmaktadır. “Denklem sisteme geleceğini söyler” türünden ifadeler, matematiksel ilişkiyi bir karar mekanizması gibi konumlandırır. Gerçekte denklem, belirlenen T ve P değerlerinde gözlemlenen denge durumunu sayısal olarak tanımlar; herhangi bir sistemi yönlendirmez. Clapeyron denklemi de faz sınırının eğimini betimler; faz geçişini başlatmaz. Betimleyici ile kurucu arasındaki bu ayrımın göz ardı edilmesi, termodinamik yasaların gerçek açıklama gücünü belirsizleştirir.

Kaçınıklık kavramı da benzer bir eleştirinin konusudur. “Gerçek gazlar ideal davranıştan kaçar” ifadesi kasıt içermez; gözlemlenen basınç-hacim sapmaları, moleküller arası kuvvetlerin yarattığı fiili etkileşimlerin sonucudur. Kaçınıklık bu sapmayı betimleyen matematiksel bir araçtır, sapmaya neden olan bir mekanizma değildir.

4.3 Hammadde ve Sanat Ayrımı Analizi

Gibbs serbest enerjisi, bileşen büyüklüklerinden hiçbirinde mevcut olmayan bir özelliği temsil eder. Enthalpi (H) ısıl içeriği, entropi (S) sistemin düzensizliğini ölçer; ancak G=HTS bağıntısındaki T ile çarpım ve çıkarma işlemi, ayrı ayrı ölçülebilir olmayan yeni bir işlevsellik ortaya çıkarır: belirli koşullar altında gerçekleştirilecek maksimum faydalı iş kapasitesini.

Bu durum Clapeyron denkleminde somut bir anlam kazanır. Denge eğrisinin eğimi (dP/dT=ΔH/TΔV) ne tek başına entalpiden, ne entropiden ne de hacim değişiminden elde edilebilir; bu büyüklüklerin sabit sıcaklık ve basınç koşulunda birbiriyle kombinasyonundan ortaya çıkan bir özelliktir. Daha çarpıcı bir örnek, suyun faz diyagramının üçlü noktasıdır: 273{,}16 K ve 611{,}73 Pa değerlerinde katı, sıvı ve gaz fazlarının eşzamanlı dengede bulunduğu bu nokta, üç ayrı Gibbs-enerji yüzeyinin kesişiminden meydana gelir. Tek tek fazların özelliklerinden bu kesişim noktasının koordinatlarını türetmek mümkün değildir.

İdeal gazdaki RTln(P/P) basınç terimi de aynı özelliği taşır: bu büyüklük, gaz sabitinin (R), mutlak sıcaklığın (T) ve standart duruma göre logaritmik basınç oranının birlikte oluşturduğu bir ifadedir. Her birim başına R=8,314 J mol1 K1 değeri sayısal bir sabit olsa da bu sabiti barındıran RTln(P/P) terimi, bireysel moleküllerin taşımadığı kolektif bir termodinamik büyüklüğü—molar bazda yapılabilecek maksimum geri dönüşümlü çalışı—temsil eder.


5. Sonuç

Gibbs serbest enerjisinin sıcaklık ve basınca bağlılığı, termodinamiğin en verimli ve uyumlu yapılarından birini barındırır. dG=SdT+VdP ifadesi; entalpi, entropi ve hacmi tek bir diferansiyel denklemde birleştirir. Bu ifadeden türetilen Gibbs-Helmholtz denklemi, sıcaklığın denge sabitleri ve serbest enerji üzerindeki etkisini ortaya koyarken; G=G+RTlnP bağıntısı, basıncın gazlar üzerindeki logaritmik etkisini tanımlar. Gerçek gazlardaki sapmalar kaçınıklık (f) kavramıyla sistematik biçimde ele alınır ve her durumda Gibbs enerjisi, kimyasal sistemin potansiyelini tek bir sayısal büyüklükte ifade eder.

Clapeyron ve Clausius-Clapeyron denklemleri, faz geçişinin eğimini, gizli ısı ve hacim değişimiyle ilişkilendirerek faz diyagramlarının matematiksel temelini oluşturur. Suyun katı-sıvı faz sınırındaki negatif eğim, bu genel şemanın içindeki olağandışı bir özel durum olarak öne çıkar. Makine öğrenimi temelli Gibbs-Helmholtz modellerinden jeofiziksel faz dengesi çalışmalarına uzanan geniş bir araştırma çerçevesinde, bu ilişkilerin deneysel ve hesaplamalı geçerliliği sürekli biçimde sınanmakta ve genişletilmektedir.

Entalpi ile entropinin sıcaklık katsayısıyla çarpımının birbirini giderek dengelediği ya da pekiştirdiği koşullar, faz geçişlerini, kimyasal dengeleri ve biyokimyasal süreçleri düzenleyen termodinamik yapıyı oluşturur. Bu yapının matematiksel tutarlılığı, iç uyumu ve pek çok farklı ölçekte geçerliliği—atomik düzeyden jeolojik derinliklere kadar—deliller ışığında değerlendirildiğinde, nihai karar okuyucunun kendi aklına ve vicdanına bırakılmaktadır.


Kaynakça

  1. Atkins, P., & de Paula, J. (2014). Atkins’ Physical Chemistry (10th ed.). Oxford University Press.
  2. Ellgen, P. (2020). Thermodynamics and Chemical Equilibrium. Creative Commons Attribution ShareAlike 4.0 License. https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Thermodynamics_and_Chemical_Equilibrium_(Ellgen)/10:_Some_Mathematical_Consequences_of_the_Fundamental_Equation/10.11:_The_Gibbs-Helmholtz_Equation
  3. Middelburg, J. J. (2024). The Gibbs free energy. In Thermodynamics and Equilibria in Earth System Sciences: An Introduction. SpringerBriefs in Earth System Sciences. Springer. https://doi.org/10.1007/978-3-031-53407-2_4
  4. Peverati, R. (2023). The Live Textbook of Physical Chemistry 1. Florida Tech. https://peverati.github.io/pchem1/GibbsFree.html
  5. Lewis, G. N. (1908). The osmotic pressure of concentrated solutions, and the laws of the perfect solution. Journal of the American Chemical Society, 30(5), 668–683.
  6. Wikipedia contributors. (2025). Fugacity. In Wikipedia, The Free Encyclopedia. https://en.wikipedia.org/wiki/Fugacity
  7. Tatum, J. B. (2021). Heat and Thermodynamics. LibreTexts. https://phys.libretexts.org/Bookshelves/Thermodynamics_and_Statistical_Mechanics/Heat_and_Thermodynamics_(Tatum)/14:_The_Clausius-Clapeyron_Equation
  8. Wikipedia contributors. (2025). Clausius–Clapeyron relation. In Wikipedia, The Free Encyclopedia. https://en.wikipedia.org/wiki/Clausius%E2%80%93Clapeyron_relation
  9. Levine, I. N. (2009). Physical Chemistry (6th ed.). McGraw-Hill.
  10. Middelburg, J. J. (2024). The Gibbs free energy. SpringerBriefs in Earth System Sciences. Springer. https://doi.org/10.1007/978-3-031-53407-2_4
  11. Mathias, P. M., & Kister, H. Z. (2013). The Gibbs–Helmholtz equation and the thermodynamic consistency of chemical absorption data. Industrial & Engineering Chemistry Research, 52(23), 7807–7817. https://doi.org/10.1021/ie202668k
  12. Mathias, P. M. (2016). The Gibbs–Helmholtz equation in chemical process technology. Industrial & Engineering Chemistry Research, 55(4), 1076–1087. https://doi.org/10.1021/acs.iecr.5b03405
  13. Sanchez Medina, A. I., Rehner, P., & Theophil, A. (2023). Gibbs–Helmholtz graph neural network: Capturing the temperature dependency of activity coefficients at infinite dilution. Digital Discovery, 2, 781–798. https://doi.org/10.1039/D2DD00142J
  14. Noor, E., Haraldsdóttir, H. S., Milo, R., & Fleming, R. M. T. (2013). Consistent estimation of Gibbs energy using component contributions. PLOS Computational Biology, 9(7), e1003098. https://doi.org/10.1371/journal.pcbi.1003098 (güncel grup-katkı metodolojisinin temeli; PMC6129446 numaralı çalışmayla güncellenmiştir)
  15. Middelburg, J. J. (2024). The Gibbs free energy. SpringerBriefs in Earth System Sciences. Springer. https://doi.org/10.1007/978-3-031-53407-2_4
  16. Elgeti, K. (2014). A simple method to calculate the temperature dependence of the Gibbs energy and chemical equilibrium constants. Journal of Chemical Education, 91(3), 396–401. https://doi.org/10.1021/ed3006874