Source author record

Wing-Huen Ip

Wing-Huen Ip appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

The 3D Direct Simulation Monte Carlo Study of Europa Gas Plume

Europa has been spotted to have water outgassing activities by the space and ground-based telescopes as well as reanalysis of the Galileo data (Roth et al. 2014; Sparks et al. 2016, 2017; Paganini et al. 2020; Jia et al. 2018; Arnold et al. 2019). However, these observations only provided limited information about plume dynamics, which is critical in understanding the eruption mechanism and preparation of future exploration. We adopt a 3D DSMC model to investigate the plume characteristics of Europa assuming supersonic expansion originated from the undersurface vent. The main goal is to understand the physical processes and structures of Europa water vapor plumes, which can play a key role on probing its undersurface vent condition and outgassing mechanism. With a parametric study of the total gas production rate and initial gas bulk velocity, the gas number density, temperature and velocity information of the outgassing plumes from the various case studies are derived. Our results show that the plume gases experience acceleration through mutual collisions and adiabatic cooling when exiting and expanding from the surface. The central part of the plume with the relatively large gas production rates (of 1029 and 1030 H2O s-1) is found to sustain thermal equilibrium and nearly continuum condition. Column density maps integrated along two different viewing angles are presented to demonstrate the importance of the projection effect on remote sensing diagnostics. Finally, the density profiles at different altitudes are provided to prepare for observations of Europa plumes including the upcoming spacecraft missions such as JUICE and Europa Clipper.

preprint2021arXiv

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years.

preprint2021arXiv

OSSOS: The Eccentricity and Inclination Distributions of the Stable Neptunian Trojans

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 $σ_{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 $σ_e$ of 0.045 and a typical sin(i) x Gaussian inclination distribution with a width $σ_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 $σ_{e-cold}$ ~ 0.02 / $σ_{i-cold}$ ~ 6$^\circ$ and $σ_{e-hot}$~ 0.05 / $σ_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 members 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.

preprint2020arXiv

Asteroid migration due to the Yarkovsky effect and the distribution of the Eos family

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $Δ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 $Δa$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $Δ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 $ω_{\rm rot}$ of asteroid depends on its size $R$ by $ω_{\rm rot}\propto R^{-1}$.

preprint2020arXiv

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections.

preprint2020arXiv

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

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)

preprint2020arXiv

Transit of asteroids across the 7/3 Kirkwood gap under the Yarkovsky effect

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.

preprint2019arXiv

A Twilight Search for Atiras, Vatiras and Co-orbital Asteroids: Preliminary Results

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.

preprint2019arXiv

Asteroid Discovery and Light Curve Extraction Using the Hough Transform -- A Rotation Period Study for Sub-Kilometer Main-Belt Asteroids

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.

preprint2019arXiv

Low activity main belt comet 133P/Elst-Pizarro: New constraints on its Albedo, Temperature and Active Mechanism from a thermophysical perspective

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.

preprint2016arXiv

Be Stars in the Open Cluster NGC 6830

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.

preprint2016arXiv

Discovery of A New Retrograde Trans-Neptunian Object: Hint of A Common Orbital Plane for Low Semi-Major Axis, High Inclination TNOs and Centaurs

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

preprint2016arXiv

Discovery of an X-ray Emitting Contact Binary System 2MASS J11201034$-$2201340

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

preprint2016arXiv

OSSOS: IV. Discovery of a dwarf planet candidate in the 9:2 resonance with Neptune

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.

preprint2016arXiv

Palomar Transient Factory and RR Lyrae: Metallicity-Light Curve Relation Based on ab-Type RR Lyrae in Kepler Field

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 ϕ_{31}$ (where $P$ is pulsational period and $ϕ_{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.

preprint2016arXiv

The Outer Solar System Origins Survey: I. Design and First-Quarter Discoveries

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.

preprint2016arXiv

The Pan-STARRS 1 Discoveries of five new Neptune Trojans

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.

preprint2015arXiv

A Quick Test on Rotation Period Clustering for the Small Members of the Koronis Family

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.

preprint2015arXiv

Asteroid lightcurves from the Palomar Transient Factory survey: Rotation periods and phase functions from sparse photometry

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.

preprint2015arXiv

Asteroid Spin-Rate Study using the Intermediate Palomar Transient Factory

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

preprint2015arXiv

Boulders on asteroid Toutatis as observed by Chang'e-2

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.

preprint2015arXiv

Orientation and Rotational Parameters of Asteroid 4179 Toutatis: New Insights from Chang'e-2's Close Flyby

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 $λ_{H}=180.2^{+0.2^\circ}_{-0.3^\circ}$ and $β_{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.

preprint2015arXiv

VI-Band Follow-Up Observations of Ultra-Long-Period Cepheid Candidates in M31

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 $μ_{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.

preprint2014arXiv

313 new asteroid rotation periods from Palomar Transient Factory observations

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.

preprint2014arXiv

A New Large Super-Fast Rotator: (335433) 2005 UW163

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.

preprint2013arXiv

Detection of Large Color Variation of Potentially Hazardous Asteroid (297274) 1996 SK

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.

preprint2013arXiv

Discovery of a New Member of the Inner Oort Cloud from The Next Generation Virgo Cluster Survey

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 $α\simeq 0.7 \pm 0.2$, with only two detections in this region this slope estimate is highly uncertain.

preprint2013arXiv

Flare-associated type III radio bursts and dynamics of the EUV jet from SDO/AIA and RHESSI observations

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.

preprint2013arXiv

Long-Term Monitoring of Comet 103P/Hartley 2

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.

preprint2013arXiv

Using the Palomar Transient Factory to Search for Ultra-Long-Period Cepheid Candidates in M31

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.

preprint2012arXiv

Discovery of Main-Belt Comet P/2006 VW139 by Pan-STARRS1

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.

preprint2009arXiv

Swifts GRB GRB071010B : outlier of the $\rm E^{src}_{peak} - E_γ$ and ${\rm E_{iso}-E^{src}_{peak}-t^{src}_{jet}}$ correlations

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_γ}$ 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 $α=-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_γ}$ and ${\rm E_{iso}-E^{src}_{peak}-t^{src}_{jet}}$ at more than the 3$σ$ level.

preprint2007arXiv

A multi band study of the optically dark GRB 051028

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