In this work, we investigate configuration formation of two inner terrestrial planets near mean motion resonance (MMRs) induced by the perturbation of a distant gas-giant for the Kepler-68 system, by conducting thousands of numerical simulations. The results show that the formation of terrestrial planets is relevant to the speed of Type I migration, the mass of planets, and the existence of giant planet. The mass and eccentricity of the giant planet may play a crucial role in shaping the final configuration of the system. The inner planet pair can be trapped in 5:3 or 7:4 MMRs if the giant planet revolves the central star with an eccentric orbit, which is similar to the observed configuration of Kepler-68. Moreover, we find that the eccentricity of the middle planet can be excited to roughly 0.2 if the giant planet is more massive than 5 $M_J$,otherwise the terrestrial planets are inclined to remain near-circular orbits. Our study may provide a likely formation scenario for the planetary systems that harbor several terrestrial planets near MMRs inside and one gas-giant exterior to them.

In this work, we investigate the size, thermal inertia, surface roughness and geometric albedo of 10 Vesta family asteroids by using the Advanced Thermophysical Model (ATPM), based on the thermal infrared data acquired by mainly NASA's Wide-field Infrared Survey Explorer (WISE). Here we show that the average thermal inertia and geometric albedo of the investigated Vesta family members are 42 $\rm J m^{-2} s^{-1/2} K^{-1}$ and 0.314, respectively, where the derived effective diameters are less than 10 km. Moreover, the family members have a relatively low roughness fraction on their surfaces. The similarity in thermal inertia and geometric albedo among the V-type Vesta family member may reveal their close connection in the origin and evolution. As the fragments of the cratering event of Vesta, the family members may have undergone similar evolution process, thereby leading to very close thermal properties. Finally, we estimate their regolith grain sizes with different volume filling factors.

The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au.

High spatial resolution observations of protoplanetary disks (PPDs) by ALMA have revealed many details that are providing interesting constraints on the disk physics as well as dust dynamics, both of which are essential for understanding planet formation. We carry out high-resolution, 2D global hydrodynamic simulations, including the effects of dust feedback, to study the stability of dusty rings. When the ring edges are relatively sharp and the dust surface density becomes comparable to the gas surface density, we find that dust feedback enhances the radial gradients of both the azimuthal velocity profile and the potential vorticity profile at the ring edges. This eventually leads to instabilities on meso-scales (spatial scales of several disk scale heights), causing dusty rings to be populated with many compact regions with highly concentrated dust densities on meso-scales. We also produce synthetic dust emission images using our simulation results and discuss the comparison between simulations and observations.

The Insight-Hard X-ray Modulation Telescope (Insight-HXMT) is a broad band X-ray and gamma-ray (1-3000 keV) astronomy satellite. The High Energy X-ray telescope (HE) is one of its three main telescopes. The main detector plane of HE is composed of 18 NaI(Tl)/CsI(Na) phoswich detectors, where NaI(Tl) serves as primary detector to measure ~ 20-250 keV photons incident from the field of view (FOV) defined by the collimators, and CsI(Na) is used as an active shield detector to NaI(Tl) by pulse shape discrimination. CsI(Na) is also used as an omnidirectional gamma-ray monitor. The HE collimators have a diverse FOV: 1.1{\deg}x 5.7{\deg} (15 units), 5.7{\deg}x 5.7{\deg} (2 units) and blocked (1 unit), thus the combined FOV of HE is about 5.7{\deg}x 5.7{\deg}. Each HE detector has a diameter of 190 mm, resulting in the total geometrical area of about 5100 cm_2. The energy resolution is ~15% at 60 keV. The timing accuracy is better than 10 {\mu}s and dead-time for each detector is less than 10 {\mu}s. HE is devoted to observe the spectra and temporal variability of X-ray sources in the 20-250 keV band either by pointing observations for known sources or scanning observations to unveil new sources, and to monitor the gamma-ray sky in 0.2-3 MeV. This paper presents the design and performance of the HE instruments. Results of the on-ground calibration experiments are also reported.

In this paper, we investigate thermophysical characteristics of near-Earth asteroid (341843) 2008 EV5, based on our improved Advanced Thermal Physical Model (ATPM) by considering the contribution of sunlight-reflection for rough surface, along with four wavebands observations from Wide-field Infrared Survey Explorer (WISE) and the radar-derived shape model. Here we derive that 2008 EV5 has a relatively low thermal inertia of $\Gamma=110 _{-12}^{+40}\rm~J m^{-2} s^{-1/2} K^{-1}$ but a high roughness fraction. The geometric albedo and effective diameter are then constrained to be $p_v=0.095_{-0.003}^{+0.016}$ and $D_{\rm eff}=431_{-33}^{+6}\rm~m$, respectively. The low thermal inertia indicates that 2008 EV5 may have undergone sufficient space weathering over secular evolution. The high roughness may have resemblances to the appearances of Bennu and Ryugu recently observed by spacecrafts, where a great number of boulders are widely distributed on the asteroid's surface. Moreover, we numerically perform 1000 backward simulations of 2008 EV5's cloned orbits within 1 Myr to explore its origin, and present a probability of $\sim6.1\%$ that the asteroid originates from the main belt. Finally, we estimate that the mean grain size of the surface ranges from 0.58 to 1.3 mm, and infer that it is unlikely to find water ice on most area of 2008 EV5, but there may exist water ice on high-latitudes near polar region.

Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHa 330, are unlikely to be induced by the candidate vortices in them.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

