Higher-Order Corrections to Algebraic Derivation of Electric Dipole-Dipole Interaction

Nucleons are fermions with intrinsic spins exhibiting dipole character. Di-pole-dipole interaction via their dipole moments is the key feature quantify-ing the short-range nucleonics interaction in two-body physics. For a pair of interacting dipoles, the energy of a pair is the quantity of interest. The same is true for chemical polar molecules. For both cases, derivation of energy almost exclusively is carried out vectorially [1]. Although uncommon the interacting energy can be derived algebraically too. For the latter Taylor, expansion is applied [2]. The given expression although appears to be correct it is incom-plete. In our report, by applying Taylor’s expansion up to the 4th order and utilizing a Computer Algebra System we formulate the missing terms. Our report highlights the impact of correcting missing terms by giving two expli-cit examples.


Introduction
Reviewing literature on the format of dipole-dipole interaction energy in three areas of physics, magnetic, nuclear, and electric shows the similarity of the formulations. Denoting dipole moments by, µ  , σ  , and m  , respectfully, these are, The constant coefficients C's are associated with the coupling constants of each case. In these expressions, r is the vector distance between the dipoles with r being its associated unit vector.

Calculation
These expressions have been derived applying vector analysis e.g. [3]. For instance, the derivation of (1) is based on the gained energy of a magnet with a magnetic dipole 2 µ  in the magnetic field of another magnet with a dipole mo- Motivated with [2] we revisited its algebraic derivation. We consider a pair of electric dipoles in a 2D space. Figure 1 describes the scenario.  where ∑' is the algebraic sum, and dist i is the distance between the pairs, namely i = 1, 2, 3 and 4. Applying this notation according to vectors shown in Figure 1 the potential energy reads, As noted, for the sake of simplicity in (4), 12 r  is replaced with r  . The first two terms and the last two terms in (4) correspond to lines 1, 2 and 3, 4, shown in Figure 1, respectively. Note also according to notation: 1 . These are conducive to dipole moments, 1 One notices all four terms in (4) are common in 1/r and since we are confined in a 2D space this means The objective of the calculation is to apply Taylor expansion, as such the latter fraction is expanded for each of the four terms of (4) according to their increments. Let's define Expanding the first term of (4) with its associated increment gives, In (6), O 3 is the 3 rd order term, its format is like the 2 nd and the 4 th term. Because of its collective null value, it is not given explicitly. Noting also 4 , etc. Equation (6) must be applied to the remaining three terms of (4).
Realizing the first two terms of (4) are negative and the last two are positive yields to the cancelation of the first term of (6) when applied to (4). The second term of (6) contains first-order partial derivatives. Increments of four terms of (4) are such that when the second term of (6) is applied to (4) gives zero. Further analysis shows in general all odd orders of (6) when applied to (4) yield zero. In other words, only even terms of (6) survive. i.e., Taylor expansion of (4) may replace (4) in terms of only even powered terms.
As such we calculate the first even powered term of Taylor expansion, this is, Equation (7) reveals the functional format of the expression. Simply put denominators are of form 1/r 3 and 1/r 5 , these are compatible with the functional formats of (1)-(3). Now each term of (7) needs to be multiplied by their asso-ciated increments. Carrying out this algebraic simplification is cumbersome, that is the reason it is carried out applying a CAS, specifically Mathematica [4]. Nonetheless, its simplified output needs to be format-ed using the defined 1 m  and 2 m  in the text. This step requires regrouping of terms and is handled by hand.
The result is identical to (3).
Successfully this algebraic procedure bypasses the vector analysis producing the standard energy associated with the dipole-dipole interaction. However, Taylor expansion should not be terminated to the 2 nd order! As mentioned earlier higher even order terms should be included. It is one of the objectives of this report to investigate the impact of the next, i.e., the 4 th . One may argue contribution of terms beyond the 2 nd order "should" be ignored. Reasoning these terms would be falling off vs. distance as 1/r 7 and 1/r 9 . Of course, the argument holds so long as separation distance is large, however, for short distances this is not true. Even more interesting, here we show the 4 th order term contains a term such as 1/r 5 , and that needs to be added to the term of the same format that comes about from the 2 nd order Taylor expansion. Here is the proof. Noting the 4 th order term of the Taylor expansion as pointed out in (6) contains 4 th order partial derivatives we form its derivatives, x y x y x y x y x y x y x y x y x y x y x y x y + with r, we see that there are three terms in proportion to 1/r 5 ; these cannot be ignored and ought to be added to the terms of the same format when the 2 nd order expansion is considered. The other terms in (8) using the r notation are of the form 1/r 7 and 1/r 9 ; depending on the situation these terms may be ignored. It needs to be noted that it is obvious that higher-order expansion terms appear to contain overlapping 1/r n terms of the lower order. E.g., a sixth-order Taylor expansion would contain a term such as 1/r 7 and that needs to be added to the 4 th order, etc. Now as it is done for the 2 nd order, terms of the (8) are to be multiplied with their associated increments and then to be factored with 1/r 5 . This is indeed a very labor-intensive algebraic procedure. This is carried out applying CAS, Mathematica. The Equation (9) is compared to (3). The second expression in the angled bracket is the correction term. Its dimension is L −1 as is the dimension of the expression in the curly bracket. The constant C 3 is kq 2 where the k is the electrostatic coupling constant and q is the charge of the dipole making the dimension (9) energy as expected.
To illustrate the impact of the correction term we consider two instances. without and with correction terms (3) yields, respectfully. 3 3 5 In (10)    In Figure 3, the blue curves are the uncorrected expressions in (10) and (11).
The red curves are the corrected expressions. It is notice-able the correction terms have no neglect-able impacts. For the case shown in Figure 2(a), the impact is large, while for Figure 2(b) it is not severe.
We shown two instances to highlight the impact of the corrections. There are other cases correlated with Figure 2, such as antiparallel-oriented dipoles. These and other cases are left to the interested reader to explore.

Conclusion
Dipole-dipole interaction conducive to either mutual paired forces or energy redundantly is utilized in various areas of physics. Derivation of these quantities exclusively is based on vector calculation. An attempt was made in the referenced article to derive the same expressions applying algebra, specifically Taylor expansion that bypasses the latter. Although it is shown the goal is "achievable", but the calculation falls short by limiting the number of the expansion terms. In our investigation explicitly we have shown higher-order terms beyond the second order ought to be considered. This stresses that although the algebraic approach is novel, it must be done with much care.