Project

General

Profile

embed_2016.log

log of a successfull simulation with gromacs 2016.3 - Daniel Bauer, 12/22/2017 01:22 PM

 
1
grompp
2
Replacing old mdp entry 'nstxtcout' by 'nstxout-compressed'
3

    
4
Back Off! I just backed up mdout.mdp to ./#mdout.mdp.1#
5

    
6
NOTE 1 [file g_membed.mdp]:
7
  The group cutoff scheme is deprecated since GROMACS 5.0 and will be
8
  removed in a future release when all interaction forms are supported for
9
  the verlet scheme. The verlet scheme already scales better, and it is
10
  compatible with GPUs and other accelerators.
11

    
12

    
13
NOTE 2 [file g_membed.mdp]:
14
  With exact cut-offs, rlist should be larger than rcoulomb and rvdw, so
15
  that there is a buffer region for particle motion between neighborsearch
16
  steps
17

    
18

    
19
NOTE 3 [file g_membed.mdp]:
20
  For energy conservation with switch/shift potentials, rlist should be 0.1
21
  to 0.3 nm larger than rvdw.
22

    
23
Setting the LD random seed to 723256370
24
Generated 97877 of the 97903 non-bonded parameter combinations
25
Generating 1-4 interactions: fudge = 1
26
Generated 64492 of the 97903 1-4 parameter combinations
27
Excluding 3 bonded neighbours molecule type 'kcsa'
28
Excluding 3 bonded neighbours molecule type 'POPC'
29
Excluding 3 bonded neighbours molecule type 'POPC'
30
Excluding 2 bonded neighbours molecule type 'SOL'
31
Excluding 1 bonded neighbours molecule type 'K'
32
Excluding 1 bonded neighbours molecule type 'CL'
33

    
34
NOTE 4 [file system.top, line 19]:
35
  System has non-zero total charge: 4.000003
36
  Total charge should normally be an integer. See
37
  http://www.gromacs.org/Documentation/Floating_Point_Arithmetic
38
  for discussion on how close it should be to an integer.
39
  
40

    
41

    
42
Velocities were taken from a Maxwell distribution at 298 K
43
Number of degrees of freedom in T-Coupling group Other is 63997.00
44
Number of degrees of freedom in T-Coupling group Water_and_ions is 63045.00
45
Number of degrees of freedom in T-Coupling group Protein is 0.00
46

    
47
WARNING 1 [file g_membed.mdp]:
48
  Can not exclude the lattice Coulomb energy between energy groups
49

    
50

    
51
NOTE 5 [file g_membed.mdp]:
52
  The sum of the two largest charge group radii (0.157881) is larger than
53
  rlist (1.200000) - rvdw (1.200000).
54
  With exact cut-offs, better performance can be obtained with
55
  cutoff-scheme = Verlet, because it does not use charge groups at all.
56

    
57
Estimate for the relative computational load of the PME mesh part: 0.07
58

    
59
There were 5 notes
60

    
61
There was 1 warning
62

    
63
Back Off! I just backed up g_membed.tpr to ./#g_membed.tpr.1#
64
turning H bonds into constraints...
65
turning H bonds into constraints...
66
turning H bonds into constraints...
67
turning H bonds into constraints...
68
turning H bonds into constraints...
69
Analysing residue names:
70
There are:   388    Protein residues
71
There are:   200      Other residues
72
There are: 10480      Water residues
73
There are:    56        Ion residues
74
Analysing Protein...
75
Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...
76
Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...
77
Largest charge group radii for Van der Waals: 0.079, 0.079 nm
78
Largest charge group radii for Coulomb:       0.079, 0.079 nm
79
Calculating fourier grid dimensions for X Y Z
80
Using a fourier grid of 72x72x80, spacing 0.108 0.108 0.120
81
This run will generate roughly 6 Mb of data
82
Copy input topology to output topology
83
Run embedding simulation
84

    
85
Running on 1 node with total 8 cores, 16 logical cores, 1 compatible GPU
86
Hardware detected:
87
  CPU info:
88
    Vendor: AMD
89
    Brand:  AMD Ryzen 7 1700X Eight-Core Processor         
90
    SIMD instructions most likely to fit this hardware: AVX2_256
91
    SIMD instructions selected at GROMACS compile time: AVX_256
92

    
93
  Hardware topology: Full, with devices
94
  GPU info:
95
    Number of GPUs detected: 1
96
    #0: NVIDIA GeForce GTX 1070, compute cap.: 6.1, ECC:  no, stat: compatible
97

    
98
Compiled SIMD instructions: AVX_256, GROMACS could use AVX2_256 on this machine, which is better.
99

    
100
Reading file g_membed.tpr, VERSION 2016.3 (single precision)
101

    
102
NOTE: GPU(s) found, but the current simulation can not use GPUs
103
      To use a GPU, set the mdp option: cutoff-scheme = Verlet
104

    
105

    
106
Membrane embedding only supports a single rank. Choosing to use only a single thread-MPI rank.
107
Using 1 MPI thread
108

    
109
1 compatible GPU detected in the system, but none will be used.
110
Consider trying GPU acceleration with the Verlet scheme!
111

    
112
Initializing membed
113
Back Off! I just backed up embed.dat to ./#embed.dat.1#
114

    
115
Setting -cpt to -1, because embedding cannot be restarted from cpt-files.
116

    
117
Select a group to embed in the membrane:
118
Group     0 (         System) has 64175 elements
119
Group     1 (        Protein) has  5879 elements
120
Group     2 (      Protein-H) has  2901 elements
121
Group     3 (        C-alpha) has   388 elements
122
Group     4 (       Backbone) has  1164 elements
123
Group     5 (      MainChain) has  1551 elements
124
Group     6 (   MainChain+Cb) has  1899 elements
125
Group     7 (    MainChain+H) has  1929 elements
126
Group     8 (      SideChain) has  3950 elements
127
Group     9 (    SideChain-H) has  1350 elements
128
Group    10 (    Prot-Masses) has  5879 elements
129
Group    11 (    non-Protein) has 58296 elements
130
Group    12 (          Other) has 26800 elements
131
Group    13 (           POPC) has 26800 elements
132
Group    14 (              K) has    28 elements
133
Group    15 (             CL) has    28 elements
134
Group    16 (          Water) has 31440 elements
135
Group    17 (            SOL) has 31440 elements
136
Group    18 (      non-Water) has 32735 elements
137
Group    19 (            Ion) has    56 elements
138
Group    20 (           POPC) has 26800 elements
139
Group    21 (              K) has    28 elements
140
Group    22 (             CL) has    28 elements
141
Group    23 ( Water_and_ions) has 31496 elements
142
Select a group: 
143
Select a group to embed Protein into (e.g. the membrane):
144
Group     0 (         System) has 64175 elements
145
Group     1 (        Protein) has  5879 elements
146
Group     2 (      Protein-H) has  2901 elements
147
Group     3 (        C-alpha) has   388 elements
148
Group     4 (       Backbone) has  1164 elements
149
Group     5 (      MainChain) has  1551 elements
150
Group     6 (   MainChain+Cb) has  1899 elements
151
Group     7 (    MainChain+H) has  1929 elements
152
Group     8 (      SideChain) has  3950 elements
153
Group     9 (    SideChain-H) has  1350 elements
154
Group    10 (    Prot-Masses) has  5879 elements
155
Group    11 (    non-Protein) has 58296 elements
156
Group    12 (          Other) has 26800 elements
157
Group    13 (           POPC) has 26800 elements
158
Group    14 (              K) has    28 elements
159
Group    15 (             CL) has    28 elements
160
Group    16 (          Water) has 31440 elements
161
Group    17 (            SOL) has 31440 elements
162
Group    18 (      non-Water) has 32735 elements
163
Group    19 (            Ion) has    56 elements
164
Group    20 (           POPC) has 26800 elements
165
Group    21 (              K) has    28 elements
166
Group    22 (             CL) has    28 elements
167
Group    23 ( Water_and_ions) has 31496 elements
168
Select a group: Embedding piece 0 with center of geometry: 3.892237 3.892237 4.296000
169

    
170

    
171
Back Off! I just backed up temp.topyAlhm6 to ./#temp.topyAlhm6.1#
172

    
173
Back Off! I just backed up embedded.top to ./#embedded.top.1#
174

    
175
NOTE: This file uses the deprecated 'group' cutoff_scheme. This will be
176
removed in a future release when 'verlet' supports all interaction forms.
177

    
178

    
179
NOTE: Thread affinity setting failed. This can cause performance degradation.
180
      If you think your settings are correct, ask on the gmx-users list.
181

    
182

    
183
Back Off! I just backed up g_membed.xtc to ./#g_membed.xtc.1#
184

    
185
Back Off! I just backed up g_membed.edr to ./#g_membed.edr.1#
186
starting mdrun 'KcsA in POPC'
187
1000 steps,      2.0 ps.
188
Selected 1: 'Protein'
189
Selected 12: 'Other'
190
The estimated area of the protein in the membrane is 18.833 nm^2
191

    
192
There are 140 lipids in the membrane part that overlaps the protein.
193
The area per lipid is 0.5924 nm^2.
194
Maximum number of lipids that will be removed is 63.
195

    
196
Will resize the protein by a factor of 0.100 in the xy plane and 1.000 in the z direction.
197
This resizing will be done with respect to the geometrical center of all protein atoms
198
that span the membrane region, i.e. z between 1.338 and 7.254
199

    
200
Will remove 0 kcsa molecules
201
Will remove 24 POPC molecules
202
Will remove 29 SOL molecules
203
Will remove 0 K molecules
204
Will remove 0 CL molecules
205

    
206
Writing final coordinates.
207

    
208
               Core t (s)   Wall t (s)        (%)
209
       Time:      319.639      319.639      100.0
210
                 (ns/day)    (hour/ns)
211
Performance:        0.541       44.350