I am trying to model a complex molecule containing gold and carbon atoms and bonds (organogold). May I know whether there are any parameters which I can use to simulate these in MD (either GROMACS or LAMMPS)?
The simple way is to look at x-ray structures of gold derivatives that contain gold-carbon bonds. The following are publications including x-ray structures of gold compounds that you can utilize the determined Au-C bonds tfor simulation calculations:
1- Hussong, M. W., Rominger, F., Krämer, P. and Straub, B. F.: "Isolation of a Non-Heteroatom-Stabilized Gold–Carbene Complex." Angew. Chem. Int. Ed. (published online 20 June 2014), DOI: 10.1002/anie.201404032
You could try out the improved ReaxFF parameters by Gyun-Taeck Bae and Christine Aikens for Au–S–C–H Systems. You could also check out other work by Aikens' group on organo-gold, depending on the size of your system perhaps you may be able to model it at the DFT level?