Advanced Electron Simulation in Copper Wire

I have a simulation in Python/CuPy of electrons in a wire that I want to translate into Rust (or Vulkan, or OpenCL, or CUDA, or Julia or fix the Python so it uses far less memory, maybe calling PyCUDA) and make a major modification. The simulation calculates the coulumbic forces between all pairs of electrons every delta-t and then after that uses the force to calculate the new positions and velocities. It also outputs a graph showing the density of electrons at each slice of the wire at each delta-t. From this I can see how fast a signal is propagating down the wire in the simulation. I have tried to allow for the propagation delay between the electrons by using the velocity to estimate where the electron would really have been back in time and calculating the force from that point but this is an approximation I would like to improve on. The major modification is to store the positions and velocities of each electron for the last N delta-t time slices. Then when calculating the forces use past positions of electrons that are at the right distance. Also want to make an adjustment in the force if it is moving forward or away (will have some constant). I am looking at a sort of a philosophical point, can what we know of as "magnetic fields" just emerge as the mathematical consequences of time-delayed electric fields and large numbers of electrons? This would mean magnetic force is not really a separate thing from electric force. Some scientists such as Parrott, Schwarzschild, Tetrode, Fokker, Hoyle, and Narlikar, Wilhelm Weber, André-Marie Ampère, John Wheeler, and Richard Feynman have said similar things. I believe that Maxwell's equations and Einstein's relativity are a way to explain how large numbers of electrons interact, but that on an individual electron level they are doing something simpler like delayed electric fields. Individual electrons would not be able to solve Maxwell's or Einstein's equations, so must use something more simple and fundamental. The larger and more complex behavior that Maxwell and Einstein explain resulting from the interactions of large numbers of simple electrons. I have a RTX A600 GPU with 10,000 cores and 48 GB of RAM. This is around the power of the most powerful computer in the world 22 years ago. This is enough power that I am able to simulate tens of thousands of individual electrons. For a tiny wire like inside a computer chip, this is enough electrons to model what is going on. My goal is to get reasonably accurate results without using magnetic fields or Maxwell's equations. My Python code now is around 400 lines. The code must make good use of my RTX A6000 GPU. It should be able to handle 50,000 electrons and 100 past positions for each. It should be clean and well documented code so that someone new to language used can understand what is going on.