A Solvable N-body Problem of Goldfish Type Featuring N2 Arbitrary Coupling Constants

A N-body problem “of goldfish type” is introduced, the Newtonian (“acceleration equal force”) equations of motion of which describe the motion of N pointlike unit-mass particles moving in the complex z-plane. The model—for arbitrary N—is solvable, namely its configuration (positions and velocities of the N “particles”) at any later time t can be obtained from its configuration at the initial time by algebraic operations. It features specific nonlinear velocity-dependent many-body forces depending on N² arbitrary (complex) coupling constants. Sufficient conditions on these constants are identified which cause the model to be isochronous—so that all its solutions are then periodic with a fixed period independent of the initial data. A variant with twice as many arbitrary coupling constants, or even more, is also identified.     

Solvable and/or Integrable and/or Linearizable N-Body Problems in Ordinary (Three-Dimensional) Space. I

Abstract

Several N -body problems in ordinary (3-dimensional) space are introduced which are characterized by Newtonian equations of motion (“acceleration equal force;” in most cases, the forces are velocity-dependent) and are amenable to exact treatment (“solvable” and/or “integrable” and/or “linearizable”). These equations of motion are always rotation-invariant, and sometimes translation-invariant as well. In many cases they are Hamiltonian, but the discussion of this aspect is postponed to a subsequent paper. We consider “few-body problems” (with, say, N =1,2,3,4,6,8,12,16,...) as well as “many-body problems” (N an arbitrary positive integer). The main focus of this paper is on various techniques to uncover such N -body problems. We do not discuss the detailed behavior of the solutions of all these problems, but we do identify several models whose motions are completely periodic or multiply periodic, and we exhibit in rather explicit form the solutions in some cases.     

Statistical properties of two-dimensional periodic Lorentz gas with infinite horizon

We study the asymptotic statistical behavior of the 2-dimensional periodic Lorentz gas with an infinite horizon. We consider a particle moving freely in the plane with elastic reflections from a periodic set of fixed convex scatterers. We assume that the initial position of the particle in the phase space is random with uniform distribution with respect to the Liouville measure of the periodic problem. We are interested in the asymptotic statistical behavior of the particle displacement in the plane as the timet goes to infinity. We assume that the particle horizon is infinite, which means that the length of free motion of the particle is unbounded. Then we show that under some natural assumptions on the free motion vector autocorrelation function, the limit distribution of the particle displacement in the plane is Gaussian, but the normalization factor is (t logt)^1/2 and nott ^1/2 as in the classical case. We find the covariance matrix of the limit distribution.     

A Hamiltonian yielding damped motion in an homogeneous magnetic field: quantum treatment

In earlier work, a Hamiltonian describing the classical motion of a particle moving in two dimensions under the combined influence of a perpendicular magnetic field and of a damping force proportional to the particle velocity, was indicated. Here we derive the quantum propagator for the Hamiltonian in different representations, one corresponding to momentum space, the other to position, and the third to a natural choice of “velocity” variables. We call attention to the following noteworthy fact: the Hamiltonian contains three parameters which do not in any way influence the motion of the position of the particle. However, at the quantum level, the propagator, even in the position representation, depends in an intricate way on these classically irrelevant parameters. This creates considerable doubt as to the validity of such a quantization procedure, as the physical results predicted differ for various Hamiltonians, all of which describe the dissipative dynamics equally well.     

PARTICULAR SOLITONS IN NONLINEAR LATTICE

One soliton of particle velocity and its amplitude (maximum particle velocity of one soliton) in Toda lattice is given analytically. It has also been known numerically that the maximum particle velocity (when the collision of two solitons reaches their maximum, we define V_n at this time as its maximum particle velocity) during the collision of two solitons moving in the same direction is equal to the difference between the amplitudes of two solitons if the difference is large enough; however, the maximum particle velocity is equal to the adding up of the amplitudes of two solitons moving in the opposite directions. The relationship between the maximum value of e^-(r_n)-1 and their initial amplitude of e^-(r_n)-1 is also given analytically in Toda lattice if the amplitudes of the two solitons are the same and their moving directions are opposite. Compared with the Boussinesq equation, there are differences between the Toda lattice equation and the Boussinesq equation for solitons during the collision.     

