We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $\sigma_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.

The goal of Project GAUSS is to return samples from the dwarf planet Ceres. Ceres is the most accessible ocean world candidate and the largest reservoir of water in the inner solar system. It shows active cryovolcanism and hydrothermal activities in recent history that resulted in minerals not found in any other planets to date except for Earth's upper crust. The possible occurrence of recent subsurface ocean on Ceres and the complex geochemistry suggest possible past habitability and even the potential for ongoing habitability. Aiming to answer a broad spectrum of questions about the origin and evolution of Ceres and its potential habitability, GAUSS will return samples from this possible ocean world for the first time. The project will address the following top-level scientific questions: 1) What is the origin of Ceres and the origin and transfer of water and other volatiles in the inner solar system? 2) What are the physical properties and internal structure of Ceres? What do they tell us about the evolutionary and aqueous alteration history of icy dwarf planets? 3) What are the astrobiological implications of Ceres? Was it habitable in the past and is it still today? 4) What are the mineralogical connections between Ceres and our current collections of primitive meteorites? GAUSS will first perform a high-resolution global remote sensing investigation, characterizing the geophysical and geochemical properties of Ceres. Candidate sampling sites will then be identified, and observation campaigns will be run for an in-depth assessment of the candidate sites. Once the sampling site is selected, a lander will be deployed on the surface to collect samples and return them to Earth in cryogenic conditions that preserves the volatile and organic composition as well as the original physical status as much as possible.

Small planets are common around late-M dwarfs and can be detected through highly precise photometry by the transit method. Planets orbiting nearby stars are particularly important as they are often the best-suited for future follow-up studies. We present observations of three nearby M-dwarfs referred to as EIC-1, EIC-2, and EIC-3, and use them to search for transits and set limits on the presence of planets. On most nights our observations are sensitive to Earth-sized transiting planets, and photometric precision is similar to or better than TESS for faint late-M dwarfs of the same magnitude (I=15 mag). We present our photometry and transit search pipeline, which utilizes simple median detrending in combination with transit least squares based transit detection (Hippke & Heller 2019).For these targets, and transiting planets between one and two Earth radii, we achieve an average transit detection probability of 60% between periods of 0.5 and 2 days, 30% between 2 and 5 days,and 10% between 5 and 10 days. These sensitivities are conservative compared to visual searches.

Many asteroids in the main belt are continuously pushed by Yarkovsky effect into regions of different mean motion resonances (MMRs) and then ejected out. They are considered as the principal source of near-Earth objects. We investigate in this paper the effects of the 7/3 MMR with Jupiter (J7/3 MMR) on the transportation of asteroids from Koronis and Eos families that reside respectively on the inner and outer side of the resonance. The fraction of asteroids that make successful crossing through the resonance and the escaping rate from the resonance are found to depend on the Yarkovsky drifting rate, the initial inclination and the migrating direction. The excitation of eccentricity and inclination due to the combined influence from both the resonance and Yarkovsky effect is discussed. Only the eccentricity can be pumped up considerably, and it is attributed mainly to the resonance. In the observational data, family members are also found in the resonance and on the opposite side of the resonance with respect to the corresponding family centre. The existence of these family members is explained using our results of numerical simulations. Finally, the replenishment of asteroids in the J7/3 MMR and the transportation of asteroids by it are discussed.

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $\Delta a$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $\Delta a$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $\Delta a$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $\omega_{\rm rot}$ of asteroid depends on its size $R$ by $\omega_{\rm rot}\propto R^{-1}$.

Near-Earth Objects (NEOs) that orbit the Sun on or within Earth's orbit are tricky to detect for Earth-based observers due to their proximity to the Sun in the sky. These small bodies hold clues to the dynamical history of the inner solar system as well as the physical evolution of planetesimals in extreme environments. Populations in this region include the Atira and Vatira asteroids, as well as Venus and Earth co-orbital asteroids. Here we present a twilight search for these small bodies, conducted using the 1.2-m Oschin Schmidt and the Zwicky Transient Facility (ZTF) camera at Palomar Observatory. The ZTF twilight survey operates at solar elongations down to $35^\circ$ with limiting magnitude of $r=19.5$. During a total of 40 evening sessions and 62 morning sessions conducted between 2018 November 15 and 2019 June 23, we detected 6 Atiras, including 2 new discoveries 2019 AQ$_3$ and 2019 LF$_6$, but no Vatiras or Earth/Venus co-orbital asteroids. NEO population models show that these new discoveries are likely only the tip of the iceberg, with the bulk of the population yet to be found. The population models also suggest that we have only detected 5--$7\%$ of the $H<20$ Atira population over the 7-month survey. Co-orbital asteroids are smaller in diameters and require deeper surveys. A systematic and efficient survey of the near-Sun region will require deeper searches and/or facilities that can operate at small solar elongations.

133P/Elst-Pizarro is the firstly recognized main-belt comet, but we still know little about its nucleus. Firstly we use mid-infrared data of Spitzer-MIPS, Spitzer-IRS and WISE to estimate its effective diameter $D_{\rm eff}=3.9^{+0.4}_{-0.3}$ km, geometric albedo $p_{\rm v}=0.074\pm0.013$ and mean Bond albedo $A_{\rm eff,B}=0.024\pm0.004$. The albedo is used to compute 133P's temperature distribution, which shows significant seasonal variation, especially polar regions, ranging from $\sim40$ to $\sim200$ K. Based on current activity observations, the maximum water gas production rate is estimated to be $\sim1.4\times10^{23}\rm~s^{-1}$, being far weaker than $\sim10^{26}\rm~s^{-1}$ of JFC 67P at similar helio-centric distance $\sim2.7$ AU, indicating a thick dust mantle on the surface to lower down the gas production rate. The diameter of the sublimation area may be $<\sim200$ m according to our model prediction. We thus propose that 133P's activity is more likely to be caused by sublimation of regional near-surface ice patch rather than homogeneous buried ice layer. Such small near-surface ice patch might be exposed by one impact event, before which 133P may be an extinct comet (or ice-rich asteroid) with ice layer buried below $\sim$40 m depth. The proposed ice patch may be located somewhere within latitude $-50\sim50^\circ$ by comparing theoretical variation of sublimation temperature to the constraints from observations. The time scale to form such a thick dust mantle is estimated to be $>100$ Myr, indicating that 133P may be more likely to be a relatively old planetesimals or a member of an old family than a recently formed fragment of some young family.

The intra-night trajectories of asteroids can be approximated by straight lines, and so are their intra-night detections. Therefore, the Hough transform, a line detecting algorithm, can be used to connect the line-up detections to find asteroids. We applied this algorithm to a high-cadence Pan-STARRS 1 (PS1) observation, which was originally designed to collect asteroid light curves for rotation period measurements (Chang et al., 2019). The algorithm recovered most of the known asteroids in the observing fields and, moreover, discovered 3574 new asteroids with magnitude mainly of 21.5 < w_p1 < 22.5 mag. This magnitude range is equivalent to sub-kilometer main-belt asteroids (MBAs), which usually lack of rotation period measurements due to their faintness. Using the light curves of the 3574 new asteroids, we obtained 122 reliable rotation periods, of which 13 are super-fast rotators (SRFs; i.e., rotation period of < 2 hr). The required cohesion to survive these SFRs range from tens to thousands of Pa, a value consistent with the known SFRs and the regolith on the Moon and Mars. The higher chance of discovering SFRs here suggests that sub-kilometer MBAs probably harbor more SFRs.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

We describe ZStreak, a semi-real-time pipeline specialized in detecting small, fast-moving near-Earth asteroids (NEAs) that is currently operating on the data from the newly-commissioned Zwicky Transient Facility (ZTF) survey. Based on a prototype originally developed by Waszczak et al. (2017) for the Palomar Transient Factory (PTF), the predecessor of ZTF, ZStreak features an improved machine-learning model that can cope with the $10\times$ data rate increment between PTF and ZTF. Since its first discovery on 2018 February 5 (2018 CL), ZTF/ZStreak has discovered $45$ confirmed new NEAs over a total of 232 observable nights until 2018 December 31. Most of the discoveries are small NEAs, with diameters less than $\sim100$ m. By analyzing the discovery circumstances, we find that objects having the first to last detection time interval under 2 hr are at risk of being lost. We will further improve real-time follow-up capabilities, and work on suppressing false positives using deep learning.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

