Phase diagram:

Changing the threshold in the Dual Energy implementation and the integration time step:

Adding a Second Order Gravity Integration Scheme:

From Nyx implementation: