Lipid Oxidation / Rancidity Kinetics

1 Introduksi

1.1 Latar Belakang

Oksidasi lipid adalah penyebab utama rancidity (ketengikan) pada minyak, lemak, dan produk pangan fatty selama penyimpanan, melalui mekanisme radical chain autoxidation (initiation → propagation → termination) yang menghasilkan hidroperoksida primer lalu carbonyl sekunder (aldehida, keton) yang flavor-active dan mungkin toxic (Frankel, 2005; Choe & Min, 2006). Pada level industri, oksidasi mengurangi shelf life produk fatty, menurunkan nilai nutrisi (kehilangan vitamin E + omega-3), dan menghasilkan off-flavor yang ditolak konsumen.

Kuantifikasi oksidasi memakai tiga metrik klasik: Peroxide Value (PV) untuk hidroperoksida primer (mEq O₂/kg, AOCS Cd 8-53), TBARS untuk aldehida sekunder (mg MDA/kg, malondialdehyde), dan Anisidine Value (AV) untuk α/β-unsaturated aldehyde (AOCS Cd 18-90). Kinetika oksidasi tipikal memiliki induction period + propagation + plateau — sigmoid shape. Akselerasi via Rancimat (110 °C, air flow) memberikan induction time $t_{\text{ind}}$ yang Arrhenius-extrapolable ke kondisi storage normal (Frankel, 2005; Farhoosh, 2007).

1.2 Tujuan Modul

Modul Lipid Oxidation / Rancidity Kinetics di SQalytics ditujukan untuk:

Memfit kinetika autocatalytic atau sigmoid untuk PV/TBARS vs time data.
Mendeteksi induction time $t_{\text{ind}}$ via tangent intersection (Farhoosh method).
Mendukung multi-temperature untuk Arrhenius fit dengan $Q_{10}$ output.
Memberikan shelf life prediction terhadap PV threshold (mis. 10 mEq/kg AOCS limit edible oils).
Audiens: mahasiswa S2 ilmu pangan/chemistry yang menjalankan shelf life study minyak, R&D antioksidan, QC industri minyak/snack.

1.3 Posisi di Antara Alternatif

Pilih Lipid Oxidation / Rancidity Kinetics untuk fit kinetika autocatalytic + induction time. Untuk Maillard non-enzymatic browning, gunakan Maillard / Browning Kinetics (paralel kinetic untuk warna). Untuk shelf life multi-quality (browning + oxidation + flavor bersama), pakai ASLT / Shelf-Life Prediction. Untuk inaktivasi mikroba termal, pakai D-value / Z-value / F-value. Untuk kuantifikasi PV tunggal dengan kalibrasi, pakai Standard Curve.

2 Metode

2.1 Dasar Teoretis

Autocatalytic kinetics (Özilgen & Özilgen, 1990) — diferensial:

$$ \frac{d[\text{PV}]}{dt} = k \cdot [\text{PV}] \cdot ([\text{PV}]_\infty - [\text{PV}]) $$

dengan rate proportional ke konsentrasi PV existing (autocatalysis) dan distance dari asymptote. Solusi:

$$ \boxed{\, [\text{PV}](t) = \frac{[\text{PV}]_\infty}{1 + e^{-k[\text{PV}]_\infty (t - t_{1/2})}} \,} $$

dengan $t_{1/2}$ inflection point. Sigmoid shape: induction period (lag) + rapid propagation + plateau.

Logistic model alternatif untuk PV/TBARS sigmoid:

$$ [\text{PV}](t) = \frac{a}{1 + \exp(-b(t - c))} $$

dengan $a$ = asymptote, $b$ = rate, $c$ = inflection time.

Induction time $t_{\text{ind}}$ — Farhoosh method (Farhoosh, 2007): perpotongan tangent dari dua fase (lag + propagation) pada kurva PV vs time. Praktis: $t_{\text{ind}}$ = waktu di mana PV mulai naik signifikan dari baseline (mis. PV > 2 × initial PV).

Rancimat (Frankel, 2005) — accelerated oxidation pada 110 °C dengan air flow. Induction time Rancimat $t_R$ correlate dengan shelf life normal via Arrhenius:

