# Autocatalysis

A single chemical reaction is said to be autocatalytic if one of the reaction products is also a catalyst for the same or a coupled reaction.[1] Such a reaction is called an autocatalytic reaction.

A set of chemical reactions can be said to be "collectively autocatalytic" if a number of those reactions produce, as reaction products, catalysts for enough of the other reactions that the entire set of chemical reactions is self-sustaining given an input of energy and food molecules (see autocatalytic set).

## Chemical reactions

A chemical reaction of two reactants and two products can be written as

${\displaystyle \alpha A+\beta B\rightleftharpoons \sigma S+\tau T}$

where the Greek letters are stoichiometric coefficients and the capital Latin letters represent chemical species. The chemical reaction proceeds in both the forward and reverse direction. This equation is easily generalized to any number of reactants, products, and reactions.

### Chemical equilibrium

In chemical equilibrium the forward and reverse reaction rates are such that each chemical species is being created at the same rate it is being destroyed. In other words, the rate of the forward reaction is equal to the rate of the reverse reaction.

${\displaystyle k_{+}[A]^{\alpha }[B]^{\beta }=k_{-}[S]^{\sigma }[T]^{\tau }\,}$

Here, the brackets indicate the concentration of the chemical species, in moles per liter, and k+ and k are rate constants.

### Far from equilibrium

Far from equilibrium, the forward and reverse reaction rates no longer balance and the concentration of reactants and products is no longer constant. For every forward reaction ${\displaystyle \alpha }$ molecules of A are destroyed. For every reverse reaction ${\displaystyle \alpha }$ molecules of A are created. In the case of an elementary reaction step the reaction order in each direction equals the molecularity, so that the rate of change in the number of moles of A is then

${\displaystyle {d \over dt}[A]=-\alpha k_{+}[A]^{\alpha }[B]^{\beta }+\alpha k_{-}[S]^{\sigma }[T]^{\tau }\,}$
${\displaystyle {d \over dt}[B]=-\beta k_{+}[A]^{\alpha }[B]^{\beta }+\beta k_{-}[S]^{\sigma }[T]^{\tau }\,}$
${\displaystyle {d \over dt}[S]=\sigma k_{+}[A]^{\alpha }[B]^{\beta }-\sigma k_{-}[S]^{\sigma }[T]^{\tau }\,}$
${\displaystyle {d \over dt}[T]=\tau k_{+}[A]^{\alpha }[B]^{\beta }-\tau k_{-}[S]^{\sigma }[T]^{\tau }\,}$

This system of equations has a single stable fixed point when the forward rates and the reverse rates are equal (when ${\displaystyle {d \over dt}=0}$ for every species). This means that the system evolves to the equilibrium state, and this is the only state to which it evolves.[2]

### Autocatalytic reactions

Sigmoid variation of product concentration in autocatalytic reactions

Autocatalytic reactions are those in which at least one of the products is a reactant. Perhaps the simplest autocatalytic reaction can be written[1]

${\displaystyle A+B\rightleftharpoons 2B}$

with the rate equations (for an elementary reaction)

${\displaystyle {d \over dt}[A]=-k_{+}[A][B]+k_{-}[B]^{2}\,}$
${\displaystyle {d \over dt}[B]=+k_{+}[A][B]-k_{-}[B]^{2}\,}$.

This reaction is one in which a molecule of species A interacts with a molecule of species B. The A molecule is converted into a B molecule. The final product consists of the original B molecule plus the B molecule created in the reaction.

The key feature of these rate equations is that they are nonlinear; the second term on the right varies as the square of the concentration of B. This feature can lead to multiple fixed points of the system, much like a quadratic equation can have multiple roots. Multiple fixed points allow for multiple states of the system. A system existing in multiple macroscopic states is more orderly (has lower entropy) than a system in a single state.

The concentrations of A and B vary in time according to[1][3]

