Next Article in Journal
Top-N Recommender Systems Using Genetic Algorithm-Based Visual-Clustering Methods
Next Article in Special Issue
Cubature Formulas of Multivariate Polynomials Arising from Symmetric Orbit Functions
Previous Article in Journal
Parity-Time Symmetry and the Toy Models of Gain-Loss Dynamics near the Real Kato’s Exceptional Points
Previous Article in Special Issue
On a Reduction Formula for a Kind of Double q-Integrals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants

by
Francesco Calogero
1,2
1
Physics Department, University of Rome “La Sapienza”, 00185 Rome, Italy
2
Istituto Nazionale di Fisica Nucleare, Sezione di Roma, 00185 Rome, Italy
Symmetry 2016, 8(7), 53; https://doi.org/10.3390/sym8070053
Submission received: 31 March 2016 / Revised: 16 June 2016 / Accepted: 17 June 2016 / Published: 24 June 2016
(This article belongs to the Special Issue Symmetry in Orthogonal Polynomials)

Abstract

:
Three new classes of N-body problems of goldfish type are identified, with N an arbitrary positive integer ( N 2 ). These models are characterized by nonlinear Newtonian (“accelerations equal forces”) equations of motion describing N equal point-particles moving in the complex z-plane. These highly nonlinear equations feature many arbitrary coupling constants, yet they can be solved by algebraic operations. Some of these N-body problems are isochronous, their generic solutions being all completely periodic with an overall period T independent of the initial data (but quite a few of these solutions are actually periodic with smaller periods T / p with p a positive integer); other models are isochronous for an open region of initial data, while the motions for other initial data are not periodic, featuring instead scattering phenomena with some of the particles incoming from, or escaping to, infinity in the remote past or future.

1. Introduction

