# Error "Unknown bond\_atomtype CC3151"

**URL:** https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173
**Category:** User discussions
**Created:** [September 12, 2023, 9:22pm UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173 "2023-09-12T21:22:20Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![silvio0101](https://avatars.discourse-cdn.com/v4/letter/s/4af34b/32.png) [@silvio0101](https://gromacs.bioexcel.eu/u/silvio0101)
#### Post date: [September 12, 2023, 9:22pm UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173/1 "2023-09-12T21:22:20Z")

</div>

GROMACS version: 2021.4  
GROMACS modification: Yes

Dear GROMACS users,  
I’m trying simulate a polymeric matrix composed by lignin + PLA system with charmm36 force field.  
I generated the GROMACS files with charmm2gmx. And there is a problem when I used the ‘gmx grompp’.

The command line below:  
gmx grompp -f minim-2.mdp -c syr\_box.gro -p topol.top o minim.tpr

When I run the command line occur this error:  
ERROR 1 [file ffbonded.itp, line 25]:  
Unknown bond\_atomtype CC3151

If somebody have any idea to solve this error? Thank you for reply

---

<div class="post-metadata">

### Author: ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [September 13, 2023, 2:38am UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173/2 "2023-09-13T02:38:13Z")

</div>

CC3151 is a standard carbohydrate atom type. Did you generate the entire force field yourself? That’s not necessary; the latest CHARMM toppar files were ported to GROMACS format last year and are available from Alex MacKerell’s website. If you aren’t use those files, you should. How did you generate the lignin topology?

---

<div class="post-metadata">

### Author: ![silvio0101](https://avatars.discourse-cdn.com/v4/letter/s/4af34b/32.png) [@silvio0101](https://gromacs.bioexcel.eu/u/silvio0101)
#### Post date: [October 3, 2023, 4:29pm UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173/3 "2023-10-03T16:29:09Z")

</div>

Thank you so much for your reply.  
I used the “Systematic parameterization of lignin for the CHARMM force field” made by Vermass et,al; 2018 and then convert into a GROMACS format using the charmm2gmx tool.

---

<div class="post-metadata">

### Author: ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)
#### Post date: [October 4, 2023, 1:16pm UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173/4 "2023-10-04T13:16:59Z")

</div>

To clarify, are you converting the entire CHARMM force field with our charmm2gmx program (which isn’t necessary, we’ve already done that for you) or some standalone topology?

---

<div class="post-metadata">

### Author: ![silvio0101](https://avatars.discourse-cdn.com/v4/letter/s/4af34b/32.png) [@silvio0101](https://gromacs.bioexcel.eu/u/silvio0101)
#### Post date: [October 4, 2023, 6:42pm UTC](https://gromacs.bioexcel.eu/t/error-unknown-bond-atomtype-cc3151/7173/5 "2023-10-04T18:42:54Z")

</div>

You are right, I see that CHARMM toppar files for lignin already exist.  
Thank you so much for your reply it was helpful.