${\displaystyle [A]={\frac {[A]_{0}+[B]_{0}}{1+{\frac {[B]_{0}}{[A]_{0}}}e^{([A]_{0}+[B]_{0})kt}}}}$

and

${\displaystyle [B]={\frac {[A]_{0}+[B]_{0}}{1+{\frac {[A]_{0}}{[B]_{0}}}e^{-([A]_{0}+[B]_{0})kt}}}}$.

The graph for these equations is a sigmoid curve (specifically a logistic function), which is typical for autocatalytic reactions: these chemical reactions proceed slowly at the start (the induction period) because there is little catalyst present, the rate of reaction increases progressively as the reaction proceeds as the amount of catalyst increases and then it again slows down as the reactant concentration decreases. If the concentration of a reactant or product in an experiment follows a sigmoid curve, the reaction may be autocatalytic.

These kinetic equations apply for example to the acid-catalyzed hydrolysis of some esters to carboxylic acids and alcohols.[3] There must be at least some acid present initially to start the catalyzed mechanism; if not the reaction must start by an alternate uncatalyzed path which is usually slower. The above equations for the catalyzed mechanism would imply that the concentration of acid product remains zero forever.[3]

## Creation of order

### Background

The second law of thermodynamics states that the disorder (entropy) of a physical or chemical system and its surroundings (a closed system) must increase with time. Systems left to themselves become increasingly random, and orderly energy of a system like uniform motion degrades eventually to the random motion of particles in a heat bath.

There are, however, many instances in which physical systems spontaneously become emergent or orderly. For example, despite the destruction they cause, hurricanes have a very orderly vortex motion when compared to the random motion of the air molecules in a closed room. Even more spectacular is the order created by chemical systems; the most dramatic being the order associated with life.

This is consistent with the Second Law, which requires that the total disorder of a system and its surroundings must increase with time. Order can be created in a system by an even greater decrease in order of the system's surroundings.[4] In the hurricane example, hurricanes are formed from unequal heating within the atmosphere. The Earth's atmosphere is then far from thermal equilibrium. The order of the Earth's atmosphere increases, but at the expense of the order of the sun. The sun is becoming more disorderly as it ages and throws off light and material to the rest of the universe. The total disorder of the sun and the earth increases despite the fact that orderly hurricanes are generated on earth.

A similar example exists for living chemical systems. The sun provides energy to green plants. The green plants are food for other living chemical systems. The energy absorbed by plants and converted into chemical energy generates a system on earth that is orderly and far from chemical equilibrium. Here, the difference from chemical equilibrium is determined by an excess of reactants over the equilibrium amount. Once again, order on earth is generated at the expense of entropy increase of the sun. The total entropy of the earth and the rest of the universe increases, consistent with the Second Law.

Some autocatalytic reactions also generate order in a system at the expense of its surroundings. For example, (clock reactions) have intermediates whose concentrations oscillate in time, corresponding to temporal order. Other reactions generate spatial separation of chemical species corresponding to spatial order. More complex reactions are involved in metabolic pathways and metabolic networks in biological systems.

The transition to order as the distance from equilibrium increases is not usually continuous. Order typically appears abruptly. The threshold between the disorder of chemical equilibrium and order is known as a phase transition. The conditions for a phase transition can be determined with the mathematical machinery of non-equilibrium thermodynamics.

### Temporal order

A chemical reaction cannot oscillate about a position of final equilibrium because the second law of thermodynamics requires that a thermodynamic system approach equilibrium and not recede from it. For a closed system at constant temperature and pressure, the Gibbs free energy must decrease continuously and not oscillate. However it is possible that the concentrations of some reaction intermediates oscillate, and also that the rate of formation of products oscillates.[5]

#### Idealized example: Lotka–Volterra equation

The Lotka–Volterra equation is isomorphic with the predator–prey model and the two-reaction autocatalytic model. In this example baboons and cheetahs are equivalent to two different chemical species X and Y in autocatalytic reactions.

