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.

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.

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.

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.

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.

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.

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.

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.

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.

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 explore conventional Neptune migration model with one additional planet of mass at 0.1-2.0 Me. This planet inhabited in the 3:2 mean motion resonance with Neptune during planet migration epoch, and then escaped from the Kuiper belt when Jovian planets parked near the present orbits. Adding this extra planet and assuming the primordial disk truncated at about 45 AU in the conventional Neptune migration model, it is able to explain the complex structure of the observed Kuiper belt better than the usual Neptune migration model did in several respects. However, numerical experiments imply that this model is a low-probability event. In addition to the low probability, two features produced by this model may be inconsistent with the observations. They are small number of low-inclination particles in the classical belt, and the production of a remnant population with near-circular and low-inclination orbit within a = 50-52 AU. According to our present study, including one extra planet in the conventional Neptune migration model as the scenario we explored here may be unsuitable because of the low probability, and the two drawbacks mentioned above, although this model can explain better several features which is hard to produce by the conventional Neptune migration model. The issues of low-probability event and the lack of low-inclination KBOs in the classical belt are interesting and may be studied further under a more realistic consideration.

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.

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.

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.

To derive physical properties of the neutron star surface with observed spectra, a realistic model spectrum of neutron star surface emission is essential. Limited by computing resources, a full computation of the radiative transfer equations without the diffusion approximation has been conducted up to date only for the case of local magnetic fields being perpendicular to the stellar surface. In this paper we report the full-computation result for an arbitrary field direction. For comparison we also compute the radiative transfer equation using the diffusion approximation. For a given effective temperature, the computed spectrum with the diffusion approximation is always softer than that of a full computation at a non-negligible level. It leads to an over-estimate of the effective temperature if the diffusion approximation spectrum is employed in the spectral fitting. Other characteristics for different magnetic field orientations, such as the beaming pattern of the two polarization modes and the structure of the atmosphere, are also discussed.

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.

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

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