Bridge mediated two-electron transfer reactions: Analysis of stepwise and concerted pathways

Size: px
Start display at page:

Download "Bridge mediated two-electron transfer reactions: Analysis of stepwise and concerted pathways"

Transcription

1 JOURAL OF CHEMICAL PHYSICS VOLUME 120 UMBER 9 1 MARCH 2004 Bridge mediated two-electron transfer reactions: Analysis of stepwise concerted pathways E. G. Petrov Bogolyubov Institute for Theoretical Physics ational Academy of Sciences of Ukraine UA Kiev Ukraine V. May Institüt für Physik Humboldt Universität zu Berlin D Berlin Germany Received 26 September 2003; accepted 4 December 2003 A theory of nonadiabatic donor (D) acceptor (A) two-electron transfer TET mediated by a single regular bridge (B) is developed. The presence of different intermediate two-electron states connecting the reactant state D BA with the product state DBA results in complex multiexponential kinetics. The conditions are discussed at which a reduction to two-exponential as well as single-exponential kinetics becomes possible. For the latter case the rate K TET is calculated which describes the bridge-mediated reaction as an effective two-electron D A transfer. In the limit of small populations of the intermediate TET states D B A DB A D BA DB A K TET is obtained as a sum of the rates K (step) TET K (sup) TET. The first rate describes stepwise TET originated by transitions of a single electron. It starts at D BA reaches DBA via the intermediate state D BA. These transitions cover contributions from sequential as well as superexchange reactions all including reduced bridge states. In contrast a specific two-electron superexchange mechanism from D BA to DBA defines K (sup) TET. An analytic dependence of K (step) TET K (sup) TET on the number of bridging units is presented different regimes of D A TET are studied American Institute of Physics. DOI: / I. ITRODUCTIO Two-electron transfer TET processes processes which incorporate even more electrons are common for enzyme regulated reactions in biological systems. For example multielectron transfer could be clearly attributed to different enzyme complex mediated oxidation-reduction reactions. TET reduction process have been reported for different types of reductases 12 oxidases 3 5 hydrogenases. 6 8 All these reactions involve numerous intermediate states proceed in a much more complex way as single-electron transfer SET processes. Protein film voltametry represents an experimental technique which enables direct observation of such electron transfer reactions In biological systems SET TET higher multielectron reactions incorporate charge motion along specific structures like peptide chains DA strs as well as complete redox chains. In the various types of fumarate reductases for example the redox chains are given either by the chain of hemes or by Fe S clusters. 2 Such structures form molecular bridges mediate the electron transfer across large spatial distances. To account for the influence of molecular bridges on electron transfer reactions represents a long-sting task of theoretical chemical physics. We consider it as a particular challenge to formulate a theory of bridge-mediated multielectron transfer reactions In recent studies on nonadiabatic bridge-mediated SET the derivation of kinetic equations could be demonstrated as well as rate constants which account for different transfer mechanisms sequential superexchange SET. As a particular result conditions could be specified for which the overall SET rate follows by additive contributions from the sequential the superexchange SET mechanism. 20 Turning back to TET reactions it has to be considered as a continuous problem to correctly characterize the different mechanisms contributing to the observed overall processes. At present there are some basic results on the description of TET in polar liquids where the transfer directly proceeds Current theoretical studies on TET reactions in biological systems concentrate on electronic structure calculations of e.g. the active centers of enzymes their dependency on the position of the transferred electrons 47 or on possible TET pathways To estimate the rate constants characterizing the TET process so far phenomenological versions of the Marcus theory 16 have been used. In particular the matrix element responsible for the concerted two-electron hopping transition has been introduced as a phenomenological parameter. 17 Such treatments of TET reactions indicates the need for a theory which is founded on nonequilibrium quantum statistics in this way gives a correct description of the concerted as well as stepwise TET pathways mediated by the molecular bridge. In the following we will present such a theoretical description of nonadiabatic long-range TET reactions. The reaction to be investigated starts with two electrons located at the donor (D). Both are transferred through a chain of bridging molecules (B) thereafter they are captured at the acceptor (A). The whole course of the reaction will be covered by a description which generalizes earlier approaches derived for the study of nonadiabatic D A SET for an over /2004/120(9)/4441/16/$ American Institute of Physics

2 4442 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May view see e.g. Refs the textbooks It is just the nonadiabatic character of the transfer process which enables us to utilize a coarse-grained approach to simplify the necessary rate equations considerably The paper is organized as follows: In the next section we set up a model for the TET process covering all two-electron configurations of the DBA system. This is continued by the description of a particular density matrix approach which enables us to derive rate equations respective rate constants for the nonadiabatic TET. The latter include expressions originated by the sequential the superexchange electron transfer mechanism. A reduction of the TET kinetics to two-exponential single-exponential ones is carried out in Sec. III. The conditions to be fulfilled to have this simplified behavior are specified. Section IV is devoted to an analysis of the bridge-length dependence of the stepwise the concerted part of the overall D A TET rate. A detailed discussion of those mechanisms originating the bridgeassisted D A TET is given in the concluding Sec. V. Computational details are displaced to the Appendix. II. COARSE-GRAIED DESCRIPTIO OF TET As it is well-known nonadiabatic SET takes place across different sites of electron localization. Due to the weak intersite electronic coupling the transfer proceeds against the background of considerably faster intrasite relaxation processes the whole mechanism is known as the sequential SET. If the SET occurs between a D A center connected by a molecular bridge long-range superexchange D A SET may represent an alternative charge transfer mechanism Since for both mentioned mechanisms the characteristic time scale t of the overall electron motion exceeds the characteristic time rel of intrasite relaxation the SET kinetics can be described by a simplified set of coarse-grained rate equations see Refs The sequential superexchange mechanisms enter via respective rate expressions. We must expect similar conditions for nonadiabatic TET reactions if the characteristic transfer time TET substantially exceeds rel. Therefore the inequality rel TET 1 represent the precondition to employ a coarse-grained description of TET processes on a time scale t rel. A. Basic model Theoretical studies on TET processes in a DA complex dissolved in a polar liquid are based on the use of three distinct two-electron configurations. The TET starts at the initial configuration D A with the two electrons at the D site. It is followed by the intermediate configuration D A the TET ends at the final configuration DA with the two electrons located at the A site This course of TET directly corresponds to the stepwise mechanism D A D A DA which has to be confronted with the concerted TET D A DA. Thus the intermediate electronic configuration D A is of basic importance for the TET process between the contacting redox centers. FIG. 1. Pathway scheme of bridge-assisted TET reactions. Thick dashed lines indicate the sequential single-electron pathways. Sole double solid lines display single-electron two-electron superexchange pathways respectively. The thin dashed lines correspond to the pathways via twofold reduced bridging states. In the present paper we consider a generalization of this TET model to the case of long-range TET mediated by a bridge of units connecting the D the A. According to a possible large number of bridge units the number of accessible two-electron configurations become also large cf. Fig. 1. Let us start with those states in which the bridge has not been reduced. To these types of states belongs the initial state of the TET reaction DD B 1 B 2 B m B A 2 the intermediate state ID B 1 B 2 B m B A 3 the final state of the reaction ADB 1 B 2 B m B A. 4 Furthermore there exist two categories of states including a singly reduced bridge B m D B 1 B 2 B m B A B ndb 1 B 2 B n B A. 6 Every category has to be subdivided into different states. If both transferred electrons in the bridge we arrive at the manifold of 2 doubly reduced bridge states note that a double occupation of a single bridge unit is not forbidden B mn D B 1 B 2 B m B n B A. 7 The possible TET routes are depicted in Fig. 1. Below however we will exclusively restrict our studies to TET pro- 5

