Again, 𝜌𝜌 � 𝑎𝑎 (𝐫𝐫), 𝜌𝜌 � 𝑏𝑏 (𝐫𝐫), and 𝐷𝐷(𝐫𝐫) are defined according to eqs 2, 3, and 5 using the newly defined 𝑓𝑓(𝑛𝑛) in eq 13. New method 4 is an improvement to new method 3 to remove the discontinuity at 𝑛𝑛 = 1. Here we simply mention two key points: (i) Equation 13 is a hyperbola that satisfies 𝑓𝑓(𝑛𝑛) ≤ min(𝑛𝑛, 2 -𝑛𝑛), which also makes 𝐷𝐷(𝐫𝐫) ≤ 𝜌𝜌(𝐫𝐫) for all points in space. (ii) Equation 13 has one adjustable parameter, which it is convenient to take as 𝑚𝑚 because this is the value of 𝑓𝑓(𝑛𝑛) at n = 1; if one sets 𝑚𝑚 = 1, method 4 reduces to method 3. Below we will discuss the effect of the choice of 𝑚𝑚 on the performance of the DC functional. With some trial and error, we determined that 𝑚𝑚 = 0.96 is the most appropriate value for practical calculations in terms of accuracy and physical interpretation.