$$ \ln t_R = \ln A_R - \frac{E_a}{R T} $$

Extrapolasi ke ambient: $t_{\text{shelf, 25°C}} = t_R \cdot Q_{10}^{(T_R - 25)/10}$ dengan $T_R$ Rancimat suhu (110 °C) dan $Q_{10}$ Rancimat 2.0–2.5 tipikal.

Peroxide Value calculation (AOCS Cd 8-53):

$$ \text{PV} (\text{mEq O}_2/\text{kg}) = \frac{(V_s - V_b) \cdot N \cdot 1000}{m} $$

dengan $V_s$ titer thiosulfate sampel (mL), $V_b$ blank, $N$ normalitas Na₂S₂O₃ (0.01 N tipikal), $m$ massa sampel (g).

TBARS (AOCS Cd 19-90):

$$ \text{TBARS} (\text{mg MDA/kg}) = A_{532} \cdot \frac{M_{\text{MDA}}}{\varepsilon \cdot l} \cdot \frac{V_{\text{dilution}}}{m} $$

dengan $\varepsilon = 1.56 \times 10^5$ L·mol⁻¹·cm⁻¹ untuk MDA-TBA adduct.

Shelf life ($\text{PV}_{\text{end}}$ threshold):

$$ t_{\text{shelf}} = c + \frac{1}{b} \ln\!\left(\frac{a - \text{PV}_{\text{end}}}{\text{PV}_{\text{end}}}\right) \quad \text{(logistic)} $$

AOCS recommended limit: edible refined oils PV ≤ 10 mEq/kg, virgin olive oil PV ≤ 20 mEq/kg, crude oils PV ≤ 5 mEq/kg.

2.2 Persamaan Inti

Autocatalytic: $[\text{PV}](t) = [\text{PV}]_\infty / (1 + e^{-k[\text{PV}]_\infty (t - t_{1/2})})$

Logistic: $[\text{PV}](t) = a / (1 + e^{-b(t-c)})$

Induction time: $t_{\text{ind}}$ = inflection tangent intersection

PV (AOCS): $(V_s - V_b) \cdot N \cdot 1000 / m$

TBARS (AOCS): $A_{532} \cdot M_{\text{MDA}} / (\varepsilon \cdot l) \cdot V/m$

Rancimat extrapolation: $t_{\text{shelf, ambient}} = t_R \cdot Q_{10}^{(T_R - T_{\text{ambient}})/10}$

2.3 Asumsi & Batas Validitas

Asumsi	Konsekuensi jika dilanggar	Cara cek di SQalytics
Sampel disimpan dengan O₂ access (atau closed-system terkontrol)	Kinetika berbeda	Catat packaging (closed/open)
Suhu konstan selama storage interval	Arrhenius fit invalid	Catat $T$ real-time
Cukup titik pada lag + propagation phase	$t_{\text{ind}}$ tidak teridentifikasi	Sample weekly + bi-weekly
Tidak ada hydrolytic rancidity overlap	TBARS bias	Cek FFA secara terpisah
Antioksidan natural di sampel tidak diabaikan	$t_{\text{ind}}$ extended	Catat tocopherol, polyphenol
Rancimat 110 °C representatif	Extrapolasi bias bila mekanisme berubah	Pakai 80–120 °C multi-level

3 Cara Kerja

3.1 Step-by-Step di SQalytics

Buka Lipid Oxidation / Rancidity Kinetics dari domain Kimia.
Muat tabel: Time_day, PV_mEq_per_kg (atau TBARS), opsional Temperature_C.
Pilih model:
- Logistic sigmoid (default untuk PV/TBARS).
- Autocatalytic.
- Zero-order (induction-free, untuk crude oil).
(Opsional) Aktifkan induction time detection (Farhoosh tangent method).
(Multi-suhu) Aktifkan Arrhenius mode + input $T$ dalam Kelvin.
Atur target PV end untuk shelf life (default 10 mEq/kg, AOCS edible oils).
Klik Run Oxidation Kinetics.
Tinjau hasil:
- Tab Fit Plot — sigmoid + tangent lines + $t_{\text{ind}}$ marker.
- Tab Kinetic Parameters — $a$, $b$, $c$, $t_{\text{ind}}$, $R^2$.
- Tab Arrhenius / $Q_{10}$ — extrapolation Rancimat → ambient.
- Tab Shelf Life Prediction — $t_{\text{shelf}}$ + 95% CI.