3 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4443 cesses where doubly reduced bridge states do not take part. In particular this would be the case for TET through short bridges where the Coulomb repulsion between the transferred electrons strongly increases the energetic position of the states B nm relative to all other TET states see also the discussion in Ref. 35. Resulting from this restriction the mechanisms of TET named stepwise TET from now on can be characterized as follows see Fig. 1. It consists of two sequential SET routes each related to a single-electron hopping transitions between nearest neighbor bridge units. The first type of a sequential SET route i.e. the sequence of reactions: D B 1 B A D B 1 B A D B 1 B 2 B A D B 1 B A D B 1 B A represents an electron transfer from the initial state D to the intermediate state I. Then the second type of sequential SET let move the second electron from the D to the A: D B 1 B A DB 1 B A DB 1 B 2 B A DB 1 B A DBA. In the scheme of state Eqs. 2 6 the overall TET is completed by the transition from the intermediate state I to the final state A. Apparently the transitions D I I A can also proceed directly via the single-electron superexchange mechanisms full lines in Fig. 1. It will be the particular aim of the following derivations to present an approach which accounts for stepwise TET as well as includes all possible superexchange mechanism. Beside the described superexchange we have also to expect the so-called concerted TET leading directly from state D to state A double solid line in Fig. 1. In order to derive rate equations which incorporate the different types of TET pathways in a first step we have to set up a respective Hamiltonian. We abbreviate all twoelectron states Eqs. 2 6 bym (MDB m IB n A) get TET Hamiltonian as H TET H 0 V with the zero-order part 8 H 0 M H M MM 9 where the H M denote the vibrational Hamiltonian which belong to the states M. The couplings transfer integrals H M between the different two-electron states are included in the interaction Hamiltonian V M 1 M H M M. 10 However only those H M are included for which the state M can be translated into one another by a singleelectron exchange. Let us introduce the eigenstates energies of the Hamiltonian H M as M M E(M M ) respectively where M denotes the quantum number of the vibrational state belonging to the electronic state M. In the case that all vibrational quantum numbers of the electronic state M are zero we will denote the respective energy by E M.Ifan energetic bias is present we assume E m E m1 Ẽ n Ẽ n1 may write for a regular bridge E m E B m1 Ẽ n Ẽ B n1. 11 There is a variety of reasons for. Its presence may result from externally applied as well as internal intermembrane electric fields acting along the bridge. However the bias can also be introduced by slightly changing the chemical structure of the bridge units or of the groups surrounding the bridge. Moreover non-condon effects are neglected the coupling matrix elements are approximated by V MM M M M 12 where M denotes the vibrational overlap integral. The transfer matrix elements M M MVˆ tr characterize the transitions between the electronic states M Eqs. 2 6 of the whole DBA system. We shall specify these matrix elements utilizing the tight binding model where the transfer operator Vˆ tr describes the single-electron transitions. It allows to express M M via the couplings V between the molecular orbitals related to the neighboring bridge sites. For instance if MD B 1 then M DB1 V D1 if MI B 1 then M IB 1V D1. Analogously one obtains M IB V A M AB V A.In the case of a regular bridge we set for all sites in the bridge M Bm B m1 M B n B V n1 mm1v B. All bridge-assisted electron transfer processes caused by the single-electron couplings are depicted in Fig. 2. The presence of two single-electron stepwise pathways can be clearly identified. The direct couplings T DI T IA are the result of a superexchange mechanism cf. the Appendix C Eqs. C1 C2. The energetic position of all considered TET states are represented in Fig. 3. B. Derivation of kinetic equations Our considerations on TET reactions will be based on the fact that inequality 1 is fulfilled. The related fast vibrational relaxation within every electronic state of the TET results in electron vibrational level broadening 18 dephasing Here we denote the broadening of the energy level E(M M )by(m M ). It has been discussed at length in Ref. 18 how to compute the electron vibrational dynamics in a case characterized by Eq. 1. In a first step one has to set up equations of motion for the electron vibrational density matrix MM tm Mt 13 where (t) is the related reduced density operator. The influence of an additional heat bath causing fast vibrational relaxation has been taken into consideration by energy relaxation dephasing rates. ext one changes to a time region which is large compared to the fast vibrational relaxation processes coarse-graining approximation. This change is

4 4444 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May FIG. 2. Energetic scheme of the DBA states participating in the TET. The single-electron transitions along the first (D BA) (D B 1 A) (D B A) (D BA ) the second (D BA ) (DB 1 A ) (DB A ) (DBA ) stepwise pathways are shown panels a b respectively. Wavy lines indicate the fast relaxation within the various electronic state. achieved in reducing the exact density matrix equations to the following approximate versions valid for the diagonal density matrix elements: t M M M M t i M MVt M M 14 the off-diagonal density matrix elements note in particular M) MM t M MVt. 15 M M The reduced set of equations of motion assumes that the time dependence of the diagonal density matrix elements is exclusively determined by the coupling to other TET states whereas relaxation processes do not contribute they are just finished at every time step. In contrast the time derivative has been removed in the equation of motion for the offdiagonal density matrix elements. This approximation results in a determination of the time dependence mainly by the dephasing rates is justified by the rapidness of the dephasing. The related dephasing rates enter the equation via (M M )E(M M )E( )i((m M ) ( )). Equations will be used to derive rate equations for the electronic state populations P M t m P MM t m MM M M t 16 where the transition rates are obtained as expansions with respect to the transfer integrals V M. In order to derive the related perturbation theory we first introduce the projector PA m A MM M M M M M M 17 the orthogonal complement Q1P. The latter projects any operator A defined in the electron vibrational state space on the part which has only off-diagonal matrix elements. A projection on the part of A including only diagonal matrix elements with respect to the states M M is

5 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4445 but in any case they only contain an even number of interactions. Taking the matrix elements of Eq. 21 we obtain for the electron vibrational state populations t P M M t M M R M M P t. 24 According to the coarse graining approximation there is vibrational equilibrium established at every electronic state at every time t rel. Therefore we can assume the following thermal equilibrium relation to be fulfilled: FIG. 3. Energetic positions E M E(M0 M ) of the states involved in the TET together with respective energy gaps. achieved by P. Furthermore we change from the matrix notation Eqs to the operator notation t PtiPL VQt QttiQRL V PtQt ote the introduction of L V (1/)V.... It allows to replace L V (t) byl V Q(t) in the first equation. The resolvent superoperator R 0 d S 0 () could be introduced in the second equation since it generates an expression as in Eq. 15 when changing to matrix elements. The newly defined time-propagation superoperator S 0 is determined by the zeroorder Liouvillian L 0 (1/)H 0... as well as by a part describing dissipation. To get the required rate equations in a next step we have to derive a master equation for the diagonal part P of the density operator. This is easily achieved by formally solving Eq. 19 as Qti1iQRL V Pt inserting the result into Eq. 18. It follows 20 t PtRPt. 21 Here R denotes the superoperator which describes any type of TET. Its expansion with respect to the transfer coupling 10 RR 2 R 4 R 6 22 enables one to arrange the different TET processes in a proper order. The various terms of the expansion can be written as P MM t/p MM texpem MEM M /k B T. 25 Changing to the total TET state population Eq. 16 wemay set see e.g. Refs P MM twem M P M t 26 where WEM M Z M 1 expem M /k B T Z M 1 M expem M /k B T 27 is the statistical weight of the electron-vibrational state M M. oting Eq. 26 we obtain from Eq. 24 the set of balancelike equations with Ṗ M t K M P t 28 K M M WE M M R M M. 29 ow noting the fact that the symbols M indicate electronic states 2 6 related to the whole DBA system we obtain see Appendix A the set of coupled rate equations A8 A16 with respective sequential superexchange rate constants. III. OVERALL TET RATES The concrete form of the balancelike equations 28 as given by the set of rate equations Eqs. A8 A16 reflects the rather complicated kinetics of the TET process also shown in scheme a of Fig. 4. Generally such kinetics has a multiexponential character for each of the 23 state populations P M (t) (MDIAB 1...B B 1...B ) R 2k i 2k 0 d2k1 0 d2k2 0 d1 22 P M tp M r1 B M (r) expk r t. 30 PL V QS 0 2k1 L V QS 0 2k2 L V QS 0 1 L V 23 The steady-state populations P M ()lim s 0 (sf M (s)) as well as the overall transfer rates K r can be deduced from the

6 4446 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May If the total bridge state population P B (t) remains small during the TET reaction the solution 31 describes TET only between the electronic states D I A. Just the inequality P B t m1 P m tp m t1 32 represents a necessary sufficient condition for the reduction of the multiexponential TET kinetics to a twoexponential process. If t the inequality 32 can be rewritten as FIG. 4. Complete kinetic scheme of the bridge-assisted TET process part a. For a small bridge-state population the kinetics is reduced to the transitions between the three electronic states D I A the related transfer rates are also shown part b. If the population of the intermediate state I becomes also small the three-state kinetics reduce to D A TET kinetics with an effective single forward single backward rate only part c. exact solution of the set A8 A16 for all Laplacetransformed populations F M (s) the solution is given in Appendix B. In most cases a numerical computation of the overall transfer rates becomes necessary. However further analytical considerations are possible whenever the intermediate states become less populated. This has been discussed at length for SET which is reduced to a single-exponential D A SET kinetics. 20 Therefore we will concentrate on the description of such a simple overall time behavior. A twoexponential TET process will be described in the following section whereas single-exponential behavior is discussed in the section afterwards. A. Two-exponential TET kinetics According to the general solution Eq. 30 of the rate equations A8 A16 two-exponential TET kinetics requires the existence of two rates which are much smaller than all other ones. Let these rates be K 1 K 2. Then the exact solution 30 reduces to P M tp M B M (1) e K 1 t B M (2) e K 2 t 31 which of course is only valid in a time domain t K 1 3 K K 22. lim sf B s s 0 m1 lim sf m s lim sf m s1. s 0 s 0 33 This relation has to be considered as a necessary condition for the appearance of two-exponential kinetics see also Ref. 18. Since all the Laplace-transformed populations F M (s) have a rather complicated form this is also valid for the inequality 33. To get analytical expressions for the two smallest transfer rates K 1 K 2 we shall start from the exact Eqs. B10 which however only determine the Laplace-transformed populations of the DBA states D I A the bridge populations are given by Eqs. B1 B2. If this reduced set of rate equations is transformed back into the time domain we get the following three kinetic equations see the scheme b of Fig. 4: Ṗ D tk f (1) k DA P D tk b (1) P I tk AD P A t Ṗ I tk (2) f k (1) b P I tk (1) f P D tk (2) b P A t 34 Ṗ A tk b (2) k AD P A tk f (2) P I tk DA P 2 t. The rates introduced here have to be identified with Laplacetransformed expressions as introduced in the Appendix B Eq. B10 k f (1) k DA R DD 0 k f (2) k b (1) R II 0 k b (2) k AD R AA 0 k f (1) R ID 0 35 k b (1) R DI 0 k f (2) R AI 0 k b (2) R IA 0. Using the Eqs. B3 B4 B6 B11 as well as the recursion formula B7 we may identify the rates of the Eqs. 34 with basic rate expressions of sequential superexchange transitions k (1) f k (seq) DI k (sup) DI k (seq) DI k D1 k I 1 /D 1 k (seq) ID k 1D k I 1 /D 1 k (2) f k (seq) IA k (sup) IA k (seq) IA k I1 k A 1 /D 2 k (seq) AI k 1I k A 1 /D 2 k (1) b k (seq) ID k (sup) ID k (sup) DI k DI k (sup) ID k ID k (2) b k (seq) AI k (sup) AI k (sup) IA k IA k (sup) AI k AI 36 37