Consider a coupled set of two autocatalytic reactions in which the concentration of one of the reactants A is much larger than its equilibrium value. In this case, the forward reaction rate is so much larger than the reverse rates that we can neglect the reverse rates.

${\displaystyle A+X\rightarrow 2X}$
${\displaystyle X+Y\rightarrow 2Y}$
${\displaystyle Y\rightarrow E}$

with the rate equations

${\displaystyle {d \over dt}[X]=k_{1}[A][X]-k_{2}[X][Y]\,}$
${\displaystyle {d \over dt}[Y]=k_{2}[X][Y]-k_{3}[Y]\,}$.

Here, we have neglected the depletion of the reactant A, since its concentration is so large. The rate constants for the three reactions are ${\displaystyle k_{1}}$, ${\displaystyle k_{2}}$, and ${\displaystyle k_{3}}$, respectively.

This system of rate equations is known as the Lotka–Volterra equation and is most closely associated with population dynamics in predator–prey relationships. This system of equations can yield oscillating concentrations of the reaction intermediates X and Y. The amplitude of the oscillations depends on the concentration of A (which decreases without oscillation). Such oscillations are a form of emergent temporal order that is not present in equilibrium.

#### Another idealized example: Brusselator

Another example of a system that demonstrates temporal order is the Brusselator (see Prigogine reference). It is characterized by the reactions

${\displaystyle A\rightarrow X}$
${\displaystyle 2X+Y\rightarrow 3X}$
${\displaystyle B+X\rightarrow Y+D}$
${\displaystyle X\rightarrow E}$

with the rate equations

${\displaystyle {d \over dt}[X]=[A]+[X]^{2}[Y]-[B][X]-[X]\,}$
${\displaystyle {d \over dt}[Y]=[B][X]-[X]^{2}[Y]\,}$

where, for convenience, the rate constants have been set to 1.

The Brusselator in the unstable regime. A=1. B=2.5. X(0)=1. Y(0)=0. The system approaches a limit cycle. For B<1+A the system is stable and approaches a fixed point.

The Brusselator has a fixed point at

${\displaystyle [X]=A\,}$
${\displaystyle [Y]={B \over A}\,}$.

The fixed point becomes unstable when

${\displaystyle B>1+A^{2}\,}$

leading to an oscillation of the system. Unlike the Lotka-Volterra equation, the oscillations of the Brusselator do not depend on the amount of reactant present initially. Instead, after sufficient time, the oscillations approach a limit cycle.[6]

### Spatial order

An idealized example of spatial spontaneous symmetry breaking is the case in which we have two boxes of material separated by a permeable membrane so that material can diffuse between the two boxes. It is assumed that identical Brusselators are in each box with nearly identical initial conditions. (see Prigogine reference)

${\displaystyle {d \over dt}[X_{1}]=[A]+[X_{1}]^{2}[Y_{1}]-[B][X_{1}]-[X_{1}]+D_{x}\left(X_{2}-X_{1}\right)\,}$
${\displaystyle {d \over dt}[Y_{1}]=[B][X_{1}]-[X_{1}]^{2}[Y_{1}]+D_{y}\left(Y_{2}-Y_{1}\right)\,}$
${\displaystyle {d \over dt}[X_{2}]=[A]+[X_{2}]^{2}[Y_{2}]-[B][X_{2}]-[X_{2}]+D_{x}\left(X_{1}-X_{2}\right)\,}$
${\displaystyle {d \over dt}[Y_{2}]=[B][X_{2}]-[X_{2}]^{2}[Y_{2}]+D_{y}\left(Y_{1}-Y_{2}\right)\,}$

Here, the numerical subscripts indicate which box the material is in. There are additional terms proportional to the diffusion coefficient D that account for the exchange of material between boxes.

If the system is initiated with the same conditions in each box, then a small fluctuation will lead to separation of materials between the two boxes. One box will have a predominance of X, and the other will have a predominance of Y.

