 disc is larger than a critical angle $I_c$ given by:

$$\cos^2 I_c = \frac{3}{5}.$$

The maximum value reached by the eccentricity is then (Innanen et al. 1997):

$$e_{\rm max}=\left( 1- \frac{5}{3} \cos^2 I_0 \right)^{1/2},
\label{emax}$$

and the time $t_{\rm evol}$ it takes to reach $e_{\rm max}$ starting from $e_0$ is (Innanen et al. 1997):

$$\frac{t_{\rm evol}}{\tau} =  
0.42  
\left( \sin^2 I_0 - \frac{2}{5}  \right)^{-1/2} 
\ln \left( \frac{e_{\rm max}}{e_0}  \right),
\label{tevol}$$

with the time $\tau$ defined as:

$$\tau=\frac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1})}
 \frac{R_o^3 M_{\star}}{a^3 M_{\rm disc}} \frac{T}{2 \pi},
\label{tau}$$

where $T$ is the orbital period of the planet and $a$ is its semi–major axis. Note that the function $ \left( \sin^2 I_0 - 2/5
\right)^{-1/2}$ decreases very sharply from infinity to $\sim 3$ as $I_0$ increases from $I_c$ to about 45$^{\circ}$ and then decreases by about 50% as $I_0$ continues to increase up to 90$^{\circ}$.

The $Z$–component of the angular momentum of the orbit, $L_z \propto
\sqrt{1-e^2} \cos I$, where $I$ is the inclination angle between the orbit and the plane of the disc, is constant. Therefore $I$ also oscillates with time and is out of phase with $e$.

Numerical simulations {#sec:simulations}
=====================

We consider a star of mass $M_{\star}=1$ M$_{\sun}$ surrounded by a disc of mass $M_{\rm disc}$
disc is larger than a critical angle $ I_c$ given by: 

 $ $ \cos^2 I_c = \frac{3}{5}.$$ 

 The maximal value pass by the eccentricity is then (Innanen et al.   1997 ): 

 $ $ e_{\rm max}=\left (1- \frac{5}{3 } \cos^2 I_0 \right)^{1/2 }, 
 \label{emax}$$ 

 and the time $ t_{\rm evol}$ it lead to pass $ e_{\rm max}$ starting from $ e_0 $ is (Innanen et al.   1997 ): 

 $ $ \frac{t_{\rm evol}}{\tau } =  
 0.42  
 \left (\sin^2 I_0 - \frac{2}{5 }   \right)^{-1/2 } 
 \ln \left (\frac{e_{\rm max}}{e_0 }   \right), 
 \label{tevol}$$ 

 with the time $ \tau$ define as: 

 $ $ \tau=\frac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1 }) } 
  \frac{R_o^3 M_{\star}}{a^3 M_{\rm disk } } \frac{T}{2 \pi }, 
 \label{tau}$$ 

 where $ T$ is the orbital period of the satellite and $ a$ is its semi – major bloc. Note that the function $ \left (\sin^2 I_0 - 2/5 
 \right)^{-1/2}$ decrease very sharply from infinity to $ \sim 3 $ as $ I_0 $ increases from $ I_c$ to approximately 45$^{\circ}$ and then decreases by about 50% as $ I_0 $ continue to increase up to 90$^{\circ}$. 

 The $ Z$–component of the angular momentum of the orbit, $ L_z \propto 
 \sqrt{1 - e^2 } \cos I$, where $ I$ is the inclination slant between the orbit and the plane of the disc, is ceaseless. Therefore $ I$ also hover with time and is out of phase with $ e$. 

 Numerical simulation { # sec: simulations } 
 = = = = = = = = = = = = = = = = = = = = = 

 We consider a star of mass $ M_{\star}=1 $   M$_{\sun}$ surrounded by a magnetic disk of mass $ M_{\rm disc}$
 didc is larger than a critlcal angle $I_c$ girwn by:

$$\cms^2 I_c = \frac{3}{5}.$$

The maximum value reached by thx ecxentrucity is then (Innanen dt al. 1997):

$$e_{\rm lax}=\left( 1- \frar{5}{3} \cos^2 I_0 \right)^{1/2},
\labxm{emax}$$

and the flme $t_{\xm evol}$ it takes to reach $a_{\rm max}$ startitg fxom $e_0$ is (Innanen et al. 1997):

$$\frac{t_{\rm evol}}{\tat} =  
0.42  
\lrfh( \sin^2 I_0 - \frac{2}{5}  \rigne)^{-1/2} 
\ln \left( \frac{e_{\rm max}}{e_0}  \right),
\label{tevom}$$

with tie time $\tau$ defoned as:

$$\tau=\frac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1})}
 \fgac{R_l^3 M_{\star}}{a^3 M_{\rm disc}} \vrac{T}{2 \pi},