Setelah hasil utama tampil di Step 4, gunakan Step 5 opsional untuk mengirim Oxidation curves atau Induction review ke Publication Graph Studio (PGS). Di sana Anda bisa lanjut mengatur typography, legend, preview Print, lalu ekspor SVG/PNG/PDF untuk kebutuhan publikasi.

3.2 Template Tabel Input + Contoh Data Sintetis

Kolom	Tipe	Wajib	Catatan
`Time_day`	numeric	✓	Storage time
`PV_mEq_per_kg`	numeric	✓	Peroxide Value (atau TBARS)
`Temperature_C`	numeric	◯	Untuk Arrhenius

Contoh data sintetis (palm olein RBD storage 60 °C accelerated):

Time_day	PV_mEq_per_kg
0	0.5
3	0.8
6	1.2
10	2.5
14	5.8
18	12.4
22	18.5
26	22.8
30	24.5
35	25.0

Data sintetis — palm olein RBD pada storage accelerated 60 °C dengan O₂ access. Lag period 0–10 days, rapid propagation 10–22 days, plateau setelah 25 days. Induction time $\approx$ 10 days, $t_{\text{shelf}}$ to PV = 10 ≈ 16 days at 60 °C → ekstrapolasi ambient 25 °C $\approx$ 5–6 bulan dengan $Q_{10}$ 2.3.

3.3 Contoh Luaran

Logistic fit:

Parameter	Value	Interpretasi
$a$ (asymptote)	25.2 mEq/kg	Maximum PV achievable
$b$ (rate)	0.38 day⁻¹	Steepness sigmoid
$c$ (inflection $t_{1/2}$)	18.5 day	Midpoint kurva
$t_{\text{ind}}$ (Farhoosh)	10.2 day	Induction time
$R^2$	0.998	Excellent fit
RMSE	0.42 mEq/kg	Residual error kecil

Shelf life prediction (target PV = 10 mEq/kg, AOCS edible oil limit):

Storage T (°C)	$t_{\text{shelf}}$ (day)	Konversi (bulan)
60 (accelerated)	16.0	0.53
35 (warm)	$\sim$ 38 (Q₁₀ 2.3)	1.27
25 (ambient)	$\sim$ 88	2.93
4 (chilled)	$\sim$ 470	15.6

Kesimpulan ringkas: "Palm olein RBD pada accelerated 60 °C menunjukkan classical lipid oxidation sigmoid: lag 0–10 hari (induction $t_{\text{ind}}$ = 10.2 hari), propagation 10–22 hari, plateau pada PV ≈ 25 mEq/kg. Logistic fit excellent ($R^2 = 0.998$). Shelf life ke AOCS limit PV = 10 mEq/kg: pada 60 °C = 16 hari, pada ambient 25 °C ≈ 88 hari (≈ 3 bulan) via Arrhenius $Q_{10}$ = 2.3 extrapolation. Konsisten dengan label palm olein komersial 6 bulan dengan margin keamanan packaging O₂-barrier. Rekomendasi: tambahkan antioksidan TBHQ atau tocopherol untuk extend $t_{\text{ind}}$; pakai metallic chelator (citric acid) untuk minyak unrefined dengan trace metal. Lanjut ke Compare Many Groups untuk uji formulasi antioksidan."