The Zwicky Transient Facility (ZTF) is a new optical time-domain survey that uses the Palomar 48-inch Schmidt telescope. A custom-built wide-field camera provides a 47 deg$^2$ field of view and 8 second readout time, yielding more than an order of magnitude improvement in survey speed relative to its predecessor survey, the Palomar Transient Factory (PTF). We describe the design and implementation of the camera and observing system. The ZTF data system at the Infrared Processing and Analysis Center provides near-real-time reduction to identify moving and varying objects. We outline the analysis pipelines, data products, and associated archive. Finally, we present on-sky performance analysis and first scientific results from commissioning and the early survey. ZTF's public alert stream will serve as a useful precursor for that of the Large Synoptic Survey Telescope.

The Zwicky Transient Facility (ZTF), a public-private enterprise, is a new time domain survey employing a dedicated camera on the Palomar 48-inch Schmidt telescope with a 47 deg$^2$ field of view and 8 second readout time. It is well positioned in the development of time domain astronomy, offering operations at 10% of the scale and style of the Large Synoptic Survey Telescope (LSST) with a single 1-m class survey telescope. The public surveys will cover the observable northern sky every three nights in g and r filters and the visible Galactic plane every night in g and r. Alerts generated by these surveys are sent in real time to brokers. A consortium of universities which provided funding ("partnership") are undertaking several boutique surveys. The combination of these surveys producing one million alerts per night allows for exploration of transient and variable astrophysical phenomena brighter than r $\sim$ 20.5 on timescales of minutes to years. We describe the primary science objectives driving ZTF including the physics of supernovae and relativistic explosions, multi-messenger astrophysics, supernova cosmology, active galactic nuclei and tidal disruption events, stellar variability, and Solar System objects.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

A class of asteroids, called large super-fast rotators (large SFRs), have rotation periods shorter than 2 hours and diameters larger than ~ 0.3 km. They pose challenges to the usual interior rubble-pile structure unless a relatively high bulk density is assumed. So far, only six large SFRs have been found. Therefore, we present a survey of asteroid rotation period using Pan-STARRS 1 telescopes during 2016 October 26 to 31 to search more large SFRs and study their properties. A total of 876 reliable rotation periods are measured, among which seven are large SFRs, thereby increasing the inventory of known large SFRs. These seven newly discovered large SFRs have diverse colors and locations in the main asteroid belt, suggesting that the taxonomic tendency and the location preference in the inner main belt of the six perviously known large SFRs could be a bias due to various observational limits. Interestingly, five out of the seven newly discovered large SFRs are mid main-belt asteroids. Considering the rare discovery rates of large SFR in the previously similar surveys (Chang et al., 2015, 2016) and the survey condition in this work, the chance of detecting a large SFR in the inner main belt seems to be relatively low. This probably suggests that the inner main belt harbors less large SFRs than the mid main belt. From our survey, we also found the drop in number appearing at f > 5 rev/day on the spin-rate distribution for the outer main-belt asteroids of D < 3 km, which was reported for the inner and mid main belt by (Chang et al., 2015, 2016).

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

The bulk composition of the Earth is dramatically carbon poor compared to that of the interstellar medium, and this phenomenon extends to the asteroid belt. To interpret this carbon deficit problem, the carbonaceous component in grains must have been converted into the gas-phase in the inner regions of protoplanetary disks prior to planetary formation. We examine the effect of carbon grain destruction on the chemical structure of disks by calculating the molecular abundances and distributions using a comprehensive chemical reaction network. When carbon grains are destroyed and the elemental abundance of the gas becomes carbon-rich, the abundances of carbon-bearing molecules, such as HCN and carbon-chain molecules, increase dramatically near the midplane, while oxygen-bearing molecules, such as H$_2$O and CO$_2$, are depleted. We compare the results of these model calculations with the solid carbon fraction in the solar system. Although we find a carbon depletion gradient, there are some quantitative discrepancies: the model shows a higher value at the position of asteroid belt and a lower value at the location of the Earth. In addition, using the obtained molecular abundances distributions, coupled with line radiative transfer calculations, we make predictions for ALMA to potentially observe the effect of carbon grain destruction in nearby protoplanetary disks. The results indicate that HCN, H$^{13}$CN and c-C$_3$H$_2$ may be good tracers.

A long-term sublimation model to explain how Phaethon could provide the Geminid stream is proposed. We find that it would take $\sim6$ Myr or more for Phaethon to lose all of its internal ice (if ever there was) in its present orbit. Thus, if the asteroid moved from the region of a 5:2 or 8:3 mean motion resonance with Jupiter to its present orbit less than $1$ Myr ago, it may have retained much of its primordial ice. The dust mantle on the sublimating body should have a thickness of at least $15$ m but the mantle could have been less than $1$ m thick $1000$ years ago. We find that the total gas production rate could have been as large as $10^{27}\rm~s^{-1}$ then, and the gas flow could have been capable of lifting dust particles of up to a few centimeters in size. Therefore, gas production during the past millennium could have been sufficient to blow away enough dust particles to explain the entire Geminid stream. For present-day Phaethon, the gas production is comparatively weak. But strong transient gas release with a rate of $\sim4.5\times10^{19}\rm~m^{-2}s^{-1}$ is expected for its south polar region when Phaethon moves from $0^\circ$ to $2^\circ$ mean anomaly near perihelion. Consequently, dust particles with radii of $<\sim260~\mu m$ can be blown away to form a dust tail. In addition, we find that the large surface temperature variation of $>600$ K near perihelion can generate sufficiently large thermal stress to cause fracture of rocks or boulders and provide an efficient mechanism to produce dust particles on the surface. The time scale for this process should be several times longer than the seasonal thermal cycle, thereby dominating the cycle of appearance of the dust tail.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet.

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

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

Machine learning techniques are widely applied in many modern optical sky surveys, e.q. Pan-STARRS1, PTF/iPTF and Subaru/Hyper Suprime-Cam survey, to reduce human intervention for data verification. In this study, we have established a machine learning based real-bogus system to reject the false detections in the Subaru/Hyper-Suprime-Cam StrategicSurvey Program (HSC-SSP) source catalog. Therefore the HSC-SSP moving object detection pipeline can operate more effectively due to the reduction of false positives. To train the real-bogus system, we use the stationary sources as the real training set and the "flagged" data as the bogus set. The training set contains 47 features, most of which are photometric measurements and shape moments generated from the HSC image reduction pipeline (hscPipe). Our system can reach a true positive rate (tpr) ~96% with a false positive rate (fpr) ~ 1% or tpr ~99% at fpr ~5%. Therefore we conclude that the stationary sources are decent real training samples, and using photometry measurements and shape moments can reject the false positives effectively.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

(144977) 2005 EC127 is an V-/A-type inner-main-belt asteroid with a diameter of 0.6 +- 0.1 km. Asteroids of this size are believed to have rubble-pile structure, and, therefore, cannot have a rotation period shorter than 2.2 hours. However, our measurements show that asteroid 2005 EC127 completes one rotation in 1.65 +- 0.01 hours with a peak-to-peak light-curve variation of ~0.5 mag. Therefore, this asteroid is identified as a large super-fast rotator. Either a rubble-pile asteroid with a bulk density of ~6 g cm^-3 or an asteroid with an internal cohesion of 47 +- 30 Pa can explain 2005 EC127. However, the scenario of high bulk density is very unlikely for asteroids. To date, only six large super-fast rotators, including 2005 EC127, have been reported, and this number is very small when compared with the much more numerous fast rotators. We also note that none of the six reported large SFRs are classified as C-type asteroids.

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.

