# Properly including bromide in GROMOS96 54a7
**URL:** https://gromacs.bioexcel.eu/t/properly-including-bromide-in-gromos96-54a7/4836
**Category:** User discussions
**Tags:** forcefield
**Created:** [September 22, 2022, 1:38pm UTC](https://gromacs.bioexcel.eu/t/properly-including-bromide-in-gromos96-54a7/4836 "2022-09-22T13:38:25Z")
**Posts on this page:** 3
**Page:** 1
### Author:  [@esmee](https://gromacs.bioexcel.eu/u/esmee)
#### Post date: [September 22, 2022, 1:38pm UTC](https://gromacs.bioexcel.eu/t/properly-including-bromide-in-gromos96-54a7/4836/1 "2022-09-22T13:38:25Z")
GROMACS version: 2021.3
GROMACS modification: Yes
Dear GROMACS users,
I am a beginner to GROMACS and have some previous work form a colleague available. For my purposes, I need to add bromide (Br-) to gromos 54a7 (I have a copy of the force field in my working directory).
The way my colleague has done this is by modifying _only_ the ions.itp-file with the following:
[moleculetype]
; molname nrexcl
BR 1
[atoms]
; id at type res nr residu name at name cg nr charge mass
1 BR 1 BR BR 1 -1 79.904
I thought the form of “at type” (i.e. BR) was a bit odd compared to the other standard ions that have some kind of charge in their entry for “at type” (e.g. CL-), so when I adapted this I put BR- instead to match the syntax, thinking the omitted minus-sign was a typo. No other files of the force field have been modified by my colleague.
When I then tried to run the energy minimization step (following the Lysozyme tutorial) i got the error message “atomtype BR- not found”.
I then changed to “at type” BR instead of BR-, as my colleague used this, which resolved the error. I am however not completely sure that using BR is the correct approach. In atomtypes.atp I can see that there is bromine (non-ionic) present, which is denoted BR. **Does this mean that when putting BR instead of BR- under “at type” in ions.itp, the parameters in e.g. ffnonbonded.itp for the neutral bromine are used also for the bromide**? If so, is this an ok approach anyway (I can for example see that values in ffnonbonded.itp differ between entries for CL and CL-, which is my reason for doubting the approach)?
If it is not an ok practice, how do I properly add BR-? Do I have to add it also to atomtypes.atp, and to add values for c6, c12 and nonbond\_params in ffnonbonded (where does one obtain such values for bromide, if so)? Are there any other files of the gromos54a7 force field that have to be modified as well?
Thank you in advance to anyone willing to sort out my confusion.
---
### Author:  [@alevilla](https://gromacs.bioexcel.eu/u/alevilla)
#### Post date: [September 30, 2022, 11:39am UTC](https://gromacs.bioexcel.eu/t/properly-including-bromide-in-gromos96-54a7/4836/3 "2022-09-30T11:39:23Z")
Hi,
> [@esmee](#):
>
> When I then tried to run the energy minimization step (following the Lysozyme tutorial) i got the error message “atomtype BR- not found”.
The reason is that you change the name from BR → BR- and BR- is not defined as atomtype in the force field files. Atom type information are storage in ffnonbonded.itp in the [atomtypes] directive and [nonbond\_params] directive (for gromos force field) and in the database atomtypes.atp.
To know more on the format see [https://manual.gromacs.org/current/reference-manual/topologies/topology-file-formats.html](https://manual.gromacs.org/current/reference-manual/topologies/topology-file-formats.html)
> [@esmee](#):
>
> In atomtypes.atp I can see that there is bromine (non-ionic) present, which is denoted BR. **Does this mean that when putting BR instead of BR- under “at type” in ions.itp, the parameters in e.g. ffnonbonded.itp for the neutral bromine are used also for the bromide**?
yes.
> [@esmee](#):
>
> If so, is this an ok approach anyway (I can for example see that values in ffnonbonded.itp differ between entries for CL and CL-, which is my reason for doubting the approach)?
BR and BR- are different, it is more appropriate to describe them using different LJ parameters.
> [@esmee](#):
>
> If it is not an ok practice, how do I properly add BR-? Do I have to add it also to atomtypes.atp, and to add values for c6, c12 and nonbond\_params in ffnonbonded (where does one obtain such values for bromide, if so)? Are there any other files of the gromos54a7 force field that have to be modified as well?
When you add a new ion, it is good to add it also in ions.itp.
\Alessandra
---
### Author:  [@esmee](https://gromacs.bioexcel.eu/u/esmee)
#### Post date: [September 30, 2022, 3:52pm UTC](https://gromacs.bioexcel.eu/t/properly-including-bromide-in-gromos96-54a7/4836/4 "2022-09-30T15:52:11Z")
Thank you so much for your detailed answer Alessandra!
So, if I understand you correctly I would have to do the following steps in order to properly add BR- and not use the LJ parameters for BR:
1. Add BR- with correct values for c6 and c12 in the [atomtypes] directive of ffnonbonded.itp.
2. Add BR- with correct values for c6 and c12 in the [nonbond\_params] directive of ffnonbonded.itp
3. Add BR- to the database atomtypes.atp
4. Add BR- to ions.itp
Step 3 and 4 seem perfectly straightforward to me, but I am not completely sure how to do step 1 and 2.
> [@alevilla](#):
>
> BR and BR- are different, it is more appropriate to describe them using different LJ parameters.
Where would one normally find the LJ parameters for an ion, in my case bromide, if one would like to add these to ffnonbonded.itp?