![]() ![]() It is meant to be a complementary tool to psfgen, which is very much optimized for building topologies for biomolecules. OPLS/AA has an “increment system” where those can be computed based on the specifics of the atom type assignment, which must be done by a human (although it can be scripted) for each particular case.TopoTools is a VMD plugin for manipulating topology information. xyz files and there is no simple way to recover that information heuristically. for example, the (partial) charge information that is crucial in many molecular force fields is missing in PDB or. of course, heuristics are dependent on the choice of force field. and sometimes, heuristic topology generation will give wrong results. there are add-on functions to apply certain heuristics to guess missing data, but the correctness of that is depending on the information provided. that cannot be, since most coordinate files (like PDB files, which are a very, very bad choice, BTW) are missing a lot of information, TopoTools (or rather VMD that it relies on) can only use the information available. The fact, that one can write scripts on top of TopoTools to, e.g., read or write LAMMPS data file (but not only those) is an additional feature, but not its main purpose and for certain, it is not meant to be a fully automatic, do-the-right-thing converted from coordinate files to topology files. that is ok under many circumstances, but it can become a problem, when trying to do something non-standard. indeed, one of the main motivations is to be different from other tools that often make certain assumption and make choices for the user. ![]() ![]() It very deliberately (like VMD) is force field agnostic, and it depends on users knowing what they are doing. Please also keep in mind, that TopoTools was not explicitly written for LAMMPS but for use with a variety of MD software packages. preparing a data file for use with ReaxFF using TopoTools is about the easiest conversion, since you don’t need to do any atom typing (just assign elements) and you don’t need any topology information (just write out the data file in atom style charge). So using this default forcefield to preapare the system and again using different force field (reaxff) to simulate the structure is not the wrong way to do? So I think I can not prepare ZnO lammps data file other than defining explicitly inside lammps input script. The VMD TopoTools plugin ignores the force field and only inserts commented out dummy entries. But here my doubt is when I am converting it to lammps data file topotools takes one of the default force field (which I dont know) which is not obviously reaxff. I have made a ZnO monolayer using vesta then converted it to lammps data file using vmd topo tools. I have a confusion of simulating ZnO monolayer using reaxff.
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |