Partial Evaluation for Scientific Computing: The Supercomputer Toolkit Experience

We describe the key role played by partial evaluation in the Supercomputer Toolkit, a parallel computing system for scientific applications that effectively exploits the vast amount of parallelism exposed by partial evaluation. The Supercomputer Toolkit parallel processor and its associated partial evaluation-based compiler have been used extensively by scientists at M.I.T., and have made possible recent results in astrophysics showing that the motion of the planets in our solar system is chaotically unstable.

Structure of glasses and melts

Wilding, M. C., Benmore, C. J. (2006). Structure and Glasses and Melts. reviews in Mineralogy and Geochemistry, 63 (1), 275-311 RAE2008

Badanie rotacji wybranych planetoid pojedynczych i podwójnych przy wykorzystaniu obserwacji fotometrycznych

Wydział Fizyki: Instytut Obserwatorium Astronomiczne

Resonant and Secular Orbital Interations

In stable solar systems, planets remain in nearly elliptical orbits around their stars. Over longer timescales, however, their orbital shapes and sizes change due to mutual gravitational perturbations. Orbits of satellites around a planet vary for the same reason. Because of their interactions, the orbits of planets and satellites today are different from what they were earlier. In order to determine their original orbits, which are critical constraints on formation theories, it is crucial to understand how orbits evolve over the age of the Solar System. Depending on their timescale, we classify orbital interactions as either short-term (orbital resonances) or long-term (secular evolution). My work involves examples of both interaction types. Resonant history of the small Neptunian satellites In satellite systems, tidal migration brings satellite orbits in and out of resonances. During a resonance passage, satellite orbits change dramatically in a very short period of time. We investigate the resonant history of the six small Neptunian moons. In this unique system, the exotic orbit of the large captured Triton (with a circular, retrograde, and highly tilted orbit) influences the resonances among the small satellites very strongly. We derive an analytical framework which can be applied to Neptune's satellites and to similar systems. Our numerical simulations explain the current orbital tilts of the small satellites as well as constrain key physical parameters of both Neptune and its moons. Secular orbital interactions during eccentricity damping Long-term periodic changes of orbital shape and orientation occur when two or more planets orbit the same star. The variations of orbital elements are superpositions of the same number of fundamental modes as the number of planets in the system. We investigate how this effect interacts with other perturbations imposed by external disturbances, such as the tides and relativistic effects. Through analytical studies of a system consisting of two planets, we find that an external perturbation exerted on one planet affects the other indirectly. We formulate a general theory for how both orbits evolve in response to an arbitrary externally-imposed slow change in eccentricity.

CCD photometry of distant comets. III - Ensemble properties of Jupiter-family comets

We describe the results of a ground-based observational "snapshot" study of Jupiter-family comets in the heliocentric range 2.29 AU less than or equal to R-h less than or equal to 5.72 AU. Results are presented based on observations from the 1m JKT on the island of La Palma. A total of 25 comets were targeted with 15 being positively detected. Broad-band VRI photometry was performed to determine dimensions, colour indices, and dust production rates in terms of the "A frho" formalism. The results for selected comets are compared with previous investigations. Ensemble properties of the Jupiter- family population have been investigated by combining the results presented here with those of Lowry et al. (1999), and Lowry & Fitzsimmons (2001). We find that the cumulative size distribution of the Jupiter-family comets can be described by a power law of the form; Sigma(> r) proportional to r(-1.6+/- 0.1). This size distribution is considerably shallower than that found for the observed Edgeworth-Kuiper belt objects, which may reflect either an intrinsic difference at small km- sizes in the belt, or the various processes affecting the nuclei of comets as their orbits evolve from the Edgeworth- Kuiper belt to the inner Solar system. Also, there would appear to be no correlation between nuclear absolute magnitude and perihelion distance. Finally, for the sample of active comets, there is a distinct correlation between absolute R band magnitude and perihelion distance, which can be explained by either a discovery bias towards brighter comets or in terms of "rubble" mantle formation.

Exploration of the Kuiper Belt by High-Precision Photometric Stellar Occultations: First Results

We report here the first detection of hectometer-size objects by the method of serendipitous stellar occultation. This method consists of recording the diffraction shadow created when an object crosses the observer's line of sight and occults the disk of a background star. One of our detections is most consistent with an object between Saturn and Uranus. The two other diffraction patterns detected are caused by Kuiper Belt objects beyond 100 AU from the Sun and hence are the farthest known objects in the solar system. These detections show that the Kuiper Belt is much more extended than previously believed and that the outer part of the disk could be composed of smaller objects than the inner part. This gives critical clues to understanding the problem of the formation of the outer planets of the solar system.

The nuclei of comets 7P/Pons-Winnecke, 14P/Wolf and 92P/Sanguin

Jupiter Family comets (JFCs) are short period comets which have recently entered the inner solar system, having previously orbited in the Kuiper Belt since the formation of the planets. We used two nights on the 3.6 m New Technology Telescope (NTT) at the European Southern Observatory, to obtain VRI photometry of three JFCs; 7P/Pons-Winnecke, 14P/Wolf and 92P/Sanguin. These were observed to be stellar in appearance. We find mean effective radii of 2.24 ± 0.02 km for 7P, 3.16 ± 0.01 km for 14P and 2.08 ± 0.01 km for 92P, assuming a geometric albedo of 0.04. From light-curves for each comet we find rotation periods of 7.53 ± 0.10 and 6.22 ± 0.05 h for 14P and 92P respectively. 7P exhibits brightness variations which imply a rotation period of 6.8 = Prot = 9.5 h. Assuming the nuclei to be ellipsoidal the measured brightness variations imply minimum axial ratios a/b of 1.3 ± 0.1 for 7P and 1.7 ± 0.1 for both 14P and 92P. This in turn implies minimum densities of 0.23 ± 0.08 g cm-3 for 7P, 0.32 ± 0.02 g cm-3 for 14P and 0.49 ± 0.06 g cm-3 for 92P. Finally, we measure colour indices of (V-R) = 0.40 ± 0.05 and (R-I) = 0.41 ± 0.06 for 7P/Pons-Winnecke, (V-R) = 0.57 ± 0.07 and (R-I) = 0.51 ± 0.06 for 14P/Wolf, and (V-R) = 0.54 ± 0.04 and (R-I) = 0.54 ± 0.04 for 92P/Sanguin.

Direct Detection of the Asteroidal YORP Effect

The Yarkovsky-O'Keefe-Radzievskii-Paddack effect is believed to alter the spin states of small bodies in the solar system. However, evidence for the effect has so far been indirect. Here we report precise optical photometric and radar observations acquired over four years of a small near-Earth asteroid (54509) 2000 PH5. We find the asteroid has been continuously increasing its rotation rate omega over this period by domega/dt = (2.0+/- 0.2) x 10-4 deg/day2. We simulated the close Earth approaches from 2001-2005, showing that gravitational torques cannot explain the observed spin-rate increase. Dynamical simulations also suggest that 2000 PH5 may reach a rotation period of ~20 seconds towards the end of its expected lifetime

Serendipitous Asteroid Lightcurve Survey Using SuperWASP

The SuperWASP project is an ultra-wide angle search for extra solar planetary transits. However, it can also serendipitously detect solar system objects, such as asteroids and comets. Each SuperWASP instrument consists of up to eight cameras, combined with high-quality peltier-cooled CCDs, which photometrically survey large numbers of stars in the magnitude range 7 15. Each camera covers a 7.8 × 7.8 degree field of view. Located on La Palma, the SuperWASP-I instrument has been observing the Northern Hemisphere with five cameras since its inauguration in April 2004. The ultra-wide angle field of view gives SuperWASP the possibility of discovering new fast moving (near to Earth) asteroids that could have been missed by other instruments. However, it provides an excellent opportunity to produce a magnitude-limited lightcurve survey of known main belt asteroids. As slow moving asteroids stay within a single SuperWASP field for several weeks, and may be seen in many fields, a survey of all objects brighter than magnitude 15 is possible. This will provide a significant increase in the total number of lightcurves available for statistical studies without the inherent bias against longer periods present in the current data sets. We present the methodology used in the automated collection of asteroid data from SuperWASP and some of the first examples of lightcurves from numbered asteroids.

An orbital period of 0.94days for the hot-Jupiter planet WASP-18b

