Open AccessArticle

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

Francesco Calogero

^1,2

Physics Department, University of Rome “La Sapienza”, 00185 Rome, Italy

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)

Download Versions Notes

Abstract

Three new classes of N-body problems of goldfish type are identified, with N an arbitrary positive integer (

N \geq 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.

Keywords:

solvable many-body problems; integrable many-body problems; integrable dynamical systems; solvable many-body problems of goldfish type

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:

\begin{matrix} {\ddot{z}}_{n} = (2 r + 1) i ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} \\ - \sum_{ℓ = 1, ℓ \neq n}^{N} [\frac{({\dot{z}}_{n} + i r ω z_{n}) ({\dot{z}}_{ℓ} + i r ω z_{ℓ})}{z_{n} - z_{ℓ}}] + f_{n} (\vec{z}, \overset{\cdot}{\vec{z}}), \end{matrix}

(1)

with the nonlinear functions

f_{n} (\vec{z}, \overset{\cdot}{\vec{z}})

appropriately defined (see below). (The original goldfish model is the special case of these equations of motion with

r = 0

and

f_{n} (\vec{z}, \overset{\cdot}{\vec{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 \geq 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} \equiv 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} \equiv ζ_{n} (τ)

that depend on the complex variable τ related to the real variable t (“time”) as follows:

τ \equiv τ (t) = \frac{exp (i ω t) - 1}{i ω} .

(2a)

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} = \frac{2 π}{|ω|} .

(2b)

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,

\dot{τ} \equiv \dot{τ} (t) = exp (i ω t),

(2c)

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} \equiv ζ_{n} (τ)

are related to the N complex coordinates

z_{n} \equiv 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),

(3a)

implying

\begin{matrix} {\dot{z}}_{n} & = & i r ω z_{n} + exp [i (r + 1) ω t] ζ_{n}^{'} \\ (3b) & ζ_{n}^{'} & = & exp [- i (r + 1) ω t] ({\dot{z}}_{n} - i r ω z_{n}), \end{matrix}

\begin{matrix} {\ddot{z}}_{n} & = & i (2 r + 1) ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} + exp [i (r + 2) ω t] ζ_{n}^{″} \\ (3c) & ζ_{n}^{″} & = & exp [- i (r + 2) ω t] [{\ddot{z}}_{n} - i (2 r + 1) ω {\dot{z}}_{n} - r (r + 1) ω^{2} z_{n}] . \end{matrix}

Above, and hereafter, the number r is required to be real and rational,

r = \frac{q}{p},

(3d)

where, of course, the numerator q and the denominator p are coprime integers, and, for definiteness, we hereafter assume that p is positive,

p \geq 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} (ζ; \underset{̲}{ζ} (τ); \vec{γ} (τ)) = \prod_{n = 1}^{N} [ζ - ζ_{n} (τ)] = ζ^{N} + \sum_{m = 1}^{N} [γ_{m} (τ) ζ^{N - m}],

(4a)

which features the N zeros

ζ_{n} (τ)

and the N coefficients

γ_{m} (τ)

; and by its counterpart,

P_{N} (z; \underset{̲}{z} (t); \vec{w} (t)) = \prod_{n = 1}^{N} [z - z_{n} (t)] = z^{N} + \sum_{m = 1}^{N} [w_{m} (t) z^{N - m}],

(4b)

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} \equiv γ_{m} (τ)

of the polynomial (4a) are of course expressed in terms of its N zeros

ζ_{n} \equiv ζ_{n} (τ)

via the formulas

γ_{m} = {(- 1)}^{m} σ_{m} (\underset{̲}{ζ}),

(5a)

σ_{m} (\underset{̲}{ζ}) = \frac{1}{m!} \sum_{n_{1}, n_{2}, . . ., n_{m} = 1}^{N}^{*} (ζ_{n_{1}} ζ_{n_{2}} \cdot \cdot \cdot ζ_{n_{m}}) .

(5b)

Above, and hereafter, the symbol

\underset{̲}{ζ}

denotes the unordered set of the N complex numbers

ζ_{n}

(and likewise

\underset{̲}{z}

is the unordered set of the N complex numbers

z_{n}

), while

\vec{γ}

(see (4a) and below) is the N-vector with components

γ_{m}

(and likewise

\vec{w}

is the N-vector with components

w_{m}

). The symbol

\sum_{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},

(5c)

for

N = 3

γ_{1} = - (ζ_{1} + ζ_{2} + ζ_{2}), γ_{2} = ζ_{1} ζ_{2} + ζ_{2} ζ_{3} + ζ_{3} ζ_{1}, γ_{3} = - ζ_{1} ζ_{2} ζ_{3},

(5d)

and so on. Let us also display the corresponding formulas for the τ-derivative of

γ_{m} \equiv γ_{m} (τ)

γ_{m}^{'} = {(- 1)}^{m} \sum_{n = 1}^{N} [ζ_{n}^{'} σ_{n, m} (ζ_{n}; \underset{̲}{ζ})],

(6a)

\begin{matrix} σ_{n, 1} (ζ_{n}; \underset{̲}{ζ}) = 1, \\ σ_{n, m} (ζ_{n}; \underset{̲}{ζ}) = \frac{1}{(m - 1)!} \sum_{n_{1}, n_{2}, . . ., n_{m - 1} = 1}^{N}^{* n *} (ζ_{n_{1}} ζ_{n_{2}} \cdot \cdot \cdot ζ_{n_{m - 1}}), \\ (6b) & m = 2, . . ., N, \end{matrix}

where the symbol

\sum_{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}^{'},

(6c)

for

N = 3

\begin{matrix} γ_{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}), \end{matrix}

and so on. Analogous formulas hold, of course, for the relations among the N coefficients

w_{m} \equiv w_{m} (t)

and the N zeros

z_{n} \equiv z_{n} (t)

of the polynomial (4b), and their time-derivatives; for instance,

w_{m} = \frac{{(- 1)}^{m}}{m!} \sum_{n_{1}, n_{2}, . . ., n_{m} = 1}^{N}^{*} (z_{n_{1}} z_{n_{2}} \cdot \cdot \cdot z_{n_{m}}) .

(7)

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} (τ),

(8a)

with τ related to t via (2) and r defined by (3d), implying

\begin{matrix} {\dot{w}}_{m} (t) = i m r ω w_{m} (t) + exp [i (m r + 1) ω t] γ_{m}^{'} (τ), \\ (8b) & γ_{m}^{'} (τ) = exp [- i (m r + 1) ω t] [{\dot{w}}_{m} (t) - i m r ω w_{m} (t)] . \end{matrix}

\begin{matrix} {\ddot{w}}_{m} (t) = i (2 m r + 1) ω {\dot{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] [{\ddot{w}}_{m} (t) \\ (8c) & - i (2 m r + 1) ω {\dot{w}}_{m} (t) - m r (m r + 1) ω^{2} w_{m} (t)] . ■ \end{matrix}

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 \leq 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; \underset{̲}{z} (t); \vec{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) \neq z_{ℓ} (t)$ if $n \neq ℓ$ (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; \underset{̲}{z} (t); \vec{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}^{″} = \sum_{ℓ = 1; ℓ \neq n}^{N} (\frac{2 ζ_{n}^{'} ζ_{ℓ}^{'}}{ζ_{n} - ζ_{ℓ}}) - [\prod_{ℓ = 1; ℓ \neq n}^{N} {(ζ_{n} - ζ_{ℓ})}^{- 1}] \sum_{m = 1}^{N} [γ_{m}^{″} {(ζ_{n})}^{N - m}] .

(9)

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}},

(10a)

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 - \frac{(1 + b_{m}) m r}{1 + m r}, b_{m} = 1 - a_{m} + \frac{2 - a_{m}}{m r},

(10b)

which—provided neither r nor

1 + m r

vanish, i.e.,

r \neq 0

and

r \neq - 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)):

\begin{matrix} {\ddot{w}}_{m} (t) & = & i (2 m r + 1) ω {\dot{w}}_{m} (t) + m r (m r + 1) ω^{2} w_{m} (t) \\ (10c) & + g_{m} {({\dot{w}}_{m} - i m r ω w_{m})}^{a_{m}} {(w_{m})}^{b_{m}} . \end{matrix}

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} \equiv 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:

\begin{matrix} {\ddot{z}}_{n} = i (2 r + 1) ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} \\ + \sum_{ℓ = 1, ℓ \neq n}^{N} [\frac{2 ({\dot{z}}_{n} - i r ω z_{n}) ({\dot{z}}_{ℓ} - i r ω z_{ℓ})}{z_{n} - z_{ℓ}}] - [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] \\ (11a) & \cdot \sum_{m = 1}^{N} [g_{m} {({\dot{w}}_{m} - i m r w_{m})}^{1 + [1 / (1 + m r)]} {(z_{n})}^{N - m}], \end{matrix}

with r an arbitrary rational number (see (3d); of course

r \neq - 1 / m

for

m = 1, 2, . . ., N

g_{m} \neq 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} = \frac{{(- 1)}^{m}}{m!} \sum_{n_{1}, n_{2}, . . ., n_{m} = 1}^{N}^{*} (z_{n_{1}} z_{n_{2}} \cdot \cdot \cdot z_{n_{m}}),

(11b)

{\dot{w}}_{m} = \frac{{(- 1)}^{m}}{(m - 1)!} \sum_{n = 1}^{N} \{{\dot{z}}_{n} [\sum_{n_{1}, n_{2}, . . ., n_{m - 1} = 1}^{N}^{* n *} (z_{n_{1}} z_{n_{2}} \cdot \cdot \cdot z_{n_{m}})]\} .

(11c)

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),

{\dot{z}}_{n} (0),

by the following formulas:

w_{m} (t) = exp (i r m ω t) γ_{m} (τ (t)), τ (t) = \frac{exp (i ω t) - 1}{i ω},

(12a)

γ_{m} (τ) = [w_{m} (0) + B_{m}^{(1)}] {(1 - \frac{τ}{{\overset{ˇ}{τ}}_{m}^{(1)}})}^{- m r} - B_{m}^{(1)},

(12b)

{\overset{ˇ}{τ}}_{m}^{(1)} = \frac{m r {\{A_{m} [w_{m} (0) + B_{m}^{(1)}]\}}^{- 1 / (m r)}}{A_{m}},

(12c)

A_{m} = \frac{g_{m} m r}{1 + m r}

(12d)

B_{m}^{(1)} = \frac{{[\dot{w} (0) - i m r ω w_{m} (0)]}^{(m r) / (1 + m r)}}{A_{m}} - w_{m} (0),

(12e)

w_{m} (0) = \frac{{(- 1)}^{m}}{m!} \sum_{n_{1}, n_{2}, . . ., n_{m} = 1}^{N}^{*} [z_{n_{1}} (0) z_{n_{2}} (0) \cdot \cdot \cdot z_{n_{m}} (0)],

(12f)

{\dot{w}}_{m} (0) = \frac{{(- 1)}^{m}}{(m - 1)!} \sum_{n = 1}^{N} \{{\dot{z}}_{n} (0) \sum_{n_{1}, n_{2}, . . ., n_{m - 1} = 1}^{N}^{* n *} [z_{n_{1}} (0) z_{n_{2}} (0) \cdot \cdot \cdot z_{n_{m - 1}} (0)]\} .

(12g)

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) + [\dot{w} (0) - i m r ω w_{m} (0)] τ .

(12h)

Let us also report the properties of these solutions, (12), when they are generated by generic initial data,

z_{n} (0),

{\dot{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) \neq z_{ℓ} (t)

n \neq ℓ

; 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},

(13)

with the positive integer K restricted as follows:

p \leq K \leq 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

\bar{m}

of the index m in its range from 1 to N such that

g_{\bar{m}} \neq 0

and

- \bar{m} r

is not a positive integer—the equality

|1 + i ω {\overset{ˇ}{τ}}_{\bar{m}}^{(1)}| = 1,

(14)

where

{\overset{ˇ}{τ}}_{\bar{m}}^{(1)}

is defined in terms of the initial data

z_{n} (0)

{\dot{z}}_{n} (0)

as above (see (12c)).

This singularity occurs at the time

t = {\overset{ˇ}{t}}^{(1)}

defined as follows:

{\overset{ˇ}{t}}^{(1)} = min [{\overset{ˇ}{t}}_{\bar{m}}^{(1)}]

(with the minimum taken over the values of the index

\bar{m}

satisfying the condition (14)), with

{\overset{ˇ}{t}}_{\bar{m}}^{(1)} = \frac{θ_{\bar{m}}^{(1)}}{ω}, {\overset{ˇ}{t}}_{\bar{m}}^{(1)} > 0,

(15a)

where

θ_{\bar{m}}^{(1)}

is defined

m o d [2 π]

as follows:

1 + i ω {\overset{ˇ}{τ}}_{\bar{m}}^{(1)} = exp (i θ_{\bar{m}}^{(1)}) .

(15b)

Note that the fact that

θ_{\bar{m}}^{(1)}

is real is implied by (14), and that certainly

0 < {\overset{ˇ}{t}}_{\bar{m}}^{(1)} < T_{0}

, hence

0 < {\overset{ˇ}{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 \neq ℓ

(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:

\begin{matrix} {\ddot{z}}_{n} = i (2 r + 1) ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} - {(- 1)}^{n} {(z_{1} - z_{2})}^{- 1} \\ \cdot \{2 ({\dot{z}}_{1} - i r z_{1}) ({\dot{z}}_{2} - i r z_{2}) \\ + (g_{1} {[{\dot{z}}_{1} + {\dot{z}}_{2} - i r ω (z_{1} + z_{2})]}^{1 + [1 / (1 + r)]} z_{n} \\ (16) & - g_{2} {[{\dot{z}}_{1} z_{2} + {\dot{z}}_{2} z_{1} - 2 i r ω z_{1} z_{2}]}^{1 + [1 / (1 + 2 r)]})\}, n = 1, 2, \end{matrix}

with r an arbitrary rational number (

r \neq 0,

r \neq - 1 / 2,

r \neq - 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:

\begin{matrix} {\ddot{z}}_{n} = i ω {\dot{z}}_{n} + \sum_{ℓ = 1, ℓ \neq n}^{N} (\frac{2 {\dot{z}}_{n} {\dot{z}}_{ℓ}}{z_{n} - z_{ℓ}}) \\ (17) & - [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] \sum_{m = 1}^{N} [g_{m} {({\dot{w}}_{m})}^{2} {(z_{n})}^{N - m}], \end{matrix}

with the various quantities defined as above (in particular

{\dot{w}}_{m}

defined in terms of

{\dot{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),

{\dot{z}}_{n} (0)

, by the following formula:

w_{m} (t) = w_{m} (0) - \frac{1}{g_{m}} ln [1 - \frac{τ (t)}{{\overset{ˇ}{τ}}_{m}^{(2)}}], {\overset{ˇ}{τ}}_{m}^{(2)} = \frac{1}{g_{m} {\dot{w}}_{m} (0)},

(18)

of course with

τ (t) = [exp (i ω t) - 1] / (i ω)

(see (2a)). Here and below, the function

ln \{1 - [τ (t) / {\overset{ˇ}{τ}}_{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,

{\dot{w}}_{m} (0)

are defined in terms of the initial data

z_{n} (0)

and

{\dot{z}}_{n} (0)

by (12f), respectively, (12g), and if

g_{m} = 0

then (18) becomes

w_{m} (t) = w_{m} (0) + {\dot{w}}_{m} (0) τ (t)

, while if

{\dot{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),

(19a)

iff the initial data satisfy the inequality

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(2)}| > 1,

(19b)

iff instead the initial data satisfy the opposite inequality,

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(2)}| < 1,

(20a)

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) - \frac{2 π i k}{g_{m}}, k = \pm 1, \pm 2, \pm 3, . . .,

(20b)

hence, in this second case, they diverge as

t \to \pm \infty

. In addition—in the intermediate, nongeneric case in which the initial data imply the equality

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(2)}| = 1

(21a)

—the coefficient

w_{m} (t)

diverges at the finite times

{\overset{ˇ}{t}}_{m}^{(2)} = {(i ω)}^{- 1} ln (1 + i ω {\overset{ˇ}{τ}}_{m}^{(2)}) mod (T_{0}),

(21b)

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 = {\overset{ˇ}{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

\begin{matrix} {\ddot{z}}_{n} = i (2 r + 1) ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} \\ + \sum_{ℓ = 1, ℓ \neq n}^{N} [\frac{2 ({\dot{z}}_{n} - i r ω z_{n}) ({\dot{z}}_{ℓ} - i r ω z_{ℓ})}{z_{n} - z_{ℓ}}] - [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] \\ (22) & \cdot \sum_{m = 1}^{N} [g_{m} {({\dot{w}}_{m} - i m r w_{m})}^{2 - \frac{1}{1 + m r}} {(w_{m})}^{- 1 + \frac{1}{m r}} {(z_{n})}^{N - m}], \end{matrix}

where ω is an arbitrary nonvanishing real number, r is an arbitrary rational number (of course if

g_{m} \neq 0

then

r \neq 0

and

r \neq 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} \equiv w_{m} (t),

respectively,

{\dot{w}}_{m} \equiv {\dot{w}}_{m} (t)

are expressed in terms of the coordinates

z_{n} \equiv z_{n} (t)

and

{\dot{z}}_{n} \equiv {\dot{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),

{\dot{z}}_{n} (0)

, by the following formulas:

w_{m} (t) = exp (i m r ω t) γ_{m} (τ (t)), τ (t) = \frac{exp (i ω t) - 1}{i ω},

(23a)

γ_{m} (τ) = {[B_{m}^{(3)}]}^{m r} {\{{[\frac{{\overset{ˇ}{τ}}_{m}^{(3)} - {\hat{τ}}_{m}^{(3)}}{{\overset{ˇ}{τ}}_{m}^{(3)} - τ}]}^{\frac{1}{m r}} - 1\}}^{- m r},

(23b)

B_{m}^{(3)} = \frac{{[{\dot{w}}_{m} (0) - i m r ω w (0)]}^{\frac{1}{1 + m r}}}{A_{m}} - {[w (0)]}^{\frac{1}{m r}}, A_{m} = \frac{g_{m} m r}{1 + m r},

(23c)

{\overset{ˇ}{τ}}_{m}^{(3)} = - \frac{{\{1 + B_{m}^{(3)} {[w (0)]}^{- \frac{1}{m r}}\}}^{- m r}}{{(A_{m})}^{1 + m r} B_{m}^{(3)}},

(23d)

{\hat{τ}}_{m}^{(3)} = {\overset{ˇ}{τ}}_{m}^{(3)} \{1 - {(1 + B_{m}^{(3)} {[w (0)]}^{- \frac{1}{m r}})}^{m r}\},

(23e)

where

w_{m} (0),

respectively,

{\dot{w}}_{m} (0),

are defined in terms of the initial data

z_{n} (0)

and

{\dot{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

τ = {\overset{ˇ}{τ}}_{m}^{(3)}

, respectively, at

τ = {\hat{τ}}_{m}^{(3)},

and that the periodicity of

γ_{m} (t) \equiv γ_{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

\tilde{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

\tilde{C}

—and the condition on the initial data determining this outcome is clearly validity of both inequalities

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(3)}| > 1, |1 + i ω {\hat{τ}}_{m}^{(3)}| > 1,

(24a)

with

{\overset{ˇ}{τ}}_{m}^{(3)}

, respectively,

{\hat{τ}}_{m}^{(3)}

defined in terms of the initial data by (23d), respectively, (23e)—then clearly

γ_{m} (t) \equiv γ_{m} (τ (t))

as a function of t is periodic with period

T_{0}

, see (2a),

γ_{m} (t + T_{0}) = γ_{m} (t),

(24b)

hence

w_{m} (t)

is periodic in t with period

T_{{\hat{p}}_{m}} = {\hat{p}}_{m} T_{0},

w_{m} (t + T_{{\hat{p}}_{m}}) = w_{m} (t),

(24c)

with

{\hat{p}}_{m} = p / m

if this number is integer; otherwise,

{\hat{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} \{\underset{\bar{m} = 1, 2, . . ., N}{MinimumCommonMultiple} [{\hat{p}}_{\bar{m}}]\},

(24d)

where the MinimumCommonMultiple must be evaluated for the values of

m = \bar{m}

such that

g_{\bar{m}} \neq 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

\tilde{Q}

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

\tilde{Q}

—requires, to begin with, a standard analysis of the structure of the Riemann surface associated to the function

φ_{m} (τ) = {\{{[\frac{{\overset{ˇ}{τ}}_{m}^{(3)} - {\hat{τ}}_{m}^{(3)}}{{\overset{ˇ}{τ}}_{m}^{(3)} - τ}]}^{\frac{1}{m r}} - 1\}}^{- m r}

(25)

of the complex variable

τ,

and consequently of the periodicity of

φ_{m} (t) \equiv φ_{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

\tilde{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:

\begin{matrix} {\ddot{z}}_{n} = i (2 r + 1) ω {\dot{z}}_{n} + r (r + 1) ω^{2} z_{n} - {(- 1)}^{n} {[z_{1} - z_{2}]}^{- 1} \\ \cdot \{2 ({\dot{z}}_{1} - i r z_{1}) ({\dot{z}}_{2} - i r z_{2}) \\ - g_{1} {[- {\dot{z}}_{1} - {\dot{z}}_{2} + i r (z_{1} + z_{2})]}^{2 - \frac{1}{1 + r}} {[- (z_{1} + z_{2})]}^{- 1 + \frac{1}{r}} z_{n} \\ (26) & - g_{2} {({\dot{z}}_{1} z_{2} + z_{1} {\dot{z}}_{2} - 2 i r z_{1} z_{2})}^{2 - \frac{1}{1 + 2 r}} {(z_{1} z_{2})}^{- 1 + \frac{1}{2 r}}\}, n = 1, 2, \end{matrix}

with r an arbitrary rational number (

r \neq 0, r \neq - 1 / 2,

r \neq - 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

\begin{matrix} {\ddot{z}}_{n} = - \frac{1}{4} ω^{2} z_{n} + \sum_{ℓ = 1, ℓ \neq n}^{N} [\frac{(2 {\dot{z}}_{n} + i ω z_{n}) (2 {\dot{z}}_{ℓ} + i ω z_{ℓ})}{2 (z_{n} - z_{ℓ})}] \\ (27a) & + g_{1} {(z_{n})}^{N - 1} [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] {[\sum_{n_{1} = 1}^{N} (z_{n_{1}})]}^{- 3}; \end{matrix}

case (ii), with

r = - 1 / 4

and

g_{m} = 0

for

m \neq 2

, when these equations of motion read

\begin{matrix} {\ddot{z}}_{n} = \frac{1}{2} i ω {\dot{z}}_{n} - \frac{3}{16} ω^{2} z_{n} + \sum_{ℓ = 1, ℓ \neq n}^{N} [\frac{(4 {\dot{z}}_{n} + i ω z_{n}) (4 {\dot{z}}_{ℓ} + i ω z_{ℓ})}{8 (z_{n} - z_{ℓ})}] \\ (27b) & - \frac{1}{2} g_{2} {(z_{n})}^{N - 2} [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] {[\sum_{n_{1}, n_{2} = 1; n_{1} \neq n_{2}}^{N} (z_{n_{1}} z_{n_{2}})]}^{- 3} . \end{matrix}

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})\}}^{\frac{1}{2 - a_{m}}},

(28a)

A_{m} = \frac{g_{m} (2 - a_{m})}{1 + b_{m}} = \frac{g_{m} m r}{1 + m r},

(28b)

B_{m} = \frac{{[γ_{m}^{'} (0)]}^{2 - a_{m}}}{A_{m}} - {[γ_{m} (0)]}^{1 + b_{m}} .

(28c)

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,

(29a)

implying via (10b)

a_{m} = 2 - \frac{m r}{1 + m r} = \frac{2 + m r}{1 + m r}

(29b)

and

γ_{m}^{'} (τ) = {\{A_{m} [γ_{m} (τ) + B_{m}^{(1)}]\}}^{1 + \frac{1}{m r}},

(29c)

with

A_{m}

, respectively,

B_{m}^{(1)}

defined by (28b) (or, equivalently, (12d)), respectively, (12e). Here, we assume, of course, that

r \neq 0

(implying

a_{m} \neq 2

) and moreover that

r \neq - 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) \equiv γ_{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

\tilde{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

τ = {\overset{ˇ}{τ}}_{m}^{(1)}

(see (12b)). In this case,

γ_{m} (t) \equiv γ_{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

{\overset{ˇ}{τ}}_{m}^{(1)}

falls on the circle

\tilde{C}

—note that

{\overset{ˇ}{τ}}_{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 ω {\overset{ˇ}{τ}}_{m}^{(1)}| = 1,

(30a)

then the function

γ_{m} (t) \equiv γ_{m} (τ (t))

diverges at the times

{\overset{ˇ}{t}}_{m}^{(1)} = {(i ω)}^{- 1} ln (1 + i ω {\overset{ˇ}{τ}}_{m}^{(1)}) mod (T_{0}) .

(30b)

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

τ = {\overset{ˇ}{τ}}_{m}^{(1)}

(see (12b)), then the evolution of

γ_{m} (t) \equiv γ_{m} (τ (t))

as function of the time t depends on the location of the branch point

{\overset{ˇ}{τ}}_{m}^{(1)}

in the complex τ-plane, whether it falls outside, inside or just on the circle

\tilde{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

\tilde{C}

is validity of the inequality

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(1)}| > 1,

(31)

and for the corresponding nongeneric initial data the function

γ_{m} (t) \equiv γ_{m} (τ (t))

is again periodic with period

T_{0}

. If instead the branch point falls inside the circle

\tilde{C}

, and the condition on the initial data for this to happen is validity of the opposite inequality

|1 + i ω {\overset{ˇ}{τ}}_{m}^{(1)}| < 1,

(32)

the periodicity of

γ_{m} (t) \equiv γ_{m} (τ (t))

gets modified: the period is then

{\tilde{T}}_{p} = p T_{0}

(see (12b) and (3d)), unless

p / m = {\overset{ˇ}{p}}_{m}

is an integer, in which case the period is

{\overset{ˇ}{T}}_{{\overset{ˇ}{p}}_{m}} = {\overset{ˇ}{p}}_{m} T_{0} .

Next, let us discuss the t-periodicity of

w (t),

related to

γ (t) \equiv γ (τ (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

{\overset{ˇ}{T}}_{{\overset{ˇ}{p}}_{m}} = {\overset{ˇ}{p}}_{m} T_{0}

, where

{\overset{ˇ}{p}}_{m} = p

unless

p / m

is an integer, in which case

{\overset{ˇ}{p}}_{m} = p / m

. We may therefore conclude that, for all generic initial data, the functions

w_{m} (t)

are periodic with period

{\overset{ˇ}{T}}_{{\overset{ˇ}{p}}_{m}}

w_{m} (t + {\overset{ˇ}{T}}_{{\overset{ˇ}{p}}_{m}}) = w_{m} (t) .

(33)

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

\tilde{T}

P_{N} (z; t + \tilde{T}) = P_{N} (z; t)

, the unordered set

\underset{̲}{z} (t)

of its N zeros

z_{n} (t)

is of course periodic with the same period

\tilde{T},

\underset{̲}{z} (t + \tilde{T}) = \underset{̲}{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

\tilde{T}

; although of course that multiple cannot exceed the number

N!

of permutations of the N elements of the unordered set

\underset{̲}{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 \leq K \leq N! p,

(34)

since the very definition of

{\overset{ˇ}{p}}_{m}

(see above) implies that p is the Minimum Common Multiple of the N parameters

{\overset{ˇ}{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},

(35a)

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) - \frac{ln [1 - g_{m} γ_{m}^{'} (0) τ]}{g_{m}},

(35b)

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

{\ddot{w}}_{m} = i ω {\dot{w}}_{m} + g {({\dot{w}}_{m})}^{2}

(36a)

is Hamiltonian, since it is implied by the standard Hamiltonian equations of motion

\dot{w} = \frac{\partial H (\vec{v}, \vec{w}; t)}{\partial v_{m}}, \dot{v} = - \frac{\partial H (\vec{v}, \vec{w};)}{\partial w_{m}},

(36b)

with the (time-dependent) Hamiltonian

H (\vec{v}, \vec{w}; t) = exp (i ω t) v^{α} exp (g α w),

(36c)

as the diligent reader will easily verify. Note the arbitrariness of the (nonvanishing!) parameter α.

However, via the two identities

{\dot{z}}_{n} = - [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] \sum_{m = 1}^{N} [\dot{w} {(z_{n})}^{N - m}],

(37a)

{\ddot{z}}_{n} = \sum_{ℓ = 1, ℓ \neq n}^{N} (\frac{2 {\dot{z}}_{n} {\dot{z}}_{ℓ}}{z_{n} - z_{ℓ}}) - [\prod_{ℓ = 1, ℓ \neq n}^{N} {(z_{n} - z_{ℓ})}^{- 1}] \sum_{m = 1}^{N} [\ddot{w} {(z_{n})}^{N - m}],

(37b)

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}},

(38a)

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}})]}^{\frac{1}{2 - a_{m}}},

(38b)

α_{m} = \frac{b_{m} + 1 + (s_{m} - 1) (2 - a_{m})}{s_{m}}, β_{m} = \frac{(s_{m} - 1) (2 - a_{m})}{s_{m}},

(38c)

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} = - \frac{1}{m r}, b_{m} \equiv b_{m}^{(3)} = - 1 + \frac{1}{m r}, a_{m} \equiv a_{m}^{(3)} = 2 - \frac{1}{1 + m r} .

(39)

The ODE (38b) then reads

η_{m}^{'} = - \frac{1}{m r} {[A_{m} (1 + B_{m}^{(3)} η_{m})]}^{1 + m r},

(40)

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 \neq 0, r \neq - 1 / m

). This ODE, (40), can then be immediately integrated, yielding

η_{m} (τ) = [\frac{1}{B_{m}^{(3)}} + η_{m} (0)] {(1 - \frac{τ}{{\overset{ˇ}{τ}}_{m}^{(3)}})}^{- \frac{1}{m r}} - \frac{1}{B_{m}^{(3)}},

(41)

B_{m}^{(3)} = \frac{{[γ^{'} (0)]}^{\frac{1}{1 + m r}}}{A_{m}} - {[γ (0)]}^{\frac{1}{m r}},

(42)

{\overset{ˇ}{τ}}_{m}^{(3)} = - \frac{{[1 + B_{m}^{(3)} η_{m} (0)]}^{- m r}}{{(A_{m})}^{1 + m r} B_{m}^{(3)}} .

(43)

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

Calogero, F. New solvable variants of the goldfish many-body problem. Studies Appl. Math. 2016. [Google Scholar] [CrossRef]
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]
Bihun, O.; Calogero, F. A new solvable many-body problem of goldfish type. J. Nonlinear Math. Phys. 2016, 23, 28–46. [Google Scholar] [CrossRef]
Bihun, O.; Calogero, F. Novel solvable many-body problems. J. Nonlinear Math. Phys. 2016, 23, 190–212. [Google Scholar] [CrossRef]
Calogero, F. A solvable N-body problem of goldfish type featuring N² arbitrary coupling constants. J. Nonlinear Math. Phys. 2016, 23, 300–305. [Google Scholar] [CrossRef]
Calogero, F. Novel isochronous N-body problems featuring N arbitrary rational coupling constants. J. Math. Phys. 2016, in press. [Google Scholar]
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]
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]
Calogero, F. Classical many-body problems amenable to exact treatments; Springer: Heidelberg, Germany, 2001; p. 749. [Google Scholar]
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]
Gómez-Ullate, D.; Sommacal, M. Periods of the goldfish many-body problem. J. Nonlinear Math. Phys. 2005, 12, 351–362. [Google Scholar] [CrossRef]
Suris, Y.B. Time discretization of F. Calogero’s “Goldfish” System. J. Nonlinear Math. Phys. 2005, 12, 633–647. [Google Scholar] [CrossRef]
Calogero, F.; Iona, S. Novel solvable extensions of the goldfish many-body model. J. Math. Phys. 2005, 46, 103515. [Google Scholar] [CrossRef]
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]
Bruschi, M.; Calogero, F. Novel solvable variants of the goldfish many-body model. J. Math. Phys. 2006, 47, 022703. [Google Scholar] [CrossRef]
Calogero, F.; Langmann, E. Goldfishing by gauge theory. J. Math. Phys. 2006, 47, 082702. [Google Scholar] [CrossRef]
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]
Calogero, F. Isochronous Systems; Oxford University Press: Oxford, UK, 2008; p. 264. [Google Scholar]
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]
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]
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]
Sommacal, M. The Transition from Regular to Irregular Motion, Explained as Travel on Riemann Surfaces. Ph.D. Thesis, SISSA, Trieste, Italy, 2005. [Google Scholar]
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]
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]
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]

© 2016 by the author; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).

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

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

Article Menu

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

Abstract

1. Introduction

2. Results

3. Proofs

4. Outlook

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI