[gmx-users] RDF converge to zero than 1, why?
Dallas Warren
Dallas.Warren at monash.edu
Thu Sep 16 00:36:42 CEST 2010
What does the system actually look like? Can you provide an image
showing how the molecules are distributed?
As already noted, it simply sounds like your surfactant has aggregated,
so at large distances away from the aggregate the probability is zero.
Catch ya,
Dr. Dallas Warren
Medicinal Chemistry and Drug Action
Monash Institute of Pharmaceutical Sciences, Monash University
381 Royal Parade, Parkville VIC 3010
dallas.warren at monash.edu
+61 3 9909 9304
---------------------------------
When the only tool you own is a hammer, every problem begins to resemble
a nail.
> -----Original Message-----
> From: gmx-users-bounces at gromacs.org [mailto:gmx-users-
> bounces at gromacs.org] On Behalf Of teklebrh at ualberta.ca
> Sent: Thursday, 16 September 2010 4:34 AM
> To: gmx-users at gromacs.org
> Subject: Re: [gmx-users] RDF converge to zero than 1, why?
>
>
> Dear Justine,
>
> > Then my guess is you don't have a PBC problem, but perhaps a
> > convergence problem. You're only analyzing 1 ns of data; I don't
> > know what implication that will have for your outcome.
>
> I have read about this problem and some people suggested that my box
> expanded a lot means I have a gas phase. I check the density It was
> fine.
>
> I have done the same kind of analysis using the same procedure for
> part of the molecule and it converged very well to 1. But for the
> entire surfactant it did not converge to one (approx to zero)...... Do
> you think this is an artifact or else.
>
> If I discard these values below 1 .... what is the implication? What
> will be then my justification for that.
>
> I will be happy as well if other people comment on this topic.
>
> WHY my RDF converge to zero not 1.
>
> Rob
>
>
>
>
>
>
>
> Quoting "Justin A. Lemkul" <jalemkul at vt.edu>:
>
> >
> >
> > teklebrh at ualberta.ca wrote:
> >>
> >> Dear Justine,
> >>
> >> I have tried every thing but the g_rdf of my surfactant molecules
> ....
> >> Surfactant = PAP
> >>
> >> PAP-PAP .. RDF looks like this..
> >> 0 0
> >> 0.002 0
> >> 0.004 0
> >> 0.006 0
> >> 0.008 0
> >> 0.01 0
> >> 0.012 0
> >> 0.014 0
> >> 0.016 0
> >> 0.018 0
> >> 0.02 0
> >> 0.022 0
> >> 0.024 0
> >> 0.026 0
> >> 0.028 0
> >> 0.03 0
> >> 0.032 0
> >> 0.034 0
> >> 0.036 0
> >> 0.038 0
> >> 0.04 0
> >> 0.042 0
> >> 0.044 0
> >> 0.046 0
> >> 0.048 0
> >> 0.05 0
> >> 0.052 0
> >> 0.054 0
> >> 0.056 0
> >> 0.058 0
> >> 0.06 0
> >> 0.062 0
> >> 0.064 0
> >> 0.066 0
> >> 0.068 0
> >> 0.07 0
> >> 0.072 0
> >> 0.074 0
> >> 0.076 0
> >> 0.078 0
> >> 0.08 0
> >> 0.082 0
> >> 0.084 0
> >> 0.086 0
> >> 0.088 0
> >> 0.09 0
> >> 0.092 0
> >> 0.094 0
> >> 0.096 0
> >> 0.098 0
> >> 0.1 0
> >> 0.102 0
> >> 0.104 0
> >> 0.106 0
> >> 0.108 0
> >> 0.11 0
> >> 0.112 0
> >> 0.114 0
> >> 0.116 0
> >> 0.118 0
> >> 0.12 0
> >> 0.122 0
> >> 0.124 0.0386708
> >> 0.126 0
> >> 0.128 0
> >> 0.13 0
> >> 0.132 0
> >> 0.134 0
> >> 0.136 0
> >> 0.138 0.0312225
> >> 0.14 0
> >> 0.142 0.0294884
> >> 0.144 0
> >> 0.146 0
> >> 0.148 0
> >> 0.15 0
> >> 0.152 0
> >> 0.154 0.050144
> >> 0.156 0
> >> 0.158 0.047637
> >> 0.16 0.0464535
> >> 0.162 0.022657
> >> 0.164 0.0221076
> >> 0.166 0.0215781
> >> 0.168 0.0632021
> >> 0.17 0.0617243
> >> 0.172 0.0803955
> >> 0.174 0.0392791
> >> 0.176 0.0767833
> >> 0.178 0.112601
> >> 0.18 0.165168
> >> 0.182 0.0718038
> >> 0.184 0.12294
> >> 0.186 0.171872
> >> 0.188 0.117764
> >> 0.19 0.0988267
> >> 0.192 0.0645195
> >> 0.194 0.205385
> >> 0.196 0.17026
> >> 0.198 0.121336
> >> 0.2 0.178384
> >> 0.202 0.145723
> >> 0.204 0.257183
> >> 0.206 0.196166
> >> 0.208 0.261133
> >> 0.21 0.350562
> >> 0.212 0.2646
> >> 0.214 0.324599
> >> 0.216 0.382335
> >> 0.218 0.450424
> >> 0.22 0.491411
> >> 0.222 0.591186
> >> 0.224 0.533271
> >> 0.226 0.651932
> >> 0.228 0.834996
> >> 0.23 1.19148
> >> 0.232 1.95537
> >> 0.234 2.38904
> >> 0.236 3.52306
> >> 0.238 4.37741
> >> 0.24 5.56415
> >> 0.242 6.89399
> >> 0.244 8.25963
> >> 0.246 10.1401
> >> 0.248 10.5862
> >> 0.25 11.5687
> >> 0.252 11.9851
> >> 0.254 12.3869
> >> 0.256 11.786
> >> 0.258 11.3537
> >> 0.26 10.8894
> >> 0.262 10.1434
> >> 0.264 9.45289
> >> 0.266 8.20198
> >> 0.268 7.95582
> >> 0.27 7.26746
> >> 0.272 7.00034
> >> 0.274 6.43903
> >> 0.276 6.16654
> >> 0.278 6.10889
> >> 0.28 6.47701
> >> 0.282 6.68453
> >> 0.284 7.12153
> >> 0.286 7.97467
> >> 0.288 8.63126
> >> 0.29 9.60143
> >> 0.292 10.5722
> >> 0.294 11.6465
> >> 0.296 12.6569
> >> 0.298 14.349
> >> 0.3 15.1626
> >> 0.302 16.1753
> >> 0.304 16.5934
> >> 0.306 17.2663
> >> 0.308 18.0833
> >> 0.31 18.3891
> >> 0.312 19.2229
> >> 0.314 19.0814
> >> 0.316 19.8116
> >> 0.318 19.5628
> >> 0.32 21.2178
> >> 0.322 20.5652
> >> 0.324 21.5411
> >> 0.326 21.5295
> >> 0.328 21.4391
> >> 0.33 22.2337
> >> 0.332 22.1936
> >> 0.334 22.6958
> >> 0.336 22.5159
> >> 0.338 23.9627
> >> 0.34 24.1445
> >> 0.342 24.5645
> >> 0.344 25.3701
> >> 0.346 25.3908
> >> 0.348 26.6904
> >> 0.35 27.0803
> >> 0.352 27.5124
> >> 0.354 30.0258
> >> 0.356 30.3648
> >> 0.358 30.5508
> >> 0.36 31.3087
> >> 0.362 32.5661
> >> 0.364 33.0299
> >> 0.366 33.3802
> >> 0.368 33.9448
> >> 0.37 35.4638
> >> 0.372 36.3684
> >> 0.374 38.6968
> >> 0.376 42.1603
> >> 0.378 46.8835
> >> 0.38 55.6397
> >> 0.382 63.6524
> >> 0.384 70.6934
> >> 0.386 72.4445
> >> 0.388 68.5321
> >> 0.39 63.1321
> >> 0.392 59.9827
> >> 0.394 59.5433
> >> 0.396 59.5346
> >> 0.398 59.6736
> >> 0.4 60.1412
> >> 0.402 63.2898
> >> 0.404 68.8691
> >>
> >> -
> >> -
> >> -
> >> -
> >> -
> >> -
> >>
> >> 3.618 1.49702
> >> 3.62 1.49392
> >> 3.622 1.49265
> >> 3.624 1.47739
> >> 3.626 1.48502
> >> 3.628 1.49045
> >> 3.63 1.47957
> >> 3.632 1.48315
> >> 3.634 1.47226
> >> 3.636 1.47088
> >> 3.638 1.47169
> >> 3.64 1.47815
> >> 3.642 1.45689
> >> 3.644 1.44958
> >> 3.646 1.46132
> >> 3.648 1.45749
> >> 3.65 1.45743
> >> 3.652 1.45316
> >> 3.654 1.44515
> >> 3.656 1.43456
> >> 3.658 1.45053
> >> 3.66 1.42725
> >> 3.662 1.42837
> >> 3.664 1.42375
> >> 3.666 1.43557
> >> 3.668 1.41623
> >> 3.67 1.41688
> >> 3.672 1.4284
> >> 3.674 1.39468
> >> 3.676 1.41252
> >> 3.678 1.41158
> >> 3.68 1.42698
> >> 3.682 1.39971
> >> 3.684 1.40453
> >> 3.686 1.40363
> >> 3.688 1.38676
> >> 3.69 1.3831
> >> 3.692 1.39389
> >> 3.694 1.38261
> >> 3.696 1.38849
> >> 3.698 1.39562
> >> 3.7 1.37748
> >> 3.702 1.3559
> >> 3.704 1.38009
> >> 3.706 1.37901
> >> 3.708 1.36219
> >> 3.71 1.36057
> >> 3.712 1.34968
> >> 3.714 1.35544
> >> 3.716 1.35059
> >> 3.718 1.34517
> >> 3.72 1.35307
> >> 3.722 1.33468
> >> 3.724 1.33697
> >> 3.726 1.32972
> >> 3.728 1.35317
> >> 3.73 1.34359
> >> 3.732 1.33036
> >> 3.734 1.32676
> >> 3.736 1.32371
> >> 3.738 1.31825
> >> 3.74 1.30929
> >> 3.742 1.31097
> >> 3.744 1.29373
> >> 3.746 1.30557
> >> 3.748 1.31186
> >> 3.75 1.30615
> >> 3.752 1.29304
> >> 3.754 1.30105
> >> 3.756 1.28575
> >> 3.758 1.29182
> >> 3.76 1.27846
> >> 3.762 1.29297
> >> 3.764 1.27195
> >> 3.766 1.27255
> >> 3.768 1.27472
> >> 3.77 1.26995
> >> 3.772 1.26761
> >> 3.774 1.24985
> >> 3.776 1.26267
> >> 3.778 1.25907
> >> 3.78 1.25731
> >> 3.782 1.24681
> >> 3.784 1.24505
> >> 3.786 1.25134
> >> 3.788 1.24811
> >> 3.79 1.24578
> >> 3.792 1.23609
> >> 3.794 1.23061
> >> 3.796 1.22387
> >> 3.798 1.23296
> >> 3.8 1.22636
> >> 3.802 1.23006
> >> 3.804 1.19768
> >> 3.806 1.22099
> >> 3.808 1.20974
> >> 3.81 1.21327
> >> 3.812 1.20263
> >> 3.814 1.19753
> >> 3.816 1.19019
> >> 3.818 1.19533
> >> 3.82 1.19154
> >> 3.822 1.19894
> >> 3.824 1.19425
> >> 3.826 1.19786
> >> 3.828 1.18763
> >> 3.83 1.1856
> >> 3.832 1.17327
> >> 3.834 1.1841
> >> 3.836 1.1705
> >> 3.838 1.16918
> >> 3.84 1.16837
> >> 3.842 1.16476
> >> 3.844 1.16393
> >> 3.846 1.15772
> >> 3.848 1.15619
> >> 3.85 1.14672
> >> 3.852 1.15073
> >> 3.854 1.14239
> >> 3.856 1.13708
> >> 3.858 1.12996
> >> 3.86 1.14463
> >> 3.862 1.11898
> >> 3.864 1.12187
> >> 3.866 1.11642
> >> 3.868 1.12246
> >> 3.87 1.12548
> >> 3.872 1.11631
> >> 3.874 1.11326
> >> 3.876 1.10855
> >> 3.878 1.11599
> >> 3.88 1.1087
> >> 3.882 1.09331
> >> 3.884 1.09874
> >> 3.886 1.09685
> >> 3.888 1.09971
> >> 3.89 1.10151
> >> 3.892 1.09531
> >> 3.894 1.08087
> >> 3.896 1.07741
> >> 3.898 1.08176
> >> 3.9 1.06993
> >> 3.902 1.07141
> >> 3.904 1.06588
> >> 3.906 1.06933
> >> 3.908 1.06166
> >> 3.91 1.0674
> >> 3.912 1.0484
> >> 3.914 1.05787
> >> 3.916 1.0477
> >> 3.918 1.04269
> >> 3.92 1.04876
> >> 3.922 1.03885
> >> 3.924 1.04079
> >> 3.926 1.04169
> >> 3.928 1.04109
> >> 3.93 1.0364
> >> 3.932 1.03547
> >> 3.934 1.01927
> >> 3.936 1.02606
> >> 3.938 1.0213
> >> 3.94 1.02079
> >> 3.942 1.00924
> >> 3.944 1.00848
> >> 3.946 1.00971
> >> 3.948 1.01064
> >> 3.95 0.996613
> >> 3.952 0.996559
> >> 3.954 0.996959
> >> 3.956 0.998992
> >> 3.958 0.98474
> >> 3.96 0.988146
> >> 3.962 0.99238
> >> 3.964 0.998688
> >> 3.966 0.982488
> >> 3.968 0.964017
> >> 3.97 0.970145
> >> 3.972 0.974564
> >> 3.974 0.966096
> >> 3.976 0.96833
> >> 3.978 0.957817
> >> 3.98 0.960695
> >> 3.982 0.948339
> >> 3.984 0.954437
> >> 3.986 0.949187
> >> 3.988 0.949588
> >> 3.99 0.932854
> >> 3.992 0.943381
> >> 3.994 0.943342
> >> 3.996 0.935859
> >> 3.998 0.925329
> >> 4 0.930935
> >> 4.002 0.930638
> >> 4.004 0.915742
> >> 4.006 0.916356
> >> 4.008 0.914661
> >> 4.01 0.913785
> >> 4.012 0.912741
> >> 4.014 0.914998
> >> 4.016 0.915678
> >> 4.018 0.915294
> >> 4.02 0.90147
> >> 4.022 0.895573
> >> 4.024 0.897313
> >> 4.026 0.881787
> >> 4.028 0.887577
> >> 4.03 0.890193
> >> 4.032 0.892145
> >> 4.034 0.88283
> >> 4.036 0.885348
> >> 4.038 0.876876
> >> 4.04 0.872154
> >> 4.042 0.876419
> >> 4.044 0.878118
> >> 4.046 0.865789
> >> 4.048 0.870954
> >> 4.05 0.868049
> >> 4.052 0.862877
> >> 4.054 0.858095
> >> 4.056 0.851822
> >> 4.058 0.855818
> >> 4.06 0.853619
> >> 4.062 0.842873
> >> 4.064 0.84844
> >> 4.066 0.845933
> >> 4.068 0.837848
> >> 4.07 0.833639
> >> 4.072 0.837754
> >> 4.074 0.830918
> >> 4.076 0.830175
> >> 4.078 0.825299
> >> 4.08 0.818966
> >> 4.082 0.820927
> >> 4.084 0.818022
> >> 4.086 0.812026
> >> 4.088 0.804889
> >> 4.09 0.815197
> >> 4.092 0.807292
> >> 4.094 0.810965
> >> 4.096 0.800183
> >> 4.098 0.793872
> >> 4.1 0.807127
> >> 4.102 0.789674
> >> 4.104 0.795408
> >> 4.106 0.798583
> >> 4.108 0.79196
> >> 4.11 0.786059
> >> 4.112 0.787969
> >> 4.114 0.775823
> >> 4.116 0.779488
> >> 4.118 0.777093
> >> 4.12 0.781073
> >> 4.122 0.771884
> >> 4.124 0.775363
> >> 4.126 0.766022
> >> 4.128 0.766019
> >> 4.13 0.763528
> >> 4.132 0.759282
> >> 4.134 0.761387
> >> 4.136 0.750792
> >> 4.138 0.755269
> >> 4.14 0.748107
> >> 4.142 0.744142
> >> 4.144 0.741016
> >> 4.146 0.746871
> >> 4.148 0.737341
> >> 4.15 0.730139
> >> 4.152 0.727265
> >> 4.154 0.742743
> >> 4.156 0.726984
> >> 4.158 0.725515
> >> 4.16 0.711536
> >> 4.162 0.731623
> >> 4.164 0.703933
> >> 4.166 0.721472
> >> 4.168 0.716057
> >> 4.17 0.720965
> >> 4.172 0.719051
> >> 4.174 0.704044
> >> 4.176 0.701346
> >> 4.178 0.696225
> >> 4.18 0.700132
> >> 4.182 0.701525
> >> 4.184 0.697064
> >> 4.186 0.700028
> >> 4.188 0.697499
> >> 4.19 0.687055
> >> 4.192 0.677398
> >> 4.194 0.682056
> >> 4.196 0.684369
> >> 4.198 0.678339
> >> 4.2 0.673119
> >> 4.202 0.67668
> >> 4.204 0.678804
> >> 4.206 0.668657
> >> 4.208 0.666873
> >> 4.21 0.661783
> >> 4.212 0.662844
> >> 4.214 0.666264
> >> 4.216 0.664936
> >> 4.218 0.659085
> >> 4.22 0.657377
> >> 4.222 0.647065
> >> 4.224 0.643763
> >> 4.226 0.65035
> >> 4.228 0.64951
> >> 4.23 0.644331
> >> 4.232 0.642256
> >> 4.234 0.638672
> >> 4.236 0.641769
> >> 4.238 0.645202
> >> 4.24 0.637489
> >> 4.242 0.628153
> >> 4.244 0.624075
> >> 4.246 0.627562
> >> 4.248 0.624409
> >> 4.25 0.622008
> >> 4.252 0.617578
> >> 4.254 0.618082
> >> 4.256 0.616705
> >> 4.258 0.620233
> >> 4.26 0.613723
> >> 4.262 0.613542
> >> 4.264 0.606722
> >> 4.266 0.601474
> >> 4.268 0.602512
> >> 4.27 0.610145
> >> 4.272 0.603124
> >> 4.274 0.598812
> >> 4.276 0.588374
> >> 4.278 0.59133
> >> 4.28 0.586432
> >> 4.282 0.587866
> >> 4.284 0.585618
> >> 4.286 0.573144
> >> 4.288 0.57663
> >> 4.29 0.578556
> >> 4.292 0.57521
> >> 4.294 0.580931
> >> 4.296 0.56578
> >> 4.298 0.561448
> >> 4.3 0.573121
> >> 4.302 0.573457
> >> 4.304 0.563928
> >> 4.306 0.560218
> >> 4.308 0.555828
> >> 4.31 0.556358
> >> 4.312 0.563511
> >> 4.314 0.557304
> >> 4.316 0.561251
> >> 4.318 0.545093
> >> 4.32 0.550265
> >> 4.322 0.540745
> >> 4.324 0.536129
> >> 4.326 0.552832
> >> 4.328 0.548235
> >> 4.33 0.537861
> >> 4.332 0.535545
> >> 4.334 0.528127
> >> 4.336 0.535767
> >> 4.338 0.529951
> >> 4.34 0.528368
> >> 4.342 0.523278
> >> 4.344 0.520704
> >> 4.346 0.527159
> >> 4.348 0.52593
> >> 4.35 0.525731
> >> 4.352 0.519521
> >> 4.354 0.513252
> >> 4.356 0.51457
> >> 4.358 0.518719
> >> 4.36 0.510497
> >> 4.362 0.509598
> >> 4.364 0.506181
> >> 4.366 0.504419
> >> 4.368 0.497795
> >> 4.37 0.500966
> >> 4.372 0.504876
> >> 4.374 0.489321
> >> 4.376 0.507235
> >> 4.378 0.493448
> >> 4.38 0.496172
> >> 4.382 0.491397
> >> 4.384 0.482856
> >> 4.386 0.491271
> >> 4.388 0.480308
> >> 4.39 0.481918
> >> 4.392 0.488769
> >> 4.394 0.475618
> >> 4.396 0.478773
> >> 4.398 0.479026
> >> 4.4 0.476978
> >> 4.402 0.470825
> >> 4.404 0.481324
> >> 4.406 0.466229
> >> 4.408 0.46433
> >> 4.41 0.463153
> >> 4.412 0.461773
> >> 4.414 0.452375
> >> 4.416 0.457859
> >> 4.418 0.460082
> >> 4.42 0.454811
> >> 4.422 0.463698
> >> 4.424 0.45372
> >> 4.426 0.451355
> >> 4.428 0.4544
> >> 4.43 0.444576
> >> 4.432 0.445647
> >> 4.434 0.440643
> >> 4.436 0.442069
> >> 4.438 0.442747
> >> 4.44 0.437023
> >> 4.442 0.434889
> >> 4.444 0.440704
> >> 4.446 0.432446
> >> 4.448 0.432088
> >> 4.45 0.428407
> >> 4.452 0.428726
> >> 4.454 0.427902
> >> 4.456 0.426622
> >> 4.458 0.420097
> >> 4.46 0.424059
> >> 4.462 0.423266
> >> 4.464 0.425177
> >> 4.466 0.418946
> >> 4.468 0.417772
> >> 4.47 0.417076
> >> 4.472 0.413366
> >> 4.474 0.405287
> >> 4.476 0.406673
> >> 4.478 0.406363
> >> 4.48 0.400713
> >> 4.482 0.403016
> >> 4.484 0.409838
> >> 4.486 0.404124
> >> 4.488 0.398131
> >> 4.49 0.400173
> >> 4.492 0.407882
> >> 4.494 0.397001
> >> 4.496 0.396174
> >> 4.498 0.397463
> >> 4.5 0.397709
> >> 4.502 0.390852
> >> 4.504 0.386567
> >> 4.506 0.38765
> >> 4.508 0.390087
> >> 4.51 0.382283
> >> 4.512 0.375047
> >> 4.514 0.377796
> >> 4.516 0.38053
> >> 4.518 0.377125
> >> 4.52 0.376932
> >> 4.522 0.37867
> >> 4.524 0.369377
> >> 4.526 0.365101
> >> 4.528 0.372542
> >> 4.53 0.370364
> >> 4.532 0.366127
> >> 4.534 0.366403
> >> 4.536 0.363976
> >> 4.538 0.369198
> >> 4.54 0.357528
> >> 4.542 0.364936
> >> 4.544 0.356557
> >> 4.546 0.354742
> >> 4.548 0.356217
> >> 4.55 0.352334
> >> 4.552 0.356217
> >> 4.554 0.349794
> >> 4.556 0.354308
> >> 4.558 0.352529
> >> 4.56 0.349079
> >> 4.562 0.346863
> >> 4.564 0.348607
> >> 4.566 0.344551
> >> 4.568 0.336456
> >> 4.57 0.344059
> >> 4.572 0.342846
> >> 4.574 0.343857
> >> 4.576 0.342621
> >> 4.578 0.337633
> >> 4.58 0.332792
> >> 4.582 0.330428
> >> 4.584 0.335745
> >> 4.586 0.331965
> >> 4.588 0.332565
> >> 4.59 0.327592
> >> 4.592 0.324987
> >> 4.594 0.325261
> >> 4.596 0.323394
> >> 4.598 0.321757
> >> 4.6 0.320289
> >> 4.602 0.324011
> >> 4.604 0.317495
> >> 4.606 0.319709
> >> 4.608 0.31597
> >> 4.61 0.309665
> >> 4.612 0.311883
> >> 4.614 0.310811
> >> 4.616 0.309345
> >> 4.618 0.309076
> >> 4.62 0.306849
> >> 4.622 0.306059
> >> 4.624 0.302499
> >> 4.626 0.305746
> >> 4.628 0.302484
> >> 4.63 0.299231
> >> 4.632 0.300029
> >> 4.634 0.298585
> >> 4.636 0.299477
> >> 4.638 0.299858
> >> 4.64 0.296264
> >> 4.642 0.294672
> >> 4.644 0.293695
> >> 4.646 0.289712
> >> 4.648 0.290992
> >> 4.65 0.290417
> >> 4.652 0.288514
> >> 4.654 0.283104
> >> 4.656 0.283552
> >> 4.658 0.28347
> >> 4.66 0.290235
> >> 4.662 0.285285
> >> 4.664 0.277959
> >> 4.666 0.279521
> >> 4.668 0.273256
> >> 4.67 0.277063
> >> 4.672 0.276483
> >> 4.674 0.268187
> >> 4.676 0.271514
> >> 4.678 0.268644
> >> 4.68 0.271048
> >> 4.682 0.26606
> >> 4.684 0.275102
> >> 4.686 0.267121
> >> 4.688 0.266818
> >> 4.69 0.26788
> >> 4.692 0.264629
> >> 4.694 0.264184
> >> 4.696 0.267243
> >> 4.698 0.259231
> >> 4.7 0.262627
> >> 4.702 0.259008
> >> 4.704 0.256346
> >> 4.706 0.25399
> >> 4.708 0.255914
> >> 4.71 0.255821
> >> 4.712 0.254408
> >> 4.714 0.250453
> >> 4.716 0.24994
> >> 4.718 0.252659
> >> 4.72 0.248023
> >> 4.722 0.250594
> >> 4.724 0.24799
> >> 4.726 0.247123
> >> 4.728 0.24647
> >> 4.73 0.240303
> >> 4.732 0.242484
> >> 4.734 0.238625
> >> 4.736 0.239147
> >> 4.738 0.233552
> >> 4.74 0.238918
> >> 4.742 0.236821
> >> 4.744 0.234888
> >> 4.746 0.236861
> >> 4.748 0.233715
> >> 4.75 0.232752
> >> 4.752 0.23089
> >> 4.754 0.229704
> >> 4.756 0.231364
> >> 4.758 0.230899
> >> 4.76 0.231003
> >> 4.762 0.22809
> >> 4.764 0.228475
> >> 4.766 0.221141
> >> 4.768 0.228355
> >> 4.77 0.227719
> >> 4.772 0.221959
> >> 4.774 0.220529
> >> 4.776 0.217598
> >> 4.778 0.221624
> >> 4.78 0.219706
> >> 4.782 0.219037
> >> 4.784 0.221091
> >> 4.786 0.215576
> >> 4.788 0.214745
> >> 4.79 0.212482
> >> 4.792 0.214316
> >> 4.794 0.213902
> >> 4.796 0.215149
> >> 4.798 0.207367
> >> 4.8 0.212391
> >> 4.802 0.212657
> >> 4.804 0.211209
> >> 4.806 0.206088
> >> 4.808 0.205848
> >> 4.81 0.20539
> >> 4.812 0.20463
> >> 4.814 0.205654
> >> 4.816 0.201655
> >> 4.818 0.199889
> >> 4.82 0.203828
> >> 4.822 0.200193
> >> 4.824 0.203232
> >> 4.826 0.197828
> >> 4.828 0.197051
> >> 4.83 0.195867
> >> 4.832 0.197971
> >> 4.834 0.192717
> >> 4.836 0.194464
> >> 4.838 0.195964
> >> 4.84 0.193304
> >> 4.842 0.190948
> >> 4.844 0.19491
> >> 4.846 0.192267
> >> 4.848 0.193662
> >> 4.85 0.185335
> >> 4.852 0.187522
> >> 4.854 0.186923
> >> 4.856 0.185904
> >> 4.858 0.184322
> >> 4.86 0.181996
> >> 4.862 0.185581
> >> 4.864 0.18334
> >> 4.866 0.179438
> >> 4.868 0.182014
> >> 4.87 0.179545
> >> 4.872 0.180184
> >> 4.874 0.179702
> >> 4.876 0.178584
> >> 4.878 0.177543
> >> 4.88 0.179362
> >> 4.882 0.175853
> >> 4.884 0.176353
> >> 4.886 0.172697
> >> 4.888 0.173782
> >> 4.89 0.176536
> >> 4.892 0.171556
> >> 4.894 0.170033
> >> 4.896 0.170748
> >> 4.898 0.167734
> >> 4.9 0.170794
> >> 4.902 0.169994
> >> 4.904 0.169182
> >> 4.906 0.16627
> >> 4.908 0.166794
> >> 4.91 0.167013
> >> 4.912 0.165671
> >> 4.914 0.166407
> >> 4.916 0.164396
> >> 4.918 0.167214
> >> 4.92 0.164345
> >> 4.922 0.166479
> >>
> >>
> >> It terminated with value not 1 but close to zero.... any idea on
> this
> >>
> >> g_rdf -f PAP_8_Box12_HEP_EM_NPTNHPR_20ns_PBC.xtc -s
> >> PAP_8_Box12_HEP_EM_NPTNHPR_20ns.tpr -b 19000 -e 20000 -o RDFPBC.xvg
> >>
> >> Try to correct for PBC and still the same....
> >>
> >
> > Then my guess is you don't have a PBC problem, but perhaps a
> > convergence problem. You're only analyzing 1 ns of data; I don't
> > know what implication that will have for your outcome.
> >
> > -Justin
> >
> >> Regards,
> >>
> >> Rob
> >>
> >>
> >>
> >>
> >> Quoting teklebrh at ualberta.ca:
> >>
> >>>
> >>>
> >>> Dear Justine,
> >>>
> >>> No, I have 20 surfactants in my system in a pure solvent and try
> >>> to see the aggregation and adsorption properties of these
> >>> molecules. I have run them for 20ns. And now start to analyze my
> >>> data.
> >>>
> >>> Just in short I want to correct for periodicity! so that to get
> >>> the accurate result for my g_dist and g_rdf.
> >>>
> >>> trjconv on the xxx.xtc.... Like
> >>>
> >>> trjconv -f x_input.xtc -o x_out.xtc -b 19000 -e 20000 -s
> >>> x_input.tpr -pbc nojump
> >>>
> >>> or I have to cluster it first.
> >>>
> >>> Rob
> >>>
> >>> Quoting "Justin A. Lemkul" <jalemkul at vt.edu>:
> >>>
> >>>>
> >>>>
> >>>> teklebrh at ualberta.ca wrote:
> >>>>> Dear Justine,
> >>>>>
> >>>>> I just checked at the archive list and found out lots of
> information
> >>>>> on how to correct for PBC! And if I understood it correctly is
> this
> >>>>> the right way to follow...
> >>>>> "If you are trying to correct for periodicity for all species in
> >>>>> the system
> >>>>> (protein and water) then a few iterations of trjconv may be
> >>>>> necessary, i.e.
> >>>>> -center (on protein), followed by -pbc nojump or -pbc mol -ur
> >>>>> compact. Using
> >>>>> trjconv is a bit hit-or-miss, and just requires a bit of playing
> to
> >>>>> get things working how you want them."
> >>>>>
> >>>>> trjconv -f xxx.xtc -o x_cluster.gro -b 19000 -e 20000 -pbc
> cluster
> >>>>>
> >>>>> grompp -f xx.mdp -c x_cluster.gro -o x_cluster.tpr
> >>>>>
> >>>>> trjconv -f xxx.xtc -o x_cluster.xtc -b 19000 -e 20000 -s
> >>>>> xx_cluster.tpr -pbc nojump
> >>>>>
> >>>>> Then I have to do anaysis on the new .tpr and .xtc...... is that
> >>>>> correct....
> >>>>>
> >>>>> Is this the right way to correct for periodicity... I took the
> idea
> >>>>> from you and Chris. I have done the RDF and distance measurement
> but
> >>>>> looks a bit off my RDF did not converge to 1. They recommend me
> to
> >>>>> do PBC correction on my system which is run for 20ns. I am doing
> the
> >>>>> analysis on the last 1ns (from 19-20 ns)
> >>>>>
> >>>>
> >>>> Without context, I have no idea if what you've done is right. The
> >>>> clustering algorithm is only really useful for the formation of
> >>>> clustered molecules (i.e. micelles), so if you're trying to apply
> it
> >>>> to something else, I don't know that it will work.
> >>>>
> >>>> -Justin
> >>>>
> >>>>> You suggestion is always helpful!
> >>>>>
> >>>>> Rob
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>>
> >>>>> Quoting "Justin A. Lemkul" <jalemkul at vt.edu>:
> >>>>>
> >>>>>>
> >>>>>>
> >>>>>> teklebrh at ualberta.ca wrote:
> >>>>>>> Dear Gromacs,
> >>>>>>>
> >>>>>>> I want to correct the periodic boundry condition before
> analyzing
> >>>>>>> my data, how do I perform that?
> >>>>>>>
> >>>>>>
> >>>>>> Please see trjconv -h, as well as any of the thousands of posts
> in
> >>>>>> the list archive related to this topic.
> >>>>>>
> >>>>>> -Justin
> >>>>>>
> >>>>>>> any suggest!
> >>>>>>>
> >>>>>>> Rob
> >>>>>>
> >>>>>> --
> >>>>>> ========================================
> >>>>>>
> >>>>>> Justin A. Lemkul
> >>>>>> Ph.D. Candidate
> >>>>>> ICTAS Doctoral Scholar
> >>>>>> MILES-IGERT Trainee
> >>>>>> Department of Biochemistry
> >>>>>> Virginia Tech
> >>>>>> Blacksburg, VA
> >>>>>> jalemkul[at]vt.edu | (540) 231-9080
> >>>>>> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin[1]
> >>>>>>
> >>>>>> ========================================
> >>>>>> --
> >>>>>> gmx-users mailing list gmx-users at gromacs.org
> >>>>>> http://lists.gromacs.org/mailman/listinfo/gmx-users[2]
> >>>>>> Please search the archive at
> >>>>>> http://www.gromacs.org/Support/Mailing_Lists/Search[3] before
> posting!
> >>>>>> Please don't post (un)subscribe requests to the list. Use the
> www
> >>>>>> interface or send it to gmx-users-request at gromacs.org.
> >>>>>> Can't post? Read
http://www.gromacs.org/Support/Mailing_Lists[4]
> >>>>>>
> >>>>>>
> >>>>>
> >>>>
> >>>> --
> >>>> ========================================
> >>>>
> >>>> Justin A. Lemkul
> >>>> Ph.D. Candidate
> >>>> ICTAS Doctoral Scholar
> >>>> MILES-IGERT Trainee
> >>>> Department of Biochemistry
> >>>> Virginia Tech
> >>>> Blacksburg, VA
> >>>> jalemkul[at]vt.edu | (540) 231-9080
> >>>> http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin[5]
> >>>>
> >>>> ========================================
> >>>> --
> >>>> gmx-users mailing list gmx-users at gromacs.org
> >>>> http://lists.gromacs.org/mailman/listinfo/gmx-users[6]
> >>>> Please search the archive at
> >>>> http://www.gromacs.org/Support/Mailing_Lists/Search[7] before
> posting!
> >>>> Please don't post (un)subscribe requests to the list. Use the www
> >>>> interface or send it to gmx-users-request at gromacs.org.
> >>>> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists[8]
> >>>>
> >>>>
> >>>
> >>>
> >>>
> >>> Links:
> >>> ------
> >>> [1] http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
> >>> [2] http://lists.gromacs.org/mailman/listinfo/gmx-users
> >>> [3] http://www.gromacs.org/Support/Mailing_Lists/Search
> >>> [4] http://www.gromacs.org/Support/Mailing_Lists
> >>> [5] http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
> >>> [6] http://lists.gromacs.org/mailman/listinfo/gmx-users
> >>> [7] http://www.gromacs.org/Support/Mailing_Lists/Search
> >>> [8] http://www.gromacs.org/Support/Mailing_Lists
> >>
> >
> > --
> > ========================================
> >
> > Justin A. Lemkul
> > Ph.D. Candidate
> > ICTAS Doctoral Scholar
> > MILES-IGERT Trainee
> > Department of Biochemistry
> > Virginia Tech
> > Blacksburg, VA
> > jalemkul[at]vt.edu | (540) 231-9080
> > http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
> >
> > ========================================
> > --
> > gmx-users mailing list gmx-users at gromacs.org
> > http://lists.gromacs.org/mailman/listinfo/gmx-users
> > Please search the archive at
> > http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> > Please don't post (un)subscribe requests to the list. Use the www
> > interface or send it to gmx-users-request at gromacs.org.
> > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> >
> >
>
> --
> gmx-users mailing list gmx-users at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-users
> Please search the archive at
> http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
> Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
More information about the gromacs.org_gmx-users
mailing list