The `hot Jupiters' that abound in lists of known extrasolar planets are thought to have formed far from their host stars, but migrate inwards through interactions with the proto-planetary disk from which they were born, or by an alternative mechanism such as planet-planet scattering. The hot Jupiters closest to their parent stars, at orbital distances of only ~0.02 astronomical units, have strong tidal interactions, and systems such as OGLE-TR-56 have been suggested as tests of tidal dissipation theory. Here we report the discovery of planet WASP-18b with an orbital period of 0.94days and a mass of ten Jupiter masses (10MJup), resulting in a tidal interaction an order of magnitude stronger than that of planet OGLE-TR-56b. Under the assumption that the tidal-dissipation parameter Q of the host star is of the order of 106, as measured for Solar System bodies and binary stars and as often applied to extrasolar planets, WASP-18b will be spiralling inwards on a timescale less than a thousandth that of the lifetime of its host star. Therefore either WASP-18 is in a rare, exceptionally short-lived state, or the tidal dissipation in this system (and possibly other hot-Jupiter systems) must be much weaker than in the Solar System.

Turning Asteroids into Comets with ELTs

It is clear that ELTs will be able to detect extremely weak outgassing from Solar system bodies via a number of different methods. Occultations will allow probing for outgassing around 20 km main-belt asteroids. Imaging can reveal dust emission rates of only milligrams/second in the inner solar system, while sublimation rates of gasses should be measurable down to gram/second levels. Suitable targets will be identified via the coming all-sky surveys, through both the classical dynamical Tisserand Invariant and long-baseline lightcurves. It is possible that using these methods, ELTs may allow the discovery of much more activity throughout the Solar system than is presently known.

A high C/O ratio and weak thermal inversion in the atmosphere of exoplanet WASP-12b

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior, as opposed to the silicate-dominated composition found on Earth; the atmosphere can also differ from those in the Solar System. The solar C/O is 0.54 (ref. 3). Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b (ref. 6) that reveals C/O>=1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapour and enhanced in methane, each by more than two orders of magnitude compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T>2,500K) of WASP-12b lacks a prominent thermal inversion (or stratosphere) and has very efficient day-night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

Physical and dynamical characterisation of low ΔV NEA (190491) 2000 FJ10

Aims: We investigated the physical properties and dynamical evolution of near-Earth asteroid (NEA) (190491) 2000 FJ10 in order to assess the suitability of this accessible NEA as a space mission target. Methods: Photometry and colour determination were carried out with the 1.54 m Kuiper Telescope (Mt Bigelow, USA) and the 10 m Southern African Large Telescope (SALT; Sutherland, South Africa) during the object's recent favourable apparition in 2011-12. During the earlier 2008 apparition, a spectrum of the object in the 6000-9000 Angstrom region was obtained with the 4.2 m William Herschel Telescope (WHT; Canary Islands, Spain). Interpretation of the observational results was aided by numerical simulations of 1000 dynamical clones of 2000 FJ10 up to 106 yr in the past and in the future. Results: The asteroid's spectrum and colours determined by our observations suggest a taxonomic classification within the S-complex although other classifications (V, D, E, M, P) cannot be ruled out. On this evidence, it is unlikely to be a primitive, relatively unaltered remnant from the early history of the solar system and thus a low priority target for robotic sample return. Our photometry placed a lower bound of 2 h to the asteroid's rotation period. Its absolute magnitude was estimated to be 21.54 ± 0.1 which, for a typical S-complex albedo, translates into a diameter of 130 ± 20 m. Our dynamical simulations show that it has likely been an Amor for the past 105 yr. Although currently not Earth-crossing, it will likely become so during the period 50-100 kyr in the future. It may have arrived from the inner or central main belt >1 Myr ago as a former member of a low-inclination S-class asteroid family. Its relatively slow rotation and large size make it a suitable destination for a human mission. We show that ballistic Earth-190491-Earth transfer trajectories with ΔV <2 km s-1 at the asteroid exist between 2052 and 2061. Based on observations made with the Southern African Large Telescope (SALT).

Main-belt Comet P/2012 T1 (PANSTARRS)

We present initial results from observations and numerical analyses aimed at characterizing the main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between 2012 October and 2013 February using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research Telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 μm that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of Q CN <1.5 × 1023 mol s-1, from which we infer a water production rate of Q_H_2O100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveals that it is dynamically linked to the ~155 Myr old Lixiaohua asteroid family. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration, and made possible by the generous financial support of the W. M. Keck Foundation, the Magellan Telescopes located at Las Campanas Observatory, Chile, and the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, e Inovação (MCTI) da República Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU).