Difference between revisions of "Example selectivity"
From MKMCXX
(→Example: Selectivity and Degree of Selectivity Control) |
|||
Line 1: | Line 1: | ||
== Introduction == | == Introduction == | ||
In this example, we show how to calculate selectivities and how to perform a degree of selectivity control analysis. | In this example, we show how to calculate selectivities and how to perform a degree of selectivity control analysis. | ||
+ | |||
+ | == Input == | ||
+ | In this example, the following input file is used: | ||
+ | |||
+ | <pre> | ||
+ | &compounds | ||
+ | A; 0; 1 | ||
+ | E; 0; 0 | ||
+ | F; 0; 0 | ||
+ | A*; 1; 0 ; 1 | ||
+ | B*; 1; 0 ; 1 | ||
+ | C*; 1; 0 ; 1 | ||
+ | D*; 1; 0 ; 1 | ||
+ | E*; 1; 0 ; 1 | ||
+ | F*; 1; 0 ; 1 | ||
+ | *; 1; 1.0 | ||
+ | &reactions | ||
+ | # m^2 amu K sigma sticking J/mol | ||
+ | HK; {A} + {*} => {A*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 | ||
+ | HK; {E} + {*} => {E*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 | ||
+ | HK; {F} + {*} => {F*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 | ||
+ | # vf vb Eaf Eab | ||
+ | AR; {A*} => {B*}; 1e13; 1e13; 90e3; 90e3 ;1 | ||
+ | AR; {B*} => {C*}; 1e13; 1e13; 80e3; 90e3 ;1 | ||
+ | AR; {B*} => {D*}; 1e13; 1e13; 70e3; 80e3 ;1 | ||
+ | AR; {C*} => {E*}; 1e13; 1e13; 70e3; 120e3 ;1 | ||
+ | AR; {D*} => {F*}; 1e13; 1e13; 70e3; 90e3 ;1 | ||
+ | &settings | ||
+ | TYPE = SEQUENCERUN | ||
+ | PRESSURE = 1 | ||
+ | REAGENTS = {A} | ||
+ | KEYCOMPONENT = {F} | ||
+ | EACT=1 | ||
+ | ORDERS=1 | ||
+ | DRC=1 | ||
+ | TDRC=0 | ||
+ | DEBUG=1 | ||
+ | USETIMESTAMP=0 | ||
+ | &runs | ||
+ | 300; 1e6; 1e-10; 1e-10 | ||
+ | 350; 1e6; 1e-10; 1e-10 | ||
+ | 400; 1e6; 1e-10; 1e-10 | ||
+ | 450; 1e5; 1e-10; 1e-10 | ||
+ | 500; 1e4; 1e-10; 1e-10 | ||
+ | 550; 1e3; 1e-10; 1e-10 | ||
+ | 600; 1e2; 1e-10; 1e-10 | ||
+ | 650; 1e1; 1e-10; 1e-10 | ||
+ | 700; 1e1; 1e-10; 1e-10 | ||
+ | 750; 1e1; 1e-10; 1e-10 | ||
+ | 800; 1e0; 1e-10; 1e-10 | ||
+ | 850; 1e0; 1e-10; 1e-10 | ||
+ | 900; 1e0; 1e-10; 1e-10 | ||
+ | 950; 1e0; 1e-10; 1e-10 | ||
+ | 1000; 1e0; 1e-10; 1e-10 | ||
+ | 1100; 1e0; 1e-10; 1e-10 | ||
+ | 1200; 1e0; 1e-10; 1e-10 | ||
+ | 1300; 1e0; 1e-10; 1e-10 | ||
+ | 1400; 1e0; 1e-10; 1e-10 | ||
+ | 1500; 1e0; 1e-10; 1e-10 | ||
+ | 1600; 1e0; 1e-10; 1e-10 | ||
+ | 1700; 1e0; 1e-10; 1e-10 | ||
+ | 1800; 1e0; 1e-10; 1e-10 | ||
+ | 1900; 1e0; 1e-10; 1e-10 | ||
+ | 2000; 1e0; 1e-10; 1e-10 | ||
+ | &selectivity | ||
+ | species_balance; {A}; {E},{F} | ||
+ | </pre> |
Revision as of 19:36, 16 July 2017
Introduction
In this example, we show how to calculate selectivities and how to perform a degree of selectivity control analysis.
Input
In this example, the following input file is used:
&compounds A; 0; 1 E; 0; 0 F; 0; 0 A*; 1; 0 ; 1 B*; 1; 0 ; 1 C*; 1; 0 ; 1 D*; 1; 0 ; 1 E*; 1; 0 ; 1 F*; 1; 0 ; 1 *; 1; 1.0 &reactions # m^2 amu K sigma sticking J/mol HK; {A} + {*} => {A*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 HK; {E} + {*} => {E*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 HK; {F} + {*} => {F*}; 1e-20; 1; 1; 1; 1; 80e3 ; 1 # vf vb Eaf Eab AR; {A*} => {B*}; 1e13; 1e13; 90e3; 90e3 ;1 AR; {B*} => {C*}; 1e13; 1e13; 80e3; 90e3 ;1 AR; {B*} => {D*}; 1e13; 1e13; 70e3; 80e3 ;1 AR; {C*} => {E*}; 1e13; 1e13; 70e3; 120e3 ;1 AR; {D*} => {F*}; 1e13; 1e13; 70e3; 90e3 ;1 &settings TYPE = SEQUENCERUN PRESSURE = 1 REAGENTS = {A} KEYCOMPONENT = {F} EACT=1 ORDERS=1 DRC=1 TDRC=0 DEBUG=1 USETIMESTAMP=0 &runs 300; 1e6; 1e-10; 1e-10 350; 1e6; 1e-10; 1e-10 400; 1e6; 1e-10; 1e-10 450; 1e5; 1e-10; 1e-10 500; 1e4; 1e-10; 1e-10 550; 1e3; 1e-10; 1e-10 600; 1e2; 1e-10; 1e-10 650; 1e1; 1e-10; 1e-10 700; 1e1; 1e-10; 1e-10 750; 1e1; 1e-10; 1e-10 800; 1e0; 1e-10; 1e-10 850; 1e0; 1e-10; 1e-10 900; 1e0; 1e-10; 1e-10 950; 1e0; 1e-10; 1e-10 1000; 1e0; 1e-10; 1e-10 1100; 1e0; 1e-10; 1e-10 1200; 1e0; 1e-10; 1e-10 1300; 1e0; 1e-10; 1e-10 1400; 1e0; 1e-10; 1e-10 1500; 1e0; 1e-10; 1e-10 1600; 1e0; 1e-10; 1e-10 1700; 1e0; 1e-10; 1e-10 1800; 1e0; 1e-10; 1e-10 1900; 1e0; 1e-10; 1e-10 2000; 1e0; 1e-10; 1e-10 &selectivity species_balance; {A}; {E},{F}