Variance spectral analysis of superconducting gravimeter (SG) decadal data (noise inclusive) is presented suggesting that the Earth tectonogenesis is based on magnification of the mass (mainly the mantle) mechanical resonance, in addition to or instead of previously hypothesized causes. Here the use of raw (gapped and unaltered) data is regarded as the criterion for a physical result validity, so data were not altered in any way. Then analogously to the atmospheric tidal forcing of global high frequency free oscillation, I propose that the Moon synodically recurring pull could likewise drive the long-periodic (12 to 120 minutes) oscillation of the Earth. To demonstrate this, I show that the daily magnitudes of mass (gravity) oscillation, as a relative measure of the non-stationary Earth kinetic energy, get synodically periodic while correlating up to 0.97 with seismic energies on the day of shallow and 3 days before deep earthquakes. The forced oscillator equations for the mantle usual viscosity and the Earth springtide and grave mode periods successfully model an identical 3 days phase. Finally, whereas reports on gravest earthquakes (of around M9.5) put the maximum coseismic displacement at around 10 m, the same equations predict the maximum displacement as 9.8 m, too. Hence, the same mechanism that causes bridges to collapse under the soldiers step marching could be making the lithosphere fail under the springtide induced magnification of mantle resonance resulting in strong earthquakes of unspecified type. If this assertion is correct, then many if not most large earthquakes could be spatially and temporally predictable.

We propose a new turbulence closure model based on the budget equations for the key second moments: turbulent kinetic and potential energies: TKE and TPE (comprising the turbulent total energy: TTE = TKE + TPE) and vertical turbulent fluxes of momentum and buoyancy (proportional to potential temperature). Besides the concept of TTE, we take into account the non-gradient correction to the traditional buoyancy flux formulation. The proposed model grants the existence of turbulence at any gradient Richardson number, Ri. Instead of its critical value separating – as usually assumed – the turbulent and the laminar regimes, it reveals a transition interval, 0.1< Ri <1, which separates two regimes of essentially different nature but both turbulent: strong turbulence at Ri<<1; and weak turbulence, capable of transporting momentum but much less efficient in transporting heat, at Ri>1. Predictions from this model are consistent with available data from atmospheric and lab experiments, direct numerical simulation (DNS) and large-eddy simulation (LES).

In this paper we very preliminarily investigate the possibility of measuring the post-Newtonian general relativistic gravitoelectric and gravitomagnetic components of the acceleration of gravity on the Earth, in continuous regime, with two absolute measurements at the equator and the south pole with superconducting gravimeters. The magnitudes of such relativistic effects are 10^-10 m s^-2 and 10^-11 m s^-2, respectively. Unfortunately, the present-day uncertainties in the Earth’s geodetic parameters which enter the classical Newtonian terms induce systematic errors 1-2 orders of magnitude larger than the relativistic ones. Moreover, a \sim 1 ngal sensitivity can be reached by the currently available superconducting gravimeters, but only for relative measurements.

We have developed a characterization of the geological evolution of the Earths atmosphere and surface in order to model the observable spectra of an Earth-like planet through its geological history. These calculations are designed to guide the interpretation of an observed spectrum of such a planet by future instruments that will characterize exoplanets. Our models focus on spectral features that either imply habitability or are required for habitability. These features are generated by H2O, CO2, CH4, O2, O3, N2O, and vegetation-like surface albedos. We chose six geological epochs to characterize. These epochs exhibit a wide range in abundance for these molecules, ranging from a CO2 rich early atmosphere, to a CO2/CH4-rich atmosphere around 2 billion years ago to a present-day atmosphere. We analyzed the spectra to quantify the strength of each important spectral feature in both the visible and thermal infrared spectral regions, and the resolutions required to unambiguously observe the features for each epoch. We find a wide range of spectral resolutions required for observing the different features. For example, H2O and O3 can be observed with relatively low resolution, while O2 and N2O require higher resolution. We also find that the inclusion of clouds in our models significantly affects both the strengths and resolutions required to observe all spectral features.

The requirements on space missions designed to study Terrestrial exoplanets are discussed. We then investigate whether the design of such a mission, specifically the Darwin nulling interferometer, can be carried out in a simplified scenario. The key element here is accepting somewhat higher levels of stellar leakage. We establish detailed requirements resulting from the scientific rationale for the mission, and calculate detailed parameters for the stellar suppression required to achieve those requirements. We do this utilizing the Darwin input catalogue. The dominating noise source for most targets in this sample is essentially constant for all targets, while the leakage diminishes with the square of the distance. This means that the stellar leakage has an effect on the integration time only for the nearby stars, while for the more distant targets its influence decreases significantly. We assess the impact of different array configurations and nulling profiles and identify the stars for which the detection efficiency can be maximized.

The Darwin and TPF-I missions are Infrared free flying interferometer missions based on nulling interferometry. Their main objective is to detect and characterize other Earth-like planets, analyze the composition of their atmospheres and their capability to sustain life, as we know it. Darwin and TPF-I are currently in study phase. A number of mission architectures of 3 and 4 free flying telescopes are evaluated on the basis of the interferometer’s response, ability to distinguish multiple planet signatures and starlight rejection capabilities. The characteristics of the new configurations are compared also to the former, more complex Bowtie baseline architectures as well as evaluated on base of their science capability.

The search for life on extrasolar planets is based on the assumption that one can screen extrasolar planets for habitability spectroscopically. The first space born instruments able to detect as well as characterize extrasolar planets, Darwin and terrestrial planet finder (TPF-I and TPF-C) are scheduled to launch before the end of the next decade. The composition of the planetary surface, atmosphere, and its temperature-pressure profile influence a detectable spectroscopic signal considerably. For future space-based missions it will be crucial to know this influence to interpret the observed signals and detect signatures of life in remotely observed atmospheres. We give an overview of biomarkers in the visible and IR range, corresponding to the TPF-C and TPF-I/DARWIN concepts, respectively. We also give an overview of the evolution of biomarkers over time and its implication for the search for life on extrasolar Earth-like planets. We show that atmospheric features on Earth can provide clues of biological activities for at least 2 billion years.

The DARWIN mission is an Infrared free flying interferometer mission based on the new technique of nulling interferometry. Its main objective is to detect and characterize other Earth-like planets, analyze the composition of their atmospheres and their capability to sustain life, as we know it. DARWIN is currently in definition phase. This PhD work that has been undertaken within the DARWIN team at the European Space Agency (ESA) addresses two crucial aspects of the mission. Firstly, a DARWIN target star list has been established that includes characteristics of the target star sample that will be critical for final mission design, such as, luminosity, distance, spectral classification, stellar variability, multiplicity, location and radius of the star. Constrains were applied as set by planet evolution theory and mission architecture. Secondly, a number of alternative mission architectures have been evaluated on the basis of interferometer response as a function of wavelength, achievable modulation efficiency, number of telescopes and starlight rejection capabilities. The study has shown that the core mission goals should be achievable with a lower level of complexity as compared to the current baseline configuration.

Impact basins identified by Mariner 10 and Messenger flyby images provide us a fossilized record of the impactor flux of asteroids on Mercury during the last stages of the early Solar System. The distribution of these basins is not uniform across the surface, and is consistent with a primordial synchronous rotation (Wieczorek et al. 2012). By analyzing the size of the impacts, we show that the distribution for asteroid diameters D < 110 km is compatible with an index power law of 1.2, a value that matches the predicted primordial distribution of the main-belt. We then derive a simple collisional model coherent with the observations, and when combining it with the secular evolution of the spin of Mercury, we are able to reproduce the present 3/2 spin-orbit resonance (about 50% of chances), as well as a primordial synchronous rotation. This result is robust with respect to variations in the dissipation and collisional models, or in the initial spin state of the planet.

We investigate baroclinic instability in flow conditions relevant to hot extrasolar planets. The instability is important for transporting and mixing heat, as well as for influencing large-scale variability on the planets. Both linear normal mode analysis and non-linear initial value calculations are carried out — focusing on the freely-evolving, adiabatic situation. Using a high-resolution general circulation model (GCM) which solves the traditional primitive equations, we show that large-scale jets similar to those observed in current GCM simulations of hot extrasolar giant planets are likely to be baroclinically unstable on a timescale of few to few tens of planetary rotations, generating cyclones and anticyclones that drive weather systems. The growth rate and scale of the most unstable mode obtained in the linear analysis are in qualitative, good agreement with the full non-linear calculations. In general, unstable jets evolve differently depending on their signs (eastward or westward), due to the change in sign of the jet curvature. For jets located at or near the equator, instability is strong at the flanks — but not at the core. Crucially, the instability is either poorly or not at all captured in simulations with low resolution and/or high artificial viscosity. Hence, the instability has not been observed or emphasized in past circulation studies of hot extrasolar planets.

In this work, we investigate the eclipse timing of the polar binary HU Aquarii that has been observed for almost two decades. Recently, Qian et al. attributed large (O-C) deviations between the eclipse ephemeris and observations to a compact system of two massive jovian companions. We improve the Keplerian, kinematic model of the Light Travel Time (LTT) effect and re-analyse the whole currently available data set. We add almost 60 new, yet unpublished, mostly precision light curves obtained using the time high-resolution photo-polarimeter OPTIMA, as well as photometric observations performed at the MONET/N, PIRATE and TCS telescopes. We determine new mid–egress times with a mean uncertainty at the level of 1 second or better. We claim that because the observations that currently exist in the literature are non-homogeneous with respect to spectral windows (ultraviolet, X-ray, visual, polarimetric mode) and the reported mid–egress measurements errors, they may introduce systematics that affect orbital fits. Indeed, we find that the published data, when taken literally, cannot be explained by any unique solution. Many qualitatively different and best-fit 2-planet configurations, including self-consistent, Newtonian N-body solutions may be able to explain the data. However, using high resolution, precision OPTIMA light curves, we find that the (O-C) deviations are best explained by the presence of a single circumbinary companion orbiting at a distance of ~4.5 AU with a small eccentricity and having ~7 Jupiter-masses. This object could be the next circumbinary planet detected from the ground, similar to the announced companions around close binaries HW Vir, NN Ser, UZ For, DP Leo or SZ Her, and planets of this type around Kepler-16, Kepler-34 and Kepler-35.

The dominant non-instrumental background source for space-based infrared observatories is the zo- diacal light. We present Spitzer Infrared Array Camera (IRAC) measurements of the zodiacal light at 3.6, 4.5, 5.8, and 8.0 {\mu}m, taken as part of the instrument calibrations. We measure the changing surface brightness levels in approximately weekly IRAC observations near the north ecliptic pole (NEP) over the period of roughly 8.5 years. This long time baseline is crucial for measuring the annual sinusoidal variation in the signal levels due to the tilt of the dust disk with respect to the ecliptic, which is the true signal of the zodiacal light. This is compared to both Cosmic Background Explorer Diffuse Infrared Background Experiment (COBE DIRBE) data and a zodiacal light model based thereon. Our data show a few percent discrepancy from the Kelsall et al. (1998) model including a potential warping of the interplanetary dust disk and a previously detected overdensity in the dust cloud directly behind the Earth in its orbit. Accurate knowledge of the zodiacal light is important for both extragalactic and Galactic astronomy including measurements of the cosmic infrared background, absolute measures of extended sources, and comparison to extrasolar interplanetary dust models. IRAC data can be used to further inform and test future zodiacal light models.

Numerical simulations of global three-dimensional (3D), self-gravitating discs with a gap opened by an embedded planet are presented. The simulations are customised to examine planetary gap stability. Previous results, obtained by Lin & Papaoizou from two-dimensional (2D) disc models, are reproduced in 3D. These include (i) the development of vortices associated with local vortensity minima at gap edges and their merging on dynamical timescales in weakly self-gravitating discs, (ii) the increased number of vortices as the strength of self-gravity is increased and their resisted merging, and (iii) suppression of the vortex instability and development of global spiral arms associated with local vortensity maxima in massive discs. The vertical structure of these disturbances are examined. In terms of the relative density perturbation, the vortex disturbance has weak vertical dependence when self-gravity is neglected. Vortices become more vertically stratified with increasing self-gravity. This effect is seen even when the unperturbed region around the planet’s orbital radius has a Toomre stability parameter ~10. The spiral modes display significant vertical structure at the gap edge, with the midplane density enhancement being several times larger than that near the upper disc boundary. However, for both instabilities the vertical Mach number is typically a few per cent,and on average vertical motions near the gap edge do not dominate horizontal motions.

