[gmx-users] GROMACS 4.5.4 keep crashing all the time

chris.neale at utoronto.ca chris.neale at utoronto.ca
Thu Sep 1 04:20:01 CEST 2011


Itamar:

We really are trying to help. I think that perhaps you don't grasp how  
difficult it is to help without being able to access the simulation  
directly. Therefore we have ideas and we ask you to do specific things  
that are going to move us toward a solution, either by finding answers  
or by ruling out possibilities.

It is actually useful information to know that "Sometimes it is the  
peptide N and H, like in the case of 981 and 982, and sometimes  
others" ... but when it seems like you don't want to provide the  
requested information my first inclination is to give up on trying to  
help.

At this point, there are a few unanswered old questions and I have  
some new questions.

1. Can you reproduce this with a water box?
2. Can you reproduce this with your protein in vacuum?
3. If neither 2 or 3, then can you step slowly from one of these  
systems toward your final system and identify the point at which the  
lincs warnings arise?
4. Do you get the warnings without Ca also, or just with Ca?
5. Can you reproduce this with the SD integrator? If you are really  
against trying this, then at least can you reproduce this with a  
single Berendsen temperature coupling group?
6. Can you reproduce this without using the reaction field? Either  
with PME or a simple cutoff?
7. Can you trace down the version of gromacs (between 4.0.7 and 4.5.4)  
where you start to see this warning?

Chris.

-- original message --

Hi Justin,

I did repeat it using gen_val and running temperature the same, with  
no effect, it is still crash. I didn't replied point #6 because the  
atoms which triggers the LINCS are different between each try.  
Sometimes it is the peptide N and H, like in the case of 981 and 982,  
and sometimes others. In addition, there is no visible difference in  
dynamics between 4.0.7 and 4.5.4 I can find.

Itamar


On 01/09/2011, at 11:14 AM, Justin A. Lemkul wrote:

>
>
> Itamar Kass wrote:
>> Hi Mark,
>> I didn't had the time to do the SD yet, but serial run end with the  
>> same results. I didn't try water only system, as this is of no  
>> interest to me, but I will simplify the system later on.
>
> Being of interest to you and being a useful diagnostic may be  
> different.  It's important to rule out different variables to arrive  
> at a solution, which I suspect is of interest to you.  You also  
> haven't addressed points #1 and #6 in Chris' message.
>
> -Justin
>
>> Cheers,
>> Itamar
>> On 01/09/2011, at 10:51 AM, Mark Abraham wrote:
>>> On 1/09/2011 10:20 AM, Itamar Kass wrote:
>>>> Hi Chris,
>>>>
>>>> Thanks for the email, I am sorry it took me some time to replay.  
>>>> I tried 4.5.4 again, now starting from a 5 ns simulations run  
>>>> using 4.0.7, and again the simulations had stopped after 1000  
>>>> LINCS error (I can extend the simulations using 4.0.7).
>>>>
>>>> I know that gromacs stopped after 1000 LINCS, but this is usually  
>>>> a sign that something bad is going on in the system.
>>> OK. Chris suggested a number of other strategies that will help  
>>> determine which aspect of 4.5.4 is behaving differently. How did  
>>> those strategies work out?
>>>
>>> Mark
>>>
>>>> Cheers,
>>>> Itamar
>>>>
>>>> On 18/08/2011, at 12:03 PM, chris.neale at utoronto.ca wrote:
>>>>
>>>>> OK, here's my last few ideas:
>>>>>
>>>>> 1. Please try to repeat this with gen_vel set to the same value  
>>>>> as your temperature coupling
>>>>>
>>>>> 2. Can you reproduce this in serial?
>>>>>
>>>>> 3. Can you reproduce this with the sd integrator?
>>>>>
>>>>> 4. Can you reproduce this with a simpler system? protein in  
>>>>> vacuum or just water or remove the ions, etc?
>>>>>
>>>>> 5. Take the output .gro from 4.0.7 that ran fine for X ns and  
>>>>> run it under 4.5.4. Do you get the same lincs warnings?
>>>>>
>>>>> 6. Also, note that you are getting warnings and the run does not  
>>>>> actually crash but just stops after too many warnings. So what  
>>>>> are atoms 981 and 982? Does their motion look different in an  
>>>>> important ways between the 4.0.7 and 4.5.4 trajectories?
>>>>>
>>>>> Chris.
>>>>>
>>>>> -- original message --
>>>>>
>>>>> Hi Chris,
>>>>>
>>>>> thanks for the advice, I have to say I tried this as well  
>>>>> without any success.
>>>>>
>>>>> Itamar
>>>>>
>>>>> On 18/08/2011, at 11:11 AM, chris.neale at utoronto.ca wrote:
>>>>>
>>>>>> run an EM with flexible water. I often find that this is the  
>>>>>> only way to get a stable system. 500 steps of steep with  
>>>>>> define=-DFLEXIBLE (or different depending on your water model I  
>>>>>> think) should be enough.
>>>>>>
>>>>>> Chris.
>>>>>>
>>>>>> Hi Chris and Justin,
>>>>>>
>>>>>> On 18/08/2011, at 9:36 AM, Justin A. Lemkul wrote:
>>>>>>
>>>>> --
>>>>> 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
>>>> -----
>>>> "In theory, there is no difference between theory and practice.  
>>>> But, in practice, there is." - Jan L.A. van de Snepscheut
>>>>
>>>> ===========================================
>>>> | Itamar Kass, Ph.D.
>>>> | Postdoctoral Research Fellow
>>>> |
>>>> | Department of Biochemistry and Molecular Biology
>>>> | Building 77 Clayton Campus
>>>> | Wellington Road
>>>> | Monash University,
>>>> | Victoria 3800
>>>> | Australia
>>>> |
>>>> | Tel: +61 3 9902 9376
>>>> | Fax: +61 3 9902 9500
>>>> | E-mail: Itamar.Kass at monash.edu
>>>> ============================================
>>>>
>>> -- 
>>> 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
>> -----
>> "In theory, there is no difference between theory and practice.  
>> But, in practice, there is." - Jan L.A. van de Snepscheut
>> ===========================================
>> | Itamar Kass, Ph.D.
>> | Postdoctoral Research Fellow
>> |
>> | Department of Biochemistry and Molecular Biology
>> | Building 77 Clayton Campus
>> | Wellington Road
>> | Monash University,
>> | Victoria 3800
>> | Australia
>> |
>> | Tel: +61 3 9902 9376
>> | Fax: +61 3 9902 9500
>> | E-mail: Itamar.Kass at monash.edu
>> ============================================
>
> -- 
> ========================================
>
> 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

-----
"In theory, there is no difference between theory and practice. But,  
in practice, there is." - Jan L.A. van de Snepscheut

===========================================
| Itamar Kass, Ph.D.
| Postdoctoral Research Fellow
|
| Department of Biochemistry and Molecular Biology
| Building 77 Clayton Campus
| Wellington Road
| Monash University,
| Victoria 3800
| Australia
|
| Tel: +61 3 9902 9376
| Fax: +61 3 9902 9500
| E-mail: Itamar.Kass at monash.edu





More information about the gromacs.org_gmx-users mailing list