7 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4447 where D 1 k I 1 k 1D 1 k 1D k I D2 38 D 2 k A 1 k 1I 1 k 1I k A D2. The main bridge length dependence of rates of transfer is originated by the function M 1 M1 DM Ds0;M 39 1 where /exp/k B T1 40 denotes the ratio of the backward the forward site-to-site rate constants. In Eq. 40 is the intersite energy bias introduced in Eq. 11. The transfer rates define two overall TET rates K a 1 d 1 a 1 d 1 2 4a 2 d 2 41 K a 1 d 1 a 1 d 1 2 4a 2 d 2 where we have introduced the following abbreviations: a 1 k f (2) k b (2) k AD a 2 k DA k f (2) d 1 k f (1) k b (1) k DA d 2 k AD k b (1). 42 Having determined the two overall TET rates we may present the other quantities entering the Eqs. 31 for the time-dependent populations. They read P D a 1k (1) (2) b d 2 k f K 1 K 2 P I K 1K 2 a 1 a 2 k b (1) d 1 d 2 k f (2) K 1 K 2 43 P A d 1k f (2) a 2 k b (1) K 1 K 2 as well as ( j12) ( B j) D 1 a 1K j1 j k (1) (2) b K j d 2 k f K j K 1 K 2 B I ( j) 1 j a 1K j k b (1) K j K j K 1 K 2 a 2k b (1) K j k f (2) d 1 d 2 K j K j K 1 K 2 44 ( B j) A 1 k (2) j1 f d 1 K j a 2 k (1) b K j. K j K 1 K 2 At this point of our discussion let us formulate the conditions for which the exact solution Eqs. B10 becomes identical with that of the kinetic equations 34. Therefore we have to consider the quantities R M (s) which have been introduced when considering the Laplace-transformed rate equations cf. Appendix B. We exp the R M (s) with respect to s assume that the approximations sr MM (s)sr MM (0) sr MM (0)sR MM (0) as well as R M (s)r M (0) sr M (0)R M (0) are fulfilled for both roots sk 1 sk 2 MDIA; R M (0)(dR M (s)/ds) s0 ]. Then one can see that the multistate DBA system shows a two-exponential TET kinetics Eq. 31 only if R M 0 M K 1 j R M 0 j A comparison of this inequality with the condition Eq. 33 for a small steady bridge population indicates their complete correspondence. B. Single-exponential D A TET kinetics So far we arrived at a description of bridge-assisted TET where the two forward (k (1) f k (2) f ) two backward (k (1) b k (2) b ) transfer rates contain additive contributions from the single-electron sequential the single-electron superexchange mechanisms. Remember that this additivity is a direct consequence of the small population of all bridging states. This two-exponential TET corresponds to an effective three-state system. If the intermediate state I carries a pronounced population the overall transfer rates K 1 K 2 note that K 1 K 2 ) contain a mixture of rate expressions caused by the single-electron stepwise the two-electron superexchange concerted mechanism. The situation changes substantially when K 1 K 2 46 which indicates that the TET process Eq. 31 is determined by two strongly separated kinetic phases a fast a slow one with the characteristic times 1 K K 1 2 respectively. If t 2 we have exp(k 2 t)1 thus the fast phase is defined by the time-dependent factor exp(k 1 t) for all three populations P D (t) P I (t) P A (t). The first phase of the TET is finished for t 1 but t 2. As a result a redistribution of the initial populations P M (0) MD took place for 2 t 1 the transient populations P M (t 1 )P M (0) have been formed. They are easily derived from Eq. 31 at exp(k 1 t)0 exp(k 2 t)1. Below we present concrete expressions for the overall transfer rates K 1 K 2 as well as for the intermediate populations P M (0) bearing in mind that the energy E I exceeds the energies E D E A cf. Fig. 3. For such a condition the following inequality is satisfied: k b (1) k f (2) k f (1) k b (2) k DA k AD. 47 ote that in line with Eq. 41 the inequality 46 is reduced to a 1 d 1 2 4a 1 d 1 a 2 d Therefore the following approximate expression for the overall transfer rates can be presented K 1 a 1 d 1 K 2 a 1d 1 a 2 d a 1 d 1 According to the inequalities 47 this can be rewritten as K 1 k b (1) k f (2) K 2 k f (1) k f (2) k b (2) k b (1) k b (1) k f (2) k DA k AD 50

8 4448 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May the transient populations take the following form: P D 01k f (1) k DA /k b (1) k f (2) P I 0k f (1) /k b (1) k f (2) P A 0k b (2) /k b (1) k f (2). 51 From these expressions it is seen that the ratio P I (0)/P A (0) is determined by the ratio of the stepwise the concerted rate constants k (1) f /k DA. ote that P I (0) gives the maximal value possible for the population of the intermediate state I cf. Fig. 4. But inequality 47 tells us also that P I (0)1 i.e. the population of state I remains very small at any time t P I t1. 52 This inequality indicates that at tk 1 1 the two-exponential kinetics 31 is reduced to a single-exponential one P D tp D 1P D e K TET t P A tp A 1e K TET t. 53 This kinetics describes a D A TET with an overall transfer rate K 2 K TET. The latter can be represented as the sum of a stepwise a concerted contribution K TET 1 TET K (step) TET K (conc) TET. 54 Since K TET k f k b contains a forward a backward part cf. scheme c in Fig. 4 K (step) TET K (conc) TET separate correspondingly. We have K (step) ( f TET K ) step K (b) step ( f K ) step k (1) (2) f k f k (1) (2) b k f K (conc) TET K conc ( f ) (b) K conc ( f K ) conc k DA K (b) conc k AD. IV. DISCUSSIO OF RESULTS K (b) step k b (2) (1) k b k (1) (2) b k f We consider the derivation of the overall TET rates the description of the conditions at which a multiexponential TET is reduced to a two-exponential as well as a singleexponential TET as the main result of the present studies. Additionally it could be shown in which way the overall transfer rates are expressed by the sequential the superexchange rate constants. When considering distant electron transfer it is of interest to describe the dependence of rate expressions on the number of bridge units. Concentrating on the singleexponential D A TET this dependence is contained in the stepwise concerted components of the overall transfer rate 54. First let us consider the stepwise channel of the D A TET. The stepwise transfer rates contain contributions from the single-electron sequential singleelectron superexchange transfer rates between the D the I as well as between the I the A electronic states of a whole DBA system. Using the definitions one can rewrite the sequential contribution as (seq) k DI(IA) k (1(2)seq) f 0 1 1(2) R 57 (seq) k ID(AI) k b0 (1(2)seq) 11R. 1 1(2) R It follows from Eq. 57 that the dependence of the sequential mechanism of SET is concentrated in the factor R1 1 /1 58 provided that 0 that the sequential decay parameters 1 2 are not too small. The quantities (1 k seq) f 0 k D1k I k I k 1D (2 k seq) f 0 k I1k A k A k 1I (1 k seq) b0 k Ik 1D k I k 1D (2 k seq) b0 k Ak 1I k A k 1I define rates of sequential transfer through a bridge with a single unit. The related decay parameters read 1 k 1Dk I 1 k I k 1D 61 2 k 1Ik A k A k 1I ote that due to a single-electron character of sequential pathway the rate constants in Eqs are identical to the hopping rate constants between the nearest sites of electron localization. Therefore e.g. the rate k D1 (k I1 ) characterizes the hopping of an electron from the state with two electrons single electron at the D to the adjacent wire unit while the rate k A (k I ) characterizes the similar process with respect to the A. To derive the single-electron superexchange contribution to the stepwise transfer rates one has to specify the superexchange rate constants k DI(ID) k IA(AI). These rate constants are defined in Appendix A Eq. A7 are proportional to the square of the superexchange couplings T DI 2 T IA 2 respectively. Here we restrict ourself to the consideration of a small intersite bridge bias where the effective transfer couplings are given by Eqs. C21 C23 Appendix C. A small change in the notation leads to the following expressions: (sup) (1(2) k DI(ID) k sup) f 0 e 1 (1) 63 (sup) (1(2) k IA(AI) k sup) b0 e 2 (1). 64 Here the rate k (sup) DI (k (sup) AI ) characterizes the distant coherent transfer of an electron from the state with two electrons at the D(A) to the state with one electron at the A(D). The rate (sup) ) describes the reverse process. Compare the k (sup) ID (k IA

9 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4449 definition of DBA electronic states 2 4 the scheme a of Fig. 4. The superexchange decay parameters V 1 2ln B 1/2 65 E D E I V 2 2ln B 1/2 66 Ẽ A Ẽ I are mainly responsible for the distant dependence of the superexchange rate. ote that in the presence of an energetic bias in the bridge an additional dependence also appears in the quantities (1(2) k sup) f 0 2 (1(2) k sup) b0 2 V D1 V A 2 E D E I FC DI(ID) V D1 V A 2 67 Ẽ A Ẽ I FC IA(AI). 68 The dependence is located in the Franck Condon factors through the driving forces E ID E IA cf. Eqs. C10 C11 characterizing the corresponding single-electron reaction. Two couples of decay parameters ( 1 1 ) ( 2 2 ) characterize the efficiency of the first the second singleelectron bridge-mediated stepwise pathway respectively. At the same time the efficiency of the concerted pathway is defined by the two-electron superexchange rate constants k DA k AD compare Eq. 56. Introducing the twoelectron superexchange decay parameter as V 2ln B 2 1/2 69 E D E A Ẽ D Ẽ A utilizing Eq. C25 a respective simplified notation of each concerted rate can be given (0) k DA(AD) k DA(AD) e (1). 70 Here the rate k DA (k AD ) characterizes the distant coherent transfer of two electrons from the state with two electrons at the D(A) to the state without extra electrons at the A(D). Such a process is possible owing to a virtual electronic configuration D BA corresponding to the state with one electron at the D center one electron at the A center. In the case of a bridge with an energy bias the factor (0) k DA(AD) 2 V D1 V A V D1 V A 2 E D E A Ẽ D Ẽ A E ID E IA FC DA(AD) 71 includes an additional dependence originated by the energy gaps E ID E IA as well as by the Franck Condon factor through the driving force E DA of the TET reaction cf. Eq. C12. In order to demonstrate the bridge-length dependence of the rate we take Jortner s form 39 for all rate constants. Let M be the single characteristic frequency accompanying the M electronic transition M be the corresponding reorganization energy. Then Jortner s form of the rate constant A7 can be written as see also Refs k M 2 V M 2 M 72 M where we have introduced M exp S M coth M k B T 1n M n M M /2 I M 2S M n M 1n M 1/2. 73 The expression contains the modified Bessel function I (z) the Bose distribution function n()exp(/k B T)1 1 we have set S M M / M M E M / M. First let us discuss the reduction of the multiexponential TET kinetics to the two-exponential one. In line with the theory proposed such a reduction is only possible at a small integral bridge population which is completely justified by the results of Fig. 5. Moreover the derived analytic results on the two-exponential kinetics cf. Eqs coincide completely cf. Fig. 5a with those obtained by a numerical solution of the basic set of Eqs. A8 A16 provided that 32 is fulfilled during the TET process. ote that the bridge population evolve in time as well cf. the numerical results depicted in Fig. 5b. But due to the small population of each bridging state this kinetics is well separated from the kinetics of the populations P D (t) P I (t) P A (t). If in the course of the TET process the bridging states B m Eq. 5 B n Eq. 6 as well as the intermediate state I Eq. 3 have a small population then the TET is reduced to a single-exponential D A TET. This statement is supported by the results shown in Fig. 6. Again despite the fact that both types of bridge populations P m (t) P ñ (t) denoted in Fig. 6b by the numerals 1 2 respectively as well as the intermediate state population P I (t) show a time dependence according to their respective transfer rates the basic populations P D (t) P A (t) display a singleexponential kinetics 53 with overall D A TET rate 54. Figures 7 10 show the D A TET rate in dependency on the number of bridging units. Special attention is put on the temperature dependence since it determines the contribution by the thermally activated stepwise transfer mechanism. Figure 7a demonstrates that the stepwise contribution to the overall transfer rate K TET Eq. 54 exceeds the concerted one so that K TET K (step) TET. Indeed the analysis of the interstate transfer rates indicates that at a given set of parameters specifying the elementary rate constants 72 the inequality k (2) f k (1) b cf. Fig. 7b is fulfilled thus the stepwise transfer rate Eq. 55 follows as K (step) TET k (1) b expe ID /k B T. 74 This expression indicates that the limiting step of the stepwise TET process is related to a single-electron D I transition. ote that a noticeable contribution in this thermally activated transition is given by the superexchange at

10 4450 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May FIG. 5. TET kinetics in a bridge with five units without energy bias for a small population P m (t) P m (t) of the two types of bridge states. Analytic results for the populations P D (t) P I (t) P A (t) according to Eqs the numerical solution of the set of Eqs. A8 A16 completely agree part a. The bridge populations related to the first 1 the second type of bridge states 2 show their own kinetics panel b. The calculations are based on the Eqs by choosing the following parameters: E D 0.29 ev E I 0.27 ev Ẽ I 0.21 ev Ẽ A 0.26 ev E0.03 ev; 1D I DI 0.6 ev I1 A IA 0.4 ev DA 0.8 ev mm1 B 0.8 ev; M cm ; V D1 V A V D1 V A 0.03 ev V B 0.04 ev. 1 4) the sequential at 4) single-electron pathway compare the dependence of both contributions plotted in Fig. 7b. At lower temperatures the contribution of the concerted mechanism increases compare Fig. 7a with Fig. 8a Fig. 8b. For instance at T150 K Fig. 8a the concerted mechanism dominates the TET at 1 2 while at T100 K the same mechanism works effectively up to 4 Fig. 8b. In contrast to the stepwise mechanism the concerted mechanism works also at small temperatures down to T0. ote that the stepwise pathway also includes a superexchange single-electron contribution. But due to the thermally activated origin of the stepwise rate Eq. 74 the superexchange single-electron contribution can only work at the presence of thermal activation. This fact reflects the specificity of the D A TET process. However in the case of single-electron D A transfer the superexchange channel contributes down to T0.) Therefore in the case of D A TET the contribution of the concerted mechanism increases at decreasing temperature. It follows from the comparison of Fig. 7b with Fig. 9 that the contribution of the single-electron superexchange FIG. 6. D A kinetics of a TET reaction in a bridge with five units without energy bias for a small bridge state population a small population of the intermediate state. Analytic results cf. Eqs the numerical solution of the set of Eqs. A8 A16 completely agree part a. The bridge populations related to the first 1 the second type of bridge states 2 as well as the population of the intermediate state I show their own kinetics panel b. The calculations are based on the E D 0.29 ev E I 0.09 ev Ẽ I 0.03 ev Ẽ A 0.3 ev E0.07 ev. All other parameters are identical to those used for Fig. 5. channel to the interstate transfer rates characterizing the stepwise pathway strongly increases at decreasing temperature. This means that at low temperatures even the stepwise pathway is used since it covers single-electron superexchange channels which however have to be thermally activated in contrast with the two-electron superexchange channel. Figure 10 displays the influence of the energy bias on the concerted the stepwise bridge-assisted superexchange contribution to the overall D A TET. At low temperatures under consideration just the single-electron the two-electron superexchange channels determine the thermally activated stepwise the concerted pathway of the electron transfer respectively. It is seen that at 1 4 the concerted contribution exceeds the stepwise one while at 4 the distance dependence of the K TET is determined by the stepwise mechanism of the D A TET. This means that at low temperatures the main contribution to the D A TET stems from the single-electron the two-electron superexchange mechanisms stepwise concerted pathways respectively. Figure 10 also shows that a change of generally affects the stepwise transfer rate. This follows from the fact that at 4 for which K TET is given by Eq. 74 the

11 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4451 FIG. 7. Bridge length dependence of the overall D A TET transfer rate K TET at room temperature as well as its stepwise (K (step) TET ) concerted (K (conc) TET ) components part a. The different contributions to the stepwise (step) transfer rate K TET single-electron transfer channels superexchange channel up to 3 4 bridge units the sequential transfer channel at 4) are shown in part b. The calculations are based on the Eqs by choosing the following parameters: E D 0.20 ev E I 0.10 ev Ẽ I 0.05 ev E0.15 ev; 1D I 1I 1.2 ev DI 0.8 ev A 1.4 ev IA 0.5 ev DA 0.2 ev mm1 B 0.5 ev; M cm ; V D1 V A V D1 V A 0.02 ev V B 0.07 ev. FIG. 8. Bridge length dependence of the overall D A TET transfer rate K TET as well as its stepwise (K (step) TET ) concerted (K (conc) TET ) components at different temperatures other parameters like those of Fig. 7. Part a T 150 K part b T100 K. kinetics of the TET process Eq. 30 is reduced to a twoexponential kinetics Eq. 31. The respective two overall transfer rates K 1 K 2 are given in Eq. 41. It has been demonstrated that such a reduction becomes possible if the main dependence of K TET on is concentrated at E ID Eq. C10 rather than in the rate constant k b (1) Eq. 36. This circumstance indicates directly that the stepwise D A TET process is originated by a thermally activated mechanism. V. COCLUSIO The present paper has been devoted to the study of nonadiabatic two-electron transfer TET in a donator (D) acceptor (A) complex mediated by a molecular bridge. Focusing on nonadiabatic reactions the intersite electron transitions take place on a time scale t rel where rel characterizes the time of intrasite vibrational relaxation. Therefore the whole analysis can be based on the coarse-grained kinetic equations A8 A16 valid for t rel. An analysis of the exact solution based on the Laplacetransformed state populations allowed to formulate the conditions Eqs for which the multiexponential FIG. 9. Bridge length dependence of the fastest interstate transfer rates k b (1) k f (2) characterizing the transitions from the intermediate state I in the three-state DBA system see scheme b in Fig. 4. The calculations have been done with the same parameters as in Fig. 7.

