PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of rsosThe Royal Society PublishingRoyal Society Open ScienceAboutBrowse by SubjectAlertsEditorial Board
 
R Soc Open Sci. 2017 November; 4(11): 170708.
Published online 2017 November 8. doi:  10.1098/rsos.170708
PMCID: PMC5717635

A DFT investigation of the blue bottle experiment: E[composite function (small circle)]half-cell analysis of autoxidation catalysed by redox indicators

Abstract

The blue bottle experiment is a collective term for autoxidation reactions catalysed by redox indicators. The reactions are characterized by their repeatable cycle of colour changes when shaken/left to stand and intricate chemical pattern formation. The blue bottle experiment is studied based on calculated solution-phase half-cell reduction potential of related reactions. Our investigation confirms that the reaction in various versions of the blue bottle experiment published to date is mainly the oxidation of an acyloin to a 1,2-dicarbonyl structure. In the light of the calculations, we also propose new non-acyloin reducing agents for the experiment. These results can help guide future experimental studies on the blue bottle experiment.

Keywords: blue bottle experiment, reduction potential, density functional theory, chemical pattern formation, reaction mechanism

1. Introduction

Autoxidation of reducing agents catalysed by redox indicators have been reported in the literature since 1946 [1,2]. The most notable reaction is the ‘blue bottle experiment’ [3], an oxidation of glucose catalysed by methylene blue under an alkaline condition [4]. For bulk reaction in a flask, methylene blue is reduced to colourless (leuco) form by the aldose sugar when left to stand and is oxidized to blue form by atmospheric oxygen when shaken. The cycle can be repeated many times before the reactants run out or the solution turns brown due to side reactions [5]. For thin-layer reaction in a Petri dish, dot and line patterns of the oxidized indicator develop over a period of time [68]. Figure 1 shows the general reaction framework and figure 2 shows patterns and colours observed in different variations of the reaction.

Figure 1.
The general framework of the blue bottle experiment.
Figure 2.
Patterns formed in a green version of the reaction (ascorbic acid is reducing agent) and colours of various dyes in the rapid version of the reaction (benzoin is reducing agent). Adapted with permission from Rajchakit and Limpanuparb [7,8]. Copyright ...

The blue bottle reaction and its analogues are popular chemical demonstrations due to their visual appeal and simplicity, and the majority of the research on this topic is published in the Journal of Chemical Education [118]. Fundamental research and non-education applications of the reactions have also been discussed elsewhere [1926].

In essence, there were attempts to test new reducing agents [7,11,16] and indicators [7,8], to elucidate the mechanism and kinetics of the reactions [5,11] and to model the pattern formation [6,1923,26]. Pattern formation may be comparable to the Belousov–Zhabotinsky [27] and the Briggs–Rauscher [28] reactions. However, the blue bottle model is relatively simple and versatile because it requires only a few reactants and many alternative reagents can produce similar results.

Despite numerous reports, the understanding of the reaction is advanced incrementally by mostly trial-and-error experiments. A number of papers reported only one new reducing agent [16] or indicator [9,10,17] or pattern formation in one specific system [6]. The reports can also be conflicting or incomplete. For example, in 2012, Anderson et al. [5] suggested that gluconate is not a major product and the reaction may produce hydrogen peroxide, but later work [26] as late as 2014 still discusses the gluconate compound as the main product; in 1974, Chen [29] reported the use of indophenol as a dye for the experiment, but, in 2016, Rajchakit & Limpanuparb [7] failed to reproduce it. Experimental reports usually mention only a structure of a dye in its solid form but do not explicitly show oxidized or reduced form(s) of the compound [7,17] and there was no experimental identification of the products in all cases except one [5].

In this first density functional theory (DFT) investigation of the blue bottle experiment, we aim to propose a theoretical framework to resolve discrepancies in the current literature and guide future experimental studies. The manuscript is structured as follows: Methodology describes reactions and computation approach; Results and discussion presents the main results based on reduction potentials, and preliminary experimental evidence, detailed computational/experimental results are given as the electronic supplementary material; Concluding remarks and future work are discussed at the end of the paper.

