Water Activity / Sorption Isotherm

1 Introduksi

1.1 Latar Belakang

Aktivitas air ($a_w$) didefinisikan sebagai rasio tekanan uap air di sampel terhadap tekanan uap air murni pada suhu yang sama: $a_w = p / p_0$. $a_w$ adalah prediktor utama stabilitas pangan karena mengontrol pertumbuhan mikroba (bakteri patogen butuh $a_w > 0.86$; jamur > 0.65), reaksi enzimatik, oksidasi lipid, dan Maillard browning (Labuza, 1968; Rahman, 2009). Hubungan antara $a_w$ dan kadar air kesetimbangan ($X$, g·g⁻¹ basis kering) pada suhu tertentu disebut kurva isoterm sorpsi dan berbentuk sigmoid Tipe II (Brunauer, Emmett, & Teller, 1938) untuk sebagian besar pangan dengan komponen makromolekul.

Dua model paling sering dipakai untuk fitting: BET (valid pada rentang $a_w$ 0.05–0.45) dan GAB (Guggenheim-Anderson-de Boer, valid 0.05–0.90) (van den Berg & Bruin, 1981). Parameter kritis adalah monolayer moisture content ($X_m$) — kadar air pada lapisan tunggal molekul air terikat di situs polar permukaan padatan. Pada $X = X_m$, laju oksidasi lipid berada di minimum global (Karel & Lund, 2003) sehingga $X_m$ menjadi target moisture untuk penyimpanan jangka panjang produk kering (susu bubuk, kopi instan, snack ekstrudat). Penentuan $X_m$ yang akurat adalah inti dari penetapan shelf life.

1.2 Tujuan Modul

Modul Water Activity / Sorption Isotherm di SQalytics ditujukan untuk:

Membaca tabel pasangan ($a_w$, moisture) dan memvisualisasi kurva sorpsi sigmoid.
Memfitkan BET dan GAB secara otomatis dengan output parameter $X_m$, $C$, dan $K$.
Membandingkan goodness-of-fit (RMSE, $R^2$, MRPE) antar dua model untuk memilih yang terbaik.
Memberikan estimasi target moisture penyimpanan ($X_m$) untuk diskusi shelf life.
Audiens: mahasiswa S2 teknologi pangan yang menjalankan studi shelf life, peneliti R&D produk kering, dan QC industri pangan dehidrasi.

1.3 Posisi di Antara Alternatif

Pilih Water Activity / Sorption Isotherm ketika data berisi pasangan ($a_w$, moisture) dari chamber humidity-controlled atau aw-meter + oven moisture. Untuk profil $a_w$ tunggal (monitoring storage), pakai Quick Stats Check. Untuk kinetika moisture loss selama drying (Page, Newton, Midilli), gunakan modul drying kinetics. Untuk kalkulasi nilai $a_w$ dari larutan garam standar, gunakan Dilution & Unit Calculator. Untuk uji statistik $X_m$ antar produk, lanjut ke Compare 2 Groups.

2 Metode

2.1 Dasar Teoretis

Model BET (Brunauer-Emmett-Teller, 1938) berasal dari teori adsorpsi gas multilayer, diadaptasi untuk sorpsi air:

$$ \frac{a_w}{(1 - a_w) X} = \frac{1}{X_m C} + \frac{C - 1}{X_m C} \cdot a_w $$

Persamaan ini linear dalam $a_w$ bila ruas kiri di-plot terhadap $a_w$. Slope dan intercept langsung memberi $X_m$ dan $C$. BET hanya valid pada $a_w < 0.45$ karena di atas itu kondensasi kapilar dan struktur multilayer mendominasi (Rahman, 2009, hal. 421).

Model GAB (Guggenheim 1966, Anderson 1946, de Boer 1953) memperluas BET dengan parameter ketiga $K$ untuk mengakomodasi multilayer pada $a_w$ tinggi:

$$ X = \frac{X_m \cdot C \cdot K \cdot a_w}{(1 - K a_w)(1 - K a_w + C K a_w)} $$

GAB valid pada $a_w$ 0.05–0.90 dan memberikan fit lebih baik untuk pangan kompleks (van den Berg & Bruin, 1981). Bila $K = 1$, GAB tereduksi ke BET. Tipikal $K = 0.7$–$0.9$ untuk pangan kering.

Parameter fisik utama:

$X_m$ — monolayer moisture (g·g⁻¹ db, target storage).
$C$ — konstanta Guggenheim, energi sorpsi monolayer ($C > 2$: sigmoid Tipe II).
$K$ (GAB only) — faktor koreksi multilayer ($K < 1$ untuk fit fisik valid).

Goodness-of-fit:

$$ \text{RMSE} = \sqrt{\frac{1}{n} \sum_{i=1}^{n} (X_i^{\text{obs}} - X_i^{\text{pred}})^2} $$ $$ R^2 = 1 - \frac{\sum (X_i^{\text{obs}} - X_i^{\text{pred}})^2}{\sum (X_i^{\text{obs}} - \bar{X}^{\text{obs}})^2} $$

Target praktis untuk model sorpsi yang dapat dipublikasikan: $R^2 > 0.99$ dan mean relative percentage error $\text{MRPE} < 10\%$ (Lomauro, Bakshi, & Labuza, 1985):

$$ \text{MRPE} = \frac{100}{n} \sum_{i=1}^{n} \left| \frac{X_i^{\text{obs}} - X_i^{\text{pred}}}{X_i^{\text{obs}}} \right| $$

2.2 Persamaan Inti

BET linear:

$$ \frac{a_w}{(1 - a_w) X} = \frac{1}{X_m C} + \frac{C - 1}{X_m C} a_w $$

GAB nonlinear:

$$ X = \frac{X_m C K a_w}{(1 - K a_w)(1 - K a_w + C K a_w)} $$

Monolayer dari BET:

$$ X_m = \frac{1}{\text{slope} + \text{intercept}} $$

Mean relative percentage error:

$$ \text{MRPE} = \frac{100}{n} \sum_{i=1}^{n} \left| \frac{X_i^{\text{obs}} - X_i^{\text{pred}}}{X_i^{\text{obs}}} \right| $$

2.3 Asumsi & Batas Validitas

Asumsi	Konsekuensi jika dilanggar	Cara cek di SQalytics
Kesetimbangan tercapai (≥ 7 hari di chamber)	Fit miring, $X_m$ underestimate	Tinjau raw — moisture stabil antar replikasi
Suhu konstan selama isoterm	Pergeseran kurva (efek termodinamika)	Catat suhu di metadata; gunakan satu suhu per fit
Tidak ada histeresis adsorpsi-desorpsi	Fit gabungan menghasilkan bias	Pisahkan dataset adsorpsi vs desorpsi
Titik $a_w$ tersebar (≥ 5–7 titik pada 0.1–0.9)	Fit GAB tidak konvergen	Modul fail-safe → fallback ke BET-only
Sampel tidak rusak (denaturasi, glass transition)	Isoterm sigmoid berubah bentuk	Cek visual kurva — bentuk Tipe II sigmoid?
Massa kering akurat (oven 105 °C)	Basis kering bias → semua $X$ off	Replikasi penentuan moisture awal

3 Cara Kerja

3.1 Step-by-Step di SQalytics

Buka Water Activity / Sorption Isotherm dari domain Kimia di sidebar.
Pada langkah Check the current aw table, pastikan tabel sorpsi aktif dikenali.
Pada langkah Confirm the aw and moisture columns, sesuaikan kolom aw (0 < $a_w$ < 1) dan Moisture (g·g⁻¹ basis kering atau %).
Pada langkah Choose the model strategy and run, pilih: BET only, GAB only, atau Compare BET vs GAB (rekomendasi).
Klik Run Water Activity / Sorption Isotherm.
Tinjau hasil: kesimpulan ringkas → kartu aw points/Best model/$X_m$ → tab Isotherm Graph → tab Model Table.

3.2 Template Tabel Input + Contoh Data Sintetis

Skema tabel input minimum:

Kolom	Tipe	Satuan	Wajib	Catatan
`aw`	numeric	0–1	✓	Aktivitas air (rasio tekanan uap)
`Moisture_eq`	numeric	g·g⁻¹ db atau %	✓	Kadar air kesetimbangan

Contoh data sintetis (9 baris — isoterm adsorpsi susu bubuk skim, 25 °C, basis kering):

aw	Moisture_eq (g·g⁻¹ db)
0.10	0.0250
0.20	0.0410
0.32	0.0570
0.43	0.0720
0.55	0.0925
0.65	0.1180
0.75	0.1580
0.85	0.2350
0.92	0.3450

