We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio

Since September 2018, LAMOST starts a new 5-year medium-resolution spectroscopic survey (MRS) using bright/gray nights. We present the scientific goals of LAMOST-MRS and propose a near optimistic strategy of the survey. A complete footprint is also provided. Not only the regular medium-resolution survey, but also a time-domain spectroscopic survey is being conducted since 2018 and will be end in 2023. According to the detailed survey plan, we expect that LAMOST-MRS can observe about 2 million stellar spectra with ~7500 and limiting magnitude of around G=15 mag. Moreover, it will also provide about 200 thousand stars with averagely 60-epoch observations and limiting magnitude of G~14 mag. These high quality spectra will give around 20 elemental abundances, rotational velocities, emission line profiles as well as precise radial velocity with uncertainty less than 1 km/s. With these data, we expect that LAMOST can effectively leverage sciences on stellar physics, e.g. exotic binary stars, detailed observation of many types of variable stars etc., planet host stars, emission nebulae, open clusters, young pre-main-sequence stars etc.

The Kepler mission has discovered thousands of exoplanets around various stars with different spectral types (M, K, G, and F) and thus different masses and effective temperatures. Previous studies have shown that the planet occurrence rate, in terms of average number of planets per star, drops with increasing stellar effective temperature (Teff). In this paper, with the final Kepler Data Release (DR25) catalog, we revisit the relation between stellar effective temperature (as well as mass) and planet occurrence, but in terms of the fraction of stars with planets and the number of planets per planetary system (i.e. planet multiplicity). We find that both the fraction of stars with planets and planet multiplicity decrease with increasing stellar temperature and mass. Specifically, about 75% late-type stars (Teff<5000 K) have Kepler-like planets with an average planet multiplicity of ~2.8; while for early-type stars (Teff>6500 K), this fraction and the average multiplicity fall down to ~35% and ~1.8, respectively. The decreasing trend in the fraction of stars with planets is very significant with ${\Delta}$AIC> 30, though the trend in planet multiplicity is somewhat tentative with ${\Delta}$AIC~5. Our results also allow us to derive the dispersion of planetary orbital inclinations in relationship with stellar effective temperature. Interestingly, it is found to be similar to the well-known trend between obliquity and stellar temperature, indicating that the two trends might have a common origin.

The radial velocity (RV) is a basic physical quantity which can be determined through Doppler shift of the spectrum of a star. The precision of RV measurement depends on the resolution of the spectrum we used and the accuracy of wavelength calibration. In this work, radial velocities of LAMOST-II medium resolution (R ~ 7500) spectra are measured for 1,594,956 spectra (each spectrum has two wavebands) through matching with templates. A set of RV standard stars are used to recalibrate the zero point of the measurement, and some reference sets with RVs derived from medium/high-resolution observations are used to evaluate the accuracy of the measurement. Comparing with reference sets, the accuracy of our measurement can get 0.0227 km s/1 with respect to radial velocities standard stars. The intrinsic precision is estimated with the multiple observations of single stars, which can achieve to 1.36 km s/1,1.08 km s/1, 0.91 km s/1 for the spectra at signal-to-noise levels of 10, 20, 50, respectively.

There is no consensus on the progenitors of Type Ia supernovae (SNe Ia) despite their importance for cosmology and chemical evolution. We address this question by using our previously published catalogs of Mg, Si, Ca, Cr, Fe, Co, and Ni abundances in dwarf galaxy satellites of the Milky Way to constrain the mass at which the white dwarf explodes during a typical SN Ia. We fit a simple bi-linear model to the evolution of [X/Fe] with [Fe/H], where X represents each of the elements mentioned above. We use the evolution of [Mg/Fe] coupled with theoretical supernova yields to isolate what fraction of the elements originated in SNe Ia. Then, we infer the [X/Fe] yield of SNe Ia for all of the elements except Mg. We compare these observationally inferred yields to recent theoretical predictions for two classes of Chandrasekhar-mass (M_Ch) SN Ia as well as sub-M_Ch SNe Ia. Most of the inferred SN Ia yields are consistent with all of the theoretical models, but [Ni/Fe] is consistent only with sub-M_Ch models. We conclude that the dominant type of SN Ia in ancient dwarf galaxies is the explosion of a sub-M_Ch white dwarf. The Milky Way and dwarf galaxies with extended star formation histories have higher [Ni/Fe] abundances, which could indicate that the dominant class of SN Ia is different for galaxies where star formation lasted for at least several Gyr.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

We present measurements of the abundances of chromium, cobalt, and nickel in 4113 red giants, including 2277 stars in globular clusters, 1820 stars in the Milky Way's dwarf satellite galaxies, and 16 field stars. We measured the abundances from mostly archival Keck/DEIMOS medium-resolution spectroscopy with a resolving power of R ~ 6500 and a wavelength range of approximately 6500-9000 A. The abundances were determined by fitting spectral regions that contain absorption lines of the elements under consideration. We used estimates of temperature, surface gravity, and metallicity that we previously determined from the same spectra. We estimated systematic error by examining the dispersion of abundances within mono-metallic globular clusters. The median uncertainties for [Cr/Fe], [Co/Fe], and [Ni/Fe] are 0.20, 0.20, and 0.13, respectively. Finally, we validated our estimations of uncertainty through duplicate measurements, and we evaluated the accuracy and precision of our measurements through comparison to high-resolution spectroscopic measurements of the same stars.

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.