As planets form they tidally interact with their natal disks. Though the tidal perturbation induced by Earth and super-Earth mass planets is generally too weak to significantly modify the structure of the disk, the interaction is potentially strong enough to cause the planets to undergo rapid type I migration. This physical process may provide a source of short-period super-Earths, though it may also pose a challenge to the emergence and retention of cores on long-period orbits with sufficient mass to evolve into gas giants. Previous numerical simulations have shown that the type I migration rate sensitively depends upon the circumstellar disk’s properties, particularly the temperature and surface density gradients. Here, we derive these structure parameters for 1) a self-consistent viscous-disk model based on a constant \alpha-prescription, 2) an irradiated disk model that takes into account heating due to the absorption of stellar photons, and 3) a layered-accretion disk model with variable \alpha-parameter. We show that in the inner viscously-heated regions of typical protostellar disks, the horseshoe and corotation torques of super-Earths can exceed their differential Lindblad torque and cause them to undergo outward migration. However, the temperature profile due to passive stellar irradiation causes type I migration to be inwards throughout much of the disk. For disks in which there is outwards migration, we show that location and the mass range of the “planet traps” depends on some uncertain assumptions adopted for these disk models. Competing physical effects may lead to dispersion in super-Earths’ mass-period distribution.

Aims. In order to search for and study the nature of the low-amplitude and long-periodic radial velocity (RV) variations of massive stars, we have been carrying out a precise RV survey for supergiants that lie near or inside the Cepheid instability strip. Methods. We have obtained high-resolution spectra of {\alpha} Per (F5 Ib) from November 2005 to September 2011 using the fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at Bohyunsan Optical Astronomy Observatory (BOAO). Results. Our measurements reveal that {\alpha} Per shows a periodic RV variation of 128 days and a semi-amplitude of 70 m/s. We find no strong correlation between RV variations and bisector velocity span (BVS), but the 128-d peak is indeed present in the BVS variations among several other significant peaks in periodogram. Conclusions. {\alpha} Per may have an exoplanet, but the combined data spanning over 20 years seem to suggest that the 128-d RV variations have not been stable on long-term scale, which is somewhat difficult to reconcile with the exoplanet explanation. We do not exclude the pulsational nature of the 128-d variations in {\alpha} Per. Although we do not find clear evidence for surface activity or rotational modulations by spots, coupled with the fact that the expected rotation period is ~ 130 days the rotational modulation seems to be the most likely cause of the RV variations. More observational data and research are needed to clearly determine the origin of RV the variations in {\alpha} Per.

Astrophysical media are turbulent and therefore reconnection should be treated in the presence of pre-existing turbulence. We consider the model of fast magnetic reconnection in Lazarian & Vishniac (1999) which predicts that the rate of reconnection is controlled by the intensity and the injection scale of turbulent motions. We provide new evidence of successful testing of the model and argue that the model presents a generic set up for astrophysical reconnection events. We study particle acceleration that takes place in volumes of driven turbulence as well turbulent volumes in the presence of large scale reconnection. We show that in the latter case the acceleration is of the first order Fermi type thus supporting the model of acceleration proposed in Gouveia dal Pino & Lazarian (2005).

Context. The discovery of about 700 extrasolar planets, so far, has lead to the first statistics concerning extrasolar planets. The presence of giant planets seems to depend on stellar metallicity and mass. For example, they are more frequent around metal-rich stars,with an exponential increase in planet occurrence rates with metallicity. Aims. We analyzed two samples of metal-poor stars (-2.0 \leq [Fe/H] \leq 0.0) to see if giant planets are indeed rare around these objects. Radial velocity datasets were obtained with two different spectrographs (HARPS and HIRES). Detection limits for these data,expressed in minimum planetary mass and period, are calculated. These produce trustworthy numbers for the planet frequency. Methods. A general Lomb Scargle (GLS) periodogram analysis was used together with a bootstrapping method to produce the detection limits. Planet frequencies were calculated based on a binomial distribution function within metallicity bins. Results. Almost all hot Jupiters and most giant planets should have been found in these data. Hot Jupiters around metal-poor stars have a frequency lower than 1.0% at one sigma. Giant planets with periods up to 1800 days, however, have a higher frequency of $f_p = 2.63^{+2.5}_{-0.8}%$. Taking into account the different metallicities of the stars, we show that giant planets appear to be very frequent $(f_p = 4.48^{+4.04}_{-1.38}%)$ around stars with [Fe/H] > -0.7, while they are rare around stars with [Fe/H] \leq -0.7 (\leq 2.36% at one sigma). Conclusions. Giant planet frequency is indeed a strong function of metallicity, even in the low-metallicity tail. However, the frequencies are most likely higher than previously thought.

We present 8 years of high-precision radial velocity (RV) data for HD 204313 from the 2.7 m Harlan J. Smith Telescope at McDonald Observatory. The star is known to have a giant planet (M sin i = 3.5 M_J) on a ~1900-day orbit, and a Neptune-mass planet at 0.2 AU. Using our own data in combination with the published CORALIE RVs of Segransan et al. (2010), we discover an outer Jovian (M sin i = 1.6 M_J) planet with P ~ 2800 days. Our orbital fit suggests the planets are in a 3:2 mean motion resonance, which would potentially affect their stability. We perform a detailed stability analysis, and verify the planets must be in resonance.

We present the preliminary results of an analysis of the sub-populations within the near-Earth asteroids, including the Atens, Apollos, Amors, and those that are considered potentially hazardous using data from the Wide-field Infrared Survey Explorer (WISE). In order to extrapolate the sample of objects detected by WISE to the greater population, we determined the survey biases for asteroids detected by the project’s automated moving object processing system (known as NEOWISE) as a function of diameter, visible albedo, and orbital elements. Using this technique, we are able to place constraints on the number of potentially hazardous asteroids (PHAs) larger than 100 m and find that there are $\sim4700\pm1450$ such objects. As expected, the Atens, Apollos, and Amors are revealed by WISE to have somewhat different albedo distributions, with the Atens being brighter than the Amors. The cumulative size distributions of the various near-Earth object (NEO) subgroups vary slightly between 100 m and 1 km. A comparison of the observed orbital elements of the various sub-populations of the NEOs with the current best model is shown.

Since 1995, more than 500 extrasolar planets have been discovered orbiting very close to their parent star, where they experience strong tidal interactions. Their orbital evolution depends on the physical mechanisms that cause tidal dissipation, and these are still not well understood. We refine the theory of the equilibrium tide in fluid bodies that are partly or entirely convective, to predict the dynamical evolution of the systems. In particular, we examine the validity of modeling the tidal dissipation by the quality factor Q, as is commonly done. We consider here the simplest case where the considered star or planet rotates uniformly, all spins are aligned, and the companion is reduced to a point-mass. The first manifestation of the tide is to distort the shape of the star or planet adiabatically along the line of centers. This generates the divergence-free velocity field of the adiabatic equilibrium tide which is decoupled from the dynamical tide. The tidal kinetic energy is dissipated into heat through turbulent friction, which is modeled here as an eddy-viscosity acting on the adiabatic tidal flow. This dissipation induces a second velocity field, the dissipative equilibrium tide, which is in quadrature with the exciting potential; it is responsible for the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which one derives characteristic evolution times. The rate at which the system evolves depends on the physical properties of tidal dissipation, and specifically on how the eddy viscosity varies with tidal frequency and on the thickness of the convective envelope for the fluid equilibrium tide. At low frequency, this tide retards by a constant time delay, whereas it lags by a constant angle when the tidal frequency exceeds the convective turnover rate.

Single metal polluted white dwarfs with no dusty disks are believed to be actively accreting metals from a circumstellar disk of gas caused by the destruction of asteroids perturbed by planetary systems. We report, for the first time, the detection of circumstellar Ca~II gas in absorption around the DAZ WD~1124-293, which lacks an infrared excess. We constrain the gas to $>$7 $R_{\rm WD}$ and $<$32000~AU, and estimate it to be at $\sim$54~R$_{\rm WD}$, well within WD~1124-293's tidal disruption radius. This detection is based on several epochs of spectroscopy around the Ca~II H and K lines ($\lambda$=3968\AA, 3933\AA) with the MIKE spectrograph on the Magellan/Clay Telescope at Las Campanas Observatory. We confirm the circumstellar nature of the gas by observing nearby sightlines and finding no evidence for gas from the local interstellar medium. Through archival data we have measured the equivalent width of the two photospheric Ca lines over a period of 11 years. We see $$ R$_{\rm \oplus}$ using the WASP survey. The presence of gas in orbit around WD~1124-293 implies that most DAZs could harbor planetary systems. Since 25-30\% of white dwarfs show metal line absorption, the dynamical process for perturbing small bodies must be robust.

The restricted (equilateral) four-body problem consists of three bodies of masses m1, m2 and m3 (called primaries) lying in a Lagrangian configuration of the three-body problem i.e., they remain fixed at the apices of an equilateral triangle in a rotating coordinate system. A massless fourth body moves under the Newtonian gravitational law due to the three primaries, as in the Restricted three-body problem (R3BP), the fourth mass does not affect the motion of the three primaries. In this paper we explore symmetric periodic orbits of the restricted four-body problem (R4BP) for the case of two equal masses where they satisfy approximately the Routh’s critical value.

This treatise outlines how a real non-systematic based Pioneer anomaly, with its implied violation (re: ‘low’ mass bodies only) of both general relativity’s weak equivalence principle and the Newtonian inverse-square law, can be successfully modelled. These theoretical hurdles and various awkward observational constraints, such as the low value of Pioneer 11’s anomaly pre-Saturn encounter, have not been convincingly modelled to date. Notwithstanding the recent trend to embrace a non-constant Sun/Earth-directed heat based explanation of this anomalous deceleration, the actual: nature, direction, and temporal and spatial variation of the Pioneer anomaly remain an open arena of research. Working backwards from the observational evidence, and rethinking: time, mass, quantum entanglement and non-locality, we hypothesise a mechanism involving a quantum mechanical energy source and a new type of ‘gravitational’ field; neither of which lie within general relativity’s domain of formulation/application. By way of a systemic conservation of energy principle, an internally inexpressible (aggregate) non-inertial energy discrepancy/uncertainty – involving a myriad of quantum (lunar/third-body residing) atomic and molecular systems moving in analog curved spacetime – is (non-locally) re-expressed externally as a (rotating) non-Euclidean spatial geometry perturbation. At a moving body each “rotating space-warp” induces sinusoidal proper acceleration and speed perturbations, as well as a path-based constant (per cycle) rate of speed shortfall relative to predictions that omit the additional effect. ‘Solutions’ of the new model extend to: the Earth flyby anomaly; solar system related large-scale anomalies in the CMB radiation data; the nature of dark energy; and how a theory of everything unification agenda is inadvertently impeding a deeper understanding of physical reality and quantum entanglement.

We study the secondary electron emissions induced by the impact of electrons on dust grains and the resulting dust charging processes in the nonequilibrium dusty plasma with power-law distributions. We derive new expressions of the secondary emitted electron flux and the dust charging currents that are generalized by the power-law q-distributions, where the nonlinear core functions are numerically studied for the nonextensive parameter q. Our numerical analyses show that the power-law q-distribution of the primary electrons has a significant effect on the secondary emitted electron flux as well as the dust charging currents, and this effect depends strongly on the ratio of the electrostatic potential energy of the primary electrons at the dust grain’s surface to the thermodynamic energy, implying that a competition in the dusty plasma between these two energies plays a crucial role in this novel effect.