Kinetika oksidasi lipid palm olein logistic + induction time + AOCS threshold — Gambar 1. Kinetika oksidasi lipid palm olein RBD pada storage accelerated 60 °C dengan akses O₂ (n=1 per timepoint, data sintetis). Peroxide value (PV) mengikuti model logistic dengan *baseline* ($pv_0$ = 0.5 mEq/kg fresh oil): PV($t$) = 0.5 + 24.7/(1 + exp(–0.38($t$ – 18.5))), $R^2$ = 0.998. Tiga fitur kunci ditandai: induction time $t_{ind}$ (Farhoosh) = 10.2 day (titik kelengkungan maksimum fase awal — ditandai bintang hijau), $t_{shelf}$ pada threshold AOCS PV = 10 mEq/kg = 17.3 day @ 60 °C (ditandai X merah), dan asimtot $a$ = 25.2 mEq/kg sebagai PV maksimum yang dapat dicapai. Ekstrapolasi Arrhenius ke *ambient* 25 °C dengan $Q_{10} \approx 2.3$ memberi shelf life ≈ 5–6 bulan.

4 Kesimpulan

4.1 Relevansi Real-World

Minyak goreng komersial — PV monitoring untuk label expired by.
Snack berbasis minyak (keripik, kacang) — TBARS untuk shelf life prediction.
Daging beku — TBARS untuk freezer burn quality marker.
Susu bubuk — oksidasi lipid susu (oxidative milk powder) selama 24 bulan.
Omega-3 supplement — induction time + PV untuk QC fish oil capsule.
Antioksidan natural (rosemary, ekstrak Singkil) — uji efficacy via extension induction time.
Pakan ternak — TBARS untuk feed quality kontrol.

4.2 Where to Go from Here

Maillard / Browning Kinetics — paralel kinetik untuk warna.
ASLT / Shelf-Life Prediction — multi-quality framework.
Compare Many Groups — uji statistik antar formulasi antioksidan.
AOCS Official Methods Cd 8-53, Cd 18-90, Cd 19-90 untuk standar metode.

⚙ Troubleshooting Cepat

PV plateau sangat awal (< 5 mEq/kg). Sample mungkin tidak access O₂ — buka kemasan atau pakai shaker.

TBARS tinggi pada Day 0. Indikasi sample sudah teroksidasi sebelumnya — gunakan fresh material.

Logistic fit gagal. Coba autocatalytic atau zero-order; cek apakah data masih dalam lag phase.

$t_{\text{ind}}$ tidak terdeteksi. Tambah titik di lag region (mis. weekly 0–14 days).

Arrhenius $Q_{10}$ > 4 atau < 1.5. Mekanisme berubah dengan suhu; rentang lebih sempit untuk Arrhenius valid.

i Riwayat Revisi

Tanggal	Revisi	Penulis
2026-05-12	Migrasi MD v2 → HTML final dengan figure publikasi + caption Elsevier-style (W1 batch malam 12 Mei)	Claude
2026-05-12	Draft v2 publikasi (KaTeX autocatalytic/logistic/Rancimat + AOCS PV/TBARS + APA Frankel/Choe/Farhoosh)	Claude

4 Referensi

Frankel, E. N. (2005). Lipid oxidation (2nd ed.). The Oily Press. https://doi.org/10.1533/9780857097927
Choe, E., & Min, D. B. (2006). Mechanisms and factors for edible oil oxidation. Comprehensive Reviews in Food Science and Food Safety, 5(4), 169–186. https://doi.org/10.1111/j.1541-4337.2006.00009.x
Farhoosh, R. (2007). The effect of operational parameters of the Rancimat method on the determination of the oxidative stability measures and shelf-life prediction of soybean oil. Journal of the American Oil Chemists' Society, 84(3), 205–209. https://doi.org/10.1007/s11746-006-1030-4
Özilgen, S., & Özilgen, M. (1990). Kinetic model of lipid oxidation in foods. Journal of Food Science, 55(2), 498–501. https://doi.org/10.1111/j.1365-2621.1990.tb06797.x
AOCS. (2017). Official methods and recommended practices of the AOCS (7th ed., Methods Cd 8-53 PV, Cd 18-90 AV, Cd 19-90 TBARS). American Oil Chemists' Society.
Shahidi, F., & Zhong, Y. (2010). Lipid oxidation and improving the oxidative stability. Chemical Society Reviews, 39(11), 4067–4079. https://doi.org/10.1039/B922183M