The LAMOST-Kepler (LK-) project was initiated to use the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) to make spectroscopic follow-up observations for the targets in the field of the Kepler mission. The Kepler field is divided into 14 subfields that are adapted to the LAMOST circular field with a diameter of 5 degrees. During the regular survey phase of LAMOST, the LK-project took data from 2012 June to 2017 June and covered all the 14 subfields at least twice. In particular, we describe in this paper the second Data Release of the LK-project, including all spectra acquired through 2015 May to 2017 June together with the first round observations of the LK-project from 2012 June to 2014 September. The LK-project now counts 227 870 spectra of 156 390 stars, among which we have derived atmospheric parameters (log g, T eff and [Fe/H]) and heliocentric radial velocity (RV) for 173 971 spectra of 126 172 stars. These parameters were obtained with the most recent version of the LAMOST Stellar Parameter Pipeline v 2.9.7. Nearly one half, namely 76 283 targets, are observed both by LAMOST and Kepler telescopes. These spectra, establishing a large spectroscopy library, will be useful for the entire astronomical community, particularly for planetary science and stellar variability on Kepler targets.

During the evolution of proto-planetary disc, photo-evaporations of both central and external stars play important roles. Considering the complicated radiation surroundings in the clusters, where the star formed, the proto-planetary discs survive in different lifetimes due to flyby events. In this paper, we mainly focus on the disc around a T Tauri star, which encounters with another main-sequence star with different temperatures in hyperbolic orbits with different peri-center distances, eccentricities and inclinations. We find the criterion for gap-opening due to photo-evaporation of central star after the flyby event. A gap is opened in the late stage of gas disc, and induce that the gap only influence the planet formation and migration limitedly. If the flyby orbit has a moderate value of peri-center distance, which weakly depends on the eccentricity and inclination, the external photo-evaporation lead to a maximum mass loss during the flyby event. Flyby stars in orbits with smaller eccentricities or larger inclinations induce larger mass loss. Adopting a simple multiple flyby models, we conclude that in open clusters, gas discs usually survive in typical lifetimes between 1 and 10 Myr, except there are many massive stars in dense open clusters. In globular clusters, discs disperse very quickly and hardly produce the gas giant planets. The fast-depleted discs are probably responsible for the null detection of giant planets in globular clusters.

Stars are commonly formed in binary systems, which provide a natural laboratory for studying planet formation in extreme conditions. In our first paper (Paper I) of a series Xie et al. (2011), we have shown that the intermediate stage - from planetesimals to planetary embryos/cores - of planet formation can proceed even in highly inclined binaries. Following Paper I, here we numerically study the late stage of terrestrial planet formation, i.e., from embryos to full planets, in binary systems of various orbital configurations. We identify an orbital alignment or anti-alignment effect; namely, although an inclined binary generally misaligns the planetary orbits with respect to the spin axis of the primary host star (i.e., causing large obliquity), it could align or anti-align the planetary orbits with respect to the binary orbit. Such an orbital (anti-)alignment effect is caused by the combination of orbital differential precession and self-damping, and it is mostly significant in cases of intermediate binary separations, i.e., a_B~40-200 AU for terrestrial planet formation around 1 AU from the primary stars. In such intermediate separation binaries, somewhat contrary to intuition, the binary companion can aid planet growth by having increased the rate of collisions, forming significantly more massive but fewer planets. In the other two ends, the companion is ether too close thus plays a violently disruptive role or too wide to have significant effect on planet formation. Future observations, which can discover more planet-bearing binary star systems and constrain their masses and 3-D orbital motions will test our numerical findings.

We constrain the intrinsic architecture of Kepler planetary systems by modeling the observed multiplicities of the transiting planets (tranets) and their transit timing variations (TTVs). We robustly determine that the fraction of Sun-like stars with Kepler-like planets, $\eta_{\rm Kepler}$, is $30\pm3\%$. Here Kepler-like planets are planets that have radii $R_{\rm p} \gtrsim R_\oplus$ and orbital periods $P<400$~days. Our result thus significantly revises previous claims that more than 50\% of Sun-like stars have such planets. Combining with the average number of Kepler planets per star ($\sim0.9$), we obtain that on average each planetary system has $3.0\pm0.3$ planets within 400 days. We also find that the dispersion in orbital inclinations of planets within a given planetary system, $\sigma_{i,k}$, is a steep function of its number of planets, $k$. This can be parameterized as $\sigma_{i,k}\propto k^\alpha$ and we find that $-4<\alpha<-2$ at 2-$\sigma$ level. Such a distribution well describes the observed multiplicities of both tranets and TTVs with no excess of single tranets. Therefore we do not find evidence supporting the so-called "Kepler dichotomy." Together with a previous study on orbital eccentricities, we now have a consistent picture: the fewer planets in a system, the hotter it is dynamically. We discuss briefly possible scenarios that lead to such a trend. Despite our Solar system not belonging to the Kepler club, it is interesting to notice that the Solar system also has three planets within 400 days and that the inclination dispersion is similar to Kepler systems of the same multiplicity.

We discover a population of short-period, Neptune-size planets sharing key similarities with hot Jupiters: both populations are preferentially hosted by metal-rich stars, and both are preferentially found in Kepler systems with single transiting planets. We use accurate LAMOST DR4 stellar parameters for main-sequence stars to study the distributions of short-period 1d < P < 10d Kepler planets as a function of host star metallicity. The radius distribution of planets around metal-rich stars is more "puffed up" as compared to that around metal-poor hosts. In two period-radius regimes, planets preferentially reside around metal-rich stars, while there are hardly any planets around metal-poor stars. One is the well-known hot Jupiters, and the other is a population of Neptune-size planets (2 R_Earth <~ R_p <~ 6 R_Earth), dubbed as "Hoptunes". Also like hot Jupiters, Hoptunes occur more frequently in systems with single transiting planets though the fraction of Hoptunes occurring in multiples is larger than that of hot Jupiters. About 1% of solar-type stars host "Hoptunes", and the frequencies of Hoptunes and hot Jupiters increase with consistent trends as a function of [Fe/H]. In the planet radius distribution, hot Jupiters and Hoptunes are separated by a "valley" at approximately Saturn size (in the range of 6 R_Earth <~ R_p <~ 10 R_Earth), and this "hot-Saturn valley" represents approximately an order-of-magnitude decrease in planet frequency compared to hot Jupiters and Hoptunes. The empirical "kinship" between Hoptunes and hot Jupiters suggests likely common processes (migration and/or formation) responsible for their existence.

