Hi everyone, I have a receptor-ligand system and from umbrella sampling at different receptor-ligand distances I could find the distances at which the ligand is most probably going to bind. So, in order to get an idea of time taken for binding, I am running 10-independent unbiased MD simulations each starting from same initial structure (corresponding to a specific receptor-ligand distance where it is most likely to bind); so in order to perform 10-independent simulations I modify the velocity generation random seeds i.e. gen-seed in mdp files, e.g. if I have trajectories like: seed_1,seed_2,…,seed_10 then I put gen-seed values= 1,2,…,10 respectively for seed_1,seed_2,…,seed_10; Is this a correct way of running these or the manual choice of random seeds that differ by only one can impact the randomness of the simulations? I am of the view that using different random seeds will always give us different starting velocities. Any help will be much appreciated, thank you

Seeds values that are consecutive integers will lead to essentially identical starting conditions. I tested this a long time ago. You need to assign dramatically different values to get proper randomness.

Thank you for your reply Dr. Lemkul. In my case the starting structure is same for the independent simulations; in order to make them random I would like to assign different sets of velocities to the initial structure; So, if consecutive integers are used as gen-seed then the velocities of the starting structure in different simulations will be very similar to each other, I hope this is what you meant by identical starting conditions; I have done a small test: generated 1000 random numbers in cpp twice using random-seed=1 and 2 respectively and then measured and found they are very much uncorrelated; so, does gromacs work in a different way? and what is then the optimal process to run multiple copies of a simulation starting from same structure (using gen-seed=-1)?

gen_seed = -1 is generally a safe way to get suitable random numbers. You can always verify what the seed was in mdout.mdp or by printing the contents of the .tpr file with gmx dump.