We search a large parameter space of the LkCa 15's disk density profile to fit its observed radial intensity profile of $^{12}$CO (J = 3-2) obtained from ALMA. The best-fit model within the parameter space has a disk mass of 0.1 $M_{\odot}$ (using an abundance ratio of $^{12}$CO/H$_2$ $=$ 1.4 $\times 10^{-4}$ in mass), an inner cavity of 45 AU in radius, an outer edge at $\sim$ 600 AU, and a disk surface density profile follows a power-law of the form $\rho_r \propto r^{-4}$. For the disk density profiles that can lead to a small reduced $\chi^2$ of goodness-of-fit, we find that there is a clear linear correlation between the disk mass and the power-law index $\gamma$ in the equation of disk density profile. This suggests that the $^{12}$CO disk of LkCa 15 is optically thick and we can fit its $^{12}$CO radial intensity profile using either a lower disk mass with a smaller $\gamma$ or a higher disk mass with a bigger $\gamma$. By comparing the $^{12}$CO channel maps of the best-fit model with disk models with higher or lower masses, we find that a disk mass of $\sim$ 0.1 $M_{\odot}$ can best reproduce the observed morphology of the $^{12}$CO channel maps. The dust continuum map at 0.87 mm of the LkCa 15 disk shows an inner cavity of the similar size of the best-fit gas model, but its out edge is at $\sim$ 200 AU, much smaller than the fitted gas disk. Such a discrepancy between the outer edges of the gas and dust disks is consistent with dust drifting and trapping models.

The MIRO instrument's remote sensing capability is integral in constraining water density, temperature and velocity fields in the coma of 67P/Churyumov-Gersimenko. Our aim is to quantify the contribution to the water density from all facets inside and outside the field-of-view (FOV) of MIRO, in both, nadir and limb geometries. This information is crucial for understanding the MIRO derived coma production rates and their relation to the nucleus characteristics, and inherent spatial resolution of the data. This study relies on a detailed 3D nucleus shape model, illumination conditions and the pointing information of the viewing geometry. With these parameters, we can evaluate the relative contribution of water density originating from facets directly inside the MIRO beam as well outside of the beam as a function of distance along the MIRO line-of-sight. We also calculate the ratio of in-beam versus out-of-beam number density. We demonstrate that despite the rather small MIRO field-of-view there is only a small fraction of molecules that originate from facets within the MIRO beam. This is true for nadir, but also translated into the limb observing geometry. The MIRO instrument cannot discriminate active from non-active regions directly from observations. This study also suggests that the beam averaged solar incidence angle, local time and mean normal vectors are not necessary related to molecules within the MIRO beam. These results also illustrate why the 1D spherical Haser model can be applied with relative success for analyzing the MIRO data (and generally any Rosetta measurements). The future possibilities of constraining gas activity distribution on the surface should use 3D codes extracting information from the MIRO spectral line shapes which contain additional information. The presented results are in fact applicable to all relevant instruments onboard Rosetta.

Our aim is to investigate early activity (July 2014) of 67P/CG with 3D coma and radiative transfer modeling of MIRO measurements, accounting for nucleus shape, illumination, and orientation of the comet. We investigate MIRO line shape information for spatial distribution of water activity on the nucleus during the onset of activity. During this period we show that MIRO line shape have enough information to clearly isolate contribution from Hapi and Inhotep independently, and compare it to the nominal case of activity from the entire illuminated surface. We also demonstrate that spectral line shapes differ from the 1D model for different viewing geometries and coma conditions relevant to this study. Specifically, line shapes are somewhat sensitive to the location of the terminator in the coma. At last, fitting the MIRO observations we show that the Imhotep region (possible extended source of H$_{2}$O due to CO$_{2}$ activity) contributes only a small fraction of the total number of water molecules into MIRO beam in the early activity. On the other hand, a strong enhancement of water activity from the Hapi region seems required to fit the MIRO line shapes. This is consistent with earlier Rosetta results. Nevertheless, within the assumption of our coma and surface boundary conditions, we cannot get a reasonable fit to all MIRO mapping observations in July 2014, which may illustrate that a more sophisticated coma model or more accurate temperature/velocity distribution is needed.

Proxima Centauri was recently discovered to host an Earth-mass planet of Proxima b, and a 215-day signal which is probably a potential planet c. In this work, we investigate the dynamical evolution of the Proxima Centauri system with the full equations of motion and semi-analytical models including relativistic and tidal effects. We adopt the modified Lagrange-Laplace secular equations to study the evolution of eccentricity of Proxima b, and find that the outcomes are consistent with those from the numerical simulations. The simulations show that relativistic effects have an influence on the evolution of eccentricities of planetary orbits, whereas tidal effects primarily affects the eccentricity of Proxima b over long timescale. Moreover, using the MEGNO (the Mean Exponential Growth factor of Nearby Orbits) technique, we place dynamical constraints on orbital parameters that result in stable or quasi-periodic motions for coplanar and non-coplanar configurations. In the coplanar case, we find that the orbit of Proxima b is stable for the semi-major axis ranging from 0.02 au to 0.1 au and the eccentricity being less than 0.4. This is where the best-fitting parameters for Proxima b exactly fall. Additional simulations show that the robust stability of this system would favor an eccentricity of Proxima b less than 0.45 and that of Proxima c below 0.65. In the non-coplanar case, we find that mutual inclinations of two planets must be lower than $50^{\circ}$ in order to provide stability. Finally, we estimate the mass of Proxima c to be $3.13~{M_ \oplus } \le {m_c} \le 70.7~{M_ \oplus }$ when $1.27~{M_ \oplus } \le {m_b} \le 1.6~{M_ \oplus }$, if ${i_{mutual}} \le {50^{\circ}}$ and $\Delta \Omega = 0^{\circ}$.

Several nearby protoplanetary disks have been observed to display large scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas over-densities and are expected to concentrate dust particles with Stokes number around unity. This process might catalyze the formation of planetesimals. Whereas recent observations showed that dust crescent are indeed regions where millimeter-size particles have abnormally high concentration relative to the gas and smaller grains, no observations have yet shown that the gas within the crescent region counter-rotates with respect to the protoplanetary disk. Here we investigate the detectability of anticyclonic features through measurement of the line-of-sight component of the gas velocity obtained with ALMA. We carry out 2D hydrodynamic simulations and 3D radiative transfer calculation of a protoplanetary disk characterized by a vortex created by the tidal interaction with a massive planet. As a case study, the disk parameters are chosen to mimic the IRS 48 system, which has the most prominent crescent observed to date. We generate synthetic ALMA observations of both the dust continuum and 12CO emission around the frequency of 345 GHz. We find that the anticyclonic features of vortex are weak but can be detected if both the source and the observational setup are properly chosen. We provide a recipe for maximizing the probability to detect such vortex features and present an analysis procedure to infer their kinematic properties.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

