We perform numerical simulations to investigate tidal evolution of two single-planet systems, that is, WASP-50 and GJ 1214 and a two-planet system CoRoT-7. The results of orbital evolution show that tidal decay and circularization may play a significant role in shaping their final orbits, which is related to the initial orbital data in the simulations. For GJ 1214 system, different cases of initial eccentricity are also considered as only an upper limit of its eccentricity (0.27) is shown, and the outcome suggests a possible maximum initial eccentricity (0.4) in the adopted dynamical model. Moreover, additional runs with alternative values of dissipation factor $Q^\prime_1$ are carried out to explore tidal evolution for GJ 1214b, and these results further indicate that the real $Q^\prime_1$ of GJ 1214b may be much larger than its typical value, which may reasonably suggest that GJ 1214b bears a present-day larger eccentricity, undergoing tidal circularization at a slow rate. For the CoRoT-7 system, tidal forces make two planets migrating towards their host star as well as producing tidal circularization, and in this process tidal effects and mutual gravitational interactions are coupled with each other. Various scenarios of the initial eccentricity of the outer planet have also been done to investigate final planetary configuration. Tidal decay arising from stellar tides may still work for each system as the eccentricity decreases to zero, and this is in association with the remaining lifetime of each planet used to predict its future.

We investigate the formation of terrestrial planets in the late stage of planetary formation using two-planet model. At that time, the protostar has formed for about 3 Myr and the gas disk has dissipated. In the model, the perturbations from Jupiter and Saturn are considered. We also consider variations of the mass of outer planet, and the initial eccentricities and inclinations of embryos and planetesimals. Our results show that, terrestrial planets are formed in 50 Myr, and the accretion rate is about 60% - 80%. In each simulation, 3 - 4 terrestrial planets are formed inside "Jupiter" with masses of $0.15 - 3.6 M_{\oplus}$. In the 0.5 - 4AU, when the eccentricities of planetesimals are excited, planetesimals are able to accrete material from wide radial direction. The plenty of water material of the terrestrial planet in the Habitable Zone may be transferred from the farther places by this mechanism. Accretion may also happen a few times between two giant planets only if the outer planet has a moderate mass and the small terrestrial planet could survive at some resonances over time scale of $10^8$ yr.

A diversity of resonance configurations may be formed under different migration of two giant planets. And the researchers show that the HD 128311 and HD 73526 planetary systems are involved in a 2:1 mean motion resonance but not in apsidal corotation, because one of the resonance argument circulates over the dynamical evolution. In this paper, we investigate potential mechanisms to form the 2:1 librating-circulating resonance configuration. In the late stage of planetary formation, scattering or colliding among planetesimals and planetary embryos can frequently occur. Hence, in our model, we consider a planetary configuration of two giants together with few terrestrial planets. We find that both colliding or scattering events at very early stage of dynamical evolution can influence the configurations trapped into resonance. A planet-planet scattering of a moderate terrestrial planet, or multiple scattering of smaller planets in a crowded planetary system can change the resonant configuration. In addition, collision or merging can alter the masses and location of the giant planets, which also play an important role in shaping the resonant configuration during the dynamical evolution. In this sense, the librating-circulating resonance configuration is more likely to form by a hybrid mechanism of scattering and collision.

We investigate the formation of terrestrial planets in the late stage of planetary formation using two-planet model. At that time, the protostar has formed for about 3 Myr and the gas disk has dissipated. In the model, the perturbations from Jupiter and Saturn are considered. We also consider variations of the mass of outer planet, and the initial eccentricities and inclinations of embryos and planetesimals. Our results show that, terrestrial planets are formed in 50 Myr, and the accretion rate is about $60% - 80%$. In each simulation, 3 - 4 terrestrial planets are formed inside "Jupiter" with masses of $0.15 - 3.6 M_{\oplus}$. In the $0.5 - 4$ AU, when the eccentricities of planetesimals are excited, planetesimals are able to accrete material from wide radial direction. The plenty of water material of the terrestrial planet in the Habitable Zone may be transferred from the farther places by this mechanism. Accretion could also happen a few times between two major planets only if the outer planet has a moderate mass and the small terrestrial planet could survive at some resonances over time scale of $10^8$ yr. In one of our simulations, com-mensurability of the orbital periods of planets is very common. Moreover, a librating-circulating 3:2 configuration of mean motion resonance is found.