The wide-field synoptic sky surveys, known as the Palomar Transient Factory (PTF) and the intermediate Palomar Transient Factory (iPTF), will accumulate a large number of known and new RR Lyrae. These RR Lyrae are good tracers to study the substructure of the Galactic halo if their distance, metallicity, and galactocentric velocity can be measured. Candidates of halo RR Lyrae can be identified from their distance and metallicity before requesting spectroscopic observations for confirmation. This is because both quantities can be obtained via their photometric light curves, because the absolute V-band magnitude for RR Lyrae is correlated with metallicity, and the metallicity can be estimated using a metallicity-light curve relation. To fully utilize the PTF and iPTF light-curve data in related future work, it is necessary to derive the metallicity-light curve relation in the native PTF/iPTF R-band photometric system. In this work, we derived such a relation using the known ab-type RR Lyrae located in the Kepler field, and it is found to be $[Fe/H]_{PTF} = -4.089 - 7.346 P + 1.280 \phi_{31}$ (where $P$ is pulsational period and $\phi_{31}$ is one of the Fourier parameters describing the shape of the light curve), with a dispersion of 0.118 dex. We tested our metallicity-light curve relation with new spectroscopic observations of a few RR Lyrae in the Kepler field, as well as several data sets available in the literature. Our tests demonstrated that the derived metallicity-light curve relation could be used to estimate metallicities for the majority of the RR Lyrae, which are in agreement with the published values.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

A sample of the LAMOST spectra of the early type M0-M3 dwarfs is compared with the Kepler observations. It is found that M dwarfs with strong chromospheric emission in $H_{\alpha}$ have large flare activity in general. The rotational periods derived from the Kepler measurements have close correlations with the sizes of the flares, the power-law distribution index and the equivalent widths of the $H_{\alpha}$ emission. A clear trend exists for higher magnetic activities being detected in faster rotating M dwarfs (rotation periods $<$ 20 day).

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

Infrared observations of the coma of 67P/Churyumov-Gerasimenko were carried out from July to September 2015, i.e., around perihelion (13 August 2015), with the high-resolution channel of the VIRTIS instrument onboard Rosetta. We present the analysis of fluorescence emission lines of H$_2$O, CO$_2$, $^{13}$CO$_2$, OCS, and CH$_4$ detected in limb sounding with the field of view at 2.7-5 km from the comet centre. Measurements are sampling outgassing from the illuminated southern hemisphere, as revealed by H$_2$O and CO$_2$ raster maps, which show anisotropic distributions, aligned along the projected rotation axis. An abrupt increase of water production is observed six days after perihelion. In the mean time, CO$_2$, CH$_4$, and OCS abundances relative to water increased by a factor of 2 to reach mean values of 32%, 0.47%, and 0.18%, respectively, averaging post-perihelion data. We interpret these changes as resulting from the erosion of volatile-poor surface layers. Sustained dust ablation due to the sublimation of water ice maintained volatile-rich layers near the surface until at least the end of the considered period, as expected for low thermal inertia surface layers. The large abundance measured for CO$_2$ should be representative of the 67P nucleus original composition, and indicates that 67P is a CO$_2$-rich comet. Comparison with abundance ratios measured in the northern hemisphere shows that seasons play an important role in comet outgassing. The low CO$_2$/H$_2$O values measured above the illuminated northern hemisphere are not original, but the result of the devolatilization of the uppermost layers.

In this work we report the detection of seven Neptune Trojans (NTs) in the Pan-STARRS 1 (PS1) survey. Five of these are new discoveries, consisting of four L4 Trojans and one L5 Trojan. Our orbital simulations show that the L5 Trojan stably librates for only several million years. This suggests that the L5 Trojan must be of recent capture origin. On the other hand, all four new L4 Trojans stably occupy the 1:1 resonance with Neptune for more than 1 Gyr. They can, therefore, be of primordial origin. Our survey simulation results show that the inclination width of the Neptune Trojan population should be between $7^{\circ}$ and $27^{\circ}$ at $>$ 95% confidence, and most likely $\sim 11^{\circ}$. In this paper, we describe the PS1 survey, the Outer Solar System pipeline, the confirming observations, and the orbital/physical properties of the new Neptune Trojans.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

In order to look for large super-fast rotators, five dedicated surveys covering ~ 188 square degree in the ecliptic plane have been carried out in R-band with ~10 min cadence using the intermediate Palomar Transient Factory in late 2014 and early 2015. Among 1029 reliable rotation periods obtained from the surveys, we discovered one new large super-fast rotator, (40511) 1999 RE88, and other 18 candidates. (40511) 1999 RE88 is an S-type inner main-belt asteroid with a diameter of D = 1.9 +- 0.3 km, which has a rotation period of P = 1.96 +- 0.01 hr and a lightcurve amplitude of ~0.1 mag. To maintain such fast rotation, an internal cohesive strength of ~780 Pa is required. Combining all known large super-fast rotators, their cohesive strengths all fall in the range of 100 to 1000 Pa of lunar regolith. However, the number of large super-fast rotators seems to be far less than the whole asteroid population. This might indicate a peculiar asteroid group for them. Although the detection efficiency for a long rotation period is greatly reduced due to our two-day observation time span, the spin-rate distributions of this work show consistent results with Chang et al. (2015) after considering the possible observational bias in our surveys. It shows a number decrease with increase of spin rate for asteroids with diameter of 3 < D < 15 km, and a number drop at spin-rate of f = 5 rev/day for asteroids with D < 3 km.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

Although the majority of Centaurs are thought to have originated in the scattered disk, with the high-inclination members coming from the Oort cloud, the origin of the high inclination component of trans-Neptunian objects (TNOs) remains uncertain. We report the discovery of a retrograde TNO, which we nickname "Niku", detected by the Pan-STARRS 1 Outer Solar System Survey. Our numerical integrations show that the orbital dynamics of Niku are very similar to that of 2008 KV$_{42}$ (Drac), with a half-life of $\sim 500$ Myr. Comparing similar high inclination TNOs and Centaurs ($q > 10$ AU, $a < 100$ AU and $i > 60^\circ$), we find that these objects exhibit a surprising clustering of ascending node, and occupy a common orbital plane. This orbital configuration has high statistical significance: 3.8-$\sigma$. An unknown mechanism is required to explain the observed clustering. This discovery may provide a pathway to investigate a possible reservoir of high-inclination objects.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

We report the discovery of 2 new Be stars, and re-identify one known Be star in the open cluster NGC 6830. Eleven H-alpha emitters were discovered using the H-alpha imaging photometry of the Palomar Transient Factory Survey. Stellar membership of the candidates was verified with photometric and kinematic information using 2MASS data and proper motions. The spectroscopic confirmation was carried out by using the Shane 3-m telescope at Lick observatory. Based on their spectral types, three H-alpha emitters were confirmed as Be stars with H-alpha equivalent widths > -10 Angstrom. Two objects were also observed by the new spectrograph SED-Machine on the Palomar 60 inch Telescope. The SED-Machine results show strong H-alpha emission lines, which are consistent with the results of the Lick observations. The high efficiency of the SED-Machine can provide rapid observations for Be stars in a comprehensive survey in the future.

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

The Karin family is a very young asteroid family created by an asteroid breakup 5.8 Myr ago. Since the members of this family probably have not experienced significant orbital or collisional evolution yet, it is possible that they still preserve properties of the original family-forming event in terms of their spin state. As we carried out a series of photometric observations of the Karin family asteroids, here we report an analysis result of lightcurves including the rotation period of eleven members. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs or smaller MBAs (diameter D<12 km), and even lower than that of larger MBAs (D>130 km). We investigated a correlation between the peak-to-peak variation magnitude reduced to zero solar phase angle and the rotation period of the eleven Karin family asteroids, and found a possible trend that elongated members have lower spin rate, and less elongated members have higher spin rate. However, this has to be confirmed by another series of future observations.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

Rotation period clustering in prograde/retrograde rotators might indicate the preliminary indication of the Slivan state in the Koronis family as a result of the YORP effect. We follow the general scenario of dispersion in semimajor axis of the asteroid family members to separate prograde and retrograde rotators in the Koronis family. From the available rotation periods obtained from PTF/iPTF, we can not found the rotation period clustering of objects with H grater than 12 mag in the Koronis family. This could be the result of the intermittent collisional process of small asteroids (D less than ~20 km) which leads to astray Yarkovsky drifting. Measurement of the pole orientations of our sample will verify our preliminary result and validate our method.

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.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.

