Implementing a new Reactive Forcefield

[AsifShah]

Dear Admin,
I am trying to implement a new MoS2 reactive forcefield file into my simulation. The method I use is as follows:

1. I send the configuration to workflow and then choose the forcefield calculator. Then, I choose some inbuilt forcefield and add an optimization block.
2. Then send it to editor with all details present.
3. Then change reactive forcefield parameters with new parameters.

But it shows error. I have attached the python file and as well as the reactive forcefield parameter file. Kindly look into the matter why simulations shows error "Segmentation fault- core dumped"

Thanks
Regards

[AsifShah]

Dear Admin,

Kindly help in this regard.

[Julian Schneider]

Dear Asif,

We will check and get back to you.

[AsifShah]

Dear Admin,
Sure. Looking forward to your response.

[AsifShah]

Dear Admin,
Any update on this?
Here is the actual paper where this ReaxFF was implemented


https://pubs.acs.org/doi/10.1021/acs.jpcc.2c01075#:~:text=Abstract-,We%20present%20a%20new%20Reactive%20Force%20Field%20(ReaxFF)%20parameter%20set,simulations%20of%20Mo%2DS%20structures.&text=2-,structure%20in%20crystallization%20simulations.,unknown%20rock%20salt%20MoS%20structure.

[AsifShah]

Dear admin,
Is there any update on this? I am not able to find a solution.

[AsifShah]

Dear Admin,
Is there any update? It is urgent.

[Anders Blom]

I am trying to be constructive here, not criticize or downplay the problem, but why not fit an MTP for MoS2 instead, if this is urgent (and anyway)? I am sure it will be more accurate and perform similarly.

Even if that takes a week month to fit (which it doesn't, more like a day or two) it's faster than waiting for us to try to debug someone else's potential which may have some mistake in the numerical values, which is almost impossible to sort out; even if we see where the problem is we don't know the correct value. Or if there is a mistake in our ReaxFF model, it will take until September to release a new version...

So in any scenario you are probably much better off fitting a new MTP, and after that you can use the same procedure for any other 2D material (for which there are no ReaxFF potentials).