We perform numerical simulations to study the Habitable zones (HZs) and dynamical structure for Earth-mass planets in multiple planetary systems. For example, in the HD 69830 system, we extensively explore the planetary configuration of three Neptune-mass companions with one massive terrestrial planet residing in 0.07 AU $\leq a \leq$ 1.20 AU, to examine the asteroid structure in this system. We underline that there are stable zones of at least $10^5$ yr for low-mass terrestrial planets locating between 0.3 and 0.5 AU, and 0.8 and 1.2 AU with final eccentricities of $e < 0.20$. Moreover, we also find that the accumulation or depletion of the asteroid belt are also shaped by orbital resonances of the outer planets, for example, the asteroidal gaps at 2:1 and 3:2 mean motion resonances (MMRs) with Planet C, and 5:2 and 1:2 MMRs with Planet D. In a dynamical sense, the proper candidate regions for the existence of the potential terrestrial planets or HZs are 0.35 AU $< a < $ 0.50 AU, and 0.80 AU $< a < $ 1.00 AU for relatively low eccentricities, which makes sense to have the possible asteroidal structure in this system.

We perform numerical simulations to study the secular orbital evolution and dynamical structure in the quintuplet planetary system 55 Cancri with the self-consistent orbital solutions by Fischer and coworkers (2008). In the simulations, we show that this system can be stable at least for $10^{8}$ yr. In addition, we extensively investigate the planetary configuration of four outer companions with one terrestrial planet in the wide region of 0.790 AU $\leq a \leq $ 5.900 AU to examine the existence of potential asteroid structure and Habitable Zones (HZs). We show that there are unstable regions for the orbits about 4:1, 3:1 and 5:2 mean motion resonances (MMRs) with the outermost planet in the system, and several stable orbits can remain at 3:2 and 1:1 MMRs, which is resemblance to the asteroidal belt in solar system. In a dynamical point, the proper candidate HZs for the existence of more potential terrestrial planets reside in the wide area between 1.0 AU and 2.3 AU for relatively low eccentricities.

We use a multi-dimensional hydrodynamics code to study the gravitational interaction between an embedded planet and a protoplanetary disk with emphasis on the generation of vortensity (Potential Vorticity or PV) through a Baroclinic Instability. We show that the generation of PV is very common and effective in non-barotropic disks through the Baroclinic Instability, especially within the coorbital region. Our results also complement previous work that non-axisymmetric Rossby-Wave Instabilities (RWIs, Lovelace et al. 1999) are likely to develop at local minima of PV distribution that are generated by the interaction between a planet and an inviscid barotropic disk. The development of RWIs results in non-axisymmetric density blobs, which exert stronger torques onto the planet when they move to the vicinity of the planet. Hence, large amplitude oscillations are introduced to the time behavior of the total torque acted on the planet by the disk. In current simulations, RWIs do not change the overall picture of inward orbital migration but cause a non-monotonic behavior to the migration speed. As a side effect, RWIs also introduce interesting structures into the disk. These structures may help the formation of Earth-like planets in the Habitable Zone or Hot Earths interior to a close-in giant planet.