Modelling detonation of heterogeneous explosives with embedded inert particles using detonation shock dynamics: Normal and divergent propagation in regular and simplified microstructure

This paper discusses the mathematical formulation of Detonation Shock Dynamics (DSD) regarding a detonation shock wave passing over a series of inert spherical particles embedded in a high-explosive material. DSD provides an efficient method for studying detonation front propagation in such materials without the necessity of simulating the combustion equations for the entire system. We derive a series of partial differential equations in a cylindrical coordinate system and a moving shock-attached coordinate system which describes the propagation of detonation about a single particle, where the detonation obeys a linear shock normal velocity-curvature (D_n–κ) DSD relation. We solve these equations numerically and observe the short-term and long-term behaviour of the detonation shock wave as it passes over the particles. We discuss the shape of the perturbed shock wave and demonstrate the periodic and convergent behaviour obtained when detonation passes over a regular, periodic array of inert spherical particles.     

Some Stability and Instability Criteria for Ideal Plane Flows

We investigate stability and instability of steady ideal plane flows for an arbitrary bounded domain. First, we obtain some general criteria for linear and nonlinear stability. Second, we find a sufficient condition for the existence of a growing mode to the linearized equation. Third, we construct a steady flow which is nonlinearly and linearly stable in the L² norm of vorticity but linearly unstable in the L² norm of velocity.     

High-temperature brownian motors: Deterministic and stochastic fluctuations of a periodic potential

The high-temperature unidirectional motion of a Brownian particle with time-dependent potential energy described by a spatially asymmetric periodic function is considered. A general formula derived for the mean velocity ν of such a motion is specified for dichotomic deterministic and Markovian stochastic processes. In both cases, ν increases linearly for low-frequencies γ of potential-energy fluctuations and reaches maxima for γ about the inverse time of diffusion by the spatial period of the potential. The behaviors of ν for large γ values are different in these cases: ν ∝ γ^?2 and ν ts γ^?1 for the deterministic and stochastic processes, respectively. It is shown that the direction of the motor motion depends on the relative lifetimes of each of the dichotomic-process states if the amplitude of the potential-energy fluctuations is fairly large in comparison with the mean value.     

Diffusive Behavior for Randomly Kicked Newtonian Particles in a Spatially Periodic Medium

We prove a central limit theorem for the momentum distribution of a particle undergoing an unbiased spatially periodic random forcing at exponentially distributed times without friction. The start is a linear Boltzmann equation for the phase space density, where the average energy of the particle grows linearly in time. Rescaling time, the momentum converges to a Brownian motion, and the position is its time-integral showing superdiffusive scaling with time t ^3/2. The analysis has two parts: (1) to show that the particle spends most of its time at high energy, where the spatial environment is practically invisible; (2) to treat the low energy incursions where the motion is dominated by the deterministic force, with potential drift but where symmetry arguments cancel the ballistic behavior.     

Features of the resonant interaction of moving neutral atoms and clusters with superlattice surface

Features of the interaction of moving neutral atoms, molecules, and clusters with a superlattice field (for example, the system of linear magnetic and electric domains) are considered. It is shown that the character of the particle motion depends on the ratio of the frequency ω₂₁ of the internal electromagnetic resonance to the bounce frequency Ω _s determined by the superlattice period, the velocity of the particle motion, and the possible moments of the particle in the ground d ₁₁ and excited d ₂₂ states. The conditions for regimes of attraction and repulsion of particles by the superlattice are considered. The preconditions for formation of a one-dimensional potential well located far from the superlattice and for stable channeling of neutral and charged particles in this well are also considered. Depending on the ratio of ω₂₁ to Ω _s , particle sorting and beam separation occur during interaction of the multicomponent beam consisting of different particles with the superlattice field.     

A convenient expression of the time-derivative,of arbitrary order k,of the zero zn(t) of a time-dependent polynomial pN(z;t) of arbitrary degree N in z,and solvable dynamical systems

Let p_N (z; t) be a (monic) time-dependent polynomial of arbitrary degree N in z, and let z_n ≡ z_n (t) be its N zeros: . In this paper we report a convenient expression of the k-th time-derivative of the zero z_n (t). This formula plays a key role in the identification of classes of solvable dynamical systems describing the motion of point-particles moving in the complex z-plane while nonlinearly interacting among themselves; one such example, featuring many arbitrary parameters, is reported, including its variation describing the motion of many particles moving in the real Cartesian xy-plane and interacting among themselves via rotation-invariant Newtonian equations of motion (”accelerations equal forces”).     