\labgm{tat}$$

qhere $T$ is tfe orbital period of tge planet and $a$ is its semi–majof axix. Note thaj fhf function $ \oeft( \fin^2 I_0 - 2/5
\right)^{-1/2}$ decreasas very sharply from lnfinmty ro $\sim 3$ as $I_0$ increasev from $I_c$ to about 45$^{\circ}$ and tken decreases by abour 50% as $I_0$ condinudw tu ihcceaae up ho 90$^{\rirc}$.

The $Z$–cojponent of rhe angular momentuk jd the orbit, $L_a \propeo
\fqrt{1-e^2} \cos I$, where $I$ is the inclination dngme between the orbit ane the plane of the didc, is conftant. Therefore $I$ also oscillates with time and iv out uf khqse wkrh $e$.

Numerical simulations {#sec:simulations}
=====================

We consyset s star of mass $M_{\star}=1$ M$_{\sun}$ xugrponded by a disz of mcas $M_{\rm disc}$
disc is larger than a critical angle by: I_c = The maximum value then et al. 1997): max}=\left( 1- \frac{5}{3} I_0 \right)^{1/2}, \label{emax}$$ and the time evol}$ it takes to reach $e_{\rm max}$ starting from $e_0$ is (Innanen et 1997): $$\frac{t_{\rm evol}}{\tau} = 0.42 \left( \sin^2 I_0 - \frac{2}{5} \right)^{-1/2} \ln \left( max}}{e_0} \label{tevol}$$ the $\tau$ defined as: $$\tau=\frac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1})} \frac{R_o^3 M_{\star}}{a^3 M_{\rm disc}} \frac{T}{2 \pi}, \label{tau}$$ where $T$ is the orbital of the planet and $a$ is its semi–major Note that the function \left( \sin^2 I_0 - 2/5 decreases sharply from to 3$ $I_0$ increases from to about 45$^{\circ}$ and then decreases by about 50% as $I_0$ continues to increase up to 90$^{\circ}$. $Z$–component of momentum of orbit, \propto \cos I$, where the inclination angle between the orbit of the disc, is constant. Therefore $I$ also with time is out of phase with $e$. simulations {#sec:simulations} ===================== We consider a star of $M_{\star}=1$ M$_{\sun}$ surrounded by a disc of mass $M_{\rm disc}$
 disc is larger than a critical Angle $I_c$ givEn by:

$$\cOs^2 I_C = \frAc{3}{5}.$$

the mAximUm value reached BY the Eccentricity is then (InnaNen et Al. 1997):

$$E_{\Rm maX}=\LeFt( 1- \fraC{5}{3} \cos^2 I_0 \riGHt)^{1/2},
\LABel{EmAx}$$

And ThE TiMe $t_{\rm EvoL}$ it takeS to reach $e_{\rM maX}$ sTarting from $e_0$ IS (INnanen et al. 1997):

$$\FraC{t_{\rm evol}}{\tau} =  
0.42  
\lEft( \Sin^2 I_0 - \frAc{2}{5}  \RigHT)^{-1/2} 
\ln \leFt( \fRac{e_{\rM max}}{e_0}  \rIGht),
\labEl{tevol}$$

wiTh THe time $\TAu$ definED As:

$$\Tau=\fRac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1})}
 \frac{r_O^3 M_{\STar}}{a^3 M_{\rm disc}} \fraC{T}{2 \pi},
\laBeL{TaU}$$

WHerE $T$ iS the orbitaL pEriod OF the plaNEt AND $A$ is ITs semi–major axIs. Note that tHE fuNction $ \LeFt( \sIN^2 I_0 - 2/5
\righT)^{-1/2}$ decrEaSEs vEry sharply fRom iNfinity to $\Sim 3$ as $I_0$ INcreaseS From $I_c$ tO about 45$^{\CirC}$ anD theN DeCrEasEs BY abOUt 50% As $I_0$ COntInues to iNcReAse up To 90$^{\ciRC}$.

tHE $Z$–coMpoNent Of the Angular momentUm oF the ORbiT, $L_z \prOpto
\sQrt{1-e^2} \CoS I$, wheRe $I$ is tHe incLiNation angle betwEen tHe orbit anD thE pLanE oF the dISc, is coNstAnt. therefoRe $I$ also OSciLlATES wIth time and is out of pHaSE WiTh $e$.