Although several S-type and P-type planets in binary systems were discovered in past years, S-type planets have not yet been found in close binaries with an orbital separation not more than 5 au. Recent studies suggest that S-type planets in close binaries may be detected through high-accuracy observations. However, nowadays planet formation theories imply that it is difficult for S-type planets in close binaries systems to form in situ. In this work, we extensively perform numerical simulations to explore scenarios of planet-planet scattering among circumbinary planets and subsequent tidal capture in various binary configurations, to examine whether the mechanism can play a part in producing such kind of planets. Our results show that this mechanism is robust. The maximum capture probability is $\sim 10\%$, which can be comparable to the tidal capture probability of hot Jupiters in single star systems. The capture probability is related to binary configurations, where a smaller eccentricity or a low mass ratio of the binary will lead to a larger probability of capture, and vice versa. Furthermore, we find that S-type planets with retrograde orbits can be naturally produced via capture process. These planets on retrograde orbits can help us distinguish in situ formation and post-capture origin for S-type planet in close binaries systems. The forthcoming missions (PLATO) will provide the opportunity and feasibility to detect such planets. Our work provides several suggestions for selecting target binaries in search for S-type planets in the near future.

In order to obtain the substantial information about the surface physics and thermal property of the target asteroid (162173) 1999 JU3, which will be visited by Hayabusa 2 in a sample return mission, with the Advanced Thermal Physical Model (ATPM) we estimate the possible thermal inertia distribution over its surface, and infer the major material composition of its surface materials. In addition, the effective diameter and geometric albedo are derived to be $D_{\rm eff}=1.13\pm0.03\rm~km$, $p_{\rm v}=0.042\pm0.003$, respectively, and the average thermal inertia is estimated to be about $(300\pm50)\rm~J\cdot m^{-2}\cdot s^{-0.5}\cdot K^{-1}$. According to the derived thermal inertia distribution, we infer that the major area on the surface of the target asteroid may be covered by loose materials, such as rock debris, sands, and so on, but few bare rocks may exist in a very small region. In this sense, the sample return mission of Hayabusa 2 is feasible, when it is performed successfully, it will certainly bring significant scientific information to the research of asteroids.

The planets with a radius $<$ 4 $R$$_\oplus$ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius $<$ 4 $R_\oplus$ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ.

Two photometric follow-up transit (primary eclipse) observations on WASP-43 b and four observations on TrES-3 b are performed using the Xuyi Near-Earth Object Survey Telescope. After differential photometry and light curve analysis, the physical parameters of the two systems are obtained and are in good match with the literature. Combining with transit data from a lot of literature, the residuals ($O-C$) of transit observations of both systems are fitted with the linear and quadratic functions. With the linear fitting, the periods and transit timing variations (TTVs) of the planets are obtained, and no obvious periodic TTV signal is found in both systems after an analysis. The maximum mass of a perturbing planet located at the 1:2 mean motion resonance (MMR) for WASP-43 b and TrES-3 b is estimated to be 1.826 and 1.504 Earth mass, respectively. By quadratic fitting, it is confirmed that WASP-43 b may have a long-term TTV which means an orbital decay. The decay rate is shown to be $\dot{P} =(-0.005248 \pm 0.001714)$~s$\cdot$yr$^{-1}$, and compared with the previous results. Based on this, the lower limit of the stellar tidal quality parameter of WASP-43 is calculated to be $Q'_{\rm *} \ge 1.5\times10^5$, and the remaining lifetimes of the planets are presented for the different $Q'_{\rm *}$ values of the two systems, correspondingly.

Recently, low frequency, broadband radio emission has been observed accompanying bright meteors by the Long Wavelength Array (LWA). The broadband spectra between 20 and 60 MHz were captured for several events, while the spectral index (dependence of flux density on frequency, with $S_\nu \propto \nu^\alpha$) was estimated to be $-4\pm1$ during the peak of meteor afterglows. Here we present a survey of meteor emission and other transient events using the Murchison Widefield Array (MWA) at 72-103 MHz. In our 322-hour survey, down to a $5\sigma$ detection threshold of 3.5 Jy/beam, no transient candidates were identified as intrinsic emission from meteors. We derived an upper limit of -3.7 (95% confidence limit) on the spectral index in our frequency range. We also report detections of other transient events, like reflected FM broadcast signals from small satellites, conclusively demonstrating the ability of the MWA to detect and track space debris on scales as small as 0.1 m in low Earth orbits.

The optical images of near-Earth asteroid 4179 Toutatis acquired by Chang'e-2 spacecraft show that Toutatis has an elongated contact binary configuration, with the contact point located along the long axis. We speculate that such configuration may have resulted from a low-speed impact between two components. In this work, we performed a series of numerical simulations and compared the results with the optical images, to examine the mechanism and better understand the formation of Toutatis. Herein we propose an scenario that an assumed separated binary precursor could undergo a close encounter with Earth, leading to an impact between the primary and secondary, and the elongation is caused by Earth's tide. The precursor is assumed to be a doubly synchronous binary with a semi-major axis of 4 $R_p$ (radius of primary) and the two components are represented as spherical cohesionless self-gravitating granular aggregates. The mutual orbits are simulated in a Monte Carlo routine to provide appropriate parameters for our \emph{N}-body simulations of impact and tidal distortion. We employ the \emph{pkdgrav} package with a soft-sphere discrete element method (SSDEM) to explore the entire scenarios. The results show that contact binary configurations are natural outcomes under this scenario, whereas the shape of the primary is almost not affected by the impact of the secondary. However, our simulations further provide an elongated contact binary configuration best-matching to the shape of Toutatis at an approaching distance $r_p$ = 1.4 $\sim$ 1.5 $R_e$ (Earth radius), indicative of a likely formation scenario for configurations of Toutatis-like elongated contact binaries.

