Siirtymätilateoria

Kemiallisen reaktion kinetiikkaa voidaan tarkastella teoreettisesti siirtymätilateorian avulla, jolloin reaktion nopeusvakio voidaan määrittää ilman empiirisiä mittauksia tai parametrejä. Tämän lisäksi siirtymätilateorialla yhdessä kokeellisesti määritettyjen nopeusvakioiden avulla voidaan laskea reaktion siirtymätilaan liittyvien termodynaamisten funktioiden arvot. Siirtymätilateoria voidaan jakaa konventionaaliseen ja termodynaamiseen lähestymistapaan.s. 159s. 177 ==Siirtymätila== Kemiallisen reaktion siirtymätila on reaktion potentiaalienergiapinnan, jota pitkin reaktiopolku (reaktiokoordinaatti) kulkee, reaktiopolun maksimi. Siirtymätila ("col") on potentiaalienergiavallin huipulla (ks. oheinen kuva). Potentiaalienergia on siirtymätilassa maksimissa reaktiokoordinaatin suuntaisesti, mutta minimissä reaktiokoordinaattiin suorakulmassa olevien koordinaattien suhteen. Yksittäistasapainoperiaatteen mukaan kemiallisen reaktion ollessa tasapainotilassa, on lähtöaineiden kulkeutuminen reaktion potentiaalienergiavallin yli oltava tarkalleen tasapainossa reaktion tuotteista lähtöaineisiin tapahtuvan kulkeutumisen kanssa. Tämän seurauksena todennäköisin reaktiopolku ja reaktion siirtymätila ovat samat molemmille reaktiosuunnille.s. 372 Siirtymätilateorian mukaan siirtymätilalle valitaan reaktiokoordinaatin suuntaisesti hyvin lyhyt reaktiopolun osuus \delta, jossa sillä on reaktiopolun normaalin suuntaisesti kaikki vibraatio-, rotaatio, ja translaatiovapausasteet. Vain reaktiopolun suuntainen yksi tangentiaalinen vibraatiomoodi (tai yksi translaatiomoodi) on erityisasemassa. Reaktion **aktivoitu kompleksi** omaa muut kaikki siirtymätilan konfiguraation kaltaiset rakenteet siirtymätilan ympäristössä.s. 60 Siirtymätilateorian mukaan reaktiosuunnan reaktionopeus riippuu Gibbsin vapaaenergian suurenemisesta reaktion edetessä lähtöaineista siirtymätilaan. Koska siirtymätila on sama molemmille reaktiosuunnille, niin erityisesti alkeisreaktiolle reaktion Gibbsin vapaaenergian muutos on yhtäsuuri kuin reaktiosuuntien Gibbsin vapaaenergioiden erotus. ==Konventionaalinen siirtymätilateoria== Oheisessa kuvassa E_0 on potentiaalienergiaero siirtymätilan alimman energiatilan ja lähtöaineiden muodostaman yhteisen alimman energiatilan välillä. Näiden tilojen tarkastelussa on nollapiste- energiat on otettu huomioon. Reaktion etenemisen kannalta katsoen E_0, joka vastaa törmäysteorian kynnysenergiaa, on pienin energia, joka tarvitaan reaktion tuotteiden syntymiseksi. Reaktion kinetiikan määrittämisen kannalta katsoen oletetaan, että siirtymätilan ja aktivoidun kompleksin konsentraatio on [C]^{\ddagger}. Aktivoitu kompleksi, joka liikkuu reaktiokoordinaatin suuntaisesti (vasemmalta oikealle) kohti reaktiotuotteita, omaa positiivisen nopeuden, \dot{x}, Tällöin aktivoidun kompleksin keskimääräinen nopeus on kineettisen kaasuteorian mukaan ::(1)\qquad \bar{\dot{x}} = \sqrt{k_B \, T \over 2 \pi \, m^{\ddagger}} Tässä m^{\ddagger} on aktivoidun kompleksin reaktiokoordinaatin suuntaiseen liikkeeseen (1D-laatikko) liittyvä massa. Aktivoidun kompleksin käyttämä keskimääräinen aika siirtymätilalle varatun reaktiokoordinaatin \delta-pituuden kulkemiseksi on t = \frac{\delta}{\bar{\dot{x}}}. Reaktion nopeudelle voidaan nyt kirjoittaa:s. 107s. 164 ::(2)\qquadreaktionopeus \, = \frac{d[C]^{\ddagger}}{dt} = \Bigg( {[C]^{\ddagger} \over \delta} \Bigg) \sqrt{k_B \,T \over 2 \pi \, m^{\ddagger}} Oletetaan reaktion olevan bimolekulaarinen: A \mathit{}+B -> C \mathit{^{\ddagger}}-> tuotteet. Aktivoidun kompleksin ja lähtöaineiden konsentraatioiden suhde Boltzmannin jakauman mukaan on: ::(3)\qquad \frac{[C]^{\ddagger}}{[A][B]} = {{Q_{C^{\ddagger}} \over L \, V} \over \Bigg({Q_A \over L \, V}\Bigg)\Bigg({Q_B \over L \, V}\Bigg)} e^{-\frac {E_0}{k_B T}} Tässä L on Avogadron luku ja V on tilavuus. Tästä ratkaisemalla [C]^{\ddagger} ja sijoittamalla edelliseen yhtälöön saadaan: ::(4)\qquadreaktionopeus \, = (L\,V) \delta^{-1} \sqrt{k_B \,T \over 2 \pi \, m^{\ddagger}} \Bigg({Q_{C^{\ddagger}} \over{Q_A \, Q_B}}\Bigg) e^{-\frac {E_0}{k_B T}}[A][B] Tästä on todettavissa, että bimolekulaarisen reaktion kokonaiskertaluku on 2 ja että alkeisreaktion kokonaiskertaluku on yhtä kuin reaktion molekulaarisuus. Konventionaalisen siirtymätilateorian mukaisesti aktivoidun kompleksin molekulaarinen kokonaisjakaumafunktio eroaa pysyvän molekyylin vastaavasta vain yhden vapausasteen osalta. Tämä liike potentiaalienergiavallin yli \delta-matkan verran muistuttaa translaatiovapauastetta. Matka \delta voidaan valita niin lyhyeksi, että potentiaalienergian voi katsoa olevan vakio. Tällöin voidaan olettaa tämän translaatiovapausasteen käyttäytyvän normaalin translaatiojakaumafunktion mukaisesti: ::(5)\qquad Q_{t,[C]^{\ddagger}} = Q_{\ddagger} \sqrt{2 \pi \, m^{\ddagger} \, k_B \, T} \Bigg(\frac{\delta}{h}\Bigg) Tässä Q_{\ddagger} sisältää kaikki muut siirtymätilan vapausasteet paitsi yhden reaktiokoordinaatin suuntaisen translaatiovapausasteen. Kun tämä jakaumafunktio sijoitetaan edelliseen reaktion nopeuslakiin, saadaan reaktion nopeusvakioksi: ::(6)\qquad k = (L \, V) \Bigg(\frac{k_B \, T}{ h} \Bigg) \frac{Q_{\ddagger}}{Q_A Q_B}e^{-\frac {E_0}{k_B \, T}} Tämä yhtälö voidaan kirjoittaa koskemaan yleisesti yhtälöä aA \mathit{}+bB \mathit{}+...-> C \mathit{^{\ddagger}} -> tuotteet: ::(7)\qquad k = (L \, V)^{a+b+...-1} \Bigg(\frac{R \, T}{L \, h} \Bigg) \frac{Q_{\ddagger}}{Q^{a}_A Q^{b}_B...}e^{-\frac {L \, E_0}{R \, T}} = \Bigg(\frac{R \, T}{L \, h}\Bigg) K^{\ddagger} Tässä R = k_B\,L on yleinen kaasuvakio ja h on Planckin vakio. Yhtälön lopussa olevaa sulkulauseketta kutsutaan taajuustekijäksi ja sen arvo ilmaisee millä taajuudella siirtymätilan rakenteen omaava aktivoitu kompleksi muuntuu reaktion tuotteiksi. Taajuustekijä riippuu vain lämpötilasta ja sen yksikkö on s-1. Yhtälössä K^{\ddagger} on tasapainovakion kaltainen vakio ja L \, E_0 on reaktion molaarinen kynnysenergia, joka todellisuudessa riippuu lähtöainemolekyylien reaktiivisen poikkipinta-alan muuttumisesta törmäysenergian funktiona. ===Konventionaalisen siirtymätilateorian oletukset=== * Klassisella mekaniikalla ei voida täysin selittää liikettä pitkin reaktiokoordinaattia. Kevyellä lähtöaineelle kuten vetyatomi pitäisi ottaa huomioon kvanttimekaaninen tunneloituminen potentiaalivallin läpi. * Aktivoidut kompleksit muodostuvat lähtöainemolekyylien keskinäisistä törmäyksistä. Niiden oletetaan olevan samanlaisia. Kun lähtöainemolekyylien oletetaan noudattavan Boltzmannin jakaumalakia, pitäisi muodostuneen aktivoidun kompleksin myös noudattaa sitä. * Kun aktivoitu kompleksi ylittää potentiaalienergiavallin, ei sen teorian mukaan ole mahdollista dissosioitua takaisin lähtöaineiksi. On kuitenkin mahdollista, että jotkut törmäyksellisesti aktivoidun kompleksin reaktiopolut risteävät potentiaalienergiavallin yli useasti. Tämän korjaamiseksi on käytetty läpäisykerrointa \kappa, jolloin nopeusvakio on: ::(8)\qquad k = \kappa \Bigg(\frac{R\,T}{L\,h}\Bigg) K^{\ddagger} Jos törmäävien lähtöainemolekyylien muodostaman aktivoidun kompleksin energia on vain hiukan enemmän kuin E_0, niin \kappa ≈ 1\. Mikäli kyseessä on lähtöainemolekyylien suurenergiset törmäykset, voi \kappa << 1\. ==Variaatiosiirtymätilateoria== Reaktiodynamiikan näkökulmasta konventionaalisella siirtymätilateorialla laskemalla saadaan liian suuria nopeusvakioita, koska teoria ei mahdollista reaktiopolun risteävän useasti siirtymätilan yli. Tämän epätarkkuuden poistamiseksi variaatiosiirtymätilateoriassa voi siirtymätilan paikkaa liikuttaa pitkin reaktiopolkua.s. 1035 Toisin sanoen variaatiosiirtymätilateorian avulla voidaan laskea reaktion nopeusvakion arvo eri kohdissa reaktiopolkua ja pienin arvo pitäisi olla lähinnä todellista arvoa kun oletetaan ettei kvanttimekaanista tunneloitumista tapahdu. Konventionaalisen teorian nopeusvakion (so. kanoninen nopeusvakio) suhde mikrokanoniseen nopeusvakioon k(E) on seuraava:s. 125 ::(9)\qquad k(T)\,=\, \frac{\int \limits_{0}^{\infty} N(E)\,e^{-\frac{E}{k_B T}}\,k(E)dE}{Q_{\text{LA}}(T)} Tässä N(E) ja Q_{\text{LA}}(T) ovat lähtöainemolekyylien energiatilojen tiheys ja molekulaariset jakaumafunktiot lämpötilassa T. Variaatiosiirtymätilateoria sisältää mikrokanonisen ja kanonisen vaihtoehdon. ==Termodynaaminen siirtymätilateoria== Termodynaamisten funktioiden käyttö molekulaaristen jakaumafunktioiden sijasta on käytännöllistä laskettaessa nopeusvakioita liuosfaasiselle reaktiolle. Tässä lähestymistavassa tavanomaista reaktion tasapainovakiota muistuttava K^{\ddagger} kuvaa siirtymätilan ja lähtöaineiden välistä tasapainotilaa, joka eroaa tavanomaisesta tasapainotilasta. Sille on kuitenkin määritelty termodynaamisia funktioita aivan kuten tasapainoreaktion tasapainovakio on kytketty standardiseen Gibbsin vapaaenergiaan. Määritelmän mukaan standardinen aktivoitumisen Gibbsin vapaaenergia on:s. 492 ::(10)\qquad \Delta G^{\ddagger \ominus} = -R\,T\,\ln [K^{\ddagger}(c^{\ominus})^{m-1}] Tässä m on reaktion molekulaarisuus ja c^{\ominus} on standarditilan konsentraatio. Aktivoitumisen standardinen entalpia on määritelty olevan: ::(11)\qquad \Delta H^{\ddagger \ominus} = R\,T^2 \frac{d\,\ln[K^{\ddagger}_{p} (p^{\ominus})^{m-1}]}{dT} = R\,T^2 \frac{d\,\ln[K^{\ddagger}(c^{\ominus})^{m-1}]}{dT} - (m-1)R\,T Tässä K^{\ddagger}_{p} on aktivoitumisen painetasapainevakio. Aktivoitumisen standardinen entalpia voidaan ilmaista myös konsentraatiotermeillä, koska ihannekaasuille pätee ::(12)\qquad K^{\ddagger}_{p}(p^{\ominus})^{m-1} = K^{\ddagger}(c^{\ominus})^{m-1} \Bigg(\frac{c^{\ominus}R\,T}{p^{\ominus}}\Bigg)^{1-m} Aktivoitumisen standardiseksi entropiaksi on määritelty: ::(13)\qquad \Delta S^{\ddagger \ominus} = \frac{\Delta H^{\ddagger \ominus}-\Delta G^{\ddagger \ominus}}{T} Konventionaalisen siirtymätilateorialla ilmaistu nopeusvakio voidaan nyt esittää edellä mainituilla termodynaamisilla funktioilla: ::(14)\qquad k=(c^{\ominus})^{1-m} \Bigg( \frac{R\,T}{L\,h}\Bigg) e^{-\frac{\Delta G^{\ddagger \ominus}}{R\,T}}=(c^{\ominus})^{1-m} \Bigg( \frac{R\,T}{L\,h}\Bigg) e^{\frac{\Delta S^{\ddagger \ominus}}{R}} e^{-\frac{\Delta H^{\ddagger \ominus}}{R\,T}} Tässä k:n yksikkö on (c^{\ominus})^{1-m} \text{s}^{-1}. Tämän ns. Eyringin yhtälön kehitti Henry Eyring vuonna 1935. Eyringin yhtälöstä voidaan todeta, että esim. kaasufaasisen bimolekulaarisen reaktion Arrheniuksen taajuustekijä on ::(15)\qquad A = \Bigg( \frac{R\,T}{L\,h\,c^{\ominus}} \Bigg)e^2 e^{\frac{\Delta S^{\ddagger \ominus}}{R}} Tämä pätee, koska \Delta H^{\ddagger \ominus} = E_a - m\,R\,T. Aktivoitumisen standardinen entropia korreloi törmäysteorian ns. steerisen tekijän vaikutusta nopeusvakioon. Aktivoitumisen standardinen entropia on negatiivisempi moniatomisten molekyylien reaktiossa kuin atomien reaktiossa. Jos aktivoitumisen standardinen entropian arvo on noin -50 \text{ JK mol}^{-1}, vastaa Arrheniuksen taajuustekijä törmäysteorian avulla laskettua arvoa. ===Termodynaaminen siirtymätilateoria liuosreaktiossa=== Konventionaalisen siirtymätilateorian soveltaminen liuosfaasireaktion kinetiikkaan ei ole mielekästä, koska molekulaaristen jakaumafunktioiden muodostaminen tällöin on työlästä. Tähän monimutkaisuuteen vaikuttaa reaktion liuotin, jolla voi olla dramaattinen vaikutus kemiallisen reaktion potentiaalienergiapinnalle. Esimerkiksi SN2-reaktiossa poolinen liuotin stabiloi lähtöaineita ja tuotteita siirtymätilaa voimakkaammin. Tästä johtuen nopeusvakion arvo liusfaasireaktiolla verrattuna vastaavaan kaasufaasireaktiolla voi olla hyvin pieni.s. 140 Sovellettaessa siirtymätilateoriaa esimerkiksi bimolekulaariseen liusreaktioon, reaktioyhtälö kirjoitetaan ::(16)\qquad \text{A + B}\mathit{\ce{<=>[K_s^{\ddagger}]}}\text{AB}^{\ddagger}\,\,\ce{->[ u^{\ddagger}]} H Tässä \text{AB}^{\ddagger} on reaktion siirtymätila. Tähän reaktioon voidaan soveltaa oletusta tasapainotilasta lähtöaineiden ja siirtymätilan välillä, jos siirtymätilan unimolekulaarinen hajoaminen lähtöaineiksi on paljon nopeampaa kuin sen hajoaminen tuotteiksi. Tämän lisäksi kaasufaasireaktiossa käytetyt konsentraatiot pitää korvata liuosfaasireaktioissa aktiivisuuksilla, joten aktivoitumisen tasapainovakiolle liuosfaasissa,K_s^{\ddagger}, voidaan kirjoittaa ::(17)\qquad K_s^{\ddagger}\,=\,\frac{k_2}{k_{uni}}\,=\,\frac{c^{\ominus}(AB)^{\ddagger}}{(A)(B)}\,=\,\frac{c^{\ominus}[AB]^{\ddagger}\gamma^{\ddagger}}{[A]\gamma_A\,[B]\gamma_B} Oletuksen mukaan k_{uni}\gg u^{\ddagger} ja k_2 on bimolekulaarinen nopeusvakio. Yhtälössä (17) \gamma aktiivisuuskerroin ottaa huomioon poikkeaman kaasufaasisesta "ideaalisesta" olosuhteesta. Reaktionopeus on ::(18)\qquad R\,=\,+\frac{d(H)}{dt}\,=\, u^{\ddagger}(AB^{\ddagger})\,=\, u^{\ddagger}K_s^{\ddagger}(A)(B)\,=\,k(A)(B) Tässä u^{\ddagger} = \frac{R\,T}{L\,h} siirtymätilateorian mukaisesti. Tästä saadaan nopeusvakiolle ::(19)\qquad k\,=\,\Bigg(\frac{R\,T}{L\,h\,c^{\ominus}}\Bigg) K_s^{\ddagger}\frac{\gamma_A \gamma_B}{\gamma^{\ddagger}} Voidaan määritellään, että k^{\ominus} = \frac{R\,T\,K_s^{\ddagger}}{L\,h\,c^{\ominus}}, jossa k^{\ominus} on k:n standardiarvo äärettämässä laimennuksessa ( \gamma=1 ). Usein standardikonsentraatio, c^{\ominus}, on valittu olemaan 1 \text{ mol dm}^{-3}. Täten yhtälö (19) voidaan kirjoittaa myös muotoon ::(20)\qquad k\,=\,k^{\ominus} \frac{\gamma_A \gamma_B}{\gamma^{\ddagger}} Tätä sanotaan Brønsted-Bjerrum-yhtälöksi. Yhtälön (20) tärkeä sovellus on sen käyttö ionien väliseen reaktioon: ::(21)\qquad \text{A}^{z_A} + \text{B}^{z_B} \longrightarrow \left( \text{AB}^{\ddagger} \right)^{z_A + z_B} \longrightarrow \text{tuotteet} Debye-Hückel-teorians. 185 mukaan aktiivisuuskerroin riippuu ionivahvuudestas. 555 ::(22)\qquad \lg\,\gamma_i\,=\,-Bz_i^2\,\sqrt{I} Yhtälössä B voidaan laskea Debye-Hückel- yhtälöistä ja sen arvo 298 K:n lämpöisessä vesiliuoksessa on 0,509 \text{ dm}^{-3/2}\text{mol}^{-1/2}. Ionivahvuus I = \frac{1}{2}m_j z_j, jossa m_j on ionisten ainesosien molaalisuus liuoksessa. Ionivahvuuden yksikkö on \text{mol dm}^{-3}. Sijoittamalla yhtälö (21) yhtälöön (19) saadaan ::(23)\qquad \lg\,k\,=\,\lg\,k^{\ominus} + \lg\,\gamma_A + \lg\,\gamma_B - \lg\,\gamma^{\ddagger}\,=\,\lg\,k^{\ominus} + 2B z_A z_B\sqrt{I} Kun olosuhteiksi otetaan 298 K:n vesilius, niin nopeusvakiosuhde on ::(24)\qquad \lg\frac{k}{k^{\ominus}}\,=\,1,018 z_A z_B\sqrt{I} Tämä yhtälö on hyödyllinen kun verrataan kaasufaasi- ja liuosfaasireaktioita keskenään ja jos vertailutila valitaan niin, että k^{\ominus} vastaa kaasufaasista nopeusvakiota. Piirrettäessä kuvaaja \lg\,(k/k^{\ominus}) versus \sqrt{I}, on kuvaajan kulmakerroin 1,018 z_A z_B. ==Paineen vaikutus nopeusvakioon== Reaktiopaineen vaikutus nopeusvakioon on normaalisti vähäinen verrattuna nopeusvakion lämpötilariippuvuuteen. Poikkeuksena tästä ovat unimolekulaariset reaktiot. Vasta verrattain suuret, yli sadan kilobaarin, reaktiopaineet vaikuttavat nopeusvakion arvoon mitattavasti. Kaasuräjähdys on käytännön esimerkki paineen vaikutuksesta. Tarkasteltaessa paineen vaikutusta nopeuvakioon, lähtökohdaksi voidaan ottaa yhtälö (7). Tämä voidaan uudelleen kirjoittaa ottamalla yhtälöstä luonnollinen logaritmi ja käyttäen aktivoitumisen vapaaenergiaa konsentraatiotasapainovakiolle (ks. yhtälö (10)): ::(25)\qquad \ln\,k\,=\,\ln \Bigg(\frac{\,R\,T}{L\,h}\Bigg)-\frac{\Delta G^{\ddagger \ominus}}{R\,T} Tässä yhtälössä ensimmäinen termi on vakio. Derivoimalla yhtälö paineen suhteen,{\partial P}\Bigg)_T\,=\,-\frac{1}{RT}\Bigg(\frac{\partial \Delta G^{\ddagger \ominus}}{\partial P}\Bigg)_T\,=\,-\frac{\Delta V^{\ddagger}}{RT} Laimeissa liuoksissa aktiivisuuskertoimilla on vain vähäinen paineriippuvuus, joten yhtälöstä (19) saadaan nopeusvakion paineriippuvuudelle yhtälö (28).}} saadaan nopeusvakion paineriippuvuudelles. 875 ::(28)\qquad \Bigg(\frac{\partial \ln\,k}{\partial P}\Bigg)_T\,\simeq\,\Bigg(\frac{\partial \ln\,K_s^{\ddagger}}{\partial P}\Bigg)_T\,=\,-\frac{\Delta V^{\ddagger \ominus}}{R\,T} Tässä \Delta V^{\ddagger \ominus} on standardinen aktivoitumisen molaarinen tilavuus. Integroitaessa yhtälö (28) ::(29)\qquad \int_{k_0}^k\,d\ln k\,=\,-\frac{\Delta V^{\ddagger \ominus}}{RT}\int_0^P dP saadaan ::(30)\qquad RT\,\ln\frac{k}{k_0}\,=\,-P\,\Delta V^{\ddagger \ominus} Tässä R = 0,082 \text{ dm}^3\text{ atm mol}^{-1}\text{ T}^{-1} ja k_0 on nopeusvakio vertailupaineessa (usein 1 atm). Tyypillisesti nopeusvakion pienen paineriippuvuuden takia mittauksissa reaktiopaineen P on oltava vähintään \text{kbar}, jotta nopeusvakio k:lla olisi mitattavissa oleva muutos. Esimerkiksi jos \Delta V^{\ddagger \ominus} = \pm 15 \text{ ml mol}^{-1}, niin reaktiopaineen ollessa 1000 atm lämpötilassa 298 K, muuttuu nopeusvakio 80%. == Huomautukset == ==Katso myös== *Reaktiokinetiikka *Arrheniuksen yhtälö *Alkeisreaktio *Boltzmannin jakauma *Molekulaarinen jakaumafunktio *Törmäysteoria *Unimolekulaarinen reaktio *RRKM == Lähteet == 1\. John W. Moore ja Ralph G. Pearson, _Kinetics and Mechanism_ , 3. painos, (1981), John Wiley & Sons, ISBN 0-471-03558-0 2\. Keith J. Laidler, _Chemical Kinetics_ , 3. painos, (1987), HarperCollinsPublisher, ISBN 0-06-043862-2 3\. H. Eyring, _J. Chem. Phys._ , vol 3, (1935) 4\. J.C. Keck, _J. Chem. Phys._ , vol 32, (1960) 5\. W.F.K. Wynne-Jones ja H. Eyring, _J. Chem. Phys._ , vol 3, (1935) 6\. Frank Wilkinson, _Chemical Kinetics and Reaction Mechanisms_ , (1980), van Nostrand Reinhold Company, ISBN 0-442-30248-7 7\. P. J. W. Debye ja E. Hückel, _Phys. Z._ , vol 24, (1923) 8\. J. N. Brønsted ja V. K. La Mer, J. Am. Chem. Soc., vol 46, (1924) 9\. Irving M. Klotz ja Robert M. Rosenberg, _Chemical Thermodynamics_ , 5. painos, (1994), John Wiley & Sons, ISBN 0-471-53439-0 10\. M. G. Evans ja M. Polanyi, Trans. Faraday Soc., vol 31, (1935) Luokka:Kemiallinen kinetiikka