We study the secondary electron emissions induced by the impact of electrons on dust grains and the resulting dust charging processes in the nonequilibrium dusty plasma with power-law distributions. We derive new expressions of the secondary emitted electron flux and the dust charging currents that are generalized by the power-law q-distributions, where the nonlinear core functions are numerically studied for the nonextensive parameter q. Our numerical analyses show that the power-law q-distribution of the primary electrons has a significant effect on the secondary emitted electron flux as well as the dust charging currents, and this effect depends strongly on the ratio of the electrostatic potential energy of the primary electrons at the dust grain’s surface to the thermodynamic energy, implying that a competition in the dusty plasma between these two energies plays a crucial role in this novel effect.

We present an investigation of spin-orbit angles for planetary system candidates reported by Kepler. By combining the rotational period $P_s$ inferred from the flux variation due to starspots and the projected rotational velocity $V\sin I_s$ and stellar radius obtained by a high resolution spectroscopy, we attempt to estimate the inclination $I_s$ of the stellar spin axis with respect to the line-of-sight. For transiting planetary systems, in which planetary orbits are edge-on seen from us, the stellar inclination $I_s$ can be a useful indicator of a spin-orbit alignment/misalignment. We newly conducted spectroscopic observations with Subaru/HDS for 15 KOI systems, whose lightcurves show periodic flux variations. After detailed analyses of their lightcurves and spectra, it turned out that some of them are binaries, or the flux variations are too coherent to be caused by starspots, probably representing ellipsoidal variations, and consequently we could constrain stellar inclinations $I_s$ for eight systems. Among them, KOI-262 and 280 are in good agreement with $I_s=90^\circ$ suggesting a spin-orbit alignment, while at least one system, KOI-261, shows a possible spin-orbit misalignment. We also obtain a small $I_s$ for KOI-1463, but the transiting companion seems to be a star rather than a planet. The results for KOI-257, 269, 367, and 974 are ambiguous, and can be explained with either misalignments or moderate differential rotation. Since our method can be applied to any system having starspots regardless of the planet size, future observations will allow for the expansion of the parameter space in which the spin-orbit relations are investigated.

The observed lifetimes of gaseous protoplanetary discs place strong constraints on gas and ice giant formation in the core accretion scenario. The approximately 10-Earth-mass solid core responsible for the attraction of the gaseous envelope has to form before gas dissipation in the protoplanetary disc is completed within 1-10 million years. Building up the core by collisions between km-sized planetesimals fails to meet this time-scale constraint, especially at wide stellar separations. Nonetheless, gas-giant planets are detected by direct imaging at wide orbital distances. In this paper, we numerically study the growth of cores by the accretion of cm-sized pebbles loosely coupled to the gas. We measure the accretion rate onto seed masses ranging from a large planetesimal to a fully grown 10-Earth-mass core and test different particle sizes. The numerical results are in good agreement with our analytic expressions, indicating the existence of two accretion regimes, one set by the azimuthal and radial particle drift for the lower seed masses and the other, for higher masses, by the velocity at the edge of the Hill sphere. In the former, the optimally accreted particle size increases with core mass, while in the latter the optimal size is centimeters, independent of core mass. We discuss the implications for rapid core growth of gas-giant and ice-giant cores. We conclude that pebble accretion can resolve the long-standing core accretion time-scale conflict. This requires a near-unity dust-to-gas ratio in the midplane, particle growth to mm and cm and the formation of massive planetesimals or low radial pressure support. The core growth time-scale is shortened by a factor 30-1,000 at 5 AU and by a factor 100-10,000 at 50 AU, compared to the gravitationally focused accretion of, respectively, low-scale-height planetesimal fragments or standard km-sized planetesimals.

Gravitational instability is one of considerable mechanisms to explain the formation of giant planets. We study the gravitational stability for the protoplanetary disks around a protostar. The temperature and Toomre’s Q-value are calculated by assuming local equilibrium between viscous heating and radiative cooling (local thermal equilibrium). We assume constant $\alpha$ viscosity and use a cooling function with realistic opacity. Then, we derive the critical surface density $\Sigma_{\rm{c}}$ that is necessary for a disk to become gravitationally unstable as a function of $r$. This critical surface density $\Sigma_{\rm c}$ is strongly affected by the temperature dependence of the opacity. At the radius $r_{\rm c}\sim 20$AU, where ices form, the value of $\Sigma_{\rm c}$ changes discontinuously by one order of magnitude. This $\Sigma_{\rm c}$ is determined only by local thermal process and criterion of gravitational instability. By comparing a given surface density profile to $\Sigma_{\rm c}$, one can discuss the gravitational instability of protoplanetary disks. As an example, we discuss the gravitational instability of two semi-analytic models for protoplanetary disks. One is the steady state accretion disk, which is realized after the viscous evolution. The other is the disk that has the same angular momentum distribution with its parent cloud core, which corresponds to the disk that has just formed. As a result, it is found that the disks tend to become gravitationally unstable for $r\ge r_{\rm c}$ because ices enable the disks to become low temperature. In the region closer to the protostar than $r_{\rm c}$, it is difficult for a typical protoplanetary disk to fragment because of the high temperature and the large Coriolis force. From this result, we conclude that the fragmentation near the central star is possible but difficult.

The Karin family is the first recognized very young asteroid family that was created by an asteroid breakup only 5.8 Myr ago. As 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 collisional event in terms of their rotational status and surface color. We have been carrying out a series of photometric observations of the Karin family asteroids, and here we report the analysis result of lightcurves including the rotation period of eleven members as well as those of an interloper asteroid: (832) Karin, (4507) 1990 FV (an interloper), (7719) 1997 GT36, (10783) 1999 RB9, (11728) Einer, (13765) Nansmith, (16706) Svojsik, (28271) 1999 CK16, (40917) 1999 TR171, (43032) 1999 VR26, (69880) 1998 SQ81, and (71031) 1999 XE68. As for four of them we estimated their absolute magnitudes H_R and the slope parameter G_R of the solar phase curves: (832) Karin, (4507) 1990 FV, (13765) Nansmith, and (69880) 1998 SQ81. Relation between the reduced peak-to-peak variation magnitude and the rotation period of the asteroids indicates that there is a general trend that elongated members have lower spin rate, and less elongated members have higher spin rate. The mean rotation rate of the Karin family members turned out to be much lower than those of NEAs (near-Earth asteroids) or smaller MBAs (main belt asteroids; diameter D130 km), suggesting the existence of post-disruption evolution process over the spin state of the older asteroids, such as the YORP effect.

We report the detection of two new planets from the Anglo-Australian Planet Search. These planets orbit two stars each previously known to host one planet. The new planet orbiting HD 142 has a period of 6005\pm427 days, and a minimum mass of 5.3M_Jup. HD142c is thus a new Jupiter analog: a gas-giant planet with a long period and low eccentricity (e = 0.21 \pm 0.07). The second planet in the HD 159868 system has a period of 352.3\pm1.3 days, and m sin i=0.73\pm0.05 M_Jup. In both of these systems, including the additional planets in the fitting process significantly reduced the eccentricity of the original planet. These systems are thus examples of how multiple-planet systems can masquerade as moderately eccentric single-planet systems.

We report the detection of two new planets from the Anglo-Australian Planet Search. These planets orbit two stars each previously known to host one planet. The new planet orbiting HD 142 has a period of 6005\pm427 days, and a minimum mass of 5.3M_Jup. HD142c is thus a new Jupiter analog: a gas-giant planet with a long period and low eccentricity (e = 0.21 \pm 0.07). The second planet in the HD 159868 system has a period of 352.3\pm1.3 days, and m sin i=0.73\pm0.05 M_Jup. In both of these systems, including the additional planets in the fitting process significantly reduced the eccentricity of the original planet. These systems are thus examples of how multiple-planet systems can masquerade as moderately eccentric single-planet systems.

We report the discovery of two new transiting planets from the WASP survey. WASP-42 b is a 0.500 +- 0.035 M_J planet orbiting a K1 star at a separation of 0.0548 +- 0.0017 AU with a period of 4.9816872 +- 0.0000073 days. The radius of WASP-42 is 1.080 +- 0.057 R_J while its equilibrium temperature is T_eq = 995 +- 34 K. We detect some evidence of a small but non-zero eccentricity of e = 0.060 +- 0.013. WASP-49 b is a 0.378 +- 0.027 M_J planet around an old G6 star. It has a period of 2.7817387 +- 5.6 x 10-6 days and a separation of 0.0379 +- 0.0011 AU. This planet is slightly bloated, having a radius of 1.115 +- 0.056 R_J and an equilibrium temperature of T_eq = 1369 +- 42 K. Both planets have been followed up intensively in photometry, in total we have obtained 5 full and one partial transit light curves of WASP-42 and 4 full and one partial light curves of WASP-49 using the Euler-Swiss, TRAPPIST and Faulkes South telescopes.

The MESSENGER spacecraft conducted its first flyby of Mercury on 14th January 2008, followed by two subsequent encounters on 6th October 2008 and 29th September 2009, prior to Mercury orbit insertion on 18th March 2011. We have reviewed MESSENGER flight telemetry and X-ray Spectrometer observations from the first two encounters, and correlate several prominent features in the data with the presence of astrophysical X-ray sources in the instrument field of view. We find that two X-ray peaks attributed in earlier work to the detection of suprathermal electrons from the Mercury magnetosphere, are likely to contain a significant number of events that are of astrophysical origin. The intensities of these two peaks cannot be explained entirely on the basis of astrophysical sources, and we support the previous suprathermal explanation but suggest that the electron fluxes derived in those studies be revised to correct for a significant astrophysical signal.

I summarize recent surveys of protoplanetary disks at millimeter wavelengths and show that the distribution of luminosity, equivalent to the mass in small dust grains, declines rapidly. This contrasts with statistics on the lifetime of disks from infrared observations and the high occurrence of planets from radial velocity and transit surveys. I suggest that these disparate results can be reconciled if most of the dust in a disk is locked up in millimeter and larger sized particles within about 2 Myr. This general statistical result on disk evolution agrees with detailed modeling of a small number of individual disks and with cosmochemical measurements of chondrule ages.

The recent Kepler discovery of KOI-152 reveals a system of three hot super-Earth candidates that are in or near a 4:2:1 mean motion resonance. It is unlikely that they formed in situ, the planets probably underwent orbital migration during the formation and evolution process. The small semimajor axes of the three planets suggest that migration stopped at the inner edge of the primordial gas disk. In this paper we focus on the influence of migration halting mechanisms, including migration “dead zones”, and inner truncation by the stellar magnetic field. We show that the stellar accretion rate, stellar magnetic field and the speed of migration in the proto-planetary disk are the main factors affecting the final configuration of KOI-152. Our simulations suggest that three planets may be around a star with low star accretion rate or with high magnetic field. On the other hand, slow type I migration, which decreases to one tenth of the linear analysis results, favors forming the configuration of KOI-152. Under such formation scenario, the planets in the system are not massive enough to open gaps in the gas disk. The upper limit of the planetary masses are estimated to be about $15,~19$, and $24 M_\oplus$, respectively. Our results are also indicative of the near Laplacian configurations that are quite common in planetary systems.

The Habitable Zone for a given star describes the range of circumstellar distances from the star within which a planet could have liquid water on its surface, which depends upon the stellar properties. Here we describe the development of the Habitable Zone concept, its application to our own Solar System, and its subsequent application to exoplanetary systems. We further apply this to planets in extreme eccentric orbits and show how they may still retain lifebearing properties depending upon the percentage of the total orbit which is spent within the Habitable Zone.

