First EURONEAR NEA discoveries from La Palma using the INT

Since 2006, the European Near Earth Asteroids Research (EURONEAR) project has been contributing to the research of near-Earth asteroids (NEAs) within a European network. One of the main aims is the amelioration of the orbits of NEAs, and starting in 2014 February we focus on the recovery of one-opposition NEAs using the Isaac Newton Telescope (INT) in La Palma in override mode. Part of this NEA recovery project, since 2014 June EURONEAR serendipitously started to discover and secure the first NEAs from La Palma and using the INT, thanks to the teamwork including amateurs and students who promptly reduce the data, report discoveries and secure new objects recovered with the INT and few other telescopes from the EURONEAR network. Five NEAs were discovered with the INT, including 2014 LU14, 2014 NL52 (one very fast rotator), 2014 OL339 (the fourth known Earth quasi-satellite), 2014 SG143 (a quite large NEA), and 2014 VP. Another very fast moving NEA was discovered but was unfortunately lost due to lack of follow-up time. Additionally, another 14 NEA candidates were identified based on two models, all being rapidly followed-up using the INT and another 11 telescopes within the EURONEAR network. They include one object discovered by Pan-STARRS, two Mars crossers, two Hungarias, one Jupiter trojan, and other few inner main belt asteroids (MBAs). Using the INT and Sierra Nevada 1.5 m for photometry, then the Gran Telescopio de Canarias for spectroscopy, we derived the very rapid rotation of 2014 NL52, then its albedo, magnitude, size, and its spectral class. Based on the total sky coverage in dark conditions, we evaluate the actual survey discovery rate using 2-m class telescopes. One NEA is possible to be discovered randomly within minimum 2.8 deg2 and maximum 5.5 deg2. These findings update our past statistics, being based on double sky coverage and taking into account the recent increase in discovery.

Exploring S-type orbits in the Pluto-Charon binary system

This work generates, through a sample of numerical simulations of the restricted three-body problem, diagrams of semimajor axis and eccentricity which defines stable and unstable zones for particles in S-type orbits around Pluto and Charon. Since we consider initial conditions with 0 <= e <= 0.99, we found several new stable regions. We also identified the nature of each one of these newly found stable regions. They are all associated to families of periodic orbits derived from the planar circular restricted three-body problem. We have shown that a possible eccentricity of the Pluto-Charon system slightly reduces, but does not destroy, any of the stable regions.

Stability regions around the components of the triple system 2001 SN263

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Near-Earth asteroid binaries in close encounters with the Earth

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A multidomain approach to asteroid families' identification

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Stable retrograde orbits around the triple system 2001 SN263

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Absolute Magnitudes of Pan-STARRS PS1 Asteroids

Absolute magnitude (H) of an asteroid is a fundamental parameter describing the size and the apparent brightness of the body. Because of its surface shape, properties and changing illumination, the brightness changes with the geometry and is described by the phase function governed by the slope parameter (G). Although many years have been spent on detailed observations of individual asteroids to provide H and G, vast majority of minor planets have H based on assumed G and due to the input photometry from multiple sources the errors of these values are unknown. We compute H of ~ 180 000 and G of few thousands asteroids observed with the Pan-STARRS PS1 telescope in well defined photometric systems. The mean photometric error is 0.04 mag. Because on average there are only 7 detections per asteroid in our sample, we employed a Monte Carlo (MC) technique to generate clones simulating all possible rotation periods, amplitudes and colors of detected asteroids. Known asteroid colors were taken from the SDSS database. We used debiased spin and amplitude distributions dependent on size, spectral class distributions of asteroids dependent on semi-major axis and starting values of G from previous works. H and G (G12 respectively) were derived by phase functions by Bowell et al. (1989) and Muinonen et al. (2010). We confirmed that there is a positive systematic offset between H based on PS1 asteroids and Minor Planet Center database up to -0.3 mag peaking at 14. Similar offset was first mentioned in the analysis of SDSS asteroids and was believed to be solved by weighting and normalizing magnitudes by observatory codes. MC shows that there is only a negligible difference between Bowell's and Muinonen's solution of H. However, Muinonen's phase function provides smaller errors on H. We also derived G and G12 for thousands of asteroids. For known spectral classes, slope parameters agree with the previous work in general, however, the standard deviation of G in our sample is twice as larger, most likely due to sparse phase curve sampling. In the near future we plan to complete the H and G determination for all PS1 asteroids (500,000) and publish H and G values online. This work was supported by NASA grant No. NNX12AR65G.

1. Astronomical photometry at wavelengths of 8.5, 10.5, and 11.6 µm. 2. Infrared emission from asteroids at wavelengths of 8.5, 10.5, and 11.6 µm