The high planetary multiplicity revealed by Kepler implies that Transit Time Variations (TTVs) are intrinsically common. The usual procedure for detecting these TTVs is biased to long-period, deep transit planets whereas most transiting planets have short periods and shallow transits. Here we introduce the Spectral Approach to TTVs technique that allows expanding the TTVs catalog towards lower TTV amplitude, shorter orbital period, and shallower transit depth. In the Spectral Approach we assume that a sinusoidal TTV exists in the data and then calculate the improvement to $\chi^2$ this model allows over that of linear ephemeris model. This enables detection of TTVs even in cases where the transits are too shallow so individual transits cannot be timed. The Spectral Approach is more sensitive due to the reduced number of free parameters in its model. Using the Spectral Approach, we: (a) detect 131 new periodic TTVs in Kepler data (an increase of ~2/3 over a previous TTV catalog); (b) Constrain the TTV periods of 34 long-period TTVs and reduce amplitude errors of known TTVs; (c) Identify cases of multi-periodic TTVs, for which absolute planetary mass determination may be possible. We further extend our analysis by using perturbation theory assuming small TTV amplitude at the detection stage, which greatly speeds up our detection (to a level of few seconds per star). Our extended TTVs sample shows no deficit of short period or low amplitude transits, in contrast to previous surveys in which the detection schemes were significantly biased against such systems.

{Context}. The HIFI instrument on the Herschel Space Observatory performed over 9100 astronomical observations, almost 900 of which were calibration observations in the course of the nearly four-year Herschel mission. The data from each observation had to be converted from raw telemetry into calibrated products and were included in the Herschel Science Archive. {Aims}. The HIFI pipeline was designed to provide robust conversion from raw telemetry into calibrated data throughout all phases of the HIFI missions. Pre-launch laboratory testing was supported as were routine mission operations. {Methods}. A modular software design allowed components to be easily added, removed, amended and/or extended as the understanding of the HIFI data developed during and after mission operations. {Results}. The HIFI pipeline processed data from all HIFI observing modes within the Herschel automated processing environment as well as within an interactive environment. The same software can be used by the general astronomical community to reprocess any standard HIFI observation. The pipeline also recorded the consistency of processing results and provided automated quality reports. Many pipeline modules were in use since the HIFI pre-launch instrument level testing. {Conclusions}. Processing in steps facilitated data analysis to discover and address instrument artefacts and uncertainties. The availability of the same pipeline components from pre-launch throughout the mission made for well-understood, tested, and stable processing. A smooth transition from one phase to the next significantly enhanced processing reliability and robustness.

The nearly circular (mean eccentricity <e>~0.06) and coplanar (mean mutual inclination <i>~3 deg) orbits of the Solar System planets motivated Kant and Laplace to put forth the hypothesis that planets are formed in disks, which has developed into the widely accepted theory of planet formation. Surprisingly, the first several hundred extrasolar planets (mostly Jovian) discovered using the Radial Velocity (RV) technique are commonly on eccentric orbits (<e> ~ 0.3). This raises a fundamental question: Are the Solar System and its formation special? The Kepler mission has found thousands of transiting planets dominated by sub-Neptunes, but most of their orbital eccentricities remain unknown. By using the precise spectroscopic host star parameters from the LAMOST observations, we measure the eccentricity distributions for a large (698) and homogeneous Kepler planet sample with transit duration statistics. Nearly half of the planets are in systems with single transiting planets (singles), while the other half are multiple-transiting planets (multiples). We find an eccentricity dichotomy: on average, Kepler singles are on eccentric orbits with <e>~0.3, while the multiples are on nearly circular (<e> = 0.04^{+0.03}_{-0.04}) and coplanar (<i> = 1.4^{+0.8}_{-1.1} deg) orbits similar to the Solar System planets. Our results are consistent with previous studies of smaller samples and individual systems. We also show that Kepler multiples and solar system objects follow a common relation <e>~(1-2)x<i> between mean eccentricities and mutual inclinations. The prevalence of circular orbits and the common relation may imply that the solar system is not so atypical in the galaxy after all.

Many exoplanets have been found in orbits close to their host stars and thus they are subject to the effects of photo-evaporation. Previous studies have shown that a large portion of exoplanets detected by the Kepler mission have been significantly eroded by photo-evaporation. In this paper, we numerically study the effects of photo-evaporation on the orbital evolution of a hypothesized moon system around a planet. We find that photo-evaporation is crucial to the stability of the moon system. Photo-evaporation can erode the atmosphere of the planet thus leading to significant mass loss. As the planet loses mass, its Hill radius shrinks and its moons increase their orbital semi-major axes and eccentricities. When some moons approach their critical semi-major axes, global instability of the moon system would be triggered, which usually ends up with two, one or even zero surviving moons. Some lost moons could escape from the moon system to become a new planet orbiting the star or run away further to become a free-floating object in the Galaxy. Given the destructive role of photo-evaporation, we speculate that exo-moons are less common for close-in planets ($<0.1$ AU), especially those around M-type stars, because they are more X-ray luminous and thus enhancing photo-evaporation. The lessons we learn in this study may be helpful for the target selection of on-going/future exomoon searching programs.

Many multiple-planet systems have been found by the Kepler transit survey and various radial velocity (RV) surveys. Kepler planets show an asymmetric feature, namely, there are small but significant deficits/excesses of planet pairs with orbital period spacing slightly narrow/wide of the exact resonance, particularly near the first order mean motion resonance (MMR), such as 2:1 and 3:2 MMR. Similarly, if not exactly the same, an asymmetric feature (pileup wide of 2:1 MMR) is also seen in RV planets, but only for massive ones. We analytically and numerically study planets' orbital evolutions near and in the MMR. We find that their orbital period ratios could be asymmetrically distributed around the MMR center regardless of dissipation. In the case of no dissipation, Kepler planets' asymmetric orbital distribution could be partly reproduced for 3:2 MMR but not for 2:1 MMR, implying that dissipation might be more important to the latter. The pileup of massive RV planets just wide of 2:1 MMR is found to be consistent with the scenario that planets formed separately then migrated toward the MMR. The location of the pileup infers a K value of 1-100 on the order of magnitude for massive planets, where K is the damping rate ratio between orbital eccentricity and semimajor axis during planet migration.

The Chinese Small Telescope ARray (CSTAR) is the first Chinese astronomical instrument placed in Antarctica. It is a group of four identical, fully automatic $14.5\,\rm{cm}$ telescopes, with an field of view (FOV) of $20\,\rm{deg^2}$ centered on the South Celestial Pole. Placed at Antarctic Dome A, CSTAR is designed to provide high-cadence photometry for site monitoring and variable sources detection. During the 2008 observing season, CSTAR has taken high-precision photometric data for 18,145 stars around the South Celestial Pole. At $i\,=\,7.5$ and $12$, the photometric precision reaches $\sim 8$ mmag and $\sim 30$ mmag with a cadence of 20s or 30s, respectively. Using robust detection method, we have found 15 stellar flares on 13 sources, including two classified variables. We have also found a linear relation between the decay times and the total durations of the stellar flares. The details of all detected flares along with their stellar properties are presented in this work.

We report the detections of a giant planet (MARVELS-7b) and a brown dwarf candidate (MARVELS-7c) around the primary star in the close binary system, HD 87646. It is the first close binary system with more than one substellar circum-primary companion discovered to the best of our knowledge. The detection of this giant planet was accomplished using the first multi-object Doppler instrument (KeckET) at the Sloan Digital Sky Survey (SDSS) telescope. Subsequent radial velocity observations using ET at Kitt Peak National Observatory, HRS at HET, the "Classic" spectrograph at the Automatic Spectroscopic Telescope at Fairborn Observatory, and MARVELS from SDSS-III confirmed this giant planet discovery and revealed the existence of a long-period brown dwarf in this binary. HD 87646 is a close binary with a separation of $\sim22$ AU between the two stars, estimated using the Hipparcos catalogue and our newly acquired AO image from PALAO on the 200-inch Hale Telescope at Palomar. The primary star in the binary, HD 87646A, has Teff = 5770$\pm$80K, log(g)=4.1$\pm$0.1 and [Fe/H] = $-0.17\pm0.08$. The derived minimum masses of the two substellar companions of HD 87646A are 12.4$\pm$0.7M$_{\rm Jup}$ and 57.0$\pm3.7$M$_{\rm Jup}$. The periods are 13.481$\pm$0.001 days and 674$\pm$4 days and the measured eccentricities are 0.05$\pm$0.02 and 0.50$\pm$0.02 respectively. Our dynamical simulations show the system is stable if the binary orbit has a large semi-major axis and a low eccentricity, which can be verified with future astrometry observations.

The Kepler mission provides a wealth of multiple transiting planet systems (MTPS). The formation and evolution of multi-planet systems are likely to be influenced by companion stars given the abundance of multi stellar systems. We study the influence of stellar companions by measuring the stellar multiplicity rate of MTPS. We select 138 bright (KP < 13.5) Kepler MTPS and search for stellar companions with AO imaging data and archival radial velocity (RV) data. We obtain new AO images for 73 MTPS. Other MTPS in the sample have archival AO imaging data from the Kepler Community Follow-up Observation Program (CFOP). From these imaging data, we detect 42 stellar companions around 35 host stars. For stellar separation 1 AU < a < 100 AU, the stellar multiplicity rate is 5.2 $\pm$ 5.0% for MTPS, which is 2.8{\sigma} lower than 21.1 $\pm$ 2.8% for the control sample, i.e., the field stars in the solar neighborhood. We identify two origins for the deficit of stellar companions within 100 AU to MTPS: (1) a suppressive planet formation, and (2) the disruption of orbital coplanarity due to stellar companions. To distinguish between the two origins, we compare the stellar multiplicity rates of MTPS and single transiting planet systems (STPS). However, current data are not sufficient for this purpose. For 100 AU < a < 2000 AU, the stellar multiplicity rates are comparable for MTPS (8.0 $\pm$ 4.0%), STPS (6.4 $\pm$ 5.8%), and the control sample (12.5 $\pm$ 2.8%).

As hundreds of gas giant planets have been discovered, we study how these planets form and evolve in different stellar environments, specifically in multiple stellar systems. In such systems, stellar companions may have a profound influence on gas giant planet formation and evolution via several dynamical effects such as truncation and perturbation. We select 84 Kepler Objects of Interest (KOIs) with gas giant planet candidates. We obtain high-angular resolution images using telescopes with adaptive optics (AO) systems. Together with the AO data, we use archival radial velocity data and dynamical analysis to constrain the presence of stellar companions. We detect 59 stellar companions around 40 KOIs for which we develop methods of testing their physical association. These methods are based on color information and galactic stellar population statistics. We find evidence of suppressive planet formation within 20 AU by comparing stellar multiplicity. The stellar multiplicity rate for planet host stars is 0$^{+5}_{-0}$\% within 20 AU. In comparison, the stellar multiplicity rate is 18\%$\pm$2\% for the control sample, i.e., field stars in the solar neighborhood. The stellar multiplicity rate for planet host stars is 34\%$\pm$8\% for separations between 20 and 200 AU, which is higher than the control sample at 12\%$\pm$2\%. Beyond 200 AU, stellar multiplicity rates are comparable between planet host stars and the control sample. We discuss the implications of the results to gas giant planet formation and evolution.

The Chinese Small Telescope ARray (CSTAR) is the first telescope facility built at Dome A, Antarctica. During the 2008 observing season, the installation provided long-baseline and high-cadence photometric observations in the i-band for 18,145 targets within 20 deg2 CSTAR field around the South Celestial Pole for the purpose of monitoring the astronomical observing quality of Dome A and detecting various types of photometric variability. Using sensitive and robust detection methods, we discover 274 potential variables from this data set, 83 of which are new discoveries. We characterize most of them, providing the periods, amplitudes and classes of variability. The catalog of all these variables is presented along with the discussion of their statistical properties.

\We present the sixth catalog of Kepler candidate planets based on nearly 4 years of high precision photometry. This catalog builds on the legacy of previous catalogs released by the Kepler project and includes 1493 new Kepler Objects of Interest (KOIs) of which 554 are planet candidates, and 131 of these candidates have best fit radii <1.5 R_earth. This brings the total number of KOIs and planet candidates to 7305 and 4173 respectively. We suspect that many of these new candidates at the low signal-to-noise limit may be false alarms created by instrumental noise, and discuss our efforts to identify such objects. We re-evaluate all previously published KOIs with orbital periods of >50 days to provide a consistently vetted sample that can be used to improve planet occurrence rate calculations. We discuss the performance of our planet detection algorithms, and the consistency of our vetting products. The full catalog is publicly available at the NASA Exoplanet Archive.

Almost half of the stellar systems in the solar neighborhood are made up of multiple stars. In multiple-star systems, planet formation is under the dynamical influence of stellar companions, and the planet occurrence rate is expected to be different from that for single stars. There have been numerous studies on the planet occurrence rate of single star systems. However, to fully understand planet formation, the planet occurrence rate in multiple-star systems needs to be addressed. In this work, we {{infer}} the planet occurrence rate in multiple-star systems by measuring the stellar multiplicity rate for planet host stars. For a sub-sample of 56 $Kepler$ planet host stars, we use adaptive optics (AO) imaging and the radial velocity (RV) technique to search for stellar companions. The combination of these two techniques results in high search completeness for stellar companions. We detect 59 visual stellar companions to 25 planet host stars with AO data. {{Three stellar companions are within 2$^{\prime\prime}$, and 27 within 6$^{\prime\prime}$. We also detect 2 possible stellar companions (KOI 5 and KOI 69) showing long-term RV acceleration.}} After correcting for a bias against planet detection in multiple-star systems due to flux contamination, we find that planet formation is suppressed in multiple-star systems with separations smaller than 1500 AU. Specifically, we find that compared to single star systems, planets in multiple-star systems occur $4.5\pm3.2$, $2.6\pm1.0$, and $1.7\pm0.5$ times less frequently when a stellar companion is present at a distance of 10, 100, and 1000 AU, respectively. This conclusion applies only to circumstellar planets; the planet occurrence rate for circumbinary planets requires further investigation.

A transiting planet exhibits sinusoidal transit-time-variations (TTVs) if perturbed by a companion near a mean-motion-resonance (MMR). We search for sinusoidal TTVs in more than 2600 Kepler candidates, using the publicly available Kepler light-curves (Q0-Q12). We find that the TTV fractions rise strikingly with the transit multiplicity. Systems where four or more planets transit enjoy roughly five times higher TTV fraction than those where a single planet transits, and about twice higher than those for doubles and triples. In contrast, models in which all transiting planets arise from similar dynamical configurations predict comparable TTV fractions among these different systems. One simple explanation for our results is that there are at least two different classes of Kepler systems, one closely packed and one more sparsely populated.

Following on from Paper I in our series (Xie 2013), we report the confirmation by Transit Timing Variations (TTVs) of a further 30 planets in 15 multiple planet systems, using the publicly available Kepler light curves (Q0-Q16). All of these fifteen pairs are near first-order Mean Motion Resonances (MMR), showing sinusoidal TTVs consistent with theoretically predicted periods, which demonstrate they are orbiting and interacting in the same systems. Although individual masses cannot be accurately extracted based only on TTVs (because of the well known degeneracy between mass and eccentricity), the measured TTV phases and amplitudes can still place relatively tight constraints on their mass ratios and upper limits on their masses, which confirm their planetary nature. Some of these systems (KOI-274, KOI-285, KOI-370 and KOI-2672) are relatively bright and thus suitable for further follow-up observations.

The planet occurrence rate for multiple stars is important in two aspects. First, almost half of stellar systems in the solar neighborhood are multiple systems. Second, the comparison of the planet occurrence rate for multiple stars to that for single stars sheds light on the influence of stellar multiplicity on planet formation and evolution. We developed a method of distinguishing planet occurrence rates for single and multiple stars. From a sample of 138 bright (KP <13.5) Kepler multi-planet candidate systems, we compared the stellar multiplicity rate of these planet host stars to that of field stars. Using dynamical stability analyses and archival Doppler measurements, we find that the stellar multiplicity rate of planet host stars is significantly lower than field stars for semi-major axes less than 20 AU, suggesting that planet formation and evolution are suppressed by the presence of a close-in companion star at these separations. The influence of stellar multiplicity at larger separations is uncertain because of search in- completeness due to a limited Doppler observation time baseline and a lack of high resolution imaging observation. We calculated the planet confidence for the sample of mutlti-planet candidates, and find that the planet confidences for KOI 82.01, KOI 115.01, KOI 282.01 and KOI 1781.02 are higher than 99.7% and thus validate the planetary nature of these four planet candidates. This sample of bright Kepler multi-planet candidates with refined stellar and orbital parameters, planet confidence estimation, and nearby stellar companion identification offers a well-characterized sample for future theoretical and observational study.

