SHORT VERSION: I have written my fixes and pair styles to allow concentration diffusion. It works nicely in serial, but it chokes in parallel. What are the common errors that people make when implementing new fixes/pair styles that I could look into to fix this?
LONGER VERSION (my apologise...): I don't think this is the right way to do it, but I got something to work (also it helped understanding how LAMMPS structure). I defined a new property using fix property/atom in the input script, here:
Then I added a new fixes and pair_style to allow the diffusion of the property here (the concentration cA):
While I would truly appreciate being criticised about code-styles and good/bad practice (please do! I am still learning), I would love to know which step I missed out? Do I need to regenerate/update the values of the ghost atoms (from what I collected here http://lammps.sandia.gov/doc/Developer.pdf
) at each time step? Bear in mind the code is still simply a proof on concept at this stage - I have not sanitise all inputs errors and make it truly robust...
Also, I have a feeling like my approach is fundamentally wrong. Some instincts tell me that it would be better to implement the new properties inside the atom styles instead of introducing them anew in each script?