SYNTHETIC Sorpsi adsorpsi susu bubuk skim pada 25 °C, ekuilibrasi 14 hari di chamber larutan garam jenuh (LiCl 0.11, MgCl₂ 0.33, K₂CO₃ 0.43, NaBr 0.58, NaCl 0.75, KCl 0.85; Greenspan, 1977). Bentuk Tipe II sigmoid khas pangan kaya laktosa+protein. Fit GAB memberi $X_m$ ≈ 0.052 g·g⁻¹ db — target moisture penyimpanan susu bubuk komersial. CSV setara tersedia di docs/assets/example-data/id/chemistry/template_isotherm.csv.

3.3 Contoh Luaran

Tabel Model Table:

Model	$X_m$ (g·g⁻¹ db)	$C$	$K$	$R^2$	RMSE	MRPE (%)
BET (0.05–0.45)	0.046	12.4	—	0.998	0.0021	4.1
GAB (0.05–0.92)	0.052	10.8	0.875	0.9985	0.0028	3.8

Interpretasi parameter (GAB):

Parameter	Nilai	Interpretasi
$X_m$	0.052 g·g⁻¹ db (5.2%)	Target moisture penyimpanan susu bubuk
$C$	10.8	Energi sorpsi kuat di monolayer (Tipe II)
$K$	0.875 (< 1)	Multilayer well-behaved, GAB physically valid
$R^2$	0.9985	Linearitas excellent
MRPE	3.8%	< 10% (Lomauro 1985 acceptance) ✓

Kesimpulan ringkas: "Fit isoterm sukses. GAB lebih superior dibanding BET pada rentang $a_w$ 0.10–0.92 (BET hanya valid sampai 0.45). Parameter GAB: $X_m = 0.052$ g·g⁻¹ db (5.2% basis kering), $C = 10.8$, $K = 0.875$ — semua secara fisik valid (Tipe II sigmoid, $K < 1$). $X_m = 5.2\%$ adalah target moisture penyimpanan susu bubuk untuk shelf life jangka panjang — di kadar air ini, laju oksidasi lipid dan Maillard berada di minimum global. Direkomendasikan: keringkan produk sampai ≤ 5.5% moisture, simpan di RH < 30% ($a_w < 0.30$) untuk margin keamanan. MRPE 3.8% memenuhi kriteria publikasi (< 10%, Lomauro et al. 1985)."

Grafik utama: scatter ($a_w$, $X$) dengan dua kurva fit — BET garis putus-putus merah valid sampai 0.45, GAB garis solid teal valid penuh — plus marker $X_m$ pada sumbu y (5.2%) dengan label dan garis horizontal/vertikal putus-putus.

Isoterm sorpsi susu bubuk skim 25 C dengan fit BET dan GAB plus marker Xm target storage — Gambar 1. Isoterm sorpsi adsorpsi susu bubuk skim pada 25 °C dengan perbandingan fit dua model klasik. Lingkaran biru menandakan data observasi (9 titik dari chamber larutan garam jenuh; Greenspan, 1977). **Garis putus-putus merah** adalah fit BET (Brunauer-Emmett-Teller, 1938) yang valid pada rentang transparan merah ($a_w$ 0.05–0.45, $X_m$ = 0.046, $R^2$ = 0.998). **Garis solid teal** adalah fit GAB (Guggenheim-Anderson-de Boer; van den Berg & Bruin, 1981) yang valid pada seluruh rentang ($X_m$ = 0.052, $C$ = 10.8, $K$ = 0.875, $R^2$ = 0.9985, MRPE = 3.8%). Bintang hijau pada ($a_w$ = 0.30, $X$ = 0.052) menandai **monolayer moisture target storage**: pada kadar air ini laju oksidasi lipid dan Maillard browning berada di minimum global (Karel & Lund, 2003) — rekomendasi target produk kering pada $a_w < 0.30$ untuk margin keamanan shelf life.

4 Kesimpulan

4.1 Relevansi Real-World

Water Activity / Sorption Isotherm adalah pondasi penetapan shelf life produk kering:

Susu bubuk ($X_m$ 5–6%) — target setelah spray drying, RH packaging ≤ 30%.
Kopi instan ($X_m$ 3–5%) — di atas $X_m$ caking dan agglomerasi terjadi.
Snack ekstrudat / keripik ($a_w < 0.30$) — di atas itu crispness hilang.
Buah kering / dendeng ($a_w$ 0.65–0.70) — batas atas untuk mencegah jamur xerofilik.
Bumbu kering / rempah ($a_w < 0.60$) — antara $X_m$ dan glass transition.
Farmasi padat — $a_w < 0.40$ untuk stabilitas API hidrolitik.