The Kepler mission has released over 4496 planetary candidates, among which 3483 planets have been confirmed as of April 2017. The statistical results of the planets show that there are two peaks around 1.5 and 2.0 in the distribution of orbital period ratios. The observations indicate that a plenty of planet pairs could have firstly been captured into mean motion resonances (MMRs) in planetary formation. Subsequently, these planets depart from exact resonant locations to be near MMRs configurations. Through type I migration, two low-mass planets have a tendency to be trapped into first-order MMRs (2:1 or 3:2 MMRs), however two scenarios of mass accretion of planets and potential outward migration play an important role in reshaping their final orbital configurations. Under the scenario of mass accretion, the planet pairs can cross 2:1 MMRs and then enter into 3:2 MMRs, especially for the inner pairs. With such formation scenario, the possibility that two planets are locked into 3:2 MMRs can increase if they are formed in a flat disk. Moreover, the outward migration can make planets have a high likelihood to be trapped into 3:2 MMRs. We perform additional runs to investigate the mass relationship for those planets in three-planet systems, and we show that two peaks near 1.5 and 2.0 for the period ratios of two planets can be easily reproduced through our formation scenario. We further show that the systems in chain resonances (e.g., 4:2:1, 3:2:1, 6:3:2 and 9:6:4 MMRs), have been observed in our simulations. This mechanism can be applicable to understand the formation of systems of Kepler-48, Kepler-53, Kepler-100, Kepler-192, Kepler-297, Kepler-399, and Kepler-450.

Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that the scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planetary configuration. These studies indicate that circumbinary planets discovered by Kepler may have experienced scattering process. In this work, we extensively investigate the scattering outcomes of circumbinary planets focusing on the effects of planet mass ratio. We find that the planetary mass ratio and the the initial relative locations of planets act as two important parameters which affect the eccentricity distribution of the surviving planets. As an application of our model, we discuss the observed orbital configurations of Kepler-34b and Kepler-413b. We first adopt the results from the disk-driven models as the initial conditions, then simulate the scattering process occurred in the late evolution stage of circumbinary planets. We show that the present orbital configurations of Kepler-34b and Kepler-413b can be well reproduced when considering two unequal-mass planet ejection model. Our work further suggests that some of the currently discovered circumbinary single-planet systems may be the survivals of original multiple-planet systems. The disk-driven migration and the scattering events occurring in the late stage both play an irreplaceable role in sculpting the final systems.

In asteroid rendezvous missions, the dynamical environment near the asteroid's surface should be made clear prior to the mission launch. However, most of the asteroids have irregular shapes, which lower the efficiency of calculating their gravitational field by adopting the traditional polyhedral method. In this work, we propose a method to partition the space near the asteroid adaptively along three spherical coordinates and use Chebyshev polynomials interpolation to represent the gravitational acceleration in each cell. Moreover, we compare four different interpolation schemes to obtain the best precision with the identical initial parameters. An error-adaptive octree division is combined to improve the interpolation precision near the surface. As an example, we take the typical irregular-shaped near-Earth asteroid 4179 Toutatis to show the advantage of this method, as a result, we show that the efficiency can be increased by hundreds to thousands times with our method. In a word, this method can be applicable to other irregular-shaped asteroids and can greatly improve the evaluation efficiency.

(349) Dembowska, a large, bright main-belt asteroid, has a fast rotation and oblique spin axis. It may have experienced partial melting and differentiation. We constrain Dembowska's thermophysical properties, e.g., thermal inertia, roughness fraction, geometric albedo and effective diameter within 3$\sigma$ uncertainty of $\Gamma=20^{+12}_{-7}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.25^{+0.60}_{-0.25}$, $p_{\rm v}=0.309^{+0.026}_{-0.038}$, and $D_{\rm eff}=155.8^{+7.5}_{-6.2}\rm~km$, by utilizing the Advanced Thermophysical Model (ATPM) to analyse four sets of thermal infrared data obtained by IRAS, AKARI, WISE and Subaru/COMICS at different epochs. In addition, by modeling the thermal lightcurve observed by WISE, we obtain the rotational phases of each dataset. These rotationally resolved data do not reveal significant variations of thermal inertia and roughness across the surface, indicating the surface of Dembowska should be covered by a dusty regolith layer with few rocks or boulders. Besides, the low thermal inertia of Dembowska show no significant difference with other asteroids larger than 100 km, indicating the dynamical lives of these large asteroids are long enough to make the surface to have sufficiently low thermal inertia. Furthermore, based on the derived surface thermophysical properties, as well as the known orbital and rotational parameters, we can simulate Dembowska's surface and subsurface temperature throughout its orbital period. The surface temperature varies from $\sim40$ K to $\sim220$ K, showing significant seasonal variation, whereas the subsurface temperature achieves equilibrium temperature about $120\sim160$ K below $30\sim50$ cm depth.

In this work, we investigate the surface thermophysical properties (thermal emissivity, thermal inertia, roughness fraction and geometric albedo) of asteroid (99942) Apophis, using the currently available thermal infrared observations of CanariCam on Gran Telescopio CANARIAS and far-infrared data by PACS of Herschel, on the basis of the Advanced thermophysical model. We show that the thermal emissivity of Apophis should be wavelength dependent from $8.70~\mu m$ to $160~\mu m$, and the maximum emissivity may arise around $20~\mu m$ similar to that of Vesta. Moreover, we further derive the thermal inertia, roughness fraction, geometric albedo and effective diameter of Apophis within a possible 1$\sigma$ scale of $\Gamma=100^{+100}_{-52}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.78\sim1.0$, $p_{\rm v}=0.286^{+0.030}_{-0.026}$, $D_{\rm eff}=378^{+19}_{-25}\rm~m$, and 3$\sigma$ scale of $\Gamma=100^{+240}_{-100}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.2\sim1.0$, $p_{\rm v}=0.286^{+0.039}_{-0.029}$, $D_{\rm eff}=378^{+27}_{-29}\rm~m$. The derived low thermal inertia but high roughness fraction may imply that Apophis could have regolith on its surface, and less regolith migration process has happened in comparison with asteroid Itokawa. Our results show that small-size asteroids could also have fine regolith on the surface, and further infer that Apophis may be delivered from the Main Belt by Yarkovsky effect.