We investigate the secular resonances for massless small bodies and Earth-like planets in several planetary systems. We further compare the results with those of Solar System. For example, in the GJ 876 planetary system, we show that the secular resonances $\nu_1$ and $\nu_2$ (respectively, resulting from the inner and outer giant planets) can excite the eccentricities of the Earth-like planets with orbits 0.21 AU $\leq a <$ 0.50 AU and eject them out of the system in a short timescale. However, in a dynamical sense, the potential zones for the existence of Earth-like planets are in the area 0.50 AU $\leq a \leq$ 1.00 AU, and there exist all stable orbits last up to $10^5$ yr with low eccentricities. For other systems, e.g., 47 UMa, we also show that the Habitable Zones for Earth-like planets are related to both secular resonances and mean motion resonances in the systems.

We perform numerical simulations to study the secular orbital evolution and dynamical structure in the HD 69830 system with the best-fit orbital solutions by Lovis and coworkers (2006). In the simulations, we show that the triplet Neptunian system can be stable at least for 2 Gyr and the stability would not be greatly influenced even if we vary the planetary masses. In addition, we employ the Laplace-Lagrange secular theory to investigate the long-term behaviors of the system, and the outcomes demonstrate that this theory can well describe the secular orbital evolution for all planets, where the secular periods and amplitudes in the eccentricities well agrees with those of the direct numerical integrations. We first reveal that the secular periods of the eccentricity $e_{1}$ and $e_{2}$ are identical about 8,300 yr. Moreover, we explore the planetary configuration of three Neptune-mass companions with one massive terrestrial planet in 0.07 AU $\leq a \leq 1.20$ AU, to examine the asteroid structure in this system. We underline that there are stable zones at least $10^{5}$ yr for low-mass terrestrial planets locating between 0.3 and 0.5 AU, and 0.8 and 1.2 AU with final low eccentricities. Still, we also find that the secular resonance $\nu_{1}$ and $\nu_{2}$ of two inner planets can excite the eccentricities of the terrestrial bodies, and the accumulation or depletion of the asteroid belt are also shaped by orbital resonances of the outer planets, i.e., 5:2 and 1:2 MMRs with Planet D... (abridged)

We have numerically explored the stable planetary geometry for the multiple systems involved in a 2:1 mean motion resonance, and herein we mainly study the HD 82943 system by employing two sets of the orbital parameters (Mayor et al. 2004; Ji et al. 2004). In the simulations, we find that all stable orbits are related to the 2:1 resonance that can help to remain the semi-major axes for two companions almost unaltered over the secular evolution for $10^{8}$ yr. In addition, we also show that there exist three possible stable configurations:(1) Type I, only $\theta_{1} \approx 0^{\circ}$, (2) Type II, $\theta_{1}\approx\theta_{2}\approx\theta_{3}\approx 0^{\circ}$ (aligned case), and (3) Type III, $\theta_{1}\approx 180^{\circ}$, $\theta_{2}\approx0^{\circ}$, $\theta_{3}\approx180^{\circ}$ (antialigned case), where two resonant arguments are $\theta_{1} = \lambda_{1} - 2\lambda _{2} + \varpi_{1}$ and $\theta_{2} = \lambda_{1} - 2\lambda_{2} + \varpi_{2}$, the relative apsidal longitudes $\theta_{3} = \varpi_{1}-\varpi_{2}=\Delta\varpi$. And we find that other 2:1 resonant systems (e.g., GJ 876) may possess one of three stable orbits in their realistic motions. Moreover, we also study the existence of the assumed terrestrial bodies at $\sim 1$ AU for HD 82943 and GJ 876 systems (see main texts).

We integrate the orbital solutions of the planets orbiting 55 Cnc. In the simulations, we find that not only three resonant arguments $\theta_{1}=\lambda_{1}-3\lambda_{2}+2\tilde\omega_{1}$, $\theta_{2}=\lambda_{1}-3\lambda_{2}+2\tilde\omega_{2}$ and $\theta_{3}=\lambda_{1}-3\lambda_{2}+(\tilde\omega_{1}+\tilde\omega_{2})$ librate respectively, but the relative apsidal longitudes $\Delta\omega$ also librates about $250^{\circ}$ for millions of years. The results imply the existence of the 3:1 resonance and the apsidal resonance for the studied system. We emphasize that the mean motion resonance and apsidal locking can act as two important mechanisms of stabilizing the system. In addition, we further investigate the secular dynamics of this system by comparing the numerical results with those given by Laplace-Lagrange secular theory.