Dalam publikasi IMRAD, kurva isoterm + tabel parameter GAB sering muncul sebagai Figure 2 + Table 2 pada artikel drying/storage stability. Pada Rencana Publikasi Singkil v5, modul ini dipakai pada T4 Tahap 1 untuk penetapan target moisture penyimpanan ekstrak kering Singkil hasil spray drying / freeze drying.

4.2 Where to Go from Here

Modul lanjutan di SQalytics:

Compare 2 Groups — uji $X_m$ adsorpsi vs desorpsi (deteksi histeresis).
X-Y Relationship — regresi residual fit untuk diagnostik outlier.
Predict One Result — prediksi moisture pada $a_w$ target storage.
Trend Graphs — kustomisasi grafik isoterm untuk publikasi (Wong palette).
Modul drying kinetics (akan datang) — kinetika pengeringan menuju target $X_m$.

Pembacaan lanjutan:

Rahman (2009), Bab 16 — review lengkap model isoterm sorpsi.
Lomauro et al. (1985) — paper klasik perbandingan 8 model isoterm, dasar MRPE acceptance.
Brunauer et al. (1938) — paper sumber BET (recommended reading).
van den Berg & Bruin (1981) — fundamentals teori sorpsi GAB.
Greenspan (1977) — referensi nilai $a_w$ larutan garam jenuh standar.

⚙ Troubleshooting Cepat

Route gagal difit (model error). Periksa apakah kolom aw dan Moisture_eq numerik. Modul perlu ≥ 5 titik untuk BET dan ≥ 7 titik untuk GAB.

GAB fit tidak konvergen. Coba mode BET only dulu — gunakan $X_m$ BET sebagai initial guess. Pastikan rentang $a_w$ mencakup region multilayer ($a_w > 0.5$).

$K > 1$ di GAB (fisik invalid). Indikasi dataset noisy atau ada outlier. Tinjau titik $a_w$ tinggi — hilangkan replikasi yang menyimpang.

Model tampak tidak cocok ($R^2 < 0.99$). Bandingkan BET dan GAB di tabel model. Cek apakah ada efek histeresis (titik adsorpsi vs desorpsi tercampur).

Nilai $a_w$ aneh (mis. 55 vs 0.55). Pastikan satuan: $a_w$ adalah rasio 0–1, bukan persen. Modul memvalidasi $a_w < 1.0$.

$X_m$ jauh dari ekspektasi literatur. Cek (i) ekuilibrasi cukup lama (≥ 7 hari), (ii) suhu chamber stabil, (iii) basis kering vs basis basah konsisten.

i Riwayat Revisi

Tanggal	Revisi	Penulis
2026-05-12	Migrasi MD v2 → HTML final dengan figure isoterm BET vs GAB + marker $X_m$ target storage + caption Elsevier-style	Claude
2026-05-12	Migrasi v1 → v2 (template publikasi + KaTeX BET/GAB/MRPE + tabel produk komersial $X_m$ + APA)	Claude
2026-05-09	Draft awal v1	Tim docs

4 Referensi

Brunauer, S., Emmett, P. H., & Teller, E. (1938). Adsorption of gases in multimolecular layers. Journal of the American Chemical Society, 60(2), 309–319. https://doi.org/10.1021/ja01269a023
van den Berg, C., & Bruin, S. (1981). Water activity and its estimation in food systems: Theoretical aspects. In L. B. Rockland & G. F. Stewart (Eds.), Water activity: Influences on food quality (pp. 1–61). Academic Press.
Labuza, T. P. (1968). Sorption phenomena in foods. Food Technology, 22(3), 263–272.
Lomauro, C. J., Bakshi, A. S., & Labuza, T. P. (1985). Evaluation of food moisture sorption isotherm equations. Part I: Fruit, vegetable and meat products. Lebensmittel-Wissenschaft & Technologie, 18(2), 111–117.
Rahman, M. S. (Ed.). (2009). Food properties handbook (2nd ed.). CRC Press. https://doi.org/10.1201/9781420003093
Karel, M., & Lund, D. B. (2003). Physical principles of food preservation (2nd ed.). Marcel Dekker.