We present the results of a search for planetary companions orbiting near hot Jupiter planet candidates (Jupiter-size candidates with orbital periods near 3 days) identified in the Kepler data through its sixth quarter of science operations. Special emphasis is given to companions between the 2:1 interior and exterior mean-motion resonances. A photometric transit search excludes companions with sizes ranging from roughly 2/3 to 5 times the size of the Earth, depending upon the noise properties of the target star. A search for dynamically induced deviations from a constant period (transit timing variations or TTVs) also shows no significant signals. In contrast, comparison studies of warm Jupiters (with slightly larger orbits) and hot Neptune-size candidates do exhibit signatures of additional companions with these same tests. These differences between hot Jupiters and other planetary systems denote a distinctly different formation or dynamical history.

We determine the fraction of F, G, and K dwarfs in the Solar Neighborhood hosting hot jupiters as measured by the California Planet Survey from the Lick and Keck planet searches. We find the rate to be 1.2\pm0.38%, which is consistent with the rate reported by Mayor et al. (2011) from the HARPS and CORALIE radial velocity surveys. These numbers are more than double the rate reported by Howard et al. (2011) for Kepler stars and the rate of Gould et al. (2006) from the OGLE-III transit search, however due to small number statistics these differences are of only marginal statistical significance. We explore some of the difficulties in estimating this rate from the existing radial velocity data sets and comparing radial velocity rates to rates from other techniques.

The population of Neptune Trojans is believed to be bigger than that of Jupiter Trojans and that of asteroids in the main belt, although only eight members of this far distant asteroid swarm have been observed up to now. Six leading Neptune Trojans around the Lagrange point L4 discovered earlier have been studied in detail, but two trailing ones found recently around the L5 point, 2004 KV18 and 2008 LC18, have not been investigated yet. In this paper, we report our investigations on the dynamical behaviors of these two new Neptune Trojans. Our calculations show that the asteroid 2004 KV18 is a temporary Neptune Trojan. Most probably, it was captured into the trailing Trojan cloud no earlier than 203kyr ago, and it will not keep this identity no later than 165kyr in future. Based on the statistics on our orbital simulations, we argue that this object is more like a scattered Kuiper belt object. On the contrary, the orbit of asteroid 2008 LC18 is much more stable. Among the clone orbits spread within the orbital uncertainties, a considerable portion of clones may survive on the L5 tadpole orbits for 4Gyr. The strong dependence of the stability on the semimajor axis and resonant angle suggests that further observations are badly needed to confine the orbit in the stable region. We also discuss the implications of the existence and dynamics of these two trailing Trojans on the Solar system history.

The population of Neptune Trojans is believed to be bigger than that of Jupiter Trojans and that of asteroids in the main belt, although only eight members of this far distant asteroid swarm have been observed up to now. Six leading Neptune Trojans around the Lagrange point L4 discovered earlier have been studied in detail, but two trailing ones found recently around the L5 point, 2004 KV18 and 2008 LC18, have not been investigated yet. In this paper, we report our investigations on the dynamical behaviors of these two new Neptune Trojans. Our calculations show that the asteroid 2004 KV18 is a temporary Neptune Trojan. Most probably, it was captured into the trailing Trojan cloud no earlier than 203kyr ago, and it will not keep this identity no later than 165kyr in future. Based on the statistics on our orbital simulations, we argue that this object is more like a scattered Kuiper belt object. On the contrary, the orbit of asteroid 2008 LC18 is much more stable. Among the clone orbits spread within the orbital uncertainties, a considerable portion of clones may survive on the L5 tadpole orbits for 4Gyr. The strong dependence of the stability on the semimajor axis and resonant angle suggests that further observations are badly needed to confine the orbit in the stable region. We also discuss the implications of the existence and dynamics of these two trailing Trojans on the Solar system history.

Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous paper, we performed three-dimensional N-body simulations of sequential collisions of aggregates composed of a number of sub-micron-sized icy particles and examined radii of gyration (and bulk densities) of the obtained aggregates. We showed that collisional compression of aggregates is not efficient and that aggregates remain fluffy. In the present study, we examine geometrical cross sections of the aggregates. Their cross sections decreases due to the compression as well as their gyration radii. It is found that a relation between the cross section and the gyration radius proposed by Okuzumi et al. is valid for the compressed aggregates. We also refine the compression model proposed in our previous paper. The refined model enables us to calculate the evolution of both gyration radii and cross sections of growing aggregates and reproduces well our numerical results of sequential aggregate collisions. The refined model can describe non-equal-mass collisions as well as equal-mass case. Although we do not take into account oblique collisions in the present study, oblique collisions would further hinder compression of aggregates.

We use the distribution of extrasolar planets in circular orbits around stars with surface convective zones detected by ground based transit searches to constrain how efficiently tides raised by the planet are dissipated on the parent star. We parameterize this efficiency as a tidal quality factor (Q*). We conclude that the population of currently known planets is inconsistent with Q*<10^7 at the 99% level. Previous studies show that values of Q* between 10^5 and 10^7 are required in order to explain the orbital circularization of main sequence low mass binary stars in clusters, suggesting that different dissipation mechanisms might be acting in the two cases, most likely due to the very different tidal forcing frequencies relative to the stellar rotation frequency occurring for star–star versus planet–star systems.

We report on the detection of infrared light from the super-Earth 55 Cnc e, based on four occultations obtained with Warm Spitzer at 4.5 microns. Our data analysis consists of a two-part process. In a first step, we perform individual analyses of each dataset and compare several baseline models to optimally account for the systematics affecting each lightcurve. We apply independent photometric correction techniques, including polynomial detrending and pixel-mapping, that yield consistent results at the 1-sigma level. In a second step, we perform a global MCMC analysis including all four datasets, that yields an occultation depth of 131+-28ppm, translating to a brightness temperature of 2360+-300 K in the IRAC-4.5 micron channel. This occultation depth suggests a low Bond albedo coupled to an inefficient heat transport from the planetary dayside to the nightside, or else possibly that the 4.5-micron observations probe atmospheric layers that are hotter than the maximum equilibrium temperature (i.e., a thermal inversion layer or a deep hot layer). The measured occultation phase and duration are consistent with a circular orbit and improves the 3-sigma upper limit on 55 Cnc e’s orbital eccentricity from 0.25 to 0.06.

Signatures of the dissipation region of collisionless magnetic reconnection are investigated by the Geotail spacecraft for the 15 May 2003 event. The energy dissipation in the rest frame of the electron’s bulk flow is considered in an approximate form D*_e, which is validated by a particle-in-cell simulation. The dissipation measure is directly evaluated from the {plasma moments}, the electric field, and the magnetic field. Using D*_e, a compact dissipation region is successfully detected in the vicinity of the possible X-point in Geotail data. The dissipation rate is 45 pWm**{-3}. The length of the dissipation region is estimated to 1–2 local ion inertial length. The Lorentz work W, the work rate by Lorentz force to plasmas, is also introduced. It is positive over the reconnection region and it has a peak around the pileup region away from the X-point. These new measures D*_e and W provide useful information to understand the reconnection structure.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m$. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

High rate sampling detectors measuring the potential difference between the main body and boom antennas of interplanetary spacecraft have been shown to be efficient means to measure the voltage pulses induced by nano dust impacts on the spacecraft body itself (see Meyer-Vernet et al, Solar Phys. 256, 463 (2009)). However, rough estimates of the free charge liberated in post impact expanding plasma cloud indicate that the cloud’s own internal electrostatic field is too weak to account for measured pulses as the ones from the TDS instrument on the STEREO spacecraft frequently exceeding 0.1 V/m$. In this paper we argue that the detected pulses are not a direct measure of the potential structure of the plasma cloud, but are rather the consequence of a transitional interruption of the photoelectron return current towards the portion of the antenna located within the expanding cloud.

We present a method to optimize absorption cells for precise wavelength calibration in the near-infrared. We apply it to design and optimize methane isotopologue cells for precision radial velocity measurements in the K band. We also describe the construction and installation of two such cells for the CSHELL spectrograph at NASA’s IRTF. We have obtained their high-resolution laboratory spectra, which we can then use in precision radial velocity measurements and which can also have other applications. In terms of obtainable RV precision methane should out-perform other proposed cells, such as the ammonia cell ($^{14}$NH$_{3}$) recently demonstrated on CRIRES/VLT. The laboratory spectra of Ammonia and the Methane cells show strong absorption features in the H band that could also be exploited for precision Doppler measurements. We present spectra and preliminary radial velocity measurements obtained during our first-light run. These initial results show that a precision down to 20-30 m s$^{-1}$ can be obtained using a wavelength interval of only 5 nm in the K band and S/N$\sim$150. This supports the prediction that a precision down to a few m s$^{-1}$ can be achieved on late M dwarfs using the new generation of NIR spectrographs, thus enabling the detection of terrestrial planets in their habitable zones. Doppler measurements in the NIR can also be used to mitigate the radial velocity jitter due to stellar activity enabling more efficient surveys on young active stars.

An analytic expression for the variation in surface and sub-surface temperature is developed for worlds whose surface pressures are nearly constant with latitude and longitude and whose atmospheres are in vapor-pressure equilibrium with the dominant surface volatiles. Such worlds include the current Pluto and Triton, and other volatile-covered Kuiper Belt Objects during some portion of their heliocentric orbit. The expressions also apply on worlds with negligible horizontal heat flow, such as asteroids. Temperature variations in volatile-covered or bare areas as a function of time is derived in terms of three thermal parameters relating to (1) the thermal wave within the substrate, (2) the energy needed to heat an isothermal volatile slab, and (3) the buffering by the latent heat needed to change the atmospheric surface pressure. For Pluto’s current surface pressure (~17 microbar), atmospheric buffering dominates over subsurface effects on diurnal timescales, and should keep the surface pressure over a Pluto day constant to within 0.2%.

We investigate the orbital stability of a putative Jovian planet reported by Ramm et. al in a compact binary \nu Octantis. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework Mechanic which we developed to run massive numerical experiments on CPU clusters. The code is illustrated on a model Hamiltonian introduced by Froeschl\’e et al. We confirm that the \nu Octantis planet could reside in a retrograde orbit, according with a hypothesis of Eberle & Cuntz. It may be present in a zone of stable motions which has a structure of the Arnold web formed due to overlapping of low-order mean motion resonances and their sub-resonances. We also re-analyzed the available radial velocity data in terms of self-consistent Newtonian N-body model. We found stable best-fit solutions that obey the observational constraints. They correspond to retrograde, strictly anti-aligned orbits of the binary and the planet. However, these solutions are confined in very small stable regions of the phase space. The presence of a real planet in the system is still questionable, because its formation would be hindered by strong dynamical perturbations.

We investigate the orbital stability of a putative Jovian planet reported by Ramm et. al in a compact binary \nu Octantis. Our numerical study makes use of a new computational Message Passing Interface (MPI) framework Mechanic which we developed to run massive numerical experiments on CPU clusters. The code is illustrated on a model Hamiltonian introduced by Froeschl\’e et al. We confirm that the \nu Octantis planet could reside in a retrograde orbit, according with a hypothesis of Eberle & Cuntz. It may be present in a zone of stable motions which has a structure of the Arnold web formed due to overlapping of low-order mean motion resonances and their sub-resonances. We also re-analyzed the available radial velocity data in terms of self-consistent Newtonian N-body model. We found stable best-fit solutions that obey the observational constraints. They correspond to retrograde, strictly anti-aligned orbits of the binary and the planet. However, these solutions are confined in very small stable regions of the phase space. The presence of a real planet in the system is still questionable, because its formation would be hindered by strong dynamical perturbations.

Recent infrared observations have demonstrated that photoevaporation driven by high-energy photons from the central star contributes to the dispersal of protoplanetary disks. Here, we show that photoevaporative winds should produce a detectable free-free continuum emission given the range of stellar ionizing photons and X-ray luminosities inferred for young sun-like stars. We point out that VLA observations of the nearby disk around TWHya might have already detected this emission at centimeter wavelengths and calculate the wind electron density and mass flow rate. We also estimate the intensities of H radio recombination lines tracing the wind and discuss which ones could be detected with current instrumentation. The detection and profiles of these recombination lines would unambiguously prove our inference of free-free emission from photoevaporating disks like TWHya. In addition, radio/millimeter data can help constraining wind parameters such as temperature and electron density that are fundamental in measuring mass flow rates.

We examine a low energy mechanism for the transfer of meteoroids between two planetary systems embedded in a star cluster using quasi-parabolic orbits of minimal energy. Using Monte Carlo simulations, we find that the exchange of meteoroids could have been significantly more efficient than previously estimated. Our study is relevant to astrobiology as it addresses whether life on Earth could have been transferred to other planetary systems in the solar system’s birth cluster and whether life on Earth could have been transferred here from beyond the solar system. In the solar system, the timescale over which solid material was delivered to the region from where it could be transferred via this mechanism likely extended to several hundred million years (as indicated by the 3.8-4.0 Ga epoch of the Late Heavy Bombardment). This timescale could have overlapped with the lifetime of the Solar birth cluster (~100-500 Myr). Therefore, we conclude that lithopanspermia is an open possibility if life had an early start. Adopting parameters from the minimum mass solar nebula, considering a range of planetesimal size distributions derived from observations of asteroids and Kuiper Belt Objects and theoretical coagulation models, and taking into account Oort Cloud formation models, the expected number of bodies with mass > 10 kg that could have been transferred between the Sun and its nearest cluster neighbor could be of the order of 1E14-3E16, with transfer timescales of 10s Myr. We estimate that of the order of 3E8 x l(km) could potentially be life-bearing, where l(km) is the depth of the Earth crust in km that was ejected as the result of the early bombardment.

Identifying liquid water on the surface of planets is a high priority, as this traditionally defines habitability. One proposed signature of oceans is specular reflection (“glint”), which increases the apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth-like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations even though we do not include specular reflection in our model. This false positive is the product of two generic properties: 1) for modest obliquities, a planet’s poles receive less orbit-averaged stellar flux than its equator, so the poles are more likely to be covered in highly reflective snow and ice, and 2) we show that reflected light from a modest-obliquity planet at crescent phases probes higher latitudes than at gibbous phases, therefore a planet’s apparent albedo will naturally increase at crescent phase. We suggest that this “latitude-albedo effect” will operate even for large obliquities: in that case the equator receives less orbit-averaged flux than the poles, and the equator is preferentially sampled at crescent phase. Using rotational and orbital color variations to map the surfaces of directly imaged planets and estimate their obliquity will therefore be a necessary pre-condition for properly interpreting their reflected phase variations. The latitude-albedo effect is a particularly convincing glint false positive for zero-obliquity planets, and such worlds are not amenable to latitudinal mapping. This effect severely limits the utility of specular reflection for detecting oceans on exoplanets.

Debris disks with extremely large infrared excesses (fractional luminosities $> 10^{-2}$) are rare. Those with ages between 30 and 130 Myr are of interest because their evolution has progressed well beyond that of protoplanetary disks (which dissipate with a timescale of order 3 Myr), yet they represent a period when dynamical models suggest that terrestrial planet building may still be progressing through large, violent collisions that could yield large amounts of debris and large infrared excesses. For example, our Moon was formed through a violent collision of two large proto-planets during this age range. We report two disks around the solar-like stars ID8 and HD 23514 in this age range where the 24 {\mu}m infrared excesses vary on timescales of a few years, even though the stars are not variable in the optical. Variations this rapid are difficult to understand if the debris is produced by collisional cascades, as it is for most debris disks. It is possible that the debris in these two systems arises in part from condensates from silicate-rich vapor produced in a series of violent collisions among relatively large bodies. If their evolution is rapid, the rate of detection of extreme excesses would indicate that major collisions may be relatively common in this age range.

The G0V dwarf TrES-2A, which is transited by a hot Jupiter, is one of the main short-cadence targets of the Kepler telescope and, therefore, among the photometrically best-studied planetary systems known today. Given the near-grazing geometry of the planetary orbit, TrES-2 offers an outstanding opportunity to search for changes in its orbital geometry. Our study focuses on the secular change in orbital inclination reported in previous studies. We present a joint analysis of the first four quarters of Kepler photometry together with the publicly available ground-based data obtained since the discovery of TrES-2b in 2006. We use a common approach based on the latest information regarding the visual companion of TrES-2A and stellar limb darkening to further refine the orbital parameters. We find that the Kepler observations rule out a secular inclination change of previously claimed order as well as variations of the transit timing, however, they also show slight indication for further variability in the inclination which remains marginally significant.

Our aim is to identify and classify mean-motion resonances (MMRs) for the coplanar circular restricted three-body problem (CR3BP) for mass ratios between 0.10 and 0.50. Our methods include the maximum Lyapunov exponent, which is used as an indicator for the location of the resonances, the Fast Fourier Transform (FFT) used for determining what kind of resonances are present, and the inspection of the orbital elements to classify the periodicity. We show that the 2:1 resonance occurs the most frequently. Among other resonances, the 3:1 resonance is the second most common, and furthermore both 3:2 and 5:3 resonances occur more often than the 4:1 resonance. Moreover, the resonances in the coplanar CR3BP are classified based on the behaviour of the orbits. We show that orbital stability is ensured for high values of resonance (i.e., high ratios) where only a single resonance is present. The resonances attained are consistent with the previously established resonances for the Solar System, i.e., specifically, in regards to the asteroid belt. Previous work employed digital filtering and Lyapunov characteristic exponents to determine stochasticity of the eccentricity, which is found to be consistent with our usage of Lyapunov exponents as an alternate approach based on varying the mass ratio instead of the eccentricity. Our results are expected to be of principal interest to future studies, including augmentations to observed or proposed resonances, of extra-solar planets in binary stellar systems.

The Kepler-11 star hosts at least six transiting super-Earth planets detected through the precision photometric observations of the KEPLER mission (Lissauer et al.). In this paper, we re-analyze the available KEPLER data, using the direct N-body approach rather than an indirect TTV method in the discovery paper. The orbital modeling in the realm of the direct approach relies on the whole data set, rather than the times of mid-transits only. Most of the results in the original paper are confirmed and extended. We constrained the mass of the outermost planet g to less than 30 Earth masses. The mutual inclinations between orbits b and c as well as between orbits d and e are determined with a good precision, in the range of [1,5] degrees. Having several solutions to four qualitative orbital models of the Kepler-11 system, we analyze its global dynamics with the help of dynamical maps. They reveal very complex structure of the phase space with narrow regions of regular motion. The dynamics are governed by a dense net of three- and four-body mean motion resonances, forming the Arnold web. Overlapping of these resonances is a main source of instability. We found that the Kepler-11 system may be long-term stable only in particular multiple resonant configurations with small relative inclinations. The mass-radius data derived for all companions reveal a clear anti-correlation between the mean density of the planets with their distance from the star. It may reflect the formation and early evolution history of the system.

The potential existence of a giant planet orbiting within a few AU of a stellar remnant has profound implications for both the survival and possible regeneration of planets during post-main sequence stellar evolution. This paper reports Hubble Space Telescope Fine Guidance Sensor and U.S. Naval Observatory relative astrometry of GD 66, a white dwarf thought to harbor a giant planet between 2 and 3 AU based on stellar pulsation arrival times. Combined with existing infrared data, the precision measurements here rule out all stellar-mass and brown dwarf companions, implying that only a planet remains plausible, if orbital motion is indeed the cause of the variations in pulsation timing.

Comet 17P/Holmes underwent a dramatic outburst in October 2007, caused by the sudden fragmentation of its nucleus and the production of a large quantity of grains scattering sunlight. We report on 90 GHz continuum observations carried out with the IRAM Plateau de Bure interferometer on 27.1 and 28.2 October 2007 UT, i.e., 4-5 days after the outburst. These observations probed the thermal radiation of large dust particles, and therefore provide the best constraints on the mass in the ejecta debris. The thermal emission of the debris was modelled and coupled to a time-dependent description of their expansion after the outburst. The analysis was performed in the Fourier plane. Visibilities were computed for the two observing dates and compared to the data to measure their velocity and mass. Optical data and 250-GHz continuum measurements published in the literature were used to further constrain the dust kinematics and size distribution. Two distinct dust components in terms of kinematic properties are identified in the data. The large-velocity component, with typical velocities V0 of 50-100 m/s for 1 mm particles, displays a steep size distribution with a size index estimated to q = -3.7 (\pm0.1), assuming a minimum grain size of 0.1 \mum. It corresponds to the fast expanding shell observed in optical images. The slowly-moving “core” component (V0 = 7-9 m/s) detected near the nucleus has a size index |q| < 3.4 and contains a higher proportion of large particles than the shell. The dust mass in the core is in the range 0.1-1 that of the shell. Using optical constants pertaining to porous grains (50% porosity) made of astronomical silicates mixed with water ice (48% in mass), the total dust mass Mdust injected by the outburst is estimated to 4-14 x 10**11 kg, corresponding to 3-9% the nucleus mass.

As the number of discovered extrasolar planets has been increasing, diversity of planetary systems requires studies of new formation scenarios. It is important to study satellite formation in circumplanetary disks, which is often viewed as analogous to formation of rocky planets in protoplanetary disks. We investigated satellite formation from satellitesimals around giant planets through N-body simulations that include gravitational interactions with a circumplanetary gas disk. Our main aim is to reproduce the observable properties of the Galilean satellites around Jupiter through numerical simulations, as previous N-body simulations have not explained the origin of the resonant configuration. We performed accretion simulations based on the work of Sasaki et al. (2010), in which an inner cavity is added to the model of Canup & Ward (2002, 2006). We found that several satellites are formed and captured in mutual mean motion resonances outside the disk inner edge and are stable after rapid disk gas dissipation, which explains the characteristics of the Galilean satellites. In addition, owing to the existence of the disk edge, a radial compositional gradient of the Galilean satellites can also be reproduced. An additional objective of this study is to discuss orbital properties of formed satellites for a wide range of conditions by considering large uncertainties in model parameters. Through numerical experiments and semianalytical arguments, we determined that if the inner edge of a disk is introduced, a Galilean-like configuration in which several satellites are captured into a 2:1 resonance outside the disk inner cavity is almost universal. In fact, such a configuration is produced even for a massive disk and rapid type I migration. This result implies the inevitability of a Galilean satellite formation in addition to providing theoretical predictions for extrasolar satellites.

Phase reddening is an effect that produces an increase of the spectral slope and variations in the strength of the absorption bands as the phase angle increases. In order to understand its effect on spectroscopic observations of asteroids, we have analyzed the visible and near-infrared spectra (0.45-2.5 \mu m) of 12 near-Earth asteroids observed at different phase angles. All these asteroids are classified as either S-complex or Q-type asteroids. In addition, we have acquired laboratory spectra of three different types of ordinary chondrites at phase angles ranging from 13\degree to 120\degree. We have found that both asteroid and meteorite spectra show an increase in band depths with increasing phase angle. The spectral slope of the ordinary chondrites spectra shows a significant increase with increasing phase angle for g > 30\degree. Variations in band centers and band area ratio (BAR) values were also found, however they seems to have no significant impact on the mineralogical analysis. Our study showed that the increase in spectral slope caused by phase reddening is comparable to certain degree of space weathering. In particular, an increase in phase angle in the range of 30\degree to 120\degree will produce a reddening of the reflectance spectra equivalent to exposure times of ~ 0.1×10^6 to 1.3×10^6 years at about 1 AU from the Sun. Furthermore, we found that under some circumstances phase reddening could lead to an ambiguous taxonomic classification of asteroids.