Two dedicated asteroid rotation-period surveys have been carried out using data taken on January 6-9 and February 20-23 of 2014 by the Intermediate Palomar Transient Factory (iPTF) in the $R$~band with $\sim 20$-min cadence. The total survey area covered 174~deg$^2$ in the ecliptic plane. Reliable rotation periods for 1,438 asteroids are obtained from a larger data set of 6,551 mostly main-belt asteroids, each with $\geq 10$~detections. Analysis of 1751, PTF based, reliable rotation periods clearly shows the "spin barrier" at $\sim 2$~hours for "rubble-pile" asteroids. We also found a new large-sized super-fast rotator, 2005 UW163 (Chang et al., 2014), and other five candidates as well. Our spin-rate distributions of asteroids with $3 < D < 15$~km shows number decrease when frequency greater than 5 rev/day, which is consistent to that of the Asteroid Light Curve Database (LCDB, Warner et al., 2009) and the result of (Masiero et al., 2009). We found the discrepancy in the spin-rate distribution between our result and (Pravec et al., 2008, update 2014-04-20) is mainly from asteroids with $\Delta m < 0.2$ mag that might be primarily due to different survey strategies. For asteroids with $D \leq 3$~km, we found a significant number drop at $f = 6$ rev/day. The YORP effect timescale for small-sized asteroid is shorter that makes more elongate objets spun up to reach their spin-rate limit and results in break-up. The K-S test suggests a possible difference in the spin-rate distributions of C- and S-type asteroids. We also find that C-type asteroids have a smaller spin-rate limit than the S-type, which agrees with the general sense that the C-type has lower bulk density than the S-type.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.

We fit 54,296 sparsely-sampled asteroid lightcurves in the Palomar Transient Factory to a combined rotation plus phase-function model. Each lightcurve consists of 20+ observations acquired in a single opposition. Using 805 asteroids in our sample that have reference periods in the literature, we find the reliability of our fitted periods is a complicated function of the period, amplitude, apparent magnitude and other attributes. Using the 805-asteroid ground-truth sample, we train an automated classifier to estimate (along with manual inspection) the validity of the remaining 53,000 fitted periods. By this method we find 9,033 of our lightcurves (of 8,300 unique asteroids) have reliable periods. Subsequent consideration of asteroids with multiple lightcurve fits indicate 4% contamination in these reliable periods. For 3,902 lightcurves with sufficient phase-angle coverage and either a reliably-fit period or low amplitude, we examine the distribution of several phase-function parameters, none of which are bimodal though all correlate with the bond albedo and with visible-band colors. Comparing the theoretical maximal spin rate of a fluid body with our amplitude versus spin-rate distribution suggests that, if held together only by self-gravity, most asteroids are in general less dense than 2 g/cm$^3$, while C types have a lower limit of between 1 and 2 g/cm$^3$, in agreement with previous density estimates. For 5-20km diameters, S types rotate faster and have lower amplitudes than C types. If both populations share the same angular momentum, this may indicate the two types' differing ability to deform under rotational stress. Lastly, we compare our absolute magnitudes and apparent-magnitude residuals to those of the Minor Planet Center's nominal $G=0.15$, rotation-neglecting model; our phase-function plus Fourier-series fitting reduces asteroid photometric RMS scatter by a factor of 3.

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 size-dependent effects of asteroids on surface regolith and collisional lifetimes suggest that small asteroids are younger than large asteroids. In this study, we performed multicolor main-belt asteroid (MBA) survey by Subaru telescope/Suprime-Cam to search for subkilometer-sized ordinary chondrite (Q-type) like MBAs. The total survey area was 1.5 deg^2 near ecliptic plane and close to the opposition. We detected 150 MBAs with 4 bands (B, V , R, I) in this survey. The range of absolute magnitude of detected asteroids was between 13 and 22 magnitude, which is equivalent to the size range of kilometer to sub-kilometer diameter in MBAs. From this observation, 75 of 150 MBAs with color uncertainty less than 0.1 were used in the spectral type analysis, and two possible Q-type aster- oids were detected. This mean that the Q-type to S-type ratio in MBAs is < 0.05. Meanwhile, the Q/S ratio in near Earth asteroids (NEAs) has been estimated to be 0.5 to 2 (Binzel et al., 2004; Dandy et al., 2003). Therefore, Q-type NEAs might be delivered from the main belt region with weathered, S-type surface into near Earth region and then obtain their Q-type, non- weathered surface after undergoing re-surfacing process there. The resur- facing mechanisms could be: 1. dispersal of surface material by tidal effect during planetary encounters (Binzel et al., 2010; Nesvorny et al., 2010), 2. the YORP spin-up induced rotational-fission (Polishook et al., 2014) or surface re-arrangement, or 3. thermal degradation (Delbo et al., 2014).

The ultra-long period Cepheids (ULPCs) are classical Cepheids with pulsation periods exceeding $\approx 80$ days. The intrinsic brightness of ULPCs are ~1 to ~3 mag brighter than their shorter period counterparts. This makes them attractive in future distance scale work to derive distances beyond the limit set by the shorter period Cepheids. We have initiated a program to search for ULPCs in M31, using the single-band data taken from the Palomar Transient Factory, and identified eight possible candidates. In this work, we presented the VI-band follow-up observations of these eight candidates. Based on our VI-band light curves of these candidates and their locations in the color-magnitude diagram and the Period-Wesenheit diagram, we verify two candidates as being truly ULPCs. The six other candidates are most likely other kinds of long-period variables. With the two confirmed M31 ULPCs, we tested the applicability of ULPCs in distance scale work by deriving the distance modulus of M31. It was found to be $\mu_{M31,ULPC}=24.30\pm0.76$ mag. The large error in the derived distance modulus, together with the large intrinsic dispersion of the Period-Wesenheit (PW) relation and the small number of ULPCs in a given host galaxy, means that the question of the suitability of ULPCs as standard candles is still open. Further work is needed to enlarge the sample of calibrating ULPCs and reduce the intrinsic dispersion of the PW relation before re-considering ULPCs as suitable distance indicators.

We present Be star candidates in the open cluster NGC\,663, identified by H$\alpha$ imaging photometry with the Palomar Transient Factory Survey, as a pilot program to investigate how the Be star phenomena, the emission spectra, extended circumstellar envelopes, and fast rotation, correlate with massive stellar evolution. Stellar membership of the candidates was verified by 2MASS magnitudes and colors, and by proper motions. We discover 4 new Be stars and exclude one known Be star from being a member due to its inconsistent proper motions. The fraction of Be stars to member stars [N(Be)/N(members)] in NGC\,663 is 3.5\%. The spectral type of the 34 Be stars in NGC\,663 shows bimodal peaks at B0--B2 and B5--B7, which is consistent with the statistics in most star clusters. Additionally, we also discover 23 emission-line stars of different types, including non-member Be stars, dwarfs, and giants.

Asteroids of size larger than 150 m generally do not have rotation periods smaller than 2.2 hours. This spin cutoff is believed to be due to the gravitationally bound rubble-pile structures of the asteroids. Rotation with periods exceeding this critical value will cause asteroid breakup. Up until now, only one object, 2001 OE84, has been found to be an exception to this spin cutoff. We report the discovery of a new super-fast rotator, (335433) 2005 UW163, spinning with a period of 1.290 hours and a lightcurve variation of $r'\sim0.8$ mag from the observations made at the P48 telescope and the P200 telescope of the Palomar Observatory. Its $H_{r'} = 17.69 \pm 0.27$ mag and multi-band colors (i.e., $g'-r' = 0.68\pm0.03$ mag, $r'-i' = 0.19\pm0.02$ mag and SDSS $i-z = -0.45$ mag) show it is a V-type asteroid with a diameter of $0.6 +0.3/-0.2$ km. This indicates (335433) 2005 UW163 is a super-fast rotator beyond the regime of the small monolithic asteroid.

A new asteroid rotation period survey have been carried out by using the Palomar Transient Factory (PTF). Twelve consecutive PTF fields, which covered an area of 87 deg$^2$ in the ecliptic plane, were observed in $R$ band with a cadence of $\sim$20 min during February 15--18, 2013. We detected 2500 known asteroids with a diameter range of 0.5 km $\leq D \leq$ 200 km. Of these, 313 objects had highly reliable rotation periods and exhibited the "spin barrier" at $\sim2$ hours. In contrast to the flat spin rate distribution of the asteroids with 3 km $\leq D \leq$ 15 km shown by Pravec et al. (2008), our results deviated somewhat from a Maxwellian distribution and showed a decrease at the spin rate greater than 5 rev/day. One super-fast-rotator candidate and two possible binary asteroids were also found in this work.

