Kinetische analysemethoden bij thermische analyse zijn rekenkundige procedures waarmee uit thermische meetdata — doorgaans afkomstig van thermogravimetrische analyse (TGA) of differentiële scanning calorimetrie (DSC) — kwantitatieve kinetiekparameters worden afgeleid: de activeringsenergie (Ea), de pre-exponentiële factor (A) en het reactiemodel f(α). Met deze parameters wordt de snelheid waarmee een thermisch geactiveerde reactie verloopt beschreven als functie van temperatuur en omzettingsgraad. De toepassingen lopen uiteen van levensduurvoorspelling van polymeren en farmaceutische stabiliteitsmodellering tot procesveiligheidsbeoordeling en optimalisatie van pyrolyse- en sintertechnologie.
De omzettingssnelheid dα/dt van een thermisch geactiveerde reactie wordt beschreven door de combinatie van de Arrhenius-vergelijking met een reactiemodelfunctie f(α):
dα/dt = A · exp(−Ea/RT) · f(α)
Hierin is α de omzettingsgraad (dimensieloos, van 0 tot 1), A de pre-exponentiële factor (s⁻¹), Ea de activeringsenergie (J/mol), R de universele gasconstante (8,314 J/mol·K) en T de absolute temperatuur (K). Bij niet-isothermische metingen (dynamische TGA) geldt T = T₀ + βt, waarbij β de verwarmingssnelheid in K/min is. Na deling door β ontstaat de niet-isothermische kinetiekvergelijking:
dα/dT = (A/β) · exp(−Ea/RT) · f(α)
De drie grootheden Ea, A en f(α) worden gezamenlijk aangeduid als de kinetische driehoek of het kinetische triplet. Een volledige kinetische karakterisering levert alle drie; voor levensduurberekeningen volstaat in de praktijk de activeringsenergie als primaire parameter.
De omzettingsgraad α wordt bij TGA gedefinieerd als het aandeel van het totale verwachte massaverlies dat al heeft plaatsgevonden:
α = (m₀ − m) / (m₀ − m∞)
Hierin is m₀ de beginmassa, m de actuele massa en m∞ de eindmassa. Bij DSC-metingen wordt α bepaald als de gedeeltelijke piekoppervlakte op het totale piekoppervlak.
Het diagram toont het isoconversieprincipe: bij dezelfde omzettingsgraad α = 0,5 wordt de karakteristieke temperatuur Tα bij vier verwarmingssnelheden afgelezen (linkerpaneel). Deze waarden worden uitgezet in de Ozawa-Flynn-Wall-plot (rechterpaneel), waarbij de helling van de regressielijn direct Ea oplevert.
Model-afhankelijke methoden gaan ervan uit dat het reactiemodel f(α) vooraf bekend is. Het meest gebruikte model in niet-isothermische kinetiek is de n-de-orde reactie: f(α) = (1−α)n. Andere veelgebruikte modellen zijn de Avrami-Erofeev-nucleatiemodellen (f(α) = n(1−α)[−ln(1−α)](n-1)/n) voor kristallisatie en sinteren, en diffusiemodellen (D1–D4) voor vaste-stofreacties waarbij diffusie de limiterende stap is.
De Kissinger-methode is de meest gebruikte model-afhankelijke benadering voor niet-isothermische kinetiek. Zij maakt gebruik van de DTG-piektemperatuur Tp bij verschillende verwarmingssnelheden β. De Kissinger-vergelijking luidt:
ln(β/Tp2) = −Ea/(R · Tp) + ln(A · R / Ea)
Door ln(β/Tp2) uit te zetten tegen 1/Tp ontstaat een rechte lijn met helling −Ea/R. Voordelen van de Kissinger-methode zijn de eenvoud en de directe afleesbaarheid van Ea en A uit één lineaire regressie. De beperking is dat de methode strikt alleen geldig is voor eerste-orde reacties (n = 1), al wordt zij in de praktijk ook voor andere orde gebruikt als benadering. Conform de ICTAC Kinetics Committee-aanbevelingen (Vyazovkin et al., 2011) dient de Kissinger-methode te worden gecombineerd met isoconversiemethoden voor een volledige kinetische evaluatie.
Model-vrije methoden — ook isoconversiemethoden of iso-conversional methods — maken geen aanname over het reactiemodel f(α). Zij bepalen Ea als functie van de omzettingsgraad α uit metingen bij meerdere verwarmingssnelheden (of temperaturen bij isothermische methoden). Dit maakt het mogelijk een activeringsenergieprofiel Ea(α) op te stellen, waarmee complexe meerstaps- of parallelle reacties worden geïdentificeerd: een constant Ea(α) wijst op een enkelvoudige reactie, een veranderend Ea(α) op een samengesteld mechanisme.
De Flynn-Wall-Ozawa-methode (FWO) is de meest toegepaste lineaire isoconversiemethode. De Doyle-benadering van de temperatuurintegraal leidt tot de FWO-vergelijking:
ln(β) = ln(A · Ea / R · g(α)) − 5,331 − 1,052 · Ea / (R · Tα)
Bij vaste α wordt ln(β) uitgezet tegen 1/Tα; de helling bedraagt −1,052 Ea/R. Uit de helling van de regressielijn volgt direct Ea. De methode is van toepassing op TGA- en DSC-data, vereist minimaal drie (bij voorkeur vier tot vijf) verwarmingssnelheden, en is genormaliseerd in ASTM E1641 en ISO 11358-2.
De Friedman-methode is een differentiële isoconversiemethode die werkt zonder de Doyle-benadering en daardoor nauwkeuriger is, maar gevoeliger voor ruis in de meetsignalen. De Friedman-vergelijking:
ln(dα/dt) = ln(A · f(α)) − Ea / (R · T)
Bij constante α worden de instantane omzettingssnelheid dα/dt en de corresponderende temperatuur T bij elke verwarmingssnelheid afgelezen. Lineaire regressie van ln(dα/dt) versus 1/T geeft de helling −Ea/R. De Friedman-methode is nauwkeuriger dan FWO maar vereist differentiatie van de experimentele data, wat ruisgevoeliger is. Zij is geïmplementeerd in alle gangbare thermische analysesoftware (NETZSCH Proteus, Mettler-Toledo STARe, TA Instruments Trios).
De Vyazovkin-methode is een niet-lineaire isoconversiemethode die geen benadering van de temperatuurintegraal gebruikt en daardoor de meest nauwkeurige Ea(α)-profielen oplevert. De methode minimaliseert numeriek de doelstelling:
Φ(Ea) = ∑i ∑j≠i [ I(Ea, Tα,i) · βj ] / [ I(Ea, Tα,j) · βi ]
waarbij de temperatuurintegraal I(Ea, Tα) numeriek wordt berekend. De Vyazovkin-methode is de aanbevolen methode van de ICTAC Kinetics Committee voor een nauwkeurig kinetisch triplet en is genormaliseerd in ISO 11358-3. De methode is beschikbaar in gespecialiseerde kinetische analysesoftware en als plug-in voor TGA-softwarepakketten.
Bij isothermische kinetiekanalyse worden metingen uitgevoerd bij meerdere vaste temperaturen Ti in plaats van bij meerdere verwarmingssnelheden. De tijd tα,i waarop omzettingsgraad α wordt bereikt bij temperatuur Ti vormt de experimentele invoer. De isothermische versie van de isoconversievergelijking:
ln(tα) = ln(g(α)/A) + Ea / (R · T)
Door ln(tα) uit te zetten tegen 1/T ontstaat een rechte lijn met helling Ea/R. Isothermische methoden hebben als voordeel dat zij geen verwarmingssnelheidseffecten kennen en thermodynamisch dichter bij evenwicht werken; als nadeel geldt dat het instellen van de gewenste temperatuur tijd kost en dat bij hoge temperaturen de reactie al begint vóór de isotherm is bereikt (het zogenoemde non-zero conversion at the start probleem). Isothermische kinetiek is de voorkeursmethode bij langetermijn-levensduurstudies (polymeerveroudering, voedselstabiliteit) conform ASTM E2550.
Bij model-afhankelijke analyse en bij de afsluitende modelidentificatie na isoconversieanalyse wordt een reactiemodel gekozen. De meest gebruikte modellen zijn:
Modelidentificatie geschiedt door de geïntegreerde functie g(α) = ∫₀α dα/f(α) te berekenen en lineaire regressie van g(α) versus de temperatuurintegraal uit te voeren. Het model met de hoogste correlatiecoëfficiënt R² en de best overeenkomende Ea met de isoconversieanalyse wordt geselecteerd.
De International Confederation for Thermal Analysis and Calorimetry (ICTAC) publiceerde in 2011 uitgebreide aanbevelingen voor thermische kinetiekanalyse (Vyazovkin et al., Thermochimica Acta 520, 1–19). De kernpunten zijn: gebruik altijd minimaal drie verwarmingssnelheden (aanbevolen: vier tot vijf); gebruik bij voorkeur isoconversiemethoden boven model-afhankelijke methoden als eerste stap; rapporteer Ea(α) als functie van α in plaats van één enkel Ea-getal; valideer het kinetische triplet door voorspelling versus experiment te vergelijken; gebruik de Friedman- of Vyazovkin-methode voor de meest nauwkeurige resultaten. De relevante normen zijn:
Kinetische analyse van thermische oxidatieve ontleding van polymeren is de meest toegepaste toepassing van TGA-kinetiek. Via multi-rate TGA in lucht worden Ea en A bepaald. Met het Arrhenius-model wordt de tijd tot een vastgesteld percentageverlies (bijv. 5% massaverlies) geëxtrapoleerd naar gebruikstemperatuur. Voor een polyamide met Ea = 180 kJ/mol kan de levensduur bij 100 °C worden berekend uit metingen bij 250–400 °C: de extrapolatie beslaat drie tot vier ordes van grootte in de tijdschaal. Conform ISO 11358-2 wordt de levensduur uitgedrukt als tend(T) = g(α)/(A · exp(−Ea/RT)).
In de farmaceutische industrie wordt thermische kinetiek gebruikt voor de houdbaarheidsmodellering van werkzame stoffen (API's) en formulaties. Gecombineerde TGA/DSC-meting — zie ook het artikel over simultane thermische analyse (STA) — bij meerdere verwarmingssnelheden geeft de ontledingsonset-temperatuur en het kinetisch triplet. Met de Arrhenius-extrapolatie naar opslagtemperatuur (25 °C, 40 °C) wordt de chemische stabiliteit onder ICH Q1A-condities voorspeld. Polymorfvormen en solvaten vertonen karakteristieke kinetische profielen, waardoor TGA-kinetiek ook wordt ingezet voor polymorfscreening.
Voor chemische procesvoering is de bepaling van de zelf-versnellende ontledingstemperatuur (SADT) en de maximale temperatuurstijging bij ongecontroleerde ontleding (adiabatische opschaling) een kritische veiligheidstaak. DSC en TGA geven de ontledingsenthalpie (ΔHdec) en de ontledingsonset-temperatuur; kinetische analyse levert de temperatuurafhankelijkheid van de ontledingssnelheid. De adiabatische opschaling van laboratoriumschaal (mg) naar productieschaal (ton) vereist correctie voor de zogenoemde phi-factor (Φ). Materialen met Ea < 80 kJ/mol en een hoge ΔHdec vereisen bijzondere procesveiligheidsmaatregelen en adiabatische calorimetrische validatie. Voor de kinetiek van de gewenste synthese-reactie zelf — onder realistische procescondities met dosering, oplosmiddel en menging — is reactiecalorimetrie de aangewezen techniek.
Biomassapyrolyse verloopt via meerdere parallelle en sequentiële reacties (cellulose, hemicellulose, lignine). Kinetische analyse via Friedman of Vyazovkin geeft het Ea(α)-profiel dat de bijdragen van de afzonderlijke componenten onderscheidt. Ea voor celluloseontleding ligt typisch in het bereik van 150–200 kJ/mol; ligninepyrolyse verloopt over een breed temperatuurbereik met Ea van 40–200 kJ/mol als functie van α. Kinetische parameters worden gebruikt voor de ontwerpsimulatie van pyrolyse-reactoren en de optimalisatie van biochar-productie.
Sinterings- en kristallisatiekinetiek van keramische poeders, glas en katalysatoren worden bestudeerd via dilatometrische data (zie dilatometrie) of via DSC-kristallisatiepieken. Avrami-Erofeev-modellen zijn bij uitstek geschikt voor deze systemen. De piektemperatuurverschuiving als functie van de verwarmingssnelheid geeft via de Kissinger-methode direct Ea voor kristallisatienucleatie en -groei.
Het Chuzhakin-rekensysteem (ook aangeduid als de Chuzhakin-methode of Chuzhakin-Borovinskaya-methode) is een modelidentificatieprocedure voor vaste-stofreacties waarbij door middel van lineaire transformaties van TGA-data het meest passende reactiemodel f(α) wordt geïdentificeerd. De methode construeert een genormaliseerd mastercurveplot: experimentele g(α)/g(0,5)-waarden worden uitgezet tegen α en vergeleken met de theoretische mastercurven van een bibliotheek van twintig tot veertig reactiemodellen. Het model waarvan de mastercurve het dichtst aansluit bij de experimentele data wordt geselecteerd. Het Chuzhakin-systeem is een verfijning van de klassieke Coats-Redfern-methode en is met name gangbaar in de Oost-Europese wetenschappelijke literatuur. In Westeuropese en Amerikaanse laboratoria worden doorgaans de ICTAC-aanbevolen isoconversiemethoden (Friedman, Vyazovkin) gecombineerd met een aansluitende modelidentificatiestap als volwaardige vervanging gebruikt.
De activeringsenergie Ea wordt in thermische kinetiekanalyse bepaald via de helling van een lineaire regressie waarbij een logaritmische grootheid wordt uitgezet tegen de reciproke temperatuur 1/T. De drie meest gebruikte formules zijn:
In alle gevallen geldt: Ea = −R × helling (gecorrigeerd voor de methodefactor bij FWO). De gasconstante R = 8,314 J/(mol·K). Een Ea van 100 kJ/mol geeft een helling van −12.029 K bij Kissinger en −12.651 K bij FWO. Praktisch: bij een verwarmingssnelheid van 10 °C/min en Ea = 150 kJ/mol verschuift de piektemperatuur Tp met ongeveer 8–12 K per verdubbeling van de verwarmingssnelheid — dit is de experimenteel waarneembare vuistregel die bevestigt dat de reactie temperatuurgeactiveerd is.
Een correcte kinetische analyse vereist een zorgvuldige experimentele opzet. De ICTAC-aanbevelingen schrijven voor: gebruik minimaal drie verwarmingssnelheden, bij voorkeur 2, 5, 10 en 20 °C/min voor een breed temperatuurbereik; houd de monstermassa klein (1–10 mg) om temperatuurgradiënten en diffusie-effecten in het monster te minimaliseren; gebruik een inerte atmosfeer (stikstof of argon, 50 ml/min) tenzij de reactie oxidatief is; voer een buoyantiecorrectie uit (blinde meting met lege pan); en kalibreer temperatuur en massa voor elke meetreeks conform ISO 11358. Voor DSC-kinetiek gelden vergelijkbare eisen: kleine monsters (2–10 mg), hermetisch gesloten kroezen voor vluchtige monsters, en verwarmingssnelheden van 1–40 °C/min. Raadpleeg ook het artikel over simultane thermische analyse (STA) wanneer TGA en DSC gelijktijdig worden ingezet voor kinetiekstudies — STA elimineert monsteronhomogeniteit als foutbron.
Model-vrije methoden (Ozawa-Flynn-Wall, Friedman, Vyazovkin) bepalen Ea(α) zonder aanname van een reactiemodel. Zij zijn geschikt als eerste stap en geven inzicht in de reactiecomplexiteit: een constant Ea(α) wijst op een enkelvoudige reactie; een variërend Ea(α) op een meerstapsmechanisme. Model-afhankelijke methoden (Kissinger, Coats-Redfern) gaan uit van een vooraf gekozen f(α) en zijn geschikt als het mechanisme bekend is of na een isoconversieanalyse als aanvullende stap voor de volledige kinetische tripletbepaling.
ICTAC beveelt minimaal drie verwarmingssnelheden aan; vier tot vijf geven robuustere resultaten. De snelheden dienen voldoende van elkaar te verschillen — een factor twee tot tien tussen opeenvolgende snelheden is gebruikelijk (bijv. 2, 5, 10, 20 °C/min). Te dicht bij elkaar gelegen snelheden geven een slechte regressie door de kleine piektemperatuurverschuiving; te ver uit elkaar kan de isoconversieaanname schenden bij complexe reacties met een veranderend mechanisme als functie van de verwarmingssnelheid.
Ja — alle isoconversiemethoden zijn direct toepasbaar op DSC-pieken: α wordt dan bepaald als de gedeeltelijke piekoppervlakte. DSC-kinetiek is bij uitstek geschikt voor kristallisatie (Avrami-analyse), uithardingsreacties van harsen (Kissinger voor exotherme pieken) en smeltkinetiek. TGA is beter geschikt voor ontledingsreacties waarbij massaverlies de primaire kwantificeerbare grootheid is. Voor reacties waarbij zowel massaverlies als een warmte-effect optreedt, geeft gecombineerde STA-kinetiek (gelijktijdige TGA en DSC op hetzelfde monster) de meest consistente resultaten, zoals beschreven in het artikel over simultane thermische analyse (STA).
Een stijgend Ea(α) bij toenemende omzettingsgraad duidt op een reactie waarbij de gemakkelijkst-reagerende fractie als eerste reageert (distributies in reactiviteit, bijv. in biomassa of polymeermengsels). Een dalend Ea(α) wijst op autocatalyse: de reactie wordt versneld door zijn eigen producten. Een Ea(α)-profiel met een minimum in het middengebied (α ≈ 0,3–0,7) is kenmerkend voor een parallelle reactie waarbij twee processen met verschillende Ea gelijktijdig verlopen. In alle gevallen is een variërend Ea(α)-profiel een signaal dat een enkelvoudig Ea-getal onvoldoende is voor een correcte levensduurvoorspelling of procesmodellering.
Kinetische analysemethoden worden vrijwel altijd toegepast op data van thermogravimetrische analyse (TGA), de primaire techniek voor massaverliesmeting als functie van temperatuur. Voor gelijktijdige meting van massaverlies en warmtestroom — wat de kinetiekanalyse aanzienlijk versterkt doordat beide tegelijkertijd beschikbaar zijn — is simultane thermische analyse (STA) de aangewezen combinatietechniek. Identificatie van de gasvormige ontledingsproducten die de kinetisch bestudeerde reactie genereren, wordt uitgevoerd via Evolved Gas Analysis (EGA): TGA-MS en TGA-FTIR. Voor de kinetiek van gewenste synthesereacties onder realistische procescondities is reactiecalorimetrie de complementaire methode. Sintering- en kristallisatiekinetiek worden ook bestudeerd via dilatometrische data, zie dilatometrie. Voor de calorimetrische karakterisering van de ontledingsenthalpie bij constante druk is differentiële scanning calorimetrie (DSC) de basistechniek, terwijl voor de verbrandingswarmte bij constant volume bomcalorimetrie wordt ingezet.
Disclaimer: De informatie in dit artikel is met zorg samengesteld op basis van wetenschappelijke bronnen en ICTAC-normering. Canidae Seal B.V. / Labvakhandel.nl is niet aansprakelijk voor de toepassing van kinetische analysemethoden in specifieke industriële of veiligheidsgerelateerde situaties. Raadpleeg voor regulatoire toepassingen de relevante normen (ISO 11358, ASTM E1641) en betrek een gekwalificeerde thermal analysis specialist.
Inloggen
Wachtwoord vergeten
Account aanmaken
Uw winkelwagen is leeg.
