Social Distancing via Coulomb’s Law

This research effort addresses the social-distancing problem. As the COVID-19 pandemic continues, we’ve learned the importance of keeping proper distance, so as to avoid (or minimize) the spread of infection. For this paper, individuals are represented as positively-charged particles, behaving in accor-dance with Coulomb’s Law. Additionally, negatively-charged stationary (non-moving) particles are positioned such that their attraction to the posi-tively-charged particles guides the movement of the positively-charged particles in a desirable fashion. During a simulation process, Coulomb’s Law guides particle behavior such that the positively-charged particles arrange themselves in a way such that their spacing is essentially optimal. Of course, these positively charged particles can be thought of as a surrogate for individuals, resulting in the optimal spacing of individuals.


Introduction
COVID-19 became a major problem in the world in the early months of 2020.
Countries have used varying approaches in addressing the problem, but the concept of "social-distancing" has become one of the most important defenses against the virus. Two major factors contribute to the spread of infection: population density [1] and transience [2]. The number of cases is disproportionately higher in densely-populated areas, such as New York City and Madrid. Additionally, these two cities also have a lot of "transience." That is, many people move through these cities, due to the transportation infrastructure, leaving these cities more exposed to possible infections. New Orleans, Louisiana is a good example of transience being a problem. New Orleans is not a large city compared to those mentioned above, but the city held its annual Mardi Gras festival health officials realized holding Mardi Gras was a mistake, given that so many people in New Orleans became infected, presumably because of the large number of people who came to New Orleans [3].
In short, social-distancing is the best defense against the virus. Because of the importance of social distancing, a good social distancing tool is desired. The main motivation of this research effort is to assist with social distancing-finding the best way to address social distancing via our vast array of mathematical tools.
These tools can be used to optimize the spacing of individuals in our collective attempt to minimize transmission of the virus. The body of literature contains a good amount of research relating to the circle-packing problem-a problem committed to optimizing placement of circles in a confined area such that the total area dedicated to circles is maximized. The research presented here slightly differs in that because particles, not circles, are placed into a unit square.
This research exploits Coulomb's Law in an attempt to optimize the placement of charged particles in a confined environment. In this context, the charged particles are surrogates for individuals. Coulomb's Law guides the charged particles as to where to go. This process is organized as a part of the meta-heuristic Simulated Annealing-a search process engineered to find near-optimal solutions for all sizes of optimization problems. The pursuance of Coulomb's Law in this endeavor is considered the unique feature of this research.
In short, the contributions of this paper are summarized as follows: • The social distancing problem is addressed, specifically via mathematical means • Coulomb's Law is employed in a novel fashion • Near-optimal results are found for problems proven optimal • The presented methodology can be employed to larger problems without loss of generality The following sections describe Circle-Packing/Particle Arrangement, Coulomb's Law, and how Coulomb's Law is used to assist in particle movement via the Simulated Annealing search heuristic. The presented methodology is also used on sample problems resulting in near-optimal solutions. General conclusions are also offered, along with recommendations for subsequent research.
It should be noted that there are several terms defined in this paper. Specifically, Table 1 defines terms related to forces and Coulomb's Law, Table 2 shows an example of this, while Table 3 defines terms related to the Simulated Annealing Process.

Circle-Packing and Arrangement of Particles
In its most general form, "circle-packing" is the practice of placing non-overlapping circles in a confined space. The circles need not be of the same diameter. The intention is to arrange the circles such that the wasted space (non-circle area) is minimized. This practice is important in the steel-stamping can only be determined via mathematical proof for a specific number of circles, and optimality has only been proven for small values of n. Optimal solutions for larger values of n remain an unresolved issue [7].
The research here differs from the above problem in that we are only interested in optimally spacing the centers of the n circles of concern. As such, we treat the centers of the circles as points on a two-dimensional plane, and wish to maximize the distance between centers of the two closest circles. Figure 1 shows optimal circle packing for n = 2. Maximizing the distance between the two closest circle centers is the objective for this research. Given that the larger values of n result in problems with unknown optimal solutions, we are motivated to provide a tool for addressing not just "large" values of n, but any particular value of n. From this point forward, these circle centers are referred to as "particles," and in the context of social distancing, we wish to maximize the distance between the two closest particles (persons).

Particle Movement
In the context of this research, there are two ways to move individuals with the intent of social distancing in an optimal fashion. First of all, individuals hereafter are referred to as "particles," which are simply points in space. It is desired that these particles be arranged such that the minimum pairwise distance between all particles be maximized. In order to accomplish this objective, the particles are moved during a simulation process. For this research effort, two different ways to move the particles are pursued. The first is by exploiting the magnetic nature of charged particles, which is referred to as "polar movement." The second way to move the particles is via "random movement."

