Chemical Sciences: A Manual for CSIR-UGC National Eligibility Test for Lectureship and JRF/Jacobi coordinates Jacobi coordinates for two-body problem; Jacobi coordinates are ${\boldsymbol {R}}={\frac {m_{1}}{M}}{\boldsymbol {x}}_{1}+{\frac {m_{2}}{M}}{\boldsymbol {x}}_{2}$ and ${\boldsymbol {r}}={\boldsymbol {x}}_{1}-{\boldsymbol {x}}_{2}$ with $M=m_{1}+m_{2}\$ . A possible set of Jacobi coordinates for four-body problem; the Jacobi coordinates are r1, r2, r3 and the center of mass R. See Cornille.

In the theory of many-particle systems, Jacobi coordinates often are used to simplify the mathematical formulation. These coordinates are particularly common in treating polyatomic molecules and chemical reactions, and in celestial mechanics. An algorithm for generating the Jacobi coordinates for N bodies may be based upon binary trees. In words, the algorithm is described as follows:

Let mj and mk be the masses of two bodies that are replaced by a new body of virtual mass M = mj + mk. The position coordinates xj and xk are replaced by their relative position rjk = xj − xk and by the vector to their center of mass Rjk = (mj qj + mkqk)/(mj + mk). The node in the binary tree corresponding to the virtual body has mj as its right child and mk as its left child. The order of children indicates the relative coordinate points from xk to xj. Repeat the above step for N − 1 bodies, that is, the N − 2 original bodies plus the new virtual body.

For the four-body problem the result is:

${\boldsymbol {r_{1}=x_{1}-x_{2}}}\ ,$ ${\boldsymbol {r_{j}}}={\frac {1}{m_{0j}}}\sum _{k=1}^{j}m_{k}{\boldsymbol {x_{k}}}\ -\ {\boldsymbol {x_{j+1}}}\ ,$ with

$m_{0j}=\sum _{k=1}^{j}\ m_{k}\ .$ The vector R is the center of mass of all the bodies:

${\boldsymbol {R}}={\frac {1}{m_{0}}}\sum _{k=1}^{N}\ m_{k}{\boldsymbol {x_{k}}}\ ;$ $m_{0}=\sum _{k=1}^{N}\ m_{k}\ .$ 