In this work, we extensively investigate the formation of near 4:2:1 mean motion resonances (MMRs) configuration by performing two sets of N-body simulations. We model the eccentricity damping, gas drag, type I and type II planetary migration of planetesimals, planetary embryos and giant planets in the first sets. For the simulations of giant planets with type II migration, the massive terrestrial planets, with a mass up to several Earth masses, are likely produced in the systems. We further show that by shepherding and/or scattering mechanisms through Jovian planet's type II migration, the terrestrial planets and giant planets in the systems can be evolved into a near 4:2:1 MMRs. Moreover, the models are applicable to the formation of Kepler-238 and 302 systems. In the second set, we study the 4:2:1 MMRs formation in the terrestrial planetary systems, where the planets undergo type I migration and eccentricity damping. By considering type I migration, $\sim$ 17.1\% of the simulations indicate that terrestrial planets are evolved into 4:2:1 MMRs. However, this probability should depend on the initial conditions of planets. Hence, we conclude that both type I and type II migration can play a crucial role in close-in terrestrial planet formation.

We use extensive global two-dimensional hydrodynamic disk gas+dust simulations with embedded planets, coupled with three dimensional radiative transfer calculations, to model the dust ring and gap structures in the HL Tau protoplanetary disk observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). We include the self-gravity of disk gas and dust components and make reasonable choices of disk parameters, assuming an already settled dust distribution and no planet migration. We can obtain quite adequate fits to the observed dust emission using three planets with masses 0.35, 0.17, and 0.26 $M_{Jup}$ at 13.1, 33.0, and 68.6 AU, respectively. Implications for the planet formation as well as the limitations of this scenario are discussed.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

This paper presents analysis of the rotational parameters of Toutatis based on the observational results from Chang'e-2's close flyby. The 3-D shape model derived from ground-based radar observation is used to calculate the 3-1-3 Euler angles at the flyby epoch, which are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$. The large amplitude of Toutatis' tumbling attitude is demonstrated to be the result of the large deviation of the angular momentum axis and the rotational axis. Two rotational periods are evaluated to be $5.38\pm0.03$ days for rotation about the long axis and $7.40\pm0.03$ days for precession of the long axis about the angular momentum vector based on Fourier analysis. These results provide a further understanding of rotational state of Toutatis.

In this work, we investigate the thermophysical properties, including thermal inertia, roughness fraction and surface grain size of OSIRIS-REx target asteroid (101955) Bennu by using a thermophysical model with the recently updated 3D radar-derived shape model (\cite[Nolan et al., 2013]{Nolan2013}) and mid-infrared observations (\cite[M$\ddot{u}$ller et al, 2012]{Muller2012}, \cite[Emery et al., 2014]{Emery2014}). We find that the asteroid bears an effective diameter of $510^{+6}_{-40}$ m, a geometric albedo of $0.047^{+0.0083}_{-0.0011}$, a roughness fraction of $0.04^{+0.26}_{-0.04}$, and thermal inertia of $240^{+440}_{-60}\rm~Jm^{-2}s^{-0.5}K^{-1}$ for our best-fit solution. The best-estimate thermal inertia suggests that fine-grained regolith may cover a large portion of Bennu's surface, where a grain size may vary from $1.3$ to $31$~mm. Our outcome suggests that Bennu is suitable for the OSIRIS-REx mission to return samples to Earth.

On 13 December 2012, Chang'e-2 completed a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 meters from the asteroid's surface. The observations show that Toutatis has an irregular surface and its shape resembles a ginger-root of a smaller lobe (head) and a larger lobe (body). Such bilobate shape is indicative of a contact binary origin for Toutatis. In addition, the high-resolution images better than 3 meters provide a number of new discoveries about this asteroid, such as an 800-meter depression at the end of the large lobe, a sharply perpendicular silhouette near the neck region, boulders, indicating that Toutatis is probably a rubble-pile asteroid. Chang'e-2 observations have significantly revealed new insights into the geological features and the formation and evolution of this asteroid. In final, we brief the future Chinese asteroid mission concept.

Boulders are ubiquitously found on the surfaces of small rocky bodies in the inner solar system and their spatial and size distributions give insight into the geological evolution and collisional history of the parent bodies. Using images acquired by the Chang'e-2 spacecraft, more than 200 boulders have been identified over the imaged area of the near-Earth asteroid Toutatis. The cumulative boulder size frequency distribution (SFD) shows a steep slope of -4.4 $\pm$ 0.1, which is indicative of a high degree of fragmentation. Similar to Itokawa, Toutatis probably has a rubble-pile structure, as most boulders on its surface cannot solely be explained by impact cratering. The significantly steeper slope for Toutatis' boulder SFD compared to Itokawa may imply a different preservation state or diverse formation scenarios. In addition, the cumulative crater SFD has been used to estimate a surface crater retention age of approximately 1.6 $\pm$ 0.3 Gyr.

The X-ray Timing and Polarization (XTP) is a mission concept for a future space borne X-ray observatory and is currently selected for early phase study. We present a new design of X-ray polarimeter based on the time projection gas chamber. The polarimeter, placed above the focal plane, has an additional rear window that allows hard X-rays to penetrate (a transmission of nearly 80% at 6 keV) through it and reach the detector on the focal plane. Such a design is to compensate the low detection efficiency of gas detectors, at a low cost of sensitivity, and can maximize the science return of multilayer hard X-ray telescopes without the risk of moving focal plane instruments. The sensitivity in terms of minimum detectable polarization, based on current instrument configuration, is expected to be 3% for a 1mCrab source given an observing time of 10^5 s. We present preliminary test results, including photoelectron tracks and modulation curves, using a test chamber and polarized X-ray sources in the lab.

