Periodic orbits and escapes in dynamical systems

We study the periodic orbits and the escapes in two different dynamical systems, namely (1) a classical system of two coupled oscillators, and (2) the Manko-Novikov metric which is a perturbation of the Kerr metric (a general relativistic system). We find their simple periodic orbits, their characteristics and their stability. Then we find their ordered and chaotic domains. As the energy goes beyond the escape energy, most chaotic orbits escape. In the first case we consider escapes to infinity, while in the second case we emphasize escapes to the central ??bumpy?? black hole. When the energy reaches its escape value, a particular family of periodic orbits reaches an infinite period and then the family disappears (the orbit escapes). As this family approaches termination it undergoes an infinity of equal period and double period bifurcations at transitions from stability to instability and vice versa. The bifurcating families continue to exist beyond the escape energy. We study the forms of the phase space for various energies, and the statistics of the chaotic and escaping orbits. The proportion of these orbits increases abruptly as the energy goes beyond the escape energy.     

磁重联电流片的参数变化对电子加速的影响

通过采用试验粒子的方法,研究了在有引导磁场Bz存在的磁重联电流片中,电子被super-Dreicer电场Ez加速后的运动特征.首先,考虑了引导磁场恒定且与电场有不同方向时对粒子加速的影响.在这种情况下,Bz方向的改变直接改变了电子的运动轨迹,使其沿着不同的路径离开电流片.在Bz和Ez同向时,高能电子的pitch-angle接近于180°.然而,当2者反向时,高能电子的pitch-angle接近0°.引导磁场的取向只是使电场有选择地对不同区域的电子进行加速,不会最终影响电子的能量分布,最终得到的能谱是普遍的幂率谱E-γ.在典型的日冕条件下, γ大约等于2.9.进一步的研究表明γ的大小依赖于引导磁场及磁重联电场的强弱,以及电流片的尺度.随后,也研究了包含多个X-点和O-点电流片中被加速粒子的运动特征.结果表明X-点和O-点的存在使得粒子被束缚在加速区并获得最大的加速,而且最终的能谱具有多幂率谱的特征.     

Analytical description of charged particle motion in a reconnecting current sheet

We have obtained an analytical solution to the equation of motion in the guiding center approximation for nonrelativistic charged particles in a reconnecting current sheet with a three-component magnetic field. Given the electric field attributable to magnetic reconnection, the solution describes stable and unstable three-dimensional particle orbits. We have found the domain of input parameters at which the motion is stable. A physical interpretation of the processes affecting the stability of the motion is given. Charge separation is shown to take place in the sheet during the motion: oppositely charged particles are localized mostly in different regions of the current sheet. A formula is derived for the particle energy in stable and unstable orbits. The results obtained by numerical and analytical methods are compared.     

Chaos-induced resistivity of collisionless magnetic reconnection in the presence of a guide field

One of the most puzzling problems in astrophysics is to understand the anomalous resistivity in collisionless magnetic reconnection that is believed extensively to be responsible for the energy release in various eruptive phenomena. The magnetic null point in the reconnecting current sheet, acting as a scattering center, can lead to chaotic motions of particles in the current sheet, which is one of the possible mechanisms for anomalous resistivity and is called chaos-induced resistivity. In many interesting cases, however, instead of the magnetic null point, there is a nonzero magnetic field perpendicular to the merging field lines, usually called the guide field, whose effect on chaos-induced resistivity has been an open problem. By use of the test particle simulation method and statistical analysis, we investigate chaos-induced resistivity in the presence of a constant guide field. The characteristics of particle motion in the reconnecting region, in particular, the chaotic behavior of particle orbits and evolving statistical features, are analyzed. The results show that as the guide field increases, the radius of the chaos region increases and the Lyapunov index decreases. However, the effective collision frequency, and hence the chaos-induced resistivity, reach their peak values when the guide field approaches half of the characteristic strength of the reconnection magnetic field. The presence of a guide field can significantly influence the chaos of the particle orbits and hence the chaos-induced resistivity in the reconnection sheet, which decides the collisionless reconnection rate. The present result is helpful for us to understand the microphysics of anomalous resistivity in collisionless reconnection with a guide field.     