The purpose of this thesis is to present new observations of thermal-infrared radiation from asteroids. Stellar photometry was performed to provide standards for comparison with the asteroid data. The details of the photometry and the data reduction are discussed in Part 1. A system of standard stars is derived for wavelengths of 8.5, 10.5 and 11.6 µm and a new calibration is adopted. Sources of error are evaluated and comparisons are made with the data of other observers.

The observations and analysis of the thermal-emission observations of asteroids are presented in Part 2. Thermal-emission lightcurve and phase effect data are considered. Special color diagrams are introduced to display the observational data. These diagrams are free of any model-dependent assumptions and show that asteroids differ in their surface properties.

On the basis of photometric models, (4) Vesta is thought to have a bolometric Bond albedo of about 0.1, an emissivity greater than 0.7 and a true radius that is close to the model value of 300^(+50)_(-30)km. Model albedos and model radii are given for asteroids 1, 2, 4, 5, 6, 7, 15, 19, 20, 27, 39, 44, 68, 80, 324 and 674. The asteroid (324) Bamberga is extremely dark with a model (~bolometric Bond) albedo in the 0.01 - 0.02 range, which is thought to be the lowest albedo yet measured for any solar-system body. The crucial question about such low-albedo asteroids is their number and the distribution of their orbits.

Solar System Small-Body Demographics with the Palomar Transient Factory Survey

Observational studies of our solar system's small-body populations (asteroids and comets) offer insight into the history of our planetary system, as these minor planets represent the left-over building blocks from its formation. The Palomar Transient Factory (PTF) survey began in 2009 as the latest wide-field sky-survey program to be conducted on the 1.2-meter Samuel Oschin telescope at Palomar Observatory. Though its main science program has been the discovery of high-energy extragalactic sources (such as supernovae), during its first five years PTF has collected nearly five million observations of over half a million unique solar system small bodies. This thesis begins to analyze this vast data set to address key population-level science topics, including: the detection rates of rare main-belt comets and small near-Earth asteroids, the spin and shape properties of asteroids as inferred from their lightcurves, the applicability of this visible light data to the interpretation of ultraviolet asteroid observations, and a comparison of the physical properties of main-belt and Jovian Trojan asteroids. Future sky-surveys would benefit from application of the analytical techniques presented herein, which include novel modeling methods and unique applications of machine-learning classification. The PTF asteroid small-body data produced in the course of this thesis work should remain a fertile source of solar system science and discovery for years to come.

Tight constraints on the existence of additional planets around HD 189733†

We report a transit timing study of the transiting exoplanetary system HD 189733. In total, we observed 10 transits in 2006 and 2008 with the 2.6-m Nordic Optical Telescope, and two transits in 2007 with the 4.2-m William Herschel Telescope. We used Markov Chain Monte Carlo simulations to derive the system parameters and their uncertainties, and our results are in a good agreement with previously published values. We performed two independent analyses of transit timing residuals to place upper mass limits on putative perturbing planets. The results show no evidence for the presence of planets down to 1 Earth mass near the 1:2 and 2:1 resonance orbits, and planets down to 2.2 Earth masses near the 3:5 and 5:3 resonance orbits with HD 189733b. These are the strongest limits to date on the presence of other planets in this system. Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by Denmark, Finland, Iceland, Norway and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. ‡

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Resonance and chaos .1. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted three-body problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincare surfaces of section with a mass ratio of 10(-3) (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. on theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Resonance and chaos: I. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted threebody problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincaré surfaces of section with a mass ratio of 10-3 (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. On theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

The Gemini planet-finding campaign: the frequency of giant planets aroud debris disk stars

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.''5 and 14.1 mag at 1'' separation. Follow-up observations of the 66 candidates with projected separation <500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known β Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a ≥5 M Jup planet beyond 80 AU, and <21% of debris disk stars have a ≥3 M Jup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly imaged planets as d 2 N/dMdavpropm α a β, where m is planet mass and a is orbital semi-major axis (with a maximum value of a max). We find that β < –0.8 and/or α > 1.7. Likewise, we find that β < –0.8 and/or a max < 200 AU. For the case where the planet frequency rises sharply with mass (α > 1.7), this occurs because all the planets detected to date have masses above 5 M Jup, but planets of lower mass could easily have been detected by our search. If we ignore the β Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that <20% of debris disk stars have a ≥3 M Jup planet beyond 10 AU, and β < –0.8 and/or α < –1.5. Likewise, β < –0.8 and/or a max < 125 AU. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation >40 AU and planet masses >3 M Jup do not carve the central holes in these disks.