Difference between revisions of "Keywords and settings new"

From MKMCXX
Jump to: navigation, search
Line 46: Line 46:
 
| 0
 
| 0
 
| 0 or 1
 
| 0 or 1
 +
|-
 +
| ABSTOL
 +
| Set the absolute tolerance for a TPD run
 +
| 1e-12
 +
| positive float
 +
|-
 +
| RELTOL
 +
| Set the relative tolerance for a TPD run
 +
| 1e-12
 +
| positive float
 +
|-
 +
| CACHE_WORKPOINT
 +
| Reuse the sequencerun results as a workpoint in additional routines like DRC.
 +
* 0 : Do not reuse. Rerun workpoint with zero perturbation.
 +
* 1 : Reuse the sequencerun results and skip the zero perturbation workpoint.
 +
| 1
 +
| 0 or 1
 +
|-
 +
| NUMDIFF
 +
| Resolution of perturbations at each side of the workpoint.
 +
| 2
 +
| int >= 1
 +
|-
 +
| ORDERSDIFF
 +
| Step size used in the linear fitting of reaction orders (fractional pressure)
 +
* Wide: 0.1    Normal: 0.01  Tight: 0.001
 +
| 0.01
 +
| float
 +
|-
 +
| EACTDIFF
 +
| Step size used in the linear fitting of the apparent activation energy (fractional temperature)
 +
* Wide: 0.001  Normal: 0.0001 Tight: 0.00001
 +
| 0.0001
 +
|-
 +
| DRCDIFF
 +
| Step size used in the linear fitting for the degree of rate control analysis (fractional k)
 +
* Wide: 0.1    Normal: 0.01  Tight: 0.001
 +
| 0.01
 +
|-
 +
| TDRCDIFF
 +
| Step size used in the linear fitting for the thermodynamic DRC analysis (absolute dG in J/mol)
 +
Wide: 10      Normal: 1    Tight: 0.1
 +
| 1
 
|-
 
|-
 
| PRECISION
 
| PRECISION
Line 86: Line 129:
 
| Whether to create plots
 
| Whether to create plots
 
| 1
 
| 1
|-
 
| ORDERSDIFF
 
| Step size used in the linear fitting of reaction orders
 
| 0.2 (Wide: 0.4; Normal: 0.2; Tight: 0.002)
 
|-
 
| EACTDIFF
 
| Step size used in the linear fitting of the apparent activation energy
 
| 0.002 (Wide: 0.02; Normal: 0.002; Tight: 0.00002)
 
|-
 
| DRCDIFF
 
| Step size used in the linear fitting for the degree of rate control analysis
 
| 0.002 (Wide: 0.02; Normal: 0.002; Tight: 0.00002)
 
|-
 
| TDRCDIFF
 
| Step size used in the linear fitting for the degree of rate control analysis
 
| 0.002 (Wide: 0.02; Normal: 0.002; Tight: 0.00002)
 
 
|-
 
|-
 
| BOOSTER
 
| BOOSTER

Revision as of 10:35, 23 March 2020

Settings block

Mkmcxx2.15.3.png

The list below is an overview of all the keywords that can be placed inside the &settings block. All keywords have a default value, so if these are not specified within the block, then the default value is used.

List of common keywords and their meaning
Keyword Explanation Default value Possible values
DEBUG Output additional debug statements.
  • 0 : no debug statements
  • 1-4 : various levels of debug statements
0 0-4 (int)
PRESSURE The total pressure in the gas phase.
  • <0 : Calculate the total pressure from the sum of the starting pressure of the compounds.
  • >0 : Set this total pressure in bar and scale starting pressures proportionately.
-1 (Take starting pressures) 0 - inf or negative (float)
List of less common keywords and their meaning
Keyword Explanation Default value Possible values
RERUN_OUTPUT Define the folder used for re-running sensitivity analysis and graphs run 'string'
RERUN_GRAPH Whether to re-run graphs from the RERUN_OUTPUT folder
  • 0 : Regular mkmcxx run
  • 1 : Re-run graphs