Orbital structure in N-body models of barred-spiral galaxies

We study the orbital structure in a series of self-consistent N -body configurations simulating rotating barred galaxies with spiral and ring structures. We perform frequency analysis in order to measure the angular and the radial frequencies of the orbits at two different time snapshots during the evolution of each N -body system. The analysis is done separately for the regular and the chaotic orbits. We thereby identify the various types of orbits, determine the shape and percentages of the orbits supporting the bar and the ring/spiral structures, and study how the latter quantities change during the secular evolution of each system. Although the frequency maps of the chaotic orbits are scattered, we can still identify concentrations around resonances. We give the distributions of frequencies of the most important populations of orbits. We explore the phase-space structure of each system using projections of the 4D surfaces of section. These are obtained via the numerical integration not only of the orbits of test particles, but also of the real N -body particles. We thus identify which domains of the phase space are preferred and which are avoided by the real particles. The chaotic orbits are found to play a major role in supporting the shape of the outer envelope of the bar as well as the rings and the spiral arms formed outside corotation.     

Stickiness effects in chaos

Stickiness is a temporary confinement of orbits in a particular region of the phase space before they diffuse to a larger region. In a system of 2-degrees of freedom there are two main types of stickiness (a) stickiness around an island of stability, which is surrounded by cantori with small holes, and (b) stickiness close to the unstable asymptotic curves of unstable periodic orbits, that extend to large distances in the chaotic sea. We consider various factors that affect the time scale of stickiness due to cantori. The overall stickiness (stickiness of the second type) is maximum near the unstable asymptotic curves. An important application of stickiness is in the outer spiral arms of strong-barred spiral galaxies. These spiral arms consist mainly of sticky chaotic orbits. Such orbits may escape to large distances, or to infinity, but because of stickiness they support the spiral arms for very long times.     

Planar periodic orbits in exterior resonances with Neptune

In the framework of the planar restricted three-body problem we study a considerable number of resonances associated to the basic dynamical features of Kuiper belt and located between 30 and 48 a.u. Our study is based on the computation of resonant periodic orbits and their stability. Stable periodic orbits are surrounded by regular librations in phase space and in such domains the capture of trans-Neptunian object is possible. All the periodic orbits found are symmetric and there is an indication of the existence of asymmetric ones only in a few cases. In the present work first, second and third order resonances are under consideration. In the planar circular case we found that most of the periodic orbits are stable. The families of periodic orbits are temporarily interrupted by collisions but they continue up to relatively large values of the Jacobi constant and highly eccentric regular motion exists for all cases. In the elliptic problem and for a particular eccentricity value of the primary bodies, the periodic orbits are isolated. The corresponding families, where they belong to, bifurcate from specific periodic orbits of the circular problem and seem to continue up to the rectilinear problem. Both stable and unstable orbits are obtained for each case. In the elliptic problem, the unstable orbits found are associated with narrow chaotic domains in phase space. The evolution of the orbits, which are located in such chaotic domains, seems to be practically regular and bounded for long time intervals.     

The Photogravitational Hill Problem: Numerical Exploration

We consider the recently introduced version of the classical Lunar Hill problem, the photogravitational Hill problem, and study it's equilibrium points and zero-velocity curves. The full network of families of periodic orbits is numerically explored, their stability is computed and critical orbits are determined. Non-periodic orbits are also computed as points on a surface of section, providing an outlook of the stability regions, chaotic motions and escape.     

Acceleration of Electrons and Protons in Reconnecting Current Sheets Including Single or Multiple X-points

Kinematic characteristics of electrons and protons in the magnetic reconnecting current sheet in the presence of a guide field are investigated. Particle trajectories are calculated for different values of the guide field by a test-particle calculation. The relationship between the final energy and the initial position has also been studied. We found that the addition of a guide field not only allows particles to get more energy and not only results in the separation of electrons and protons, but also causes the reconnecting electric field to selectively accelerate electrons and protons for different initial positions. The energy spectrum eventually obtained is the common power-law spectrum, and as the guide field increases, the index for the spectrum of electrons decreases rapidly. However, for a weak background magnetic field, proton spectra are not very sensitive to the guide field; but for a strong background field, the dependence of the spectrum index is similar to the electron spectrum. Meanwhile, kinematic characteristics of the accelerated particles in the current sheet including multiple X-points and O-points were also investigated. The result indicates that the existence of the multiple X- and O-points helps particles trapped in the accelerating region to gain more energy, and yields the double or multiple power-law feature.     

Generalized adiabatic theory applied to the magnetotail current sheet

We use the generalized first adiabatic invariant, an extension of the magnetic moment for regions of large field gradients, to treat particles in the magnetotail current sheet. The equations of motion can be expressed in terms of drift parameters which vary slowly and smoothly at the drift rate, not at the gyration rate. The analysis leads to boundaries in phase space which form a generalized loss cone and separate particles drifting into and out of the layer from particles trapped within the layer. These boundaries can be used in the moment integrals for densities and currents when the drifting particles differ in temperature, or in other properties, from the trapped population, as has been suggested by observations. We give examples of how different kinds of particle orbits contribute to the spatial profiles of density and current and thus to the field structure of the current sheet. We find that the parallel pressure of the drifting particles must exceed the transverse pressure for self-consistent solutions to exist, and based on this result, we give examples of fully self-consistent solutions using bi-Maxwellian ion and Maxwellian electron distributions. We give a proof, using generalized adiabatic theory, of Cowley's (1978a) theorem that particles trapped in the current layer experience zero net drift.Paper dedicated to Professor Hannes Alfvén on the occasion of his 80th birthday, 30 May 1988.     

Particle Acceleration by a Time-Varying Electric Field in Merging Magnetic Fields

Traditional models for particle acceleration by magnetic reconnection in solar flares assumed a constant electric field in a steady reconnecting magnetic field. Although this assumption may be justified during the gradual phase of flares, the situation is different during the impulsive phase. Observed rapid variations in flare emissions imply that reconnection is non-steady and a time-varying electric field is present in a reconnecting current sheet. This paper describes exploratory calculations of charged particle orbits in an oscillating electric field present either at a neutral plane or a neutral line of two-dimensional magnetic field. A simple analytical model makes it possible to explain the effects of particle trapping and resonant acceleration previously noted by Petkaki and MacKinnon in a numerical simulation. As an application, electron acceleration to X-ray generating energies in impulsive solar flares is discussed within the context of the model.     

Global Dynamics in Self-Consistent Models of Elliptical Galaxies

We compare two different N-body models simulating elliptical galaxies. Namely, the first model is a non-rotating triaxial N-body equilibrium model with smooth center, called SC model. The second model, called CM model, is derived from the SC by inserting a central mass in it, so that all possible differences between the two models are due to the effect of the central mass. The central mass is assumed to be mainly due to a massive central black hole of mass about 1% of the total mass of the galaxy. By using the fundamental frequency analysis, the two systems are thoroughly investigated as regards the types of orbits described either by test particles, or by the real particles of the systems at all the energy levels. A comparison between the orbits of test particles and the orbits of real particles at various energy levels is made on the rotation number plane. We find that extensive stable regions of phase space, detected by test particles remain empty, i.e. these regions are not occupied by real particles, while many real particles move in unstable regions of phase space describing chaotic orbits. We run self-consistently the two models for more than a Hubble time. During this run, in spite of the noise due to small variations of the potential, the SC model maintains (within a small uncertainly) the number of particles moving on orbits of each particular type. In contrast, the CM model is unstable, due to the large amount of mass in chaotic motion caused by the central mass. This system undergoes a secular evolution towards an equilibrium state. During this evolution it is gradually self-organized by converting chaotic orbits to ordered orbits mainly of the short axis tube type approaching an oblate spheroidal equilibrium. This is clearly demonstrated in terms of the fundamental frequencies of the orbits on the rotation number plane and the time evolution of the triaxiality index.     

