1
# Here we define a trivial molecule containing only one particle.
5
# atomID molID atomType charge x y z
8
$atom:C $mol @atom:C 0.0 0.0 0.0 0.0
11
write_once("Data Masses") {
15
write_once("In Settings") {
16
# If for some reason there are multiple chaperones present,
17
# I assume that they interact repulsively (hence, L=0)
19
# i j pairStyle eps sig K L
21
pair_coeff @atom:C @atom:C lj/charmm/coul/charmm/inter 1.0 3.0 1 0
23
# Optional: define the atoms in the "chaperonins" group:
24
# (Defining a group for the chaperone makes it easy to immobilize it later.)
26
group chaperones type @atom:C
29
# Specify which pair_styles, and atom styles work well with
30
# this model. (Again this can be overridden later.)
32
write_once("In Init") {
35
pair_style hybrid lj/charmm/coul/charmm/inter es4k4l maxmax 8.0 9.0
40
# We have not specified how this particle interacts with other particles
41
# besides itself. Later on you must do this.