2. Methodology

2.1. Preliminary consideration

Figure 1 shows that there are three main groups of reactions in the blue bottle experiment: oxygen reduction reactions (ORRs), oxidation/reduction of redox indicators and oxidation of reducing agents. It is natural to characterize these redox reactions in terms of standard half-cell potential, Ehalf-cell in aqueous solution at 298.15 K. Because Ehalf-cell is a ‘per electron’ quantity, it conveniently allows quick comparison and helps with our prediction whether a compound can possibly be used in a blue bottle reaction. By considering the potentials, it is equivalent to the consideration of Gibbs energy. A reaction is spontaneous provided that the cell potential made by combination of reduction potentials of two half-reactions,

Ecell=EcathodeEanode
2.1

is positive. In other words, a necessary but not sufficient condition for the combination of ORR, oxidation/reduction of dyes and oxidation of reducing agents to make up a blue bottle experiment is

EORR>Edye>Ereducingagent
2.2

as shown in figure 1. All the discussions that follow use the same potential comparison process as a thinking framework.

We include representative compounds reported in the blue bottle literature and possible reagents to explore alternative redox indicators/reducing agents and to gain mechanistic insight of the reaction. Table 1 lists the oxidized and reduced structures of all compounds in this study. If oxidized/reduced form(s) of the compounds are not explicitly mentioned in the literature, we do our best to propose them.

Table 1.
Ereduction of all compounds in this study calculated at B3LYP/6-311++G** and SMD solvation model. For complete half-reactions in acid and base, refer to the electronic supplementary material (calculations.xlsx). The calculated difference in energy ...

2.2. Grouping of oxygen reduction reactions, dyes and reducing agents

Reactions reported in the current literature and our proposal for dyes and reagents are studied as follows (structurally similar compounds are grouped together):

  1. The ORRs are trivial but instead of using the literature values [30,31], the calculations were completed to obtain reference values for comparison purpose.
  2. The dye oxidation/reduction reactions are grouped into:
    • (i) heterocyclic azine: oxazine, thiazine and pyrazine,
    • (ii) triarylmethane, and
    • (iii) other common redox dyes.
  3. The oxidation of reducing agents are grouped into
    • (i) alcohol [right harpoon over left harpoon] aldehyde/ketone
      •  simple substrates
      •  acyloin including phenols, pyranoses and lactones
      •  α-hydroxy carboxylic acids,
    • (ii) aldehyde [right harpoon over left harpoon] carboxylic acid,
    • (iii) sulfur-containing amino acids, and
    • (iv) other common reducing agents.

2.3. Computational details

Different procedures [3250] exist for computation of accurate Ehalf-cell especially in one-electron case and/or families of structurally similar compounds [3236,42,4446,48]. Since all species in our study are simply closed-shell singlet organic molecules, however, a gas-phase DFT optimization followed by a free energy of solvation calculation has been proven successful in many cases [40,43,4748]. This approach was also included and tested in recent reviews [4951].

Gas-phase geometries were obtained at B3LYP/6-311++G** level and were confirmed to be a minimum point on the potential energy surface by frequency calculation. Solvation was treated by SMD model [52] on the gas-phase structure. Some compounds in our study have a number of rotamers and diastereomers. We try to use the lowest energy structure as a representative. However, the difference due to these stereoisomers is expected to be small (1 kcal mole−1 of electron is approximately 0.04 V). All output files are provided in the electronic supplementary material. (Additional calculation at B3LYP/6-31G* (gas phase) and MP2/cc-pVTZ (solution phase) were also completed on selected compounds for the preparation of initial structures for B3LYP/6-311++G** and for benchmarking, respectively.)

