# How to use Clay.ff with gromacs for molecular dynamics simulation of mica

**URL:** <https://gromacs.bioexcel.eu/t/how-to-use-clay-ff-with-gromacs-for-molecular-dynamics-simulation-of-mica/10881>\
**Category:** User discussions\
**Tags:** forcefield\
**Created:** [December 22, 2024, 2:15pm UTC](https://gromacs.bioexcel.eu/t/how-to-use-clay-ff-with-gromacs-for-molecular-dynamics-simulation-of-mica/10881 "2024-12-22T14:15:09Z")\
**Posts on this page:** 3\
**Page:** 1

<div class="post-metadata">

**Author:** ![Rui](https://avatars.discourse-cdn.com/v4/letter/r/ba8739/32.png) [@Rui](https://gromacs.bioexcel.eu/u/Rui)\
**Post date:** [December 22, 2024, 2:15pm UTC](https://gromacs.bioexcel.eu/t/how-to-use-clay-ff-with-gromacs-for-molecular-dynamics-simulation-of-mica/10881/1 "2024-12-22T14:15:10Z")

</div>

Hi ,  
I want to study the adsorption of reagents on the muscovite surface.  
I am trying to use Clay.ff with Gromacs.  
Clay.ff downloaded form the link: ([GitHub - thomasunderwood/ClayFF: A GROMACS implementation of the ClayFF force field](http://github.com/thomasunderwood/ClayFF))  
But I can not find the bonding parameters in ffbonded.itp for all atom of muscovite.  
Or if there is any other force field for muscovite and how to creat the topology file of muscovite?  
I have been stucked in this problem for many days. Can any help me? Thanks very much

---

<div class="post-metadata">

**Author:** ![askforarun](https://avatars.discourse-cdn.com/v4/letter/a/c5a1d2/32.png) [@askforarun](https://gromacs.bioexcel.eu/u/askforarun)\
**Post date:** [January 2, 2025, 11:12pm UTC](https://gromacs.bioexcel.eu/t/how-to-use-clay-ff-with-gromacs-for-molecular-dynamics-simulation-of-mica/10881/2 "2025-01-02T23:12:35Z")

</div>

In your model is the clay rigid or flexible ?

---

<div class="post-metadata">

**Author:** ![Rui](https://avatars.discourse-cdn.com/v4/letter/r/ba8739/32.png) [@Rui](https://gromacs.bioexcel.eu/u/Rui)\
**Post date:** [January 8, 2025, 2:19am UTC](https://gromacs.bioexcel.eu/t/how-to-use-clay-ff-with-gromacs-for-molecular-dynamics-simulation-of-mica/10881/3 "2025-01-08T02:19:20Z")

</div>

Dear askforarun

Thank you very much for your response. In my model, the clay is frozen, does this mean that the bond parameters of the clay are not particularly important?