TY - JOUR
T1 - Diffusive dynamics on potential energy surfaces
T2 - Nonequilibrium CO binding to heme proteins
AU - Agmon, Noam
AU - Rabinovich, Savely
PY - 1992
Y1 - 1992
N2 - Theory and practice of reaction dynamics on two-dimensional potential energy surfaces is investigated. Nonequilibrium multidimensional barrier crossing, occurring when the initial density is located near the ridgeline separating reactants and products, is treated by solving the time-dependent Smoluchowski equation as a function of diffusion anisotropy. For a locally separable potential, and slow relaxation in the perpendicular mode, the problem reduces to a one-dimensional Smoluchowski equation with a sink term. It may be further approximated as a first-order reaction with a time-dependent rate coefficient. These approximations are compared with exact two-dimensional propagations on a potential surface representing CO binding to α heme. The intermediate-time power-law decay of the survival probability is analyzed with the aid of the above approximations. The power also shows some kind of critical behavior near the isotropic diffusion limit, where these approximations are no longer valid. For fast relaxation, a nonmonotonic survival probability is observed. The long time decay of the survival probability is governed by the equilibrium rate coefficient. We calculate its anisotropy dependence and compare it with two asymptotic expansions for the lowest eigenvalue of the Smoluchowski operator-for the one-dimensional sink-Smoluchowski operator and the fully two-dimensional operator. The difference between the fast relaxation limit of these one- and two-dimensional equations may provide a quantitative explanation for previous problems in extrapolating to high temperatures using the Agmon-Hopfield model. The implications of these results to heme protein dynamics are discussed and new experiments are proposed.
AB - Theory and practice of reaction dynamics on two-dimensional potential energy surfaces is investigated. Nonequilibrium multidimensional barrier crossing, occurring when the initial density is located near the ridgeline separating reactants and products, is treated by solving the time-dependent Smoluchowski equation as a function of diffusion anisotropy. For a locally separable potential, and slow relaxation in the perpendicular mode, the problem reduces to a one-dimensional Smoluchowski equation with a sink term. It may be further approximated as a first-order reaction with a time-dependent rate coefficient. These approximations are compared with exact two-dimensional propagations on a potential surface representing CO binding to α heme. The intermediate-time power-law decay of the survival probability is analyzed with the aid of the above approximations. The power also shows some kind of critical behavior near the isotropic diffusion limit, where these approximations are no longer valid. For fast relaxation, a nonmonotonic survival probability is observed. The long time decay of the survival probability is governed by the equilibrium rate coefficient. We calculate its anisotropy dependence and compare it with two asymptotic expansions for the lowest eigenvalue of the Smoluchowski operator-for the one-dimensional sink-Smoluchowski operator and the fully two-dimensional operator. The difference between the fast relaxation limit of these one- and two-dimensional equations may provide a quantitative explanation for previous problems in extrapolating to high temperatures using the Agmon-Hopfield model. The implications of these results to heme protein dynamics are discussed and new experiments are proposed.
UR - http://www.scopus.com/inward/record.url?scp=36449002985&partnerID=8YFLogxK
U2 - 10.1063/1.463500
DO - 10.1063/1.463500
M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???
AN - SCOPUS:36449002985
SN - 0021-9606
VL - 97
SP - 7270
EP - 7286
JO - The Journal of Chemical Physics
JF - The Journal of Chemical Physics
IS - 10
ER -