Theory of Gravitational-Inertial Field of Universe

In this paper the basic proposition is a generalization of the metric tensor by introduction of an inertial field tensor satisfying ?_ig_lm ? g_lm;i ≠ 0. On the basis of variational equations a system of more general covariant equations of gravitational-inertial field is obtained. In Einstein's approximation these equations reduce to the field equations of Einstein. The solution of fundamental problems of generl taheory of relativity by means of the new equations give the same results as Einstein's equations. However application of these equations to the cosmologic problem leads to following results: 1. All Galaxies in the Universe (actually all bodies if gravitational attraction is not considered) “disperse” from each other according to Hubble's law. Thus contrary to Friedmann's theory (according to which the “expansion of Universe” began from the singular state with an infinite velocity) the velocity of “dispersion” of bodies begins from the zero value and in the limit tends to the velocity of light. 2. The “dispertion” of bodies represents a free motion in the inertial field and Hubble's law represents a law of motion of free bodies in the inertial field - the law of inertia. All critical systems (with Schwarzschild radius) are specific because they exist in maximal inertial and gravitational potentials. The Universe represents a critical system, it exists under the Schwarzschild radius. In the high-potential inertial and gravitational fields the material mass in a static state or in the process of motion with decelleration is subject to an inertial and gravitational “annihilation”. Under the maximal value of inertial and gravitational potentials (= c²) the material mass is completely “evaporated” transforming into a radiation mass. The latter is concentrated in the “horizon” of the critical system. All critical systems –“black holes”- represent geon systems, i.e., the local formations of gravitational-electromagnetic radiations, held together by their own gravitational and inertial fields. The Universe, being a critical system, is “wrapped” in a geon crown. The Universe is in a state of dynamical equilibrium. Near the external part of its boundary surface a transformation of matter into electromagnetic-gravitational-neutrineal energy (geon mass) takes place. Inside the Universe, in the galaxies takes place the synthesis of matter from geon mass, penetrating from the external part of the world (from geon crown) by means of a tunneling mechanism. The geon system may be considered as a natural entire cybernetic system.     

Chaotic diffusion and 1/f-noise of particles in two-dimensional solids

We investigate the classical nonlinear dynamics of a particle moving conservatively in a two-dimensional periodic potential. The particle exhibits diffusive motion in the absence of random forces. In a broad range of energies above the potential barrier, the diffusion process is anomalously accelerated and associated with 1/f-noise in the power spectrum of velocity fluctuations. The analysis of Poincaré surfaces of section and the distribution of free paths indicate that the phenomenon is caused by a trapping of orbits in a self-similar hierarchy of nested cantori. We describe a statistical theory for this mechanism in terms of a renewal process and a random walk on a hierarchical lattice.Work supported by Deutsche Forschungsgemeinschaft     

Possibility of electrical conduction in filled bands

It is shown that a moving neutral particle interacting with electrons may cause an “electron drag” within a filled band. The calculation uses perturbation theory and periodic boundary conditions and is based on the one-electron model. WithN being the number and ¯v the average velocity of the electrons, one finds that for largeN the electronic velocity sumN¯v induced by the motion of the neutral particle is independent ofN, i.e. of the size of the system. The lowest-order contributions toN¯v that do not necessarily vanish are seen to be those of second order in the interaction potential. These second-order contributions are studied. In a simple one-dimensional model they are found to be, in fact, not necessarily zero and to be proportional to the velocity of the neutral particle. An order-of-magnitude formula forN¯v is derived for this case. The calculation suggests that mobile neutral particles may act as charge carriers, their effective charge possibly being much smaller than the elementary charge. In real systems, neutral particles which interact with electrons might be represented by phonons and excitons.     

Approximate indexing of icosahedral polyhedra with fibonacci numbers

In general, indexing faces of icosahedral face-forms requires irrational numbers. However, for many practical purposes an approximate indexing based on triplets of integer numbers can be used. Two possible approaches called, respectively, “Fibonacci Matrix Methods” (FMM) and the “Linear Combination Method” (LCM) are described. FMM relies on the use of “auxiliary” matrices F_n, F² _n, F³ _n and F⁴ _n which have Fibonacci numbers as their elements. These matrices allow good approximation of the results usually obtained using the standard five-fold rotation matrices which are typical of icosahedral symmetry. LCM is based on the use of a classical crystallographic rule i.e. the so-called “Goldschmidt Complication Law” which is just a particular case of linear combination of triplets of face indices, with integers as coefficients. The occurrence of large integer indices is remarked.     

On the semiclassical approach to the ionization process during sputtering

The semiclassical “rate equation” approach to the ionization process during sputtering is shown to be the correct one in the limit of high temperatures. Specifically we show that simple “ rate equation” represents a high temperature asymptotic form of the equation of motion for the occupation number 〈n_a(t)〉 of relevant valence level of the sputtered particle.     

Cascades of Particles Moving at Finite Velocity in Hyperbolic Spaces

A branching process of particles moving at finite velocity over the geodesic lines of the hyperbolic space (Poincaré half-plane and Poincaré disk) is examined. Each particle can split into two particles only once at Poisson spaced times and deviates orthogonally when splitted. At time t, after N(t) Poisson events, there are N(t)+1 particles moving along different geodesic lines. We are able to obtain the exact expression of the mean hyperbolic distance of the center of mass of the cloud of particles. We derive such mean hyperbolic distance from two different and independent ways and we study the behavior of the relevant expression as t increases and for different values of the parameters c (hyperbolic velocity of motion) and λ (rate of reproduction). The mean hyperbolic distance of each moving particle is also examined and a useful representation, as the distance of a randomly stopped particle moving over the main geodesic line, is presented.     

振动筛系统的两类余维三分岔与非常规混沌演化

建立了振动筛系统的动力学模型,推导出了其周期运动的Poincaré 映射,基于Poincaré 映射方法着重研究了系统Flip-Hopf-Hopf余维三分岔、三次强共振条件下的Hopf-Hopf余维三分岔以及三种非常规的混沌演化过程.研究结果表明,此两类余维三分岔点附近的动力学行为变得更加复杂和新颖,在分岔点附近出现了三角形吸引子、3T²环面分岔以及“五角星型”、“轮胎型”概周期吸引子,揭示了环面爆破、环面倍化以及T²环面分岔向混沌演化的过程,这些结果对于振动筛系统的动力学优化设计提供了理论参考. 关键词： 余维三分岔 非常规混沌演化 T²环面分岔')" href="#">T²环面分岔     

Superdiffusion and stable laws

The superdiffusion equation with a fractional Laplacian Δ ^α/2 in _N-dimensional space describes the asymptotic (t→∞) behavior of a generalized Poisson process with the range (discontinuity) distribution density ∼|x|−α−1. The solutions of this equation belong to a class of spherically symmetric stable distributions. The main properties of these solutions are given together with their representations in the form of integrals and series and the results of numerical calculations. It is shown that allowance for the finite velocity of free particle motion for α>1 merely amounts to a reduction in the diffusion coefficient with the form of the distribution remaining stable. For α<1 the situation changes radically: the expansion velocity of the diffusion packet exceeds the velocity of free particle motion and the superdiffusion equation becomes physically meaningless. Zh. éksp. Teor. Fiz. 115, 1411–1425 (April 1999)     

Propagation and stability of kinks in driven and damped nonlinear klein-gordon chains1

We consider the propagation of kinks in an elastic chain in a bistable or multistable potential under the action of a driving force [M. Büttiker and H. Thomas,Phys. Rev. A 37:235 (1988)]. Each element of the chain is subject to a damping force proportional to its velocity. We show that both the propagation velocity of the kinks as a function of the driving field, and the kink width as a function of propagation velocity, are determined by characteristic functions which depend only on the form of the potential. These functions can be found by considering a single particle moving in the upside-down potential of the chain. The general properties of these functions are studied and illustrated by several examples. The stability of these driven kinks is discussed. Interestingly, we find in addition to the expected discrete localized eigenmodes a two-dimensional continuum of oscillatory modes with a localized envelope.This paper will appear in a forthcoming issue of theJournal of Statistical Physics.