]> www.hackdaworld.org Git - physik/posic.git/blob - config.h
support bonds for more than one atom
[physik/posic.git] / 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