Trisurf Monte Carlo simulator
Samo Penic
2016-12-06 0634b97e436e480b203a22e1113c473d3ae4d281
commit | author | age
a6b1b5 1 ####### Vesicle definitions ###########
SP 2 # nshell is a number of divisions of dipyramid
0634b9 3 <<<<<<< HEAD
e98482 4 nshell=17
0634b9 5 =======
033407 6 nshell=10
0634b9 7 >>>>>>> nirgov
e86357 8 # dmax is the max. bond length (in units l_min)
304510 9 dmax=1.7
ea1cce 10 # dmin_interspecies in the min. dist. between different vertex species (in units l_min)
fe24d2 11 dmin_interspecies=1.2
e86357 12 # bending rigidity of the membrane (in units kT)
2dcc65 13 xk0=10.0
e86357 14 # max step size (in units l_min)
a6b1b5 15 stepsize=0.15
e86357 16
M 17 # Pressure calculations
18 # (pswitch=1: calc. p*dV energy contribution)
45c708 19 pswitch = 0
2dcc65 20 # pressure difference: p_inside - p_outside (in units kT/l_min^3):
e98482 21 pressure=10.0
a6b1b5 22
267433 23 #Constant volume constraint (0 disable constant volume, 1 enable wiht additional vertex move, 2 enable with epsvol)
37791b 24 constvolswitch=0
dd5aca 25 constvolprecision=1e-14
9166cb 26
c0ae90 27 #Constant area constraint (0 disable constant area, 2 enable constant area with epsarea)
8db203 28 constareaswitch=0
9166cb 29
624f81 30 ####### Polymer (brush) definitions ###########
a2db52 31 # npoly is a number of polymers attached to npoly distinct vertices on vesicle
e98482 32 npoly=800
48bb92 33 # nmono is a number of monomers in each polymer
e98482 34 nmono=6
48bb92 35 # Spring constant between monomers of the polymer
M 36 k_spring=800
e98482 37 #set to 1 if half of the polymeres are inside the vesicle
SP 38 internal_poly=1
a6b1b5 39
624f81 40 ####### Filament (inside the vesicle) definitions ###########
58230a 41 # nfil is a number of filaments inside the vesicle
8db203 42 nfil=0
58230a 43 # nfono is a number of monomers in each filament
8db203 44 nfono=3
b30f45 45 # Persistence lenght of the filaments (in units l_min)
5c64e2 46 xi=100
58230a 47
M 48 ####### Nucleus (inside the vesicle) ###########
49 # Radius of an impenetrable hard sphere inside the vesicle
f06f5f 50 R_nucleus=0
bfa6c4 51 R_nucleusX=0
SP 52 R_nucleusY=0
53 R_nucleusZ=0
a6b1b5 54 #######  Cell definitions ############
SP 55 nxmax=60
56 nymax=60
57 nzmax=60
58
59
60 ####### Program Control ############
d7a113 61 #how many MC sweeps between subsequent records of states to disk
a5e543 62 #200000
032273 63 mcsweeps=2000
d7a113 64 #how many initial mcsweeps*inititer MC sweeps before recording to disk?
a5e543 65 #2
a17e91 66 inititer=10
d7a113 67 #how many records do you want on the disk iteration are there in a run?
a5e543 68 #10000
M 69 iterations=100
d7a113 70
a6b1b5 71
dc77e8 72 ###### Spherical harmonics ###########
5a3862 73 # If 0 then spherical harmonics are not calculated at all.
c60a49 74 spherical_harmonics_coefficients=21
dc77e8 75
a6b1b5 76 #shut up if we are using cluster!!!
SP 77 quiet=false
5c64e2 78  
SP 79
f74313 80
SP 81 #what type of multiprocessing? (*none, smp, cluster, distributed, cuda, auto)
82 #currently only none makes sense.
83 multiprocessing=none
84 #how many cores are allowed to process in SMP?
85 smp_cores=2
86 #how many nodes in cluster?
87 cluster_nodes=50
88 #max number of processes in distributed (voluntary) environment
89 distributed_processes=50
90 #cuda options???
e5858f 91
SP 92
93 #NirGov branch only!
94 #number of vertices with spontaneous curvature (integer)
250de4 95 number_of_vertices_with_c0=100
e5858f 96 #spontaneous curvature (float)
032273 97 c0=0.5
SP 98 #energy of attraction of vertices with spontaneous curvature (float, positive value for attraction)
250de4 99 w=10.0
SP 100 #direct force on vesicles with spontaneous curvature (float)
101 F=2.0