Related papers: An algorithm for coalescence of classical objects …
Building accurate and predictive models of the underlying mechanisms of celestial motion has inspired fundamental developments in theoretical physics. Candidate theories seek to explain observations and predict future positions of planets,…
The law of centripetal force governing the motion of celestial bodies in eccentric conic sections, has been established and thoroughly investigated by Sir Isaac Newton in his Principia Mathematica. Yet its profound implications on the…
We describe the Reversibility Error Method (REM) and its applications to planetary dynamics. REM is based on the time-reversibility analysis of the phase-space trajectories of conservative Hamiltonian systems. The round-off errors break the…
Langevin simulation provides an effective way to study collisional effects in beams by reducing the six-dimensional Fokker-Planck equation to a group of stochastic ordinary differential equations. These resulting equations usually have…
When particles on a line collide, they may coalesce into one. Such systems arise in the voter model, where boundaries between opinion clusters perform coalescing random walks, and in reaction-diffusion theory, where diffusing particles…
An algorithm is proposed for finding numerical solutions of a kinetic equation that describes an infinite system of point articles placed in $\mathbb{R}^d (d \geq 1)$. The particles perform random jumps with pair wise repulsion, in the…
In this work, the relativistic phenomena of Lorentz contraction and time dilation are derived using a modified distance formula appropriate for discrete space. This new distance formula is different than Pythagoras's theorem but converges…
Towards the end of nineteenth century, Celestial Mechanics provided the most powerful tools to test Newtonian gravity in the solar system, and led also to the discovery of chaos in modern science. Nowadays, in light of general relativity,…
A novel view for the emergence of chaos in Lorenz-like systems is presented. For such purpose, the Lorenz problem is reformulated in a classical mechanical form and it turns out to be equivalent to the problem of a damped and forced one…
In finite-size scaling analyses of Monte Carlo simulations of second-order phase transitions one often needs an extended temperature range around the critical point. By combining the parallel tempering algorithm with cluster updates and an…
The first integrals of the Kepler problem are used to compute preliminary orbits starting from two short observed arcs of a celestial body, which may be obtained either by optical or radar observations. We write polynomial equations for…
In a quest towards an intelligent decision-making machine, the ability to make plausible predictions is the central pillar of its intelligence. A predicting algorithm's central idea is to understand the governing physical rules and make…
We develop a new algorithm for the Brownian dynamics of soft matter systems that evolves time by spatially correlated Monte Carlo moves. The algorithm uses vector wavelets as its basic moves and produces hydrodynamics in the low Reynolds…
Collisional fragmentation is a ubiquitous phenomenon arising in a variety of astrophysical systems, from asteroid belts to debris and protoplanetary disks. Numerical studies of fragmentation typically rely on discretizing the size…
There are many occasions when one does not have complete information in order to classify objects into different classes, and yet it is important to do the best one can since other decisions depend on that. In astronomy, especially…
Molecular simulations of many particles which move rather according to a brownian than a newtonian type of dynamics, nevertheless, can be performed by means of a "velocity-Verlet-like" algorithm. The derivation of this algorithm requires…
Linear stationary reaction-convection-diffusion equations with Dirichlet boundary conditions are approximated using a simple finite difference method corresponding to central differences and the addition of a high-order stabilization term…
The hidden shift problem is a natural place to look for new separations between classical and quantum models of computation. One advantage of this problem is its flexibility, since it can be defined for a whole range of functions and a…
We demonstrate when and how an entire left-infinite orbit of an underlying dynamical system or observations from such left-infinite orbits can be uniquely represented by a pair of elements in a different space, a phenomenon which we call…
In the past, Kepler painstakingly derived laws of planetary motion using difficult to understand and hard to follow techniques. In 1843 William Hamilton created and described the quaternions, which extend the complex numbers and can easily…