Your most recent question sounds like a VMD question. I don’t know enough about VMD to answer it.
Fortunately, perhaps I don’t need to. The issue has nothing to do with PSF files. Neither LAMMPS, not moltemplate read or care about PSF files. You can run your simulations without them. (I have only ever needed PSF files to view my molecules in VMD.)
However if you are adding ions, then you will need to inform LAMMPS and moltemplate about them. The simplest way to do that is to modify your “system.lt” file this way:
---- “system.lt” file: ----
import “spce.lt”
import “ions.lt”
wat = new SPCE [260]
na = new NaIon [3]
cl = new ClIon [3]
---- end of “system.lt” file ----
(Again, “3” and “260” are the number of ions and water molecules in this example.)
— ions.lt —
For a working example of the “ions.lt” file, see this link:
http://moltemplate.org/examples/waterSPCE+Na+Cl/
…or if that link is broken, try this one:
https://github.com/jewettaij/moltemplate/tree/master/examples/all_atom/force_field_explicit_parameters/waterSPCE%2BNa%2BCl
The “ions.lt” file contains additional pair_coeff commands to define Lennard Jones parameters for these atoms. Later on, as your system gets more complicated, you may decide it is more convenient to use a canned force field that already contains definitions for all these atom types by replacing the “import ions.lt” line with one of these choices:
import “oplsaa.lt” #(OPLS boss4.8 force field)
-or-
import “gaff.lt” #(AMBER force field)
#(…and rename your @atom types accordingly)
Cheers
Andrew
More force fields are planned. Right now, there are two scripts distributed with moltemplate which have been written to add more force field support (tinkerparm2lt.py and emcprm2lt.py) these scripts can convert some of the other force fields distributed with TINKER and EMC into moltemplate format. But force field conversion can be arduous work, and some force fields contain terms which are not yet implemented in lammps.