## Real examples

Real examples of clock reactions are the Belousov–Zhabotinsky reaction (BZ reaction), the Briggs–Rauscher reaction, the Bray–Liebhafsky reaction and the iodine clock reaction. These are oscillatory reactions, and the concentration of products and reactants can be approximated in terms of damped oscillations.

The best-known reaction, the BZ reaction, can be created with a mixture of potassium bromate ${\displaystyle {\ce {(KBrO3)}}}$, malonic acid ${\displaystyle {\ce {(CH2(COOH)2)}}}$, and manganese sulfate ${\displaystyle {\ce {(MnSO4)}}}$ prepared in a heated solution with sulfuric acid ${\displaystyle {\ce {(H2SO4)}}}$ as solvent.[7]

## Optics example

Another autocatalytic system is one driven by light coupled to photo-polymerization reactions. In a process termed optical autocatalysis, positive feedback is created between light intensity and photo-polymerization rate, via polymerization-induced increases in the refractive index. Light's preference to occupy regions of higher refractive index results in leakage of light into regions of higher molecular weight, thereby amplifying the photo-chemical reaction. The positive feedback may be expressed as:[8]

${\displaystyle {\text{polymerization rate}}\to {\text{molecular weight}}/{\text{refractive index}}\to {\text{intensity}}}$

Noting that photo-polymerization rate is proportional to intensity[9] and that refractive index is proportional to molecular weight,[10] the positive feedback between intensity and photo-polymerization establishes the auto-catalytic behavior. Optical auto-catalysis has been shown to result on spontaneous pattern formation in photopolymers.[11][12][13] Hosein and co-workers discovered that optical autocatalysis can also occur in photoreactive polymer blends, and that the process can induce binary phase morphologies with the same pattern as the light profile.[8] The light undergoes optical modulation instability, spontaneous dividing into a multitude of optical filaments, and the polymer system thereby forms filaments within the blend structure.[8] The result is a new system that couples optical autocatalytic behavior to spinodal decomposition.

## Biological example

It is known that an important metabolic cycle, glycolysis, displays temporal order.[14] Glycolysis consists of the degradation of one molecule of glucose and the overall production of two molecules of ATP. The process is therefore of great importance to the energetics of living cells. The global glycolysis reaction involves glucose, ADP, NAD, pyruvate, ATP, and NADH.

${\displaystyle {\ce {glucose{}+2ADP{}+2P_{\mathit {i}}{}+2NAD->2(pyruvate){}+2ATP{}+2NADH}}}$.

The details of the process are quite involved, however, a section of the process is autocatalyzed by phosphofructokinase (PFK). This portion of the process is responsible for oscillations in the pathway that lead to the process oscillating between an active and an inactive form. Thus, the autocatalytic reaction can modulate the process.

## Shape tailoring of thin layers

It is possible to use the results from an autocatalytic reaction coupled with reaction–diffusion system theory to tailor the design of a thin layer. The autocatalytic process allows controlling the nonlinear behavior of the oxidation front, which is used to establish the initial geometry needed to generate the arbitrary final geometry.[15] It has been successfully done in the wet oxidation of ${\displaystyle Al_{x}Ga_{1-x}As}$ to obtain arbitrary shaped layers of ${\displaystyle AlO_{x}}$.

## Phase transitions

The initial amounts of reactants determine the distance from a chemical equilibrium of the system. The greater the initial concentrations the further the system is from equilibrium. As the initial concentration increases, an abrupt change in order occurs. This abrupt change is known as phase transition. At the phase transition, fluctuations in macroscopic quantities, such as chemical concentrations, increase as the system oscillates between the more ordered state (lower entropy, such as ice) and the more disordered state (higher entropy, such as liquid water). Also, at the phase transition, macroscopic equations, such as the rate equations, fail. Rate equations can be derived from microscopic considerations. The derivations typically rely on a mean field theory approximation to microscopic dynamical equations. Mean field theory breaks down in the presence of large fluctuations (see Mean field theory article for a discussion). Therefore, since large fluctuations occur in the neighborhood of a phase transition, macroscopic equations, such as rate equations, fail. As the initial concentration increases further, the system settles into an ordered state in which fluctuations are again small. (see Prigogine reference)