In this work, we investigate the thermophysical properties of OSIRIS-REx target asteroid (101955) Bennu (hereafter, Bennu), where thermal inertia plays an important role in understanding the nature of the asteroid's surface, and will definitely provide substantial information for the sampling return mission. Using a thermophysical model incorporating the recently updated 3D radar-derived shape model \citep{Nolan2013} and mid-infrared observations of Spitzer-PUI, Spitzer-IRAC, Herschel/PACS and ESO VLT/VISIR \citep{Muller2012,Emery2014}, we derive the surface thermophysical properties of Bennu. The asteroid has an effective diameter of $510^{+6}_{-40}$ m, a geometry albedo of $0.047^{+0.0083}_{-0.0011}$, a roughness fraction of $0.04^{+0.26}_{-0.04}$, and thermal inertia of $240^{+440}_{-60}\rm~Jm^{-2}s^{-0.5}K^{-1}$ for a best-fit solution at 1$\sigma$ level. The best-estimate thermal inertia indicates that fine-grained regolith may cover a large area of Bennu's surface, with a grain size that may range from $1.3$ to $31$~mm, and our outcome further supports that Bennu would be a suitable target for the OSIRIS-REx mission to return samples from the asteroid to Earth.

Binary stars are predicted to have an important role in the evolution of globular clusters, so we obtained binary fractions for 35 globular clusters that were imaged in the F606W and F814W with the ACS on the Hubble Space Telescope. When compared to the values of prior efforts (Sollima et al. 2007; Milone et al. 2012), we find significant discrepancies, despite each group correcting for contamination effects and having performed the appropriate reliability tests. The most reliable binary fractions are obtained when restricting the binary fraction to q > 0.5. Our analysis indicates that the range of the binary fractions is nearly an order of magnitude for the lowest dynamical ages, suggesting that there is a broad distribution in the binary fraction at globular cluster formation. Dynamical effects also appears to decrease the core binary fractions by a factor of two over a Hubble time, but this is a weak relationship. We confirm a correlation from previous work that the binary fraction within the core radius decreases with cluster age, indicating that younger clusters formed with higher binary fractions. The strong radial gradient in the binary fraction with cluster radius appears to be a consequence of dynamical interactions. It is often not present in dynamically young clusters but nearly always present in dynamically old clusters.

In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth asteroid 4179 Toutatis, which was closely observed by Chang'e-2 at a distance of $770\pm120~$ meters from the asteroid's surface during the outbound flyby \citep{Huang2013} on 13 December 2012. A sequence of high-resolution images was acquired during the flyby mission. In combination with ground-based radar observations collected over the last two decades, we analyze these flyby images and determine the orientation of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be $-20.1^\circ\pm1^\circ$, $27.6^\circ\pm1^\circ$ and $42.2^\circ\pm1^\circ$, respectively. The least-squares method is utilized to determine the rotational parameters and spin state of Toutatis. The characteristics of the spin-state parameters and angular momentum variations are extensively studied using numerical simulations, which confirm those reported by \citet{Takahashi2013}. The large amplitude of Toutatis' precession is assumed to be responsible for its tumbling attitude as observed from Earth. Toutatis' angular momentum orientation is determined to be described by $\lambda_{H}=180.2^{+0.2^\circ}_{-0.3^\circ}$ and $\beta_{H}=-54.75^{+0.15^\circ}_{-0.10^\circ}$, implying that it has remained nearly unchanged for two decades. Furthermore, using Fourier analysis to explore the change in the orientation of Toutatis' axes, we reveal that the two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results of the former investigation. Hence, our investigation provides a clear understanding of the state of the rotational dynamics of Toutatis.

The Kepler mission has released ~4229 transiting planet candidates. There are approximately 222 candidate systems with three planets. Among them, the period ratios of planet pairs near 1.5 and 2.0 reveal that two peaks exist for which the proportions of the candidate systems are ~7.0% and 18.0%, respectively. In this work, we study the formation of mean motion resonance (MMR) systems, particularly for the planetary configurations near 3:2 and 2:1 MMRs, and we concentrate on the interplay between the resonant configuration and the combination of stellar accretion rate, stellar magnetic field, speed of migration and additional planets. We perform more than 1000 runs by assuming a system with a solar-like star and three surrounding planets. From the statistical results, we find that under the formation scenario, the proportions near 1.5 and 2.0 can reach 14.5% and 26.0%, respectively. In addition, $\dot M=0.1\times 10^{-8} ~M_\odot ~{\rm yr^{-1}}$ is propitious toward the formation of 3:2 resonance, whereas $\dot M=2\times 10^{-8} ~M_\odot ~{\rm yr^{-1}}$ contributes to the formation of 2:1 resonance. The speed-reduction factor of type I migration $f_1\geq 0.3$ facilitates 3:2 MMRs, whereas $f_1\geq 0.1$ facilitates 2:1 MMRs. If additional planets are present in orbits within the innermost or beyond the outermost planet in a three-planet system, 3:2:1 MMRs can be formed, but the original systems trapped in 4:2:1 MMRs are not affected by the supposed planets. In summary, we conclude that this formation scenario will provide a likely explanation for Kepler candidates involved in 2:1 and 3:2 MMRs.

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.