P/2011 S1 (Gibbs) is an outer solar system comet or active Centaur with a similar orbit to that of the famous 29P/Schwassmann-Wachmann 1. P/2011 S1 (Gibbs) has been observed by the Pan-STARRS 1 (PS1) sky survey from 2010 to 2012. The resulting data allow us to perform multi-color studies of the nucleus and coma of the comet. Analysis of PS1 images reveals that P/2011 S1 (Gibbs) has a small nucleus $< 4$ km radius, with colors $g_{P1}-r_{P1} = 0.5 \pm 0.02$, $r_{P1}-i_{P1} = 0.12 \pm 0.02$ and $i_{P1}-z_{P1} = 0.46 \pm 0.03$. The comet remained active from 2010 to 2012, with a model-dependent mass-loss rate of $\sim100$ kg s$^{-1}$. The mass-loss rate per unit surface area of P/2011 S1 (Gibbs) is as high as that of 29P/Schwassmann-Wachmann 1, making it one of the most active Centaurs. The mass-loss rate also varies with time from $\sim 40$ kg s$^{-1}$ to 150 kg s$^{-1}$. Due to its rather circular orbit, we propose that P/2011 S1 (Gibbs) has 29P/Schwassmann-Wachmann 1-like outbursts that control the outgassing rate. The results indicate that it may have a similar surface composition to that of 29P/Schwassmann-Wachmann 1. Our numerical simulations show that the future orbital evolution of P/2011 S1 (Gibbs) is more similar to that of the main population of Centaurs than to that of 29P/Schwassmann-Wachmann 1. The results also demonstrate that P/2011 S1 (Gibbs) is dynamically unstable and can only remain near its current orbit for roughly a thousand years.

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.

Low-inclination Near-Earth Asteroid (297274) 1996 SK, which is also classified as a Potentially Hazardous Asteroid, has a highly eccentric orbit. It was studied by multi-wavelength photometry within the framework of an NEA color survey at the Lulin Observatory. We report here the finding of large color variation across the surface of (297274) 1996 SK within one asteroidal rotation period of $4.656\pm0.122$ hours and classify it as a S-type asteroid according to its average colors of $B-V=0.767\pm0.033$, $V-R=0.482\pm0.021$, $V-I=0.801\pm0.025$ and the corresponding relative reflectance spectrum. It might be indicative of differential space weathering effect or compositional inhomogeneity of the surface materials.

Ultra-long-period Cepheids (ULPCs) are important in distance-scale studies due to their potential for determining distance beyond ~100 Mpc. We performed a comprehensive search for ULPCs in M31, a local benchmark to calibrate the distance ladders. We use data from the Palomar Transient Factory (PTF), which has imaged M31 using a 1.2-m telescope equipped with a ~7.26 deg2 field-of-view (FOV) camera, usually with daily sampling, since the beginning of 2010. The large FOV, together with the regular monitoring, enables us to probe ULPCs in the bulge, disk, and even out to the halo of M31. Using a difference imaging analysis technique, we found and characterized 3 promising ULPC candidates based on their luminosities, amplitudes and Fourier parameters. The mean absolute magnitude for these 3 ULPC candidates, calibrated with latest M31 distance, is M_R=-6.47mag. Two out of the 3 ULPC candidates have been reported in literature, however their published periods from Magnier et al. are about half of the periods we found in this work. The third ULPC candidate is a new discovery. We studied 5 other candidates and determined that they are probably Mira-like or ultra-long-period variables, but not ULPCs.