Regular and chaotic orbits in a self-consistent triaxial stellar system with slow figure rotation

We created a self-consistent triaxial stellar system through the cold disipationless collapse of 100,000 particles whose evolution was followed with a multipolar code. The resulting system rotates slowly even though its total angular momentum is zero, i.e., it offers an example of figure rotation. The potential of the system was subsequently approximated with interpolating formulae yielding a smooth potential stationary in the rotating frame. The Lyapunov exponents could then be computed for a randomly selected sample of 3,472 of the bodies that make up the system, allowing the recognition of regular and partially and fully chaotic orbits. The regular orbits were Fourier analyzed and classified using their locations on the frequency map. A comparison with a similar non-rotating model showed that the fraction of chaotic orbits is slightly but significantly enhanced in the rotating model; alternatively, there are no significant differences between the corresponding fractions neither of partially and fully chaotic orbits nor of long axis tubes, short axis tubes, boxes and boxlets among the regular orbits. This is a reasonable result because the rotation causes a breaking of the symmetry that may increase chaotic effects, but the rotation velocity is probably too small to produce any other significant differences. The increase in the fraction of chaotic orbits in the rotating system seems to be due mainly to the effect of the Coriolis force, rather than the centrifugal force, in good agreement with the results of other investigations.     

Regular and irregular periodic orbits

We distinguish between regular orbits, that bifurcate from the main families of periodic orbits (those that exist also in the unperturbed case) and irregular periodic orbits, that are independent of the above. The genuine irregular families cannot be made to join the regular families by changing some parameters. We present evidence that all irregular families appear inside lobes formed by the asymptotic curves of the unstable periodic orbits. We study in particular a dynamical system of two degrees of freedom, that is symmetric with respect to the x-axis, and has also a triple resonance in its unperturbed form. The distribution of the periodic orbits (points on a Poincaré surface of section) shows some conspicuous lines composed of points of different multiplicities. The regular periodic orbits along these lines belong to Farey trees. But there are also lines composed mainly of irregular orbits. These are images of the x-axis in the map defined on the Poincaré surface of section. Higher order iterations of this map , close to the unstable triple periodic orbit, produce lines that are close to the asymptotic curves of this unstable orbit. The homoclinic tangle, formed by these asymptotic curves, contains many regular orbits, that were generated by bifurcation from the central orbit, but were trapped inside the tangle as the perturbation increased. We found some stable periodic orbits inside the homoclinic tangle, both regular and irregular. This proves that the homoclinic tangle is not completely chaotic, but contains gaps (islands of stability) filled with KAM curves.     

Phase-bunching and other non-linear processes occurring in gyroresonant wave-particle interactions

We consider the movement of individual electrons in a magnetized plasma in which a monochromatic wave is propagating in the whistler mode. We derive simple expressions which give the displacement of the electrons as a function of time, the phase angle that their velocity vector makes with the magnetic component of the wave, their pitch angle and energy changes. A useful formula is obtained which gives the velocity range over which particles remain trapped inside the wave, as a function of the wave intensity and of the initial phase angle of the particle. It is shown that even strictly resonant particles can escape from the wave when their initial phase angle is very small. From the derived expressions, it is possible to compute the phase-bunching effect which occurs approximately at one trapping wavelength behind the leading edge of the interaction region. We deduce also the total amount of energy which is taken from (or given to) the wave by magnetospheric electrons in both cases of naturally existing or artificially injected particles. It is shown that these non-linear amplification processes can lead to very large VLF amplitude in the magnetosphere.     

Spatial Structure of Regular and Chaotic Orbits in A Self-Consistent Triaxial Stellar System