In this paper, a 3D convex shape model of (175706) 1996 FG3, which consists of 2040 triangle facets and 1022 vertices, is derived from the known lightcurves. The best-fit orientation of the asteroid's spin axis is determined to be $\lambda =237.7^\circ$ and $\beta=-83.8^{\circ}$ considering the observation uncertainties, and its rotation period is $\sim$ 3.5935 h . Using the derived shape model, we adopt the so-called advanced thermophysical model (ATPM) to fit three published sets of mid-infrared observations of 1996 FG3 \citep{Wolters2011,Walsh2012}, so as to evaluate its surface properties. Assuming the primary and the secondary bear identical shape, albedo, thermal inertia and surface roughness, the best-fit parameters are obtained from the observations. The geometric albedo and effective diameter of the asteroid are reckoned to be $p_{\rm v}=0.045\pm0.002$, $D_{\rm eff}=1.69^{+0.05}_{-0.02}$ km. The diameters of the primary and secondary are determined to be $D_{1}=1.63^{+0.04}_{-0.03}$ km and $D_{2}=0.45^{+0.04}_{-0.03}$ km, respectively. The surface thermal inertia $\Gamma$ is derived to be a low value of $80\pm40\rm Jm^{-2}s^{-0.5}K^{-1}$ with a roughness fraction $f_{\rm R}$ of $0.8^{+0.2}_{-0.4}$. This indicates that the primary possibly has a regolith layer on its surface, which is likely to be covered by a mixture of dust, fragmentary rocky debris and sand. The minimum regolith depth is estimated to be $5\sim20\rm mm$ from the simulations of subsurface temperature distribution, indicating that 1996 FG3 could be a very suitable target for a sample return mission.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

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.

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.

Binary stars are thought to be a controlling factor in globular cluster evolution, since they can heat the environmental stars by converting their binding energy to kinetic energy during dynamical interactions. Through such interaction, the binaries determine the time until core collapse. To test predictions of this model, we have determined binary fractions for 35 clusters. Here we present our methodology with a representative globular cluster NGC 4590. We use HST archival ACS data in the F606W and F814W bands and apply PSF-fitting photometry to obtain high quality color-magnitude diagrams. We formulate the star superposition effect as a Poisson probability distribution function, with parameters optimized through Monte-Carlo simulations. A model-independent binary fraction of (6.2 +- 0.3)% is obtained by counting stars that extend to the red side of the residual color distribution after accounting for the photometric errors and the star superposition effect. A model-dependent binary fraction is obtained by constructing models with a known binary fraction and an assumed binary mass-ratio distribution function. This leads to a binary fraction range of 6.8% to 10.8%, depending on the assumed shape to the binary mass ratio distribution, with the best fit occurring for a binary distribution that favors low mass ratios (and higher binary fractions). We also represent the method for radial analysis of the binary fraction in the representative case of NGC 6981, which shows a decreasing trend for the binary fraction towards the outside, consistent with theoretical predictions for the dynamical effect on the binary fraction.

The recent Kepler discovery of KOI-152 reveals a system of three hot super-Earth candidates that are in or near a 4:2:1 mean motion resonance. It is unlikely that they formed in situ, the planets probably underwent orbital migration during the formation and evolution process. The small semimajor axes of the three planets suggest that migration stopped at the inner edge of the primordial gas disk. In this paper we focus on the influence of migration halting mechanisms, including migration "dead zones", and inner truncation by the stellar magnetic field. We show that the stellar accretion rate, stellar magnetic field and the speed of migration in the proto-planetary disk are the main factors affecting the final configuration of KOI-152. Our simulations suggest that three planets may be around a star with low star accretion rate or with high magnetic field. On the other hand, slow type I migration, which decreases to one tenth of the linear analysis results, favors forming the configuration of KOI-152. Under such formation scenario, the planets in the system are not massive enough to open gaps in the gas disk. The upper limit of the planetary masses are estimated to be about $15,~19$, and $24 M_\oplus$, respectively. Our results are also indicative of the near Laplacian configurations that are quite common in planetary systems.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

In terms of energy resolution, temporal response to burst events, and thermal stability, lanthanum bromide doped with Ce is a much better choice than the traditional NaI(Tl) scintillator for hard X-ray astronomy. We present the test results of a phoswich detector with a diameter of 101.6 mm consisting of 6 mm thick LaBr3:Ce and 40 mm thick NaI(Tl), which is the largest one of this type reported so far. The measured energy resolution is 10.6% at 60 keV, varying inversely proportional to the square root of the energy, and the energy nonlinearity is found to be less than 1%, as good as those of smaller phoswiches. The coupled scintillators and phototube also show excellent uniformity across the detecting surface, with a deviation of 0.7% on the pulse amplitude produced by 60 keV gamma-rays. Thanks to the large ratio of light decay times of NaI(Tl) and LaBr3:Ce, 250 ns vs. 16 ns, pulse shape discrimination is much easier for this combination than for NaI(Tl)/CsI(Na). As the light decay time of LaBr3:Ce is about 15 times faster than that of NaI(Tl), this phoswich is more suitable for detection of bright, transient sources such as gamma-ray bursts and soft gamma-ray repeaters. The internal activity of lanthanum produces a count rate of about 6 counts/s at 37.5 keV in the detector. This peak could be used for in-flight spectral calibration and gain correction.

In this work, we extensively investigate the terrestrial planetary formation for the inclined planetary systems (considering the OGLE-2006-BLG-109L system as prototype) in the late stage. In the simulations, we show that the occurrence of terrestrial planets is quite common, in the final assembly stage. Moreover, we find that 40% of the runs finally occupy one planet in the habitable zone (HZ). On the other hand, the numerical results also indicate that the inner region of the planetesimal disk, ranging from $\sim 0.1$ to 0.3 AU, plays an important role in building up terrestrial planets. By examining all simulations, we note that the survivals are located either between 0.1$\sim$1.0 AU or beyond 7 AU, or at the 1:1 mean motion resonance of OGLE-2006-BLG-109Lb at $\sim$2.20 AU. The outcomes suggest that it may exist moderate possibility for the inclined systems to harbor terrestrial planets, even planets in the HZs.

The large number of exoplanets found to orbit their host stars in very close orbits have significantly advanced our understanding of the planetary formation process. It is now widely accepted that such short-period planets cannot have formed {\em in situ}, but rather must have migrated to their current orbits from a formation location much farther from their host star. In the late stages of planetary formation, once the gas in the proto-planetary disk has dissipated and migration has halted, gas-giants orbiting in the inner disk regions will excite planetesimals and planetary embryos, resulting in an increased rate of orbital crossings and large impacts. We present the results of dynamical simulations for planetesimal evolution in this later stage of planet formation. We find that a mechanism is revealed by which the collision-merger of planetary embryos can kick terrestrial planets directly into orbits extremely close to their parent stars.

