I just got this email from someone I never heard of. Any thoughts? (I thought some of the videos were pretty.) Jim Propp ---------- Forwarded message --------- From: Wolfhard Hövel <Wolfhard.Hoevel@th-nuernberg.de> Date: Sun, Aug 9, 2020 at 4:30 AM Subject: Discrete dynamical system, interacting particles with mass, spin and charge To: <jamespropp@gmail.com> Dear Professor Propp, I would like to propose an algorithm that simulates spontaneously interacting particles with mass, spin and charge, using only position and displacement vectors. Both the momentum and the energy conservation law are always fulfilled. You can see some results of the approach on YouTube. You can find the YouTube videos „Experiments with vectors“ under the following link: https://www.youtube.com/playlist?list=PLBda4QvOnxQHVsqa3l2ixuqD7xlRlYWVx Also you can view various examples interactively on the web page of Ulrich Schwebinghaus: http://www.fraktalwelt.de/myhome/selforg_2.htm I have attached a description of the calculation model below. There I tried to describe the algorithm “Experiments with Vectors” as far as possible with words. I would be happy if you could take a look at them. Sincerely yours, Wolfhard Hövel D 92268 Etzelwang Gerhardsberg 7 Email: wolfhard.hoevel@th-nuernberg.de Telephone: +491735737624 Experiments with vectors A digital computer is a wonderful tool for performing mathematical experiments. I would like to introduce an algorithm that spontaneously generates particles that move and interact in Euclidean space. This only happens with the help of position and displacement vectors. First of all I have to define some rules according to which the algorithm should run. These rules are not entirely random and when experimenting quickly show whether the rules produce something interesting. I believe that the proposed algorithm actually gives amazing results. The program describes a discrete dynamical system. In a loop calculation, the states of the system change in equidistant time steps. With each iteration, the states are always recalculated according to specified procedures. A general rule is introduced that a quantity G is made up of the integer multiples of a constant base quantity G0 (G = n G0, n is a natural number). All quantities are finite. Now I can define an important quantity in this way: time. Each cycle in the calculation loop specifies the shortest time step dt0. The time is then granular and corresponds to the number n of iterations t = n dt0. You also need objects with which something happens. Let's start with a point. The point has a certain position in space, which is determined by the position vector r (vectors are printed in bold). This position vector points from a freely chosen zero vector to the specified point. The point is by definition always in motion, I achieve this with a displacement vector dr, which is added to the position of the point with every iteration in the calculation loop. Since dr has by definition a constant length, the step width of the point is constant and is equal to the magnitude of the displacement vector |dr| = dr0. The point now moves straight ahead and has the constant speed v = dr / dt0. According to Newton's laws, this is a physical statement. After a while the point disappears from the screen and that would be a boring event. I have to constrain the movement of the point to a defined area. I achieve this by providing the point with a second point as a partner. These two points are supposed to dance around each other. For this, the dancers need certain rules. The first rule is: move with the velocities v and -v respectively exactly in opposite direction on two parallel lines with the distance s. The second rule is: If your distance becomes longer than the elongation e, turn around and run on two new parallel lines with the distance s. This distance s and the direction of rotation remaining always the same. A third rule should take into account the principle that all quantities can only change step by step, including the parameters s and e. They are integer multiples of the step width dr0. The two points form a system whose center of mass remains constant. The angular momentum of the points is constant and given by v and s. In addition, a maximum distance is set, which the dancing points should not exceed if possible. They don't quite succeed because they move step by step. The limit is not sharp, it can be exceeded by a maximum of dr0. There is still something missing to make the algorithm physical: it is the mass. So I assign a smallest possible mass m0 to each point. This assumption must not appear artificial or forced. A common center of mass can also be specified for a point cloud. Now the constant momentum for a point is given by p = m0 v and the kinetic energy by E = ½ m0 v². The described pair of points ("bivector") is the most important basic building block of the calculation model. It is the smallest possible particle of the simulation. From a geometric point of view, a parallelogram is spanned in space by the distance vector r12 = r2 - r1 between the two points and the displacement vector dr, whose area is equal to s · dr0. The parallelogram is an oriented plane and has a sense of rotation. The bivector exists in every dimension > 1. The bouncing points describe a circular disc, the edge of which is not very sharp. For example, if e >> s the points are oscillating in an almost star shape around the circle with the diameter s and for the special case e = s the points are moving unsharp on a circle with the diameter e. The bivector contains two points and therefore its mass is 2 m0. How can the bivector be represented in a program? At the beginning of each iteration, it is checked whether the distance between the points exeeds e. If this is not the case, the points can continue to move freely in a straight line. If the permitted limit is exceeded, the points must be reflected so that their distance is reduced again. The rules given above must be observed. From a geometric point of view, this is easy to illustrate. The distance between the points is halved and a circle with the diameter s is drawn around the center found, which is also the center of mass of the bivector. Then, starting from the two points, tangents are placed on the circle, whereby the sense of rotation must be the same as before. The points continue to move on these newly found tangents. One can imagine that the two points are connected by a chain with the length e. As long as the chain is still loose, the points can move freely. However, if the chain is suddenly tensioned after a step, the points must reverse abruptly. The chain can be spanned by a maximum of one chain link. This methode of reflection has the advantage that only internal quantities of the bivector are used and no impulse is exchanged to the outside. The bivectors can be moved freely in the space. This fact will become important later. Now we add more points. Each point must be at least a member of one bivector to give it a momentum. However, one point can be a partner of several bivectors and then the situation quickly becomes confusing. Therefore, it makes sense to represent the relationship between the points using a graph of graph theory. The points are the vertices and the bivectors are the edges of the graph. If the point cloud consists of N points, a maximum of N(N-1)/2 bivectors can be inserted between the points. If a point is involved in several bivectors, there is a strong coupling between these linked bivectors. This should be illustrated using a simple example. Suppose we connect three points P1, P2 and P3 with two bivectors with the pairing 1-2 and 2-3. With each iteration, each bivector is treated as before, regardless of the others. The associated displacement vectors dr12 and dr23 result from this. Finally, the displacements are superposed. Point P1 shifts by dr12, point P2 by (-dr12 + dr23) and point P3 by -dr23. Imagine that dancer P2 has to perform the steps of two dance pairs at the same time! The two dance couples will have to dance completely new figures, they are closely linked. This example describes another possible particle. Since it contains three points, its mass is 3 m0. If you look at the energy situation, it gets a little more complicated. Two displacement vectors -dr12 and dr23 are added to point P2. In the case that both vectors have the same direction, the point P2 has the speed 2v. Its kinetic energy is therefore 2 m0 v². On the other hand, if the two displacement vectors are exactly opposite, the kinetic energy is zero. So that the theorem of energy conservation remains valid for the simulation model, the energy of point P2 could now be interpreted as potential energy. As the directions of the displacement vectors will change, the potential energy is converted back into kinetic energy. Nothing is lost, the number of displacement vectors is constant. This interplay is also known from a tensioned and relaxed metal spring. If you experiment with the calculation model described, you will experience a surprise. Under certain conditions, self-organization creates stable structures in the space. These attractors can change spontaneously, they can attract and repel, they behave like particles with mass and obey the momentum conservation law. To do this, however, the algorithm still needs to be expanded. A weak interaction must be added. Let's start again with a point cloud. The points are linked by bivectors. When points are members of several bivectors, there are strongly coupled sets of points, these are our particles. Several isolated graphs and thus particles can be defined. The particles do not yet interact, their graphs are isolated networks in space. How can I couple the isolated particles? To do this, I assign a plus or minus sign to some points. The additional rule is: if the points have opposite signs, the points move iteratively radially towards each other with the step size drL, if they are like-signed, they repel each other. These shifts are added to the other displacements. The momentum conservation law is still observed, the sum of the momentum vectors remains equal to zero. With drL you can control the weak coupling between the particles, often I set drL = dr0 / 150. The algorithm can be programmed on the basis of real or integer numbers. Further information on this particle simulation can be found on the Internet under the keywords „search opus4 hoevel“ . Appendix You can make a bivector geometrically clear with a drawing. Draw two points P1 and P2 16 cm apart. Connect the two points with a straight line and halve them. The point M found in this way is the center of mass of the pair of points. Draw two circles around M with the radius s/2 = 2 cm and e/2 = 4 cm. Place two parallel tangent lines starting from points P1 and P2 to the s-circle. Starting from point P1, move P1 step by step with the step width dr0 = 2 cm along the tangent. The point now jumps to the e-circle and runs through it. As soon as P1 has left the e-circle, a new tangent line from P1 is placed to the s-circle, keeping the direction of rotation, and the game is repeated. The point P1 is now captured and circles the s-circle. Point P2 is always moved in the opposite direction at the same time. This system is a bivector, it is characterized by a constant angular momentum (~ s) and an elongation e. The bivector can be moved freely in space, its position is only determined by the two points P1 and P2.