Coulomb's Law and Polar Particle Movement
Magnetic particles are either positively charged or negatively charged. Particles of the same charge repel each other, while particles of differing charges are attracted to each other. Figure 2 shows this: (a) and (b) show repulsion, while (c) shows attraction.
Coulomb's Law quantifies the force between the particles according to the Coulomb's Law, which mathematically is as follows: The first term of this (1/4πε 0 ) is essentially a constant, used to convert the result into units of Newtons (N). This constant is not used here, as we are only interested in relative forces. The value of q i and q j are the charges of two particles, measured in Coulombs (C). The value of d ij is the Euclidean distance between particles i and j. As such, the force between two charged particles reduces to the following: We can state that the force between two charged particles is the product of their charges and inversely proportional to the squared distance between them, where the squared distance is as follows: With a basic understanding of Coulomb's Law, we now need to understand how the forces of several charged particles affect a single particle. For this research, only positive particles move, while negative particles remain stationarythe forces imposed by the negative particles do, however, affect the movement of the positively-charged particles. Table 1 provides some relevant terms: It is intended to determine the net force exerted on all positive particles. To do this, we combine the forces exerted on all positive particles in terms of horizontal and vertical components. Consider Figure 3, detailing a specific particle at the point (7,11). It is our intent to sum all forces imposed on this specific particle to determine the net force acting on this particle.  The find the net force exerted on some specific positive particle i + , the following two equations are used: It is important to note that the first component of the above two formulae relates to the force between the particle i + and other positively-charged particles, while the second component relates to the force between positively-charged particle i + and the negatively-charged particles. It is also important to note that the (n + /n − ) component is used to make sure the total positive force (in Coulombs) is the same as the total negative force (in Coulombs). For this example, there are four positively-charged particles, and one negatively-charged particle. As such, the (n + /n − ) multiplier assures the negatively-charged particle has a charge of four Coulombs, the same as the sum of the positively-charged particles.
The magnitude of the resultant net force is as follows: The direction, or heading, of this resultant net force is as follows: Using Equations (4) and (5) above, and summing the two right-most columns of Table 2, we have ( )   The importance of Figure 4 is that the resultant force (both magnitude and direction) presents the location to where the particle is inclined to move. Hereafter, this type of particle movement is referred to as "polar movement." Again, it is important to note that only positively-charged particles are subject to movement, so resultant forces are only determined for the positively-charged particles.

Random Particle Movement
Unfortunately, polar movement is not always possible, because resultant forces may take us outside the realm of feasible solutions (outside the unit square). Because in polar movement, only the positively-charged particles move, and their nature to repel each other will motivate these particles outside the realm of feasible solutions. Figure 5 shows an example of this. Figure 5 shows us that the positively-charged particle desires to move outside the search space, rendering an infeasible solution. Furthermore, the heading of this particle will never change, essentially putting our search at a "dead end." In order to combat this problem, all movable particles are treated as "neutral" particles, and they are moved a very small amount according to randomly-selected magnitudes and headings. This approach is referred to as random movement, and is illustrated in Figure 6.
This research effort utilizes both polar movement and random movement in our attempt to find optimal spacing.

Search Methodology
This section describes the methodology used for finding the best spacing of particles according to the movement rules described above.

Initialization
Prior to describing the search process, Table 3 provides definitions, some of which are a slight repeat of the definitions from Table 1.
Values for the following search parameters are chosen: Iter, k B , CR (< 1.0), T 1 and T F (where T F < T 1 ), and T = T 1 . The iteration index, q, is set to 1. The number of movable particles for spacing (n) is chosen, and for each particle i, x and y coordinates are assigned the following: This is the initial solution. The objective function of this initial solution is as follows: The values of Z C , Z T and Z B are set equal to the objective function value. In layman's terms, we wish to arrange n particles in the unit square such that the distance between the closest two particles is maximized. The mathematics of this model shows a scenario where the objective function is both nonlinear and discontinuous, meaning conventional solvers are unlikely to guarantee an optimal solution, especially for larger problems. This is the reason a search heuristic is pursued.

Movement
Prior to movement, the current solution is copied to the temporary solution.
Movement is either in the form of polar movement or random movement, both described in Section 3. The probability of polar-movement is τ, implying that the probability of random movement is (1 − τ).
In polar movement, only the positively-charged particles move, while the negatively-charged particles remain stationary. All positively charged particles simultaneously move in the direction of their heading ( ) i h + at a distance of ( ) + i F . In random movement, the n-particles move to a nearby, randomly-selected location. Regardless of the type of move, infeasible moves are not executed. The result of all movements is the temporary solution.

