[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