Is the Jacobian matrix calculation process correctly on this energy function?
That is not an equation to be solved; it is an objective function to be minimized.  They apparently describe how they optimize it in Section 3.3.3 of the paper.  I suggest trying to understand that section.