We present HST ultraviolet spectroscopy of the white dwarfs PG0843+516, PG1015+161, SDSS1228+1040, and GALEX1931+0117, which accrete circumstellar planetary debris formed from the destruction of asteroids. Combined with optical data, a minimum of five and a maximum of eleven different metals are detected in their photospheres. With metal sinking time scales of only a few days, these stars are in accretion/diffusion equilibrium, and the photospheric abundances closely reflect those of the circumstellar material. We find C/Si ratios that are consistent with that of the bulk Earth, corroborating the rocky nature of the debris. Their C/O values are also very similar to those of bulk Earth, implying that the planetary debris is dominated by Mg and Fe silicates. The abundances found for the debris at the four white dwarfs show substantial diversity, comparable at least to that seen across different meteorite classes in the solar system. PG0843+516 exhibits significant over-abundances of Fe and Ni, as well as of S and Cr, which suggests the accretion of material that has undergone melting, and possibly differentiation. PG1015+161 stands out by having the lowest Si abundance relative to all other detected elements. The Al/Ca ratio of the planetary debris around different white dwarfs is remarkably similar. This is analogous to the nearly constant abundance ratio of these two refractory lithophile elements found among most bodies in the solar system. Based on the detection of all major elements of the circumstellar debris, we calculate accretion rates of ~1.7e8g/s ~1.5e9g/s. We detect additional circumstellar absorption in the SiIV 1394,1403 doublet in PG0843+516 and SDSS1228+1040, reminiscent to similar high-ionisation lines seen in white dwarfs in cataclysmic variables. We suspect that these lines originate in hot gas close to the white dwarf, well within the sublimation radius.

We use models of thermal evolution and XUV-driven mass loss to explore the composition and history of low-mass low-density transiting planets. We investigate the Kepler-11 system in detail and provide estimates of both the current and past planetary compositions. We find that a H/He atmosphere on Kepler-11b is highly vulnerable to mass loss. By comparing to formation models, we show that in situ formation of the system is unlikely. Instead we propose that it is a water-rich system of sub-Neptunes that migrated from beyond the snow line. For the broader population of observed planets, we show that there is a threshold in bulk planet density and incident flux above which no low-mass transiting planets have been observed. We suggest that this threshold is due to the instability of H/He atmospheres to XUV-driven mass loss. Importantly, we find that this flux-density threshold is well reproduced by our thermal evolution/contraction models that incorporate a standard mass loss prescription. Treating the planets’ contraction history is essential because the planets have significantly larger radii during the early era of high XUV fluxes. Over time low mass planets with H/He envelopes can be transformed into water-dominated worlds with steam atmospheres or rocky super-Earths. Finally, we use this threshold to provide likely minimum masses and radial velocity amplitudes for the general population of Kepler candidates. Likewise, we use this threshold to provide constraints on the maximum radii of low-mass planets found by radial velocity surveys.

A re-evaluation of time-averaged accretion rates at DBZ-type white dwarfs points to historical, time-averaged rates significantly higher than the currently observed episodes at their DAZ counterparts. The difference between the ongoing, instantaneous accretion rates witnessed at DAZ white dwarfs, which often exceed 1e8 g/s, and those inferred over the past 1e5-1e6 yr for the DBZ stars can be a few orders of magnitude, and therefore must result from high-rate episodes of tens to hundreds of years so they remain undetected to date. This paper explores the likelihood that such brief, intense accretion episodes of gas-phase material can account for existing data. For reasonable assumptions about the circumstellar gas, accretion rates approaching or exceeding 1e15 g/s are possible, similar to rates observed in quiescent cataclysmic variables, and potentially detectable with future x-ray missions or wide-field monitoring facilities. Gaseous debris that is prone to such rapid accretion may be abundant immediately following a tidal disruption event via collisions and sublimation, or if additional bodies impinge upon an extant disk. Particulate disk matter accretes at or near the Poynting-Robertson drag rate for long periods between gas-producing events, consistent with rates inferred for dusty DAZ white dwarfs. In this picture, warm DAZ stars without infrared excesses have rates consistent with accretion from particulate disks that remain undetected. This overall picture has implications for quasi-steady state models of accretion and the derived chemical composition of asteroidal debris in DBZ white dwarfs.

We investigate analytically and numerically a nonlinear modification of the magnetospheric plasma density under the action of the ponderomotive force induced by ULF traveling waves, using the nonlinear stationary force balance equation. This equation is applied to both the dipole and dayside magnetosphere having one and two minima of the geomagnetic field near the magnetospheric boundary. The separate and joint actions of the ponderomotive, centrifugal, and gravitational forces on the density distribution are shown.

Trojan asteroids are minor planets that share the orbit of a planet about the Sun and librate around the L4 or L5 Lagrangian points of stability. Although only three Mars Trojans have been discovered, models suggest that at least ten times this number should exist with diameters >= 1 km. We derive a model that constrains optimal sky search areas and present a strategy for the most efficient use of telescope survey time that maximizes the probability of detecting Mars Trojans. We show that the Gaia space mission could detect any Mars Trojans larger than 1 km in diameter, provided the relative motion perpendicular to Gaia’s CCD array is less than 0.40 arcsec per second.

Debris disks should not be completely gas-free, since there is second generation gas from outgassing of planetesimals and dust grains via sublimation, photodesorption, or collisions, generating a system of dust-to-gas ratio close to unity, where hydrodynamics cannot be ignored. A clumping instability exists in this configuration, that has been hitherto explored only in one-dimensional, incompressible models. We performed 2D numerical compressible models of a disk with comparable amounts of gas and dust to study the growth and development of this instability. Our model solves the momentum equation for the gas and dust, together with energy and continuity equations. We uncover that the backreaction of the drag force from the gas onto the dust shepherds rings, similar to those observed in debris disks and usually attributed to the presence of hypothetical undetected planets. We also uncover that the eccentricity of these rings, usually presented as convincing evidence for the presence of a planet, can actually be simply explained by a standing wave propagating along the ring. The rings support a spectrum of oscillations, with one particular mode representing epicyclic motion. The apparent eccentricity matches the eccentricity in observed systems. This suggests that the planet possibility, though thrilling, is not necessarily required to explain these systems.

We present the first results from our Red Optical Planet Survey (ROPS) to search for low mass planets orbiting late type dwarfs (M5.5V – M9V) in their habitable zones (HZ). Our observations, with the red arm of the MIKE spectrograph (0.5 – 0.9 microns) at the 6.5 m Magellan Clay telescope at Las Campanas Observatory indicate that >= 92 per cent of the flux lies beyond 0.7 microns. We use a novel approach that is essentially a hybrid of the simultaneous iodine and ThAr methods for determining precision radial velocities. We apply least squares deconvolution to obtain a single high S/N ratio stellar line for each spectrum and cross correlate against the simultaneously observed telluric line profile, which we derive in the same way. Utilising the 0.62 – 0.90 micron region, we have achieved an r.m.s. precision of 10 m/s for an M5.5V spectral type star with spectral S/N ~160 on 5 minute timescales. By M8V spectral type, a precision of ~30 m/s at S/N = 25 is suggested, although more observations are needed. An assessment of our errors and scatter in the radial velocity points hints at the presence of stellar radial velocity variations. Of our sample of 7 stars, 2 show radial velocity signals at 6-sigma and 10-sigma of the cross correlation uncertainties. If the signals are planetary in origin, our findings are consistent with estimates of Neptune mass planets that predict a frequency of 13 – 27 per cent for early M dwarfs.Our current analysis indicates the we can achieve a sensitivity that is equivalent to the amplitude induced by a 6 M_Earth planet orbiting in the habitable zone. Based on simulations, we estimate that <10 M_Earth habitable zone planets will be detected in a new stellar mass regime, with <=20 epochs of observations.

We present the first results from our Red Optical Planet Survey (ROPS) to search for low mass planets orbiting late type dwarfs (M5.5V – M9V) in their habitable zones (HZ). Our observations, with the red arm of the MIKE spectrograph (0.5 – 0.9 microns) at the 6.5 m Magellan Clay telescope at Las Campanas Observatory indicate that >= 92 per cent of the flux lies beyond 0.7 microns. We use a novel approach that is essentially a hybrid of the simultaneous iodine and ThAr methods for determining precision radial velocities. We apply least squares deconvolution to obtain a single high S/N ratio stellar line for each spectrum and cross correlate against the simultaneously observed telluric line profile, which we derive in the same way. Utilising the 0.62 – 0.90 micron region, we have achieved an r.m.s. precision of 10 m/s for an M5.5V spectral type star with spectral S/N ~160 on 5 minute timescales. By M8V spectral type, a precision of ~30 m/s at S/N = 25 is suggested, although more observations are needed. An assessment of our errors and scatter in the radial velocity points hints at the presence of stellar radial velocity variations. Of our sample of 7 stars, 2 show radial velocity signals at 6-sigma and 10-sigma of the cross correlation uncertainties. If the signals are planetary in origin, our findings are consistent with estimates of Neptune mass planets that predict a frequency of 13 – 27 per cent for early M dwarfs.Our current analysis indicates the we can achieve a sensitivity that is equivalent to the amplitude induced by a 6 M_Earth planet orbiting in the habitable zone. Based on simulations, we estimate that <10 M_Earth habitable zone planets will be detected in a new stellar mass regime, with <=20 epochs of observations.

John Goodricke and Edward Pigott, working in York, England, between 1781 and 1786, determined the periods of variation of eclipsing binaries such as Algol and Beta Lyrae and speculated that the eclipses of Algol might be caused by a “dark body,” perhaps even a planet. They also determined the periods of variation of the first two known Cepheid variables, the stars whose period-luminosity relation today enables astronomers to determine distances to distant galaxies. Goodricke holds special interest because he was completely deaf and because he died at the age of 21. The lives and work of these two astronomers are described.

As part of NASA’s mission to explore habitable planets orbiting nearby stars, this paper explores the detection and characterization capabilities of a 4-m space telescope plus 50-m starshade located at the Earth-Sun L2 point, a.k.a. the New Worlds Observer (NWO). Our calculations include the true spectral types and distribution of stars on the sky, an iterative target selection protocol designed to maximize efficiency based on prior detections, and realistic mission constraints. We carry out both analytical calculations and simulated observing runs for a wide range in exozodiacal background levels ({\epsilon} = 1 – 100 times the local zodi brightness) and overall prevalence of Earth-like terrestrial planets ({\eta}\oplus = 0.1 – 1). We find that even without any return visits, the NWO baseline architecture (IWA = 65 mas, limiting FPB = 4\times10-11) can achieve a 95% probability of detecting and spectrally characterizing at least one habitable Earth-like planet, and an expectation value of ~3 planets found, within the mission lifetime and {\Delta}V budgets, even in the worst-case scenario ({\eta}\oplus = 0.1 and {\epsilon} = 100 zodis for every target). This achievement requires about one year of integration time spread over the 5 year mission, leaving the remainder of the telescope time for UV-NIR General Astrophysics. Cost and technical feasibility considerations point to a “sweet spot” in starshade design near a 50-m starshade effective diameter, with 12 or 16 petals, at a distance of 70,000-100,000 km from the telescope.

We have identified and quantified semi-major axis drifts in Near-Earth Asteroids (NEAs) by performing orbital fits to optical and radar astrometry of all numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the most reliable and strongest drift rates. Our selection criteria include a Yarkovsky sensitivity metric that quantifies the detectability of semi-major axis drift in any given data set, a signal-to-noise metric, and orbital coverage requirements. In 42 cases, the observed drifts (~10^-3 AU/Myr) agree well with numerical estimates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is the dominant non-gravitational process affecting these orbits, and allows us to derive constraints on asteroid physical properties. In 12 cases, the drifts exceed nominal Yarkovsky predictions, which could be due to inaccuracies in our knowledge of physical properties, faulty astrometry, or modeling errors. If these high rates cannot be ruled out by further observations or improvements in modeling, they would be indicative of the presence of an additional non-gravitational force, such as that resulting from a loss of mass of order a kilogram per second. We define the Yarkovsky efficiency f_Y as the ratio of the change in orbital energy to incident solar radiation energy, and we find that typical Yarkovsky efficiencies are ~10^-5.

