Terapan Matematika dalam Kesetimbangan Kimia SMA dan OSN

Dalam pokok bahasan kesetimbangan kimia, matematika bukan sekadar alat hitung, melainkan bahasa untuk merumuskan hubungan antarspesi dan memprediksi arah pergeseran reaksi.

Berikut lima terapan matematika utamanya, lengkap dengan contoh soal nyata dan penyelesaiannya.

1

Aljabar & Persamaan Kuadrat

Aljabar

Tabel ICE (Initial, Change, Equilibrium) sering menghasilkan persamaan polinomial. Jika nilai perubahan $x$ tidak dapat diabaikan, rumus kuadrat dipakai untuk mencari $x$ secara eksak.

Contoh Soal

Reaksi $\mathrm{N_2O_4(g) \rightleftharpoons 2NO_2(g)}$ memiliki $K_c = 0,21$ pada 350 K. Mula-mula hanya ada $\mathrm{N_2O_4}$ dengan konsentrasi $0.100\ \mathrm{M}$. Hitung konsentrasi tiap spesi saat setimbang.

Penyelesaian

Dari tabel ICE: $[\mathrm{N_2O_4}] = 0,100 - x$ dan $[\mathrm{NO_2}] = 2x$.

$$ \begin{aligned} K_c &= \frac{(2x)^2}{0,100 - x} \\ 0.21 &= \frac{4x^2}{0,100 - x} \\ 4x^2 &= 0,21(0,100 - x) \\ 4x^2 + 0,21x - 0,021 &= 0 \\ x &= \frac{-0,21 + \sqrt{(0,21)^2 - 4(4)(-0,021)}}{2(4)} \\ x &= \frac{-0,21 + \sqrt{0,0441 + 0,336}}{8} \\ x &= \frac{-0,21 + \sqrt{0,3801}}{8} \\ x &= \frac{-0,21 + 0,6165}{8} \\ x &\approx 0,0508\ \mathrm{M} \end{aligned} $$

Hasil

$[\mathrm{N_2O_4}] = 0,100 - 0,0508 = 0,0492\ \mathrm{M}$ dan $[\mathrm{NO_2}] = 2(0,0508) = 0,1016\ \mathrm{M}$.

Verifikasi: $(0,1016)^2 / 0,0492 \approx 0,21$ ✓

2

Logaritma & Fungsi Eksponensial

Logaritma

Hubungan energi bebas Gibbs dan tetapan kesetimbangan dirumuskan secara logaritmik: $\Delta G^\circ = -RT\ln K$. Bentuk ini memungkinkan perhitungan $K$ pada berbagai suhu.

Contoh Soal

Suatu reaksi memiliki $\Delta G^\circ = -32,9\ \mathrm{kJ/mol}$ pada 298 K. Tentukan nilai tetapan kesetimbangannya.

Penyelesaian

$$ \begin{aligned} \Delta G^\circ &= -RT\ln K \\ \ln K &= -\frac{\Delta G^\circ}{RT} \\ \ln K &= -\frac{-32\,900}{(8,314)(298)} \\ \ln K &= \frac{32\,900}{2\,477,6} \\ \ln K &\approx 13,28 \\ K &= e^{13,28} \\ K &\approx 5,85 \times 10^{5} \end{aligned} $$

Hasil

$K \approx 5,9 \times 10^{5}$, reaksi berlangsung sangat jauh ke kanan (produk dominan).

3

Kalkulus (Turunan & Optimasi)

Kalkulus

Keadaan setimbang tercapai saat energi bebas Gibbs minimum, yaitu $\dfrac{dG}{d\xi}=0$. Dari syarat ini pula lahir relasi $\Delta G^\circ=-RT\ln K$. Secara praktis, turunan juga memberi sensitivitas $K$ terhadap suhu.

Contoh Soal

Suatu reaksi endotermik memiliki $\Delta H^\circ = +58\ \mathrm{kJ/mol}$. Perkirakan berapa persen $K$ berubah untuk setiap kenaikan 1 K di sekitar suhu 500 K.

Penyelesaian

Gunakan bentuk diferensial persamaan van't Hoff:

$$ \begin{aligned} \frac{d\ln K}{dT} &= \frac{\Delta H^\circ}{RT^2} \\ \frac{dK/K}{dT} &= \frac{58\,000}{(8,314)(500)^2} \\ &= \frac{58\,000}{2\,078\,500} \\ &\approx 2,79 \times 10^{-2}\ \mathrm{K^{-1}} \end{aligned} $$

Hasil

$K$ naik sekitar 2,79% per kenaikan 1 K di sekitar 500 K. Inilah wujud matematis dari "reaksi endotermik bergeser ke kanan saat suhu dinaikkan".

4

Grafik & Metode Numerik

Grafik / Numerik

Plot van't Hoff ($\ln K$ terhadap $1/T$) memiliki kemiringan $-\Delta H^\circ/R$, sehingga $\Delta H^\circ$ bisa ditentukan dari data dua suhu. Untuk sistem rumit yang tak bisa diselesaikan analitik, dipakai metode iterasi seperti Newton–Raphson.

Contoh Soal

Dari percobaan diperoleh $K = 1,0\times10^{-3}$ pada 500 K dan $K = 5,0\times10^{-2}$ pada 700 K. Tentukan $\Delta H^\circ$ reaksi.

Penyelesaian

$$ \begin{aligned} \ln\frac{K_2}{K_1} &= -\frac{\Delta H^\circ}{R}\left(\frac{1}{T_2} - \frac{1}{T_1}\right) \\ \ln\left(\frac{5,0\times10^{-2}}{1,0\times10^{-3}}\right) &= -\frac{\Delta H^\circ}{8,314}\left(\frac{1}{700} - \frac{1}{500}\right) \\ \ln(50) &= -\frac{\Delta H^\circ}{8,314}\left(1,4286\times10^{-3} - 2,0\times10^{-3}\right) \\ 3,912 &= -\frac{\Delta H^\circ}{8,314}\left(-5,714\times10^{-4}\right) \\ 3,912 &= \Delta H^\circ \times 6,873\times10^{-5} \\ \Delta H^\circ &= \frac{3,912}{6,873\times10^{-5}} \\ \Delta H^\circ &\approx 5,69\times10^{4}\ \mathrm{J/mol} \end{aligned} $$

Hasil

$\Delta H^\circ \approx +56,9\ \mathrm{kJ/mol}$, reaksi bersifat endotermik, konsisten dengan $K$ yang naik saat suhu ditingkatkan.

5

Perbandingan & Aritmetika Dasar

Aritmetika

Perhitungan $K_c$ dan $K_p$ melibatkan operasi pangkat dan akar sesuai koefisien stoikiometri, serta konversi satuan lewat hubungan $K_p = K_c(RT)^{\Delta n}$. Neraca massa dan muatan pada titrasi juga berupa sistem aritmetika yang harus dipenuhi bersamaan.

Contoh Soal

Untuk reaksi $\mathrm{N_2(g) + 3H_2(g) \rightleftharpoons 2NH_3(g)}$ diketahui $K_c = 0,060$ pada 500 K. Hitung $K_p$ reaksi tersebut ($R = 0,0821\ \mathrm{L\,atm\,mol^{-1}K^{-1}}$).

Penyelesaian

Hitung dulu perubahan jumlah mol gas:

$$ \begin{aligned} \Delta n &= 2 - (1 + 3) = -2 \\ K_p &= K_c (RT)^{\Delta n} \\ K_p &= 0,060 \times (0.0821 \times 500)^{-2} \\ &= 0,060 \times (41,05)^{-2} \\ &= \frac{0,060}{1685,1} \\ &\approx 3,56 \times 10^{-5} \end{aligned} $$

Hasil

$K_p \approx 3,6 \times 10^{-5}$. Nilainya jauh lebih kecil dari $K_c$ karena $\Delta n < 0$ (jumlah mol gas berkurang).

6

Fungsi Eksponensial, Asimtot & Waktu Relaksasi

Bentuk Kurva

Kurva konsentrasi–waktu pada reaksi reversibel tidak berhenti tajam di garis setimbang, melainkan mendekati secara asimtotik. Selisih konsentrasi dari nilai setimbang meluruh mengikuti fungsi eksponensial negatif, dengan tetapan peluruhan sama dengan jumlah tetapan laju maju dan balik.

[A] pereaksi [B] hasil garis asimtot

Contoh Soal

Reaksi reversibel orde satu \(A \rightleftharpoons B\) memiliki \(k_+ = 0,20\ \mathrm{s^{-1}}\) dan \(k_- = 0,05\ \mathrm{s^{-1}}\). Mula-mula \([A]_0 = 1,00\ \mathrm{M}\) dan \([B]_0 = 0\). Tentukan \([A]_{\text{eq}}\), waktu relaksasi \(\tau\), dan \([A]\) pada \(t = 5\ \mathrm{s}\).

Penyelesaian

Konsentrasi setimbang \(A\) untuk sistem ini adalah:

$$ \begin{aligned} [A]_{\text{eq}} &= \frac{k_-}{k_+ + k_-}[A]_0 \\ &= \frac{0,05}{0,20 + 0,05}(1,00) \\ &= \frac{0,05}{0,25}(1,00) \\ &= 0,200\ \mathrm{M} \end{aligned} $$

Waktu relaksasi adalah kebalikan dari total tetapan laju:

$$ \begin{aligned} \tau &= \frac{1}{k_+ + k_-} \\ &= \frac{1}{0,20 + 0,05} \\ &= \frac{1}{0,25} \\ &= 4,00\ \mathrm{s} \end{aligned} $$

Untuk \(t = 5\ \mathrm{s}\), gunakan persamaan peluruhan eksponensial:

$$ \begin{aligned} [A] - [A]_{\text{eq}} &= \left([A]_0 - [A]_{\text{eq}}\right)e^{-(k_+ + k_-)t} \\ [A] - 0,200 &= (1,00 - 0,200)\,e^{-(0,25)(5)} \\ [A] - 0,200 &= 0,800\,e^{-1,25} \\ [A] - 0,200 &= 0,800 \times 0,2865 \\ [A] - 0,200 &= 0,2292 \\ [A] &= 0,4292\ \mathrm{M} \end{aligned} $$

Hasil

\([A]_{\text{eq}} = 0,200\ \mathrm{M}\), \(\tau = 4,00\ \mathrm{s}\), dan pada \(t = 5\ \mathrm{s}\) diperoleh \([A] \approx 0,429\ \mathrm{M}\). Setelah satu \(\tau\) (4 s), deviasi dari kesetimbangan sudah menyusut menjadi sekitar \(37\%\) dari deviasi awal, inilah “sidik jari” matematis dari pendekatan asimtotik pada kurva kesetimbangan.

Kurva \([A]\) turun melengkung mendekati \(0,200\ \mathrm{M}\) sebagai asimtot horizontal; kemiringannya makin landai seiring waktu. Plot \(\ln([A]-[A]_{\text{eq}})\) terhadap \(t\) akan menghasilkan garis lurus dengan kemiringan \(-(k_+ + k_-)\), sehingga parameter kinetik bisa diekstrak lewat regresi linear.