config.h

   1 /*
   2  * configuration file
   3  *
   4  */
   5
   6 #include "moldyn.h"
   7
   8 // simulation volume
   9
  10 #define LCNTX           31
  11 #define LCNTY           31
  12 #define LCNTZ           31
  13
  14 // initial lattice
  15
  16 #define INIT_SI         // plain silicon
  17 //#define INIT_C                // plain diamond
  18 //#define INIT_3CSIC    // 3C-SiC
  19
  20 // berendsen thermostat
  21
  22 //#define T_SCALE_TAU   10.0
  23 #define T_SCALE_TAU     100.0
  24 #define P_SCALE_TAU     (0.01/(100*GPA))
  25
  26 // prerun
  27
  28 #define PRERUN          500
  29 #define PRE_TAU         1.0
  30
  31 // insertrun
  32
  33 #define INS_RUNS        600
  34 #define INS_ATOMS       10
  35
  36 #define INS_CARBON      /* comment for silicon */
  37
  38 #ifdef INS_CARBON
  39  #define INS_TYPE       C
  40  #define INS_MASS       M_C
  41  #define INS_BRAND      1
  42 #else
  43  #define INS_TYPE       SI
  44  #define INS_MASS       M_SI
  45  #define INS_BRAND      0
  46 #endif
  47
  48 #define INS_RAND        // random nsert
  49 //#define INS_HEXA      // hexagonal interstitial position
  50 //#define INS_TETRA     // tetrahedral interstitial position
  51 //#define INS_110DB     // 110 dumbbell interstitial position
  52 //#define INS_001DB     // 001 dumbbell interstitial position
  53 //#define INS_USER      // user defined insertion location
  54  #define INS_UX         -0.25
  55  #define INS_UY         -0.25
  56  #define INS_UZ         +0.25
  57
  58 //#define INS_ATTR      ATOM_ATTR_HB                            // heat bath
  59 #define INS_ATTR        ATOM_ATTR_HB|ATOM_ATTR_VB               // + vis bonds
  60 //#define INS_ATTR      ATOM_ATTR_HB|ATOM_ATTR_VB|ATOM_ATTR_FP  // + fixed pos
  61
  62 //#define INS_LENX      (12*ALBE_LC_SIC)
  63 //#define INS_LENY      (12*ALBE_LC_SIC)
  64 //#define INS_LENZ      (12*ALBE_LC_SIC)
  65 #define INS_LENX        (9*ALBE_LC_SI)
  66 #define INS_LENY        (9*ALBE_LC_SI)
  67 #define INS_LENZ        (9*ALBE_LC_SI)
  68 #define INS_OFFSET      (ALBE_LC_SI/8.0)
  69 #define INS_DYNAMIC_LEN // with p-ctrl we have a varying volume
  70
  71 #define INS_R_C         1.5
  72 #define INS_DELTA_TC    1.0
  73 #define INS_DELTA_PC    (1.0*BAR)
  74 #define INS_RELAX       100
  75 #define INS_TAU         1.0
  76
  77 // postrun
  78
  79 #define POST_RUNS       430
  80 #define POST_DELTA_TC   1.0
  81 #define POST_DELTA_PC   (1.0*BAR)
  82 #define POST_DT         -1.0
  83 #define POST_RELAX      100
  84 #define POST_TAU        1.0
  85
  86 // log
  87
  88 #define R_TITLE         "Carbon into silicon"
  89 #define LOG_E           10
  90 #define LOG_T           10
  91 #define LOG_P           10
  92 #define LOG_V           10
  93 #define LOG_S           2000
  94 #define LOG_A           2000
  95
  96 #define AVG_SKIP        0