The existence of chaotic behavior for the geodesics of the test particles orbiting compact objects is a subject of much current research. Some years ago, Gu\'eron and Letelier [Phys. Rev. E \textbf{66}, 046611 (2002)] reported the existence of chaotic behavior for the geodesics of the test particles orbiting compact objects like black holes induced by specific values of the quadrupolar deformation of the source using as models the Erez--Rosen solution and the Kerr black hole deformed by an internal multipole term. In this work, we are interesting in the study of the dynamic behavior of geodesics around astrophysical objects with intrinsic quadrupolar deformation or nonisotropic stresses, which induces nonvanishing quadrupolar deformation for the nonrotating limit. For our purpose, we use the Tomimatsu-Sato spacetime [Phys. Rev. Lett. \textbf{29} 1344 (1972)] and its arbitrary deformed generalization obtained as the particular vacuum case of the five parametric solution of Manko et al [Phys. Rev. D 62, 044048 (2000)], characterizing the geodesic dynamics throughout the Poincar\'e sections method. In contrast to the results by Gu\'eron and Letelier we find chaotic motion for oblate deformations instead of prolate deformations. It opens the possibility that the particles forming the accretion disk around a large variety of different astrophysical bodies (nonprolate, e.g., neutron stars) could exhibit chaotic dynamics. We also conjecture that the existence of an arbitrary deformation parameter is necessary for the existence of chaotic dynamics.Comment: 7 pages, 5 figure

In the present paper, using the first-order approximation of the n-body Lagrangian (derived on the basis of the post-Newtonian gravitational theory of Einstein, Infeld, and Hoffman), we explicitly write down the equations of motion for the planar circular restricted three-body problem. Additionally, with some simplified assumptions, we obtain two formulas for estimating the values of the mass-distance and velocity-speed of light ratios appropriate for a given post-Newtonian approximation. We show that the formulas derived in the present study, lead to a good numerical conservation of the Jacobi constant and allow for an approximate equivalence between the Lagrangian and Hamiltonian approaches at the same post-Newtonian order. Accordingly, the dynamics of the system is analyzed in terms of the Poincar sections method and Lyapunov exponents, finding that for specific values of the Jacobi constant the dynamics can be either chaotic or regular. Our results suggest that the chaoticity of the post-Newtonian system is slightly in-creased in comparison with its Newtonian counterpart.

In 1998, Shibata and Sasaki [Phys. Rev. D 58, 104011 (1998)] presented an approximate analytical formula for the radius of the innermost stable circular orbit (ISCO) of a neutral test particle around a massive, rotating and deformed source. In the present paper, we generalize their expression by including the magnetic dipole moment. We show that our approximate analytical formulas are accurate enough by comparing them with the six-parametric exact solution calculated by Pachón et. al. [Phys. Rev. D 73, 104038 (2006)] along with the numerical data presented by Berti and Stergioulas [MNRAS 350, 1416[MNRAS 350, (2004] for realistic neutron stars. As a main result, we find that in general, the radius at ISCO exhibits a decreasing behavior with increasing magnetic field. However, for magnetic fields below 100GT the variation of the radius at ISCO is negligible and hence the non-magnetized approximate expression can be used. In addition, we derive approximate analytical formulas for angular velocity, energy and angular momentum of the test particle at ISCO.