Recently, a new technique to identify many-body problems solvable by algebraic operations has been introduced [1,2], and several examples of such models have been discussed [1,2,3,4,5,6]. In the present paper, three new classes of such models are introduced and discussed (the impatient reader may immediately glance at these findings reported in the next section). We term these models “of goldfish type” because their Newtonian (“accelerations equal forces”) equations of motion read as follows:
z ¨ n = 2 r + 1 i ω z ˙ n + r r + 1 ω 2 z n - = 1 , n N z ˙ n + i r ω z n z ˙ + i r ω z z n - z + f n z , z · ,
with the nonlinear functions f n z , z · appropriately defined (see below). (The original goldfish model is the special case of these equations of motion with r = 0 and f n z , z · = 0 ; after its first identification as a solvable model [7], and its tentative recognition as a “goldfish” [8], this N-body problem and some of its extensions have been investigated in several publications (see, for instance, [9,10,11,12,13,14,15,16,17,18,19,20]).
Notation 1.1. 
Above, and hereafter, N is an arbitrary positive integer ( N 2 ); indices such as n , m , run over the positive integers from 1 to N (unless otherwise indicated: see, for instance, the restriction on the values of in (1)); the N complex coordinates z n z n t identify the positions of N points moving in the complex z-plane as functions of the time t; below, we also introduce other complex variables ζ n ζ n τ that depend on the complex variable τ related to the real variable t (“time”) as follows:
τ τ t = exp i ω t - 1 i ω .
Above, and hereafter, i is the imaginary unit ( i 2 = - 1 ), ω is an arbitrary (nonvanishing) real constant, so that τ t vanishes at the initial time t = 0 and is periodic in time with period T 0 ,
τ 0 = 0 , τ t + T 0 = τ t , T 0 = 2 π ω .
Above, and hereafter, we adopt the standard notation according to which superimposed dots denote differentiations with respect to the time variable t, so that, for instance,
τ ˙ τ ˙ t = exp i ω t ,
while primes appended to a function indicate differentiations with respect to the argument of that function (see, for instance, (3b) and (3c), where the relevant variable is of course τ). The auxiliary complex coordinates ζ n ζ n τ are related to the N complex coordinates z n z n t —which are the main protagonists of our N-body problems, as they identify the positions of N equal unit-mass point-particles moving in the complex z-plane according to the Newtonian (“accelerations equal forces”) equations of motion of type (1)—by the following relations:
z n t = exp i r ω t ζ n τ , ζ n τ = exp - i r ω t z n t ,
implying
z ˙ n = i r ω z n + exp i r + 1 ω t ζ n (3b) ζ n = exp - i r + 1 ω t z ˙ n - i r ω z n ,
z ¨ n = i 2 r + 1 ω z ˙ n + r r + 1 ω 2 z n + exp i r + 2 ω t ζ n (3c) ζ n = exp - i r + 2 ω t z ¨ n - i 2 r + 1 ω z ˙ n - r r + 1 ω 2 z n .
Above, and hereafter, the number r is required to be real and rational,
r = q p ,
where, of course, the numerator q and the denominator p are coprime integers, and, for definiteness, we hereafter assume that p is positive, p 1 .
Finally, note that above we often omitted indicating explicitly the argument of functions, and we will do so throughout whenever this can be done without causing confusion.
In this paper, a key role is played by the following monic τ-dependent polynomial of degree N in the (complex) variable ζ,
p N ζ ; ζ ̲ τ ; γ τ = n = 1 N ζ - ζ n τ = ζ N + m = 1 N γ m τ ζ N - m ,
which features the N zeros ζ n τ and the N coefficients γ m τ ; and by its counterpart,
P N z ; z ̲ t ; w t = n = 1 N z - z n t = z N + m = 1 N w m t z N - m ,
featuring the N zeros z n t and the N coefficients w m t (see Notation 1.1).
Notation 1.2. 
The N coefficients γ m γ m τ of the polynomial (4a) are of course expressed in terms of its N zeros ζ n ζ n τ via the formulas
γ m = - 1 m σ m ζ ̲ ,
σ m ζ ̲ = 1 m ! n 1 , n 2 , . . . , n m = 1 N * ζ n 1 ζ n 2 · · · ζ n m .
Above, and hereafter, the symbol ζ ̲ denotes the unordered set of the N complex numbers ζ n (and likewise z ̲ is the unordered set of the N complex numbers z n ), while γ (see (4a) and below) is the N-vector with components γ m (and likewise w is the N-vector with components w m ). The symbol n 1 , n 2 , . . . , n m N * denotes the sum from 1 to N over the m indices n 1 , n 2 , . . . , n m with the restriction that these indices be all different among themselves, so that for N = 2 ,
γ 1 = - ζ 1 + ζ 2 , γ 2 = ζ 1 ζ 2 ,
for N = 3 ,
γ 1 = - ζ 1 + ζ 2 + ζ 2 , γ 2 = ζ 1 ζ 2 + ζ 2 ζ 3 + ζ 3 ζ 1 , γ 3 = - ζ 1 ζ 2 ζ 3 ,
and so on. Let us also display the corresponding formulas for the τ-derivative of γ m γ m τ ,
γ m = - 1 m n = 1 N ζ n σ n , m ζ n ; ζ ̲ ,
σ n , 1 ζ n ; ζ ̲ = 1 , σ n , m ζ n ; ζ ̲ = 1 m - 1 ! n 1 , n 2 , . . . , n m - 1 = 1 N * n * ζ n 1 ζ n 2 · · · ζ n m - 1 , (6b) m = 2 , . . . , N ,
where the symbol n 1 , n 2 , . . . , n m N * n * denotes the sum from 1 to N over the m indices n 1 , n 2 , . . . , n m with the restriction that these indices be all different among themselves and all different from the index n, so that for N = 2 ,
γ 1 = - ζ 1 + ζ 2 , γ 2 = ζ 1 ζ 2 + ζ 1 ζ 2 ,
for N = 3 ,
γ 1 = - ζ 1 + ζ 2 + ζ 2 , γ 2 = ζ 1 ζ 2 + ζ 3 + ζ 2 ζ 3 + ζ 1 + ζ 3 ζ 1 + ζ 2 , (6d) γ 3 = - ζ 1 ζ 2 ζ 3 + ζ 2 ζ 3 ζ 1 + ζ 3 ζ 1 ζ 2 ,
and so on. Analogous formulas hold, of course, for the relations among the N coefficients w m w m t and the N zeros z n z n t of the polynomial (4b), and their time-derivatives; for instance,
w m = - 1 m m ! n 1 , n 2 , . . . , n m = 1 N * z n 1 z n 2 · · · z n m .
Let us also report the relations—implied by (3a)—among the N coefficients w m t of the polynomial (4b) and the N coefficients γ m τ of the polynomial (4a)
w m t = exp i m r ω t γ m τ ,
with τ related to t via (2) and r defined by (3d), implying
w ˙ m t = i m r ω w m t + exp i m r + 1 ω t γ m τ , (8b) γ m τ = exp - i m r + 1 ω t w ˙ m t - i m r ω w m t .
w ¨ m t = i 2 m r + 1 ω w ˙ m t + m r m r + 1 ω 2 w m t + exp i m r + 2 ω t γ m τ , γ m τ = exp - i m r + 2 ω t w ¨ m t (8c) - i 2 m r + 1 ω w ˙ m t - m r m r + 1 ω 2 w m t .
Remark 1.1. 
Two comments on the problem to determine the N zeros ζ n , respectively, z n of a monic polynomial of degree N in ζ, respectively, z from its N coefficients γ m , respectively, w m (see (4a), respectively, (4b)).
(i)
Of course, the assignment of the N coefficients of a polynomial defines uniquely the corresponding unordered set of its N zeros, but generally it only allows to compute explicitly these N zeros for N 4 .
(ii)
Moreover—and quite relevantly in our context (see below)—if a polynomial features a dependence on an additional variable (as, for instance, the dependence of the polynomial P N z ; z ̲ t ; w t on the real variable t ("time") (see (4b)), then the unordered character of the set of its N zeros z n t is generally only relevant at one value of time, say at the "initial" time t = 0 , since, at other values of time, the ordering gets generally determined by the natural requirement that the functions z n t evolve continuously over time. This prescription then fixes, for all time, the ordering of the zeros z n t —i.e., the assignment of the value n of its index to each zero z n t —as long as the coefficients w m t evolve themselves continuously and unambiguously over time and moreover no "collision" of two or more zeros occurs over the time evolution, i.e., for all time z n t z t if n (since clearly such collisions imply a loss of identity of the coinciding zeros). However, this identification requires an analysis of the time evolution of the N zeros z n t not only in the complex z-plane, but in fact over the N-sheeted Riemann surface associated to the N roots of the polynomial P N z ; z ̲ t ; w t , and/or over the evolution of each coefficient w m t if its time evolution takes itself place on a Riemann surface (as it indeed happens in the cases discussed below). ■
A key formula for the identification and investigation of solvable Newtonian N-body problems is the following relation [1] among the τ-evolutions of the N zeros and the N coefficients of the monic polynomial p N ζ ; τ (see 4a)) of degree N in its argument ζ and depending on the extra variable τ:
ζ n = = 1 ; n N 2 ζ n ζ ζ n - ζ - = 1 ; n N ζ n - ζ - 1 m = 1 N γ m ζ n N - m .
In the present paper—which is a sequel to [6], so that we dispense below from some of the remarks reported there about the significance of the solvable many-body models treated herein—we focus on the N-body models that obtain via (9) for the coordinates ζ n τ —and especially for the corresponding coordinates z n t (see (3))—when the coefficients γ m τ evolve according to the following system of decoupled nonlinear Ordinary Differential Equations (ODEs):
γ m = g m γ m a m γ m b m ,
where the N "coupling constants" g m are arbitrary (possibly complex), while the two (sets of) parameters a m and b m are instead hereafter required to be real rational numbers and to be related to each other (and to the parameter r , (see (3d)) by the single relation:
a m = 2 - 1 + b m m r 1 + m r , b m = 1 - a m + 2 - a m m r ,
which—provided neither r nor 1 + m r vanish, i.e., r 0 and r - 1 / m for m in its range from 1 to N, as we now assume (the case r = 0 shall be discussed below separately)—is necessary and sufficient to guarantee that the corresponding system of ODEs satisfied by the coefficients w m t be autonomous, reading as follows (as the diligent reader will easily verify via (10a) and (8)):
w ¨ m t = i 2 m r + 1 ω w ˙ m t + m r m r + 1 ω 2 w m t (10c) + g m w ˙ m - i m r ω w m a m w m b m .
In fact, we will focus below only on certain specific assignments of the parameters r , a m , b m which allow the explicit solution of the ODEs (10a) in terms of elementary functions. The corresponding solvable N-body models satisfied by the coordinates z n z n t are displayed—and their properties discussed—in the following Section 2, with the corresponding proofs provided in Section 3, while the special case with a m = 2 and b m = - 1 (and with g m an arbitrary rational number) is treated in [6].

2. Results

The Newtonian equations of motion of the first class of N-body models treated in this paper read as follows:
z ¨ n = i 2 r + 1 ω z ˙ n + r r + 1 ω 2 z n + = 1 , n N 2 z ˙ n - i r ω z n z ˙ - i r ω z z n - z - = 1 , n N z n - z - 1 (11a) · m = 1 N g m w ˙ m - i m r w m 1 + 1 / 1 + m r z n N - m ,
with r an arbitrary rational number (see (3d); of course r - 1 / m for m = 1 , 2 , . . . , N if g m 0 ), ω an arbitrary nonvanishing real number, the N coupling constants g m arbitrary complex numbers (not all vanishing) and (see Notation 1.2)
w m = - 1 m m ! n 1 , n 2 , . . . , n m = 1 N * z n 1 z n 2 · · · z n m ,
w ˙ m = - 1 m m - 1 ! n = 1 N z ˙ n n 1 , n 2 , . . . , n m - 1 = 1 N * n * z n 1 z n 2 · · · z n m .
We then assert that this N-body model is solvable by algebraic operations, its remarkable properties being detailed by the following.
Proposition 2.1. 
The N complex coordinates z n t providing the solution at time t of the initial-values problem of the Newtonian equations of motion (11) are the N zeros of the monic polynomial (4b), the N coefficients w m t of which are given, in terms of the initial data z n 0 , z ˙ n 0 , by the following formulas:
w m t = exp i r m ω t γ m τ t , τ t = exp i ω t - 1 i ω ,
γ m τ = w m 0 + B m 1 1 - τ τ ˇ m 1 - m r - B m 1 ,
τ ˇ m 1 = m r A m w m 0 + B m 1 - 1 / m r A m ,
A m = g m m r 1 + m r
B m 1 = w ˙ 0 - i m r ω w m 0 m r / 1 + m r A m - w m 0 ,
w m 0 = - 1 m m ! n 1 , n 2 , . . . , n m = 1 N * z n 1 0   z n 2 0 · · · z n m 0 ,
w ˙ m 0 = - 1 m m - 1 ! n = 1 N z ˙ n 0 n 1 , n 2 , . . . , n m - 1 = 1 N * n * z n 1 0   z n 2 0 · · · z n m - 1 0 .
These functions w m t are, of course, defined by continuity in t from their initial values w m 0 (i.e., this prescription identifies the determination of all the rational roots appearing in the above formulas), and note that, if g m = 0 , the formula (12b) must be replaced by
γ m τ = w m 0 + w ˙ 0 - i m r ω w m 0 τ .
Let us also report the properties of these solutions, (12), when they are generated by generic initial data, z n 0 , z ˙ n 0 , excluding the nongeneric initial data identified in Proposition 2.2 for which the system of evolution equations (11) runs into a singularity at a finite time.
These solutions are all nonsingular for all time, remaining in a finite region—the size of which depends on the initial data—of the complex z-plane and featuring no particle collisions, i.e., for all time t, z n t z t if n ; and they are all completely periodic with a period T which is an integer multiple of the basic period T 0 (see (2b)),
z n t + T = z n t , T = K T 0 ,
with the positive integer K restricted as follows: p K p N ! , where, of course, p is the denominator of the rational number r (see (3d)). ■
Note the arbitrariness of the real (nonvanishing) rational number r, and especially of the N complex coupling constants g m . In addition, note that Proposition 2.1 implies that all generic solutions of the N-body model characterized by the Newtonian equations of motion (11a)—excluding the nongeneric solutions which are singular (see below Proposition 2.2)—are completely periodic with the same period T M A X = p N ! T 0 . However, there are lots of solutions that are completely periodic with periods which are integer submultiples of T M A X . The detailed identification of these solutions and their periods is a nontrivial matter, as shown, for instance, by the discussion of this phenomenology in the paper [11]—that treats the "periodic goldfish model" (for this terminology, see [8]), which is in fact characterized by the same equations of motions (11a), but with all coupling constants vanishing, g m = 0 —and by the detailed investigation of the structure of the Riemann surfaces associated with other analogous many-body models [21,22,23,24,25].
Proposition 2.2. 
The solutions of the Newtonian equations of motion (11) may feature singularities via two phenomena, both of which correspond to nongeneric initial data.
The first phenomenon is characterized by initial data satisfying—for at least one value m ¯ of the index m in its range from 1 to N such that g m ¯ 0 and - m ¯ r is not a positive integer—the equality
1 + i ω τ ˇ m ¯ 1 = 1 ,
where τ ˇ m ¯ 1 is defined in terms of the initial data z n 0 , z ˙ n 0 as above (see (12c)).
This singularity occurs at the time t = t ˇ 1 defined as follows: t ˇ 1 = min t ˇ m ¯ 1 (with the minimum taken over the values of the index m ¯ satisfying the condition (14)), with
t ˇ m ¯ 1 = θ m ¯ 1 ω , t ˇ m ¯ 1 > 0 ,
where θ m ¯ 1 is defined m o d 2 π as follows:
1 + i ω τ ˇ m ¯ 1 = exp i θ m ¯ 1 .
Note that the fact that θ m ¯ 1 is real is implied by (14), and that certainly 0 < t ˇ m ¯ 1 < T 0 , hence 0 < t ˇ 1 < T 0 .
The second phenomenon causing the equations of motion (11a) to hit a singularity is the occurrence of a collision of two (or possibly more) particles at some time t = t c ; so that, for some indices n , with n (of course both in their interval from 1 to N) there holds the equality z n t c = z t c , causing the term z n t - z t appearing in the denominators in the right-hand side of (11a) to vanish. This phenomenon corresponds to the fact that two of the zeros of the polynomial (4b) with (12) coincide, a fact that clearly only happens for nongeneric initial data (in the complex z-plane), although the condition on the initial data that would cause this phenomenon to happen—and the time at which this phenomenon would happen—can be computed explicitly only for small values of N (and even then the result is not very enlightening). ■
Let us complete our discussion of the first class of solvable N-body problems by displaying the equations of motion (11) in the (simplest) N = 2 case:
z ¨ n = i 2 r + 1 ω z ˙ n + r r + 1 ω 2 z n - - 1 n z 1 - z 2 - 1 · 2 z ˙ 1 - i r z 1 z ˙ 2 - i r z 2 + g 1 z ˙ 1 + z ˙ 2 - i r ω z 1 + z 2 1 + 1 / 1 + r z n (16) - g 2 z ˙ 1 z 2 + z ˙ 2 z 1 - 2 i r ω z 1 z 2 1 + 1 / 1 + 2 r , n = 1 , 2 ,
with r an arbitrary rational number ( r 0 , r - 1 / 2 , r - 1 ), and g 1 , g 2 two arbitrary complex numbers (not both vanishing). Note that, for r = - 2 / 3 , both exponents on the right-hand side of these ODEs are integers, since then 1 / 1 + r = 3 , 1 / 1 + 2 r = - 3 , and this is as well the case for r = - 3 / 4 , since then 1 / 1 + r = 4 , 1 / 1 + 2 r = - 2 .
Let us end our treatment of the first class of solvable N-body problems characterized by the Newtonian equations of motion (11) by pointing out that, for r = - 2 / 3 and N = 3 , the exponents on the right-hand side of (11a) are all integers, indeed for r = - 2 / 3 , the three exponents 2 + m r / 1 + m r take, for m = 1 , 2, 3, the three values 4, - 2 , 0.
The Newtonian equations of motion of the second class of solvable N-body models read as follows:
z ¨ n = i ω z ˙ n + = 1 , n N 2 z ˙ n z ˙ z n - z (17) - = 1 , n N z n - z - 1 m = 1 N g m w ˙ m 2 z n N - m ,
with the various quantities defined as above (in particular w ˙ m defined in terms of z ˙ n and z n by (11c) with Notation 1.2). Note that this might be considered the special case of the first class of models (see above) with r = 0 , a m = 2 and b m = 0 , which was previously excluded because it requires a special treatment (see Section 3).
The solvable character of this N-body system is demonstrated by the following
Proposition 2.3. 
The N complex coordinates z n t providing the solution at time t of the initial-values problem of the Newtonian equations of motion (17) are the N zeros of the monic polynomial (4b), the N coefficients w m t of which are given, in terms of the initial data z n 0 , z ˙ n 0 , by the following formula:
w m t = w m 0 - 1 g m ln 1 - τ t τ ˇ m 2 , τ ˇ m 2 = 1 g m w ˙ m 0 ,
of course with τ t = exp i ω t - 1 / i ω (see (2a)). Here and below, the function ln 1 - τ t / τ ˇ m 2 is defined by continuity in t from its vanishing value at t = 0 where τ t = 0 , and, of course, w m 0 , respectively, w ˙ m 0 are defined in terms of the initial data z n 0 and z ˙ n 0 by (12f), respectively, (12g), and if g m = 0 then (18) becomes w m t = w m 0 + w ˙ m 0 τ t , while if w ˙ m 0 = 0 it yields w m t = w m 0 . ■
Clearly, these coefficients w m t are completely periodic in the time t with period T 0 , see (2a),
w m t + T 0 = w m t ,
iff the initial data satisfy the inequality
1 + i ω τ ˇ m 2 > 1 ,
iff instead the initial data satisfy the opposite inequality,
1 + i ω τ ˇ m 2 < 1 ,
they are periodic in t except for a constant shift (independent of the initial data!) over each period T 0 , so that
w m t + k T 0 = w m t - 2 π i k g m , k = ± 1 , ± 2 , ± 3 , . . . ,
hence, in this second case, they diverge as t ± . In addition—in the intermediate, nongeneric case in which the initial data imply the equality
1 + i ω τ ˇ m 2 = 1
—the coefficient w m t diverges at the finite times
t ˇ m 2 = i ω - 1 ln 1 + i ω τ ˇ m 2 mod T 0 ,
which are, of course, real thanks to (21a).
Correspondingly, the particles coordinates z n t —being the N zeros of the polynomial (4b) with the coefficients w m t (see (18))—are periodic with period K T 0 K being a positive integer in the range from 1 to N ! (see below Remark 3.1)—iff the initial data satisfy the inequality (19b), while iff instead the initial data satisfy the opposite inequality (20a) at least one of the particle coordinates z n t comes from or escapes to infinity in the remote past and future: see, for instance, the relevant discussion in Appendix G (“Asymptotic behavior of the zeros of a polynomial whose coefficients diverge exponentially”) of the book [9]. In addition, of course, if the (nongeneric) initial data imply validity of the equality (21a), the equations of motion run into a singularity at t = t ˇ m 2 (see (21b)). Other nongeneric initial data causing the equations of motion (17) to run into a singularity at a finite time are those leading to particle collisions. Note that generally the nongeneric initial data causing singularities are also those that separate the regions of initial data associated to different behaviors of the model, including the emergence of the higher periodicities associated to values of K larger than unity as well as the periodic and nonperiodic cases.
We do not display explicitly the equations of motion of this second model in the simplest N = 2 case because they can be immediately obtained by setting r = 0 in those of the first model (see (16)).
We conclude our report on the properties of the second class of solvable N-body problems with the following remark (proven in Section 3):
Remark 2.1. 
The system of Newtonian equations of motion (17) is Hamiltonian, albeit with a time-dependent Hamiltonian. ■
The third class of solvable N-body problems is characterized by the Newtonian equations of motion
z ¨ n = i 2 r + 1 ω z ˙ n + r r + 1 ω 2 z n + = 1 , n N 2 z ˙ n - i r ω z n z ˙ - i r ω z z n - z - = 1 , n N z n - z - 1 (22) · m = 1 N g m w ˙ m - i m r w m 2 - 1 1 + m r w m - 1 + 1 m r z n N - m ,
where ω is an arbitrary nonvanishing real number, r is an arbitrary rational number (of course if g m 0 then r 0 and r 1 / m for the integer m in its range from 1 to N), the N coupling constants g m are arbitrary complex numbers (not all vanishing), and the quantities w m w m t , respectively, w ˙ m w ˙ m t are expressed in terms of the coordinates z n z n t and z ˙ n z ˙ n t via (11b), respectively (11c).
The solvable character of this N-body problem is demonstrated by the following
Proposition 2.4. 
The N complex coordinates z n t providing the solution at time t of the initial-values problem of the Newtonian equations of motion (22) are the N zeros of the monic polynomial (4b), the N coefficients w m t of which are given, in terms of the initial data z n 0 , z ˙ n 0 , by the following formulas:
w m t = exp i m r ω t γ m τ t , τ t = exp i ω t - 1 i ω ,
γ m τ = B m 3 m r τ ˇ m 3 - τ ^ m 3 τ ˇ m 3 - τ 1 m r - 1 - m r ,
B m 3 = w ˙ m 0 - i m r ω w 0 1 1 + m r A m - w 0 1 m r , A m = g m m r 1 + m r ,
τ ˇ m 3 = - 1 + B m 3 w 0 - 1 m r - m r A m 1 + m r B m 3 ,
τ ^ m 3 = τ ˇ m 3 1 - 1 + B m 3 w 0 - 1 m r m r ,
where w m 0 , respectively, w ˙ m 0 , are defined in terms of the initial data z n 0 and z ˙ n 0 by (12f), respectively, (12g), and of course the determinations of the rational roots are implied by the requirement that these formulas be valid at the initial time t = 0 and thereafter by continuity in t , and again if g m = 0 , the expression (23b) must be replaced by (12h). ■
It is plain that, if neither m r nor 1 / m r are integers, for generic initial data the function γ m τ features, in the complex τ-plane, two rational branch points at τ = τ ˇ m 3 , respectively, at τ = τ ^ m 3 , and that the periodicity of γ m t γ m τ t as function of t (see (23b)), is determined by the location in the complex τ-plane of these two branch points with respect to the circle C ˜ centered at τ = i / ω and with radius 1 / ω on which rotates the point τ t = exp i ω t - 1 / i ω as function of the time t. If both branch points are located outside the circle C ˜ —and the condition on the initial data determining this outcome is clearly validity of both inequalities
1 + i ω τ ˇ m 3 > 1 , 1 + i ω τ ^ m 3 > 1 ,
with τ ˇ m 3 , respectively, τ ^ m 3 defined in terms of the initial data by (23d), respectively, (23e)—then clearly γ m t γ m τ t as a function of t is periodic with period T 0 , see (2a),
γ m t + T 0 = γ m t ,
hence w m t is periodic in t with period T p ^ m = p ^ m T 0 ,
w m t + T p ^ m = w m t ,
with p ^ m = p / m if this number is integer; otherwise, p ^ m = p (where of course p is the denominator of r: see (23a) and (3d)). In addition, as a consequence—if the inequalities (24a) are valid for all values of m in the range from 1 to N—then the polynomial (4b) is periodic in t with period
T = T 0 MinimumCommonMultiple m ¯ = 1 , 2 , . . . , N p ^ m ¯ ,
where the MinimumCommonMultiple must be evaluated for the values of m = m ¯ such that g m ¯ 0 ; hence (see below Remark 3.1), the N coordinates z n t are periodic with period Q T , where Q is a positive integer in the range from 1 to N ! .
It is also plain that the solutions z n t are periodic in t with a period which is a positive integer multiple Q ˜ of T 0 even if the initial data imply instead that some of the inequalities (24a) are reversed, but, in these cases, the determination of the outcome—in particular, of the value of Q ˜ —requires, to begin with, a standard analysis of the structure of the Riemann surface associated to the function
φ m τ = τ ˇ m 3 - τ ^ m 3 τ ˇ m 3 - τ 1 m r - 1 - m r
of the complex variable τ , and consequently of the periodicity of φ m t φ m τ t as a function of the real variable t ("time") when τ t = exp i ω t - 1 / i ω so that τ t travels on the circle C ˜ with radius 1 / ω centered at τ = i / ω ; and then an analysis of the resulting periodicity of the solutions z n t analogous to that made in Section 3 for the first class of solvable models. We leave this task to the interested reader.
Let us complete this discussion of the third class of solvable N-body problems by displaying the equations of motion (22) in the (simplest) N = 2 case:
z ¨ n = i 2 r + 1 ω z ˙ n + r r + 1 ω 2 z n - - 1 n z 1 - z 2 - 1 · 2 z ˙ 1 - i r z 1 z ˙ 2 - i r z 2 - g 1 - z ˙ 1 - z ˙ 2 + i r z 1 + z 2 2 - 1 1 + r - z 1 + z 2 - 1 + 1 r z n (26) - g 2 z ˙ 1 z 2 + z 1 z ˙ 2 - 2 i r z 1 z 2 2 - 1 1 + 2 r z 1 z 2 - 1 + 1 2 r , n = 1 , 2 ,
with r an arbitrary rational number ( r 0 , r - 1 / 2 , r - 1 ), and g 1 , g 2 two arbitrary complex numbers (not both vanishing).
In addition, let us end this discussion of the third class of solvable N-body problems characterized by the Newtonian equations of motion (22) by also displaying these equations in the following two special cases:
case (i), with r = - 1 / 2 and g m = 0 for m > 1 , when these equations read
z ¨ n = - 1 4 ω 2 z n + = 1 , n N 2 z ˙ n + i ω z n 2 z ˙ + i ω z 2 z n - z (27a) + g 1 z n N - 1 = 1 , n N z n - z - 1 n 1 = 1 N z n 1 - 3 ;
case (ii), with r = - 1 / 4 and g m = 0 for m 2 , when these equations of motion read
z ¨ n = 1 2 i ω z ˙ n - 3 16 ω 2 z n + = 1 , n N 4 z ˙ n + i ω z n 4 z ˙ + i ω z 8 z n - z (27b) - 1 2 g 2 z n N - 2 = 1 , n N z n - z - 1 n 1 , n 2 = 1 ; n 1 n 2 N z n 1 z n 2 - 3 .

3. Proofs

In this Section 3 we prove the results reported without their proofs in the preceding Section 2.
The first task is to integrate once the ODE (10a). This is an easy task, yielding
γ m τ = A m γ m τ 1 + b m + B m 1 2 - a m ,
A m = g m 2 - a m 1 + b m = g m m r 1 + m r ,
B m = γ m 0 2 - a m A m - γ m 0 1 + b m .
Note that, above and hereafter, we exclude from consideration the special case with a m = 2 and b m = - 1 , which is treated in [6].
The next integration can be performed in terms of elementary functions only for special assignments (satisfying the restriction (10b)) of the parameters a m and b m , to which we restrict attention in the present paper.
Our first assignment is
b m = 0 ,
implying via (10b)
a m = 2 - m r 1 + m r = 2 + m r 1 + m r
and
γ m τ = A m γ m τ + B m 1 1 + 1 m r ,
with A m , respectively, B m 1 defined by (28b) (or, equivalently, (12d)), respectively, (12e). Here, we assume, of course, that r 0 (implying a m 2 ) and moreover that r - 1 / m with m = 1 , 2 , . . . , N (implying that a m is a finite rational number for all values of m in its range from 1 to N).
This ODE (29c) can now be easily integrated, yielding (12b) with (12c).
These developments clearly prove the first part of Proposition 2.1.
To prove the second part of Proposition 2.1 we ascertain, to begin with, the periodicity properties as functions of the time variable t of the coefficients γ m t γ m τ t (see (12b)). The starting point is the observation that τ t —see (2) or (12a)—is a periodic function of t with period T 0 , rotating in the complex τ-plane on the circle C ˜ centered at the point i / ω and having radius 1 / ω . Hence, any holomorphic function of τ is as well periodic in t with period T 0 ; this clearly is (for all values of m in its range from 1 to N) the case of the functions γ m τ (see (12b)), if r is a negative integer. If instead r is a positive integer, the functions γ m τ are meromorphic in τ, featuring a polar singularity at τ = τ ˇ m 1 (see (12b)). In this case, γ m t γ m τ t is again generally periodic in t with period T 0 , γ m t + T 0 = γ m t , but for the nongeneric assignments of the initial data such that τ ˇ m 1 falls on the circle C ˜ —note that τ ˇ m 1 does depend on the initial data (see (12c) and (12e)), and that the condition for this to happen is validity of the equality
1 + i ω τ ˇ m 1 = 1 ,
then the function γ m t γ m τ t diverges at the times
t ˇ m 1 = i ω - 1 ln 1 + i ω τ ˇ m 1 mod T 0 .
Finally, if r is rational but not integer, i.e., p > 1 (see (3b)), and m is not an integer multiple of p, γ m τ features a rational branch point at τ = τ ˇ m 1 (see (12b)), then the evolution of γ m t γ m τ t as function of the time t depends on the location of the branch point τ ˇ m 1 in the complex τ-plane, whether it falls outside, inside or just on the circle C ˜ . The latter case requires that the initial data satisfy the condition (30a), implying again that they are not generic. While clearly the condition that the branch point be located outside the circle C ˜ is validity of the inequality
1 + i ω τ ˇ m 1 > 1 ,
and for the corresponding nongeneric initial data the function γ m t γ m τ t is again periodic with period T 0 . If instead the branch point falls inside the circle C ˜ , and the condition on the initial data for this to happen is validity of the opposite inequality
1 + i ω τ ˇ m 1 < 1 ,
the periodicity of γ m t γ m τ t gets modified: the period is then T ˜ p = p T 0 (see (12b) and (3d)), unless p / m = p ˇ m is an integer, in which case the period is T ˇ p ˇ m = p ˇ m T 0 .
Next, let us discuss the t-periodicity of w t , related to γ t γ τ t by (8a) or (12a). We then note that the prefactor exp i r m ω t in (8a) or (12a) is also periodic in t (see (3d)), with the same period T ˇ p ˇ m = p ˇ m T 0 , where p ˇ m = p unless p / m is an integer, in which case p ˇ m = p / m . We may therefore conclude that, for all generic initial data, the functions w m t are periodic with period T ˇ p ˇ m ,
w m t + T ˇ p ˇ m = w m t .
Next, let us discuss the periodicity of the zeros z n t of the polynomial (4b) with coefficients w m t periodic as indicated just above. In this context, the following Remark 3.1 is relevant.
Remark 3.1. 
If a time-dependent polynomial P N z ; t , of degree N in z, is time-periodic with period T ˜ , P N z ; t + T ˜ = P N z ; t , the unordered set z ̲ t of its N zeros z n t is of course periodic with the same period T ˜ , z ̲ t + T ˜ = z ̲ t (since after a period the polynomial is unchanged); however, due to the possibility that these zeros, as it were, "exchange their places" over their time evolution, the period of each individual zero z n t , considered as a continuous function of time, may be a positive integer multiple of T ˜ ; although of course that multiple cannot exceed the number N ! of permutations of the N elements of the unordered set z ̲ t (for a detailed discussion of this phenomenology in analogous many-body contexts see [11,21,22,23,24,25]). ■
We can therefore conclude that the N zeros z n t of the monic polynomial (4b), the N coefficients w m t of which are periodic as described above, are also periodic, for generic initial data, with periods K T 0 , with the positive integer K restricted as follows:
p K N ! p ,
since the very definition of p ˇ m (see above) implies that p is the Minimum Common Multiple of the N parameters p ˇ m .
This ends the proofs of the findings reported in Propositions 2.1 and 2.2.
The proofs of the other two Propositions reported in Section 2 follow, but below we omit the proofs of some aspects of the results reported in Section 2, which are analogous to those detailed above in the formulations and proofs of Propositions 2.1 and 2.2.
Hence, the only aspect relevant to prove Proposition 2.3 that we do report is the derivation of the formula (18). The starting point is the ODE
γ m = g m γ m 2 ,
which is the special case of (10a) with r = 0 , a m = 2 and b m = 0 . The integration of this ODE is a trivial task, yielding
γ m τ = γ m 0 - ln 1 - g m γ m 0 τ g m ,
and the fact that this implies (18) is an immediate consequence of (8) with r = 0 .
As for the proof of Remark 2.1, it is based on the observation that the decoupled nonlinear system of N ODEs
w ¨ m = i ω w ˙ m + g w ˙ m 2
is Hamiltonian, since it is implied by the standard Hamiltonian equations of motion
w ˙ = H v , w ; t v m , v ˙ = - H v , w ; w m ,
with the (time-dependent) Hamiltonian
H v , w ; t = exp i ω t v α exp g α w ,
as the diligent reader will easily verify. Note the arbitrariness of the (nonvanishing!) parameter α.
However, via the two identities
z ˙ n = - = 1 , n N z n - z - 1 m = 1 N w ˙ z n N - m ,
z ¨ n = = 1 , n N 2 z ˙ n z ˙ z n - z - = 1 , n N z n - z - 1 m = 1 N w ¨ z n N - m ,
linking—as proven in [1]—the time evolutions of the N zeros z n t and of the N coefficients w m t of the monic time-dependent polynomial (4b), it is easily seen that the set of ODEs (36a) correspond precisely to the equations of motion (17). This implies that the equations of motion (17) are also Hamiltonian, since the N coordinates z n are linked to the N canonical coordinates w m by the "point" transformations—not involving the canonical momenta v m —which relate the coefficients and the zeros of a polynomial (see (7) or (11b)).
Let us finally deal with the third class of N-body models, obtained by identifying another set of parameters r, a m , b m allowing the ODEs (28) to be explicitly integrated in terms of elementary functions. To this end, we introduce the auxiliary functions
η m τ = γ m τ s m , γ m τ = η m τ 1 / s m ,
with the option to assign the parameters s m at our convenience (see below). We then see, from (28) with (10b), that these auxiliary functions satisfy the ODEs
η m = s m A m η α m + B m η β m 1 2 - a m ,
α m = b m + 1 + s m - 1 2 - a m s m , β m = s m - 1 2 - a m s m ,
with A m , respectively, B m defined by (28b), respectively, (28c). Thus, we set α m = 0 and β m = 1 , implying (after a bit of elementary algebra)
s m = - 1 m r , b m b m 3 = - 1 + 1 m r , a m a m 3 = 2 - 1 1 + m r .
The ODE (38b) then reads
η m = - 1 m r A m 1 + B m 3 η m 1 + m r ,
with A m , respectively, B m 3 defined by (28b) respectively (28c) with (39), while the assignments of the parameters r and g m are still free (of course r 0 , r - 1 / m ). This ODE, (40), can then be immediately integrated, yielding
η m τ = 1 B m 3 + η m 0 1 - τ τ ˇ m 3 - 1 m r - 1 B m 3 ,
B m 3 = γ 0 1 1 + m r A m - γ 0 1 m r ,
τ ˇ m 3 = - 1 + B m 3 η m 0 - m r A m 1 + m r B m 3 .
Hence, finally, via (38a) with (39) and (8), we arrive at the formulas (23), thereby proving Proposition 2.4.

4. Outlook

In this Section 4 we outline tersely possible extensions of the results obtained in this paper. A tentative list of such possible developments is provided in the last Section of [6]: the first of those reported there is essentially implemented in the present paper, while the others are still open. Also open is the possibility to investigate a combination of the three classes of models treated above, and of the class of models treated in [6], by assuming that, for different values of the parameter m, some coefficients γ m evolve according to the three different solvable subcases of the ODEs (28) discussed above in Section 3 or to the case discussed in [6]. In addition, of course, also open is the investigation of the hierarchies of solvable models associated with each of those mentioned herein via the notion of generations of monic polynomials such that the coefficients of the polynomials of the next generation coincide with the zeros of a polynomial of the current generation [2].
Other N-body problems—more general than those treated in [6] and in the present paper—can be investigated by analogous techniques to those employed above, while allowing each γ m t to satisfy a solvable second-order ODE still belonging to the class (10) but forsaking the restriction on the parameters a m and b m necessary and sufficient to allow the solution of these ODEs in terms of elementary functions (for instance, a natural extension might include special functions of elliptic and hyperelliptic type).

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Calogero, F. New solvable variants of the goldfish many-body problem. Studies Appl. Math. 2016. [Google Scholar] [CrossRef]
  2. Bihun, O.; Calogero, F. Generations of monic polynomials such that the coefficients of each polynomial of the next generation coincide with the zeros of a polynomial of the current generation, and new solvable many-body problems. Lett. Math. Phys. 2016, 106, 1011–1031. [Google Scholar] [CrossRef]
  3. Bihun, O.; Calogero, F. A new solvable many-body problem of goldfish type. J. Nonlinear Math. Phys. 2016, 23, 28–46. [Google Scholar] [CrossRef]
  4. Bihun, O.; Calogero, F. Novel solvable many-body problems. J. Nonlinear Math. Phys. 2016, 23, 190–212. [Google Scholar] [CrossRef]
  5. Calogero, F. A solvable N-body problem of goldfish type featuring N2 arbitrary coupling constants. J. Nonlinear Math. Phys. 2016, 23, 300–305. [Google Scholar] [CrossRef]
  6. Calogero, F. Novel isochronous N-body problems featuring N arbitrary rational coupling constants. J. Math. Phys. 2016, in press. [Google Scholar]
  7. Calogero, F. Motion of poles and zeros of special solutions of nonlinear and linear partial differential equations, and related “solvable” many-body problems. II Nuovo Cimento B 1978, 43, 177–241. [Google Scholar] [CrossRef]
  8. Calogero, F. The neatest many-body problem amenable to exact treatments (a “goldfish”?). Phys. D Nonlinear Phenom. 2001, 152–153, 78–84. [Google Scholar] [CrossRef]
  9. Calogero, F. Classical many-body problems amenable to exact treatments; Springer: Heidelberg, Germany, 2001; p. 749. [Google Scholar]
  10. Nucci, M.C. Calogero’s ’goldfish’ is indeed a school of free particles. J. Phys. A Math. Gen. 2004, 37, 11391–11400. [Google Scholar] [CrossRef]
  11. Gómez-Ullate, D.; Sommacal, M. Periods of the goldfish many-body problem. J. Nonlinear Math. Phys. 2005, 12, 351–362. [Google Scholar] [CrossRef]
  12. Suris, Y.B. Time discretization of F. Calogero’s “Goldfish” System. J. Nonlinear Math. Phys. 2005, 12, 633–647. [Google Scholar] [CrossRef]
  13. Calogero, F.; Iona, S. Novel solvable extensions of the goldfish many-body model. J. Math. Phys. 2005, 46, 103515. [Google Scholar] [CrossRef]
  14. Guillot, A. The Painlevé property for quasihomogeneous systems and a many-body problem in the plane. Comm. Math. Phys. 2005, 256, 181–194. [Google Scholar] [CrossRef]
  15. Bruschi, M.; Calogero, F. Novel solvable variants of the goldfish many-body model. J. Math. Phys. 2006, 47, 022703. [Google Scholar] [CrossRef]
  16. Calogero, F.; Langmann, E. Goldfishing by gauge theory. J. Math. Phys. 2006, 47, 082702. [Google Scholar] [CrossRef]
  17. Arlind, J.; Bordemann, M.; Hoppe, J.; Lee, C. Goldfish geodesics and Hamiltonian reduction of matrix dynamics. Lett. Math. Phys. 2008, 84, 89–98. [Google Scholar] [CrossRef]
  18. Calogero, F. Isochronous Systems; Oxford University Press: Oxford, UK, 2008; p. 264. [Google Scholar]
  19. Bihun, O.; Calogero, F. Solvable many-body models of goldfish type with one-, two- and three-body forces. SIGMA 2013, 9, 059. [Google Scholar] [CrossRef]
  20. Jairuk, U.; Yoo-Kong, S.M.; Tanasittikosol, M. The Lagrangian structure of Calogero’s goldfish model. Theor. Math. Phys. 2015, 183, 665–683. [Google Scholar] [CrossRef]
  21. Calogero, F.; Gómez-Ullate, D.; Santini, P.M.; Sommacal, M. The transition from regular to irregular motions, explained as travel on Riemann surfaces. J. Phys. A Math. Gen. 2005, 38, 8873–8896. [Google Scholar] [CrossRef]
  22. Sommacal, M. The Transition from Regular to Irregular Motion, Explained as Travel on Riemann Surfaces. Ph.D. Thesis, SISSA, Trieste, Italy, 2005. [Google Scholar]
  23. Grinevich, P.G.; Santini, P.M. Newtonian dynamics in the plane corresponding to straight and cyclic motions on the hyperelliptic curve: Ergodicity, isochrony and fractals. Phys. D Nonlinear Phenom. 2007, 232, 22–32. [Google Scholar] [CrossRef]
  24. Calogero, F.; Gómez-Ullate, D.; Santini, P.M.; Sommacal, M. Towards a theory of chaos explained as travel on Riemann surfaces. J. Phys. A. Math. Theor. 2009, 42, 015205. [Google Scholar] [CrossRef]
  25. Gómez-Ullate, D.; Santini, P.M.; Sommacal, M.; Calogero, F. Understanding complex dynamics by means of an associated Riemann surface. Phys. D Nonlinear Phenom. 2012, 241, 1291–1305. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Calogero, F. Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants. Symmetry 2016, 8, 53. https://doi.org/10.3390/sym8070053

AMA Style

Calogero F. Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants. Symmetry. 2016; 8(7):53. https://doi.org/10.3390/sym8070053

Chicago/Turabian Style

Calogero, Francesco. 2016. "Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants" Symmetry 8, no. 7: 53. https://doi.org/10.3390/sym8070053

APA Style

Calogero, F. (2016). Three New Classes of Solvable N-Body Problems of Goldfish Type with Many Arbitrary Coupling Constants. Symmetry, 8(7), 53. https://doi.org/10.3390/sym8070053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop