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 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.