0 0 or 1
ABSTOL Set the absolute tolerance for a TPD run 1e-12 positive float
RELTOL Set the relative tolerance for a TPD run 1e-12 positive float
CACHE_WORKPOINT Reuse the sequencerun results as a workpoint in additional routines like DRC.
  • 0 : Do not reuse. Rerun workpoint with zero perturbation.
  • 1 : Reuse the sequencerun results and skip the zero perturbation workpoint.
1 0 or 1
NUMDIFF Resolution of perturbations at each side of the workpoint. 2 int >= 1
ORDERSDIFF Step size used in the linear fitting of reaction orders (fractional pressure)
  • Wide: 0.1 Normal: 0.01 Tight: 0.001
0.01 float
EACTDIFF Step size used in the linear fitting of the apparent activation energy (fractional temperature)
  • Wide: 0.001 Normal: 0.0001 Tight: 0.00001
0.0001
DRCDIFF Step size used in the linear fitting for the degree of rate control analysis (fractional k)
  • Wide: 0.1 Normal: 0.01 Tight: 0.001
0.01
TDRCDIFF Step size used in the linear fitting for the thermodynamic DRC analysis (absolute dG in J/mol)

Wide: 10 Normal: 1 Tight: 0.1

1
PRECISION Amount of significant digits to use in output. 10
USETIMESTAMP Whether to place the output in a new folder which name is based on the current time stamp. If 0, then all output is placed inside a 'RUN' folder. 1
ORDERS Whether to calculate the reaction orders. (You need to set KEYCOMPONENTS for this) 1
EACT Whether to calculate the apparent activation energy. (You need to set KEYCOMPONENTS for this) 1
DRC Whether to perform a degree of rate control analysis. (You need to set KEYCOMPONENTS for this) 1
TDRC Whether to perform a thermodynamic degree of rate control analysis. (You need to set KEYCOMPONENTS for this) 1
REAGENTS List of reagents for which the reaction orders will be calculated NULL
KEYCOMPONENTS List of compounds on which the reaction orders, apparent activation energy, DRC, DSC and/or TDRC analysis should be based NULL
GNUPLOT Whether to output GNUPlot-style graphs 0
MAKEPLOTS Whether to create plots 1
BOOSTER Default multiplier used to speed-up reaction rates; sometimes leads to faster convergence towards the steady-state solution 1.0
SOLVERTYPE Type of integration method to use for solving the system of ordinary differential equations 1.0 (BDF)
SOLSTOPTIME Specificies when the solver should re-evaluate the Jacobian matrix. A value of 1 means at every time step. 1.0
SOLMAXSTEP Maximum number of internal steps the solver is allowed to take. 5000
SOLTESTFAIL Maximum number of test failures before the solver gives up 70
SOLCONVFAIL Maximum number of convergence failures before the solver gives up 100

Graphs block

Mkmcxx 2 3 0.png

Inside the graphs block, you can set the colors used in the non-GNUPLOT graphs for specific components. On each line, you place the compound between curly brackets followed by the RGB color code. For example:

&graphs
# fix colors for particular compounds
{A*} #E74C3D
{B*} #F29C1F
{C*} #287FB9
{*}  #15A086

Selectivity block

Mkmcxx 2 3 0.png

Inside the selectivity block, you can specify mole balances on which basis you calculate selectivity and degree of selectivity control graphs. You need to specify a name, a key component and one or more products (typically more than one, else the concept of selectivity is rather trivial). The name will be used in the generation of the corresponding graph file, so please use a safe name (i.e. no spaces and no special characters!)

&selectivity
species_balance; {A}; {E},{F}

In the above example, the name of the mole balance block is "species_balance", the key component is A and the products of interest are E and F. A detailed example on how to use this block is explained here.

Stoichiometry

Often, the stoichiometry of the key component and the products is not simply a 1:1 ratio. To account for this difference, you need to put the stoichiometric coefficient in front the compound. The stoichiometric coefficient is simply the number of key components which need to be consumed to produce one product. To illustrate this, below an example for CO hydrogenation towards C1-C3 hydrocarbons is provided.

&selectivity
carbon_balance; {CO}; {CH4}, 2{CH2CH2}, 2{CH3CH3}, 3{CH3CH2CH3}, 3{CH2CHCH3}