To whom could lend me some help,
A a beginner I am first creating the simulation box as a multi-component system. I downloaded the Force-Fields from the ATb server, for which I have successfully simulated several MD and BD runs. The procedure I have followed is:
- Generate LAMMPS input scripts using Moltemplate and compiling the “.lt” files extracted from ATb for the FF and the molecule topology.
- From the generated “data” file I manually produce the molecule “.txt” file to be compiled along with the “system.in.init” and “system.in.settings” files.
- I follow up the commands that I provide in the “run.in” file here.
I have successfully loaded the molecule I name as “trans-OMCA.txt” that has 22 atoms. However, when trying the same for the “PGAC.txt” molecule that holds 98 atoms, I get the following error:
ERROR: on proc 0: Unexpected line in molecule file while looking for the next section: 98 20
I have double checked that the “PGAC.txt” follows the same pattern as the “t-OMCA.txt” molecule succesfully compiled. I have already tried the following, but nothing works, but one try gave me a hint on what the problem is:
- Fix the topology data amounts while doing “create_box” so that this includes what is necessary for this “PGAC.txt” molecule, but any # seems to work.
- I wondered if the 98 atoms could be too much for the molecule command to be able to compile, thus I tried to fix 97 atoms instead of 98, and commented the Coords, Types, and Charges; for the leftover molecule #98, and the molecule command strangely compiled this modified “PGAC.txt” file, but without placing the Types for each atom. This results in not being able to used the tagged molecule to be inserted in the simulation box with the “create_atoms” command with which I previously inserted the “t-OMCA” molecules successfully.
The latter test I performed gave me the hint on the possibility that 98 toms may be too much for a molecule to hold. Thus, my fina try was to use a similar molecule, but with 62 atoms (CTAB.txt) to see if the amount of atoms is really affecting, or if there was something I was missing. The result was that this “CTAB.txt” molecule was compiled correctly. Being quite similar to the “PGAC.txt”, but with 62 atoms instead of 98, made me notice that I am not missing anything in terms if the necessary topology data included
However, I would really like to know if there is something I can do to insert replicas of this “PGAC.txt” molecule within my simulating system, given that this “molecule” command is limiting the atom count, and this “CTAB.txt” is not enough for my purposes.
Regards
jfernando
P.S. I tried to upload the files I mentioned before, but it seem that is not a feature available for new users. In case there is a reply from someone who requests the files as usual, and somehow the system enables me to upload them for not being a new user anymore, I will certainly do it.