# Topology file for metal ion **URL:** \ **Category:** User discussions\ **Created:** [February 4, 2024, 2:50pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218 "2024-02-04T14:50:03Z")\ **Posts on this page:** 9\ **Page:** 1
**Author:** ![Nightwolf](https://avatars.discourse-cdn.com/v4/letter/n/779978/32.png) [@Nightwolf](https://gromacs.bioexcel.eu/u/Nightwolf)\ **Post date:** [February 4, 2024, 2:50pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/1 "2024-02-04T14:50:04Z")
GROMACS version: 2023 GROMACS modification: No Hello everyone, just asking is it possible to make a topology file for a metal ion like Na+? ---
**Author:** ![milosz.wieczor](https://avatars.discourse-cdn.com/v4/letter/m/b19c9b/32.png) [@milosz.wieczor](https://gromacs.bioexcel.eu/u/milosz.wieczor)\ **Post date:** [February 4, 2024, 6:00pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/2 "2024-02-04T18:00:03Z")
Sure, here it is, fully self-contained: > [defaults] > 1 2 yes 0.5 0.8333 > > [atomtypes] > Na 11 22.99 0.0000 A 0.243934000000 0.36585 > > [moleculetype] > Na 1 > > [atoms] > 1 Na 1 Na Na 1 1.0 22.9898 > > [system] > sodium > > [molecules] > Na 1 But I guess you want to add it to an existing system? Or have Na+ mixed with something else? ---
**Author:** ![Nightwolf](https://avatars.discourse-cdn.com/v4/letter/n/779978/32.png) [@Nightwolf](https://gromacs.bioexcel.eu/u/Nightwolf)\ **Post date:** [February 4, 2024, 6:13pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/3 "2024-02-04T18:13:01Z")
Thanks. Yes I want to attach it as a ligand to an existing molecule making it one. ---
**Author:** ![milosz.wieczor](https://avatars.discourse-cdn.com/v4/letter/m/b19c9b/32.png) [@milosz.wieczor](https://gromacs.bioexcel.eu/u/milosz.wieczor)\ **Post date:** [February 4, 2024, 6:28pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/4 "2024-02-04T18:28:51Z")
Most likely you don’t want covalent bonds with sodium though, right? Except maybe for some extreme organosodium compounds, it’s going to be a non-covalent sodium-binding site, and then it’s enough to have a sodium ion defined in the topology. Normally `ions.itp` is automatically embedded in the topology, so you can just add `Na 1` to the `[molecules]` section to account for that; the rest will rely on a correct placement of the ion in the structure file. If you indeed mean a covalent bond, then it gets more complicated. ---
**Author:** ![Nightwolf](https://avatars.discourse-cdn.com/v4/letter/n/779978/32.png) [@Nightwolf](https://gromacs.bioexcel.eu/u/Nightwolf)\ **Post date:** [February 4, 2024, 6:35pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/5 "2024-02-04T18:35:57Z")
No its not a covalent bond but rather an ionic one. I just want to couple the ion with the molecule and simulate them together as a whole. ---
**Author:** ![milosz.wieczor](https://avatars.discourse-cdn.com/v4/letter/m/b19c9b/32.png) [@milosz.wieczor](https://gromacs.bioexcel.eu/u/milosz.wieczor)\ **Post date:** [February 4, 2024, 11:08pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/6 "2024-02-04T23:08:16Z")
Then, as I said, if the interaction is strong enough, it should be sufficient to have them both in a bound state at the beginning of the simulation, without any additional bonds specified in the topology. Of course, interactions with ions are never perfectly modeled in additive (non-polarizable) force fields so it’s an open question whether a naive approach will work for you; if it’s somehow experimentally known that the two should stay together and in the simulation they don’t, one can try out different sodium parameter sets, or simply add a restraint. This very much depends on the ultimate purpose of your simulation. ---
**Author:** ![Nightwolf](https://avatars.discourse-cdn.com/v4/letter/n/779978/32.png) [@Nightwolf](https://gromacs.bioexcel.eu/u/Nightwolf)\ **Post date:** [February 5, 2024, 10:16am UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/7 "2024-02-05T10:16:17Z")
The system am simulating is supposed to contain a partial dissociation between the molecule and the Na counterion. i.e some molecules will be dissociated and others will have the sodium attached the hypothesis is that this will answer some anomalies being recorded for some macroscopic property of the system. ---
**Author:** ![milosz.wieczor](https://avatars.discourse-cdn.com/v4/letter/m/b19c9b/32.png) [@milosz.wieczor](https://gromacs.bioexcel.eu/u/milosz.wieczor)\ **Post date:** [February 5, 2024, 12:06pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/8 "2024-02-05T12:06:08Z")
Right, then in theory you’re good to go with a plain non-bonded sodium (that is, no extra additions), just make sure that you’re using a sodium parametrization that was validated with a similar type of interaction - could be papers on crown ethers etc. There are many choices around, but most of them focus on reproducing co-solvent properties rather than binding to ligands. ---
**Author:** ![Nightwolf](https://avatars.discourse-cdn.com/v4/letter/n/779978/32.png) [@Nightwolf](https://gromacs.bioexcel.eu/u/Nightwolf)\ **Post date:** [February 5, 2024, 1:13pm UTC](https://gromacs.bioexcel.eu/t/topology-file-for-metal-ion/8218/9 "2024-02-05T13:13:44Z")
Ok, thanks. Really appreciate the assistance