(Abridged) We have numerically explored the stable planetary geometry for the multiple systems involved in a 2:1 mean motion resonance, and herein we mainly concentrate on the study of the HD 82943 system by employing two sets of the orbital parameters (Mayor et al. 2004). We find that all stable orbits are related to the 2:1 commensurability for $10^{7}$ yr, and the apsidal phase-locking between two orbits can further enhance the stability for this system. For HD 82943, there exist three possible stable configurations:(1) Type I, only $\theta_{1} \approx 0^{\circ}$,(2) Type II, $\theta_{1}\approx\theta_{2}\approx\theta_{3}\approx 0^{\circ}$ (aligned case), and (3) Type III, $\theta_{1}\approx 180^{\circ}$, $\theta_{2}\approx0^{\circ}$, $\theta_{3}\approx180^{\circ}$ (antialigned case), here the lowest eccentricity-type mean motion resonant arguments are $\theta_{1} = \lambda_{1} - 2\lambda_{2} + \varpi_{1}$ and $\theta_{2} = \lambda_{1} - 2\lambda_{2} + \varpi_{2}$, the relative apsidal longitudes $\theta_{3} = \varpi_{1}-\varpi_{2}=\Delta\varpi$. In addition, we also propose a semi-analytical model to study $e_{i}-\Delta\varpi$ Hamiltonian contours. With the updated fit, we examine the dependence of the stability of this system on the orbital parameters. Moreover, we numerically show that the assumed terrestrial bodies cannot survive near the habitable zones for HD 82943 and low-mass planets can be dynamically habitable in the GJ 876 system at $\sim 1$ AU in the numerical surveys.

We perform numerical simulations to explore the dynamical evolution of the HD 82943 planetary system. By simulating diverse planetary configurations, we find two mechanisms of stabilizing the system: the 2:1 mean motion resonance between the two planets can act as the first mechanism for all stable orbits. The second mechanism is a dynamical antialignment of the apsidal lines of the orbiting planets, which implies that the difference of the periastron longitudes $\theta_{3}$ librates about $180^{\circ}$ in the simulations. We also use a semi-analytical model to explain the numerical results for the system under study.

We carry out numerical simulations to explore the dynamical evolution of the HD 82943 and HD 37124 planetary systems,which both have two Jupiter-like planets. By simulating various planetary configurations in the neighborhood of the fitting orbits, we find three mechanisms to maintain the stability of these systems: For HD 82943,we find that the 2:1 mean motion resonance can act as the first mechanism for all the stable orbits. The second mechanism is the alignment of the periastron of the two planets of HD 82943 system. In the paper,we show one case is simultaneously maintained by the two mechanisms. Additionally,we also use the corresponding analytical models successfully to explain the different numerical results for the system. The third mechanism is the Kozai resonance which takes place in the mutual highly orbits of HD 37124. In the simulations,we discover that the argument of periastron $\omega$ of the inner planet librates about $90^{\circ}$ or $270^{\circ}$ for the whole time span. The Kozai mechanism can explain the stable configuration of large eccentricity of the inner planet.

We have performed simulations to investigate the dynamics of the M dwarf star GJ 876 in an attempt to reveal any stabilizing mechanism for sustaining the system.We simulated different coplanar and noncoplanar configurations of two-planet systems and other cases.From the simulations,we found that the 2 :1 mean-motion resonance between two planets can act as an effective mechanism for maintaining the stability of the system.This result is explained by a proposed analytical model.Using this model,we studied the region of motion of the inner planet by varying the parameters of the system,and we detected that the analytical results are well consistent with the numerical simulations.