From: winardi, erik (
Date: Tue Aug 30 2011 - 10:44:49 CDT

Dear Dr. Kohlmeyer,  Thank you for your respond, it really helpful.  I will explain what I am currently doing.  I ran propene oxidation  containt 1 propene and 100 O2 = 209 total atoms at Temperature of 2500K using NVT ensemble. I am using reaxFF to do the simulation and come up with several molecules like CO and H2O2.  I am trying to get the CO(2 atoms) only and H2O2(4 atoms) only from the molecule data result without taking all the other atoms.  What I did is, I picked the atoms in VMD and wrote down all the atom index for those atom molecule and make a new lammps input file manually for these 2 molecules (CO and H2O2). We want to try to run and break the molecule using reaxff and measure the energy for that both molecules.  I couldn't create the new input file for both atoms. I just wondering how to write the correct atom input for lammps so I could write a converter for write the input data for lammps using the information I got from the VMD by picking the atom (Coordinate, Atom Index, etc).  What I can think of is to delete all the molecules in the box and leave the molecule of interest, but it will take a lot of time to pick more about 200 atoms in VMD visualization panel just to delete the atoms. Is there any better way to do this ? I hope I explain it clear enough. Thank so much,   Erik    ________________________________ From: Axel Kohlmeyer <> To: "winardi, erik" <> Cc: "" <> Sent: Monday, August 29, 2011 1:20 PM Subject: Re: vmd-l: Lammps Output data format On Mon, 2011-08-29 at 10:20 -0700, winardi, erik wrote: > Dear VMD users, dear erik, > I have questions about using the lammps.trj output from VMD. > 1. I tried to use the lammps.trj format datafile for lammps data > input, I am guessing that I can't use the output of lammps.trj for > lammps input file. I understand that VMD can use the topo > writelammpsdata command to provide the initial structure data to > lammps. I just curious when do we need to provide this lammps.trj data > in lammps? i don't really understand your question, which seems to be more appropriate for the lammps-users mailing list in the first place. do you want to use a lammps trajectory dump as input for the read_data command? this is not possible. > 2. I went trough the topotools website but I couldn't find how the > data format for topo writelammpsdata is written. I am trying to make a > converter for different data format to be converted to lammps data > format. It would be great if I could get the format of writting the > initial lammps molecule data. again, this question would be better asked on the lammps-users mailing list. the format of the lammps data file is documented in great detail in the lammps manual. look up the documentation of the read_data command. what kind of a conversion are you looking for. for the most part the topotools plugin aims to provide already the infrastructure to write correctly formatted lammps data files and VMD provides (through molfile/babel) access to a large variety of input file formats. the tricky part in between is to generate/reconstruct the pieces of information that are needed for a proper topology (i.e. lammps data) file, that are missing from the many existing coordinate file formats. this is something that is very system specific and thus it is impossible to write a simple converted. the (incomplete) topotools tutorial on my home page gives a few hints at which steps are missing and how to generate the information for a couple of rather simple systems. if you can give me a more specific example of the kind of problem that you want to address, i'd be happy to make further suggestions or even provide some script templates, if they are easy enough to do. i would also encourage you to reuse existing tools as much as possible instead of writing your own tool from scratch, unless there are very, very good reasons to do so. thanks,     axel. > > Many Thanks, > Erik -- Dr. Axel Kohlmeyer Institute for Computational Molecular Science Temple University, Philadelphia PA, USA.