GAIA will provide observations of some multiple asteroid and dwarf systems. These observations are a way to determine and improve the quantification of dynamical parameters, such as the masses and the gravity fields, in these multiple systems. Here we investigate this problem in the cases of Pluto’s and Eugenia’s system. We simulate observations reproducing an approximate planning of the GAIA observations for both systems, as well as the New Horizons observations of Pluto. We have developed a numerical model reproducing the specific behavior of multiple asteroid system around the Sun and fit it to the simulated observations using least-square method, giving the uncertainties on the fitted parameters. We found that GAIA will improve significantly the precision of Pluto’s and Charon’s mass, as well as Petit Prince’s orbital elements and Eugenia’s polar oblateness.

On various occasions, several authors suggested to model bodily tides with superposition of two symmetrical bulges. One bulge is always aimed at the secondary, and thus implements the instantaneous reaction of the primary’s shape and potential to the tide-rising gravitational pull exerted on it by the secondary. This portion of the tide is called “adiabatic tide” (Zahn 1966a,b) or “elastic tide” (Ferraz Mello 2012; Krasinsky 2006). The second bulge is assumed to align orthogonally to the direction to the tide-raising secondary. So this bulge is set to implement the entire nonelastic portion of the primary’s deformation. This, second bulge is called “dissipative tide” (Zahn 1966a,b; Krasinsky 2006) or “creep tide” (Ferraz Mello 2012). We demonstrate that the two-bulge method is not a separate approximation, but ensues directly from the Fourier expansion of a linear tidal theory equipped with an arbitrary rheological model involving a departure from elasticity. While less economical mathematically, the two-bulge approach has a good illustrative power, and may be employed (like, e.g., in Remus et al. 2012a,b) on a par with a more concise method of complex amplitudes.

We present continuum high resolution Submillimeter Array (SMA) observations of the transition disk object RX J1633.9-2442, which is located in the Ophiuchus molecular cloud and has recently been identified as a likely site of ongoing giant planet formation. The observations were taken at 340 GHz (880 micron) with the SMA in its most extended configuration, resulting in an angular resolution of 0.3″ (35 AU at the distance of the target). We find that the disk is highly inclined (i ~50 deg) and has an inner cavity ~25 AU in radius, which is clearly resolved by our observations. We simultaneously model the entire optical to millimeter wavelength spectral energy distribution (SED) and SMA visibilities of RX J1633.9-2442 in order to constrain the structure of its disk. We find that an empty cavity ~25 AU in radius is inconsistent with the excess emission observed at 12, 22, and 24 micron. Instead, the mid-IR excess can be modeled by either a narrow, optically thick ring at ~10 AU or an optically thin region extending from ~7 AU to ~25 AU. The inner disk (r < 5 AU) is mostly depleted of small dust grains as attested by the lack of detectable near-IR excess. We also present deep Keck aperture masking observations in the near-IR, which rule out the presence of a companion up to 500 times fainter than the primary star (in K-band) for projected separations in the 5-20 AU range. We argue that the complex structure of the RX J1633.9-2442 disk is best explained by multiple planets embedded within the disk. We also suggest that the properties and incidence of objects such as RX J1633.9-2442, T Cha, and LkCa 15 (and those of the companions recently identified to these two latter objects) are most consistent with the runaway gas accretion phase of the core accretion model, when giant planets gain their envelopes and suddenly become massive enough to open wide gaps in the disk.

It has been suggested that the transition between magnetorotationally active and dead zones in protoplanetary disks should be prone to the excitation of vortices via Rossby wave instability (RWI). However, the only numerical evidence for this has come from alpha disk models, where the magnetic field evolution is not followed, and the effect of turbulence is parametrized by Laplacian viscosity. We aim to establish the phenomenology of the flow in the transition in 3D resistive-magnetohydrodynamical models. We model the transition by a sharp jump in resistivity, as expected in the inner dead zone boundary, using the Pencil Code to simulate the flow. We find that vortices are readily excited in the dead side of the transition. We measure the mass accretion rate finding similar levels of Reynolds stress at the dead and active zones, at the $\alpha\approx 10^{-2}$ level. The vortex sits in a pressure maximum and does not migrate, surviving until the end of the simulation. A pressure maximum in the active zone also triggers the RWI. The magnetized vortex that results should be disrupted by parasitical magneto-elliptic instabilities, yet it subsists in high resolution. This suggests that either the parasitic modes are still numerically damped, or that the RWI supplies vorticity faster than they can destroy it. We conclude that the resistive transition between the active and dead zones in the inner regions of protoplanetary disks, if sharp enough, can indeed excite vortices via RWI. Our results lend credence to previous works that relied on the alpha-disk approximation, and caution against the use of overly reduced azimuthal coverage on modeling this transition.

We present the 1.9–4.2um spectra of the five bright (L<11.2) young stars associated with silhouette disks with moderate to high inclination angle of 39–80deg in the M42 and M43 regions. The water ice absorption is seen toward d121-1925 and d216-0939, while the spectra of d182-316, d183-405, and d218-354 show no water ice feature around 3.1um within the detection limits. By comparing the water ice features toward nearby stars, we find that the water ice absorption toward d121-1925 and d216-0939 most likely originates from the foreground material and the surrounding disk, respectively. The angle of the disk inclination is found to be mainly responsible for the difference of the optical depth of the water ice among the five young stars. Our results suggest that there is a critical inclination angle between 65deg and 75deg for the circumstellar disk where the water ice absorption becomes strong. The average density at the disk surface of d216-0939 was found to be 6.38×10^(-18) g cm^(-3). The water ice absorption band in the d216-0939 disk is remarkable in that the maximum optical depth of the water ice band is at a longer wavelength than detected before. It indicates that the primary carrier of the feature is purely crystallized water ice at the surface of the d216-0939 disk with characteristic size of ~0.8um, which suggests grain growth. This is the first direct detection of purely crystallized water ice in a silhouette disk.

The Damocloids are a group of unusual asteroids, recently enrolling a new member of 2010 EJ104. The dynamical evolution for the Damocloids may uncover a connection passage from the Main Belt, the Kuiper Belt and the scattered disk beyond. According to our simulations, two regions may be considered as possible origin of the Damocloids: the scattered disk, or a part of Oort cloud which will be perturbed to a transient region locating between 700 AU and 1000 AU. Based on the potential origin, the Damocloids can be classified into two types, with relation to their semi-major axes, and about 65.5% Damocloids is classified into type I which mainly originate from Oort cloud. Whether the Damocloids is inactive nuclei of Halley Family Comets may rely on their origin.

The disk around the Herbig Ae star HD\,169142 was imaged and resolved at 18.8 and 24.5\,$\mu$m using Subaru/COMICS. We interpret the observations using a 2D radiative transfer model and find evidence for the presence of a large gap. The MIR images trace dust that emits at the onset of the strong rise in the spectral energy distribution (SED) at 20\,$\mu$m, therefore are very sensitive to the location and characteristics of the inner wall of the outer disk and its dust. We determine the location of the wall to be 23$^{+3}_{-5}$\,AU from the star. An extra component of hot dust must exist close to the star. We find that a hydrostatic optically thick inner disk does not produce enough flux in the NIR and an optically thin geometrically thick component is our solution to fit the SED. Considering the recent findings of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that such disk structures may be common in group I sources. Classification as group I should be considered a support for classification as a transitional disk, though improved imaging surveys are needed to support this speculation.

Motivated by recent discoveries of low-density super-Earths with short orbital periods, we have investigated in-situ accretion of H-He atmospheres on rocky bodies embedded in dissipating warm disks, by simulating quasi-static evolution of atmospheres that connect to the ambient disk. We have found that the atmospheric evolution has two distinctly different outcomes, depending on the rocky body’s mass: While the atmospheres on massive rocky bodies undergo runaway disk-gas accretion, those on light rocky bodies undergo significant erosion during disk dispersal. In the atmospheric erosion, the heat content of the rocky body that was previously neglected plays an important role. We have also realized that the atmospheric mass is rather sensitive to disk temperature in the mass range of interest in this study. Our theory is applied to recently-detected super-Earths orbiting Kepler-11 to examine the possibility that the planets are rock-dominated ones with relatively thick H-He atmospheres. The application suggests that the in-situ formation of the relatively thick H-He atmospheres inferred by structure modeling is possible only under restricted conditions; namely, relatively slow disk dissipation and/or cool environments. This study demonstrates that low-density super-Earths provide important clues to understanding of planetary accretion and disk evolution.

We consider a coplanar system comprised of a massive central body (a star), a less massive secondary (a planet) on a circular orbit, and a test particle on a bound orbit exterior to that of the secondary. The gravitational pull exerted on the test particle by the secondary acts as a small perturbation, wherefore the trajectory of the particle can be described as an ellipse of a precessing perihelion. While the apsidal motion is defined overwhelmingly by the Newtonian portion of the secondary’s gravity, the post-Newtonian portion, too, brings its tiny input. We explore whether this input may be of any astrophysical relevance in the next few decades. We demonstrate that the overall post-Newtonian input of the secondary’s gravity can be split into two parts. One can be expressed via the orbital angular momentum of the secondary, another via its orbital radius. Despite some moderately large numerical factors showing up in the expressions for these two parts, the resulting post-Newtonian contributions from the secondary’s gravity into the apsidal motion of the test particle turn out to be small enough to be neglected in the near-future measurements.

The chance that a planetary system will interact with another member of its host star’s nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multiplanet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form on wide orbits in situ, then there will be discernible differences in the proper motion distributions of a sample of wide-orbit planets compared with a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but auto-ionization is likely to remain the dominant source of free-floating planets.

We present seven new transiting hot Jupiters from the WASP-South survey. The planets are all typical hot Jupiters orbiting stars from F4 to K0 with magnitudes of V = 10.3 to 12.5. The orbital periods are all in the range 3.9–4.6 d, the planetary masses range from 0.4–2.3 Mjup and the radii from 1.1–1.4 Mjup. In line with known hot Jupiters, the planetary densities range from Jupiter-like to inflated (rho = 0.13–1.07 rho_jup). We use the increasing numbers of known hot Jupiters to investigate the distribution of their orbital periods and the 3–4-d “pile-up”.

Fomalhaut is a young, nearby star that is suspected to harbor an infant planetary system, interspersed with one or more belts of dusty debris. We present far-infrared images obtained with the Herschel Space Observatory with an angular resolution between 5.7 and 36.7 arcsec at wavelengths between 70 and 500 micrometer. The images show the main debris belt in great detail. Even at high spatial resolution, the belt appears smooth. The region in between the belt and the central star is not devoid of material; thermal emission is observed here as well. Also at the location of the star, excess emission is detected. We use a dynamical model together with radiative-transfer tools to derive the parameters of the debris disk. We include detailed models of the interaction of the dust grains with radiation, for both the radiation pressure and the temperature determination. Comparing these models to the spatially resolved temperature information contained in the images allows us to place strong constraints on the presence of grains that will be blown out of the system by radiation pressure. We use this to derive the dynamical parameters of the system. The appearance of the belt points towards a remarkably active system in which dust grains are produced at a very high rate by a collisional cascade in a narrow region filled with dynamically excited planetesimals. Dust particles with sizes below the blow-out size are abundantly present. The equivalent of 2000 one-km-sized comets are destroyed every day, out of a cometary reservoir amounting to 110 Earth masses. From comparison of their scattering and thermal properties, we find evidence that the dust grains are fluffy aggregates, which indicates a cometary origin. The excess emission at the location of the star may be produced by hot dust with a range of temperatures, but may also be due to gaseous free-free emission from a stellar wind.

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles’ aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles’ aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

Analyzing exoplanets detected by radial velocity or transit observations, we determine the multiplicity of exoplanet host stars in order to study the influence of a stellar companion on the properties of planet candidates. Matching the host stars of exoplanet candidates detected by radial velocity or transit observations with online multiplicity catalogs in addition to a literature search, 57 exoplanet host stars are identified having a stellar companion. The resulting multiplicity rate of at least 12 percent for exoplanet host stars is about four times smaller than the multiplicity of solar like stars in general. The mass and the number of planets in stellar multiple systems depend on the separation between their host star and its nearest stellar companion, e.g. the planetary mass decreases with an increasing stellar separation. We present an updated overview of exoplanet candidates in stellar multiple systems, including 15 new systems (compared to the latest summary from 2009).