12 4452 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 E. G. Petrov V. May FIG. 10. Low temperature bridge length dependence of the overall D A TET rate K TET at different values of the bridge internal energy bias. The calculations have been done with the same parameters as in Fig. 7. integral population of the bridging states remains small less then 10 2 ) during the whole TET process. If additionally the population of the intermediate state I Eq. 3 with one electron at the D one at the A) remains also small then the two-exponential kinetics turns to a single-exponential one cf. Eq. 53 characterized by a single overall D A TET transfer rate K TET. Certain conditions have been derived cf. Eqs at which the transfer rate K TET can be given by a sum of the stepwise the concerted transfer rates Eq. 54. Within the model of a regular bridge we have presented analytic expressions for both mentioned transfer rates Eqs as well as Eqs The stepwise transfer rate characterizes the single-electron pathways of D A TET. It contains two contributions related to the sequential superexchange mechanism of bridge-assisted singleelectron transfer between the state D both electrons at the D) the state I one electron at the D another electron at the A) as well as single-electron transfer between the state I the state A both electrons at the A). Since the single-electron pathways follow from thermal activation the stepwise contributions to K TET are strongly suppressed with decreasing temperature. In contrast the concerted transfer rate which is originated by a specific twoelectron superexchange coupling between the state with two electrons at the D the state with two electrons at the A) can act at low temperatures down to T0. Therefore the contribution of the concerted mechanism to the D A TET increases with decreasing temperature. The model of a regular bridge has been used to derive analytical results on the dependence of the stepwise as well as concerted contributions to the overall D A TET rate. The Coulomb interaction between the transferred electrons located within the bridge as well as at the D A induces a shift of the site energies. This may also happen for the bridge-site reorganization energies mm1 which can differ from B. In the case of a single-electron D A transfer this problem has been already discussed in e.g. Ref. 40. If however among all possible bridging states only those with a single electron in the bridge are involved in the TET process the general form of the nonadiabatic rate constants 72 see also the more general expression A7 does not change. In particular the square of the superexchange couplings are further defined by Eqs. C1 C3. Only the energy gaps Eqs. C4 C12 as well as the reorganization energies have to be modified. In the present paper however we assumed a strong screening of the D A redox centers by surrounding polar groups as it is often the case in enzymes. Therefore any Coulomb interaction of the transferred electrons with these centers could be neglected. Studying TET in a bridge with an energy bias a deviation from the stard exponential law exp(1) characterizing the decrease of the superexchange rate with the number of bridge units has been obtained. This deviation could be explained by the change of the energy gap for every bridging unit. The key result of the present paper is the specification of the stepwise concerted contributions to the overall transfer rate which characterizes the nonadiabatic D A TET process mediated by a molecular bridge. This has been done in the framework of a tight binding model where the electronic energies intersite couplings enter as phenomenological parameters. Such a model allowed us to derive the analytic dependence of both contributions on the bridge length number of bridging units ). ote that the examination of the bridge-length dependence of the overall transfer rate is the most direct way to specify the stepwise concerted electron transfer routes in a concrete system. However even for a SET process such an experimental examination has been done only recently see discussion in Refs Therefore a detailed test of bridge-length effects on nonadiabatic TET reactions by comparing with experimental data has to be postponed to the future. Another way to distinguish the stepwise concerted routes is an examination of the ph dependence of the overall transfer rate characterizing the TET reactions in enzymes. Recently we used this way to analyze the reduction of micothione reductase by ADPH. 41 It has been shown that such a reduction takes place along the concerted route. Probably further progress in the theoretical description of bridge-assisted TET could be achieved with the combination of semiphenomenological analytic expressions containing the precise bridge-length dependence of transfer rates quantum-chemical estimations of electronic energies of the various site-to-site superexchange electronic couplings. In the case of SET just such methods allowed to specify single-electron pathways between redox-centers of number protein macromolecules ACKOWLEDGMETS This work is supported in part by the Fund of Fundamental Research of ASU the Volkswagen-Stiftung Germany priority area Intra- Intermolecular Electron Transfer. APPEDIX A: DERIVATIO OF KIETIC EQUATIOS The concrete form of the kinetic equations follows from Eq. 28 if the expansion 22 is introduced into Eq. 29.

13 J. Chem. Phys. Vol. 120 o. 9 1 March 2004 Two-electron transfer reactions 4453 The first term R 2 of this expansion corresponds to the stard second-order Born approximation. It leads to rate constants which describe the sequential pathway of TET. Higher-order terms result in superexchange single-electron two-electron rate constants. A similar result has been presented earlier in Ref. 18 for the derivation of SET rate equations. Therefore we do not repeat details of the derivation. We only present the superexchange couplings for the transitions between the two-electron states D I I A as well as D A. Moreover we will only consider pure electronic energy gaps E M in Fig. 3 these gaps are shown for a bridge without energy bias. Let M be the donor state D. Then the Born approximation of the quantities K D Eq. 29 follows from the general expression 29 as K D K (Born) D D k D1 1 k 1D. A1 Here k D1 k 1D exp(e 1 E D )/k B T k 1D 2 D 1 V 11 D D 2 WE1 1 LE1 1 ED D. A2 The quantity L(x) denotes the Lorentzian-type broadening of the function. Employing the Condon approximation Eq. 12 following the approach given in Ref. 18 we can write the rate constant k 1D in the stard form: k 1D 2 V 1D 2 FC 1D. A3 Here (FC) 1D is the Franck Condon factor for the 1 D electron transition. In the considered case of nonadiabatic electron transfer any rate constant characterizing electron hopping between adjacent sites has a similar form. Higher-order expansion terms in Eq. 22 result in two effects the renormalization of those rate constants derived in the Born approximation the appearance of superexchange rate constants. In the case of nonadiabatic electron transfer the rate renormalization is of less importance will be neglected. The superexchange contributions if related to the D are of the single-electron type rate constants k DI k ID ) of the two-electron concerted type rate constants k DA k AD ). Therefore we get K D K (Born) (sup D K 1) (sup D K 2) D A4 where the first term is given by Eq. A1 while (sup K 1) D D k DI I k ID A5 (sup K 2) D D k DA I k AD. A6 In Eqs. A5 A6 the superexchange rate constants take the form k M 2 T M 2 FC M A7 where expressions of the square of single-electron twoelectron superexchange couplings defining the D A TET are given in Appendix C. Substituting Eqs. A4 A6 in Eq. 28 we get Ṗ D tq D k DA P D tk 1D P 1 tk ID P I t k AD P A t Ṗ 1 tq 1 P 1 tk D1 P D tp 2 t Ṗ m tqp m tp m1 tp m1 t Ṗ tq P tk I P I tp 1 t Ṗ I tq I (1) q I (2) P I tk I P tk DI P D t k 1I P 1tk AI P A t Ṗ 1 tq 1P 1 tk I1 P I tp 2 t Ṗ ñ tqp ñ tp ñ1 tp ñ1 t Ṗ Ñ tq P Ñ tk A P A tp Ñ1 t Ṗ A tq A k AD P A tk A P Ñ tk IA P I t A8 A9 A10 A11 A12 A13 A14 A15 k DA P D t. A16 ote that in the Eqs. A10 A14 mn run over Furthermore we used the abbreviation m B m ñb n so that we can identify P m (t) with P Bm (t) P ñ (t) with P B n(t). The Eqs. A8 A16 have to be complemented by the following abbreviations for the different rate expressions: q D k D1 k DI q 1 k 1D q k I q 1k 1I q k A q A k A k AI A17 q (1) I k I k ID q (2) I k I1 k IA q. Here k mm1 k m1m (m12...1) are the rates characterizing the forward backward jumps of a single-electron between neighboring bridge units. APPEDIX B: EXACT SOLUTIO OF THE LAPLACE TRASFORMED SET OF RATE EQUATIOS The solution of the Eqs. A8 A16 can be constructed by changing to the Laplace transform F M (s) 0 exp(st)p M (t) dt of all electronic populations. oting the initial condition P M (0) MD the solution of the Laplace-transformed rate equations A9 A11 A13 A15 governing the bridge populations P m P ñ are obtained as F m s 1 D 1 s k D1 m1 C mf D s as k I m C 1 mf I s F ñ s 1 D 2 s k I1 n1 C nf I s k A n C 1nF A s. The following abbreviations have been introduced: B1 B2

12.2 MARCUS THEORY 1 (12.22)

12.2 MARCUS THEORY 1 (12.22) Andrei Tokmakoff, MIT Department of Chemistry, 3/5/8 1-6 1. MARCUS THEORY 1 The displaced harmonic oscillator (DHO) formalism and the Energy Gap Hamiltonian have been used extensively in describing charge

More information

Quantum Master Equations for the Electron Transfer Problem

Quantum Master Equations for the Electron Transfer Problem 20/01/2010 Quantum Master Equations for the Electron Transfer Problem Seminarvortrag Dekohaerenz und Dissipation in Quantensystemen Antonio A. Gentile The general transport problem in micro/mesoscopic

More information

Charge and Energy Transfer Dynamits in Molecular Systems

Charge and Energy Transfer Dynamits in Molecular Systems Volkhard May, Oliver Kühn Charge and Energy Transfer Dynamits in Molecular Systems Second, Revised and Enlarged Edition WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA Contents 1 Introduction 19 2 Electronic

More information

Lecture 12. Electron Transport in Molecular Wires Possible Mechanisms

Lecture 12. Electron Transport in Molecular Wires Possible Mechanisms Lecture 12. Electron Transport in Molecular Wires Possible Mechanisms In Lecture 11, we have discussed energy diagrams of one-dimensional molecular wires. Here we will focus on electron transport mechanisms

More information

Bound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws

Bound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws PHYSICAL REVIEW B VOLUME 55, NUMBER 8 15 FEBRUARY 1997-II Bound states of two particles confined to parallel two-dimensional layers and interacting via dipole-dipole or dipole-charge laws V. I. Yudson

More information

5.74 Introductory Quantum Mechanics II

5.74 Introductory Quantum Mechanics II MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,

More information

Heyrovský Sq. 2, Prague 6, Czech Republic. Na Slovance 2, Prague 8, Czech Republic. 1. Introduction

Heyrovský Sq. 2, Prague 6, Czech Republic. Na Slovance 2, Prague 8, Czech Republic. 1. Introduction Materials Science-Poland, Vol. 27, No. 3, 2009 Nonradiative electron and energy transfer. Explicit estimation of the influence of coherent and dephasing processes in a vibrational bath on electronic dynamics

More information

11.1. FÖRSTER RESONANCE ENERGY TRANSFER

11.1. FÖRSTER RESONANCE ENERGY TRANSFER 11-1 11.1. FÖRSTER RESONANCE ENERGY TRANSFER Förster resonance energy transfer (FRET) refers to the nonradiative transfer of an electronic excitation from a donor molecule to an acceptor molecule: D *

More information

(Excerpt from S. Ji, Molecular Theory of the Living Cell: Concepts, Molecular Mechanisms, and Biomedical Applications, Springer, New York, 2012)

(Excerpt from S. Ji, Molecular Theory of the Living Cell: Concepts, Molecular Mechanisms, and Biomedical Applications, Springer, New York, 2012) 2.2 The Franck-Condon Principle (FCP) 2.2.1 FCP and Born-Oppenheimer Approximation The Franck-Condon Principle originated in molecular spectroscopy in 1925 when J. Franck proposed (and later Condon provided

More information

Theoretical Models for Chemical Kinetics

Theoretical Models for Chemical Kinetics Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular

More information

Electrons in a periodic potential

Electrons in a periodic potential Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x

More information

Laser Induced Control of Condensed Phase Electron Transfer

Laser Induced Control of Condensed Phase Electron Transfer Laser Induced Control of Condensed Phase Electron Transfer Rob D. Coalson, Dept. of Chemistry, Univ. of Pittsburgh Yuri Dakhnovskii, Dept. of Physics, Univ. of Wyoming Deborah G. Evans, Dept. of Chemistry,

More information

Marcus Theory for Electron Transfer a short introduction

Marcus Theory for Electron Transfer a short introduction Marcus Theory for Electron Transfer a short introduction Minoia Andrea MPIP - Journal Club -Mainz - January 29, 2008 1 Contents 1 Intro 1 2 History and Concepts 2 2.1 Frank-Condon principle applied to

More information

Introduction to Heisenberg model. Javier Junquera

Introduction to Heisenberg model. Javier Junquera Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics

More information

How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics

How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics Reaction rated are fractions of a second for fireworks to explode. Reaction Rates takes years for a metal

More information

arxiv:cond-mat/ v1 17 Aug 1994

arxiv:cond-mat/ v1 17 Aug 1994 Universality in the One-Dimensional Self-Organized Critical Forest-Fire Model Barbara Drossel, Siegfried Clar, and Franz Schwabl Institut für Theoretische Physik, arxiv:cond-mat/9408046v1 17 Aug 1994 Physik-Department

More information

10/26/2010. An Example of a Polar Reaction: Addition of H 2 O to Ethylene. to Ethylene

10/26/2010. An Example of a Polar Reaction: Addition of H 2 O to Ethylene. to Ethylene 6.5 An Example of a Polar Reaction: Addition of H 2 O to Ethylene Addition of water to ethylene Typical polar process Acid catalyzed addition reaction (Electophilic addition reaction) Polar Reaction All

More information

arxiv: v1 [cond-mat.stat-mech] 6 Mar 2008

arxiv: v1 [cond-mat.stat-mech] 6 Mar 2008 CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey

More information

Analysis of the ultrafast dynamics of the silver trimer upon photodetachment

Analysis of the ultrafast dynamics of the silver trimer upon photodetachment J. Phys. B: At. Mol. Opt. Phys. 29 (1996) L545 L549. Printed in the UK LETTER TO THE EDITOR Analysis of the ultrafast dynamics of the silver trimer upon photodetachment H O Jeschke, M E Garcia and K H

More information

Supplementary Figure S1 Definition of the wave vector components: Parallel and perpendicular wave vector of the exciton and of the emitted photons.

Supplementary Figure S1 Definition of the wave vector components: Parallel and perpendicular wave vector of the exciton and of the emitted photons. Supplementary Figure S1 Definition of the wave vector components: Parallel and perpendicular wave vector of the exciton and of the emitted photons. Supplementary Figure S2 The calculated temperature dependence

More information

A theoretical study of the single-molecule transistor

A theoretical study of the single-molecule transistor A theoretical study of the single-molecule transistor B. C. Friesen Department of Physics, Oklahoma Baptist University, Shawnee, OK 74804 J. K. Ingersent Department of Physics, University of Florida, Gainesville,

More information

1 Adsorption of NO 2 on Pd(100) Juan M. Lorenzi, Sebastian Matera, and Karsten Reuter,

1 Adsorption of NO 2 on Pd(100) Juan M. Lorenzi, Sebastian Matera, and Karsten Reuter, Supporting information: Synergistic inhibition of oxide formation in oxidation catalysis: A first-principles kinetic Monte Carlo study of NO+CO oxidation at Pd(100) Juan M. Lorenzi, Sebastian Matera, and

More information

Density functional calculation of nuclear magnetic resonance chemical shifts

Density functional calculation of nuclear magnetic resonance chemical shifts Density functional calculation of nuclear magnetic resonance chemical shifts Christoph van Wüllen Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany Received 22 August

More information

Identical Particles. Bosons and Fermions

Identical Particles. Bosons and Fermions Identical Particles Bosons and Fermions In Quantum Mechanics there is no difference between particles and fields. The objects which we refer to as fields in classical physics (electromagnetic field, field

More information

Dynamical principles in biological processes

Dynamical principles in biological processes Proc. Natl. Acad. Sci. USA Vol. 95, pp. 358 362, February 998 Chemistry Dynamical principles in biological processes E. W. SCHLAG, S.H.LIN, R.WEINKAUF, AND P. M. RENTZEPIS Institute for Physical Theoretical

More information

Decoherence and the Classical Limit

Decoherence and the Classical Limit Chapter 26 Decoherence and the Classical Limit 26.1 Introduction Classical mechanics deals with objects which have a precise location and move in a deterministic way as a function of time. By contrast,

More information

Competition between face-centered cubic and icosahedral cluster structures

Competition between face-centered cubic and icosahedral cluster structures Competition between face-centered cubic and icosahedral cluster structures R. S. Berry University of Chicago, Chicago, IL 60637, USA B. M. Smyrnov, and A. Yu. Strizhev High-Temperatures Institute, Russian

More information

The energy speed accuracy trade-off in sensory adaptation

The energy speed accuracy trade-off in sensory adaptation SUPPLEMENTARY INFORMATION DOI: 1.138/NPHYS2276 The energy speed accuracy trade-off in sensory adaptation 1 The effective potential H(m) In the main text, we introduced a set of F a and F m functions to

More information

Phonons and lattice dynamics

Phonons and lattice dynamics Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed

More information

Chemical Kinetics. Goal. Objectives

Chemical Kinetics. Goal. Objectives 4 Chemical Kinetics To understand the physical Goal factors that determine reaction rates. bjectives After this chapter, you should be able to: describe the factors that determine reaction rates. identify

More information

Pauli Deformation APPENDIX Y

Pauli Deformation APPENDIX Y APPENDIX Y Two molecules, when isolated say at infinite distance, are independent and the wave function of the total system might be taken as a product of the wave functions for the individual molecules.

More information

Theoretical Concepts of Spin-Orbit Splitting

Theoretical Concepts of Spin-Orbit Splitting Chapter 9 Theoretical Concepts of Spin-Orbit Splitting 9.1 Free-electron model In order to understand the basic origin of spin-orbit coupling at the surface of a crystal, it is a natural starting point

More information

Rate of Heating and Cooling

Rate of Heating and Cooling Rate of Heating and Cooling 35 T [ o C] Example: Heating and cooling of Water E 30 Cooling S 25 Heating exponential decay 20 0 100 200 300 400 t [sec] Newton s Law of Cooling T S > T E : System S cools

More information

The reaction whose rate constant we are to find is the forward reaction in the following equilibrium. NH + 4 (aq) + OH (aq) K b.

The reaction whose rate constant we are to find is the forward reaction in the following equilibrium. NH + 4 (aq) + OH (aq) K b. THE RATES OF CHEMICAL REACTIONS 425 E22.3a The reaction for which pk a is 9.25 is NH + 4 aq + H 2Ol NH 3 aq + H 3 O + aq. The reaction whose rate constant we are to find is the forward reaction in the

More information

Quantum Transport through Coulomb-Blockade Systems

Quantum Transport through Coulomb-Blockade Systems Quantum Transport through Coulomb-Blockade Systems Björn Kubala Institut für Theoretische Physik III Ruhr-Universität Bochum COQUSY6 p.1 Overview Motivation Single-electron box/transistor Coupled single-electron

More information

Matthias Lütgens, Frank Friedriszik, and Stefan Lochbrunner* 1 Concentration dependent CARS and Raman spectra of acetic acid in carbon tetrachloride

Matthias Lütgens, Frank Friedriszik, and Stefan Lochbrunner* 1 Concentration dependent CARS and Raman spectra of acetic acid in carbon tetrachloride Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 SUPPORTING INFORMATION Direct observation of the cyclic dimer in liquid acetic

More information

Part One: Reaction Rates. 1. Rates of chemical reactions. (how fast products are formed and/or reactants are used up)

Part One: Reaction Rates. 1. Rates of chemical reactions. (how fast products are formed and/or reactants are used up) A. Chemical Kinetics deals with: CHAPTER 13: RATES OF REACTION Part One: Reaction Rates 1. Rates of chemical reactions. (how fast products are formed and/or reactants are used up) 2. Mechanisms of chemical

More information

Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS

Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS Experimental Kinetics and Gas Phase Reactions Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS Professor Angelo R. Rossi http://homepages.uconn.edu/rossi Department of Chemistry, Room

More information

Quantum impurities in a bosonic bath

Quantum impurities in a bosonic bath Ralf Bulla Institut für Theoretische Physik Universität zu Köln 27.11.2008 contents introduction quantum impurity systems numerical renormalization group bosonic NRG spin-boson model bosonic single-impurity

More information

Diffusion influence on Michaelis Menten kinetics

Diffusion influence on Michaelis Menten kinetics JOURNAL OF CHEMICAL PHYSICS VOLUME 5, NUMBER 3 5 JULY 200 Diffusion influence on Michaelis Menten kinetics Hyojoon Kim, Mino Yang, Myung-Un Choi, and Kook Joe Shin a) School of Chemistry, Seoul National

More information

The quantum speed limit

The quantum speed limit The quantum speed limit Vittorio Giovannetti a,sethlloyd a,b, and Lorenzo Maccone a a Research Laboratory of Electronics b Department of Mechanical Engineering Massachusetts Institute of Technology 77

More information

F(t) equilibrium under H 0

F(t) equilibrium under H 0 Physics 17b: Statistical Mechanics Linear Response Theory Useful references are Callen and Greene [1], and Chandler [], chapter 16. Task To calculate the change in a measurement B t) due to the application