We extract Transit Timing Variation (TTV) signals for 12 pairs of transiting planet candidates that are near first-order Mean Motion Resonances (MMR), using publicly available Kepler light curves (Q0-Q14). These pairs show significant sinusoidal TTVs with theoretically predicted periods, which demonstrate these planet candidates are orbiting and interacting in the same system. Although individual masses cannot be accurately extracted based only on TTVs because of the well known degeneracy between mass and eccentricity, TTV phases and amplitudes can still place upper limits on the masses of the candidates, confirming their planetary nature. Furthermore, the mass ratios of these planet pairs can be relatively tight constrained using these TTVs. The planetary pair in Kepler-82 (KOI-880) seems to have a particularly high mass ratio and density ratio, which might indicate very different internal compositions of these two planets. Some of these newly confirmed planets are also near MMR with other candidates in the system, forming unique resonance chains, e.g., Kepler-80 (KOI-500).

The latitudinal distributions of the yearly mean rotation rates measured respectively by Suzuki in 1998 and 2012 and Pulkkinen $\&$ Tuominen in 1998 are utilized to investigate internal-cycle variation of solar differential rotation. The rotation rate at the solar Equator seems to decrease since cycle 10 onwards. The coefficient $B$ of solar differential rotation, which represents the latitudinal gradient of rotation, is found smaller in the several years after the minimum of a solar cycle than in the several years after the maximum time of the cycle, and it peaks several years after the maximum time of the solar cycle. The internal-cycle variation of the solar rotation rates looks similar in profile to that of the coefficient $B$. A new explanation is proposed to address such a solar-cycle related variation of the solar rotation rates. Weak magnetic fields may more effectively reflect differentiation at low latitudes with high rotation rates than at high latitudes with low rotation rates, and strong magnetic fields may more effectively repress differentiation at relatively low latitudes than at high latitudes. The internal-cycle variation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation.

Solar-cycle related variation of differential rotation is investigated through analyzing the rotation rates of magnetic fields, distributed along latitudes and varying with time at the time interval of August 1976 to April 2008. More pronounced differentiation of rotation rates is found to appear at the ascending part of a Schwabe cycle than at the descending part on an average. The coefficient $B$ in the standard form of differential rotation, which represents the latitudinal gradient of rotation, may be divided into three parts within a Schwabe cycle. Part one spans from the start to the $4^{th}$ year of a Schwabe cycle, within which the absolute $B$ is approximately a constant or slightly fluctuates. Part two spans from the $4^{th}$ to the $7^{th}$ year, within which the absolute $B$ decreases. Part three spans from the $7^{th}$ year to the end, within which the absolute $B$ increases. Strong magnetic fields repress differentiation of rotation rates, so that rotation rates show less pronounced differentiation, but weak magnetic fields seem to just reflect differentiation of rotation rates. The solar-cycle related variation of solar differential rotation is inferred to the result of both the latitudinal migration of the surface torsional pattern and the repression of strong magnetic activity to differentiation of rotation rates.

Compared to the giant planets in the solar system, exoplanets have many remarkable properties such as the prevalence of giant planets on eccentric orbits and the presence of hot Jupiters. Planet-planet scattering (PPS) between giant planets is a possible mechanism in interpreting above and other observed properties. If the observed giant planet architectures are indeed the outcomes of PPS, such drastic dynamical process must affect their primordial moon systems. In this Letter, we discuss the effect of the PPS on the survival of their regular moons. From the viewpoint of observations, some preliminary conclusions are drawn from the simulations. 1. PPS is a destructive process to the moon systems, single planets on eccentric orbits are not the ideal moon-search targets. 2. If hot Jupiters formed through PPS, their original moons have little chance to survive. 3. Planets in multiple systems with small eccentricities are more likely holding their primordial moons. 4. Compared to the lower-mass planets, the massive ones in multiple systems may not be the preferred moon-search targets if the system underwent a PPS history.

The metallic gas associated with the Beta Pic debris disk is not believed to be primordial, but arises from the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are exceptionally overabundant compared to other elements, by some 400 times. We study the origin of this enrichment under two opposing hypothesis, preferential production, where the gas is produced with the observed unusual abundance (as may happen if gas is produced by photo-desorption from C/O-rich icy grains), and preferential depletion, where the gas evolves to the observed state from an original solar abundance (if outgassing occurs under high-speed collisions) under a number of dynamical processes. We include ... ... We find ... ...

According to the core accretion theory, circumbinary embryos can form only beyond a critical semimajor axis (CSMA). However, due to the relatively high density of solid materials in the inner disk, significant amount of small planetesimals must exist in the inner zone when embryos were forming outside this CSMA. So embryos migration induced by the planetesimal swarm is possible after the gas disk depletion. Through numerical simulations, we found (i) the scattering-driven inward migration of embryos is robust, planets can form in the habitable zone if we adopt a mass distribution of MMSN-like disk; (ii) the total mass of the planetesimals in the inner region and continuous embryo-embryo scattering are two key factors that cause significant embryo migrations; (iii) the scattering-driven migration of embryos is a natural water-deliver mechanism. We propose that planet detections should focus on the close binary with its habitable zone near CSMA.

Most planet pairs in the Kepler data that have measured transit time variations (TTV) are near first-order mean-motion resonances. We derive analytical formulae for their TTV signals. We separate planet eccentricity into free and forced parts, where the forced part is purely due to the planets' proximity to resonance. This separation yields simple analytical formulae. The phase of the TTV depends sensitively on the presence of free eccentricity: if the free eccentricity vanishes, the TTV will be in phase with the longitude of conjunctions. This effect is easily detectable in current TTV data. The amplitude of the TTV depends on planet mass and free eccentricity, and it determines planet mass uniquely only when the free eccentricity is sufficiently small. We proceed to analyze the TTV signals of six short period Kepler pairs. We find that three of these pairs (Kepler-18,24,25) have TTV phase consistent with zero. The other three (Kepler-23,28,32) have small TTV phases, but ones that are distinctly non-zero. We deduce that the free eccentricities of the planets are small, < 0.01, but not always vanishing. Furthermore, as a consequence of this, we deduce that the true masses of the planets are fairly accurately determined by the TTV amplitudes, within a factor < 2. The smallness of the free eccentricities suggests that the planets have experienced substantial dissipation. This is consistent with the hypothesis that the observed pile-up of Kepler pairs near mean-motion resonances is caused by resonant repulsion. But the fact that some of the planets have non-vanishing free eccentricity suggests that after resonant repulsion occurred there was a subsequent phase in the planets' evolution when their eccentricities were modestly excited, perhaps by interplanetary interactions.

With the increasing number of detected exoplanet samples, the statistical properties of planetary systems have become much clearer. In this review, we summarize the major statistics that have been revealed mainly by radial velocity and transiting observations, and try to interpret them within the scope of the classical core-accretion scenario of planet formation, especially in the formation of different orbital architectures for planetary systems around main sequence stars. Based on the different possible formation routes for different planet systems, we tentatively classify them into three major catalogs: hot Jupiter systems, standard systems and distant giant planet systems. The standard system can be further categorized into three sub-types under different circumstances: solar-like systems, hot Super-Earth systems, sub-giant planet systems. We also review the planet detection and formation in binary systems as well as planets in star clusters.

Long-term variations of solar differential rotation and sunspot activity are investigated through re-analyzing the data on parameters of the differential rotation law obtained by Makarov, Tlatov, and Callebaut (1997), Javaraiah, Bertello, and Ulrich (2005a, b), and Javaraiah et al. (2009). Our results indicate that the solar surface rotation rate at the Equator (indicated by the A parameter of the standard solar rotation law) shows a secular decrease since cycle 12 onwards, given by about $1\,-\,1.5\times10^{-3}$($deg\ day^{-1} year^{-1}$). The B parameter of the standard differential rotation law seems to also show a secular decrease since cycle 12 onwards, but of weak statistical significance. The rotation rate averaged on latitudes ($0^{o}\,--\,40^{o}$) does not show a secular trend of statistical significance. Moreover, the average sunspot area shows a secular increase of statistical significance since cycle 12 onwards, while a negative correlation is found between the level of sunspot activity (indicated by the average sunspot area) and the solar equatorial rotation in the long run.

We investigate the effects of gas-disk gravity on the planetesimal dynamics in inclined binary systems, where the circumprimary disk plane is tilted by a significant angle ($i_B$) with respect to the binary disk plane. Our focus is on the Lidov-Kozai mechanism and the evolution of planetesimal eccentricity and inclination. Using both analytical and numerical methods, we find that, on one hand, the disk gravity generally narrows down the Kozai-on region, i.e., the Lidov-Kozai effect can be suppressed in certain parts of (or even the whole of) the disk, depending on various parameters. In the Kozai-off region, planetesimals would move on orbits close to the mid-plane of gas-disk, with the relative angle ($i^{'}$) following a small amplitude periodical oscillation. On the other hand, when we include the effects of disk gravity, we find that the Lidov-Kozai effect can operate even at arbitrarily low inclinations ($i_B$), although lower $i_B$ leads to a smaller Kozai-on region. Furthermore, in the Kozai-on region, most planetesimals' eccentricities can be excited to extremely high values ($\sim 1$), and such extreme high eccentricities usually accompany orbital flipping, i.e., planetesimal orbit flips back and forth between anterograde and retrograde. Once a planetesimal reaches very high orbital eccentricity, gas drag damping will shrink the planetesimal orbit, forming a "hot planetesimal" on a near circular orbit very close to the primary star. Such a mechanism, if replacing the planetesimals and gas drag damping with Jupiters and tidal damping respectively, may lead to frequent production of hot-Jupiters.

The daily sunspot numbers of the whole disk as well as the northern and southern hemispheres from January 1, 1945 to December 31, 2010 are used to investigate the temporal variation of the rotational cycle length through the continuous wavelet transformation analysis method. The auto-correlation function analysis of daily hemispheric sunspot numbers shows that the southern hemisphere rotates faster than the northern hemisphere. The results obtained from the wavelet transformation analysis are: there exists no direct relationship between the variation trend of the rotational cycle length and the variation trend of solar activity in the two hemispheres; the rotational cycle length of both hemispheres has no significant period appearing at the 11 years, but has significant period of about 7.6 years. Analysis concerning the solar cycle dependence of the rotational cycle length shows that in the whole disk and the northern hemisphere acceleration seems to appear before the minimum time of solar activity. Furthermore, the cross-correlation study indicates that the rotational cycle length of the two hemispheres has different phases, and the rotational cycle length of the whole disk as well as the northern and southern hemispheres also has phase shifts with the corresponding solar activity. What's more, the temporal variation of North-South (N-S) asymmetry of the rotational cycle length is also studied; it displays the same variation trend as the N-S asymmetry of solar activity in a solar cycle as well as in the considered time interval, and it has two significant periods of 7.7 and 17.5 years. Moreover, the N-S asymmetry of the rotational cycle length and the N-S asymmetry of solar activity are highly correlated. It's inferred that the northern hemisphere should rotate faster at the beginning of solar cycle 24.

Most detected planet-bearing binaries are in wide orbits, for which a high inclination, $i_B$, between the binary orbital plane and the plane of the planetary disk around the primary is likely to be common. In this paper, we investigate the intermediate stages - from planetesimals to planetary embryos/cores - of planet formation in such highly inclined cases. Our focus is on the effects of gas drag on the planetesimals' orbital evolution, in particular on the evolution of the planetesimals' semimajor axis distribution and their mutual relative velocities. We first demonstrate that a non-evolving axisymmetric disk model is a good approximation for studying the effects of gas drag on a planetesimal in the highly inclined case ($30^\circ <i_B<150^\circ$). We then find that gas drag plays a crucial role, and the results can be generally divided into two categories, i.e., the Kozai-on regime and the Kozai-off regime, depending on the specific value of $i_B$. For both regimes, a robust outcome over a wide range of parameters is that, planetesimals migrate/jump inwards and pile up, leading to a severely truncated and dense planetesimal disk around the primary. In this compact and dense disk, collision rates are high but relative velocities are low, providing conditions which are favorable for planetesimal growth, and potentially allow for the subsequent formation of planets.

We present a new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

We present ground-based optical observations of the September 2009 and January 2010 transits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a midtransit time of T_c [HJD] = 2455099.196 +- 0.026, which is about 1 sigma away from the previously predicted time. We observed the January 2010 event from 9 different locations, with most phases of the transit being observed by at least 3 different teams. We determine a midtransit time of T_c [HJD] = 2455210.6502 +- 0.0064, which is within 1.3 sigma of the time derived from a Spitzer observation of the same event.

The standard model of planet formation considers an initial phase in which planetesimals form from a dust disk, followed by a phase of mutual planetesimal-planetesimal collisions, leading eventually to the formation of planetary embryos. However, there is a potential transition phase (which we call the "snowball phase"), between the formation of the first planetesimals and the onset of mutual collisions amongst them, which has often been either ignored or underestimated in previous studies. In this snowball phase, isolated planetesimals move on Keplerian orbits and grow solely via the direct accretion of sub-cm sized dust entrained with the gas in the protoplanetary disk. Using a simplified model in which planetesimals are progressively produced from the dust, we consider the expected sizes to which the planetesimals can grow before mutual collisions commence and derive the dependence of this size on a number of critical parameters, including the degree of disk turbulence, the planetesimal size at birth and the rate of planetesimal creation. For systems in which turbulence is weak and the planetesimals are created at a low rate and with relatively small birth size, we show that the snowball growth phase can be very important, allowing planetesimals to grow by a factor of 10^6 in mass before mutual collisions take over. In such cases, the snowball growth phase can be the dominant mode to transfer mass from the dust to planetesimals. Moreover, such growth can take place within the typical lifetime of a protoplanetary gas disk. A noteworthy result is that ... ...(see the paper). For the specific case of close binaries such as Alpha Centauri ... ... (see the paper). From a more general perspective, these preliminary results suggest that an efficient snowball growth phase provides a large amount of "room at the bottom" for theories of planet formation.

Stellar perturbations affect planet-formation in binary systems. Recent studies show that the planet-formation stage of mutual accretion of km-sized planetesimals is most sensitive to binary effects. In this paper, the condition for planetesimal accretion is investigated around Alpha CenB, which is believed to be an ideal candidate for detection of an Earth-like planet in or near its habitable zone(0.5-0.9 AU). A simplified scaling method is developed to estimate the accretion timescale of the planetesimals embedded in a protoplanetary disk. Twenty-four cases with different binary inclinations(i_B=0, 0.1, 1.0, and 10 deg), gas densities(0.3,1,and 3 times of the Minimum Mass of Solar Nebula, MMSN hereafter), and with and without gas depletion, are simulated. We find: (1)...(2)...(3)...(4)...(see the paper for details). In other words, our results suggest that formation of Earth-like planets through accretion of km-sized planetesimals is possible in Alpha CenB, while formation of gaseous giant planets is not favorable.

Recent observations show that planet can reside in close binary systems with stellar separation of only about 20 AU. However, planet formation in such close binary systems is a challenge to current theory. One of the major theoretical problems occurs in the intermediate stage-planetesimals accretion into planetary embryos-during which the companion's perturbations can stir up the relative velocites(dV) of planetesimals and thus slow down or even cease their growth. However, all previous studies assumed a 2-dimentional (2D) disk and a coplanar binary orbit. Extending previous studies by including a 3D gas disk and an inclined binary orbit with small relative inclination of i_B=0.1-5 deg, we numerically investigate the conditions for planetesimal accretion at 1-2 AU, an extension of the habitable zone(1-1.3 AU), around alpha Centauri A in this paper. Inclusion of the binary inclination leads to: (1) differential orbital phasing is realized in the 3D space, and thus different-sized bodies are separated from each other; (2) total impact rate becomes lower, and impacts mainly occur between similar-sized bodies; (3) accretion is more favored, but the balance between accretion and erosion remains uncertain, and the "possible accretion region" extends up to 2AU when assuming an optimistic Q*(critical specific energy that leads to catastrophic fragmentation); and (4) impact velocities (dV) are significantly reduced but still much larger than their escape velocities, which infers that planetesimals grow by means of type II runaway mode. As a conclusion, inclusion of a small binary inclination is a promising mechanism that favors accretion, opening a possibility that planet formation in close binary systems can go through the difficult stage of planetesimals accretion into planetary embryos.

Currently, one of major problems concerning planet formation theory in close binary systems is, the strong perturbation from the companion star can increase relative velocities ($\triangle V$) of planetesimals around the primary and thus hinder their growth. According to previous studies, while gas drag can reduce the $\triangle V$ between bodies of the same sizes by forcing orbital alignment to planetesimals, it increases the $\triangle V$ among bodies of different sizes. In this paper, focusing on the $\gamma$ Cephei binary system, we propose a mechanism that can overcome this difficulty. We show that in a dissipating gas disk (with a typical dissipating timescale of $\sim 10^5-10^6$ years), all the planetesimals eventually converge towards the same forced orbits regardless of their sizes, leading to much lower impact velocities among them. These $\triangle V$ decrease processes progressively increase net mass accretion and even trigger runaway growth for large bodies (radius $>15$ km). The effect of size distribution of planetesimals is discussed, and found to be one of the dominant factors that determine the outcome of collisional evolution. Anyway, it can be concluded that by including the gas dissipation in the early stage of disk evolution, the conditions for planetesimal accretion become much better, and the process from planetesimal to planet-embryo can be carried out in close binary systems like $\gamma$ Cephei.