[gmx-users] trjconv, -pbc whole, and Max OS X?
Emma Falck
efalck at cc.hut.fi
Fri Nov 18 17:12:56 CET 2005
I just installed gromacs 3.2.1 to my system. Here trjconv with the option
-pbc whole works exactly as I would expect it to, ie, leaves the molecules
whole/intact. I am using exactly the same syntax as I used with the 3.3
version.
Emma
On Fri, 18 Nov 2005, Emma Falck wrote:
> Date: Fri, 18 Nov 2005 00:45:34 +0200 (EET)
> From: Emma Falck <efalck at cc.hut.fi>
> Reply-To: Discussion list for GROMACS users <gmx-users at gromacs.org>
> To: Discussion list for GROMACS users <gmx-users at gromacs.org>
> Subject: Re: [gmx-users] trjconv, -pbc whole, and Max OS X?
>
> Erik,
>
> I was wrong about the first frame. When I look at my .trr file
> using ngmx, then all lipids in all frames are whole, as they should. It is
> only when I look at my .gro file -- either with ngmx or vmd -- when I see
> broken lipids in most of the frames. This .gro file I have generated with
> trjconv using the -pbc whole option. This does not seem to restrict itself to
> Mac OS X -- I checked with RedHat, and ran into the same problem.
>
> Emma
>
> On Thu, 17 Nov 2005, Erik Lindahl wrote:
>
>> Date: Thu, 17 Nov 2005 23:22:26 +0100
>> From: Erik Lindahl <lindahl at sbc.su.se>
>> Reply-To: Discussion list for GROMACS users <gmx-users at gromacs.org>
>> To: Discussion list for GROMACS users <gmx-users at gromacs.org>
>> Subject: Re: [gmx-users] trjconv, -pbc whole, and Max OS X?
>>
>> Emma,
>>
>> Are you sure your individual molecules (single lipids) are broken into
>> parts?
>>
>> mdrun never writes broken molecules, but the periodic image is determined
>> so the first atom of each molecule is inside the box and no periodicity
>> jumps occur inside a molecule.
>>
>> There are a couple of other -pbc options in case you want to force all
>> atoms inside the box (break molecules) or to keep a membrane more intact.
>>
>> Cheers,
>>
>> Erik
>>
>> On Nov 17, 2005, at 10:58 PM, Emma Falck wrote:
>>
>>>
>>> The first frame seems to contain broken lipids, too. I thought trjconv
>>> would make the lipids whole irrespective of this, but apparently this is
>>> not the case. What should I do?
>>>
>>>> What about the initial structure of the run?
>>>
>>>> On 11/17/05, Emma Falck <efalck at cc.hut.fi> wrote:
>>>>>
>>>>> Hi,
>>>>>
>>>>> My problem is with the trjconv command. I am trying to convert a
>>>>> trajectory
>>>>> in the .trr format to a trajectory in the .gro format, such that the
>>>>> molecules are whole in the resulting trajectory. To do so, I use trjconv
>>>>> with the option '-pbc whole'. Even so, my molecules -- both lipids and
>>>>> water molecules -- are broken in the resulting .gro file.
>>>>>
>>>>> I recently installed Gromacs 3.3 to my PowerBook, that runs OS X 10.4.
>>>>> I did not compile Gromacs from the source, but used the OS X package
>>>>> provided
>>>>> on the gromacs pages.
>>>>>
>>>>> Has anyone else had related/similar problems with trjconv, in
>>> particular
>>>>> with
>>>>> OS X?
>>>>>
>>>>> Best regards,
>>>>>
>>>>> Emma
>>>>>
>>>>> --
>>>>> Dr. Emma Falck
>>>>> Beckman Institute, UIUC
>>>
>>>> --
>>>> Viswanadham Sridhara,
>>>> Graduate Research Assistant,
>>>> Old Dominion University, "VIRGINIA".
>>> _______________________________________________
>>> gmx-users mailing list
>>> gmx-users at gromacs.org
>>> http://www.gromacs.org/mailman/listinfo/gmx-users
>>> Please don't post (un)subscribe requests to the list. Use the www
>>> interface or send it to gmx-users-request at gromacs.org.
>>>
>>
>> -----------------------------------------------------------
>> Erik Lindahl <lindahl at sbc.su.se>
>> Assistant Professor, Stockholm Bioinformatics Center
>> Stockholm University, SE 106 91 Stockholm
>> Phone: +46 8 5537 8564 Fax: +46 8 5537 8214
>>
>>
>>
>>
>> _______________________________________________
>> gmx-users mailing list
>> gmx-users at gromacs.org
>> http://www.gromacs.org/mailman/listinfo/gmx-users
>> Please don't post (un)subscribe requests to the list. Use thewww interface
>> or send it to gmx-users-request at gromacs.org.
>>
> _______________________________________________
> gmx-users mailing list
> gmx-users at gromacs.org
> http://www.gromacs.org/mailman/listinfo/gmx-users
> Please don't post (un)subscribe requests to the list. Use the www interface
> or send it to gmx-users-request at gromacs.org.
>
More information about the gromacs.org_gmx-users
mailing list