More information

ON FRACTIONAL RELAXATION

ON FRACTIONAL RELAXATION Fractals, Vol. 11, Supplementary Issue (February 2003) 251 257 c World Scientific Publishing Company ON FRACTIONAL RELAXATION R. HILFER ICA-1, Universität Stuttgart Pfaffenwaldring 27, 70569 Stuttgart,

More information

First-Principles Calculation of Exchange Interactions

First-Principles Calculation of Exchange Interactions Chapter 2 First-Principles Calculation of Exchange Interactions Before introducing the first-principles methods for the calculation of exchange interactions in magnetic systems we will briefly review two

More information

σ Bonded ligands: Transition Metal Alkyls and Hydrides

σ Bonded ligands: Transition Metal Alkyls and Hydrides σ Bonded ligands: Transition Metal Alkyls and Hydrides Simplest of organo-transitionmetal species Rare until and understanding of their stability in the 60 s and 70 s Metal alkyls can be considered as

More information

Electronic structure of correlated electron systems. Lecture 2

Electronic structure of correlated electron systems. Lecture 2 Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No

More information

arxiv: v1 [physics.chem-ph] 19 Dec 2018

arxiv: v1 [physics.chem-ph] 19 Dec 2018 arxiv:82.07969v [physics.chem-ph] 9 Dec 208 ransition state theory in systems of interacting particles Doriano Brogioli, Fabio La Mantia Universität Bremen, Energiespeicher und Energiewandlersysteme, Bibliothekstraße,

More information

NMR Dynamics and Relaxation

NMR Dynamics and Relaxation NMR Dynamics and Relaxation Günter Hempel MLU Halle, Institut für Physik, FG Festkörper-NMR 1 Introduction: Relaxation Two basic magnetic relaxation processes: Longitudinal relaxation: T 1 Relaxation Return

More information

Nuclear Effects in Electron Capture into Highly Charged Heavy Ions

Nuclear Effects in Electron Capture into Highly Charged Heavy Ions Nuclear Effects in Electron Capture into Highly Charged Heavy Ions W. Scheid 1,A.Pálffy 2,Z.Harman 2, C. Kozhuharov 3, and C. Brandau 3 1 Institut für Theoretische Physik der Justus-Liebig-Universität

More information

Chapter 2 Energy Transfer Review

Chapter 2 Energy Transfer Review Chapter 2 Energy Transfer Review In this chapter, we discuss the basic concepts of excitation energy transfer, making the distinction between radiative and nonradiative, and giving a brief overview on

More information

Bimolecular processes

Bimolecular processes Bimolecular processes Electron transfer *A + B A + + B - *A + B A - + B + EA IP *EA *IP LUMO An excited state is a better oxidant and a better reductant than the ground state HOMO X X* Kinetic of electron

More information

Transport of single molecules along the periodic parallel lattices with coupling

Transport of single molecules along the periodic parallel lattices with coupling THE JOURNAL OF CHEMICAL PHYSICS 124 204901 2006 Transport of single molecules along the periodic parallel lattices with coupling Evgeny B. Stukalin The James Franck Institute The University of Chicago

More information

Susana F. Huelga. Dephasing Assisted Transport: Quantum Networks and Biomolecules. University of Hertfordshire. Collaboration: Imperial College London

Susana F. Huelga. Dephasing Assisted Transport: Quantum Networks and Biomolecules. University of Hertfordshire. Collaboration: Imperial College London IQIS2008, Camerino (Italy), October 26th 2008 Dephasing Assisted Transport: Quantum Networks and Biomolecules Susana F. Huelga University of Hertfordshire Collaboration: Imperial College London Work supported

More information

arxiv:cond-mat/ v1 1 Dec 1999

arxiv:cond-mat/ v1 1 Dec 1999 Impurity relaxation mechanism for dynamic magnetization reversal in a single domain grain Vladimir L. Safonov and H. Neal Bertram Center for Magnetic Recording Research, University of California San arxiv:cond-mat/9912014v1

More information

Kinetic equation approach to the problem of rectification in asymmetric molecular structures

Kinetic equation approach to the problem of rectification in asymmetric molecular structures Kinetic equation approach to the problem of rectification in asymmetric molecular structures Kamil Walczak Institute of Physics, Adam Mickiewicz University Umultowska 85, 6-64 Poznań, Poland Transport

More information

Lecture 27. Transition States and Enzyme Catalysis

Lecture 27. Transition States and Enzyme Catalysis Lecture 27 Transition States and Enzyme Catalysis Reading for Today: Chapter 15 sections B and C Chapter 16 next two lectures 4/8/16 1 Pop Question 9 Binding data for your thesis protein (YTP), binding

More information

arxiv: v1 [cond-mat.dis-nn] 31 Aug 2011

arxiv: v1 [cond-mat.dis-nn] 31 Aug 2011 Suppression of the virtual Anderson transition in a narrow impurity band of doped quantum well structures. N.V. Agrinskaya, V.I. Kozub, and D.S. Poloskin Ioffe Physical-Technical Institute of the Russian

More information

The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is

The dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is 1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles

More information

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics

Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle

More information

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Q. Ma, NASA/GISS R. H. Tipping, Univ. of Alabama C. Boulet, Univ. Paris-Sud, France Statements The presentation is not a

More information

Coherent and incoherent motion in a one-dimensional lattice

Coherent and incoherent motion in a one-dimensional lattice PHYSICAL REVIEW B VOLUME 62, UMBER 6 Coherent and incoherent motion in a one-dimensional lattice AUGUST 2-II M. A. Palenberg and R. J. Silbey Department of Chemistry and Center for Materials Science and

More information

A determinantal formula for the GOE Tracy-Widom distribution

A determinantal formula for the GOE Tracy-Widom distribution A determinantal formula for the GOE Tracy-Widom distribution Patrik L. Ferrari and Herbert Spohn Technische Universität München Zentrum Mathematik and Physik Department e-mails: ferrari@ma.tum.de, spohn@ma.tum.de

More information

Triplet state diffusion in organometallic and organic semiconductors

Triplet state diffusion in organometallic and organic semiconductors Triplet state diffusion in organometallic and organic semiconductors Prof. Anna Köhler Experimental Physik II University of Bayreuth Germany From materials properties To device applications Organic semiconductors

More information

α particles, β particles, and γ rays. Measurements of the energy of the nuclear

α particles, β particles, and γ rays. Measurements of the energy of the nuclear .101 Applied Nuclear Physics (Fall 004) Lecture (1/1/04) Nuclear ecays References: W. E. Meyerhof, Elements of Nuclear Physics (McGraw-Hill, New York, 1967), Chap 4. A nucleus in an excited state is unstable

More information

Applications of Quantum Theory to Some Simple Systems

Applications of Quantum Theory to Some Simple Systems Applications of Quantum Theory to Some Simple Systems Arbitrariness in the value of total energy. We will use classical mechanics, and for simplicity of the discussion, consider a particle of mass m moving

