From 60c2b879f26c5f008af0df7d7f159d33077225ac Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@gmail.com> Date: Sat, 09 Mar 2019 19:47:16 +0000 Subject: [PATCH] Cleanup of stretching energy plugin --- src/tape | 40 +++++++++++++++++++++++++++++++++------- 1 files changed, 33 insertions(+), 7 deletions(-) diff --git a/src/tape b/src/tape index 0784fac..42e22ea 100644 --- a/src/tape +++ b/src/tape @@ -12,7 +12,7 @@ # Pressure calculations # (pswitch=1: calc. p*dV energy contribution) -pswitch = 0 +pswitch = 1 # pressure difference: p_inside - p_outside (in units kT/l_min^3): pressure=-10.0 @@ -23,28 +23,35 @@ #Constant area constraint (0 disable constant area, 2 enable constant area with epsarea) constareaswitch=0 + +#Stretching +stretchswitch=1 +xkA0=1.0 + ####### Polymer (brush) definitions ########### # npoly is a number of polymers attached to npoly distinct vertices on vesicle -npoly=3 +npoly=0 # nmono is a number of monomers in each polymer nmono=10 # Spring constant between monomers of the polymer k_spring=800 +#set to 1 if half of the polymeres are inside the vesicle +internal_poly=1 ####### Filament (inside the vesicle) definitions ########### # nfil is a number of filaments inside the vesicle -nfil=2 +nfil=0 # nfono is a number of monomers in each filament -nfono=300 +nfono=3 # Persistence lenght of the filaments (in units l_min) xi=100 ####### Nucleus (inside the vesicle) ########### # Radius of an impenetrable hard sphere inside the vesicle R_nucleus=0 -R_nucleusX=6.0 -R_nucleusY=12.0 -R_nucleusZ=6.0 +R_nucleusX=0 +R_nucleusY=0 +R_nucleusZ=0 ####### Cell definitions ############ nxmax=60 nymax=60 @@ -82,3 +89,22 @@ #max number of processes in distributed (voluntary) environment distributed_processes=50 #cuda options??? + + +#number of vertices with spontaneous curvature (integer) +number_of_vertices_with_c0=100 +#c0/2 is spontaneous curvature. c0 is used as (c1+c1-c0)^2 in energy term (float) +c0=0.5 +#energy of attraction of vertices with spontaneous curvature (float, positive value for attraction) +w=10.0 +#direct force on vesicles with spontaneous curvature (float) +F=2.0 + + + +#plane confinement +plane_confinement_switch=1 +#final plane distance (float in lmin) +plane_d=10 +#plane to vesicle repulsion force while closing +plane_F=10 -- Gitblit v1.9.3