Monte Carlo Commands

Monte Carlo Commands

Change parameters for Monte Carlo-based calculations. Energy values should be given in the current working unit of energy in the program.

mcaccept

Syntax:

void mcaccept(move,  
 emax); 
string  move;
double  emax;

Sets the energy difference emax for the movetype move above which moves will be rejected.

For example:

mcaccept("translate", 0.0);

requests that only translation moves that lower the overall energy will be accepted.

mcaccept("insert", 200.0);

requests that insertion moves will be accepted provided the total energy does not rise more than 200.0 units.

mcallow

Syntax:

void mcallow(move,  
 allow); 
string  move;
bool  allow;

Indicates whether to allow moves of the specified type in Monte Carlo calculations.

For example:

mcallow("translate", FALSE);

disallows translation moves in Monte Carlo calculations.

mcallow("rotate", "true");

allows rotation moves in Monte Carlo calculations.

mcmaxstep

Syntax:

void mcmaxstep(move,  
 mcmaxstep); 
string  move;
double  mcmaxstep;

Sets the maximal step size for the move type move.

For example:

mcmaxstep("translate", 5.0);

sets the maximum translation displacement to 5 Å.

mcmaxstep("rotate", 30.0);

sets the maximum rotation to 30 degrees.

mcntrials

Syntax:

void mcntrials(move,  
 n); 
string  move;
double  n;

Sets the number of times n that the move type move should be attempted in each cycle.

For example:

mcntrials("insert", 50);

requests that there will be 50 insertion attempts per cycle per molecule type.

printmc

Syntax:

void printmc(); 

Prints the current Monte Carlo parameters.

For example:

printmc();