More information

Theoretical formulation for electron transfer coupled to multiple protons: Application to amidinium carboxylate interfaces

Theoretical formulation for electron transfer coupled to multiple protons: Application to amidinium carboxylate interfaces JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 1 1 JULY 21 Theoretical formulation for electron transfer coupled to multiple protons: Application to amidinium carboxylate interfaces Ivan Rostov and Sharon

More information

Classical and quantum simulation of dissipative quantum many-body systems

Classical and quantum simulation of dissipative quantum many-body systems 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 0 20 32 Classical and quantum simulation of dissipative quantum many-body systems

More information

1. Introduction to Chemical Kinetics

1. Introduction to Chemical Kinetics 1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism

More information

Chem1B General Chemistry II Exam 1 Summer Read all questions carefully make sure that you answer the question that is being asked.

Chem1B General Chemistry II Exam 1 Summer Read all questions carefully make sure that you answer the question that is being asked. ChemB General Chemistry II Exam Summer 20 Name: KEY GSI: Write your name on all pages of the exam. Read all questions carefully make sure that you answer the question that is being asked. Write neatly

More information

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS

NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS NUMERICAL METHODS FOR QUANTUM IMPURITY MODELS http://www.staff.science.uu.nl/~mitch003/nrg.html March 2015 Anrew Mitchell Utrecht University Quantum impurity problems Part 1: Quantum impurity problems

More information

Multi-Scale Modeling from First Principles

Multi-Scale Modeling from First Principles m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations

More information

A population is modeled by the differential equation

A population is modeled by the differential equation Math 2, Winter 2016 Weekly Homework #8 Solutions 9.1.9. A population is modeled by the differential equation dt = 1.2 P 1 P ). 4200 a) For what values of P is the population increasing? P is increasing

More information

RPA in infinite systems

RPA in infinite systems RPA in infinite systems Translational invariance leads to conservation of the total momentum, in other words excited states with different total momentum don t mix So polarization propagator diagonal in

More information

Semiclassical molecular dynamics simulations of intramolecular proton transfer in photoexcited 2-2 -hydroxyphenyl oxazole

Semiclassical molecular dynamics simulations of intramolecular proton transfer in photoexcited 2-2 -hydroxyphenyl oxazole JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 21 1 DECEMBER 2000 Semiclassical molecular dynamics simulations of intramolecular proton transfer in photoexcited 2-2 -hydroxyphenyl oxazole Victor Guallar,

More information

Lecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions

Lecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion

More information

Reaction Thermodynamics

Reaction Thermodynamics Reaction Thermodynamics Thermodynamics reflects the degree to which a reaction is favored or disfavored Recall: G = Gibbs free energy = the energy available to do work ΔG = change in G of the system as

More information

Energetics of Electron Transfer Reactions

Energetics of Electron Transfer Reactions Energetics of Electron Transfer Reactions Dmitry Matyushov Arizona State University MIT, November 17, 24 Problems Tunneling between localized states E D A LUMO E HOMO < E> Instantaneous energy gap E becomes

More information

Overview of Kinetics

Overview of Kinetics Overview of Kinetics [P] t = ν = k[s] Velocity of reaction Conc. of reactant(s) Rate of reaction M/sec Rate constant sec -1, M -1 sec -1 1 st order reaction-rate depends on concentration of one reactant

More information

Quantum control of dissipative systems. 1 Density operators and mixed quantum states

Quantum control of dissipative systems. 1 Density operators and mixed quantum states Quantum control of dissipative systems S. G. Schirmer and A. I. Solomon Quantum Processes Group, The Open University Milton Keynes, MK7 6AA, United Kingdom S.G.Schirmer@open.ac.uk, A.I.Solomon@open.ac.uk

More information

CHEM Core Chemistry 3. Inorganic Reaction Mechanisms

CHEM Core Chemistry 3. Inorganic Reaction Mechanisms CHEM3012 - Core Chemistry 3 Inorganic Reaction Mechanisms 5. Mechanisms of electron transfer between metal ions This section of the course is concerned with the mechanisms of electron transfer reactions,

More information

Enduring Understandings & Essential Knowledge for AP Chemistry

Enduring Understandings & Essential Knowledge for AP Chemistry Enduring Understandings & Essential Knowledge for AP Chemistry Big Idea 1: The chemical elements are fundamental building materials of matter, and all matter can be understood in terms of arrangements

More information

3.320 Lecture 23 (5/3/05)

3.320 Lecture 23 (5/3/05) 3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials

More information

Diffusion in Dilute Alloys

Diffusion in Dilute Alloys Chapter 3 Diffusion in Dilute Alloys Our discussion of the atomistic mechanisms for diffusion has been confined to the situations where the diffusing species is chemically identical to the host atom. We

More information

Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture

Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture Non-Equilibrium Kinetics and Transport Processes in a Hypersonic Flow of CO 2 /CO/O 2 /C/O Mixture E.V. Kustova, E.A. Nagnibeda, Yu.D. Shevelev and N.G. Syzranova Department of Mathematics and Mechanics,

More information

Splitting of a Cooper pair by a pair of Majorana bound states

Splitting of a Cooper pair by a pair of Majorana bound states Chapter 7 Splitting of a Cooper pair by a pair of Majorana bound states 7.1 Introduction Majorana bound states are coherent superpositions of electron and hole excitations of zero energy, trapped in the

More information

Lecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics

Lecture 6. Fermion Pairing. WS2010/11: Introduction to Nuclear and Particle Physics Lecture 6 Fermion Pairing WS2010/11: Introduction to Nuclear and Particle Physics Experimental indications for Cooper-Pairing Solid state physics: Pairing of electrons near the Fermi surface with antiparallel

More information

Program for the rest of the course

Program for the rest of the course Program for the rest of the course 16.4 Enzyme kinetics 17.4 Metabolic Control Analysis 19.4. Exercise session 5 23.4. Metabolic Control Analysis, cont. 24.4 Recap 27.4 Exercise session 6 etabolic Modelling

More information

Supporting Information

Supporting Information Supporting Information Fluorination of Metal Phthalocyanines: Single-Crystal Growth, Efficient N-Channel Organic Field-Effect Transistors, and Structure- Property Relationships Hui Jiang 1*, Jun Ye 2,

More information

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić

Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824

More information

Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign

Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Understanding of correlated materials is mostly phenomenological FN- DMC

More information

ELECTRONS IN A PERIODIC POTENTIAL AND ENERGY BANDS IN SOLIDS-2

ELECTRONS IN A PERIODIC POTENTIAL AND ENERGY BANDS IN SOLIDS-2 ELECTRONS IN A PERIODIC POTENTIAL AND ENERGY BANDS IN SOLIDS-2 ENERGY BANDS IN A SOLID : A FORMAL APPROACH SCHROEDINGER'S EQUATION FOR A PERIODIC POTENTIAL * Electrons motion in a crystal will now be considered

More information

Berry s phase under the Dzyaloshinskii-Moriya interaction

Berry s phase under the Dzyaloshinskii-Moriya interaction PHYSICAL REVIEW A 77, 8 erry s phase under the Dzyaloshinskii-Moriya interaction M. K. Kwan, Zeynep ilhan Gurkan, and L. C. Kwek, ational Institute of Education, anyang Technological University, anyang

More information

5.74 Introductory Quantum Mechanics II

5.74 Introductory Quantum Mechanics II MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. p. 10-0 10..

More information

STM spectroscopy (STS)

STM spectroscopy (STS) STM spectroscopy (STS) di dv 4 e ( E ev, r) ( E ) M S F T F Basic concepts of STS. With the feedback circuit open the variation of the tunneling current due to the application of a small oscillating voltage

More information

Binding Theory Equations for Affinity and Kinetics Analysis

Binding Theory Equations for Affinity and Kinetics Analysis Technology Note #101 Binding Theory Equations for Affinity and Kinetics Analysis This technology note summarizes important equations underlying the theory of binding of solute analytes to surface-tethered

More information

Major Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1

Major Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1 Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response

More information

Classical Theory of Harmonic Crystals

Classical Theory of Harmonic Crystals Classical Theory of Harmonic Crystals HARMONIC APPROXIMATION The Hamiltonian of the crystal is expressed in terms of the kinetic energies of atoms and the potential energy. In calculating the potential

More information

SMALL BIPOLARONS IN BORON CARBIDES: PAIR BREAKING IN SEMICLASSICAL HOPPING* David Emin Sandia National Laboratories Albuquerque, NM , USA

SMALL BIPOLARONS IN BORON CARBIDES: PAIR BREAKING IN SEMICLASSICAL HOPPING* David Emin Sandia National Laboratories Albuquerque, NM , USA SMALL BIPOLARONS IN BORON CARBIDES: PAIR BREAKING IN SEMICLASSICAL HOPPING* David Emin Sandia National Laboratories Albuquerque, NM 87185-1421, USA 1. INTRODUCTION A pair of charge carriers can be bound

More information

Lecture 4: Hartree-Fock Theory

Lecture 4: Hartree-Fock Theory Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism

More information

S.K. Saikin May 22, Lecture 13

S.K. Saikin May 22, Lecture 13 S.K. Saikin May, 007 13 Decoherence I Lecture 13 A physical qubit is never isolated from its environment completely. As a trivial example, as in the case of a solid state qubit implementation, the physical

More information