Description Usage Arguments Details Value References See Also Examples
Computes the non-adiabatic activation energy ( ΔEa ) involved in the calculation of the Landau-Zener and Marcus charge transfer (CT) rates using the following expression:
ΔEa = (λ + ΔE0 + ΔEField)^2/(4λ)
1 | activationFreeEnergy(lambda, dE0 = 0, dEField = 0)
|
lambda |
a scalar, a vector, a matrix or an array containing the total reorganization energies (in eV) used to calculate the activation energy. |
dE0 |
a scalar, a vector, a matrix or an array containing the site energy differences (in eV) used to calculate the activation energy. By default self-exchange reactions are considered ( |
dEField |
a scalar, a vector, a matrix or an array containing an additional contribution to the site energy differences due to an external electric field (in eV). By default no electric field is applied ( |
The arguments of these function can be scalars, vectors, matrices or arrays. Mixing scalar values with vectors, matrices or arrays is allowed but in all other cases the arguments must have the same dimensions and lengths.
Depending on the dimension of the objects passed to the function a scalar, a vector, a matrix or an array containing the activation energies (in eV) is returned.
H. Oberhofer and J. Blumberger, Physical Chemistry Chemical Physics, 14, 13846–13852, 2012
B. S. Brunschwig, J. Logan, M. D. Newton and N. Sutin, Journal of the American Chemical Society, 102, 5798, 1980
M. D. Newton and N. Sutin, Annual Review of Physical Chemistry, 35, 437, 1984
M. D. Newton, Chemical Reviews, 91, 767, 1991.
A. Nitzan, Chemical Dynamics in Condensed Phases, Oxford University Press, 2006
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 | ## Produce a map of the non-adiabatic free energy barrier and
## two other maps of the non-adiabatic free energy barrier corrected by
## the deviation between the adiabatic and non-adiabatic potential energy curves at:
## - the crossing point (dEad = dEna - abs(J))
## - both the crossing point and the minimas of the curves (dEad = dEna - Delta)
## We assume that the reorganization energies are the following:
lambdaI <- 0.14 # internal reorganization energy in eV
lambdaS <- 36E-3 # external reorganization energy in eV
lambda <- lambdaI + lambdaS
N <- 301
J <- seq( 0 , 65,length.out=N)*1E-3 # eV
dE <- seq(-0.5,0.5,length.out=N) # eV
G <- expand.grid(J, dE)
J <- G[,1]
dE <- G[,2]
dEna <- activationFreeEnergy(lambda = lambda, dE0 = dE, dEField = 0)
dEad <- dEna - adiabaticCorrection(J = J, lambda = lambda, dE0 = dE, dEField = 0)
dEadCross <- dEna - J
dEna <- matrix(dEna , nrow = N, ncol = N)
dEad <- matrix(dEad , nrow = N, ncol = N)
dEadCross <- matrix(dEadCross, nrow = N, ncol = N)
addAxis <- function(bottom = TRUE, left = FALSE, above = FALSE, right = FALSE){
useless <- lapply(1:4,axis, labels=FALSE)
if(bottom) axis(1, labels = TRUE)
if(left ) axis(2, labels = TRUE)
if(above ) axis(3, labels = TRUE)
if(right ) axis(4, labels = TRUE)
if(bottom) mtext(side=1,line=1.2, text=expression( abs(J)/eV), cex=par("cex"))
if(left ) mtext(side=2,line=1.2, text=expression(Delta*E/eV), cex=par("cex"))
if(right ) mtext(side=4,line=1.2, text=expression(Delta*E/eV), cex=par("cex"))
box()
}
layout(matrix(1:3, ncol=3))
par(cex=2, lwd=1.5, pty="s", mgp=c(1.1,0.1,0), tck=0.02, mar=rep(0.7,4), oma=rep(2,4))
contour(unique(J), unique(dE), dEna , levels = seq(-0.1,0.6,0.1),
xaxt="n", yaxt="n", labcex=2)
addAxis(TRUE, TRUE, FALSE, FALSE)
title(expression(Delta*E[na]), line=1)
contour(unique(J), unique(dE), dEadCross, levels = seq(-0.1,0.6,0.1),
xaxt="n", yaxt="n", labcex=2)
addAxis(TRUE, FALSE, FALSE, FALSE)
title(expression(Delta*E[na]-J), line=1)
contour(unique(J), unique(dE), dEad , levels = seq(-0.1,0.6,0.1),
xaxt="n", yaxt="n", labcex=2)
addAxis(TRUE, FALSE, FALSE, TRUE)
title(expression(Delta*E[na]-Delta), line=1)
|
Loading required package: parallel
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.