All calculations were performed using the Q-Chem 4.4 developer version [53]. The half-cell reduction potential was directly calculated from these equations:

For a half-cell reaction:

ΔGsoln=nFE.
2.3

For a chemical structure [49]:

Gsoln=HgasTSgas+ΔGsolv,
2.4

Hgas=ε0(gas)+Hcorr(gas),
2.5

where ΔG is the standard Gibbs energy of the half reaction, n is the number of electrons in the half reaction, F is the Faraday constant 96 485 C mol−1, G is the standard Gibbs energy, H is the standard enthalpy, T is 298.15 K, S is the standard entropy, ΔGsolv is the free energy of solvation from SMD and standard state correction, ε0 is the electronic energy (EB3LYP or EMP2 = EHF + EMP2 correlation as applicable) and Hcorr is the total enthalpy correction to ε0.

2.4. Standard state, reference potential and deviations

In solution phase, 1 M reference state is used with the exception of water where 55.34 M is used [5458]. The correction for these are 3.02 mhartree and 3.80 mhartree, respectively. We do not use reference potential and ignore electrons in ΔG calculation.

Most variations of the blue bottle experiments take place in a alkaline solution with an exception of ascorbic acid system. The reactions are therefore considered in acid and alkaline conditions separately. Hydronium ion ([H3O+] = 1 M, pH = 0) and water are used in the calculation instead of hydrogen ion (H+) for reactions under acidic condition. Similarly, hydroxide ion ([OH] = 1 M, pH = 14) is used for reactions under alkaline condition. Examples are provided in table 1 to illustrate reactions under acidic and alkaline conditions.

Reduction potentials of many redox reactions are pH dependent due to deviation from standard condition. The deviation in terms of Gibbs energy (RT ln Q) is expressed in the last term of the Nernst equation,

E=EORTnFlnQ,
2.6

where Q is the reaction quotient. Since pH is not exactly 0 or 14 and the concentration of reactants are generally lower than 1 M, the ln Q consideration may be employed for detail analysis, especially when Ecell is close to zero.

2.5. Benchmarking

The mean unsigned errors for solution-phase and gas-phase EO of 52 selected reactions obtained at B3LYP/6-311++G** and MP2/cc-pVTZ are 0.86 V and 1.10 V, respectively. Figure 3 shows satisfactory linear relationships between EO obtained by the two methods (high R2 value but slope values slightly greater than unity). These benchmarking results confirm that B3LYP/6-311++G** yields acceptable results at a relatively small computational cost [59].

Figure 3.
Comparison of reduction potentials of reactions at B3LYP/6-311++G** and MP2/cc-pVTZ levels.

3. Results and discussion

3.1. Reduction potentials

Table 1 shows half-cell reduction potentials of all possible reactions in the blue bottle experiment. For comparison purpose, figure 4 shows the reduction potentials in acidic and alkaline conditions for the three groups of compounds, respectively. The first set of data is four- and two-electron oxygen reduction potentials on the first and fourth column of figure 4. The values are far from the literature values [30,31] but the trend that acidic potentials are higher than alkaline potentials is still preserved. The second set of data is the dye reduction potentials. Figure 4 shows that all dyes including indophenol may be conveniently oxidized by either two- or four-election ORR. The third set of data is the reducing agent reduction potentials. (Refer to Methodology section for grouping of reactions. Some compounds, for example, glucose, belong to two groups of reactions.)

Figure 4.
Reduction potentials of reactions in acidic and alkaline conditions calculated at B3LYP/6-311++G** level.

In general, the reduction potentials are quite similar for reactions in the same group under the same condition, and acidic potentials are higher than alkaline potentials. As reduction potentials of dyes and reducing agents are overlapping in figure 4, there are combinations of dyes and reducing agent that may or may not work in the blue bottle experiment.

3.2. Thermodynamic considerations

Positive cell potential from combination of the compounds can be found in figure 4 if the half-cell potentials decrease from left to right. For example, in alkaline condition, the classical blue bottle experiment may proceed via two- or four-electron oxygen reduction (approx. 4.7 V), with methylene blue as a catalyst (3.70 V) and glucose as a reducing agent (3.30 V). In acidic condition, the green version of the blue bottle experiment may proceed via two- or four-electron oxygen reduction (6.5–6.8 V), with methylene blue as a catalyst (5.18 V) and ascorbic acid as a reducing agent (5.10 V).

As our reduction potential considerations here are thermodynamic, the negative prediction (non-spontaneity for large negative value of Ecell) should be valid but the positive prediction (combination of dye and reducing agent make a blue bottle reaction for positive or close to zero value of Ecell) requires further verifications. To produce repeatable cycle of colour change, the rate of reduction of dye by reducing agent must be slower than the oxidation of dye by oxygen [14] and the direct oxidation of reducing agent by oxygen [5] should be minimal compared with the dye-catalysed reaction. Additional catalysts similar to the green version of the experiment [8,16] may be needed to make the reaction occur but it is beyond the scope of this study.

3.3. Implications and experimental confirmation

The following findings are made based on the calculated results and information in the literature. To support our claims, preliminary experiments to test some reducing agents were also carried out (see the electronic supplementary material).

An external file that holds a picture, illustration, etc.
Object name is rsos170708-i8.jpg

  1. Acyloin/enediol structure is necessary for the blue bottle reaction [5,7,8,11,16] and the main product is 1,2-diketone.
    •  Replacement of dextrose by simple aldehydes, e.g. benzldehyde, and alcohols, e.g. ethanol, does not result in repeatable colour change. The solution is blue and is not reduced to colourless over time.
    •  Replacement of dextrose by acetoin and dihydroxyacetone in the blue bottle experiment yields repeatable colour change and chemical patterns [18]. Phenol-derivatives such as catechol reacts with oxygen rapidly under alkaline conditions without a dye to produce a dark-coloured solution.
    •  GC/MS analysis of rapid blue bottle experiment confirms that benzil is a product of the reaction.
  2. The mechanism of reduction proceeds via enediolate formation [5] (deprotonation at carbon attached to OH).
    •  It is possible for ascorbic acid to deprotonate without the use of base and it is the only reducing agent for the experiment in acidic condition.
    •  Replacement of dextrose by α-hydroxy carboxylic acids such as citric acid and tartaric acid does not result in repeatable colour change. The solution is blue and it is not reduced to colourless over time. For citric acid, the tertiary alcohol group cannot be oxidized to ketone. For, tartaric acid, the negative charge upon deprotonation of –COOH group may make it difficult for the oxidation of the adjacent hydroxyl group.
      An external file that holds a picture, illustration, etc.
Object name is rsos170708-i9.jpg
    •  Replacement of dextrose by α-hydroxy lactones, for example, glucono delta-lactone did not change the colour of the dye from blue to colourless.
  3. Alternative reducing agents such as amino acids may also be used in the blue bottle experiment.
    •  Replacement of dextrose by amino acid and various food products yields repeatable colour change but some of the reactions are slow [18].
    •  In an iodine clock experiment [60], cysteine is also used successfully in place of ascorbic acid.

4. Concluding remarks and future work

Half-cell reduction potentials of oxygen, redox indicators and reducing agents have been investigated using DFT calculations. The results help us better understand the blue bottle reaction and guide us to focus the experiments only on a certain number of representative compounds and only for reaction that lead to a positive cell potential. The use of alternative reducing agents can help avoid side reactions that make the solution brown after a number of cycles [5] and increase solubility of the reducing agent in water which is a known issue for benzoin [7]. Possible future computational investigation includes prediction of pKa [39,56] and stability of intermediates and activated complexes of reduction reactions and prediction of the colour of redox dyes [61].

Supplementary Material

calculations.xlsx:

Supplementary Material

output.zip:

Supplementary Material

experiment.pdf:

Supplementary Material

omp2.txt:

Supplementary Material

ob3lyp.txt:

Acknowledgements

We thank Dr Andrew T. B. Gilbert for his guidance in Q-CHEM installation, Ms Punchalee Montriwat, Mr Rattha Noorat, and MUIC Science Division technicians for support in preliminary experimental work, Dr Yuthana Tantirungrotechai for helpful suggestions, and MUIC's Information Technology Section for the computing facility.

Data accessibility

The data supporting the findings of this study are available in the electronic supplementary material.

  •  Energies and graphs are summarized in calculations.xlsx.
  •  Preliminary experimental results are in experiment.pdf.
  •  Optimized geometries are in ob3lyp.txt and omp2.txt.
  •  All Q-Chem output files are in output.zip

Authors' contributions

T.L., the principal investigator and corresponding author, directed the research and wrote the manuscript. P.R. and C.A. assisted with calculation and experiment.

Competing interests

We declare we have no competing interests.

Funding

T.L. received funding from the Mahidol University International College (MUIC 001/2016) and the Institute for the Promotion of Teaching Science and Technology (IPST 08/2557).

References

1. Michaelis L. 1946. Demonstration of the existence of a free radical. J. Chem. Educ. 23, 317 (doi:10.1021/ed023p317)
2. Michaelis L. 1947. To the Editor. J. Chem. Educ. 24, 149 (doi:10.1021/ed024p149.3)
3. Campbell J. 1963. Kinetics—early and often. J. Chem. Educ. 40, 578–583. (doi:10.1021/ed040p578)
4. Dutton F. 1960. Methylene blue-reduction and oxidation. J. Chem. Educ. 37, A799 (doi:10.1021/ed037pA799.1)
5. Anderson L, Wittkopp SM, Painter CJ, Liegel JJ, Schreiner R, Bell JA, Shakhashiri BZ 2012. what is happening when the blue bottle bleaches: an investigation of the methylene blue-catalyzed air oxidation of glucose. J. Chem. Educ. 89, 1425–1431. (doi:10.1021/ed200511d)
6. Adamcikova L, Sevcik P 1998. The blue bottle experiment: simple demonstration of self-organization. J. Chem. Educ. 75, 1580 (doi:10.1021/ed075p1580)
7. Rajchakit U, Limpanuparb T 2016. Rapid blue bottle experiment: autoxidation of benzoin catalyzed by redox indicators. J. Chem. Educ. 93, 1490–1494. (doi:10.1021/acs.jchemed.6b00018)
8. Rajchakit U, Limpanuparb T 2016. Greening the traffic light: air oxidation of vitamin C catalyzed by indicators. J. Chem. Educ. 93, 1486–1489. (doi:10.1021/acs.jchemed.5b00630)
9. Chen PS. 1970. Resazurin-reduction and oxidation. J. Chem. Educ. 47, A335 (doi:10.1021/ed047pA335)
10. Chen PS. 1970. Oxidation and reduction of indigo carmine aided by benzoin. J. Chem. Educ. 47, A335 (doi:10.1021/ed047pA335)
11. Cook AG, Tolliver RM, Williams JE 1994. The blue bottle experiment revisited: how blue? how sweet? J. Chem. Educ. 71, 160–161. (doi:10.1021/ed071p160)
12. Vandaveer WR IV, Mosher M 1997. The blue bottle revisited. J. Chem. Educ. 74, 402 (doi:10.1021/ed074p402)
13. Hile L. 1998. The blue bottle revisited. J. Chem. Educ. 75, 1067 (doi:10.1021/ed075p1067.1)
14. Engerer SC, Cook AG 1999. The blue bottle reaction as a general chemistry experiment on reaction mechanisms. J. Chem. Educ. 76, 1519–1520. (doi:10.1021/ed076p1519)
15. Jonnalagadda S, Gollapalli N 2000. Kinetics of reduction of toluidine blue with sulfite—kinetic salt effect in elucidation of mechanism. J. Chem. Educ. 77, 506–509. (doi:10.1021/ed077p506)
16. Wellman WE, Noble ME, Tom AH 2003. Greening the blue bottle. J. Chem. Educ. 80, 537–540. (doi:10.1021/ed080p537)
17. Staiger FA, Peterson JP, Campbell DJ 2015. Variations on the ‘blue-bottle’ demonstration using food items that contain FD&C Blue# 1. J. Chem. Educ. 92, 1684–1686. (doi:10.1021/acs.jchemed.5b00190)
18. Limpanuparb T, Areekul C, Montriwat P, Rajchakit U 2017. Blue bottle experiment: learning chemistry without knowing the chemicals. J. Chem. Educ. 94, 730–737. (doi:10.1021/acs.jchemed.6b00844)
19. Adamcikova L, Sevcik P 1997. The blue bottle experiment and pattern formation in this system. Z. Naturforsch A 52, 650–654. (doi:10.1515/zna-1997-8-918)
20. Adamcikova L, Pavlikovia K, Secik P 1999. The methylene blue-d-glucose-O2 system: oxidation of d-glucose by methylene blue in the presence and the absence of oxygen. Int. J. Chem. Kinet. 31, 463–468. (doi:10.1002/(SICI)1097-4601(1999)31:6<463::AID-KIN8>3.0.CO;2-2)
21. Bees M, Pons A, Sorensen PG, Sagues F 2001. Chemoconvection: a chemically driven hydrodynamic instability. J. Chem. Phys. 114, 1932–1943. (doi:10.1063/1.1333757)
22. Pons A, Sagues F, Bees M, Sorensen PG 2000. Pattern formation in the methylene-blue-glucose system. J. Phys. Chem. B 104, 2251–2259. (doi:10.1021/jp9935788)
23. Pons A, Sagues F, Bees M, Sorensen PG 2002. Quantitative analysis of chemoconvection patterns in the methylene-blue-glucose system. J. Phys. Chem. B 106, 7252–7259. (doi:10.1021/jp014346n)
24. Jang NY, Won K 2014. New pressure-activated compartmented oxygen indicator for intelligent food packaging. Int. J. Food Sci. Technol. 49, 650–654. (doi:10.1111/ijfs.12310)
25. Schafer P, van de Linde S, Lehmann J, Sauer M, Doose S 2013. Methylene blue- and thiol-based oxygen depletion for super-resolution imaging. Anal. Chem. 85, 3393–3400. (doi:10.1021/ac400035k) [PubMed]
26. Kollner T, Rossi M, Broer F, Boeck T 2014. Chemical convection in the methylene-blue-glucose system: optimal perturbations and three-dimensional simulations. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 90, 053004 (doi:10.1103/PhysRevE.90.053004) [PubMed]
27. Winfree AT. 1984. The prehistory of the Belousov-Zhabotinsky oscillator. J. Chem. Educ. 61, 661 (doi:10.1021/ed061p661)
28. Briggs TS, Rauscher WC 1973. An oscillating iodine clock. J. Chem. Educ. 50, 496 (doi:10.1021/ed050p496)
29. Chen PS. 1974. Entertaining and educational chemical demonstrations. Camarillo, CA: Chemical Elements Publishing Company.
30. Song C, Zhang J 2008. Electrocatalytic oxygen reduction reaction. In PEM fuel cell electrocatalysts and catalyst layers (ed. Zhang J, editor. ), pp. 89–134. Berlin, Germany: Springer; (doi:10.1007/978-1-84800-936-3_2)
31. Liu K, Kattel S, Mao V, Wang G 2016. Electrochemical and computational study of oxygen reduction reaction on nonprecious transition metal/nitrogen doped carbon nanofibers in acid medium. J. Phys. Chem. C 120, 1586–1596. (doi:10.1021/acs.jpcc.5b10334)
32. Beheshti A, Norouzi P, Ganjali M 2012. A simple and robust model for predicting the reduction potential of quinones family; electrophilicity index effect. Int. J. Electrochem. Sci. 7, 4811–4821. (http://www.electrochemsci.org/papers/vol7/7054811.pdf)
33. Benassi R, Ferrarini P, Fontanesi C, Benedetti L, Paolucci F 2004. Calculated electron affinities and redox E 0 values of polypyridinic derivatives. J. Electroanal. Chem. 564, 231–237. (doi:10.1016/j.jelechem.2003.10.029)
34. Bushnell EA, Boyd RJ 2015. Assessment of several DFT functionals in calculation of the reduction potentials for Ni--, Pd--, and Pt--Bis-ethylene-1, 2-dithiolene and-diselenolene complexes. J. Phys. Chem. A 119, 911–918. (doi:10.1021/jp511350z) [PubMed]
35. Fareghi-Alamdari R, Zandi F, Keshavarz MH 2015. A new model for prediction of one electron reduction potential of nitroaryl compounds. Z. Anorg. Allg. Chem. 641, 2641–2648. (doi:10.1002/zaac.201500688)
36. Frontana C, Vazquez-Mayagoitia A, Garza J, Vargas R, Gonzalez I 2006. Substituent effect on a family of quinones in aprotic solvents: an experimental and theoretical approach. J. Phys. Chem. A 110, 9411–9419. (doi:10.1021/jp060836+) [PubMed]
37. Kelly CP, Cramer CJ, Truhlar DG 2005. SM6: a density functional theory continuum solvation model for calculating aqueous solvation free energies of neutrals, ions, and solute-water clusters. J. Chem. Theory Comput. 1, 1133–1152. (doi:10.1021/ct050164b) [PubMed]
38. Kelly CP, Cramer CJ, Truhlar DG 2007. Single-ion solvation free energies and the normal hydrogen electrode potential in methanol, acetonitrile, and dimethyl sulfoxide. J. Phys. Chem. B 111, 408–422. (doi:10.1021/jp065403l) [PMC free article] [PubMed]
39. Li J, Fisher CL, Chen JL, Bashford D, Noodleman L 1996. Calculation of redox potentials and pKa values of hydrated transition metal cations by a combined density functional and continuum dielectric theory. Inorg. Chem. 35, 4694–4702. (doi:10.1021/ic951428f)
40. Lister SG, Reynolds CA, Richards WG 1992. Theoretical calculation of electrode potentials: electron-withdrawing compounds. Int. J. Quantum Chem. 41, 293–310. (doi:10.1002/qua.560410206)
41. Llano J, Eriksson LA 2002. First principles electrochemistry: electrons and protons reacting as independent ions. J. Chem. phys. 117, 10 193–10 206. (doi:10.1063/1.1516786) [PubMed]
42. Mohammad-Shiri H, Ghaemi MARS, Akbari-Sehat A 2011. Computational and electrochemical studies on the redox reaction of dopamine in aqueous solution. Int. J. Electrochem. Sci. 6, 317–336. (http://www.electrochemsci.org/papers/vol6/6020317.pdf)
43. Namazian M, Norouzi P, Ranjbar R 2003. Prediction of electrode potentials of some quinone derivatives in acetonitrile. J. Mol. Struct. THEOCHEM 625, 235–241. (doi:10.1016/S0166-1280(03)00070-8)
44. Nematollahi D, Taherpour A, Jameh-Bozorghi S, Mansouri A, Dadpou B 2010. Experimental and computational study on electrochemical oxidation of catechols. Int. J. Electrochem. Sci. 5, 867–879. (http://www.electrochemsci.org/papers/vol5/5060867.pdf)
45. Paukku Y, Hill G 2011. Theoretical determination of one-electron redox potentials for DNA bases, base pairs, and stacks. J. Phys. Chem. A 115, 4804–4810. (doi:10.1021/jp201281t) [PubMed]
46. Phillips KL, Sandler SI, Chiu PC 2011. A method to calculate the one-electron reduction potentials for nitroaromatic compounds based on gas-phase quantum mechanics. J. Comput. Chem. 32, 226–239. (doi:10.1002/jcc.21608) [PubMed]
47. Reynolds CA. 1995. Density functional calculation of quinone electrode potentials. Int. J. Quantum Chem. 56, 677–687. (doi:10.1002/qua.560560605)
48. Roy LE, Jakubikova E, Guthrie MG, Batista ER 2009. Calculation of one-electron redox potentials revisited. Is it possible to calculate accurate potentials with density functional methods? J. Phys. Chem. A 113, 6745–6750. (doi:10.1021/jp811388w) [PubMed]
49. Ho J, Klamt A, Coote ML 2010. Comment on the correct use of continuum solvent models. J. Phys. Chem. A 114, 13 442–13 444. (doi:10.1021/jp107136j) [PubMed]
50. Ho J, Coote M, Cramer C, Truhlar D 2013. Theoretical calculation of reduction potentials. In Organic electrochemistry. (ed. Hammerich O and Speiser B), pp. 229–259. Boca Raton, FL: CRC Press; (doi:10.1201/b19122-8)
51. Marenich AV, Ho J, Coote ML, Cramer CJ, Truhlar DG 2014. Computational electrochemistry: prediction of liquid-phase reduction potentials. Phys. Chem. Chem. Phys. 16, 15 068–15 106. (doi:10.1039/C4CP01572J) [PubMed]
52. Marenich AV, Cramer CJ, Truhlar DG 2009. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396. (doi:10.1021/jp810292n) [PubMed]
53. Shao Y, et al. 2015. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol. Phys. 113, 184–215. (doi:10.1080/00268976.2014.952696)
54. Bryantsev VS, Diallo MS, Goddard WA III 2008. Calculation of solvation free energies of charged solutes using mixed cluster/continuum models. J. Phys. Chem. B 112, 9709–9719. (doi:10.1021/jp802665d) [PubMed]
55. Ho J, Ertem MZ 2016. Calculating free energy changes in continuum solvation models. J. Phys. Chem. B 120, 1319–1329. (doi:10.1021/acs.jpcb.6b00164) [PubMed]
56. Ho J, Coote ML 2010. A universal approach for continuum solvent pKa calculations: are we there yet? Theor. Chem. Acc. 125, 3 (doi:10.1007/s00214-009-0667-0)
57. Winget P, Cramer CJ, Truhlar DG 2004. Computation of equilibrium oxidation and reduction potentials for reversible and dissociative electron-transfer reactions in solution. Theor. Chem. Acc. Theory, Comput. Model. (Theor. Chim. Acta) 112, 217–227. (doi:10.1007/s00214-004-0577-0)
58. Camaioni DM, Schwerdtfeger CA 2005. Comment on ‘Accurate experimental values for the free energies of hydration of H+, OH, and H3O+. J. Phys. Chem. A 109, 10 795–10 797. (doi:10.1021/jp054088k) [PubMed]
59. Wiberg KB. 2004. Basis set effects on calculated geometries: 6-311++G** vs. aug-cc-pVDZ. J. Comput. Chem. 25, 1342–1346. (doi:10.1002/jcc.20058) [PubMed]
60. Wright SW, Phil R 2002. The vitamin C clock reaction. J. Chem. Educ. 79, 41 (doi:10.1021/ed079p41)
61. Fleming S, Mills A, Tuttle T 2011. Predicting the UV–vis spectra of oxazine dyes. Beilstein J. Org. Chem. 7, 432 (doi:10.3762/bjoc.7.56) [PMC free article] [PubMed]

Articles from Royal Society Open Science are provided here courtesy of The Royal Society