banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Redox potentials and for benzoquinone from density functional theory based molecular dynamics
Rent this article for


Image of FIG. 1.
FIG. 1.

Decomposition of the reduction of quinone to hydroquinone in reduction (horizontal arrows) and protonation half reactions (vertical arrows). The quinone, semiquinone, and hydroquinone are indicated by Q, HQ, and , respectively. Diagonal arrows correspond to (de)hydrogenation (proton coupled electron transfer). Note that the usual dot indicating that a species is open shell (a radical) has been suppressed.

Image of FIG. 2.
FIG. 2.

Accumulative averages of vertical energy gaps as a function of MD time for the four oxidation reactions of Fig. 1: (a) ; (b) ; (c) ; and (d) . The three curves for each reaction are the result for three different PESs as defined by the coupling parameter [see Eq. (9)]. In the order of decreasing value of the energy gap the coupling parameter values are (reactant), (mixed state), and (product). Each curve is followed by the overall time average of the vertical energy gap in eV (see also Table III).

Image of FIG. 3.
FIG. 3.

Accumulative averages of vertical energy gaps as a function of MD time for the deprotonation and dehydrogenation reactions of Fig. 1. Labeling is continued from Fig. 2 (see also Table III): (e) ; (f) ; (g) ; (h) ; (i) ; (j) . For reactions (e), (f), (i), and (j) only the results for the two end states and the half way mixed state are given similar to the reactions (a)–(d) in Fig. 2. For reactions (g) and (h) a more closely spaced set of coupling parameter values has been used, namely, , for reaction (h) an additional state at . The gap energies are strictly monotonously decreasing in [see plot of charging curves for reactions (g) and (h) in Fig. 5]. Energies at the end of the accumulative average plots are the final values used in the calculation of the free energies.

Image of FIG. 4.
FIG. 4.

Free energy diagram for quinone reduction. The numbers alongside the arrows give the calculated reaction free energies (units in eV) in the reverse direction, i.e., hydroquinone oxidation to quinone. The free energy changes from experimental (Ref. 10) NHE reduction potentials and are given in parentheses. The numbers in bold are the free energy changes of four full reactions. The three red (blue) arrows indicate three alternative pathways for the , respectively reactions. The energies marking the arrows are the corresponding DFTMD reaction free energies which according to Hess’s law should be the same for each reaction.

Image of FIG. 5.
FIG. 5.

Correlation between vertical energy gap and solvation structure for half reactions [panels (a) and (d)], [(b) and (e)], and [(c) and (f)]. Black squares in (a)–(c) give the vertical energy gaps as determined by averaging over MD trajectories at the indicated values of the coupling parameter . Red circles are the corresponding coordination numbers obtained by integrating over the first solvation shell of the rdf between quinone O atoms and solvent H atoms (denoted by ). The rdfs at the end points (black solid curves), (blue dashed-dotted curves), and one intermediate point (red dashed curves) are given in panels (d)–(f) on the right. The blue dashed line in (a) is the linear fit of the vertical energy gap as a function of coupling parameter with the regression coefficient .


Generic image for table
Table I.

Force constants and structural parameters for the restraining potentials defined in Eq. (19). , , and are the force constants for bonds, bends, and torsions, respectively, in atomic unit. , , and are the corresponding equilibrium values in atomic unit for bond lengths and radians for angles and dihedrals. For weak acids . For strong acids is applied to the OH bond only (for both O–H bonds are restrained).

Generic image for table
Table II.

Vibrational frequencies , vibrational temperatures , moments of inertia , and rotational temperatures used in the correction for restraining potentials and zero point motion (see Sec. II E). The molecule is the hydronium ion with one proton replaced by a dummy of the same mass but without charge. Similarly is a quinone with the acid proton replaced by a dummy proton. The three vibrational frequencies given for are the normal modes of the dummy atom with all of the quinone anion atoms fixed.

Generic image for table
Table III.

DFTMD free energy changes , variances of energy gaps , and reorganization energies for the half reactions of Fig. 1. Subscripts LR and TI indicate the numerical method used to calculate free energy changes. LR stands for linear response (end point) approximation [Eq. (13)] and TI for thermodynamic integration using additional intermediate coupling parameter states (the number of integration points can be deduced from Figs. 2 and 3). Subscripts 0 and 1 denote the value of the coupling parameter corresponding to begin and end states [see Eq. (9)]. The reorganization energy has been calculated according to Eq. (15) allowing for nonlinear solvent response. All simulations have been performed in a cubic periodic cell box with edge except the results of redox reactions in parentheses, which were calculated in a smaller box . Energies are in eV and variances are in .

Generic image for table
Table IV.

Full reactions for the determination of the oxidation free energy versus NHE [(a)-(d)], Brønsted acidity [(e)–(h)], and dehydrogenation free energies [(i) and (j)] of the half reactions of Table III. The data marked repeat the free energies of the column of Table III. are the experimental free energies as collected in Ref. 10 (same data are given again in Fig. 4). is the free energy cost of deprotonating a hydronium ion estimated by aligning simulation and experimental free energies. The average obtained in this way from the acid dissociation reactions [(e)–(h)] have been used to compute the oxidation free energies and acidities denoted by (see Sec. III C for more details). Figures in parentheses include the corrections for the restraining potentials and zero point motion (see Sec. II E). All energies are in eV. The mean deviation from experiment of the restraint corrected oxidation free energies is −0.52 eV.

Generic image for table
Table V.

Quinone energetics in vacuum as computed in this work using BLYP compared with the results taken from the literature. The DFT functional used in Refs. 43, 48, and 49 is B3LYP and B3PW91 in Ref. 97. The unit of energy is eV.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Redox potentials and pKa for benzoquinone from density functional theory based molecular dynamics