The source of hot gas in elliptical galaxies is thought to be due to stellar mass loss, with contributions from supernova events and possibly from infall from a surrounding environment. This picture predicts supersolar values for the metallicity of the gas toward the inner part of the galaxy, which can be tested by measuring the gas phase abundances. We use high-quality data for 10 nearby early-type galaxy from XMM-Newton, featuring both the EPIC and the Reflection Grating Spectrometer, where the strongest emission lines are detected with little blending; some Chandra data are also used. We find excellent consistency in the elemental abundances between the different XMM instruments and good consistency with Chandra. Differences in abundances with aperture size and model complexity are examined, but large differences rarely occur. For a two-temperature thermal model plus a point source contribution, the median Fe and O abundances are 0.86 and 0.44 of the Solar value, while Si and Mg abundances are similar to that for Fe. This is similar to stellar abundances for these galaxies but supernovae were expected to enhance the gas phase abundances considerably, which is not observed.

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) through a Baroclinic Instability and subsequent development of Rossby-Wave-Instability (RWI). It is found that the generation of potential vorticity is very common and effective in non-barotropic disks through the Baroclinic Instability, especially within the coorbital region. Our results also complement previous studies by \citet{KLL03} that non-axisymmetric RWIs are likely to develop at local minima of potential vorticity distribution that are generated by the interaction between a planet and a inviscid barotropic disk. This second instability appears to be very common and robust, regardless of the equation of state, initial density distribution, and rotational law of the disk. The development of RWIs results in non-axisymmetric density blobs, which exert stronger torques onto the planet when they travel in the vicinity of the planet. As a result of that, large amplitude oscillations are introduced to the time behavior of the total torque acted on the planet by the disk. In our current simulations, RWIs do not change the overall picture of inward orbital migration but bring in 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.

(Abridged)We numerically investigated the dynamical architecture of 47 UMa with the planetary configuration of the best-fit orbital solutions by Fischer et al. We systematically studied the existence of Earth-like planets in the region 0.05 AU $\leq a \leq 2.0$ AU for 47 UMa with numerical simulations, and we also explored the packed planetary geometry and Trojan planets in the system. In the simulations, we found that "hot Earths" at 0.05 AU $\leq a < $ 0.4 AU can dynamically survive at least for 1 Myr. The Earth-like planets can eventually remain in the system for 10 Myr in areas involved in the mean motion resonances (MMR) (e.g., 3:2 MMR) with the inner companion. Moreover, we showed that the 2:1 and 3:1 resonances are on the fringe of stability, while the 5:2 MMR is unstable. Additionally, the 2:1 MMR marks out a remarkable boundary between chaotic and regular motions, inside, most of the orbits can survive, outside, they are mostly lost in the orbital evolution. In a dynamical sense, the most likely candidate for habitable environment is Earth-like planets with orbits in the ranges 0.8 AU $\leq a < 1.0$ AU and 1.0 AU $ < a < 1.30$ AU (except several unstable cases) with relatively low eccentricities. The Trojan planets with low eccentricities and inclinations can secularly last at the triangular equilibrium points of two massive planets. Hence, the 47 UMa planetary system may be a close analog to our solar system.

We investigated the apsidal motion for the multi-planet systems. In the simulations, we found that the two planets of HD 37124, HD 12661, 47 Uma and HD 82943 separately undergo apsidal alignment or antialignment. But the companions of GJ 876 and $\upsilon$ And are only in apsidal lock about $0^{\circ}$. Moreover, we obtained the criteria with Laplace-Lagrange secular theory to discern whether a pair of planets for a certain system are in libration or circulation.

To date, more than 100 giant Jupiter-like planets have been discovered in Doppler surveys of solar-type stars. In this paper, we perform simulations to investigate three systems: GJ 876, HD 82943 and 55 Cnc. The former two systems both have a pair of planets in the 2:1 Mean Motion Resonance (MMR), while the inner two companions of the later is close to 3:1 MMR. By integrating hundreds of the planetary orbits of three systems for million years, we find that for GJ 876 and HD 82943, the critical argument $\lambda_{1} - 2\lambda_{2} + \varpi_{1}$ and $\lambda_{1} - 2\lambda_{2} + \varpi_{2}$ librate about $0^{\circ}$ or $180^{\circ}$, indicating 2:1 MMR can play an important role in stabilizing the motion of the planets so that they are protected from frequent close encounters. As for 55 Cnc, we further show the three resonant arguments for 3:1 MMR execute librations for millions of years respectively, which reveals the evidence of the resonance for this system. Additionally, we should emphasize another vital mechanism is the apsidal phase-locking between a couple of planets for a certain system. For GJ 876 and HD 82943, we discover the relative apsidal longitudes $\varpi_{1} - \varpi_{2}$ move about $0^{\circ}$ or $180^{\circ}$, respectively; but for 55 Cnc, we find that there exists an asymmetric apsidal libration between two inner planets. Finally, we made a brief discussion about the Habitable Zones in the exoplanetary systems.

We present results from the analysis of X-ray power density spectra and coherence when GRS 1915+105 is in soft states. We use three data sets that belong to $\mu$, $\phi$ and $\delta$ classes in Belloni et al (2000). We find that the power density spectra appear to be complex, with several features between 0.01 and 10 Hz. The coherence deviates from unity above a characteristic frequency. We discuss our results in different models. The corona size in the sphere-disk model implied by this break frequency is on the order of 10$^4$ GM/c$^2$, which is unphysical. Our results are more consistent with the prediction of the model of a planar corona sustained by magnetic flares, in which the characteristic frequency is associated with the longest time-scale of an individual flare, which is about eight seconds.