NumerIcal siMUlAtIOns {#sec:siMuLatIons}
=====================

wE ConsiDer a STaR of mass $M_{\Star}=1$ M$_{\sUN}$ sUrRounded By A disc oF mAss $m_{\rm Disc}$
 disc is larger than a cri tical angl e $I_ c$giv en  by: 

$$ \cos^2 I_c = \ f rac{ 3}{5}.$$

The maximumvalue  r e ache d  b y the  eccent r ic i t y i sth en(I n na nen e t a l. 1997 ):

$$e_{\ rmma x}=\left( 1- \f rac{5}{3}\co s^2 I_0 \rig ht) ^{1/2} ,\la b el{em ax} $$

a nd the time $ t_{\rm ev ol } $ it t a kes tor e ac h $e _{\rm max}$ start i ng from $e_0$ is(Innan en et a l.199 7):

$$\fr ac {t_{\ r m evol} } {\ t a u } =  
0.42  
\lef t( \sin^2 I _ 0 -  \frac {2 }{5 }   \rig ht)^{ -1 / 2}
\ln \left(  \fr ac{e_{\rm  max}} { e_0}  \ r ight),\label {te vol }$$wi th  th et ime $\ tau $  de fined as :
$ $\tau =\fr a c { ( 1+n) (1- \eta ^{-n+ 2})}{(-n+2)(- 1+\ eta^ { -n- 1})} \fra c{R_ o^ 3 M_{ \star} }{a^3  M _{\rm disc}} \f rac{ T}{2 \pi} ,
\ la bel {t au}$$ 
where  $T $ i s the o rbitalp eri od o f  t he planet and $a$is i ts  semi–ma jor ax i s.  N o te thatth e f unct i o n $ \ left (  \ sin^2 I_ 0 - 2/ 5 
\ ri ght)^{- 1/ 2}$ de cr eas esverys harp ly fro m infini ty to $\sim 3$ as $I _ 0$ increasesf ro m $I _ c$ t o a bout 45$^{\ circ } $ an d th e ndec r eases  by a bo u t5 0% as $I_0$ continu es  to in creas e up to 90$^{ \circ}$.
T h e  $Z$–com pone n to f the angularmomen tum of the orbit, $ L_z \ propto
\ sqrt{1-e^ 2 }  \cos I$ , w her e $ I$i s  t he inclinatio n angl ebetween  th e orbit  an d t hepla ne  of the d isc, isco ns ta nt . T heref o re $I$ a ls o o sc ill atesw ith ti me an d is  o ut ofphase w i th $ e$.
N um eric alsi mulat ions {#s ec:simu lations}=== = ==== == == ======= ==

We consid er  a star of  m ass  $M_{\ s t ar}=1$ M $_{\sun}$ surrounded by a discofmass$M_{ \rm disc} $
 disc_is larger_than a critical angle_$I_c$ given_by:

$$\cos^2_I_c =_\frac{3}{5}.$$

The_maximum value reached_by the eccentricity_is then (Innanen et_al. 1997):

$$e_{\rm max}=\left( 1-_\frac{5}{3}_\cos^2 I_0 \right)^{1/2},
\label{emax}$$

and the time $t_{\rm evol}$ it takes to reach $e_{\rm max}$ starting_from_$e_0$ is_(Innanen_et_al. 1997):

$$\frac{t_{\rm evol}}{\tau} =  
0.42_ 
\left( \sin^2 I_0 -_\frac{2}{5} _\right)^{-1/2} 
\ln \left( \frac{e_{\rm max}}{e_0}  \right),
\label{tevol}$$

with the_time_$\tau$ defined as:

$$\tau=\frac{(1+n)(1-\eta^{-n+2})}{(-n+2)(-1+\eta^{-n-1})}
_\frac{R_o^3 M_{\star}}{a^3 M_{\rm disc}} \frac{T}{2 \pi},
\label{tau}$$

where $T$ is the_orbital period of the planet and_$a$ is its_semi–major_axis._Note that the function_$ \left( \sin^2 I_0 - 2/5
\right)^{-1/2}$_decreases very sharply from infinity to_$\sim 3$ as $I_0$ increases from $I_c$_to about 45$^{\circ}$ and then decreases_by about 50% as $I_0$_continues to_increase up to 90$^{\circ}$.

The $Z$–component_of the angular_momentum of_the orbit, $L_z_\propto
\sqrt{1-e^2} \cos I$, where $I$ is_the inclination angle_between the orbit and the plane_of_the disc, is_constant._Therefore_$I$ also_oscillates with time_and_is out_of_phase with $e$.

Numerical simulations {#sec:simulations}
=====================

We consider_a_star of mass $M_{\star}=1$ M$_{\sun}$ surrounded by a_disc of mass $M_{\rm_disc}$