We report the discovery of 2010 GB$_{174}$, a likely new member of the Inner Oort Cloud (IOC). 2010 GB$_{174}$ is one of 91 Trans Neptunian Objects (TNOs) and Centaurs discovered in a 76 deg$^2$ contiguous region imaged as part of the Next Generation Virgo Cluster Survey (NGVS) --- a moderate ecliptic latitude survey reaching a mean limiting magnitude of $g^\prime \simeq 25.5$ --- using MegaPrime on the 3.6m Canada France Hawaii Telescope. 2010 GB$_{174}$ is found to have an orbit with semi-major axis $a\simeq350.8$ AU, inclination $i \simeq 21.6^\circ$ and pericentre $q\sim48.5$ AU. This is the second largest perihelion distance among known solar system objects. Based on the sky coverage and depth of the NGVS, we estimate the number of IOC members with sizes larger than 300 km ($H_V \le 6.2$ mag) to be $\simeq 11\,000$. A comparison of the detection rate from the NGVS and the PDSSS (a characterized survey that `re-discovered' the IOC object Sedna) gives, for an assumed a power-law LF for IOC objects, a slope of $\alpha \simeq 0.7 \pm 0.2$, with only two detections in this region this slope estimate is highly uncertain.

We reported the monitoring results on spectrophotometry, photometry and imaging of comet 103P/Hartley 2 obtained at Lulin (1m), Calar Alto (2.2m) and Beijing Astronomical (2.16m) Observatory from April to December 2010. We found that a dust feature at sunward direction was detected starting from the end of September until the beginning of December (our last observation from the Lulin and Calar Alto observatory). Two distinct sunward jet features in the processed images were observed on October 11 and after October 29 until November 2. In parallel, the CN images reveal two asymmetrical jet features which are nearly perpendicular to the Sun-nucleus direction and this asymmetrical features implies that the comet was in a nearly side-on view in late-October and early-November. Additional to the jet features, the average result of the C2-to-CN production rate ratio ranges from 0.7 to 1.5 which places 103P/Hartley 2 as being of typical cometary chemistry. We found that the $r_h$ dependence for the dust production rate, Afrho (5,000 km), is -3.75+-0.45 before perihelion and is -3.44 +-1.20 during post-perihelion period. We detected the higher dust reddening is around the optocenter and getting bluer outward along the sunward jet feature and concluded that the former one, higher dust reddening, could be associated with strong jet activity and the latter one, the lowering of the reddening, might imply the optical properties changed or could be associated with outburst. The average dust color did not appear to vary significantly as the comet passed through perihelion.

The quasi-Hilda object, 212P/2000YN30 with a cometary-like orbit, was found to display a dust tail structure between January and March, 2009. From orbital calculations, it is shown that this object could have been an active comet in its past history before being transported to the current orbital configuration in quasi-stable 3:2 resonance with Jupiter.

We present a detailed description of the interrelation between the Type III radio bursts and energetic phenomena associated with the flare activities in Active region AR 11158 at 07:58 UT on 2011, Feb. 15. The timing of the Type-III radio burst measured by the radio wave experiment on the Wind/WAVE and an array of ground-based radio telescopes, coincided with an EUV jet and hard X-ray emission observed by SDO/AIA and RHESSI., respectively. There is clear evidence that the EUV jet shares the same source region as the hard X-ray emission. The temperature of the jet, as determined by multiwavelength measurements of AIA, suggests that type III emission is associated with hot, 7 MK, plasma at the jet's footpoint.

The Voyager flyby observations revealed that a very broad doughnut shaped distribution of the hydrogen atoms existed in the Saturnian magnetosphere. Recent Cassini observations confirmed the local-time asymmetry but also showed the hydrogen cloud density increases with decreasing distance to Saturn. The origin of the atomic hydrogen cloud has been debated ever since. Therefore, we have carried out a global investigation of the atomic hydrogen cloud taking into account all possible sources: 1) the Saturnian atmosphere, 2) the H2 atmosphere of main rings, 3) Enceladus H2O and OH torus, 4) Titan H2 torus and 5) the atomic hydrogen directly escaping from Titan. We show that the H ejection velocity and angle distribution are modified by collisions of the hot H, produced by electron-impact dissociation of H2, with the ambient atmospheric H2 and H. This in turn affects the morphology of the escaping hydrogen as does the morphology of the ionospheric electron distribution. That Saturn atmosphere is an important source is suggested by the fact that the H cloud peaks well below the ring plane, a feature that, so far, we can not reproduce by the dissociation of the ring H2 atmosphere or other proposed sources. Our simulations show that H directly escaping from Titan is a major contribution in the outer magnetosphere. The morphology of Titan H torus, shaped by the solar radiation pressure and the Saturnian oblateness, can account for the local time asymmetry near Titan orbit. Dissociation of H2O and OH in the Enceladus torus contributes inside ~5 RS, but dissociation of Titan H2 torus does not due to the significant energy released. The total number of H observed by Cassini inside 5 RS: our modeling results suggest ~20% from dissociation in the Enceladus torus, ~10% from dissociation of ring H2 atmosphere, and ~50% from Titan H torus implying that ~20% comes from the Saturnian atmosphere.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

We present the study on GRB 071112C X-ray and optical light curves. In these two wavelength ranges, we have found different temporal properties. The R-band light curve showed an initial rise followed by a single power-law decay, while the X-ray light curve was described by a single power-law decay plus a flare-like feature. Our analysis shows that the observed temporal evolution cannot be described by the external shock model in which the X-ray and optical emission are produced by the same emission mechanism. No significant color changes in multi-band light curves and a reasonable value of the initial Lorentz factor ({\Gamma}0 = 275 \pm 20) in a uniform ISM support the afterglow onset scenario as the correct interpretation for the early R-band rise. The result suggests the optical flux is dominated by afterglow. Our further investigations show that the X-ray flux could be created by an additional feature related to energy injection and X-ray afterglow. Different theoretical interpretations indicate the additional feature in X-ray can be explained by either late internal dissipation or local inverse-Compton scattering in the external shock.

The {\it Fermi} Large Area Telescope has recently discovered two giant gamma-ray bubbles which extend north and south of the Galactic center with diameters and heights of the order of $H\sim 10$ kpc. We suggest that the periodic star capture processes by the Galactic supermassive black hole Sgr A$^*$, with a capture rate of $\tau_{\rm cap}^{-1}\sim 3\times 10^{-5}$ yr$^{-1}$ and an energy release of $W\sim 3\times 10^{52}$ erg per capture, can result in hot plasma injecting into the Galactic halo at a wind velocity of $u\sim 10^8$ cm s$^{-1}$. The periodic injection of hot plasma can produce a series of shocks. Energetic protons in the bubble are re-accelerated when they interact with these shocks. We show that for energy larger than $E> 10^{15}$ eV, the acceleration process can be better described by the stochastic second-order Fermi acceleration. We propose that hadronic cosmic rays (CRs) within the ``knee'' of the observed CR spectrum are produced by Galactic supernova remnants distributed in the Galactic disk. Re-acceleration of these particles in the Fermi Bubble produces CRs beyond the knee. With a mean CR diffusion coefficient in this energy range in the bubble $D_B\sim 3\times 10^{30}$ cm$^2$ s$^{-1}$, we can reproduce the spectral index of the spectrum beyond the knee and within. The conversion efficiency from shock energy of the bubble into CR energy is about 10\%. This model provides a natural explanation of the observed CR flux, spectral indices, and matching of spectra at the knee.

Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10 kpc in diameter. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr A*, with a capture rate $<10^{-5}$ yr$^{-1}$ and energy release $\sim 10^{52}$ erg per one capture can produce shocks in the halo, which accelerate electrons to the energy ~ 1 TeV. These electrons generate radio emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic and the galactic soft photons. Estimates of the diffusion coefficient from the observed gamma-ray flux explains consistently the necessary maximum energy of electrons and sharp edges of the bubble.

Fermi LAT has discovered two extended gamma-ray bubbles above and below the galactic plane. We propose that their origin is due to the energy release in the Galactic center (GC) as a result of quasi-periodic star accretion onto the central black hole. Shocks generated by these processes propagate into the Galactic halo and accelerate particles there. We show that electrons accelerated up to ~10 TeV may be responsible for the observed gamma-ray emission of the bubbles as a result of inverse Compton (IC) scattering on the relic photons. We also suggest that the Bubble could generate the flux of CR protons at energies > 10^15 eV because the shocks in the Bubble have much larger length scales and longer lifetimes in comparison with those in SNRs. This may explain the the CR spectrum above the knee.

Fermi has discovered two giant gamma-ray-emitting bubbles that extend nearly 10kpc in diameter north and south of the galactic center (GC). The existence of the bubbles was first evidenced in X-rays detected by ROSAT and later WMAP detected an excess of radio signals at the location of the gamma-ray bubbles. We propose that periodic star capture processes by the galactic supermassive black hole, Sgr A$^*$, with a capture rate $3\times 10^{-5}$yr$^{-1}$ and energy release $\sim 3\times 10^{52}$erg per capture can produce very hot plasma $\sim 10$keV with a wind velocity $\sim 10^8$cm/s injected into the halo and heat up the halo gas to $\sim 1$keV, which produces thermal X-rays. The periodic injection of hot plasma can produce shocks in the halo and accelerate electrons to $\sim$TeV, which produce radio emission via synchrotron radiation, and gamma-rays via inverse Compton scattering with the relic and the galactic soft photons.

We report here the discovery of an new near-equal mass Trans-Neptunian Binaries (TNBs) L5c02 and the and the putative detection of a second TNB (L4k12) among the year two and three detections of the Canada-France-Eclipic Plane Survey (CFEPS). These new binaries (internal designation L4k12 and L5c02) have moderate separations of 0.4" and 0.6" respectively. The follow-up observation confirmed the binarity of L5c02, but L4k12 are still lack of more followup observations. L4k12 has a heliocentric orbital inclination of ~ 35?degree, marking this system as having the highest heliocentric orbital inclination among known near-equal mass binaries. Both systems are members of the classical main Kuiper belt population. Based on the sample of objects searched we determine that the fraction of near-equal mass wide binaries with separations > 0.4" is 1.5% to 20% in the cold main classical Kuiper belt and, if our detection the binarity L4k12 holds, 3% to 43% in the hot main classical objects are binary. In this manuscript we describe our detection process, the sample of objects surveyed, our confirmation observations.

We discuss the origin of continuum and line X-ray emission observed in the direction the Galactic Center. We predict a significant flux of de-excitation gamma-ray lines in this direction, which can be produced by subrelativistic protons generated by accretion processes.

The origin and properties of the source of positrons annihilating in the Galactic Center is still a mystery. One of the criterion, which may discriminate between different mechanisms of positron production there, is the positron energy injection. Beacom and Yueksel (2006) suggested a method to estimate this energy from the ratio of the 511 keV line to the MeV in-flight annihilation fluxes. From the COMPTEL data they derived that the maximum injection energy of positron should be about several MeV that cut down significantly a class of models of positron origin in the GC assuming that positrons lose their energy by Coulomb collisions only. However, observations show that the strength of magnetic field in the GC is much higher than in other parts of the Galaxy, and it may range there from 100 $\mu$G to several mG. In these conditions, synchrotron losses of positrons are significant that extends the range of acceptable values of positron injection energy. We show that if positrons injection in the GC is non-stationary and magnetic field is higher than 0.4 mG both radio and gamma-ray restrictions permit their energy to be higher than several GeV.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

We present multi-band results for GRB071010B based on \Swift, \Suzaku, and ground-based optical observations. This burst is an ideal target to evaluate the robustness of the ${\rm E^{src}_{peak}-E_{iso}}$ and ${\rm E^{src}_{peak}-E_{\gamma}}$ relations, whose studies have been in stagnation due to the lack of the combined estimation of $\rm E^{src}_{peak}$ and long term optical monitoring. The joint prompt spectral fitting using \Swift/Burst Alert Telescope and \Suzaku/Wide-band All sky Monitor data yielded the spectral peak energy as E$^{src}_{peak}$ of $86.5^{+6.4}_{-6.3}$ keV and E$_{iso}$ of $2.25^{+0.19}_{-0.16}\times10^{52}$ erg with $z=0.947$. The optical afterglow light curve is well fitted by a simple power law with temporal index $\alpha=-0.60\pm0.02$. The lower limit of temporal break in the optical light curve is 9.8 days. Our multi-wavelength analysis reveals that GRB071010B follows ${\rm E^{src}_{peak}-E_{iso}}$ but violates the ${\rm E^{src}_{peak}-E_{\gamma}}$ and ${\rm E_{iso}-E^{src}_{peak}-t^{src}_{jet}}$ at more than the 3$\sigma$ level.

Aims. The accretion of stars onto the central supermassive black hole at the center of the Milky Way is predicted to generate large fluxes of subrelativistic ions in the Galactic center region. We analyze the intensity, shape and spatial distribution of de-excitation gamma-ray lines produced by nuclear interactions of these energetic particles with the ambient medium. Methods. We first estimate the amount and mean kinetic energy of particles released from the central black hole during star disruption. We then calculate from a kinetic equation the energy and spatial distributions of these particles in the Galactic center region. These particle distributions are then used to derive the characteristics of the main nuclear interaction gamma-ray lines. Results. Because the time period of star capture by the supermassive black hole is expected to be shorter than the lifetime of the ejected fast particles against Coulomb losses, the gamma-ray emission is predicted to be stationary. We find that the nuclear de-excitation lines should be emitted from a region of maximum 5$^\circ$ angular radius. The total gamma-ray line flux below 8 MeV is calculated to be $\approx10^{-4}$ photons cm$^{-2}$ s$^{-1}$. The most promising lines for detection are those at 4.44 and $\sim$6.2 MeV, with a predicted flux in each line of $\approx$$10^{-5}$ photons cm$^{-2}$ s$^{-1}$. Unfortunately, it is unlikely that this emission can be detected with the INTEGRAL observatory. But the predicted line intensities appear to be within reach of future gamma-ray space instruments. A future detection of de-excitation gamma-ray lines from the Galactic center region would provide unique information on the high-energy processes induced by the central supermassive black hole and the physical conditions of the emitting region.

Here we present initial results of our search for extensive halos around the planetary nebulae (PNe) in our Galaxy based on imaging data from the Southern H-Alpha Sky Survey Atlas (SHASSA). A threshold surface brightness in H$\alpha$ emission was used to help identify the spatially extended features of the PNe. We investigated a sample that included 10 PNe, the large majority of which were found to be surrounded by extensive outer halos. The formation of these halos might be closely related to the AGB phase mass loss and/or the interaction of the stellar outflows with the interstellar medium. Most of these outer halos are nearly spherical. Close investigation of some specific objects, such as He 2-111 and NGC 3242, indicate a kinematic age on the order of 10$^{5}$ yrs. The mass loss history can be traced back to the late AGB phase of the evolution of the progenitors. Two objects form the sample (He 2-111, NGC 3242), were found to have outer halos with fragmented arcs that are apparently the result of their interaction with the interstellar medium.

We study the ``normal'' decay phase of the X-ray afterglows of gamma-ray bursts (GRBs), which follows the shallow decay phase, using the events simultaneously observed in the R-band. The classical external shock model -- in which neither the delayed energy injection nor time-dependency of shock micro-physics is considered -- shows that the decay indices of the X-ray and R-band light curves, $\alpha_{\rm X}$ and $\alpha_{\rm O}$, obey a certain relation, and that in particular, $\alpha_{\rm O}-\alpha_{\rm X}$ should be larger than -1/4 unless the ambient density increases with the distance from the central engine. For our selected 14 samples, we have found that 4 events violate the limit at more than the 3$\sigma$ level, so that a fraction of events are outliers of the classical external shock model at the ``normal'' decay phase.

Observations were made of the optical afterglow of GRB 051028 with the Lulin observatory's 1.0 m telescope and the WIDGET robotic telescope system. R band photometric data points were obtained on 2005 October 28 (UT), or 0.095-0.180 days after the burst. There is a possible plateau in the optical light curve around 0.1 days after the burst; the light curve resembles optically bright afterglows (e.g. GRB 041006, GRB 050319, GRB060605) in shape of the light curve but not in brightness. The brightness of the GRB 051028 afterglow is 3 magnitudes fainter than that of one of the dark events, GRB 020124. Optically dark GRBs have been attributed to dust extinction within the host galaxy or high redshift. However, the spectrum analysis of the X-rays implies that there is no significant absorption by the host galaxy. Furthermore, according to the model theoretical calculation of the Ly$\alpha$ absorption to find the limit of GRB 051028's redshift, the expected $R$ band absorption is not high enough to explain the darkness of the afterglow. While the present results disfavor either the high-redshift hypothesis or the high extinction scenario for optically dark bursts, they are consistent with the possibility that the brightness of the optical afterglow, intrinsically dark.

We retrieve the H_2O expansion velocity in a number of comets, using the 18-cm line shapes of the OH radical observed with the Nan\c{c}ay radio telescope. The H_2O velocity is derived from the large base of a trapezium fitted to the observed spectra. This method, which was previously applied to 9 comets, is now extended to 30 further comets. This allows us to study the evolution of their water molecule outflow velocity over a large range of heliocentric distances and gas production rates. Our analysis confirms and extends previous analyses. The retrieved expansion velocities increases with increasing gas production rates and decreasing heliocentric distances. Heuristic laws are proposed, which could be used for the interpretation of observations of cometary molecules and as a touchstone for hydrodynamical models. The expansion velocities retrieved from 18 cm line shapes are larger than those obtained from millimetric observations of parent molecules with smaller fields of view, which demonstrates the acceleration of the gas with cometocentric distance. Our results are in reasonable quantitative agreement with current hydrodynamical models of cometary atmospheres.

Observations of the optical afterglow of GRB 041006 with the Kiso Observatory 1.05 m Schmidt telescope, the Lulin Observatory 1.0 m telescope and the Xinglong Observatory 0.6 m telescope. Three-bands (B, V and R) of photometric data points were obtained on 2004 October 6, 0.025-0.329 days after the burst. These very early multi band light curves imply the existence of a color dependent plateau phase. The B-band light curve shows a clear plateau at around 0.03 days after the burst. The R band light curve shows the hint of a plateau, or a possible slope change, at around 0.1 days after the burst. The overall behavior of these multi-band light curves may be interpreted in terms of the sum of two separate components, one showing a monotonic decay the other exhibiting a rising and a falling phase, as described by the standard afterglow model.

Aims. Understand the shape and implications of the multiband light curve of GRB 050408, an X-ray rich (XRR) burst. Methods. We present a multiband optical light curve, covering the time from the onset of the gamma-ray event to several months after, when we only detect the host galaxy. Together with X-ray, millimetre and radio observations we compile what, to our knowledge, is the most complete multiband coverage of an XRR burst afterglow to date. Results. The optical and X-ray light curve is characterised by an early flattening and an intense bump peaking around 6 days after the burst onset. We explain the former by an off-axis viewed jet, in agreement with the predictions made for XRR by some models, and the latter with an energy injection equivalent in intensity to the initial shock. The analysis of the spectral flux distribution reveals an extinction compatible with a low chemical enrichment surrounding the burst. Together with the detection of an underlying starburst host galaxy we can strengthen the link between XRR and classical long-duration bursts.

Multi-wavelength B, V, R, I observations of the optical afterglow of GRB 050319 were performed by the 1.05-m telescope at Kiso Observatory and the 1.0-m telescope at Lulin Observatory from 1.31 hours to 9.92 hours after the burst. Our R band lightcurves, combined with other published data, can be described by the smooth broken power-law function, with $\alpha_1$ = -0.84 $\pm$0.02 to $\alpha_2$ = -0.48$\pm$0.03, 0.04 days after the GRB. The optical lightcurves are characterized by shallow decays-- as was also observed in the X-rays-- which may have a similar origin, related to energy injection. However, our observations indicate that there is still a puzzle concerning the chromatic breaks in the R band lightcurve (at 0.04 days) and the X-ray lightcurve (at 0.004 days) that remains to be solved.

Even with the renaissance in gamma-ray burst (GRB) research fostered by the Swift satellite, few bursts have both contemporaneous observations at long wavelengths and exquisite observations at later times across the electromagnetic spectrum. We present here contemporaneous imaging with the KAIT robotic optical telescope, dense optical sampling with Lulin, and supplemented with infrared data from PAIRITEL and radio to gamma-ray data from the literature. For the first time, we can test the constancy of microphysical parameters in the internal-external shock paradigm and carefully trace the flow of energy from the GRB to the surrounding medium. KAIT data taken <~1 minute after the start of GRB051111 and coinciding with the fading gamma-ray tail of the prompt emission indicate a smooth re-injection of energy into the shock. No color change is apparent in observations beginning ~1.5 minutes after the GRB and lasting for the first hour after the burst. There are achromatic flux modulations about the best-fit model at late (t~10^4 s) times, possibly due to variations in the external density. We find that the host-galaxy extinction is well fit by a curve similar to that ofthe Small Magellanic Cloud. Low visual extinction, A_V~0.2 mag, combined with high column densities determined from the X-ray and optical spectroscopy (N_H> 10^21 cm^-2), indicate a low dust-to-metals ratio and a possible over-abundance of the light metals. An apparent small ratio of total to selective extinction (R_V~2) argues against dust destruction by the GRB. Time constancy of both the IR/optical/UV spectral energy distribution and the soft X-ray absorption suggests that the absorbing material is not local to the GRB.

The stochastic acceleration of subrelativistic electrons from a background plasma is studied in order to find a possible explanation of the hard X-ray (HXR) emission detected from the Coma cluster. We calculate the necessary energy supply as a function of the plasma temperature and of the electron energy. We show that, for the same value of the HXR flux, the energy supply changes gradually from its high value (when emitting particle are non-thermal) to lower values (when the electrons are thermal). The kinetic equations we use include terms describing particle thermalization as well as momentum diffusion due to the Fermi II acceleration. We show that the temporal evolution of the particle distribution function has, at its final stationary stage, a rather specific form: it cannot be described by simple exponential or power-law expressions. A broad transfer region is formed by Coulomb collisions at energies between the Maxwellian and power-law parts of the distribution function. In this region the radiative lifetime of a single electron differs greatly from the lifetime of the distribution function as a whole. For a plasma temperature of 8 keV, the particles emitting bremsstrahlung at 20-80 keV lie in this quasi-thermal regime. We show that the energy supply required by quasi-thermal electrons to produce the observed HXR flux from Coma is one or two orders of magnitude smaller than the value derived from the assumption of a nonthermal origin of the emitting particles. This result may solve the problem of rapid cluster overheating by nonthermal electrons. We finally predict the change in Coma's SZ effect caused by the distortions of the Maxwellian electron spectrum, and we show that evidence for acceleration of subrelativistic electrons can be derived from detailed spectral measurements.

The young planetary nebulae play an important role in stellar evolution when intermediate- to low-mass stars (0.8 $\sim$ 8 M$_\odot$) evolve from the proto-planetary nebulae phase to the planetary nebulae phase. Many young planetary nebulae display distinct bipolar structures as they evolve away from the proto-planetary nebulae phase. One possible cause of their bipolarity could be due to a binary origin of its energy source. Here we report our detailed investigation of the young planetary nebula, Hubble 12, which is well-known for its extended hourglass-like envelope. We present evidence with time-series photometric observations the existence of an eclipsing binary at the center of Hubble 12. Low-resolution spectra of the central source show, on the other hand, absorption features such as CN, G-band & Mg b${\arcsec}$, which can be suggestive of a low-mass nature of the secondary component.

The 1-m telescope at Lulin Observatory and the 0.76-m Katzman Automatic Imaging Telescope at Lick Observatory were used to observe the optical afterglow of the short-duration (1.2--1.5 s) gamma-ray burst (GRB) 040924. This object has a soft high-energy spectrum, thus making it an exceptional case, perhaps actually belonging to the short-duration tail of the long-duration GRBs. Our data, combined with other reported measurements, show that the early R-band light curve can be described by two power laws with index alpha = -0.7 (at t = 16-50 min) and alpha = -1.06 (at later times). The rather small difference in the spectral indices can be more easily explained by an afterglow model invoking a cooling break rather than a jet break.

In 2004, we established a Japan-Taiwan-China collaboration for GRB study in the East-Asian region. This serves as a valuable addition to the world-wide optical and infrared follow-up network, because the East-Asia region would otherwise be blank. We have been carrying out imaging and spectroscopic follow-up observations at Lulin (Taiwan), Kiso (Japan), WIDGET (Japan) and Xinglong (China). From Xinglong and Kiso, we can locate candidates and obtain early time spectra for afterglows. While WIDGET provides early time observations before the burst, the high-time resolution for multi-band light curves can be obtained at Lulin. With the data from these sites, we can obtain detailed information about the light curve and redshift of GRBs, which are important to understand the mechanism of the afterglows. Up to March 2005, ten follow-up observations have been provided by this East-Asia cooperation. Two optical afterglows were detected, GRB 040924 and GRB 041006. The results of the two detected afterglows are reported in this article.

The Lulin GRB program, using the Lulin One-meter Telescope (LOT) in Taiwan started in July 2003. Its scientific aims are to discover optical counterparts of XRFs and short and long GRBs, then to quickly observe them in multiple bands. Thirteen follow-up observations were provided by LOT between July 2003 and Feb. 2005. One host galaxy was found at GRB 031203. Two optical afterglows were detected for GRB 040924 and GRB 041006. In addition, the optical observations of GRB 031203 and a discussion of the non-detection of the optical afterglow of GRB 031203 are also reported in this article.

We investigate the evolution and number distribution of radio halos in galaxy clusters. Without re-acceleration or regeneration, the relativistic electrons responsible for the diffuse radio emission will lose their energy via inverse-Compton and synchrotron losses in a rather short time, and radio halos will have lifetimes $\sim$ 0.1 Gyr. Radio halos could last for $\sim$ Gyr if a significant level of re-acceleration is involved. The lifetimes of radio halos would be comparable with the cosmological time if the radio-emitting electrons are mainly the secondary electrons generated by pion decay following proton-proton collisions between cosmic-ray protons and the thermal intra-cluster medium within the galaxy clusters. Adopting both observational and theoretical constraints for the formation of radio halos, we calculate the formation rates and the comoving number density of radio halos in the hierarchical clustering scheme. Comparing with observations, we find that the lifetimes of radio halos are $\sim$ Gyr. Our results indicate that a significant level of re-acceleration is necessary for the observed radio halos and the secondary electrons may not be a dominant origin for radio halos.

We investigate theoretical models for the radio halo and hard X-ray (HXR) excess in the Coma galaxy cluster. Time-independent and time-dependent re-acceleration models for relativistic electrons have been carried out to study the formation of the radio halo and HXR excess. In these models, the relativistic electrons are injected by merger shocks and re-accelerated by ensuing violent turbulence. The effects of different Mach numbers of the merger shocks on the radio and HXR excess emission are also investigated. We adopt 6 $\mu$G as the central magnetic field and reproduce the observed radio spectra via the synchrotron emission. We also obtain a central "plateau" in the radio spectral-index distribution, which have been observed in radio emission distribution. Our models can also produce the observed HXR excess emission via the inverse Compton scattering of the cosmic microwave background photons. We find that only the merger shocks with the Mach numbers around 1.6--2 can produce results in agreement with both the radio and HXR emission in the Coma cluster.

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

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

The role of cosmic rays (CR) in the formation and evolution of clusters of galaxies has been much debated. It may well be related to other fundamental questions, such as the mechanism which heats and virializes the intracluster medium (ICM), and the frequency at which the ICM is shocked. There is now compelling evidence both from the cluster soft excess (CSE) and the `hard-tail' emissions at energies above 10 keV, that many clusters are luminous sources of inverse-Compton (IC) emission. This is the first direct measurement of cluster CR: the technique is free from our uncertainties in the ICM magnetic field, and is not limited to the small subset of clusters which exhibit radio halos. The CSE emitting electrons fall within a crucial decade of energy where they have the least spectral evolution, and where most of the CR pressure resides. However their survival times do not date them back to the relic CR population. By using the CSE data of the Coma cluster, we demonstrate that the CR are energetically as important as the thermal ICM: the two components are in pressure equiparition. Thus, contrary to previous expectations, CR are a dominant component of the ICM, and their origin and effects should be explored. The best-fit CR spectral index is in agreement with the Galactic value.