We created a triaxial stellar system through the cold dissipationless collapse of 100,000 particles whose evolution was followed with a multipolar code. Once an equilibrium system had been obtained, the multipolar expansion was freezed and smoothed in order to get a stationary smooth potential. The resulting model was self-consistent and the orbits and Lyapunov exponents could then be computed for a randomly selected sample of 3472 of the bodies that make up the system. More than half of the orbits (52.7 % ) turned out to be chaotic. Regular orbits were then classified using the frequency analysis automatic code of Carpintero and Aguilar (1998, MNRAS 298(1), 1–21). We present plots of the distributions of the different kinds of orbits projected on the symmetry planes of the system. We distinguish chaotic orbits with only one non-zero Lyapunov exponent from those with two non-zero exponents and show that their spatial distributions differ, that of the former being more similar to the one of the regular orbits. Most of the regular orbits are boxes and boxlets, but the minor axis tubes play an important role filling in the wasp waists of the boxes and helping to give a lentil shape to the system. We see no problem in building stable triaxial models with substantial amounts of chaotic orbits; the difficulties found by other authors may be due not to a physical cause but to a limitation of Schwarzschild’s method.     

Energy spectra of particles accelerated in a reconnecting current sheet with the guiding magnetic field

Electron and proton acceleration by a super-Dreicer electric field is investigated in the non-neutral reconnecting current sheet (RCS) with a non-zero longitudinal component of the magnetic field ('guiding field'). The guiding field is assumed parallel to the direction of electric field and constant within an RCS. The other two magnetic field components, transverse and tangential, are considered to vary with distances from the X null point of an RCS. The proton and electron energy spectra are calculated numerically from a motion equation using the test particle approach for model RCSs with constant and variable densities. In the presence of a strong or moderate guiding field, protons were found fully or partially separated from electrons at ejection from an RCS into the opposite, 'electron' and 'proton', semiplanes. In the case of a weak guiding field, both protons and electrons are ejected symmetrically in equal proportions as neutral beams. The particles ejected from an RCS with a very weak or very strong guiding field have power-law energy spectra with spectral indices of about 1.5 for protons and 2.0 for electrons. For a moderate guiding field, the energy spectra of electrons ejected into the opposite semiplanes are mixed, i.e. in the 'electron-dominated' semiplane power-law energy spectra for electrons and thermal-like for protons, while in the 'proton' semiplane they are symmetrically mirrored.     

Chaos in Barred Galaxy Models

We investigate the regular and chaotic motion in a model potential found using the recent developments of the Inverse Problem of Dynamics. The potential describes the motion in the central parts of a barred galaxy. In the absence of rotation chaotic motion is observed when the perturbation strength is near the escape perturbation for a fixed value of the energy. In the rotating cases one observes that the area of chaotic motion on the surface of section decreases as the angular velocity Ω increases and finally all orbits become regular. The character of motion is also checked by computing the Liapunov characteristic exponents in all cases. This revised version was published online in July 2006 with corrections to the Cover Date.     

Particle distributions and X-ray spectra in single or multiple solar current sheets

The acceleration of charged particles in a site of magnetic reconnection is analysed by detailed numerical simulations. Single or multiple encounters of the particles with Harris-type reconnecting current sheets (RCSs) are modelled as an overall stochastic process taking place within an active region. RCS physical parameters are selected in a parameter space relevant to solar flares. Initially, the charged particles form a thermal (Maxwellian) distribution corresponding to coronal temperature  ≃2 × 10⁶ K  . Our main goal is to investigate how the acceleration process changes the shape of the particles' kinetic energy distribution. The evolution of the kinetic energy distribution, calculated numerically after one encounter of the particles with a single RCS, is found to be in good agreement with our previously published analytical formulae. In the case of consecutive encounters, we find that the kinetic distribution tends to converge to a practically invariant form after a relatively small number of encounters. We construct a discrete stochastic process that reproduces the numerical distributions and we provide a theoretical interpretation of the asymptotic convergence of the energy distribution. We finally compute the theoretical X-ray spectra that would be emitted by the simulated particles in a thick target model of radiation.