Studying the two-body problem associated to an anisotropic Schwarzschild-type field, Mioc et al. (2003) did not succeed in proving the existence or non-existence of periodic orbits. Here we answer this question in the affirmative. To do this, we start from two basic facts: (1) the potential generates a strong force in Gordon’s sense; (2) the vector field of the problem exhibits the symmetries \(S_{i},i=1,7\), which form, along with the identity, an Abelian group of order 8 with three generators of order 2. Resorting to \(S_{2}\) and \(S_{3}\), in connection with variational methods (particularly the classical lower-semicontinuity method), we prove the existence of infinitely many \(S_{2}\)- or \(S_{3}\)-symmetric periodic solutions. The symmetries \(S_{2}\) and \(S_{3}\) constitute an indicator of the robustness of the classical isotropic Schwarzschild-type system to perturbations (as the anisotropy may be considered).
Authors
Vasile Mioc
Astronomical Institute of the Romanian Academy, Romania
Mira Cristiana Anisiu
Tiberiu Popoviciu Institute of Numerical Analysis of the Romanian Academy, Romania
Michael Barbosu
SUNY Brockport, Department of Mathematics, Brockport, NY 14420, U.S.A.,
[1] H. Amann (1990) Ordinary Differential Equations: An Introduction to Nonlinear Analysis Walter de Gruyer Berlin, New York Google Scholar
[2] A. Ambrosetti V. Coti Zelati (1993) Periodic Solutions of Singular Lagrangian Systems. Progresses in Nonlinear Differential Equations and their Applications, No. 10 Birkhäuser Boston Google Scholar
[3] M.-C. Anisiu (1998) Methods of Nonlinear Analysis Applied to Celestial Mechanics Cluj University Press Cluj-Napoca Google Scholar
[4] M. L. Bertotti (1991)‘Forced oscillations of singular dynamical systems with an application to the restricted three-body problem’J. Diff. Eq.93102–14110.1016/0022-0396(91)90023-3 ArticleGoogle Scholar
[5] Casasayas, J. and Llibre, J.: 1984, ‘Qualitative analysis of the anisotropic Kepler problem’,Mem. Amer. Math. Soc., Vol. 52, No. 312, AMS, Providence, RI.
[6] Chenciner, A.: 2002, ‘Action minimizing periodic orbits in the Newtoniann-body problem’,Celestial Mechanics(Evanston, IL, 1999),Contemporary Mathematics292, Amer. Math. Soc., Providence, RI, pp. 71–90.
[7] A. Chenciner R. Montgomery (2000)‘A remarkable periodic solution of the three-body problem in the case of equal masses’Ann. Math.152881–901 Google Scholar
[8] Coti Zelati, V.: 1994,Introduction to Variational Methods and Singular Lagrangian Systems, School and Workshop on Variational and Local Methods in the Study of Hamiltonian Systems, International Centre for Theoretical Physics, Trieste, Italy, 10–28 October 1994, SMR 779/4.
[9] S. Craig F. N. Diacu E. A. Lacomba E. Perez (1999)‘On the anisotropic Manev problem’J. Math. Phys.401359–137510.1063/1.532807 ArticleGoogle Scholar
[10] R. L. Devaney (1978)‘Collision orbits in the anisotropic Kepler problem’Invent. Math.45221–25110.1007/BF01403170 ArticleGoogle Scholar
[11] Devaney R. L. 1981, ‘Singularities in classical mechanical systems’, inErgodic Theory and Dynamical Systems, Vol. 1, Birkhäuser, Boston, pp. 211–333.
[12] F. N. Diacu (1996)‘Near-collision dynamics for particle systems with quasihomogeneous potentials’J. Diff. Eq.12858–7710.1006/jdeq.1996.0089 ArticleGoogle Scholar
[13] F. Diacu M. Santoprete (2001)‘Nonintegrability and chaos in the anisotropic Manev problem’Physica D15639–5210.1016/S0167-2789(01)00248-2, ArticleGoogle Scholar
[14] Diacu, F. and Santoprete, M.: 2002, ‘On the global dynamics of the anisotropic Manev problem’, http://arXiv.org/abs/nlin/?0208012
[15] W. B. Gordon (1975)‘Conservative dynamical systems involving strong forces’Trans. Amer. Math. Soc.204113–135 Google Scholar
[16] M. C. Gutzwiller (1971)‘Periodic orbits and classical quantization conditions’J. Math. Phys.12343–35810.1063/1.1665596, ArticleGoogle Scholar
[17] M. C. Gutzwiller (1973)‘The anisotropic Kepler problem in two dimensions’J. Math. Phys.14139–15210.1063/1.1666164 ArticleGoogle Scholar
[18] M. C. Gutzwiller (1977)‘Bernoulli sequences and trajectories in the anisotropic Kepler problem’J. Math. Phys.18806–82310.1063/1.523310, ArticleGoogle Scholar
[19] M. Hénon C. Heiles (1964)‘The applicability of the third integral of motion: some numerical experiments’Astron. J.6973–7910.1086/109234 ArticleGoogle Scholar
[20] R. McGehee (1973)‘A stable manifold theorem for degenerate fixed points with applications to celestial mechanics’J. Diff. Eq.1470–8810.1016/0022-0396(73)90077-6, ArticleGoogle Scholar
[21] R. McGehee (1974)‘Triple collision in the collinear three-body problem’Invent. Math.27191–22710.1007/BF01390175, ArticleGoogle Scholar
[22] V. Mioc (2002)‘Symmetries in the Schwarzschild problem’Baltic Astron.11393–407, Google Scholar
[23] V. Mioc E. Pérez-Chavela M. Stavinschi (2003)‘The anisotropic Schwarzschild-type problem, main features’Celest. Mech. Dyn. Astron.8681–[24]10.1023/A:1023648616687, ArticleGoogle Scholar
[25] V. Mioc E. Radu (1992)‘Orbits in an anisotropic radiation field’Astron. Nachr.313353–357, Google Scholar
[26] R. Palais (1979)‘The principle of symmetric criticality’Comm. Math. Phys.6919–3010.1007/BF01941322, ArticleGoogle Scholar
[27] H. Poincaré (1896)‘Sur les solutions périodiques et Ie principe de la moindre action’C. R. Acad. Sci. Paris123915–918, Google Scholar
[28] M. Santoprete (2002)‘Symmetric periodic solutions of the anisotropic Manev problem’J. Math. Phys.433207–321910.1063/1.1469670, ArticleGoogle Scholar
[29] W. C. Saslaw (1978)‘Motion around a source whose luminosity changes’Astrophys. J.226240–25210.1086/156603, ArticleGoogle Scholar
[30] Schwarzschild, K.: 1916, ‘Über das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie’,Sitzber. Preuss. Akad. Wiss. Berlin, 189–196.
[31] C. Stoica V. Mioc (1997)‘The Schwarzschild problem in astrophysics’Astrophys. Space Sci.249161–17310.1023/A:1000347014891, ArticleGoogle Scholar
[32] M. Struwe (1996) Variational Methods. Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems Springer-Verlag Berlin, Google Scholar
[33] L. Tonelli (1915)‘Sur une méthode directe du calcul des variations’Rend. Circ. Mat. Palermo39233–263, Google Scholar
[34] J. P. Vinti (1972)‘Possible effects of anisotropy ofGon celestial orbits’Celest. Mech.6198–20710.1007/BF01227781, ArticleGoogle Scholar
[35] C. Will (1971)‘Relativistic gravity in the solar system. II. Anisotropy in the Newtonian gravitational constant’Astrophys. J.169141–15510.1086/151125, ArticleGoogle Scholar
2005MiocAnisiuBSymmetric
SYMMETRIC PERIODIC ORBITS IN THE ANISOTROPIC SCHWARZSCHILD-TYPE PROBLEM
VASILE MIOC ^(1){ }^{1}, MIRA-CRISTIANA ANISIU ^(2){ }^{2} and MICHAEL BARBOSU ^(3){ }^{3}^(1){ }^{1} Astronomical Institute of the Romanian Academy, Str. Cuţitul de Argint 5, RO-040557 Bucharest, Romania, e-mail: vmioc@aira.astro.ro^(2){ }^{2} T. Popoviciu Institute of Numerical Analysis of the Romanian Academy, P.O. Box 68, 400110 Cluj-Napoca, Romania, e-mail: mira@math.ubbcluj.ro^(3){ }^{3} SUNY Brockport, Department of Mathematics, Brockport, NY 14420, U.S.A., e-mail: mbarbosu@brockport.edu
(Received: 31 March 2004; revised: 16 July 2004; accepted: 29 July 2004)
Abstract
Studying the two-body problem associated to an anisotropic Schwarzschild-type field, Mioc et al. (2003) did not succeed in proving the existence or non-existence of periodic orbits. Here we answer this question in the affirmative. To do this, we start from two basic facts: (1) the potential generates a strong force in Gordon's sense; (2) the vector field of the problem exhibits the symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7}, which form, along with the identity, an Abelian group of order 8 with three generators of order 2 . Resorting to S_(2)S_{2} and S_(3)S_{3}, in connection with variational methods (particularly the classical lower-semicontinuity method), we prove the existence of infinitely many S_(2)S_{2} - or S_(3)S_{3}-symmetric periodic solutions. The symmetries S_(2)S_{2} and S_(3)S_{3} constitute an indicator of the robustness of the classical isotropic Schwarzschild-type system to perturbations (as the anisotropy may be considered).
Astronomy provides a lot of concrete situations that can be tackled via anisotropic mathematical models. The anisotropy of the gravitational constant was discussed by many authors (e.g. Will, 1971, Vinti, 1972). The two-dimensional galactic models also join this class of problems. We further mention: motions around a luminous accretion disk, far orbits around binary stars, orbits around stars with unequal luminosity over surface (pulsars, stars with spots), orbits in proto-stellar systems (and in the proto-solar, too), motion of bodies (from dust to satellites) around planets, in the field of radiation re-emitted by these ones (Saslaw, 1978; Mioc and Radu, 1992). Even the celebrated model of Hénon and Heiles (1964) involves anisotropy.
Gutzwiller (1971, 1973, 1977) defined the anisotropic Kepler problem (namely the anisotropic two-body problem associated to the Newtonian potential) with an essential goal: to identify links between classical and
quantum mechanics. Devaney (1978,1981)(1978,1981) and Casasayas and Llibre (1984) went deeper into this problem.
According to a suggestion formulated by Diacu (1996), the anisotropic Manev problem (associated to a classic potential of the form alpha//r+beta//r^(2)\alpha / r+\beta / r^{2} ) was tackled with a more ambitious purpose: to find connections between classical, quantum, and relativistic mechanics. Important results in this problem were obtained by Craig et al. (1999), Diacu and Santoprete (2001, 2002), Santoprete (2002).
Combining Gutzwiller's anisotropy with Schwarzschild's (1916) potential (of the form alpha//r+beta//r^(3)\alpha / r+\beta / r^{3} ), Mioc et al. (2003; hereafter Paper I) considered the anisotropic Schwarzschild-type problem. They used the powerful tools of the theory of dynamical systems to depict the main features of the global flow.
However, Paper I did not offer an answer to a crucial question for all dynamical systems: does the model admit periodic solutions? In this paper we solve this problem in the affirmative.
In Section 2 we recall the basic equations of the problem in configurationmomentum coordinates. We show that the anisotropic Schwarzschild-type potential generates a strong force in Gordon's (1975) sense. We also show that the corresponding vector field benefits of seven symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7}, which form, along with the identity, an Abelian group of order 8 with three generators of order 2.
In Section 3 we expose some basic notions related to Sobolev spaces, which are the natural frame for finding periodic solutions by variational methods. To get certain families of periodic orbits, we resort to the symmetries S_(2),S_(3)S_{2}, S_{3} in connection with variational methods. The subspaces of symmetric periodic paths Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} will be divided into homotopy classes by their winding number (rotation index). The action integral A_(T)A_{T} will be the functional whose critical points will be the periodic solutions of the Schwarzschild-type problem.
Section 4 contains auxiliary results, which allow us to prove the existence of critical points on some subsets with symmetry and topological constraints, which will be critical points for A_(T)A_{T} on the whole space.
Section 5 presents the main results of our endeavours. Following the methods used by Gordon (1975), Ambrosetti and Coti Zelati (1993), and Coti Zelati (1994) for general strong force fields, and by Diacu and Santoprete (2002) for the anisotropic Manev problem, we resort to the classical lowersemicontinuity method (Tonelli, 1915) to find a minimizer of the action integral in each class, avoiding the collision-type or escape-type minimizers. We prove that, for any pre-assigned period, the action integral has a critical point in each homotopy class with nonnull winding number in Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3}. This implies the existence of infinitely many S_(i)S_{i}-symmetric ( i=2,3i=2,3 ) periodic solutions of the anisotropic Schwarzschild-type problem.
Some remarks are to be formulated here. Even if the idea of using variational methods to find periodic orbits in the planar three-body problem for a potential force of the type 1//r^(n)1 / r^{n} with n >= 2n \geqslant 2 goes as far back as the end of the 19th century (Poincaré, 1896), effective results appeared only relatively recently. The most celebrated result in this context was obtained by Chenciner and Montgomery (2000), who connected symmetries to variational principles to find a new periodic solution of the three-body problem with equal masses. As it was mentioned by Anisiu (1998) for the Manev problem, and it is obviously true for the Schwarzschild one, in these cases (unlike in the Newtonian case), the force field is 'strong' (according to Gordon's (1975) definition), which makes variational methods easier to apply. Interesting results of this type were obtained by Bertotti (1991) for the restricted three-body problem.
Another remark concerns the existence of S_(i)S_{i}-symmetric ( i=2,3i=2,3 ) periodic solutions. Such solutions exist in the isotropic Schwarzschild-type problem (see Stoica and Mioc, 1997; Mioc, 2002). Their persistence (even deformed) in the anisotropic case (regarded as a perturbation of the isotropic case) makes the symmetries S_(2),S_(3)S_{2}, S_{3} constitute an indicator of the robustness of the system to perturbations.
We tried to expose in some detail the necessary mathematical tools, in connection with which we must not forget two issues: (a) almost all such mathematical methods were born from and intended to tackle concrete astronomical situations; (b) our present results add new, important features to the dynamics of the anisotropic Schwarzschild-type problem.
2. Basic Equations and Properties
The 2D2 D anisotropic Schwarzschild-type problem is described by the two-degrees-of-freedom system of ODE
where mu > 0\mu>0 and b > 0b>0 are parameters.
We mention that the Lagrangian L(q,p)=|p|^(2)//2+W(q)L(\mathbf{q}, \mathbf{p})=|\mathbf{p}|^{2} / 2+W(\mathbf{q}) of the anisotropic Schwarzschild-type problem (1) is
and has the property that L(q,p) > 0L(\mathbf{q}, \mathbf{p})>0.
Equations (1) define the motion of a unit-mass particle with respect to another unit-mass particle in an anisotropic plane, namely a plane in which the attraction forces act differently in every direction. The force function W(q)W(\mathbf{q}) characterizes the anisotropy of the plane as a function of the parameter mu\mu. For mu > 1\mu>1 the attraction is the strongest in the q_(1)q_{1}-direction and the weakest in the q_(2)q_{2}-direction; for mu < 1\mu<1 the situation is reversed. For mu=1\mu=1 we retrieve the classical Schwarzschild-type two-body problem, whose global flow was fully depicted by Stoica and Mioc (1997). We shall consider, without loss of generality, that mu > 1\mu>1.
One sees that the Hamiltonian (2) is the sum of the kinetic (K(p(t))=|p(t)|^(2)//2)\left(K(\mathbf{p}(t))=|\mathbf{p}(t)|^{2} / 2\right) and potential ( -W(q(t))-W(\mathbf{q}(t)) ) energies. It provides the first integral of energy
{:(5)H(q(t)","p(t))= tilde(h)","t inR:}\begin{equation*}
H(\mathbf{q}(t), \mathbf{p}(t))=\tilde{h}, t \in \mathbb{R} \tag{5}
\end{equation*}
where tilde(h)\tilde{h} stands for the energy constant.
Notice that, unlike in the classical Schwarzschild model, the force derived from the potential function W is not central; the anisotropy of the plane destroys the rotational invariance. Consequently, the angular momentum C(t)=p(t)xxq(t)C(t)=\mathbf{p}(t) \times \mathbf{q}(t) is not conserved; it does not provide a first integral.
A basic property of system (1) with WW given by (3) is that it satisfies the strong force condition.
DEFINITION 1 (Gordon, 1975). A potential function W:R^(2)\\{(0,0)}rarrRW: \mathbb{R}^{2} \backslash\{(0,0)\} \rightarrow \mathbb{R} generates a strong force if there exist a neighbourhood NN of (0,0)(0,0) and a C^(2)C^{2}-class function U:N\\{(0,0)}rarrRU: N \backslash\{(0,0)\} \rightarrow \mathbb{R} such that
(i) U(q_(1),q_(2))rarr-ooU\left(q_{1}, q_{2}\right) \rightarrow-\infty as (q_(1),q_(2))rarr(0,0)\left(q_{1}, q_{2}\right) \rightarrow(0,0);
(ii) W(q_(1),q_(2)) >= (del U//delq_(1))^(2)+(del U//delq_(2))^(2)=|grad U|^(2)quadW\left(q_{1}, q_{2}\right) \geqslant\left(\partial U / \partial q_{1}\right)^{2}+\left(\partial U / \partial q_{2}\right)^{2}=|\nabla U|^{2} \quad for all quad(q_(2),q_(2))quad\quad\left(q_{2}, q_{2}\right) \quad in N\\{(0,0)}N \backslash\{(0,0)\}.
Remark 2. If there exist the constants c,R > 0c, R>0 such that the potential function W:R^(2)\\{(0,0)}rarrRW: \mathbb{R}^{2} \backslash\{(0,0)\} \rightarrow \mathbb{R} verifies the inequality
then WW generates a strong force, with U(q)=sqrtcln(|q|)U(\mathbf{q})=\sqrt{c} \ln (|\mathbf{q}|) and N={q:|q| < R}N=\{\mathbf{q}:|\mathbf{q}|<R\}. Condition (6), which is easier to be verified, appears as the definition of strong forces in many papers.
THEOREM 3. The potential function WW given by (3) generates a strong force.
Proof. From the expression (3) of the potential function WW and from mu > 1\mu>1 it is obvious that W(q_(1),q_(2)) >= b(muq_(1)^(2)+q_(2)^(2))^(-3//2) >= bmu^(-3//2)(q_(1)^(2)+q_(2)^(2))^(-3//2) >= bmu^(-3//2)|q|^(2)W\left(q_{1}, q_{2}\right) \geqslant b\left(\mu q_{1}^{2}+q_{2}^{2}\right)^{-3 / 2} \geqslant b \mu^{-3 / 2}\left(q_{1}^{2}+q_{2}^{2}\right)^{-3 / 2} \geqslant b \mu^{-3 / 2} |\mathbf{q}|^{2} for 0 < |q| < 10<|\mathbf{q}|<1, hence WW satisfies (6) with c=bmu^(-3//2)c=b \mu^{-3 / 2} and R=1R=1.
By (2) and (3), the motion Equations (1) explicitly read.
An important property of these equations can be stated as
THEOREM 4. The vector field (7) benefits of seven symmetries S_(i)=S_(i)(q_(1),q_(2),p_(1),p_(2),t),i= bar(1,7)S_{i}= S_{i}\left(q_{1}, q_{2}, p_{1}, p_{2}, t\right), i=\overline{1,7}, as follows:
The set G={I}uu{S_(i)∣i= bar(1,7)}G=\{I\} \cup\left\{S_{i} \mid i=\overline{1,7}\right\}, endowed with the usual composition law "。", forms a symmetric Abelian group with an idempotent structure, isomorphic to Z_(2)xxZ_(2)xxZ_(2)\mathbb{Z}_{2} \times \mathbb{Z}_{2} \times \mathbb{Z}_{2}, where II denotes the identity. This group owns seven proper subgroups isomorphic to Klein's group.
Proof. It is easy to check that Equations (7) are invariant to the transformations (8). As regards the Abelian group structure of GG, it suffices to construct and examine the following composition table:
Since every element is its own inverse, the idempotent structure is obvious.
Observe that, among the symmetries (8), only three are independent. Consider, for instance, that these ones are S_(1),S_(2),S_(3)S_{1}, S_{2}, S_{3}. The relations S_(1)@S_(2)=S_(4),S_(1)@S_(3)=S_(5),S_(2)@S_(3)=S_(6),S_(1)@S_(2)@S_(3)=S_(7)S_{1} \circ S_{2}=S_{4}, S_{1} \circ S_{3}=S_{5}, S_{2} \circ S_{3}=S_{6}, S_{1} \circ S_{2} \circ S_{3}=S_{7} are immediate. Every other three independent symmetries generate the remaining four ones. This means that GG is an Abelian group of order eight with three generators of order two. By the Fundamental Theorem of Abelian Groups, GG is isomorphic to Z_(2)xxZ_(2)xxZ_(2)\mathbb{Z}_{2} \times \mathbb{Z}_{2} \times \mathbb{Z}_{2}.
As to the proper subgroups of GG, let us denote them by G_(ijk)={I,S_(i):}G_{i j k}=\left\{I, S_{i}\right., {:S_(j),S_(k)∣i!=j!=k!=i,S_(i)@S_(j)=S_(k)}\left.S_{j}, S_{k} \mid i \neq j \neq k \neq i, S_{i} \circ S_{j}=S_{k}\right\}. We can immediately check that the only sets {i,j,k}\{i, j, k\} that fulfil this condition lead to the subgroups tilde(G)_(124), tilde(G)_(135), tilde(G)_(167), tilde(G)_(236), tilde(G)_(257), tilde(G)_(347), tilde(G)_(456)\tilde{G}_{124}, \tilde{G}_{135}, \tilde{G}_{167}, \tilde{G}_{236}, \tilde{G}_{257}, \tilde{G}_{347}, \tilde{G}_{456}. All these subgroups are of order four with two generators of order two, hence isomorphic to Klein's subgroup. This completes the proof.
In Paper I we have shown that Theorem 4 also holds for the motion equations expressed in collision-blow-up or infinity-blow-up McGehee-type coordinates (McGehee, 1973, 1974). The respective groups of 4 symmetries, G_(0)G_{0} and G_(oo)G_{\infty}, are isomorphic to GG. This is not a trivial result, because the phase spaces corresponding to G_(0)G_{0} and to G_(oo)G_{\infty} contain the supplementary boundary manifolds of collision and infinity, respectively.
As mentioned in the introductory section, in Paper I the important question about the existence or not of periodic orbits remained open. This question will be answered in the affirmative in what follows, where the symmetries (8) and the strong force property of the potential will play a premier role.
3. The Functional Background
In order to get certain families of periodic orbits, we shall use some symmetry and topological constraints in connection with a variational principle, finding periodic orbits not as minimizers, but as extremal values of the action integral. To this end we shall resort to results concerning periodic solutions of fixed period for symmetric, singular, Lagrangian systems (see Gordon 1975; Ambrosetti and Coti Zelati, 1993; Coti Zelati, 1994), as the system (7) is.
We first need some notations. Given a number T > 0T>0, let us denote by C^(oo)([0,T],R^(2))C^{\infty}\left([0, T], \mathbb{R}^{2}\right) the space of TT-periodic C^(oo)C^{\infty} cycles (loops) f:[0,T]rarrR^(2)f:[0, T] \rightarrow \mathbb{R}^{2}. We define the inner products
and let ||f||_(L^(2))=(int_(0)^(T)|f|^(2)(d)t)^(1//2)\|f\|_{L^{2}}=\left(\int_{0}^{T}|f|^{2} \mathrm{~d} t\right)^{1 / 2} and ||f||_(H^(1))=(int_(0)^(T)(|(f^(˙))|^(2)+|f|^(2))dt)^(1//2)\|f\|_{H^{1}}=\left(\int_{0}^{T}\left(|\dot{f}|^{2}+|f|^{2}\right) \mathrm{d} t\right)^{1 / 2} be the corresponding norms. Then the completion of C^(oo)([0,T],R^(2))C^{\infty}\left([0, T], \mathbb{R}^{2}\right) with respect to ||*||_(L^(2))\|\cdot\|_{L^{2}} is denoted by L^(2)L^{2} and is the space of square integrable functions. The completion with respect to ||*||_(H^(1))\|\cdot\|_{H^{1}} is denoted by H^(1)H^{1} and is the Sobolev space of all absolutely continuous TT-periodic functions that have L^(2)L^{2} derivatives defined almost everywhere. The space H^(1)H^{1} is compactly embedded in the space of continuous functions on [0,T],C^(0)([0,T],R^(2))[0, T], C^{0}\left([0, T], \mathbb{R}^{2}\right) with ||f||=max{|f(t)|:t in[0,T]}\|f\|=\max \{|f(t)|: t \in [0, T]\} (Gordon, 1975).
The Schwarzschild potential (3) has a singularity at the origin of R^(2)\mathbb{R}^{2}, hence we shall denote by
{:(11)Lambda={f inH^(1)∣f(t)!=(0,0)quad" for all "t in[0,T]}:}\begin{equation*}
\Lambda=\left\{f \in H^{1} \mid f(t) \neq(0,0) \quad \text { for all } t \in[0, T]\right\} \tag{11}
\end{equation*}
the open subset of the cycles in H^(1)H^{1} which do not pass through the origin (noncollisional cycles). For the noncollisional cycles in Lambda\Lambda it makes sense to define the angle function vartheta_(f)inC^(0)([0,T],R)\vartheta_{f} \in C^{0}([0, T], \mathbb{R}) by
cos vartheta_(f)(t)=f_(1)(t)//sqrt(f_(1)^(2)(t)+f_(2)^(2)(t)),quad sin vartheta_(f)(t)=f_(2)(t)//sqrt(f_(1)^(2)(t)+f_(2)^(2)(t))\cos \vartheta_{f}(t)=f_{1}(t) / \sqrt{f_{1}^{2}(t)+f_{2}^{2}(t)}, \quad \sin \vartheta_{f}(t)=f_{2}(t) / \sqrt{f_{1}^{2}(t)+f_{2}^{2}(t)}
which measures the angle between the positive q_(1)q_{1}-axis and the vector f(t)f(t) in the mathematically positive direction. The rotation index (or winding number) w(f)w(f) will represent the growth of the angle function during a period, measured in units of full rotations of ff, that is
The rotation index is always an integer which shows how many times the continuous cycle f:[0,T]rarrR^(2)\\{(0,0)}f:[0, T] \rightarrow \mathbb{R}^{2} \backslash\{(0,0)\} 'winds around' the origin, having positive values for counterclockwise rotations and negative for clockwise ones. Identifying S^(1)S^{1} with R//[0,T]\mathbb{R} /[0, T] one has w(f)=deg(f)w(f)=\operatorname{deg}(f), where deg(f)\operatorname{deg}(f) is the degree of the circle map F:S^(1)rarrS^(1),F(t)=f(t)//|f(t)|F: S^{1} \rightarrow S^{1}, F(t)=f(t) /|f(t)| (Amann, 1990). It follows
where Lambda_(k)={f in Lambda∣w(f)=k}\Lambda_{k}=\{f \in \Lambda \mid w(f)=k\}. We mention that Lambda_(0)\Lambda_{0} contains the loops which are homothetic in R^(2)\\{(0,0)}\mathbb{R}^{2} \backslash\{(0,0)\} to a point.
We shall use some of the 'natural' symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7} of the system (7).
Let us denote by Sigma_(i),i= bar(1,7)\Sigma_{i}, i=\overline{1,7} the subsets of H^(1)H^{1} formed by S_(i)S_{i}-symmetric cycles, namely those which satisfy S_(i)(f(t))=f(t)S_{i}(f(t))=f(t). It is clear that each Sigma_(i)\Sigma_{i} is a subspace of H^(1)H^{1}. In the sequel we shall provide orthogonal decompositions of H^(1)H^{1} in terms of its subspaces Sigma_(i)\Sigma_{i} with i in{2,3}i \in\{2,3\} and i in{1,7}i \in\{1,7\}, respectively.
LEMMA 5. The subspaces Sigma_(i)\Sigma_{i} with i in{1,2,3,7}i \in\{1,2,3,7\} are closed, weakly closed, and complete with respect to ||*||_(H^(1))\|\cdot\|_{H^{1}}, hence they are Sobolev spaces. Moreover,
Proof. Consider a function f=(f_(1),f_(2))f=\left(f_{1}, f_{2}\right); it is well known that f_(1)f_{1} and f_(2)f_{2} can be written as sums of an even absolutely continuous function and an odd one, namely f_(j)=f_(j,e)+f_(j,o)f_{j}=f_{j, e}+f_{j, o}, where f_(j,e)(t)=(f_(j)(t)+f_(j)(-t))//2f_{j, e}(t)=\left(f_{j}(t)+f_{j}(-t)\right) / 2 and f_(j,o)(t)=(f_(j)(t)-:}{:f_(j)(-t))//2,j=1,2f_{j, o}(t)=\left(f_{j}(t)-\right. \left.f_{j}(-t)\right) / 2, j=1,2. By virtue of this fact, we can write ff as the sum of an S_(2)S_{2}-symmetric function, f_(S_(2))=(f_(1,e),f_(2,o))f_{S_{2}}=\left(f_{1, e}, f_{2, o}\right) and an S_(3)S_{3}-symmetric one, f_(S_(3))=(f_(1,o),f_(2,e))f_{S_{3}}= \left(f_{1, o}, f_{2, e}\right).
Now, let us consider an element f inSigma_(2)f \in \Sigma_{2}. Then (:f,g:)_(H^(1))=0\langle f, g\rangle_{H^{1}}=0 for every g inSigma_(3)g \in \Sigma_{3}. This is due to the fact that, by (9) and (10),
where the second integrand is an odd function, whereas the first one is an odd function almost everywhere. Thus the above inner product is zero for every g inSigma_(3)g \in \Sigma_{3}.
Denote by Sigma_(2)^(_|_)={g inSigma_(2)∣(:f,g:)_(H_(1))=0,AA_(g)inSigma_(2)}\Sigma_{2}^{\perp}=\left\{g \in \Sigma_{2} \mid\langle f, g\rangle_{H_{1}}=0, \forall_{g} \in \Sigma_{2}\right\} the space orthogonal to Sigma_(2)\Sigma_{2}. It is easy to see that Sigma_(2)^(_|_)\Sigma_{2}^{\perp} is closed and that Sigma_(3)subSigma_(2)^(_|_)\Sigma_{3} \subset \Sigma_{2}^{\perp}. To prove that Sigma_(3)=Sigma_(2)^(_|_)\Sigma_{3}=\Sigma_{2}^{\perp}, suppose that there exists varphi inSigma_(2)^(_|_)\varphi \in \Sigma_{2}^{\perp} such that varphi!=0\varphi \neq 0 and varphi!inSigma_(3)\varphi \notin \Sigma_{3}. Then write varphi=varphis_(2)+varphis_(3)\varphi=\varphi s_{2}+\varphi s_{3} and compute (:varphis_(2),varphi:)_(H_(1))=(:varphis_(2),varphis_(2)+varphis_(3):)_(H^(1))=(:varphis_(2),varphis_(2):)_(H^(1))=||varphis_(2)||_(H^(1))=delta > 0\left\langle\varphi s_{2}, \varphi\right\rangle_{H_{1}}=\left\langle\varphi s_{2}, \varphi s_{2}+\varphi s_{3}\right\rangle_{H^{1}}=\left\langle\varphi s_{2}, \varphi s_{2}\right\rangle_{H^{1}}= \left\|\varphi s_{2}\right\|_{H^{1}}=\delta>0. But this contradicts the hypothesis that varphi inSigma_(2)^(_|_)\varphi \in \Sigma_{2}^{\perp}, therefore Sigma_(3)=Sigma_(2)^(_|_)\Sigma_{3}=\Sigma_{2}^{\perp}. So Sigma_(3)\Sigma_{3} and, consequently, Sigma_(2)\Sigma_{2} are closed and such that H^(1)=Sigma_(2)o+Sigma_(3)H^{1}=\Sigma_{2} \oplus \Sigma_{3}. In addition, since H^(1)H^{1} is complete, Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} are complete. Lastly, Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} are weakly closed because they are norm-closed subspaces.
The statements for Sigma_(1)\Sigma_{1} and Sigma_(7)\Sigma_{7} can be proved similarly, by writing f=f_(S_(1))+f_(S_(7))f=f_{S_{1}}+f_{S_{7}}, with f_(S_(1))=(f_(1,e),f_(2,e))f_{S_{1}}=\left(f_{1, e}, f_{2, e}\right) and f_(S_(7))=(f_(1,o),f_(2,o))f_{S_{7}}=\left(f_{1, o}, f_{2, o}\right). This completes the proof.
Consider the sets tilde(Sigma)_(i)=Sigma_(i)nn Lambda\tilde{\Sigma}_{i}=\Sigma_{i} \cap \Lambda of symmetric noncollisional cycles, which are open submanifolds of the spaces Sigma_(i)\Sigma_{i}. We shall say that a cycle in tilde(Sigma)_(i)\tilde{\Sigma}_{i} is of class Lambda_(k),k inZ\Lambda_{k}, k \in \mathbb{Z}, if its winding number about the origin of the coordinate system is kk, namely if it performs kk loops around the origin. We remind that kk is positive for counterclockwise rotations and negative else. The family (Lambda_(k))_(k inZ)\left(\Lambda_{k}\right)_{k \in \mathbb{Z}} provides a partition of tilde(Sigma)_(i)\tilde{\Sigma}_{i} into homotopy classes (components).
We shall describe the geometric properties of the cycles in Sigma_(i),i= bar(1,7)\Sigma_{i}, i=\overline{1,7}. A cycle ff is in Sigma_(1)\Sigma_{1} if and only if f_(j)(-t)=f_(j)(t),j=1,2f_{j}(-t)=f_{j}(t), j=1,2, hence the cycles in tilde(Sigma)_(1)\tilde{\Sigma}_{1} will have null winding number ( Sigma subLambda_(0)\Sigma \subset \Lambda_{0} ). The cycles in Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} have mirror symmetry with respect to the q_(1)q_{1}, respectively to the q_(2)q_{2}-axis. Those in Sigma_(4)\Sigma_{4} and Sigma_(5)\Sigma_{5} are lying on the q_(1)q_{1}, respectively q_(2)q_{2}-axis, hence the cycles in tilde(Sigma)_(4)\tilde{\Sigma}_{4} and tilde(Sigma)_(5)\tilde{\Sigma}_{5} will have again null winding number. The cycles in Sigma_(6)\Sigma_{6} reduce themselves to the single point (0,0)(0,0). The cycles in Sigma_(7)\Sigma_{7} are symmetric with respect to the origin (0,0)(0,0) and all of them pass through (0,0)(0,0), hence they are all collisional.
We are interested in noncollisional families of cycles with nonnull winding numbers; Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} are the only subspaces of H^(1)H^{1} among the seven subspaces corresponding to the natural symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7}, which satisfy those requirements.
Remark 6. To be more clear for a nonmathematician reader, we summarize what we have done so far. We fixed a period TT and we defined the Sobolev space H^(1)H^{1} of all absolutely continuous TT-periodic functions. We defined the subsets Sigma_(i)\Sigma_{i} of H^(1)H^{1} formed by S_(i)S_{i}-symmetric cycles (characterized by the symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7}, of the motion Equation (7)). Moreover, we have shown that the couples (Sigma_(2),Sigma_(3)),(Sigma_(1),Sigma_(7))\left(\Sigma_{2}, \Sigma_{3}\right),\left(\Sigma_{1}, \Sigma_{7}\right) cover - via direct sum - the whole space H^(1)H^{1}. Being interested only in cycles that are noncollisional or nonescape type, or do not represent quasiperiodic orbits, we showed that only Sigma_(2)\Sigma_{2} and Sigma_(3)\Sigma_{3} fulfil these conditions. To continue our mathematical endeavours, we divided every subset tilde(Sigma)_(i)\tilde{\Sigma}_{i} in homotopy classes via the winding number kk; this will be useful further down.
We shall define now the action integral, whose extremal values will provide symmetric periodic orbits. The action integral A_(T):Lambda rarrRA_{T}: \Lambda \rightarrow \mathbb{R} between the instants 0 and TT, along a cycle ff whose Euclidean coordinate representation is q=(q_(1),q_(2))\mathbf{q}=\left(q_{1}, q_{2}\right), has the expression
{:(14)A_(T)(f)=int_(0)^(T)L(q(t)","p(t))dt:}\begin{equation*}
A_{T}(f)=\int_{0}^{T} L(\mathbf{q}(t), \mathbf{p}(t)) \mathrm{d} t \tag{14}
\end{equation*}
with the positive Lagrangian function given by (4). The appearance of the quadratic terms p_(1)^(2)=q^(˙)_(1)^(2)p_{1}^{2}=\dot{q}_{1}^{2} and p_(2)^(2)=q^(˙)_(2)^(2)p_{2}^{2}=\dot{q}_{2}^{2} in the expression of the Lagrangian function makes the Sobolev space H^(1)H^{1} adequate for this problem. The Schwarzschild potential belongs to a class of potentials (with W(q) >= 0W(q) \geqslant 0 for q!=(0,0)\mathbf{q} \neq(0,0) and W(q)rarr0W(\mathbf{q}) \rightarrow 0 for |q|rarr oo|\mathbf{q}| \rightarrow \infty ) for which it can be shown (Coti Zelati, 1994) that i n f_(f in Lambda)A_(T)(f)=0\inf _{f \in \Lambda} A_{T}(f)=0, and all minimizing sequences are unbounded.
In order to obtain periodic solutions of (7) we are forced to minimize the functional A_(T)A_{T} on subsets Lambda_(k)\Lambda_{k} of Lambda\Lambda, chosen by using symmetry and topological constraints. After selecting an adequate subset, we shall use a direct method of the calculus of variation, i.e. the lower-semicontinuity method (e.g., Struwe 1996) and get a minimizer in that subset, which will be proved to be an extremal value of the functional A_(T)A_{T}. Finally we show that the extremal values, which belong to the Sobolev space H^(1)H^{1}, are regular enough to constitute classical periodic solutions of (7).
Remark 7. To sketch in physical terms the minimization of the action, recall that the Lagrangian of our problem represents the sum of two positive terms: the kinetic energy KK and the force function WW (the negative of the potential energy). Also recall that KK and WW are not independent each other, they being related by the energy integral (in which the constant of energy must be negative, provided the positiveness of WW ). Since K,W > 0K, W>0, any minimization of their sum involves the minimization of both KK and WW; both push the trajectory away from the field-generating centre. But the limit imposed by the fixed energy-level and by the fixed value of TT stops the orbit expansion to a finite value, which can lead to a periodic orbit.
4. Auxiliary Results
The first statement establishes the connection between the solutions of the Schwarzschild-type problem (7) and the extremals (critical points) of the functional A_(T)A_{T} given by (14).
PROPOSITION 8. The set of noncollisional cycles Lambda\Lambda given by (11) is an open subset of H^(1)H^{1} and the functional A_(T)A_{T} is in the class C^(1)(Lambda,R)C^{1}(\Lambda, \mathbb{R}) with
if a cycle ff with coordinate expression q=(q_(1),q_(2))in Lambda\mathbf{q}=\left(q_{1}, q_{2}\right) \in \Lambda is a critical point of A_(T)A_{T} on Lambda\Lambda, then ff is a classical periodic solution of (7).
Proof. The first affirmation is standard (see Struwe, 1990). Let us consider f in Lambdaf \in \Lambda which is a critical point of A_(T)A_{T} on Lambda\Lambda. Then ff is continuous, as well as grad W(f(t))\nabla W(f(t)). We take the scalar product of
grad W(f(t))=(d)/((d)t)int_(0)^(t)grad W(f(s))ds\nabla W(f(t))=\frac{\mathrm{d}}{\mathrm{~d} t} \int_{0}^{t} \nabla W(f(s)) \mathrm{d} s
with h inH^(1),h(0)!=0h \in H^{1}, h(0) \neq 0, then integrate on [0,T][0, T] and obtain
Since f in Lambdaf \in \Lambda we obtain via Sobolev embedding that f^(˙)inC^(0)\dot{f} \in C^{0}, and from (17) that f^(˙)inC^(1)\dot{f} \in C^{1}. We integrate (16) by parts and get
(:f^(˙)(T),h(T):)-(:f^(˙)(0),(0):)=int_(0)^(T)(:f^(¨)(t)-grad W(f(t)),h(t):)dt\langle\dot{f}(T), h(T)\rangle-\langle\dot{f}(0),(0)\rangle=\int_{0}^{T}\langle\ddot{f}(t)-\nabla W(f(t)), h(t)\rangle \mathrm{d} t
The right-hand side is null, and since h(0)!=0h(0) \neq 0 it follows. f^(˙)(T)=f^(˙)(0)\dot{f}(T)=\dot{f}(0).
To an element f in Lambdaf \in \Lambda we associate a curve f_(U)f_{U}, with UU given in Remark 2:
{:(18)f_(U)(t)=(f(t)","U(f(t)))inR^(3)","0 <= t <= T:}\begin{equation*}
f_{U}(t)=(f(t), U(f(t))) \in \mathbb{R}^{3}, 0 \leqslant t \leqslant T \tag{18}
\end{equation*}
The next two lemmas contain results from Gordon's (1975) paper and rely on the fact that the force is strong.
LEMMA 9. For each f in Lambdaf \in \Lambda, the following inequality holds:
Lemma 9 will allow us to prove the corresponding of Gordon's geometrical lemma for the anisotropic Schwarzschild potential.
LEMMA 10. The functional A_(T)A_{T} has the property that, for any a > 0a>0, there exists delta=delta(a) > 0\delta=\delta(a)>0 such that if q inH^(1)q \in H^{1} and A_(T)(q) <= aA_{T}(q) \leqslant a, then |q(t)| >= delta|q(t)| \geqslant \delta for all t inRt \in \mathbb{R}, i.e. {q inH^(1):A_(T)(q) <= a}\left\{q \in H^{1}: A_{T}(q) \leqslant a\right\} is bounded away from zero.
Proof. Let us suppose that there is no such delta\delta; it means that for each n inNn \in \mathbb{N} there exist t_(n)^(**)in[0,T)t_{n}^{*} \in[0, T) and f_(n)inH^(1)f_{n} \in H^{1} for which A_(T)(f_(n)) <= aA_{T}\left(f_{n}\right) \leqslant a and i n f{|f_(n)(t)|:t in[0,T]}=|f_(n)(t_(n)^(**))| < 1//n\inf \left\{\left|f_{n}(t)\right|: t \in[0, T]\right\}=\left|f_{n}\left(t_{n}^{*}\right)\right|<1 / n. We have f_(n)inC^(0)f_{n} \in C^{0}. If ||f_(n)||rarr0\left\|f_{n}\right\| \rightarrow 0 it follows that int_(0)^(T)W(f_(n)(t))dt rarr oo\int_{0}^{T} W\left(f_{n}(t)\right) d t \rightarrow \infty as n rarr oon \rightarrow \infty, hence A_(T)(f_(n))A_{T}\left(f_{n}\right) is unbounded, contradiction. If ||f_(n)||↛0\left\|f_{n}\right\| \nrightarrow 0, there exist epsi > 0\varepsilon>0 and a subsequence of f_(n)f_{n}, denoted also f_(n)f_{n}, and t_(n)in[0,T]t_{n} \in[0, T] such that |f_(n)(t_(n))|=epsi,AA_(n)inN\left|f_{n}\left(t_{n}\right)\right|=\varepsilon, \forall_{n} \in \mathbb{N}. The arc length of f_(nU)f_{n U} between t_(n)t_{n} and t_(n)^(**)t_{n}^{*} is greater than the length of the corresponding straight line, hence
arc" length "(f_(nU)) >= |U(t_(n),f_(n)(t))-U(t_(n)^(**),f_(n)(t_(n)^(**)))|rarr oo" as "n rarr oo.\operatorname{arc} \text { length }\left(f_{n U}\right) \geqslant\left|U\left(t_{n}, f_{n}(t)\right)-U\left(t_{n}^{*}, f_{n}\left(t_{n}^{*}\right)\right)\right| \rightarrow \infty \text { as } n \rightarrow \infty .
In view of Lemma 9, this contradicts the boundedness of A_(T)(f_(n))A_{T}\left(f_{n}\right).
Remark 11. As stated by Gordon(1975), the loops in Lambda\\Lambda_(0)\Lambda \backslash \Lambda_{0} cannot be continuously moved off to infinity without either passing through (0,0)(0,0) or having its arc length become infinite (for every c_(1)c_{1}, there exists a compact subset K_(c_(1))K_{c_{1}} of R^(2)\mathbb{R}^{2} which contains every smooth cycle which is homotopic to ff in R^(2)\\{(0,0)}\mathbb{R}^{2} \backslash\{(0,0)\} and has arc length {: <= c_(1))\left.\leqslant c_{1}\right).
The next result is a special case of Palais' (1979) principle of symmetric criticality, as it was presented by Chenciner (2002).
Let us consider an orthogonal (by isometries) representation rho\rho of a finite group GG in the real Hilbert space H^(1)H^{1} such that, for any gamma\gamma in GG,
We denote by H_(rho)^(1)H_{\rho}^{1} the linear subspace of H^(1)H^{1} formed by the elements which are invariant under the representation rho\rho, and by A_(T,rho)A_{T, \rho} the restriction of the action A_(T)A_{T} to H_(rho)^(1)H_{\rho}^{1}. The principle of symmetric criticality asserts that a critical point
for the restriction of the action on the subspace H_(rho)^(1)H_{\rho}^{1} is a critical point of A_(T)A_{T} on the whole loop space H^(1)H^{1}.
PROPOSITION 12. Any critical point of A_(T,rho)A_{T, \rho} is a critical point of A_(T)A_{T}.
Proof. Using the rho\rho-invariance of A_(T)A_{T} given in (19), we have that, for any gamma in G\gamma \in G
Each f inH_(rho)^(1)f \in H_{\rho}^{1} satisfies dA_(T)(f)[rho(gamma)^(-1)*h]=dA_(T)(f)[h]\mathrm{d} A_{T}(f)\left[\rho(\gamma)^{-1} \cdot h\right]=\mathrm{d} A_{T}(f)[h]. The H^(1)H^{1}-gradient gradA_(T)(f)quad\nabla A_{T}(f) \quad will satisfy quad(:rho(gamma)*gradA_(T)(f),h:)=(:gradA_(T)(f),rho(gamma)^(-1)*h:)=(:gradA_(T)(f),h:)\quad\left\langle\rho(\gamma) \cdot \nabla A_{T}(f), h\right\rangle=\left\langle\nabla A_{T}(f), \rho(\gamma)^{-1} \cdot h\right\rangle= \left\langle\nabla A_{T}(f), h\right\rangle. This means that gradA_(T)(f)\nabla A_{T}(f) belongs to H_(rho)^(1)H_{\rho}^{1}, hence each critical point of A_(T,rho)A_{T, \rho} is a critical point of A_(T)A_{T}.
Remark 13. The symmetry S_(2)S_{2} is given by the representation s_(2)s_{2} of the group Z_(2)\mathbb{Z}_{2} acting as
In these cases, a critical point of the action restricted to the invariant subspace with respect to the symmetry will be a critical point of A_(T)A_{T} on H^(1)H^{1}.
The fact that the elements for which the action is bounded are bounded away from zero prevents the critical points from being collisional solutions of the system (7). Another property of the action, namely its coercivity, avoids the critical points at infinity. A real functional AA defined on a Hilbert space with the norm ||*||\|\cdot\| is coercive if A(x_(n))rarr+ooA\left(x_{n}\right) \rightarrow+\infty for all sequences x_(n)x_{n} such that ||x||rarr+oo\|x\| \rightarrow+\infty.
To end, we still need to recall some definitions and issues. Let XX be a topological space, and consider Psi:X rarrR\Psi: X \rightarrow \mathbb{R}. Then Psi\Psi is lower-semicontinuous if and only if Psi^(-1)(-oo,a]\Psi^{-1}(-\infty, a] is closed for every a inRa \in \mathbb{R}, in which case Psi\Psi is bounded from below and reaches its infimum on every compact subset of XX. If XX is a Hausdorff space (thus compact subsets are necessarily closed), the following result (known as Weierstrass' theorem) holds:
PROPOSITION 14. Consider a real-valued function Psi:X rarrR\Psi: X \rightarrow \mathbb{R}, where XX is a Hausdorff space, such that Psi^(-1)(-oo,a]\Psi^{-1}(-\infty, a] is compact for every a inRa \in \mathbb{R}. Then Psi\Psi is lower-semicontinuous bounded from below and reaches its infimum on XX.
Remark 15. In physical terms, we have shown that the symmetric solutions that lead to critical points of the action do not encounter either collision or escape, thus being periodic orbits of the problem.
Now we have at our disposal all the necessary mathematical instruments and results to prove the existence of periodic orbits in the anisotropic Schwarzschild-type problem.
5. Main Results
We state now the central result of our paper.
THEOREM 16. For any T > 0T>0 and any k in Z\\{0}k \in Z \backslash\{0\}, there exists at least one S_(i)S_{i}-symmetric ( i=2,3i=2,3 ) periodic solution of the anisotropic Schwarzschild-type problem, with period TT and winding number kk.
Proof. Let i in{2,3}i \in\{2,3\} be fixed and XX a component of tilde(Sigma)_(i)\tilde{\Sigma}_{i} that consists of non-simple cycles. Endow XX with the weak topology it inherits from Sigma_(i)\Sigma_{i}. We intend to apply Proposition 14 with Psi=A_(T)\Psi=A_{T}, thus we shall show that Y_(a)=X nnA_(T)^(-1)(-oo,a]Y_{a}=X \cap A_{T}^{-1}(-\infty, a] is bounded and weakly closed in Sigma_(i)\Sigma_{i}, hence weakly compact.
Let f inY_(a)f \in Y_{a}, hence f in Xf \in X and 0 <= A_(T)(f) <= a0 \leqslant A_{T}(f) \leqslant a. From Lemma 9 it follows that the elements of Y_(a)Y_{a} are bounded in arc length by the same constant, and from Lemma 10 that they are bounded away from (0,0)(0,0). The elements of XX being tied to (0,0)(0,0), from Remark 11 it follows that XX is bounded in C^(0)C^{0} norm, i.e. there exists c > 0c>0 such that ||f|| <= c\|f\| \leqslant c for all f in Xf \in X. Let f inY_(a)f \in Y_{a}; we have
||f||_(H^(1))^(2)=||f||_(L^(2))^(2)+||f^(˙)||_(L^(2))^(2) <= ||f||^(2)+2A_(T)(f) <= c^(2)+2a,\|f\|_{H^{1}}^{2}=\|f\|_{L^{2}}^{2}+\|\dot{f}\|_{L^{2}}^{2} \leqslant\|f\|^{2}+2 A_{T}(f) \leqslant c^{2}+2 a,
hence Y_(a)Y_{a} is bounded with respect to ||*||_(H^(1))\|\cdot\|_{H^{1}}.
To show that Y_(a)Y_{a} is weakly closed, we consider a sequence f_(n)in X nnA_(T)^(-1)(-oo,a]f_{n} \in X \cap A_{T}^{-1}(-\infty, a] which converges weakly to a cycle f inH^(1)f \in H^{1}. The subspaces Sigma_(i)\Sigma_{i} being weakly closed (as it was proved in Lemma 5), we have that f inSigma_(i)f \in \Sigma_{i}. As stated above, the cycles f_(n)f_{n} are bounded in arc length and bounded away from (0,0)(0,0). The weak convergence in H^(1)H^{1} implies that lim_(n rarr oo)||f_(n)-f||=0\lim _{n \rightarrow \infty}\left\|f_{n}-f\right\|=0 with ||*||\|\cdot\| the C^(0)C^{0} norm. Because f_(n),n inNf_{n}, n \in \mathbb{N}, are bounded away from (0,0)(0,0), there exists delta > 0\delta>0 such that |f_(n)(t)| >= delta\left|f_{n}(t)\right| \geqslant \delta, for each t in[0,T]t \in[0, T] and n inNn \in \mathbb{N}. Then
and making n rarr oon \rightarrow \infty it follows that |f(t)| >= delta,t in[0,T]|f(t)| \geqslant \delta, t \in[0, T]. Therefore f in tilde(Sigma)_(i)f \in \tilde{\Sigma}_{i} and it is in the same component of tilde(Sigma)_(i)\tilde{\Sigma}_{i} as f_(n),n inNf_{n}, n \in \mathbb{N}, hence f in Xf \in X.
Because of the fact that |f_(n)(t)| >= delta,|f(t)| >= delta\left|f_{n}(t)\right| \geqslant \delta,|f(t)| \geqslant \delta for t in[0,T]t \in[0, T], it follows that we may apply Fatou's lemma to obtain
int_(0)^(T)W(f(t))dt=int_(0)^(T)l i m i n f W(f_(n)(t))dt <= l i m i n fint_(0)^(T)W(f_(n)(t))dt\int_{0}^{T} W(f(t)) \mathrm{d} t=\int_{0}^{T} \liminf W\left(f_{n}(t)\right) \mathrm{d} t \leqslant \liminf \int_{0}^{T} W\left(f_{n}(t)\right) \mathrm{d} t
The H^(1)H^{1} norm is weakly sequentially lower-semicontinuous (see Struwe, 1996), thus
||f^(˙)_(L^(2))^(2)=||f||_(H^(1))^(2)-||f||_(L^(2))^(2) <= l i m i n f||f_(n)||_(H^(1))^(2)-||f||_(L^(2))^(2)=l i m i n f||f_(n)^(˙)||_(L^(2))^(2),\left\|\dot{f}_{L^{2}}^{2}=\right\| f\left\|_{H^{1}}^{2}-\right\| f\left\|_{L^{2}}^{2} \leqslant \liminf \right\| f_{n}\left\|_{H^{1}}^{2}-\right\| f\left\|_{L^{2}}^{2}=\liminf \right\| \dot{f_{n}} \|_{L^{2}}^{2},
where the last equality holds because (f_(n))\left(f_{n}\right) converges strongly in L^(2)L^{2}. We evaluate now A_(T)(f)=||f||_(L^(2))^(2)//2+int_(0)^(T)W(f(t))dtA_{T}(f)=\|f\|_{L^{2}}^{2} / 2+\int_{0}^{T} W(f(t)) \mathrm{d} t, with WW given by (3):
A_(T)(f) <= l i m i n f||(f_(n)^(˙))||_(L^(2))^(2)+l i m i n fint_(0)^(T)W(f_(n)(t))dt <= l i m i n fA_(T)(f_(n)) <= a,A_{T}(f) \leqslant \liminf \left\|\dot{f_{n}}\right\|_{L^{2}}^{2}+\liminf \int_{0}^{T} W\left(f_{n}(t)\right) \mathrm{d} t \leqslant \liminf A_{T}\left(f_{n}\right) \leqslant a,
hence f inY_(a)f \in Y_{a}.
Proposition 14 implies that A_(T)A_{T} attains its infimum on XX; as a consequence of Palais' principle of symmetric criticality, and using Proposition 6, any bar(f)\bar{f} for which A_(T)( bar(f))=min{A_(T)(g):g in X}A_{T}(\bar{f})=\min \left\{A_{T}(g): g \in X\right\} is a classical periodic solution of the system (7).
Remark 17. Even if it does not appear explicitly in the proof of Theorem 16, the winding number kk is present in XX.
Remark 18. It is known that for each periodic solution of an autonomous system (as the anisotropic Schwarzschild-problem is) there exists a minimal period (see for example Amann, 1990). Let us consider a periodic solution ff of period TT and winding number kk, and let tau=T//m\tau=T / m be its minimal period (which leads to another value of the winding number). Applying one of the existence theorems for the period tau//2\tau / 2 we obtain the existence of a tau//2\tau / 2 periodic solution f_(1)f_{1}, which is of course also periodic of period TT, and surely different from ff (which has the minimal period tau\tau ). Continuing this process, we obtain an infinite set of distinct TT-periodic solutions (each one featured by its own winding number).
6. Conclusions
To summarize, here are some of the properties of the anisotropic Schwarzs-child-type problem revealed by applying abstract mathematical results:
6.1. The vector field associated to this anisotropic problem, expressed in configuration-momentum coordinates, exhibits seven symmetries S_(i),i= bar(1,7)S_{i}, i=\overline{1,7}, which, along with the identity, form a groups isomorphic to Z_(2)xxZ_(2)xxZ_(2)\mathbb{Z}_{2} \times \mathbb{Z}_{2} \times \mathbb{Z}_{2}.
6.2. Using variational methods, we point out the existence of S_(i)S_{i}-symmetric ( i=2,3i=2,3 ) periodic orbits that may have any assigned period TT and winding number k inZ\\{0}k \in \mathbb{Z} \backslash\{0\}.
6.3. Consider the anisotropic Schwarzschild-type problem to be a perturbation of the isotropic case (Stoica and Mioc, 1997) via the anisotropy parameter mu > 1\mu>1. Observe that anisotropy, no matter how large its size, deforms the S_(i)S_{i}-symmetric ( i=2,3i=2,3 ) periodic orbits of the isotropic problem (whose symmetries were pointed out by (Mioc, 2002)). but does not destroy them. This makes the symmetries S_(i)(i=2,3)S_{i}(i=2,3) constitute an indicator of the robustness of the system to perturbations.
These results add new, important features to the dynamics of the anisotropic Schwarzschild-type problem.
Acknowledgements
The authors are grateful to Professors George Contopoulos and Martin C. Gutzwiller for many suggestions intended to improve the paper.
References
Amann, H.: 1990, Ordinary Differential Equations: An Introduction to Nonlinear Analysis, Walter de Gruyer, Berlin, New York.
Ambrosetti, A. and Coti Zelati, V.: 1993, Periodic Solutions of Singular Lagrangian Systems, Progresses in Nonlinear Differential Equations and their Applications, No. 10, Birkhäuser, Boston.
Anisiu, M.-C.: 1998, Methods of Nonlinear Analysis Applied to Celestial Mechanics, Cluj University Press, Cluj-Napoca (Romanian).
Bertotti, M. L.: 1991, 'Forced oscillations of singular dynamical systems with an application to the restricted three-body problem', J. Diff. Eq. 93, 102-141.
Casasayas, J. and Llibre, J.: 1984, 'Qualitative analysis of the anisotropic Kepler problem', Mem. Amer. Math. Soc., Vol. 52, No. 312, AMS, Providence, RI.
Chenciner, A.: 2002, 'Action minimizing periodic orbits in the Newtonian nn-body problem', Celestial Mechanics (Evanston, IL, 1999), Contemporary Mathematics 292, Amer. Math. Soc., Providence, RI, pp. 71-90.
Chenciner, A. and Montgomery, R.: 2000, 'A remarkable periodic solution of the three-body problem in the case of equal masses', Ann. Math. 152, 881-901.
Coti Zelati, V.: 1994, Introduction to Variational Methods and Singular Lagrangian Systems, School and Workshop on Variational and Local Methods in the Study of Hamiltonian Systems, International Centre for Theoretical Physics, Trieste, Italy, 10-28 October 1994, SMR 779/4.
Craig, S., Diacu, F. N., Lacomba, E. A. and Perez, E.: 1999, 'On the anisotropic Manev problem', J. Math. Phys. 40, 1359-1375.
Devaney, R. L.: 1978, 'Collision orbits in the anisotropic Kepler problem', Invent. Math. 45, 221-251.
Devaney, R. L.: 1981, 'Singularities in classical mechanical systems', in Ergodic Theory and Dynamical Systems, Vol. 1, Birkhäuser, Boston, pp. 211-333.
Diacu, F. N.: 1996, 'Near-collision dynamics for particle systems with quasihomogeneous potentials', J. Diff. Eq. 128, 58-77.
Diacu, F. and Santoprete, M.: 2001, 'Nonintegrability and chaos in the anisotropic Manev problem', Physica D 156, 39-52.
Diacu, F. and Santoprete, M.: 2002, 'On the global dynamics of the anisotropic Manev problem', http://arXiv.org/abs/nlin/?0208012
Gordon, W. B.: 1975, 'Conservative dynamical systems involving strong forces', Trans. Amer. Math. Soc. 204, 113-135.
Gutzwiller, M. C.: 1971, 'Periodic orbits and classical quantization conditions', J. Math. Phys. 12, 343-358.
Gutzwiller, M. C.: 1973, 'The anisotropic Kepler problem in two dimensions', J. Math. Phys. 14, 139-152.
Gutzwiller, M. C.: 1977, 'Bernoulli sequences and trajectories in the anisotropic Kepler problem', J. Math. Phys. 18, 806-823.
Hénon, M. and Heiles, C.: 1964, 'The applicability of the third integral of motion: some numerical experiments', Astron. J. 69, 73-79.
McGehee, R.: 1973, 'A stable manifold theorem for degenerate fixed points with applications to celestial mechanics', J. Diff. Eq. 14, 70-88.
McGehee, R.: 1974, 'Triple collision in the collinear three-body problem', Invent. Math. 27, 191-227.
Mioc, V.: 2002, 'Symmetries in the Schwarzschild problem', Baltic Astron. 11, 393-407.
Mioc, V., Pérez-Chavela, E. and Stavinschi, M.: 2003, 'The anisotropic Schwarzschild-type problem, main features', Celest. Mech. Dyn. Astron. 86, 81-106 (Paper I).
Mioc, V. and Radu, E.: 1992, 'Orbits in an anisotropic radiation field', Astron. Nachr. 313, 353-357.
Palais, R.: 1979, 'The principle of symmetric criticality', Comm. Math. Phys. 69, 19-30.
Poincaré, H.: 1896, 'Sur les solutions périodiques et Ie principe de la moindre action', C. R. Acad. Sci. Paris 123, 915-918.
Santoprete, M.: 2002, 'Symmetric periodic solutions of the anisotropic Manev problem', J. Math. Phys. 43, 3207-3219.
Saslaw, W. C.: 1978, 'Motion around a source whose luminosity changes', Astrophys. J. 226, 240-252.
Schwarzschild, K.: 1916, 'Über das Gravitationsfeld eines Massenpunktes nach der Einsteinschen Theorie', Sitzber. Preuss. Akad. Wiss. Berlin, 189-196.
Stoica, C. and Mioc, V.: 1997, 'The Schwarzschild problem in astrophysics', Astrophys. Space Sci. 249, 161-173.
Struwe, M.: 1996, Variational Methods. Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems, Springer-Verlag, Berlin.
Tonelli, L.: 1915, 'Sur une méthode directe du calcul des variations', Rend. Circ. Mat. Palermo 39, 233-263.
Vinti, J. P.: 1972, 'Possible effects of anisotropy of GG on celestial orbits', Celest. Mech. 6, 198-207.
Will, C.: 1971, 'Relativistic gravity in the solar system. II. Anisotropy in the Newtonian gravitational constant', Astrophys. J. 169, 141-155.