## Asymmetric autocatalysis

Asymmetric autocatalysis occurs when the reaction product is chiral and thus acts as a chiral catalyst for its own production. Reactions of this type, such as the Soai reaction, have the property that they can amplify a very small enantiomeric excess into a large one. This has been proposed as an important step in the origin of biological homochirality.[16]

## Role in origin of life

In 1995 Stuart Kauffman proposed that life initially arose as autocatalytic chemical networks.[17]

British ethologist Richard Dawkins wrote about autocatalysis as a potential explanation for abiogenesis in his 2004 book The Ancestor's Tale. He cites experiments performed by Julius Rebek and his colleagues at the Scripps Research Institute in California in which they combined amino adenosine and pentafluorophenyl ester with the autocatalyst amino adenosine triacid ester (AATE). One system from the experiment contained variants of AATE which catalyzed the synthesis of themselves. This experiment demonstrated the possibility that autocatalysts could exhibit competition within a population of entities with heredity, which could be interpreted as a rudimentary form of natural selection, and that certain environmental changes (such as irradiation) could alter the chemical structure of some of these self-replicating molecules (an analog for mutation) in such ways that could either boost or interfere with its ability to react, thus boosting or interfering with its ability to replicate and spread in the population.[18]

Autocatalysis plays a major role in the processes of life. Two researchers who have emphasized its role in the origins of life are Robert Ulanowicz [19] and Stuart Kauffman.[20]

Autocatalysis occurs in the initial transcripts of rRNA. The introns are capable of excising themselves by the process of two nucleophilic transesterification reactions. The RNA able to do this is sometimes referred to as a ribozyme. Additionally, the citric acid cycle is an autocatalytic cycle run in reverse.

Ultimately, biological metabolism itself can be seen as a vast autocatalytic set, in that all of the molecular constituents of a biological cell are produced by reactions involving this same set of molecules.

## References

1. ^ a b c Steinfeld J.I., Francisco J.S. and Hase W.L. Chemical Kinetics and Dynamics (2nd ed., Prentice-Hall 1999) p.151-2 ISBN 0-13-737123-3
2. ^ Ross, John; Garcia-Colin, Leopoldo S. (March 1989). "Thermodynamics of chemical systems far from equilibrium". The Journal of Physical Chemistry. 93 (5): 2091–2092. doi:10.1021/j100342a075.
3. ^ a b c Moore J.W. and Pearson R.G. Kinetics and Mechanism (John Wiley 1981) p.26 ISBN 0-471-03558-0
4. ^ Ilya Prigogine (1980). From Being to Becoming: Time and Complexity in the Physical Sciences. San Francisco: W. H. Freeman. ISBN 978-0-7167-1107-0.
5. ^ Espenson, J.H. Chemical Kinetics and Reaction Mechanisms (2nd ed., McGraw-Hill 2002) p.190 ISBN 0-07-288362-6
6. ^ "Archived copy" (PDF). Archived from the original (PDF) on 2008-12-17. Retrieved 2015-10-15.CS1 maint: archived copy as title (link) Dynamics of the Brusselator
7. ^ Peterson, Gabriel. "The Belousov-Zhabotinsky Reaction". Archived from the original on December 31, 2012.
8. ^ a b c Biria, Saeid; Malley, Phillip P. A.; Kahan, Tara F.; Hosein, Ian D. (2016-11-15). "Optical Autocatalysis Establishes Novel Spatial Dynamics in Phase Separation of Polymer Blends during Photocuring". ACS Macro Letters. 5 (11): 1237–1241. doi:10.1021/acsmacrolett.6b00659.
9. ^ Decker, Christian (1998-02-01). "The use of UV irradiation in polymerization". Polymer International. 45 (2): 133–141. doi:10.1002/(SICI)1097-0126(199802)45:2<133::AID-PI969>3.0.CO;2-F.
10. ^ Askadskii, A. A. (1990). "Influence of crosslinking density on the properties of polymer networks". Polymer Science U.S.S.R. 32 (10): 2061–2069. doi:10.1016/0032-3950(90)90361-9.
11. ^ Burgess, Ian B.; Shimmell, Whitney E.; Saravanamuttu, Kalaichelvi (2007-04-01). "Spontaneous Pattern Formation Due to Modulation Instability of Incoherent White Light in a Photopolymerizable Medium". Journal of the American Chemical Society. 129 (15): 4738–4746. doi:10.1021/ja068967b. ISSN 0002-7863. PMID 17378567.
12. ^ Basker, Dinesh K.; Brook, Michael A.; Saravanamuttu, Kalaichelvi (2015). "Spontaneous Emergence of Nonlinear Light Waves and Self-Inscribed Waveguide Microstructure during the Cationic Polymerization of Epoxides". The Journal of Physical Chemistry C. 119 (35): 20606–20617. doi:10.1021/acs.jpcc.5b07117.
13. ^ Biria, Saeid; Malley, Philip P. A.; Kahan, Tara F.; Hosein, Ian D. (2016-03-03). "Tunable Nonlinear Optical Pattern Formation and Microstructure in Cross-Linking Acrylate Systems during Free-Radical Polymerization". The Journal of Physical Chemistry C. 120 (8): 4517–4528. doi:10.1021/acs.jpcc.5b11377. ISSN 1932-7447.
14. ^ G. Nicolis and Ilya Prigogine (1977). Self-Organization in Nonequilibrium Systems. New York: John Wiley and Sons. ISBN 978-0-471-02401-9.
15. ^ Alfaro-Bittner, K.; Rojas, R.G.; Lafleur, G.; Calvez, S.; Almuneau, G.; Clerc, M.G.; Barbay, S. (22 April 2019). "Modeling the Lateral Wet Oxidation of into Arbitrary Mesa Geometries". Physical Review Applied. 11 (4): 044067. doi:10.1103/PhysRevApplied.11.044067.
16. ^ Soai K, Sato I, Shibata T (2001). "Asymmetric autocatalysis and the origin of chiral homogeneity in organic compounds". The Chemical Record. 1 (4): 321–32. doi:10.1002/tcr.1017. PMID 11893072.
17. ^ Stuart Kauffman (1995). At Home in the Universe: The Search for the Laws of Self-Organization and Complexity. Oxford University Press. ISBN 978-0-19-509599-9.
18. ^ Rebeck, Julius (July 1994). "Synthetic Self-Replicating Molecules". Scientific American: 48–55.
19. ^ Ecology, the Ascendent Perspective", Robert Ulanowicz, Columbia Univ. Press 1997.
20. ^ Investigations, Stuart Kauffman.
21. ^ Kovacs KA, Grof P, Burai L, Riedel M (2004). "Revising the Mechanism of the Permanganate/Oxalate Reaction". J. Phys. Chem. A. 108 (50): 11026–11031. Bibcode:2004JPCA..10811026K. doi:10.1021/jp047061u.
22. ^ Yin, Xi; Shi, Miao; Wu, Jianbo; Pan, Yung-Tin; Gray, Danielle L.; Bertke, Jeffery A.; Yang, Hong (11 September 2017). "Quantitative Analysis of Different Formation Modes of Platinum Nanocrystals Controlled by Ligand Chemistry". Nano Letters. 17 (10): 6146–6150. Bibcode:2017NanoL..17.6146Y. doi:10.1021/acs.nanolett.7b02751. PMID 28873317.