Periodic Orbits in the Photogravitational Restricted Problem When the Primaries Are Triaxial Rigid Bodies

We have studied periodic orbits generated by Lagrangian solutions of the restricted three-body problem when both the primaries are triaxial rigid bodies and source of radiation pressure. We have determined periodic orbits for different values of ′ ′ ′ 1 1 2 2 , , , , , , and μ h A A A A P P (h is energy constant; μ is mass ratio of the two primaries; 2 , , ′ ′ 1 1 2 and A A A A are parameters of triaxial rigid bodies and ′ and P P are radiation parameters). These orbits have been determined by giving displacements along the tangent and normal at the mobile co-ordinates as defined in our papers (Mittal et al. [1]-[3]). These orbits have been drawn by using the predictor-corrector method. We have also studied the effect of triaxial bodies and source of radiation pressure on the periodic orbits by taking fixed value of μ.


Introduction
This paper is the extension of our papers, Mittal et al. [1]- [3].Charlier [4] and Plummer [5] studied the existence of two families of small periodic motions in the neighborhood of Lagrangian solutions in the plane circular restricted three-body problem, with different values of the parameter μ.Riabov [6] investigated periodic motions analytically.Szebehely [7] has described the results on the periodic motions of circular restricted three-body problem.Deprit and Henrard [8] gave more results on periodic motions in their paper.Markeev and Sokolsky [9] worked on the small periodic motions generated by Lagrangian solutions for all values of μ and for small values of h for which the conditions of holomorphic integral theorem are valid.Hadjidemetriou [10] has discussed the continuation of periodic orbits from the restricted to the general three-body problem.Karimov and Sokolsky [11] have studied the periodic motions generated by Lagrangian solutions of the circular restricted three-body problem by using mobile co-ordinates and by taking displacements along the tangent and the normal.Aggarwal et al. [12] have discussed the non-linear stability of the triangular libration point when both the primaries are radiated axe symmetric rigid bodies in the presence of third and forth order resonance.Abouelmagd et al. [13] have studied the periodic structure of the restricted three-body problem considering the effect of the zonal harmonics J 2 and J 4 for the more massive body.They showed that the triangular points in the restricted three-body problem have long or short periodic orbits in the range 0 ≤ µ < µ c .Perdios et al. [14] have studied the equilibrium points and related periodic motions in the restricted three-body problem with angular velocity and radiation effects.Jain and Aggarwal [15] have studied the stability and existence of non-collinear libration points in restricted three-body problem with stokes drag effect when smaller primary is an oblate spheroid.
The celestial bodies are in general axis-symmetric bodies, so its shape should be taken into account as well.The replacement of mass point by rigid-body is quite important because of its wide applications.The re-entry of artificial satellite has shown the importance of periodic orbits.
That is why, we have thought of studying, in this paper, the periodic orbits generated by Lagrangian solutions of the restricted three-body problem when both the primaries are triaxial rigid bodies and source of radiation pressure.We determine the periodic orbits by giving displacements at the mobile co-ordinates along the tangent and normal.We have also determined family of periodic orbits by fixing μ (mass ratio of the two primaries) and changing the values of 1 2 1 2 , , , A A A A ′ ′ (parameters of the triaxial rigid bodies), P and P ′ (the radiation parame- ters) and varying h (energy constant).We have also studied the effect of triaxial parameters and the radiation pressure on the energy constant (h).
Most of the authors have not taken into account the effect of the solar radiation pressure in the motion of the third body whereas we have taken both the primaries as radiating triaxial rigid bodies.Besides taking both the primaries as triaxial rigid bodies and the source of radiation, we have used mobile-coordinates and given the displacement along the normal and the tangent to the orbit which has wider applications in space dynamics.We have drawn the periodic orbits by using the predictor-corrector method which is given in detail in our papers [1]- [3].

Equations of Motion
Following the procedure of our papers [1]- [3], we consider three masses Let us consider a synodic system of coordinates O(XYZ), initially coincident with the inertial system O(XYZ), rotating with the angular velocity n about Z-axis; (the Z-axis is coincident with Z-axis).We choose unit of mass such that 1 2 1, m m + = the unit of distance 1 AB = and unit of time is so chosen that 1 G = .Using the dimensionless variables, we find the Lagrangian function L and the eq- uations of motion of the infinitesimal mass in the restricted three-body problem when both the primaries are radiating triaxial rigid bodies in the synodic co-ordinate system (Figure 1).
Equations of motion with Lagrangian function L are given by Radiation pressure due to the bigger primary Gravitational force due to the bigger primary P = , Radiation pressure due to the smaller primary Gravitational force due to the smaller primary a, b, c = length of the semi axes of the triaxial body of mass 1 m , , , a b c ′ ′ ′ = length of the semi axes of the triaxial body of mass 2 m , R = dimensional distance between the primaries, ( ) ( ) U = constant to be so chosen such that h (energy constant) will vanish at 4 L (libration point).The coordinates of 4 L (libration) are where 1 1 2 Equations of motion can also be written as where The Jacobi integral is
The corresponding differential equations are given by the system of Equation ( 2) with Jacobi integral given by (3).We consider the solutions of Equations ( 2) for which C is zero.If we consider the solutions of Equation (2) given by ( 4) for some fixed parameters value p then there may exist another solution given by ( 5) with another parameter value ( ) , , , , , , , We have , , , , , , , , ( )

