Monolayer protected Gold NP

Thanks for the advice Aidan, I checked them.
Anyway, my gold core is still exploding at the beginning of minimization. The gold core was modeled with the create_atoms command, filling a region with a lattice constant of 4.08 (lattice constant for Gold clusters); Force field parameters were taken from here http://link.springer.com/article/10.1007%2Fs10853-012-6356-8.

Now, the problem is: the gold core. Coating chains are behaving as “normal”, oscillating around their positions, but the core is exploding like is part of another system.

Suggestions on how to fix that or how to put a constraint for my gold atoms?

Thanks

Michele

Thanks for the advice Aidan, I checked them.
Anyway, my gold core is still exploding at the beginning of minimization.
The gold core was modeled with the create_atoms command, filling a region
with a lattice constant of 4.08 (lattice constant for Gold clusters); Force
field parameters were taken from here
http://link.springer.com/article/10.1007%2Fs10853-012-6356-8.

but those parameters look very different from what you have in your data file.

axel.