We report the results of the transit timing variation (TTV) analysis of the extra-solar planet Qatar-1b using thirty eight light curves. Our analysis combines thirty five previously available transit light curves with three new transits observed by us between June 2016 and September 2016 using the 2-m Himalayan Chandra Telescope (HCT) at the Indian Astronomical Observatory (Hanle, India). From these transit data, the physical and orbital parameters of the Qatar-1 system are determined. In addition to this, the ephemeris for the orbital period and mid-transit time are refined to investigate the possible TTV. We find that the null-TTV model provides the better fit to the (O-C) data. This indicates that there is no evidence for TTVs to confirm the presence of additional planets in the Qatar-1 system. The use of the 3.6-m Devasthal Optical Telescope (DOT) operated by the Aryabhatta Research Institute of Observational Sciences (ARIES, Nainital, India) could improve the photometric precision to examine the signature of TTVs in this system with a greater accuracy than in the present work.

We present a comprehensive analysis (photometric and kinematical) of poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia~DR2 database. By determining the membership probabilities of stars, we identified 624 most probable members with membership probability higher than $50\%$ by using proper motion and parallax data taken from Gaia~DR2. The mean proper motion of the cluster is obtained as $\mu_{x}=-8.83\pm0.01$ and $\mu_{y}=1.49\pm0.006$ mas yr$^{-1}$. We find the normal interstellar extinction towards the cluster region. The radial distribution of members provides a cluster radius of 7.75 arcmin (5.63 pc). The estimated age of $1600\pm180$ Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as $1.46\pm0.18$ within the mass range 0.75$-$2.0 $M_\odot$, which is in fair agreement with Salpeter's value (x=1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates $(A, D)$ by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are: $(A_{conv}, D_{conv})$ = (96$^{\textrm{o}}$.27 $\pm$ 0$^{\textrm{o}}$.10, 13$^{\textrm{o}}$.14 $\pm$ 0$^{\textrm{o}}$.27) $\&$ $(A_\circ, D_\circ)$ = (100$^{\textrm{o}}$.282 $\pm$ 0$^{\textrm{o}}$.10, 9$^{\textrm{o}}$.577 $\pm$ 0$^{\textrm{o}}$.323) respectively. We also computed the Velocity Ellipsoid Parameters (VEPs), matrix elements ($\mu_{ij}$), direction cosines ($l_j$, $m_j$, $n_j$) and the Galactic longitude of the vertex ($l_2$).

