Analytic and numerical treatment of motion of dust grain particle around triangular equilibrium points with post-AGB binary star and disc

This paper investigates the motion around the triangular equilibrium points, of a passively gravitating dust particle in the gravitational field of a low-mass post-AGB binary system, surrounded by circumbinary disc. The two bodies of the binary are modeled as a triaxial star and a radiating-oblate star. Due to small deviation of disc stars on circular orbits, we have assumed that the Coriolis and centrifugal forces of the stars are slightly perturbed. The triangular equilibrium points of the particle are found. These points are defined by, triaxiality of the primary star, oblateness and radiation of the secondary one and the gravitational potential from the disc mass. Further, when the disc mass increases, the particle moves nearer to the stars and farther away from the disc. In general, these equilibrium points are linearly stable when μ < μ_C; where μ is the mass ratio and μ_C is the critical mass function, defined by the parameters of the system. The effects of each of these parameters on the size of the stability region are stated, and the periodic motion around the stable points is examined. It is seen that the orbits are ellipses, and the orientation, eccentricities, lengths of the semi-major and semi-minor axes are influenced by the parameters of the problem. In particular, for our numerical linear stability analysis, we have taken an extremely depleted pulsating star, IRAS 11472-0800 as the post-AGB triaxial star, with a weakly-radiating young white dwarf star; G29-38 as the secondary. For this system, the stability result of the triangular points comes out different. Here, μ_C < μ throughout the entire range of the mass ratio and the critical mass function. Hence, the triangular equilibrium points are unstable. The stability of the orbits is tested using the Poincaré surfaces of section (Pss). The region of stability is controlled by the introduced parameters and the Jacobi constant.