Evaluation
Once the (previously-described) move is complete, the objective function value (Z T ) is determined via Equation (10). A determination needs to be made as to whether or not to accept the temporary solution as the current solution. If the objective function value of the temporary solution is better than the objective function value of the current solution-that is, if Z T > Z C , the current solution is replaced by the temporary solution.
Otherwise, further evaluation is needed before accepting the temporary solution as the current solution. The inferiority of the temporary solution compared to the current solution is represented as δ, and is determined as follows: The probability of accepting the inferior temporary solution as the current solution occurs with the following probability: A uniformly-distributed random number of the interval (0, 1) is generated. If this value is less than P(A), the inferior temporary solution is accepted as the  (14) and (15)

Search Status
If q is less than the specified number of iterations per temperature level (Iter), the value of q is incremented by one (q = q + 1), and control if the search returns to that described in Section 4.2. Otherwise, the current Temperature (T) is multiplied by the cooling rate to set the new current temperature (T = T * CR), and q is reset to 1. If the new resultant temperature is greater than T F , control of the search returns to that described by Section 4.2, and q is reset to 1. Otherwise, the search is over and the best solution is returned.

Search in Pseudocode
For a better understanding of the search process, details are provided in pseudocode form.

Experimentation
The methodology presented above was carried out via agent-based modeling software. Agent based modeling is an extension of traditional object-oriented programming (such as Java or C++) in that the modeling environment is a dynamic system-the entire system is perpetually in a state of change [11]. The specific programming environment was NetLogo [12] [13].

Experiment in NetLogo
The methodology can be used for any value of n. As stated earlier, NetLogo was the programming software employed. The images provided via Figure 7 and Figure 8 are images from the NetLogo environment. Figure 7 details the randomly-generated initial solution for n = 18. The particles to be arranged are positively-charged, and shown via a blue "+" symbol. Additionally, the size of the particle is proportional to the net force acting on the particle, and the heading of the particle is also shown via the arrow at the end of the "+" symbol.  particles was 0.3000, where the proven-optimal distance was found to be 0.3005 [14]. Figure 8 shows what the NetLogo program found to be the best overall solution. Here, the positively-charged particles are shown via black circles, and the negatively-charged particles are not shown. The maximized distance between the two closest particles is also shown.

Computational Experience
In finding the best solutions for values of n ≥ 2, the search settings shown in Table 3 had to be adjusted for each unique value of n. For the trivial problems involving optimal spacing of 2, 3 and 4 particles (n = 2, 3, 4), the search best responded by eliminating the negative charge in the unit square. This makes sense, because for these values of n, the positively-charged particles are motivated to get as far away from each other as possible, regardless of other circumstances.
Staring with five or more particles (n ≥ 5), the negative force in the central portion of the unit square becomes more important. This is because for larger quantities of particles, some of them need to be attracted toward the central part of the unit square for optimal spacing. Adjusting the negative force from problem to problem was required to obtain desirable results.
For all problems, T 1 was set to 40, while T F was set to 1. The cooling ratio, CR was set to 0.9995, and the number of Iter increased in proportion to n. For all problems, the Boltzmann Constant, k B , was set such that a temporary solution 5% inferior to the current solution had a 5% probability of acceptance when T = T 1 . Table 4 shows the results of the methodology presented here as compared to results proven to be optimal for several problems. For the problems pursued, the presented methodology averaged 99.89% as compared to proven optimal. The worst-performance was for n = 15, which was 99.27% as compared to optimal.

Concluding Comments
Methodology has been provided to address social distancing via Coulomb's Law and a variant of the circle-packing problem. Results are encouraging. As the COVID-19 virus and its social impact continue to evolve, social distancing will likely be with us for a while, and what we've learned from the effects of social distancing will likely be applied to future crises of a similar nature.
The presented research most certainly can be taken further, particularly with regard to scale. While the virus continues to force us into social distancing, sporting events will have limited capacity, providing employment for the presented methodology. The same can be said for worship services. In these larger-scale situations; however, something other than a unit square will be required, given the differing shapes of large venues.
Another opportunity for subsequent research applies to the use of Coulomb's Law, and more effective ways of using the law to find desirable spacing. Scientific principles of other types could also be investigated in attempt to achieve desirable spacing.
In summary, the findings of this work are as follows: • Social-distancing in an optimal fashion can be addressed via an applied mathematical approach • Near-optimal results were found when compared to existing problems with known optima • A combination of Coulomb's Law and random movement was needed to obtain desirable solutions • The presented methodology can be scaled to larger problems without loss of generality • The work exploited a novel application of Coulomb's Law It is possible that other scientific principles (other than Coulomb's Law) would be employed to address this problem.