Free energy coupling coordinates

Fig. 2. Schematic plots outlining outer-shell free energy-reaction coordinate profiles for the redox couple O + e R on the basis of the hypothetical two-step charging process (Sect. 3.2) [40b]. The y axis is (a) the ionic free energy and (b) the electrochemical free energy (i.e. including free energy of reacting electron), such that the electrochemical driving force, AG° = F(E - E°), equals zero. The arrowed pathways OT S and OTS represent hypothetical charging processes by which the transition state, T, is formed from the reactant.

Here, 7 runs over all simulations and k, I run over all bins. These equations can be solved iteratively, assuming an initial set oi fj (e.g., fj =1), then calculating p°i from Eq. (34) and updating Ihe fj by Eq. (35), and so on, until thep°i no longer vary, i.e., the two equations are self-consistent. Erom the p°i = P(qt, sp and Eq. (27), one then obtains the free energy of each bin center (q, sp. Error estimates are also obtained [46]. The method can be applied to a one-dimensional reaction coordinate or generalized to more than two dimensions and to cases in which simulations are run at several different temperatures [46]. It also applies when the reaction coordinates are alchemical coupling coordinates (see below and Ref. 47). 

The left (solid) parabolic curve represents the oxidized state, the right one, the reduced state. Let us assume that the system is initially at the oxidized state (left curve). When the interaction metal-reaction species is small, the electronic coupling between is small and the system may oscillate many times on the left parabolic curve (ox) before it is transferred to the curve on the right (red). On the other hand, if the interaction is strong, the free energy should no longer be represented by the two solid curves in the intermediate region of the reaction coordinate, but rather, by the dashed... 

Fig. 1. Schematic one-dimensional cross section through the Gibbs free energy surface G(R) of a spin-state transition system along the totally symmetric stretching coordinate. The situation for three characteristic temperatures is shown (B = barrier height, ZPE = zero-point energy, 28 = asymmetry parameter, J = electronic coupling parameter, AG° = Gh — GJ...

Recently, much attention has been paid to the investigation of the role of this interaction in relation to the calculations for adiabatic reactions. For steady-state nonadiabatic reactions where the initial thermal equilibrium is not disturbed by the reaction, the coupling constants describing the interaction with the thermal bath do not enter explicitly into the expressions for the transition probabilities. The role of the thermal bath in this case is reduced to that the activation factor is determined by the free energy in the transitional configuration, and for the calculation of the transition probabilities, it is sufficient to know the free energy surfaces of the system as functions of the coordinates of the reactive modes. 

Fig. 2.5. Possible applications of a coupling parameter, A, in free energy calculations, (a) and (b) correspond, respectively, to simple and coupled modifications of torsional degrees of freedom, involved in the study of conformational equilibria (c) represents an intramolecular, end-to-end reaction coordinate that may be used, for instance, to model the folding of a short peptide (d) symbolizes the alteration of selected nonbonded interactions to estimate relative free energies, in the spirit of site-directed mutagenesis experiments (e) is a simple distance separating chemical species that can be employed in potential of mean force (PMF) calculations and (f) corresponds to the annihilation of selected nonbonded interactions for the estimation of e.g., free energies of solvation. In the examples (a), (b), and (e), the coupling parameter, A, is not independent of the Cartesian coordinates, x. Appropriate metric tensor correction should be considered through a relevant transformation into generalized coordinates...

Rao and Singh32 calculated relative solvation free energies for alkyl- and tetra-alkylammonium ions using same conditions as described before for neutral molecules (except, atomic partial charges were not scaled for ions). The values obtained with coordinate coupling were in better agreement with... 

In the iterative approach using a constant bias F, the free energy barrier is reduced each successive iteration and therefore produces complete sampling of important configurations along the coupling parameter (i.e., reaction coordinates). Furthermore, more complicated umbrella potentials (e.g. see Equation 21) may also be applied with the iterative procedure. 

Reaction Free Energy Profiles Using Free Energy Perturbation and Coordinate Coupling Methodologies Analysis of the Dihydrofolate Reductase Catalytic Mechanism... 

We carried out some measurements some years ago in order to measure the transfer coefficient in DMF for this couple at a mercury electrode [Hush, N. S. Dyke, J. M. J. Electroanal. Interfac. Electrochem. 1974, 53, 253]. The deviation from one-half for the value so obtained was consistent with reasonable values for the free energy of activation of ligand exchange in the inner-coordination spheres. 

Fig. 1. The Marcus parabolic free energy surfaces corresponding to the reactant electronic state of the system (DA) and to the product electronic state of the system (D A ) cross (become resonant) at the transition state. The curves which cross are computed with zero electronic tunneling interaction and are known as the diabatic curves, and include the Born-Oppenheimer potential energy of the molecular system plus the environmental polarization free energy as a function of the reaction coordinate. Due to the finite electronic coupling between the reactant and charge separated states, a fraction k l of the molecular systems passing through the transition state region will cross over onto the product surface this electronically controlled fraction k l thus enters directly as a factor into the electron transfer rate constant...

---