We present twelve new transit light curves of the hot-Jupiter TrES-3b observed during $2012-2018$ to probe the transit timing variation (TTV). By combining the mid-transit times determined from these twelve transit data with those re-estimated through uniform procedure from seventy one transit data available in the literature, we derive new linear ephemeris and obtain the timing residuals that suggest the possibility of TTV in TrES-3 system. However, the frequency analysis shows that the possible TTV is unlikely to be periodic, indicating the absence of an additional body in this system. To explore the other possible origins of TTV, the orbital decay and apsidal precession ephemeris models are fitted to the transit time data. We find decay rate of TrES-3b to be $\bf \dot{P_q}= -4.1 \pm 3.1$ $ms$ ${yr}^{-1}$ and the corresponding estimated modified tidal quality factor of ${Q}^{'}_{\ast}$ $\sim 1.11 \times {10}^{5}$ is consistent with the theoretically predicted values for the stars hosting the hot-Jupiters. The shift in the transit arrival time of TrES-3b after eleven years is expected to be ${T}_{shift}\sim 69.55 \ s$, which is consistent with the $RMS$ of the timing residuals. Besides, we find that the apsidal precession ephemeris model is statistically less probable than the other considered ephemeris models. It is also discussed that despite the linear ephemeris model appears to be the most plausible model to represent the transit time data, the possibility of the orbital decay cannot be completely ruled out in TrES-3 system. In order to confirm this, further high-precision and high-cadence follow-up observation of transits of TrES-3b would be important.

The kinematical parameters, spatial shape and structure of the open cluster IC 2391 and the associated stellar stream are studied here using Gaia-DR2 (GDR2) astrometry data. The apex positions are determined for the open cluster IC 2391 (data taken from Cantat-Gaudin et al.) and for the kinematical streams stars mentioned in Montes et al. using both convergent point and AD-diagram methods. The values of apex coordinates identified. The results are in good agreement with the previously calculated values. The positions of the stars in the disk and the spatial dispersion velocities are determined. The paths of cluster and associated stream are traced in the disk by orbit calculation back in time to their places of formation. A possible genetic relationship between the cluster and the stream has been detected. The approximation of the spatial and kinematical shape of the stream and the cluster is made. According to this study, even though currently the cluster and the stream seem to have spatial difference in their locations but they appear to have formed in the same region of the Galactic disk.

The resonant perturbations from planets are able to halt the drag-induced migration, and capture the inwardly drifting planetesimals into mean motion resonances. The equilibrium eccentricity of planetesimals in resonances, and the minimum size of planetesimal that can trigger resonance trapping, have been analyzed and formulated. However, the analytical works based on the assumption that the disk is axisymmetric, which is violated by the asymmetric structures developed by planets. We perform long-term 2D hydrodynamic simulations to study the dynamics of planetesimals in the $j:(j+1)$ first-order exterior resonances, and reexamine the theoretical expressions. We find the expression of equilibrium eccentricity underestimates the values for resonances with $j < 5$, in particular the 1:2 resonance that the underestimation can be $30 - 40\%$. Within the parameter space we explored, we find the equilibrium eccentricity and the minimum size are reduced in an asymmetric disk. The amount of discrepancy in eccentricity depends on the degree of asymmetric structures. For cases of Earth-sized planets, where the disk is less disturbed, the planetesimal's eccentricity can reach to the values predicted by our modified expression. For gaseous planets, however, the eccentricity can be $0.01 - 0.02$ smaller in value. We find the minimum size is 10 times smaller, and the factor seems to be independent of the planet's mass. The influences of asymmetric profiles on the eccentricity and the minimum size could affect the outcome of collisions between resonant and nonresonant planetesimals, and the amount of planetesimals migrated into the planet's feeding zone.

A machine learning technique with two-dimension convolutional neural network is proposed for detecting exoplanet transits. To test this new method, five different types of deep learning models with or without folding are constructed and studied. The light curves of the Kepler Data Release 25 are employed as the input of these models. The accuracy, reliability, and completeness are determined and their performances are compared. These results indicate that a combination of two-dimension convolutional neural network with folding would be an excellent choice for the future transit analysis.

An investigation on the possible dynamical models of the core galaxy NGC 1399 is performed. Because early-type galaxies are likely to be formed through merging events, remnant rings are considered in the modeling process. A numerical survey over three parameters is employed to obtain the best-fit models that are completely consistent with observations. It is found that the inner slope of dark matter profile is a cuspy one for this core galaxy. The existence of remnant rings in best-fit models indicates a merging history. The remnant ring explains the flatten surface brightness, and thus could be the physical counterpart of the core structure of NGC 1399.

Using astrometric techniques developed by Anderson et al., we determine proper motions (PMs) in 14.60 arcmin X 16.53 arcmin area of the kinematically "thick-disk" globular cluster M12. The cluster's proximity and sparse nature makes it a suitable target for ground-based telescopes. Archive images with time gap of 11.1 years were observed with wide-field imager (WFI) mosaic camera mounted on ESO 2.2 m telescope. The median value of PM error in both components is 0.7 mas/yr for the stars having V less than or equal to 20 mag. PMs are used to determine membership probabilities and to separate field stars from the cluster sample. In electronic form, a membership catalog of 3725 stars with precise coordinates, PMs, BVRI photometry is being provided. One of the possible applications of the catalog was shown by gathering the membership information of the variable stars, blue stragglers and X-ray sources reported earlier in the cluster's region.

A study of cluster characteristics and internal kinematical structure of the middle-aged Pleiades open star cluster is presented. The individual star apexes and various cluster kinematical parameters including the velocity ellipsoid parameters are determined using both Hipparcos and Gaia data. Modern astrometric parameters were taken from the Gaia Data Release 1 (DR1) in combination with the Radial Velocity Experiment Fifth Data Release (DR5). The necessary set of parameters including parallaxes, proper motions and radial velocities are used for n=17 stars from Gaia DR1+RAVE DR5 and for n=19 stars from the Hipparcos catalog using SIMBAD data base. Single stars are used to improve accuracy by eliminating orbital movements. RAVE DR5 measurements were taken only for the stars with the radial velocity errors not exceeding $2$~km/s. For the Pleiades stars taken from Gaia, we found mean heliocentric distance as $136.8 \pm 6.4$~pc, and the apex position is calculated as: $A_{CP}=92^\circ.52\pm1^\circ.72$, $D_{CP}=-42^\circ.28\pm2^\circ.56$ by the convergent point method and $A_0=95^\circ.59\pm2^\circ.30$ and $D_0=-50^\circ.90\pm2^\circ.04$ using AD-diagram method (n=17 in both cases). The results are compared with those obtained historically before the Gaia mission era.

With a high value of heliocentric radial velocity, a retrograde orbit, and being suspected to have an extragalactic origin, NGC 3201 is an interesting globular cluster for kinematical studies. Our purpose is to calculate the relative proper motions (PMs) and membership probability for the stars in the wide region of globular cluster NGC 3201. Proper motion based membership probabilities are used to isolate the cluster sample from the field stars. The membership catalogue will help address the question of chemical inhomogeneity in the cluster. Archive CCD data taken with a wide-field imager (WFI) mounted on the ESO 2.2m telescope are reduced using the high-precision astrometric software developed by Anderson et al. for the WFI images. The epoch gap between the two observational runs is $\sim$14.3 years. To standardize the $BVI$ photometry, Stetson's secondary standard stars are used. The CCD data with an epoch gap of $\sim$14.3 years enables us to decontaminate the cluster stars from field stars efficiently. The median precision of PMs is better than $\sim$0.8 mas~yr$^{-1}$ for stars having $V<$18 mag that increases up to $\sim$1.5 mas~yr$^{-1}$ for stars with $18<V<20$ mag. Kinematic membership probabilities are calculated using proper motions for stars brighter than $V\sim$20 mag. An electronic catalogue of positions, relative PMs, $BVI$ magnitudes and membership probabilities in $\sim$19.7$\times$17 arcmin$^2$ region of NGC 3201 is presented. We use our membership catalogue to identify probable cluster members among the known variables and $X$-ray sources in the direction of NGC 3201.

Three-dimensional equi-potential surfaces of a galactic system with supermassive binary black holes are discussed herein. The conditions of topological transitions for the important surfaces, i.e. Roche Lobes and Jiang-Yeh Lobe, are studied in this paper. In addition, the mathematical properties of the Jacobi surfaces are investigated analytically. Finally, a numerical procedure for determining the regions of the Roche Lobes and Jiang-Yeh Lobe is suggested.

In this paper, we present our results for the photometric and kinematical studies of old open cluster NGC 188. We determined various astrophysical parameters like limited radius, core and tidal radii, distance, luminosity and mass functions, total mass, relaxation time etc. for the cluster using 2MASS catalog. We obtained the cluster's distance from the Sun as 1721+/-41 pc and log (age)= 9.85+/-0.05 at Solar metallicity. The relaxation time of the cluster is smaller than the estimated cluster age which suggests that the cluster is dynamically relaxed. Our results agree with the values mentioned in the literature. We also determined the cluster's apex coordinates as (281.88 deg, -44.76 deg) using AD-diagram method. Other kinematical parameters like space velocity components, cluster center and elements of Solar motion etc. have also been computed.

Motivated by the previously reported high orbital decay rate of the planet WASP-43b, eight newly transit light curves are obtained and presented. Together with other data in literature, we perform a self-consistent timing analysis with data covering a timescale of 1849 epochs. The results give an orbital decay rate dP/dt = -0.02890795\pm 0.00772547 sec/year, which is one order smaller than previous values. This slow decay rate corresponds to a normally assumed theoretical value of stellar tidal dissipation factor. In addition, through the frequency analysis, the transit timing variations presented here are unlikely to be periodic, but could be signals of a slow orbital decay.

N-body simulations of galactic collisions are employed to investigate the formation of elliptical rings in disk galaxies. The relative inclination between disk and dwarf galaxies is studied with a fine step of five degrees. It is confirmed that the eccentricity of elliptical ring is linearly proportional to the inclination angle. Deriving from the simulational results, an analytic formula which expresses the eccentricity as a function of time and inclination angle is obtained. This formula shall be useful for the interpretations of the observations of ring systems, and therefore reveals the merging histories of galaxies.

Employing the data from orbital periods and masses of extra-solar planets in 166 multiple planetary systems, the period-ratio and mass-ratio of adjacent planet pairs are studied. The correlation between the period-ratio and mass-ratio is confirmed and found to have a correlation coefficient of 0.5303 with a 99% confidence interval (0.3807, 0.6528). A comparison with the distribution of synthetic samples from a Monte Carlo simulation reveals the imprint of planet-planet interactions on the formation of adjacent planet pairs in multiple planetary systems.

The existence and uniqueness of equilibrium points, including Lagrange Points and Jiang-Yeh Points, of a galactic system with supermassive binary black holes embedded in a centrally cuspy galactic halo are investigated herein. Differing from the previous results of non-cuspy galactic profiles that Jiang-Yeh Points only exist under a particular condition, it is found here that the Lagrange Points, L2, L3, L4 and L5, Jiang-Yeh Points, JY1 and JY2, exist under general conditions. The stability analysis shows that L2, L3, JY1 and JY2 are unstable. However, L4 and L5 are only unstable when the galactic total mass is smaller than a critical mass; otherwise they become neutrally stable centers. These results will be important for further studies on the cores of early-type galaxies.

The dynamics of galactic systems with central binary black holes is studied. The model is a modification from the restricted three body problem, in which a galactic potential is added as an external potential. Considering the case with an equal mass binary black holes, the conditions of existence of equilibrium points, including Lagrange Points and additional new equilibrium points, i.e. Jiang-Yeh Points, are investigated. A critical mass is discovered to be fundamentally important. That is, Jiang-Yeh Points exist if and only if the galactic mass is larger than the critical mass. The stability analysis is performed for all equilibrium points. The results that Jiang-Yeh Points are unstable could lead to the core formation in the centers of galaxies.

Five newly observed transit light curves of the TrES-3 planetary system are presented. Together with other light curve data from literature, 23 transit light curves in total, which cover an overall timescale of 911 epochs, have been analyzed through a standard procedure. From these observational data, the system's orbital parameters are determined and possible transit timing variations are investigated. Given that a null transit-timing-variation produces a fit with reduced chi^2=1.52, our results agree with previous work, that transit timing variations might not exist in these data. However, a 1-frequency oscillating transit-timing-variation model, giving a fit with a reduced chi^2=0.93, does possess a statistically higher probability. It is, thus, concluded that future observations and dynamical simulations for this planetary system will be very important.

Axis-symmetric ring systems can be identified from the new catalog of collisional ring galaxies in Madore et al. (2009). These are O-type-like collisional ring galaxies. Head-on collisions by dwarf galaxies moving along the symmetric axis were performed through N-body simulations to address their origins. It was found that the simulations with smaller initial relative velocities between two galaxies, or the cases with heavier dwarf galaxies, could produce rings with higher density contrasts. There are more than one generation of rings in one collision and the lifetime of any generation of rings is about one dynamical time. It was concluded that head-on penetrations could explain these O-type-like ring galaxies identified from the new catalog in Madore et al. (2009), and the simulated rings resembling the observational O-type-like collisional rings are those at the early stage of one of the ring-generations.

This paper studied the structures of debris discs, focusing on the conditions that can form an axisymmetric-looking outer disc from systems with inner clumps. The main conclusion was that as long as the dominated dust grains are smaller than the blowout size, it is easy to form an axisymmetric-looking outer debris disc, which is part of a quasi-steady state of the whole system. This quasi-steady state is established through the balance between grain generations and a continuous out-going grain flow. Assuming there is an event that starts planetesimal collisions and the corresponding grain generations, this balance can be approached in a few thousand years. This result suggested that a quasi-steady-state picture could solve the possible mass budget problem of Vega's outer debris disc.

Due to the angular momentum exchange with planetesimals, Neptune might have migrated outward to the current position, and captured many Kuiper belt objects (KBOs) into resonances. We set up a semi-analytic model to simulate the outward migration of Neptune, and the processes of resonant captures. Our model can naturally explain Neptune's currently observed semi-major axis and eccentricity. The results show that the current population ratio between 3:2 and 2:1 is mainly due to the original density distribution of KBOs, which might be related to drag-induced inward migrations of proto-KBOs.

Through the international collaborators, we recently established a network of existing and working meter-class telescopes to look for planetary transit events. As a first step, we focus on the TrES3 system, and conclude that there could be some level of transit timing variations.

Employing a catalog of 175 extrasolar planets (exoplanets) detected by the Doppler-shift method, we constructed the independent and coupled mass-period functions. It is the first time in this field that the selection effect is considered in the coupled mass-period functions. Our results are consistent with those in Tabachnik and Tremaine (2002) with the major differences that we obtain a flatter mass function but a steeper period function. Moreover, our coupled mass-period functions show that about 2.5 percent of stars would have a planet with mass between Earth Mass and Neptune Mass, and about 3 percent of stars would have a planet with mass between Neptune Mass and Jupiter Mass.

Motivated by the observations on the intra-cluster light (ICL) and inter-galactic stellar populations, n-body simulations are used to model the galactic merging events as a goal to investigate the production and distribution of gravitational unbound populations (GUPs). Both the parabolic and hyperbolic mergers are considered and each category includes six models with different relative orientations between two galaxies. Our results show that there are more (about a factor of two) GUP after a hyperbolic merging event than after a parabolic one. In general, depending on the relative orientation and also the relative velocity of the two galaxies in a merging pair, a head-on collision of a galaxy pair would only make a tiny fraction (less than one percent) of the initial stellar mass become luminous GUP but a considerable fraction (eight to fourteen percent) of the dark matter become dark GUP.

Using the period and mass data of two hundred and seventy-nine extrasolar planets, we have constructed a coupled period-mass function through the non-parametric approach. This analytic expression of the coupled period-mass function has been obtained for the first time in this field. Moreover, due to a moderate period-mass correlation, the shapes of mass/period functions vary as a function of period/mass. These results of mass and period functions give way to two important implications: (1) the deficit of massive close-in planets is confirmed, and (2) the more massive planets have larger ranges of possible semi-major axes. These interesting statistical results will provide important clues into the theories of planetary formation.

This study employed grain dynamic models to examine the density distribution of debris discs, and discussed the effects of the collisional time-intervals of asteroidal bodies, the maximum grain sizes, and the chemical compositions of the dust grains of the models, in order to find out whether a steady out-moving flow with an 1/R profile could be formed. The results showed that a model with new grains every 100 years, a smaller maximum grain size, and a composition C400 has the best fit to the 1/R profile because: (1) the grains have larger values of beta on average,therefore, they can be blown out easily; (2) the new grains are generated frequently enough to replace those have been blown out. With the above two conditions, some other models can have a steady out-moving flow with an approximate 1/R profile. However, those models in which new grains are generated every 1000 years have density distributions far from the profile of a continuous out-moving flow. Moreover, the analysis on the signatures of planets in debris discs showed that there are no indications when a planet is in a continuous out-moving flow, however, the signatures are obvious in a debris disc with long-lived grains.

Through Spitzer Space Telescope's observations, Su et al. (2005) show that the Vega debris disc is dominated by grains which are small enough to be blown out by radiation pressure. This implies the lifetime of Vega debris disc's grains is relatively short, about 1000 years, and a continuous dust production is necessary to maintain the observed debris disc. However, Krivov et al. (2006)'s theoretical calculations show that the Vega debris disc is dominated by 10 micro-meter grains, which would be in bound orbits and thus long-lived, provided that the disc is in a steady state. In order to solve the above contradiction, through dynamical simulations, we determine the grains' orbital evolutions and density profiles and seek a model of size distribution which can reproduce the observed surface brightness. Our results show that a self-consistent dynamical model with a 1/R disc density profile can be constructed when the grains have a power-law size distribution. Moreover, both types of models, dominated by short-lived and long-lived grains, are consistent with the observational data.

Using the period and mass data of two hundred and seventy-nine extrasolar planets, we have constructed a coupled period-mass function through the non-parametric approach. This analytic expression of the coupled period-mass function has been obtained for the first time in this field. Moreover, due to a moderate period-mass correlation, the shapes of mass/period functions vary as a function of period/mass. These results of mass and period functions give way to two important implications: (1) the deficit of massive close-in planets is confirmed, and (2) the more massive planets have larger ranges of possible semi-major axes. These interesting statistical results will provide important clues into the theories of planetary formation.

We have performed smoothed particle hydrodynamics (SPH) simulations to study the response of the central kiloparsec region of a gaseous disk to the imposition of nonaxisymmetric bar potentials. The model galaxies are composed of the three axisymmetric components (halo, disk, and bulge) and a non-axisymmetric bar. These components are assumed to be invariant in time in the frame corotating with the bar. The potential of spherical $\gamma$-models of Dehnen is adopted for the bulge component whose density varies as $r^{-\gamma}$ near the center and $r^{-4}$ at larger radiiand hence, possesses a central density core for $\gamma = 0$ and cusps for $\gamma > 0$. Since the central mass concentration of the model galaxies increases with the cusp parameter $\gamma$, we have examined here the effect of the central mass concentration by varying the cusp parameter $\gamma$ on the mechanism responsible for the formation of the symmetric two-armed nuclear spirals in barred galaxies. Our simulations show that the symmetric two-armed nuclear spirals are formed by hydrodynamic spiral shocks driven by the gravitational torque of the bar for the models with $\gamma = 0$ and 0.5. On the other hand, the symmetric two-armed nuclear spirals in the models with $\gamma=1$ and 1.5 are explained by gas density waves. Thus, we conclude that the mechanism responsible for the formation of the symmetric two-armed nuclear spirals in barred galaxies changes from the hydrodynamic shocks to the gas density waves when the central mass concentration increases from $\gamma = 0$ to 1.5.

We use the usual method of Schwarzschild to construct self-consistent solutions for the triaxial de Zeeuw & Carollo (1996) models with central density cusps. ZC96 models are triaxial generalisations of spherical $\gamma$-models of Dehnen whose densities vary as $r^{-\gamma}$ near the center and $r^{-4}$ at large radii and hence, possess a central density core for $\gamma=0$ and cusps for $\gamma > 0$. We consider four triaxial models from ZC96, two prolate triaxials: $(p, q) = (0.65, 0.60)$ with $\gamma = 1.0$ and 1.5, and two oblate triaxials: $(p, q) = (0.95, 0.60)$ with $\gamma = 1.0$ and 1.5. We compute 4500 orbits in each model for time periods of $10^{5} T_{D}$. We find that a large fraction of the orbits in each model are stochastic by means of their nonzero Liapunov exponents. The stochastic orbits in each model can sustain regular shapes for $\sim 10^{3} T_{D}$ or longer, which suggests that they diffuse slowly through their allowed phase-space. Except for the oblate triaxial models with $\gamma =1.0$, our attempts to construct self-consistent solutions employing only the regular orbits fail for the remaining three models. However, the self-consistent solutions are found to exist for all models when the stochastic and regular orbits are treated in the same way because the mixing-time, $\sim10^{4} T_{D}$, is shorter than the integration time, $10^{5} T_{D}$. Moreover, the ``fully-mixed'' solutions can also be constructed for all models when the stochastic orbits are fully mixed at 15 lowest energy shells. Thus, we conclude that the self-consistent solutions exist for our selected prolate and oblate triaxial models with $\gamma = 1.0$ and 1.5.

In addition to fitting the data of 233 extra-solar planets with power laws, we construct a correlated mass-period distribution function of extrasolar planets, as the first time in this field. The algorithm to generate a pair of positively correlated beta-distributed random variables is introduced and used for the construction of correlated distribution functions. We investigate the mass-period correlations of extrasolar planets both in the linear and logarithm spaces, determine the confidence intervals of the correlation coefficients, and confirm that there is a positive mass-period correlation for the extrasolar planets. In addition to the paucity of massive close-in planets, which makes the main contribution on this correlation, there are other fine structures for the data in the mass-period plane.

It is known that the discs are detected for some of the extra-solar planetary systems. It is also likely that there was a disc mixing with planets and small bodies while our Solar System was forming. From our recent results, we conclude that the discs play two roles: the gravity makes planetary systems more chaotic and the drag makes planetary systems more resonant.

We propose a new mechanism of drag-induced resonant capture, which can explain the resonant Kuiper Belt Objects in a natural way. A review and comparison with the traditional mechanism of sweeping capture by the migrating Neptune will be given.

The restricted three body problem is well-known and very important for dynamics of binary, multiple stars and also planetary systems. We extend the classical version of this problem to the situation that there are some external forces from the belt. We find that both the equilibrium points and solution curves become quite different from the classical case. We also determine the values of Lyapunov Exponent for some important orbits.

More than 100 extrasolar planets have been discovered since 1990s. Different from the solar system, these planets' orbital eccentricities cover a huge range from 0 to 0.7. Incidently, the first Kuiper Belt Object was discovered in 1992. Thus, an interesting and important question will be whether extrasolar planetary systems could have structures like Kuiper Belt or asteroid belt. We investigate the stability of these planetary systems with different orbital eccentricities by the similar procedures in Rabl & Dvorak (1988) and Holman & Wiegert (1999). We claim that most extrasolar planetary systems can have their own belts at the outer regions. However, we find that the orbits with high--eccentricity is very powerful in depletion of these populations.

The observational properties of brown dwarfs pose challenges to the theory of star formation. Because their mass is much smaller than the typical Jeans mass of interstellar clouds, brown dwarfs are most likely formed through secondary fragmentation processes, rather than through the direct collapse of a molecular cloud core. In order to prevent substantial post-formation mass accretion, young brown dwarfs must leave the high density formation regions in which they form. We propose here that brown dwarfs are formed in the optically thin outer regions of circumbinary disks. Through post-formation dynamical interaction with their host binary stars, young brown dwarfs are either scattered to large distance or removed, with modest speed, from their cradles.

Motivated by Papadakis (2005a, 2005b), we study a Chermnykh-like problem, in which an additional gravitational potential from the belt is included. In addition to the usual five equilibrium points (three collinear and two triangular points), there are some new equilibrium points for this system. We studied the conditions for the existence of these new equilibrium points both analytically and numerically.

Following Papadakis (2005)'s numerical exploration of the Chermnykh's problem, we here study a Chermnykh-like problem motivated by the astrophysical applications. We find that both the equilibrium points and solution curves become quite different from the ones of the classical planar restricted three-body problem. In addition to the usual Lagrangian points, there are new equilibrium points in our system. We also calculate the Lyapunov Exponents for some example orbits. We conclude that it seems there are more chaotic orbits for the system when there is a belt to interact with.

We use both the conventional and more recently developed methods of cluster analysis to study the data of extra-solar planets. Using the data set with planetary mass M, orbital period P, and orbital eccentricity e, we investigate the possible clustering in the ln M, ln P, ln P-ln M, e, and ln P-e spaces. There are two main implications: (1) mass distribution is continuous and (2) orbital population could be classified into three clusters, which correspond to the exoplanets in the regimes of tidal, on-going tidal and disc interaction, respectively.

We study the migration and resonant capture of planetesimals in a planetary system consisting of a gaseous disc analogous to the primordial solar nebula and a Neptune-like planet. Using a simple treatment of the drag force we find that planetesimals are mainly trapped in the 3:2 and 2:1 resonances and that the resonant populations are correlated with the gaseous drag strength in a sense that the 3:2 resonant population increases with the stronger gaseous drag, but the 2:1 resonant population does not. Since planetesimals can lead to the formation of larger bodies similar to asteroids and Kuiper Belt Objects, the gaseous drag can play an important role in the configuration of a planetary system.

It has been an interesting question that why there are one-third of Kuiper Belt Objects (KBOs) trapped into the 3:2 resonance but, in contrast, only several KBOs are claimed to be associated with the 2:1 resonance. In a model proposed by Zhou et al. (2002), the stochastic outward migration of the Neptune could reduce the number of particles in the 2:1 resonance and thus the objects in the 3:2 resonance become more distinct. As a complementary study, we investigate the effect of proto-stellar discs on the resonance capture. Our results show that the gaseous drag of a proto-stellar disc can trap KBOs into the 3:2 resonance rather easily. In addition, no objects are captured into the 2:1 resonance in our simulation.

Following Tancredi, Sanchez and Roig (2001)'s criteria of chaos, two ways of setting initial velocities are used in the numerical surveys to explore the possible chaotic and regular orbits for the disc-star-planet systems. We find that the chaotic boundary does not depend much on the disc mass for Type I initial condition, but can change a lot for different disc masses for Type II initial condition. A few sample orbits are further studied. Both Poincare surface of section and the Lyapounov Exponent Indicator are calculated and they are consistent with each other. We also find that the influence from the disc can change the locations of equilibrium points and the orbital behaviors for both types of initial conditions. Because the chaotic orbits are less likely to become the stable resonant orbits, we conclude that the proto-stellar disc shall play important roles for the capture and depletion histories of resonant orbits of both Asteroid Belt and Kuiper Belt during the formation of Solar System.

The observational properties of brown dwarfs pose challenges to the theory of star formation. Because their mass is much smaller than the typical Jeans mass of interstellar clouds, brown dwarfs are most likely formed through secondary fragmentation processes, rather than through the direct collapse of a molecular cloud core. In order to prevent substantial post-formation mass accretion, young brown dwarfs must leave the high density formation regions in which they form. We propose here that brown dwarfs are formed in the circumbinary disks. Through post-formation dynamical interaction with their host binary stars, young brown dwarfs are either scattered to large distance or removed, with modest speed, from their cradles.

The orbital evolution and stability of planetary systems with interaction from the belts is studied using the standard phase-plane analysis. In addition to the fixed point which corresponds to the Keplerian orbit, there are other fixed points around the inner and outer edges of the belt. Our results show that for the planets, the probability to move stably around the inner edge is larger than the one to move around the outer edge. It is also interesting that there is a limit cycle of semi-attractor for a particular case. Applying our results to the Solar System, we find that our results could provide a natural mechanism to do the orbit rearrangement for the larger Kuiper Belt Objects and thus successfully explain the absence of these objects beyond 50 AU.

It has been shown that there is a possible mass-period correlation for extrasolar planets from the current observational data and this correlation is, in fact, related to the absence of massive close-in planets, which are strongly influenced by the tidal interaction with the central star. We confirm that the model in P\"atzold & Rauer (2002) is a good approximation for the explanation of the absence of massive close-in planets. We thus further determine the minimum possible semimajor axis for these planets to be detected during their lifetime and also study their migration time scale at different semimajor axes by the calculations of tidal interaction. We conclude that the mass-period correlation at the time when these planets were just formed was less tight than it is now observed if these orbital migrations are taken into account.

It was proposed by Blitz et al. (1999) that High Velocity Clouds (HVCs) are remnants of Local Group formation and the average distance of HVCs is 1 Mpc, which is the result of a simple dynamical calculation leading to match the observed HVCs distributions. However, in this paper, we clearly show that fitting the observed HVCs distributions by a dynamical calculation cannot provide any constraints on the average distance of HVCs. With our choices of initial conditions, the observational results in Wakker & van Woerden (1991) are produced in our simulations for the models of both Galactic and extragalactic origins. Moreover, because Zwaan & Briggs (2000) reported that they failed to locate any extragalactic counterparts of the Local Group HVCs in a blind HI 21-cm survey of extragalactic groups, we propose to use ``remnants of galactic disc formation'' as the modification for the picture of ``remnants of galaxy group formation'' in Blitz et al. (1999) and thus reduce the possible average distances of HVCs to be about or less than a few hundred kpc.

The dynamical systems of planet-belt interaction are studied by the fixed-point analysis and the bifurcation of solutions on the parameter space is discussed. For most cases, our analytical and numerical results show that the locations of fixed points are determined by the parameters and these fixed points are either structurally stable or unstable. In addition to that, there are two special fixed points: the one on the inner edge of the belt is asymptotically stable and the one on the outer edge of the belt is unstable. This is consistent with the observational picture of Asteroid Belt between the Mars and Jupiter: the Mars is moving stablely close to the inner edge but the Jupiter is quite far from the outer edge.

A simple dynamical model is employed to study the possible orbital evolution of scattered planets and phase plane analysis is used to classify the parameter space and solutions. Our results reconfirm that there is always an increase in eccentricity when the planet was scattered to migrate outward when the initial eccentricity is zero. Applying our study on the Solar System and considering the existence of the Kuiper Belt, this conclusion implies that Neptune was dynamically coupled with the Kuiper Belt in the early phase of the Solar System, which is consistent with the simulational model in Thommes, Duncan & Levison (1999).

The bright F8 V solar-type star upsilon Andromedae has recently been reported to have a system of three planets of Jovian masses. In order to investigate the orbital stability and mutual gravitational interactions among these extrasolar planets, both forward and backward integrations from the latest observed orbital elements for all three planets' orbits have been performed under the coplanar assumption. We reconfirm that the middle and the outer planet have strong interaction leading to large time variations in the eccentricities of these planets, which was shown by the previous studies. However, we discuss the validity of the ignorance of the innermost planet. We argue that this planetary system is likely to be stable and oscillate around current orbital elements since it was formed. We suggest that one possible way to produce these orbital elements: the innermost planet has very low eccentricity but the outermost planet has high eccentricity could be the interaction with the protostellar disc.

The size distribution is determined directly from the current known Kuiper Belt Objects (KBOs). We found that there is a peak for the size distribution around 230 km in diameter. For the objects larger than 230 km, the KBOs populate as $N(s) \propto s^{-4}$ as we see in Jewitt et al. (1998). For the objects smaller than 230 km, we found $N(s) \propto s^{3}$. This result from the current observational data is new and very different from the conventional concept of the KBO size distribution. Though the true size distribution might change from what we have after more KBOs be discovered, we argue that the existence of a peak around 150 km to 250 km is likely.

Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N-body code is used to construct the first models of the Milky Way - Sgr Dwarf system in which both the Milky Way and the Sgr Dwarf are represented by full N-body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Dwarf's orbit that enable us to explore a wider parameter space than is accessible to the N-body models. We conclude that the extant data on the Dwarf are compatible with a wide range of orbital histories. At one extreme the Dwarf initially possesses 10^{11} Solar Mass and starts from a Galactocentric distance 200 kpc. At the other extreme the Dwarf starts with 10^9 Solar Mass and Galactocentric distance 60 kpc, similar to its present apocentric distance. In all cases the Dwarf is initially dark-matter dominated and the current velocity dispersion of the Dwarf's dark matter is tightly constrained to be 21 km/s. This number is probably compatible with the smaller measured dispersion of the Dwarf's stars because of (a) the dynamical difference between dark and luminous matter, and (b) velocity anisotropy.

N-body simulations show that when infall reorientates the outer parts of a galactic halo by several degrees per Gyr, a self-gravitating disk that is embedded in the halo develops an integral-sign warp that is comparable in amplitude to observed warps. Studies of angular-momentum acquisition suggest that the required rate of halo reorientation is realistic for galaxies like the Milky Way.

We study the persistence of warps in galactic discs in the presence of massive halos. A disc is approximated by a set of massive rings, while a halo is represented by a conventional n-body simulation. We confirm the conclusion of Nelson & Tremaine (1995) that a halo responds strongly to an embedded precessing disc. This response invalidates the approximations made by in the derivation of classical `modified tilt' modes. We show that the response of the halo causes the line of nodes of a disc that starts from a modified tilt mode to wind up within a few dynamical times. We explain this finding in terms of the probable spectrum of true normal modes of a combined disc-halo system.