, x x h A A A A P P y y h A A A A P P x x h A A A A P P y y h A
and We consider p ∆ and ξ as small quantities of the same order.Then have the following variational equ- ations: with the integral constructed from Equation ( 3), retaining the first order terms only, we get .
The modulus of momentary velocity on the orbit is defined as ( ) ( ) V t Q t x y = = + .We assume that (5) is not corresponding to the equilibrium state, i.e., ( ) 0 V t ≠ and we further assume that ( ) 0 V t ≠ on the whole orbit.Therefore, x and y become the mobile co-ordinates.We will, now, use the mobile coordinate system to draw the periodic orbits by resolving one of the axis along the velocity vector

Y t r t
V t = = the unit vector which is normal to the orbit, i.e., it is orthogonal to the vector s(t).

X t s t
V t = = the unit vector which is tangent to the orbit and is the last column of the matrix S.
So, we have It can be easily verified that ( ) 1, 0 and 0,1 s s r s S s e It may be noted that, Now, we may further define , the last line of S −1 .
We write α , the vector of local coordinates, in the new coordinate system as follows: ( ) ( ) where N is displacement along the normal to the orbit and M is displacement along the tangent to the orbit.
Then, the new coordinates are given by ( ) Substituting these values into the integral (7a), we have Equation ( 10) can be solved for M as ( ) ( ) Equations of motion (2) for the new coordinates are -1 2 or where ( ) 0 and 0 and N F and N F are the same as in our papers [1]- [3].
, therefore the equations of motion in normal and tangent co-ordinates can be written as ( ) and ( ) N can also be written as (

N W y W xy W x W y W x nV xy xy W V n xy xy xy xy nV xy xy x y V N W
Thus, we have derived the equation in N which possesses the remarkable property that the normal coordi- nate (N) is independent of the tangent coordinate (M).Moreover, instead of the differential equation of the second order (13), we can use the first order differential Equation (11).If the investigated motion ( 4) is periodic, then the matrix S(t) can be taken as periodic and Equations ( 12) and ( 13) will have the periodic coefficients at 0 p ∆ ≡ .

Periodic Orbits
For determining the periodic orbits, the required equations of motion and the variational equations are given as: 2 , 0 0 2 where ( ) Z t is the matrix of solutions of a homogeneous system with initial condition ( ) 2 are the solutions of Cauchy problem (iv) with (vi) of ( 14).The order of the above system is thirty-four.
After solving the above equations of motion (i) and (ii), the variational Equations (iii)-(vi) of ( 14) and applying the predictor-corrector method, we have determined the periodic orbits.
The above analysis is summed up in Table 1.By taking both the primaries as radiating triaxial rigid bodies, the difference in the behavior of the values of h is obvious.

Conclusions
Karimov and Sokolsky [11] have studied periodic orbits in the restricted three body problem by giving the displacements along the normal and the tangent to the orbit at the mobile co-ordinates.They have taken both the primaries as point masses while in this paper besides taking both the primaries as triaxial rigid bodies, we have also taken both the primaries as source of radiation pressure as well.In this paper, we have again determined   five periodic orbits in a family for fixed value of the mass parameter μ, the triaxial parameters 1 2 1 2 , , , A A A A ′ ′ and the radiation parameters P and P′ with varying energy constant h.
We have observed the following effects on the periodic orbits and on the energy constant h due to triaxial rigid bodies and radiation pressure if we compare it with the results of Karimov and Sokolsky [11] and our papers [1]- [3]. 1) The energy constant h increases in a family (for Figure 2 and Figure 3) then it decreases (for Figures 4-6) for fixed trixial parameters 1 2 1 2 , , , A A A A ′ ′ and radiation parameters P and P′ .2) As we increase the radiation parameters P and P′ , the energy constant h increases whereas the periodic or- bits shrink a little.3) The periodic orbits go away from the libration point 4 L as we increase triaxial parameters 1 2 1 2 , , , A A A A ′ ′ and radiation parameters P and P′ whereas energy constant h decrease.We have investigated the family up to the member which touches the point 4 L .It is observed that the families of periodic orbits in Karimov and Sokolsky [11] terminate at both the triangular equilibrium points simultaneously, while in our case these families are non-symmetrical, so they may continue.

m
revolve with angular velocity n (say) in circular orbits without rotation about the centre of mass O. Let there be an infinitesimal mass 3 m which is moving in the plane of motion of 1 2 and m m and is being influenced by their motion but not influencing them.Let the line joining 1 2 and m m be taken as X-axis and O their centre of mass as origin.Let the line through O and perpendicular to OX, and lying in the plane of motion of 1 2 and m m being Y-axis.

Figure 1 . 1
Figure 1.Configuration of the restricted three-body problem with masses m 1 and m 2 as radiating triaxial rigid bodies.
coordinate system, we consider the transition matrix S as follows: Consider the first column of S as the equations with zero initial conditions, i.e.,

Figure 2 3 4 5 6
Triaxial body Figure Triaxial body Figure Triaxial body Figure Triaxial body Figure Values

Table 1 .
The summary analysis of the periodic orbits.