High-magnification microlensing events provide an important channel to detect planets. Perturbations near the peak of a high-magnification event can be produced either by a planet or a binary companion. It is known that central perturbations induced by both types of companions can be generally distinguished due to the basically different magnification pattern around caustics. In this paper, we present a case of central perturbations for which it is difficult to distinguish the planetary and binary interpretations. The peak of a lensing light curve affected by this perturbation appears to be blunt and flat. For a planetary case, this perturbation occurs when the source trajectory passes the negative perturbation region behind the back end of an arrowhead-shaped central caustic. For a binary case, a similar perturbation occurs for a source trajectory passing through the negative perturbation region between two cusps of an astroid-shaped caustic. We demonstrate the degeneracy for 2 high-magnification events of OGLE-2011-BLG-0526 and OGLE-2011-BLG-0950/MOA-2011-BLG-336. For OGLE-2011-BLG-0526, the $\chi^2$ difference between the planetary and binary model is $\sim$ 3, implying that the degeneracy is very severe. For OGLE-2011-BLG-0950/MOA-2011-BLG-336, the stellar binary model is formally excluded with $\Delta \chi^2 \sim$ 105 and the planetary model is preferred. However, it is difficult to claim a planet discovery because systematic residuals of data from the planetary model are larger than the difference between the planetary and binary models. Considering that 2 events observed during a single season suffer from such a degeneracy, it is expected that central perturbations experiencing this type of degeneracy is common.

High-magnification microlensing events provide an important channel to detect planets. Perturbations near the peak of a high-magnification event can be produced either by a planet or a binary companion. It is known that central perturbations induced by both types of companions can be generally distinguished due to the basically different magnification pattern around caustics. In this paper, we present a case of central perturbations for which it is difficult to distinguish the planetary and binary interpretations. The peak of a lensing light curve affected by this perturbation appears to be blunt and flat. For a planetary case, this perturbation occurs when the source trajectory passes the negative perturbation region behind the back end of an arrowhead-shaped central caustic. For a binary case, a similar perturbation occurs for a source trajectory passing through the negative perturbation region between two cusps of an astroid-shaped caustic. We demonstrate the degeneracy for 2 high-magnification events of OGLE-2011-BLG-0526 and OGLE-2011-BLG-0950/MOA-2011-BLG-336. For OGLE-2011-BLG-0526, the $\chi^2$ difference between the planetary and binary model is $\sim$ 3, implying that the degeneracy is very severe. For OGLE-2011-BLG-0950/MOA-2011-BLG-336, the stellar binary model is formally excluded with $\Delta \chi^2 \sim$ 105 and the planetary model is preferred. However, it is difficult to claim a planet discovery because systematic residuals of data from the planetary model are larger than the difference between the planetary and binary models. Considering that 2 events observed during a single season suffer from such a degeneracy, it is expected that central perturbations experiencing this type of degeneracy is common.

We investigate the stability of prograde versus retrograde planets in circular binary systems using numerical simulations. We show that retrograde planets are stable up to distances closer to the perturber than prograde planets. We develop an analytical model to compute the prograde and retrograde mean motion resonances’ locations and separatrices. We show that instability is due to single resonance forcing, or caused by nearby resonances’ overlap. We validate our results regarding the role of single resonances and resonances’ overlap on orbit stability, by computing surfaces of section of the CR3BP. We conclude that the observed enhanced stability of retrograde planets with respect to prograde planets is due to essential differences between the phase-space topology of retrograde versus prograde resonances (at p/q mean motion ratio, prograde resonance is of order p – q while retrograde resonance is of order p + q).

We present results from a simple general circulation model (GCM) of a GJ 1214b-like super-Earth exoplanet. The dynamical core of our model is a scaled-up version of a shallow atmosphere, terrestrial planet GCM that has previously been used for Mars and therefore employs different boundary conditions and physical processes than downsized gas giant models. We assume the planet is tidally locked and has the observed characteristics of GJ-1214b [Charbonneau et al. 2009] for surface mass, surface radius, orbital period, and surface gravitational acceleration. We assume the atmosphere is composed entirely of water vapor. We assume the planet has a surface (i.e., a density discontinuity at depth), which will provide a mechanical drag and affect the radiative balance at the bottom boundary. We assume a gray atmosphere in the IR. We find that a westerly jet is present aloft at the equator and that the longitude of maximum temperature is shifted eastward of the substellar point. A wavenumber-1 feature is present in the equatorial vertical velocity field, indicative of a standing Kelvin and/or Rossby wave. As such, the circulation does not exhibit a cellular structure as on terrestrial Solar System planets; rather, air parcels move up and down on closed horizontal circuits and always return to the same location. The flow at the midlatitudes and poles has both an easterly component and a component that flows poleward along the evening terminator and equatorward along the morning terminator. Temperature inversions exist in the boundary layer and aloft, which are a result of variations in horizontal heat transport. The flow within the boundary layer is more convergent than in the atmosphere aloft. The surface pressure is higher at the poles than the equator. (abridged)

For more than four decades, scientists have been trying to find an answer to one of the most fundamental questions in paleoclimatology, the `faint young Sun problem’. For the early Earth, models of stellar evolution predict a solar energy input to the climate system which is about 25% lower than today. This would result in a completely frozen world over the first two billion years in the history of our planet, if all other parameters controlling Earth’s climate had been the same. Yet there is ample evidence for the presence of liquid surface water and even life in the Archean (3.8 to 2.5 billion years before present), so some effect (or effects) must have been compensating for the faint young Sun. A wide range of possible solutions have been suggested and explored during the last four decades, with most studies focusing on higher concentrations of atmospheric greenhouse gases like carbon dioxide, methane or ammonia. All of these solutions present considerable difficulties, however, so the faint young Sun problem cannot be regarded as solved. Here I review research on the subject, including the latest suggestions for solutions of the faint young Sun problem and recent geochemical constraints on the composition of Earth’s early atmosphere. Furthermore, I will outline the most promising directions for future research. In particular I would argue that both improved geochemical constraints on the state of the Archean climate system and numerical experiments with state-of-the-art climate models are required to finally assess what kept the oceans on the Archean Earth from freezing over completely.

In the context of planet formation, anticyclonic vortices have recently received lots of attention for the role they can play in planetesimals formation. Radial migration of intermediate size solids toward the central star may prevent their growth to larger solid grains. On the other hand, vortices can trap the dust and accelerate this growth, counteracting fast radial transport. Multiple effects have been shown to affect this scenario, such as vortex migration or decay. The aim of this paper is to study the formation of vortices by the Rossby wave instability and their long term evolution in a full three dimensional protoplanetary disc. We use a robust numerical scheme combined with adaptive mesh refinement in cylindrical coordinates, allowing to affordably compute long term 3D evolutions. We consider a full disc stratified both radially and vertically that is prone to formation of vortices by the Rossby wave instability. We show that the 3D Rossby vortices grow and survive over hundreds of years without migration. The localized overdensity which initiated the instability and vortex formation survives the growth of the Rossby wave instability for very long times. When the vortices are no longer sustained by the Rossby wave instability, their shape changes toward more elliptical vortices. This allows them to survive shear-driven destruction, but they may be prone to elliptical instability and slow decay. When the conditions for growing Rossby wave-related instabilities are maintained in the disc, large-scale vortices can survive over very long timescales and may be able to concentrate solids.

Despite being all roughly of solar composition, primitive meteorites (chondrites) present a diversity in their chemical, isotopic and petrographic properties, and in particular a first-order dichotomy between carbonaceous and non-carbonaceous chondrites. We investigate here analytically the dynamics of their components (chondrules, refractory inclusions, metal/sulfide and matrix grains) in protoplanetary disks prior to their incorporation in chondrite parent bodies. We find the dynamics of the solids, subject to gas drag, to be essentially controlled by the “gas-solid decoupling parameter” $S\equiv \textrm{St}/\alpha$, the ratio of the dimensionless stopping time to the turbulence parameter. The decoupling of the solid particles relative to the gas is significant when $S$ exceeds unity. $S$ is expected to increase with time and heliocentric distance. On the basis of (i) abundance of refractory inclusions (ii) proportion of matrix (iii) lithophile element abundances and (iv) oxygen isotopic composition of chondrules, we propose that non-matrix chondritic components had $S1$ when the other chondrites accreted. This suggests that accretion of carbonaceous chondrites predated on average that of the other chondrites and that refractory inclusions are genetically related to their host carbonaceous chondrites.

In this paper, we present a model for determination of a weakly time dependent effective recombination coefficient for the perturbed terrestrial ionospheric D-region plasma. We study consequences of a class M1.0 X-ray solar flare, recorded by GOES-15 satellite on February 18, 2011 between 14:00 UT and 14:15 UT, by analyzing the amplitude and phase real time variations of very low frequency (VLF) radio waves emitted by transmitter DHO (located in Germany) at frequency 23.4 kHz and recorded by the AWESOME receiver in Belgrade (Serbia). Our analysis is limited to ionospheric perturbations localized at altitudes around 70 km where the dominant electron gain and electron loss processes are the photo-ionization and recombination respectively.

In this paper, we study the influence of solar flares on electron concentration in the terrestrial ionospheric D-region by analyzing the amplitude and phase time variations of very low frequency (VLF) radio waves emitted by DHO transmitter (Germany) and recorded by the AWESOME receiver in Belgrade (Serbia) in real time. The rise of photo-ionization rate in the ionospheric D-region is a typical consequence of solar flare activity as recorded by GOES-15 satellite for the event on March 24, 2011 between 12:01 UT and 12:11 UT. At altitudes around 70 km, the photo-ionization and recombination are the dominant electron gain and electron loss processes, respectively. We analyze the relative contribution of each of these two processes in the resulting electron concentration variation in perturbed ionosphere.

We study large trans-neptunian objects (TNOs) using stellar occultations. We derive precise astrometric predictions for stellar occultations by Eris, Haumea, Ixion, Makemake, Orcus, Quaoar, Sedna, Varuna, 2002 TX300, and 2003 AZ84 for 2011-2015. We construct local astrometric catalogs of stars in the UCAC2 (Second US Naval Observatory CCD Astrograph Catalog) frame covering the sky path of these objects. For that purpose, during 2007-2009, we carried out an observational program at the ESO2p2/WFI (2.2 m Max-Planck ESO telescope with the Wide Field Imager) instrument. Astrometric catalogs with proper motions were produced for each TNO, containing more than 5.35 million stars covering the sky paths with 30′ width in declination. The magnitude completeness is about R = 19 with a limit of about R = 21. We predicted 2717 stellar occultation candidates for all targets. Ephemeris offsets with about from 50 mas to 100 mas precision were applied to each TNO orbit to improve the predictions. They were obtained during 2007-2010 from a parallel observational campaign carried out with from 0.6 m to 2.2 m in size telescopes. This extends our previous work for the Pluto system to large TNOs, using the same observational and astrometric procedures. The obtained astrometric catalogs are useful for follow-up programs at small to large telescopes used to improve the candidate star positions and TNO ephemeris. They also furnish valuable photometric information for the field stars. For each TNO, updates on the ephemeris offsets and candidate star positions (geometric conditions of predictions and finding charts) are made available in the web by the group.

It is very difficult to find archival images of solar system objects. While regular archive searches can find images at a fixed location, they cannot find images of moving targets. Archival images have become increasingly useful to galactic and stellar astronomers the last few years but, until now, solar system researchers have been at a disadvantage in this respect. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for images containing that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, several ESO instruments and the SDSS for a total of 6.5 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos

It is very difficult to find archival images of solar system objects. While regular archive searches can find images at a fixed location, they cannot find images of moving targets. Archival images have become increasingly useful to galactic and stellar astronomers the last few years but, until now, solar system researchers have been at a disadvantage in this respect. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for images containing that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, and several ESO instruments for a total of 1.6 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos