## 1 INTRODUCTION

In the last two decades, studies related to exoplanets have become an important field in modern astrophysics, as observational instrumentations have become more sensitive and efficient for detecting signals related to this type of planet. While there are more than 100 cataloged exoplanets orbiting single star, the studies of such planets, including detection, characterization, and evolution, are still difficult. In addition, recent discoveries confirm the presence of an exoplanet around a binary companion (Borucki et al. 2010). Therefore, many stars in the solar neighborhood are a member of binary or multiple systems, which will force us to modify our knowledge regarding the dynamics of the planets in binary star systems. The dynamics of planets in binary systems are not easy to formulate because of the second star perturbations, which require researchers to neglect other perturbations like star radiation. In a hierarchical triple system, the motion of a planet will be constrained to revolve around the center of the binary star as a two-body problem.

Systems with multiple stars and planetary components draw the attention of researchers with their rich dynamics that can carry information about planetary formation and evolution.

Raghavan et al. (2010) found approximately half of the stars in the Sun neighborhood in multiple systems. At about a few tens of AU, these local systems have planets like those hosted by single stars, presenting interactions between stellar and planetary dynamics. Marzari et al. (2005) focused on the planetary systems subjected to planet–planet scattering around a binary star system. They studied the orbits of the unejected planets. Triple-body systems in which a planet revolves around a binary star system can be treated as hierarchical stable systems, and the motion can be described by two related systems. Each system describes a two-body problem: the motion of the inner binary stars and the motion of the planet with the center of mass of the inner system. The hierarchical system dynamics, especially the long period or secular terms, attracted the attention of researchers such as Harrington (1968), Heppenheimer (1978), Ford et al. (2000), Lee & Peale (2003), Farago & Laskar (2010), Leung & Lee (2013), Naoz et al. (2013a, 2013b), Li et al. (2014), and Liu et al. (2015). As any gravitational perturbations in the binary system will strongly affect the motion of the planet, many studies introduced the dynamic stability of planets, especially the perturbation due to a secondary star (Boss 2006; Haghighipour 2006; Kley & Nelson 2008; Thébault et al. 2009; Giuppone et al. 2011). Carvalho et al. (2016) introduced triple system secular dynamics assuming a center star like the Sun and a planet like Jupiter all perturbed by another brown dwarf star.

In this work, the long period dynamics of the exoplanet will be formulated taking into consideration different rates of mass loss from the second inner companion using a hierarchical triple system to formulate the equations.

## 2 PROBLEM MODEL

In the computations, we utilized a triple body system where an outer body or exoplanet of mass, *m*_{2}, orbited the center of mass of two inner bodies or binary stars of mass (*m*_{0} and *m*_{1}). Each system was described by two perturbed Keplerian bodies.

Using Jacobi coordinates, ${\overrightarrow{r}}_{1}$ and ${\overrightarrow{r}}_{2}$ is the position vector of *m*_{1} to *m*_{0}, and ${\overrightarrow{r}}_{2}$ is the position vector of *m*_{2} to the center of mass of the binary stars, respectively (see Fig. 1). In addition to well-known Keplerian elements, the semi-major axis *a _{j}*, the eccentricity

*e*, the inclination

_{j}*i*, the right ascension of the ascending node

_{j}*Ω*, the argument of periastron

_{j}*ω*, and the mean anomaly

_{j}*M*, the true anomaly

_{j}*f*is used, where

_{j}*j*=1, 2 for the inner orbit (

*j*=1) and outer orbit (

*j*=2).

The Hamiltonian of the system is determined as (Harrington 1968; Ford et al. 2000, 2004; Takeda et al. 2008):(1)

where *g* is the gravitational constant, *α*=*a*_{1}/*a*_{2},*p*_{1} are Legendre polynomials, is the angle between *r*_{1} and *r*_{2}, and

The Delaunay’s elements are defined as: *l _{j}* is the mean of the anomalies of object, the arguments of periastron,

*g*and the longitudes of the ascending nodes,

_{j}*h*and conjugate momenta,

_{j}

The Hamiltonian is given by Michaely & Perets (2014) as:

where

with

and cosφ can be written using orbital elements as in Rahoma (2014) and Rahoma & Deleflie (2014)(3)

where the non-vanishing coefficients are

## 3 CALCULATIONS OF THE LONG TERMS DYNAMICS

In most cases, it is better to present the movement over the long term, especially when discussing exoplanets. Therefore, the short period terms must be excluded by averaging the Hamiltonian (Eq. (2)) of the system over *l*_{1} and *l*_{2}, where the conjugate momenta, *L*_{1} and *L*_{1}, were included in the unperturbed terms to provide a smooth motion over the long period .

The standard average definition is $\tilde{Q}=\frac{1}{2\pi}{\displaystyle {\int}_{0}^{2\pi}Q\text{\hspace{0.17em}}dl}$, and the normalized Hamiltonian is given as(4)(5)(6)

where

and the non-vanishing coefficients are(7)

The non-vanishing coefficients ${\alpha}_{jsk}\left({e}_{1},\text{\hspace{0.17em}}{e}_{2}\right)$ and ${\beta}_{jsk}\left({i}_{1},\text{\hspace{0.17em}}{i}_{2}\right)$ can be found in Appendix I.

Now the equations of motion can be written by using the Hamiltonian equations. For a system (*q _{i}*,

*Q*),

_{i}*q*are the coordinates and

_{i}*Q*are the momenta with Hamiltonian

_{i}*F*(

*q*,

_{i}*Q*,

_{i}*t*).

The equations of motion can be derived as:(8)

where *t* reffered to the time

## 4 STAR MASS LOSS

The key process in star studies, generally and especially for massive stars, is the mass loss; that is, the slow change in the mass of the star with respect to its orbital period. This must be understood well along with its effect on the surrounding environment. This will then guarantee complete success in astrophysical applications. All stars are undergoing mass loss with the amount depending on the star's size and its evolution stage. For instance, the Sun is losing 7 billion tons per hour, while supergiant stars emit a larger amount of energy.

In this study, the long period dynamics with mass loss will be treated by assuming the mass loss occurs isotropically. In addition, the mass loss rate will be varied depending on stellar evolutionary stage and the mechanism of mass loss. For more information about mass loss refer to Rahoma et al. (2009) and the references therein.

To simplify the dynamics, the mass loss from mass m will be assumed to be of the form(9)

where α is a constant.

## 5 NUMERICAL SIMULATIONS

In this section, the simulated effect of mass loss on the exoplanet is described. The second companion in the binary star system is chosen to loss mass over time at different rates. The orbital element variations constructed from the Hamiltonian equations are investigated numerically.

We suppose that a hierarchical triple system has the intiial condition of *a*_{1}=3 AU, *e*_{1}=0.23, *i*_{1}=67°, *ω*_{1}=179° and *Ω*_{1}=120° and *a*_{2}=40 AU, *e*_{2}=0.5, *i*_{2}=45°, *ω*_{2}=120° and *Ω*_{2}=15°, where the mass of *m*_{0} and *m*_{1} are chosen as 15 and 5 times the mass of the Sun, respectively, and the mass of the planet, *m*_{2}, is 0.01 M_{⊙}. The integration is run over a long period of 5 × 105 years.

## 6 RESULTS, COMMENTS AND DISCUSSION

Figs. 2-6 show the variation of the orbital elements as a function of time for a planet orbiting a binary star system in which the second companion member loses mass. The initial conditions were chosen empirically. Smooth curves were presented as expected; therefore, the average technique is used to eliminate the short period terms.Fig. 3Fig. 4Fig. 5

In Figs. 2-6, the orbital element evaluation is presented where the black line represents the case with no mass loss, the red line represents the case with mass loss = 10^{-6} M_{⊙}/yr, the green line represents the case with mass loss = 2 × 10^{-6} M_{⊙}/yr, the blue line represents the case with mass loss = 3 × 10^{-6} M_{⊙}/yr, and the fuchsia line represents the case with mass loss = 4 × 10^{-6} M_{⊙}/yr. It is noted that: 1) as the value of mass loss increases, the behavior of the orbital element evolution changes and becomes distant from the behavior of the black line; 2) even if all behaviors are periodic, the amplitude changes with the value of mass loss; 3) in the first approximately 15 to 20 × 10^{4} years, all behaviors nearly coincide; but after that time, the behaviors are completely different; and 4) in all figures, the secular term is imposed on the long period terms, where there is a slope in the curves.

## 7 CONCLUSIONS

Supported by observations, many exoplanets are orbiting binary star systems, with nearly 70 exoplanets found in binary stellar systems (Schwarz et al. 2015). This work studied the effects of stellar mass loss from one star in the inner binary system on its companion planetary orbit. It was found that the mass loss causes the orbit of the planet to spiral outward.

Using triple hierarchical system dynamics, the Hamiltonian of the system was constructed up to the second order with the relative semimajor axis of the binary star to the semimajor axis of the planet as a small parameter. Previous work regarding the effect of stellar mass loss on planetary orbits showed that the semimajor axis and mass production remained nearly constant (Veras et al. 2011; Mustill & Villaver 2012) unlike the results of this study, which adopted mass changes with time. The behavior of the orbital elements over a long time with different rates of mass loss was presented. Smooth curves over a long time resulted; As the short period terms are averaged. The trajectory presented an oscillatory behavior along the time where the oscillation period amplitude decreased with an increasing value of the outer companion mass loss. The accuracy could be improved by raising the order of expansion. The mass change models could be imposed on the equations.