We numerically study the first-principle dynamics and thermostatistics of a d-dimensional classical inertial Heisenberg ferromagnetic model ( d = 1 , 2 , 3 ) with interactions decaying with the distance r i j as 1 / r i j α ( α ≥ 0 ), where the limit α = 0 ( α → ∞ ) corresponds to infinite-range (nearest-neighbour) interactions, and the ratio α / d > 1 ( 0 ≤ α / d ≤ 1 ) characterizes the short-ranged (long-ranged) regime. By means of first-principle molecular dynamics we study: (i) The scaling with the system size N of the maximum Lyapunov exponent λ in the form λ ∼ N - κ , where κ ( α / d ) depends only on the ratio α / d ; (ii) The time-averaged single-particle angular momenta probability distributions for a typical case in the long-range regime 0 ≤ α / d ≤ 1 (which turns out to be well fitted by q-Gaussians), and (iii) The time-averaged single-particle energies probability distributions for a typical case in the long-range regime 0 ≤ α / d ≤ 1 (which turns out to be well fitted by q-exponentials). Through the Lyapunov exponents we observe an intriguing, and possibly size-dependent, persistence of the non-Boltzmannian behavior even in the α / d > 1 regime. The universality that we observe for the probability distributions with regard to the ratio α / d makes this model similar to the α -XY and α -Fermi-Pasta-Ulam Hamiltonian models as well as to asymptotically scale-invariant growing networks.
Keywords: complex Hamiltonian systems; long-ranged-interacting thermostatistics; lyapunov exponents; nonextensive statistical mechanics.