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.

The asteroid belt is an open window on the history of the Solar System, as it preserves records of both its formation process and its secular evolution. The progenitors of the present-day asteroids formed in the Solar Nebula almost contemporary to the giant planets. The actual process producing the first generation of asteroids is uncertain, strongly depending on the physical characteristics of the Solar Nebula, and the different scenarios produce very diverse initial size-frequency distributions. In this work we investigate the implications of the formation of Jupiter, plausibly the first giant planet to form, on the evolution of the primordial asteroid belt. The formation of Jupiter triggered a short but intense period of primordial bombardment, previously unaccounted for, which caused an early phase of enhanced collisional evolution in the asteroid belt. Our results indicate that this Jovian Early Bombardment caused the erosion or the disruption of bodies smaller than a threshold size, which strongly depends on the size-frequency distribution of the primordial planetesimals. If the asteroid belt was dominated by planetesimals less than 100 km in diameter, the primordial bombardment would have caused the erosion of bodies smaller than 200 km in diameter. If the asteroid belt was instead dominated by larger planetesimals, the bombardment would have resulted in the destruction of bodies as big as 500 km.

We present transmission spectra of the hot Jupiter HD 189733b taken with the Space Telescope Imaging Spectrograph aboard HST. The spectra cover the wavelength range 5808-6380 Ang with a resolving power of R=5000. We detect absorption from the NaI doublet within the exoplanet’s atmosphere at the 9 sigma confidence level within a 5 Ang band (absorption depth 0.09 +/- 0.01%) and use the data to measure the doublet’s spectral absorption profile. We detect only the narrow cores of the doublet. The narrowness of the feature could be due to an obscuring high-altitude haze of an unknown composition or a significantly sub-solar NaI abundance hiding the line wings beneath a H2 Rayleigh signature. We compare the spectral absorption profile over 5.5 scale heights with model spectral absorption profiles and constrain the temperature at different atmospheric regions, allowing us to construct a vertical temperature profile. We identify two temperature regimes; a 1280 +/- 240 K region derived from the NaI doublet line wings corresponding to altitudes below ~ 500 km, and a 2800 +/- 400 K region derived from the NaI doublet line cores corresponding to altitudes from ~ 500-4000 km. The zero altitude is defined by the white-light radius of Rp/Rstar=0.15628 +/- 0.00009. The temperature rises with altitude, which is likely evidence of a thermosphere. The absolute pressure scale depends on the species responsible for the Rayleigh signature and its abundance. We discuss a plausible scenario for this species, a high-altitude silicate haze, and the atmospheric temperature-pressure profile that results. In this case, the high altitude temperature rise for HD 189733b occurs at pressures of 10^-5 to 10^-8 bar.

Magnetic fields of low-mass stars and planets are thought to originate from self-excited dynamo action in their convective interiors. Observations reveal a variety of field topologies ranging from large-scale, axial dipole to more structured magnetic fields. In this article, we investigate more than 70 three-dimensional, self-consistent dynamo models obtained by direct numerical simulations. The control parameters, the aspect ratio and the mechanical boundary conditions have been varied to build up this sample of models. Both, strongly dipolar and multipolar models have been obtained. We show that these dynamo regimes can in general be distinguished by the ratio of a typical convective length scale to the Rossby radius. Models with a predominantly dipolar magnetic field were obtained, if the convective length scale is at least an order of magnitude larger than the Rossby radius. Moreover, we highlight the role of the strong shear associated with the geostrophic zonal flow for models with stress-free boundary conditions. In this case, the above transition disappears and is replaced by a region of bistability for which dipolar and multipolar dynamos co-exist. We interpret our results in terms of dynamo eigenmodes using the so-called test-field method. We can thus show that models in the dipolar regime are characterized by an isolated ’single mode’. Competing overtones become significant as the boundary to multipolar dynamos is approached. We discuss how these findings relate to previous models and to observations.

To understand the influence of additional wide stellar companions on planet formation, it is necessary to determine the fraction of multiple stellar systems amongst the known extrasolar planet population. We target recently discovered radial velocity exoplanetary systems observable from the northern hemisphere and with sufficiently high proper motion to detect stellar companions via direct imaging. We utilize the Calar Alto 2.2m telescope in combination with its lucky imaging camera AstraLux. 71 planet host stars have been observed so far, yielding one new low-mass (0.239 \pm 0.022M\odot) stellar companion, 4.5 arcsec (227AU of projected separation) northeast of the planet host star HD185269, detected via astrometry with AstraLux. We also present follow-up astrometry on three previously discovered stellar companions, showing for the first time common proper motion of the 0.5 arcsec companion to HD126614. Additionally, we determined the achieved detection limits for all targets, which allows us to characterize the detection space of possible further companions of these stars.

We present the spectra of 60 asteroids, including 47 V-types observed during the first phase of the Adler V-Type Asteroid (AVAST) Survey. SDSS photometry was used to select candidate V-type asteroids for follow up by nature of their very blue i – z color. 47 of the 61 observed candidates were positively classified as V-type asteroids, while an additional six show indications of a 0.9 {\mu}m feature consistent with V-type spectra, but not sufficient for formal classification. Four asteroids were found to be S-type, all of which had i – z values very near the adopted AVAST selection criteria of i – z \leq -0.2, including one candidate observed well outside the cut (at a mean i – z of -0.11). Three A-type asteroids were also identified. Six V-type asteroids were identified beyond the 3:1 mean motion resonance with Jupiter, and six more were found with low (< 5 deg) inclination, placing these asteroids outside of the normal dynamical range of classic Vestoids, and are suggestive of a non-Vesta origin for at least some of the population.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

We present an extension of the multi-moment advection scheme (Minoshima et al., 2011, J. Comput. Phys.) to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

We present the smoothness of the mid-infrared sky from observations by the Japanese infrared astronomical satellite AKARI. AKARI monitored the north ecliptic pole (NEP) during its cold phase with nine wavebands covering from 2.4 to 24 \mu m, out of which six mid-infrared bands were used in this study. A simple sinusoidal fit to the seasonal variation of the sky brightness shows that the mid-infrared brightness towards the NEP is not affected by small-scale features of the interplanetary dust cloud. We applied power spectrum analysis to the images in order to search for the fluctuation of the sky brightness. Observed fluctuation is explained by fluctuation of photon noise, shot noise of faint sources, and Galactic cirrus. The fluctuations at a few arcminutes scales at short mid-infrared wavelengths (7, 9, and 11 \mum) are largely caused by the diffuse Galactic light of the interstellar dust cirrus. At long mid-infrared wavelengths (15, 18, and 24 \mum), photon noise is the dominant source of fluctuation over the scale from arcseconds to a few arcminutes. The residual fluctuation power at 200″ after removing these contributions is at most 1.08 \pm 0.22 nW m^-2 sr^-1 or 0.05% of the brightness at 24 \mum and at least 0.52 \pm 0.13 nW m^-2 sr^-1 or 0.02% at 18 \mum. We conclude that the upper limit of the fluctuation in the zodiacal light is 0.02% of the sky brightness.

Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy stored in the magnetic field by permitting a change in the magnetic topology. In this article we present a review of the current state of understanding of magnetic reconnection. We discuss theoretical results regarding the formation of current sheets in complex 3D magnetic fields, and describe the fundamental differences between reconnection in two and three dimensions. We go on to outline recent developments in modelling of reconnection with kinetic theory, as well as in the MHD framework where a number of new 3D reconnection regimes have been identified. We discuss evidence from observations and simulations of solar system plasmas that support this theory, and summarise some prominent locations in which this new reconnection theory is relevant in astrophysical plasmas.

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u,g,r,i,z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Af{\rho} formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(< H) \propto 10(0.49\pm0.05)H for H < 18, with evidence of a much shallower fit N(< H) \propto 10(0.19\pm0.03)H for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 \pm 0.05 in g – r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.

Context. Mapping the brightness distribution of exoplanets is the next frontier for exoplanet infrared photometry studies. For tidally-locked hot Jupiters that transit and are eclipsed by their host star with non-zero impact parameter, the first steps are now possible. Aims. The aim is to use eclipse scanning from occultation ingress/egress and phase curve measurements to constrain exoplanet large-scale brightness structure. Methods. We use archived Spitzer/IRAC 8 {\mu}m data of HD189733 in a global MCMC procedure encompassing six transits, eight secondary eclipses, and a phase curve in a two-step analysis. The first step derives the planet-star system parameters. The second step investigates the structure found in eclipse scanning, using the previous planet-star system parameter derivation as Gaussian priors. Results. We find a 5-sigma deviation from the expected occultation ingress/egress shape for a uniform brightness disk, and demonstrate that this is dominated by large-scale brightness structure and not an occultation timing offset due to a non-zero eccentricity. Our analysis yields a 2D brightness temperature distribution showing a large-scale asymmetric hot spot whose finer structure is limited by the data quality and planet orbit geometry. We also present an improved upper limit for eccentricity, e<0.0081 (95% confidence). Conclusions. Reanalysis of archived HD 189733 data revealed brightness structure by using global analysis that mitigated systematics. Future eclipse scanning observations of the same exoplanet at other wavelengths will probe different atmosphere layers, ultimately generating a large-scale 3D map.

Planetary transits provide a unique opportunity to investigate the surface distributions of star spots. Our aim is to determine if, with continuous observation (such as the data that will be provided by the Kepler mission), we can in addition measure the rate of drift of the spot belts. We begin by simulating magnetic cycles suitable for the Sun and more active stars, incorporating both flux emergence and surface transport. This provides the radial magnetic field distribution on the stellar surface as a function of time. We then model the transit of a planet whose orbital axis is misaligned with the stellar rotation axis. Such a planet could occult spots at a range of latitudes. This allows us to complete the forward modelling of the shape of the transit lightcurve. We then attempt the inverse problem of recovering spot locations from the transit alone. From this we determine if transit lightcurves can be used to measure spot belt locations as a function of time. We find that for low-activity stars such as the Sun, the 3.5 year Kepler window is insufficient to determine this drift rate. For more active stars, it may be difficult to distinguish subtle differences in the nature of flux emergence, such as the degree of overlap of the “butterfly wings”. The rate and direction of drift of the spot belts can however be determined for these stars. This would provide a critical test of dynamo theory.

Spitzer IRAC observations of 15 metal-polluted white dwarfs reveal infrared excesses in the spectral energy distributions of HE 0110-5630, GD 61, and HE 1349-2305. All three of these stars have helium-dominated atmospheres, and their infrared emissions are consistent with warm dust produced by the tidal destruction of (minor) planetary bodies. This study brings the number of metal-polluted, helium and hydrogen atmosphere white dwarfs surveyed with IRAC to 53 and 38 respectively. It also nearly doubles the number of metal-polluted helium-rich white dwarfs found to have closely orbiting dust by Spitzer. From the increased statistics for both atmospheric types with circumstellar dust, we derive a typical disk lifetime of log[t_{disk} (yr)] = 5.6+-1.1 (ranging from 3*10^4 – 5*10^6 yr). This assumes a relatively constant rate of accretion over the timescale where dust persists, which is uncertain. We find that the fraction of highly metal-polluted helium-rich white dwarfs that have an infrared excess detected by Spitzer is only 23 per cent, compared to 48 per cent for metal-polluted hydrogen-rich white dwarfs, and we conclude from this difference that the typical lifetime of dusty disks is somewhat shorter than the diffusion time scales of helium-rich white dwarf. We also find evidence for higher time-averaged accretion rates onto helium-rich stars compared to the instantaneous accretion rates onto hydrogen-rich stars; this is an indication that our picture of evolved star-planetary system interactions is incomplete. We discuss some speculative scenarios that can explain the observations.

It has been recently suggested that the magneto-rotational instability (MRI) is a limiting case of the magneto-elliptic instability (MEI). This limit is obtained for horizontal modes in the presence of rotation and an external vertical magnetic field, when the aspect ratio of the elliptic streamlines tends to infinite. In this paper we unveil the link between these previously unconnected mechanisms, explaining both the MEI and the MRI as different manifestations of the same Magneto-Elliptic-Rotational Instability (MERI). The growth rates are found and the influence of the magnetic and rotational effects is explained, in particular the effect of the magnetic field on the range of negative Rossby numbers at which the horizontal instability is excited. Furthermore, we show how the horizontal rotational MEI in the rotating shear flow limit links to the MRI by the use of the local shearing box model, typically used in the study of accretion discs. In such limit the growth rates of the two instability types coincide for any power-type background angular velocity radial profile with negative exponent corresponding to the value of the Rossby number of the rotating shear flow. The MRI requirement for instability is that the background angular velocity profile is a decreasing function of the distance from the centre of the disk which corresponds to the horizontal rotational MEI requirement of negative Rossby numbers. Finally a physical interpretation of the horizontal instability, based on a balance between the strain, the Lorentz force and the Coriolis force is given.

We present Herschel PACS photometry of 18 Plutinos and determine sizes and albedos for these objects using thermal modeling. We analyze our results for correlations, draw conclusions on the Plutino size distribution, and compare to earlier results. Flux densities are derived from PACS mini scan-maps using specialized data reduction and photometry methods. In order to improve the quality of our results, we combine our PACS data with existing Spitzer MIPS data where possible, and refine existing absolute magnitudes for the targets. The physical characterization of our sample is done using a thermal model. Uncertainties of the physical parameters are derived using customized Monte Carlo methods. The correlation analysis is performed using a bootstrap Spearman rank analysis. We find the sizes of our Plutinos to range from 150 to 730 km and geometric albedos to vary between 0.04 and 0.28. The average albedo of the sample is 0.08 \pm 0.03, which is comparable to the mean albedo of Centaurs, Jupiter Family comets and other Trans-Neptunian Objects. We were able to calibrate the Plutino size scale for the first time and find the cumulative Plutino size distribution to be best fit using a cumulative power law with q = 2 at sizes ranging from 120-400 km and q = 3 at larger sizes. We revise the bulk density of 1999 TC36 and find a density of 0.64 (+0.15/-0.11) g cm-3. On the basis of a modified Spearman rank analysis technique our Plutino sample appears to be biased with respect to object size but unbiased with respect to albedo. Furthermore, we find biases based on geometrical aspects and color in our sample. There is qualitative evidence that icy Plutinos have higher albedos than the average of the sample.

The planet-star interaction is manifested in many ways. It was found out that a close-in exoplanet causes small but measurable variability in the cores of a few lines in the spectra of several stars which corresponds to the orbital period of the exoplanet. Stars with and without exoplanets may have different properties. The main goal of our study is to search for influence which exoplanets might have on atmospheres of their host stars. Unlike the previous studies, we do not study changes in the spectrum of a host star or differences between stars with and without exoplanets. We aim to study a large number of stars with exoplanets, current level of their chromospheric activity and look for a possible correlation with the exoplanetary properties. To analyse the chromospheric activity of stars we exploit our own (2.2m ESO/MPG telescope) and publicly available archival spectra (Keck Observatory Archive), measure the equivalent widths of the cores of Ca II H and K lines and use them as a tracer of their activity. Subsequently, we search for their dependence on the orbital parameters and mass of the exoplanet. We found a statistically significant evidence that the equivalent width of the Ca II K line emission and log R’_{HK} activity parameter of the host star varies with the semi-major axis and mass of the exoplanet. Stars with T_eff <= 5500 K having exoplanets with semi-major axis a <= 0.15 AU (P_orb <= 20 days) have a broad range of Ca II K emissions and much stronger emission in general than stars at similar temperatures but with higher values of semi-major axes. Ca II K emission of cold stars (T_eff <= 5500 K) with close-in exoplanets (a <= 0.15 AU) is also more pronounced for more massive exoplanets. The overall level of the chromospheric activity of stars may be affected by their close-in exoplanets. Stars with massive close-in exoplanets may be more active.

The f(R)-gravitational theory with torsion is considered for one family of leptons; it is found that the torsion tensor gives rise to interactions having the structure of the weak forces while the intrinsic non-linearity of the f(R) function provides an energy-dependent coupling: in this way, torsional f(R) gravity naturally generates both structure and strength of the electroweak interactions among leptons. This implies that the weak interactions among the lepton fields could be addressed as a geometric effect due to the interactions among spinors induced by the presence of torsion in the most general f(R) gravity. Phenomenological considerations are addressed.

The stability of the motion of the planet satellites is considered in the model of the general three-body problem (Sun-planet-satellite). “Sundman surfaces” are constructed, by means of which the concept “Sundman stability” is formulated. The comparison of the Sundman stability with the results of Golubev’s c2h method and with the Hill’s classical stability in the restricted three-body problem is performed. The constructed Sundman stability regions in the plane of the parameters “energy – moment of momentum” coincide with the analogous regions obtained by Golubev’s method, with the value (c2h)cr. The construction of the Sundman surfaces in the three-dimensional space of the specially selected coordinates xyR is carried out by means of the exact Sundman inequality in the general three-body problem. The determination of the singular points of surfaces, the regions of the possible motion and Sundman stability analysis are implemented. It is shown that the singular points of the Sundman surfaces in the coordinate space xyR lie in different planes. Sundman stability of all known natural satellites of planets is investigated. It is shown that a number of the natural satellites, that are stable according to Hill and also some satellites that are stable according to Golubev’s method are unstable in the sense of Sundman stability.

We examine the observational consequences of partial gaps being opened by planets in protoplanetary disks. We model the disk using a static alpha-disk model with detailed radiative transfer, parametrizing the shape and size of the partially cleared gaps based on the results of hydrodynamic simulations. Shadowing and illumination by stellar irradiation at the surface of the gap leads to increased contrast as the gap trough is deepened by shadowing and cooling and the far gap wall is puffed up by illumination and heating. In calculating observables, we find that multiple scattering is important and derive an approximation to include these effects. A gap produced by a 200 M_Earth (70 M_Earth) planet at 10 AU can lower/raise the midplane temperature of the disk by up to ~-25/+29% (~-11/+19%) by shadowing in the gap trough and illumination on the far shoulder of the gap. At the distance of Taurus, this gap would be resolvable with ~0.01″ angular resolution. The gap contrast is most significant in scattered light and at thermal continuum wavelengths characteristic of the surface temperature, reducing or raising the surface brightness by up to order of magnitude. Since gaps sizes are correlated to planet mass, this is a promising way of finding and determining the masses of planets embedded in protoplanetary disks.

We report the detection of thermal emission from the hot Jupiter WASP-3b in the KS band, using a newly developed guiding scheme for the WIRC instrument at the Palomar Hale 200in telescope. Our new guiding scheme has improved the telescope guiding precision by a factor of ~5-7, significantly reducing the correlated systematics in the measured light curves. This results in the detection of a secondary eclipse with depth of 0.181%\pm0.020% (9-{\sigma}) – a significant improvement in WIRC’s photometric precision and a demonstration of the capability of Palomar/WIRC to produce high quality measurements of exoplanetary atmospheres. Our measured eclipse depth cannot be explained by model atmospheres with heat redistribution but favor a pure radiative equilibrium case with no redistribution across the surface of the planet. Our measurement also gives an eclipse phase center of 0.5045\pm0.0020, corresponding to an ecos{\omega} of 0.0070\pm0.0032. This result is consistent with a circular orbit, although it also suggests the planet’s orbit might be slightly eccentric. The possible non-zero eccentricity provides insight into the tidal circularization process of the star-planet system, but also might have been caused by a second low-mass planet in the system, as suggested by a previous transit timing variation study. More secondary eclipse observations, especially at multiple wavelengths, are necessary to determine the temperature-pressure profile of the planetary atmosphere and shed light on its orbital eccentricity.

Using semi-analytical, one-dimensional models, we elucidate the influence of scattering and absorption on the degree of Ohmic dissipation in hot Jovian atmospheres. With the assumption of Saha equilibrium, the variation in temperature is the main driver of the variations in the electrical conductivity, induced current and Ohmic power dissipated. Atmospheres possessing temperature inversions tend to dissipate most of the Ohmic power superficially, at high altitudes, whereas those without temperature inversions are capable of greater dissipation deeper down. Scattering in the optical range of wavelengths tends to cool the lower atmosphere, thus reducing the degree of dissipation at depth. Purely absorbing cloud decks (in the infrared), of a finite extent in height, allow for localized reductions in dissipation and may reverse a temperature inversion if they are dense and thick enough, thus greatly enhancing the dissipation at depth. If Ohmic dissipation is the mechanism for inflating hot Jupiters, then variations in the atmospheric opacity (which may be interpreted as arising from variations in metallicity and cloud/haze properties) and magnetic field strength naturally produce a scatter in the measured radii at a given strength of irradiation. Future work will determine if these effects are dominant over evolutionary effects, which also contribute a scatter to the measured radii.

The recent discovery of the first Neptune Trojan at the planet’s trailing (L5) Lagrange point, 2008 LC18, offers an opportunity to confirm the formation mechanism of a member of this important tracer population for the Solar system’s dynamical history. We tested the stability of 2008 LC18’s orbit through a detailed dynamical study, using test particles spread across the \pm3{\sigma} range of orbital uncertainties in a, e, i and {\Omega}. This showed that the wide uncertainties of the published orbit span regions of both extreme dynamical instability, with lifetimes 1 Gyr lifetimes). The stability of 2008 LC18’s clones is greatly dependent on their semi-major axis and only weakly correlated with their eccentricity. Test particles on orbits with an initial semi-major axis less than 29.91 AU have dynamical half-lives shorter than 100 Myr; in contrast, particles with an initial semi-major axis greater than 29.91 AU exhibit such strong dynamical stability that almost all are retained over the 1 Gyr of our simulations. More observations of this object are necessary to improve the orbit. If 2008 LC18 is in the unstable region, then our simulations imply that it is either a temporary Trojan capture, or a representative of a slowly decaying Trojan population (like its sibling the L4 Neptunian Trojan 2001 QR322), and that it may not be primordial. Alternatively, if the orbit falls into the larger, stable region, then 2008 LC18 is a primordial member of the highly stable and highly inclined component of the Neptune Trojan population, joining 2005 TN53 and 2007 VL305. We attempted to recover 2008 LC18 using the 2.3 m telescope at Siding Spring Observatory to provide this astrometry, but were unsuccessful due to the high stellar density of its current sky location near the galactic centre. The recovery of this object will require a telescope in the 8m class.

We explore the change in Jupiter’s normalized axial moment of inertia (NMOI) assuming that Jupiter undergoes core erosion. It is found that Jupiter’s contraction combined with an erosion of 20 M_Earth from a primordial core of 30 M_Earth can change Jupiter’s NMOI over time significantly. It is shown that Jupiter’s NMOI could have changed from ~0.235 to ~0.264 throughout its evolution. We find that a NMOI value of ~0.235 as suggested by dynamical models (Ward & Canup, 2006, ApJ, 640, L91) could, in principle, be consistent with Jupiter’s primordial internal structure. Low NMOI values, however, persist only for the first ~ 10^6 years of Jupiter’s evolution. Re-evaluation of dynamical stability models as well as more sophisticated evolution models of Jupiter with core erosion seem to be required in order to provide more robust estimates for Jupiter’s primordial NMOI.

One of the least understood properties of comets is the compositional structure of their nuclei, which can either be homogeneous or heterogeneous. The nucleus structure can be conveniently studied at millimeter wavelengths, using velocity-resolved spectral time series of the emission lines, obtained simultaneously for multiple molecules as the body rotates. Using this technique, we investigated the sources of CH3OH and HCN in comet 103P/Hartley 2, the target of NASA’s EPOXI mission, which had an exceptionally favorable apparition in late 2010. Our monitoring at IRAM 30-m shows short-term variability of the spectral lines caused by nucleus rotation. The varying production rates generate changes in brightness by a factor of 5 for HCN and by a factor of 2 for CH3OH, and they are remarkably well correlated in time. With the addition of the velocity information from the line profiles, we identify the main sources of outgassing: two jets, oppositely directed in a radial sense, and icy grains, injected into the coma primarily through one of the jets. The mixing ratio of CH3OH and HCN is dramatically different in the two jets, which evidently shows large-scale chemical heterogeneity of the nucleus. We propose a network of identities linking the two jets with morphological features reported elsewhere, and postulate that the chemical heterogeneity may result from thermal evolution. The model-dependent average production rates are 2.10×10**26 molec/s for CH3OH and 1.25×10**25 molec/s for HCN and their ratio is typical of comets. The rotational temperature from CH3OH varied strongly, presumably due to nucleus rotation, with the average value 47 K.

Extrasolar planets and belts of debris orbiting post-main-sequence single stars may become unbound as the evolving star loses mass. In multiple star systems, the presence or co-evolution of the additional stars can significantly complicate the prospects for orbital excitation and escape. Here, we investigate the dynamical consequences of multi-phasic, nonlinear mass loss and establish a criterion for a system of any stellar multiplicity to retain a planet whose orbit surrounds all of the parent stars. For single stars which become white dwarfs, this criterion can be combined with the Chandrasekhar Limit to establish the maximum allowable mass loss rate for planet retention. We then apply the criterion to circumbinary planets in evolving binary systems over the entire stellar mass phase space. Through about 10^5 stellar evolutionary track realizations, we characterize planetary ejection prospects as a function of binary separation, stellar mass and metallicity. This investigation reveals that planets residing at just a few tens of AU from a central concentration of stars are susceptible to escape in a wide variety of multiple systems. Further, planets are significantly more susceptible to ejection from multiple star systems than from single star systems for a given system mass. For system masses greater than about 2 Solar masses, multiple star systems represent the greater source of free-floating planets.

In their model for the origin of the parent bodies of iron meteorites, Bottke et al proposed differentiated planetesimals that were formed in the region of 1-2 AU during the first 1.5 Myr, as the parent bodies, and suggested that these objects and their fragments were scattered into the asteroid belt as a result of interactions with planetary embryos. Although viable, this model does not include the effect of a giant planet that might have existed or been growing in the outer regions. We present the results of a concept study where we have examined the effect of a planetary body in the orbit of Jupiter on the early scattering of planetesimals from terrestrial region into the asteroid belt. We integrated the orbits of a large battery of planetesimals in a disk of planetary embryos, and studied their evolutions for different values of the mass of the planet. Results indicate that when the mass of the planet is smaller than 10 Earth-masses, its effects on the interactions among planetesimals and planetary embryos is negligible. However, when the planet mass is between 10 and 50 Earth-masses, simulations point to a transitional regime with ~50 Earth-mass being the value for which the perturbing effect of the planet can no longer be ignored. Simulations also show that further increase of the mass of the planet strongly reduces the efficiency of the scattering of planetesimals from the terrestrial planet region into the asteroid belt. We present the results of our simulations and discuss their possible implications for the time of giant planet formation.

The recent global fit of short baseline neutrino oscillation data favors the presence of one (or more) sterile neutrino state which leads to new mass splitting \Delta m^2 ~1 eV^2. We consider the effect of this new states on the evolution of neutrinos from the dark matter annihilation inside the Sun. We show that neutrinos with energy E_\nu > 100 GeV undergo resonant active-sterile oscillation which depletes the flux of neutrinos arriving at the Earth. As an example of this effect, we present the oscillation probabilities for the case of monochromatic neutrinos from the direct annihilation of dark matter particles to neutrinos and the depletion due to the presence of sterile neutrinos.

Aims. We study a subset of the planetary population characterized both by HARPS and Kepler surveys. We compare the statistical properties of planets in systems with m.sin i >5-10 M_Earth and R>2 R_Earth. If we assume that the underlying population has the same characteristics, the different detection sensitivity to the orbital inclination relative to the line of sight allows us to probe the planets’ mutual inclination. Methods. We considered the frequency of systems with one, two and three planets as dictated by HARPS data. We used Kepler’s planetary period and host mass and radii distributions (corrected from detection bias) to model planetary systems in a simple yet physically plausible way. We then varied the mutual inclination between planets in a system according to different prescriptions (completely aligned, Rayleigh distributions and isotropic) and compared the transit frequencies with one, two or three planets with those measured by Kepler. Results. The results show that the two datasets are compatible, a remarkable result especially because there are no tunable knobs other than the assumed inclination distribution. For m.sin i cutoffs of 7-10 M_Earth, which are those expected to correspond to the radius cutoff of 2 R_Earth, we conclude that the results are better described by a Rayleigh distribution with mode of 1 deg or smaller. We show that the best-fit scenario only becomes a Rayleigh distribution with mode of 5 deg if we assume a rather extreme mass-radius relationship for the planetary population. Conclusions. These results have important consequences for our understanding of the role of several proposed formation and evolution mechanisms. They confirm that planets are likely to have been formed in a disk and show that most planetary systems evolve quietly without strong angular momentum exchanges (abridged).

We obtained J-, H- and JH-band photometry of known extrasolar planet transiting systems at the 2.1-m Kitt Peak National Observatory Telescope using the FLAMINGOS infrared camera between October 2008 and October 2011. From the derived lightcurves we have extracted the mid-transit times, transit depths and transit durations for these events. The precise mid-transit times obtained help improve the orbital periods and also constrain transit-time variations of the systems. For most cases the published system parameters successfully accounted for our observed lightcurves, but in some instances we derive improved planetary radii and orbital periods. We complemented our 2.1-m infrared observations using CCD z’-band and B-band photometry (plus two Hydrogen Alpha filter observations) obtained with the Kitt Peak Visitor’s Center telescope, and with four H-band transits observed in October 2007 with the NSO’s 1.6-m McMath-Pierce Solar Telescope. The principal highlights of our results are: 1) our ensemble of J-band planetary radii agree with optical radii, with the best-fit relation being: (Rp/R*)J = 0.0017 + 0.979 (Rp/R*)optical, 2) We observe star spot crossings during the transit of WASP-11/HAT-P-10, 3) we detect star spot crossings by HAT-P-11b (Kepler-3b), thus confirming that the magnetic evolution of the stellar active regions can be monitored even after the Kepler mission has ended, and 4) we confirm a grazing transit for HAT-P-27/WASP-40. In total we present 57 individual transits of 32 known exoplanet systems.

In June 2010, we confirmed the existence of a giant planet in the disk of the young star Beta Pictoris, located between 8 AU and 15 AU from the star. This young planet offers the rare opportunity to monitor a large fraction of the orbit using the imaging technique over a reasonably short timescale. Using the NAOS-CONICA adaptive-optics instrument (NACO) at the Very Large Telescope (VLT), we obtained repeated follow-up images of the Bpic system in the Ks and L’ filters at four new epochs in 2010 and 2011. Complementing these data with previous measurements, we conduct a homogeneous analysis, which covers more than eight yrs, to accurately monitor the Bpic b position relative to the star. On the basis of the evolution of the planet’s relative position with time, we derive the best-fit orbital solutions for our measurements. More reliable results are found with a Markov-chain Monte Carlo approach. The solutions favor a low-eccentricity orbit e < 0.17, with semi-major axis in the range 8–9 AU corresponding to orbital periods of 17–21 yrs. Our solutions favor a highly inclined solution with a peak around i=88.5+-1.7 deg, and a longitude of ascending node tightly constrained at Omega = -147.5+-1.5 deg. These results indicate that the orbital plane of the planet is likely to be above the midplane of the main disk, and compatible with the warp component of the disk being tilted between 3.5 deg and 4.0 deg. This suggests that the planet plays a key role in the origin of the inner warped-disk morphology of the Bpic disk. Finally, these orbital parameters are consistent with the hypothesis that the planet is responsible for the transit-like event observed in November 1981, and also linked to the cometary activity observed in the Bpic system.

We generalize the problem of the semi-gray model to cases in which a non-negligible fraction of the stellar radiation falls on the long-wavelength range, and/or that the planetary long-wavelength emission penetrates into the transparent short wavelength domain of the absorption. Second, applying the most general assumptions and independently of any particular properties of an absorber, we show that the greenhouse effect saturates and any Earth-like planet has a maximal temperature which depends on the type of and distance to its main-sequence star, its albedo and the primary atmospheric components which determine the cutoff frequency below which the atmosphere is optically thick. For example, a hypothetical convection-less planet similar to Venus, that is optically thin in the visible, could have at most a surface temperature of 1200-1300K irrespective of the nature of the greenhouse gas. We show that two primary mechanisms are responsible for the saturation of the runaway greenhouse effect, depending on the value of the wavelength above which the atmosphere becomes optically thick. Unless this wavelength is small and resides in the optical region, saturation is achieved by radiating the thermal flux of the planet through the short wavelength tail of the thermal distribution. This has the observational implication, the radiation from such a planet should be skewed towards the NIR. Otherwise, saturation takes place by radiating through windows in the FIR.

Context. The {\beta} Pictoris system with its debris disk and a massive giant planet orbiting at \simeq 9 AU represents an ideal laboratory to study giant planet formation and evolution as well as planet-disk interactions. {\beta} Pic b can also help testing brightness-mass relations at young ages. Other planets, yet undetected, may of course be present in the system. Aims. We aim at putting direct constrains on the mass of {\beta} Pic b and at searching for additional jovian planets on orbits closer than typically 2 AU. Methods. We use high precision Harps data collected over 8 years since 2003 to measure and analyse {\beta} Pic radial velocities. Results. We show that the true mass of {\beta} Pic b is less than 10, 12, 15.5, 20 and 25 MJup if orbiting respectively at 8, 9, 10, 11 and 12 AU. This is the first direct constraint on the mass of an imaged planet. The upper mass found is well in the range predicted by brightness-mass relations provided by current “hot start” models. We also exclude the presence of giant planets more massive than 2.5 MJup with periods less than 100 days (hot Jupiters), more massive than 9 MJup for periods in the range 100-500 days. In the 500-1000 day range, the detection limit is in the brown dwarf domain. Beyond the intrinsic interest for {\beta} Pic, these results show the possibilities of precise RV measurements of early type, rapidly rotating stars.

Doppler spectroscopy has uncovered or confirmed all the known planets orbiting nearby stars. Two different approaches are used to obtain precision Doppler measurements at optical wavelengths. The first approach is the gas cell method, which is based on the least-squares matching of the absorption spectrum of Iodine over-imposed to the spectrum of the star. The second method relies on the construction of a stabilized spectrograph calibrated in wavelength with an externally fed calibration source. The most precise stabilized spectrometer in operation is HARPS, operated by ESO in La Silla/Chile. In the case of HARPS, the Doppler measurement is obtained using the so–called Cross-Correlation Function technique (CCF). It consists of multiplying the stellar spectrum with a binary mask and finding the minimum of such product as a function of the stellar Doppler shift. Such mask is weighted to account for the different depths of the stellar lines. It is known that CCF is suboptimal in exploiting the Doppler information in the stellar spectrum. Here, we describe an algorithm to obtain precision RV measurements based on least squares matching of each observation to a high signal-to-noise ratio template. Such algorithm is implemented in our software called HARPS-TERRA (Template Enhanced Radial velocity Re-analysis Application). We show that, compared to CCF, template matching provides a significant improvement in accuracy, specially when applied to M dwarfs. We conclude that other stabilized spectrographs should use a similar approach to achieve the sub \ms precision required to detect potentially habitable worlds around nearby stars.

M dwarfs have been found to often have super-Earth planets with short orbital periods. Such stars are thus preferential targets in searches for rocky or ocean planets in the solar neighbourhood. In a recent paper (Bonfils et al. 2011), we announced the discovery of respectively 1 and 2 low mass planets around the M1.5V stars Gl433 and Gl667C. We found those planets with the HARPS spectrograph on the ESO~3.6-m telescope at La Silla Observatory, from observations obtained during the Guaranteed Time Observing program of that instrument. We have obtained additional HARPS observations of those two stars, for a total of respectively 67 and 179 Radial Velocity measurements for Gl433 and Gl667C, and present here an orbital analysis of those extended data sets and our main conclusion about both planetary systems. One of the three planets, Gl667Cc, has a mass of only M2.sin(i)~4.25 M_earth and orbits in the central habitable zone of its host star. It receives just 10% less stellar energy from Gl667C than the Earth receives from the Sun. However planet evolution in habitable zone can be very different if the host star is a M dwarf or a solar-like star, without necessarily questioning the presence of water. The two other planets, Gl433b and Gl667Cb, both have M2.sin(i) of ~5.5 M_earth and periods of ~7 days. The Radial Velocity measurements of both stars contain longer time scale signals, which we fit as longer period Keplerians. For Gl433 that signal probably originates in a Magnetic Cycle, while a longer time span will be needed to conclude for Gl667C. The metallicity of Gl433 is close to solar, while Gl667C is metal poor with [Fe/H] ~ -0.6. This reinforces the recent conclusion that the occurence of Super-Earth planets does not strongly correlate with stellar metallicity.

We present a reanalysis of the relationship between asteroid albedo and polarization properties using the albedos derived from the Wide-field Infrared Survey Explorer. We find that the function that best describes this relation is a three-dimensional linear fit in the space of log(albedo)-log(polarization slope)-log(minimum polarization). When projected to two dimensions the parameters of the fit are consistent with those found in previous work. We also define p* as the quantity of maximal polarization variation when compared with albedo and present the best fitting albedo-p* relation. Some asteroid taxonomic types stand out in this three-dimensional space, notably the E, B, and M Tholen types, while others cluster in clumps coincident with the S- and C-complex bodies. We note that both low albedo and small (D<30 km) asteroids are under-represented in the polarimetric sample, and we encourage future polarimetric surveys to focus on these bodies.

The Habitable Zone Gallery (www.hzgallery.org) is a new service to the exoplanet community which provides Habitable Zone (HZ) information for each of the exoplanetary systems with known planetary orbital parameters. The service includes a sortable table with information on the percentage of orbital phase spent within the HZ, planetary effective temperatures, and other basic planetary properties. In addition to the table, we also plot the period and eccentricity of the planets with respect to their time spent in the HZ. The service includes a gallery of known systems which plot the orbits and the location of the HZ with respect to those orbits. Also provided are animations which aid in orbit visualization and provide the changing effective temperature for those planets in eccentric orbits. Here we describe the science motivation, the under-lying calculations, and the structure of the web site.

In recent years several planets have been discovered at wide orbits (>100 AU) around their host stars. Theoretical studies encounter difficulties in explaining their formation and origin. Here we propose a novel scenario for the production of planetary systems at such orbits, through the dynamical recapture of free floating planets (FFPs) in dispersing stellar clusters. This process is a natural extension of the recently suggested scenario for the formation of wide stellar binaries. We use N-body simulations of dispersing clusters with 10-1000 and f_FFP=0.5-2 to study this process. We find that planets are captured into wide orbits, ~100-10^6 AU, and a thermal eccentricity distribution. Typically, 3-6x(f_FFP/1) % of all stars capture a planetary companion (f_FFP is the number of FFP per star). The planetary capture efficiency is comparable to that of capture-formed stellar-binaries, and shows a similar dependence on the cluster size and structure. The capture efficiency is almost independent of the specific planetary mass; planets as well as sub-stellar companions of any mass can be captured, where the capture efficiency decreases with increasing cluster size. For a given cluster size the capture efficiency increases with the host/primary mass. More than one planet can be captured around the same host, and planets can be captured into binary systems. We also expect planets to be captured into pre-existing planetary systems as well as around compact objects, if these formed early enough before the cluster dispersal. In particular, stellar black holes have a high capture efficiency (>50 % and 5-10x(f_FFP/1) % for capture of stars and planetary companions, respectively) due to their large mass. Finally, although rare, two FFPs or brown dwarfs can become bound and form a FFP-binary system with no stellar host through this process.

Recently the Cassini Ultraviolet Imaging Spectrograph has revealed the presence of N2 Vegard-Kaplan band emissions in Titan’s dayglow limb observation. We present model calculations for the production of various N2 triplet states in the upper atmosphere of Titan. The Analytical Yield Spectra technique is used to calculate steady state photoelectron fluxes in Titan’s atmosphere, which are in agreement with those observed by the Cassini’s CAPS instrument. Considering direct electron impact excitation, inter-state cascading, and quenching effects, the population of different levels of N2 triplet states are calculated under statistical equilibrium. Densities of all vibrational levels of each triplet state and volume production rates for various triplet states are calculated in the model. Vertically integrated overhead intensities for the same date and lighting conditions as the reported by UVIS observations for N2 VK, 1P, 2P, Wu-Benesch, and Reverse First Positive bands of N2 are found to be 132, 114, 19, 22, and 22 R, respectively. Overhead intensities are calculated for each vibrational transition of all the triplet band emissions of N2, which span a wider spectrum of wavelengths from ultraviolet to infrared. The calculated limb intensities of total and prominent transitions of VK band are presented. The model limb intensity of VK emission within the 150-190 nm wavelength region is in good agreement with the Cassini UVIS observed limb profile. An assessment of the impact of solar EUV flux on the N2 triplet band emission intensity has been made by using three different solar flux models, viz., Solar EUV Experiment, SOLAR2000 model of Tobiska (2004), and HEUVAC model of Richards et al, (2006). The calculated N2 VK band intensity at the peak of limb intensity due to S2K and HEUVAC solar flux models is a factor of 1.2 and 0.9, respectively, of that obtained using SEE solar EUV flux.

We explore the correlation between an asteroid’s taxonomy and photometric phase curve using the H, G12 photometric phase function, with the shape of the phase function described by the single parameter G12. We explore the usability of G12 in taxonomic classification for individual objects, asteroid families, and dynamical groups. We conclude that the mean values of G12 for the considered taxonomic complexes are statistically different, and also discuss the overall shape of the G12 distribution for each taxonomic complex. Based on the values of G12 for about half a million asteroids, we compute the probabilities of C, S, and X complex membership for each asteroid. For an individual asteroid, these probabilities are rather evenly distributed over all of the complexes, thus preventing meaningful classification. We then present and discuss the G12 distributions for asteroid families, and predict the taxonomic complex preponderance for asteroid families given the distribution of G12 in each family. For certain asteroid families, the probabilistic prediction of taxonomic complex preponderance can clearly be made. The Nysa-Polana family shows two distinct regions in the proper element space with different G12 values dominating in each region. We conclude that the G12-based probabilistic distribution of taxonomic complexes through the main belt agrees with the general view of C complex asteroid proportion increasing towards the outer belt. We conclude that the G12 photometric parameter cannot be used in determining taxonomic complex for individual asteroids, but it can be utilized in the statistical treatment of asteroid families and different regions of the main asteroid belt.

(Abridged) We access the expected EKB dust disk properties by modeling. We treat the debiased population of the known transneptunian objects (TNOs) as parent bodies and generate the dust with our collisional code. The resulting dust distributions are modified to take into account the influence of gravitational scattering and resonance trapping by planets on migrating dust grains as well as the effect of sublimation. A difficulty is that the amount and distribution of dust are largely determined by sub-kilometer-sized bodies. These are directly unobservable, and their properties cannot be accessed by collisional modeling, because objects larger than 10…60m in the present-day EKB are not in a collisional equilibrium. To place additional constraints, we use in-situ measurements of the New Horizons spacecraft within 20AU. We show that the TNO population has to have a break in the size distribution at s10AU)~10^-6 and their fractional luminosity is f_d~10^-7. Planets and sublimation are found to have little effect on dust impact fluxes and dust thermal emission. The spectral energy distribution of an EKB analog, as would be seen from 10pc distance, peaks at wavelengths of 40…50\mum at F~0.5mJy, which is less than 1% of the photospheric flux at those wavelengths. Therefore, exact EKB analogs cannot be detected with present-day instruments such as Herschel/PACS.

This study determined the contribution of Martian topography and the density jump at the Mohorovicic discontinuity (M) to the gravity in a quadratic approximation. It also resolved the problem of determining the possible depths of compensation for topography harmonics of various degrees and orders. It shows that all the topography compensation is within the depth range from 0 to 1400 km. Different topographic irregularities are most likely to be compensated at depths that correspond to the upper crust (d =(4.5 +- 3.7)km), crust-mantle transition layer (d = (78 +- 24)km), lithospheric boundary (d = (200 +- 34)km), upper-middle mantle transition layer (d = (400 +- 70)km), or middle-lower mantle transition layer (d =(1120 +- 180)km). This paper presents the lateral distributions of compensation masses for these depths and the respective maps. According to calculations, stresses in the Martian crust and mantle may be as high as 10^8 Pa. This paper shows that the topographic anomalies of the Tharsis volcanic plateau and the symmetric formation in the eastern hemisphere could have originated and be dynamically maintained by two plumes of melted mantle substance enriched with fluids; these plumes may have their origin at the boundary of the lower mantle.

The eigenvalues and eigenfunctions of a linear {\alpha}^{2}-dynamo have been computed for different spatial distributions of an isotropic \alpha-effect. Oscillatory solutions are obtained when \alpha exhibits a sign change in the radial direction. The time-dependent solutions arise at so called exceptional points where two stationary modes merge and continue as an oscillatory eigenfunction with conjugate complex eigenvalues. The close proximity of oscillatory and non-oscillatory solutions may serve as the basic ingredient for reversal models that describe abrupt polarity switches of a dipole induced by noise. Whereas the presence of an inner core with different magnetic diffusivity has remarkable little impact on the character of the dominating dynamo eigenmodes, the introduction of equatorial symmetry breaking considerably changes the geometric character of the solutions. Around the dynamo threshold the leading modes correspond to hemispherical dynamos even when the symmetry breaking is small. This behavior can be explained by the approximate dipole-quadrupole degeneration for the unperturbed problem. More complicated scenarios may occur in case of more realistic anisotropies of \alpha- and \beta-effect or through non-linearities caused by the back-reaction of the magnetic field (magnetic quenching).

The eigenvalues and eigenfunctions of a linear {\alpha}^{2}-dynamo have been computed for different spatial distributions of an isotropic \alpha-effect. Oscillatory solutions are obtained when \alpha exhibits a sign change in the radial direction. The time-dependent solutions arise at so called exceptional points where two stationary modes merge and continue as an oscillatory eigenfunction with conjugate complex eigenvalues. The close proximity of oscillatory and non-oscillatory solutions may serve as the basic ingredient for reversal models that describe abrupt polarity switches of a dipole induced by noise. Whereas the presence of an inner core with different magnetic diffusivity has remarkable little impact on the character of the dominating dynamo eigenmodes, the introduction of equatorial symmetry breaking considerably changes the geometric character of the solutions. Around the dynamo threshold the leading modes correspond to hemispherical dynamos even when the symmetry breaking is small. This behavior can be explained by the approximate dipole-quadrupole degeneration for the unperturbed problem. More complicated scenarios may occur in case of more realistic anisotropies of \alpha- and \beta-effect or through non-linearities caused by the back-reaction of the magnetic field (magnetic quenching).

I briefly review the current state of the modelling of photospheric activity based on the high-precision optical light curves obtained with MOST, CoRoT, and Kepler. These models can be used to search for active longitudes where activity is preferentially concentrated, estimate the amplitude of stellar differential rotation, and look for short-term activity cycles as, e.g., in the case of CoRoT-2. In the case of a late-type star accompanied by a transiting hot Jupiter, the small light modulations observed during transits when a dark spot is occulted by the disc of the planet are also briefly considered. They can be used to derive information on individual active regions as well as on stellar rotation and the spin-orbit alignment of the system.

Due to the importance of collisions and impacts in early phases of the evolution of the planetary system, it is interesting to estimate the heating of a solid target due to an impact in it . A physically simple calculation of the temperature to which a solid target heats up after the impact of a projectile with mass $m$ and speed $v$ is performed,and possibilities for the application of this result in planetology are pointed out.

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

Relations between the length of a sunspot cycle and the average temperature in the same and the next cycle are calculated for a number of meteorological stations in Norway and in the North Atlantic region. No significant trend is found between the length of a cycle and the average temperature in the same cycle, but a significant negative trend is found between the length of a cycle and the temperature in the next cycle. This provides a tool to predict an average temperature decrease of at least 1.0 “C from solar cycle 23 to 24 for the stations and areas analyzed. We find for the Norwegian local stations investigated that 25-56% of the temperature increase the last 150 years may be attributed to the Sun. For 3 North Atlantic stations we get 63-72% solar contribution. This points to the Atlantic currents as reinforcing a solar signal.

We derive the first secondary eclipse map of an exoplanet, HD 189733b, based on Spitzer IRAC 8 micron data. We develop two complementary techniques for deriving the two dimensional planet intensity: regularized slice mapping and spherical harmonic mapping. Both techniques give similar derived intensity maps for the infrared day-side flux of the planet, while the spherical harmonic method can be extended to include phase variation data which better constrain the map. The longitudinal offset of the day-side hot spot is consistent with that found in prior studies, strengthening the claim of super-rotating winds, and eliminating the possibility of phase variations being caused by stellar variability. The latitude of the hot-spot is within 12.5 deg (68% confidence) of the planet’s equator, confirming the predictions of general circulation models for hot Jupiters and indicative of a small planet obliquity.

The torques exerted by a locally isothermal disk on an embedded planet lead to rapid inward migration. Recent work has shown that modeling the thermodynamics without the assumption of local isothermality reveals regions where the net torque on an embedded planet is positive, leading to outward migration of the planet. When a region with negative torque lies directly exterior to this, planets in the inner region migrate outwards and planets in the outer region migrate inwards, converging where the torque is zero. We incorporate the torques from an evolving non-isothermal disk into an N-body simulation to examine the behavior of planets or planetary embryos interacting in the convergence zone. We find that mutual interactions do not eject objects from the convergence zone. Small numbers of objects in a laminar disk settle into near resonant orbits that remain stable over the 10 Myr periods that we examine. However, either or both increasing the number of planets or including a correlated, stochastic force to represent turbulence drives orbit crossings and mergers in the convergence zone. These processes can build gas giant cores with masses of order ten Earth masses from sub-Earth mass embryos in 2-3 Myr.

The torques exerted by a locally isothermal disk on an embedded planet lead to rapid inward migration. Recent work has shown that modeling the thermodynamics without the assumption of local isothermality reveals regions where the net torque on an embedded planet is positive, leading to outward migration of the planet. When a region with negative torque lies directly exterior to this, planets in the inner region migrate outwards and planets in the outer region migrate inwards, converging where the torque is zero. We incorporate the torques from an evolving non-isothermal disk into an N-body simulation to examine the behavior of planets or planetary embryos interacting in the convergence zone. We find that mutual interactions do not eject objects from the convergence zone. Small numbers of objects in a laminar disk settle into near resonant orbits that remain stable over the 10 Myr periods that we examine. However, either or both increasing the number of planets or including a correlated, stochastic force to represent turbulence drives orbit crossings and mergers in the convergence zone. These processes can build gas giant cores with masses of order ten Earth masses from sub-Earth mass embryos in 2-3 Myr.

The Kilodegree Extremely Little Telescope (KELT) project is a survey for new transiting planets around bright stars. KELT-South is a small-aperture, wide-field automated telescope located at Sutherland, South Africa. The telescope surveys a set of 26 degree by 26 degree fields around the southern sky, and targets stars in the range of 8 < V < 10 mag, searching for transits by Hot Jupiters. This paper describes the KELT-South system hardware and software and discusses the quality of the observations. We show that KELT-South is able to achieve the necessary photometric precision to detect transits of Hot Jupiters around solar-type main-sequence stars.

Multi-wavelength observations of protoplanetary disks in the sub-millimeter continuum have measured spectral indices values which are significantly lower than what is found in the diffuse interstellar medium. Under the assumption that mm-wave emission of disks is mostly optically thin, these data have been generally interpreted as evidence for the presence of mm/cm-sized pebbles in the disk outer regions. In this work we investigate the effect of possible local optically thick regions on the mm-wave emission of protoplanetary disks without mm/cm-sized grains. A significant local increase of the optical depth in the disk can be caused by the concentration of solid particles, as predicted to result from a variety of proposed physical mechanisms. We calculate the filling factors and implied overdensities these optically thick regions would need to significantly affect the millimeter fluxes of disks, and we discuss their plausibility. We find that optically thick regions characterized by relatively small filling factors can reproduce the mm-data of young disks without requesting emission from mm/cm-sized pebbles. However, these optically thick regions require dust overdensities much larger than what predicted by any of the physical processes proposed in the literature to drive the concentration of solids. We find that only for the most massive disks it is possible and plausible to imagine that the presence of optically thick regions in the disk is responsible for the low measured values of the mm spectral index. For the majority of the disk population, optically thin emission from a population of large mm-sized grains remains the most plausible explanation. The results of this analysis further strengthen the scenario for which the measured low spectral indices of protoplanetary disks at mm wavelengths are due to the presence of large mm/cm-sized pebbles in the disk outer regions.

We present an analysis of the KOI-13.01 candidate exoplanet system included in the September 2011 Kepler data release. The host star is a known and relatively bright $(m_{\rm KP} = 9.95)$ visual binary with a separation significantly smaller (0.8 arcsec) than the size of a Kepler pixel (4 arcsec per pixel). The Kepler light curve shows both primary and secondary eclipses, as well as significant out-of-eclipse light curve variations. We confirm that the transit occurs round the brighter of the two stars. We model the relative contributions from (i) thermal emission from the companion, (ii) planetary reflected light, (iii) Doppler beaming, and (iv) ellipsoidal variations in the host-star arising from the tidal distortion of the host star by its companion. Our analysis, based on the light curve alone, enables us to constrain the mass of the KOI-13.01 companion to be $M_{\rm C} = 8.3 \pm 1.25M_{\rm J}$ and thus demonstrates that the transiting companion is a planet (rather than a brown dwarf which was recently proposed by \cite{b7}). The high temperature of the host star (Spectral Type A5-7V, $T_{\rm eff} = 8511-8020$ K), combined with the proximity of its companion KOI-13.01, may make it one of the hottest exoplanets known, with a detectable thermal contribution to the light curve even in the Kepler optical passband. However, the single passband of the Kepler light curve does not enable us to unambiguously distinguish between the thermal and reflected components of the planetary emission. Infrared observations may help to break the degeneracy, while radial velocity follow-up with $\sigma \sim$ 100 m s$^{-1}$ precision should confirm the mass of the planet.

Due to their low mass and luminosity, M dwarfs are ideal targets if one hopes to find low-mass planets similar to Earth by using the radial velocity (RV) method. However, stellar magnetic cycles could add noise or even mimic the RV signal of a long-period companion. Following our previous work that studied the correlation between activity cycles and long-term RV variations for K dwarfs we now expand that research to the lower-end of the main sequence. Our objective is to detect any correlations between long-term activity variations and the observed RV of a sample of M dwarfs. We used a sample of 27 M-dwarfs with a median observational timespan of 5.9 years. The cross-correlation function (CCF) with its parameters RV, bisector inverse slope (BIS), full-width-at-half- maximum (FWHM) and contrast have been computed from the HARPS spectrum. The activity index have been derived using the Na I D doublet. These parameters were compared with the activity level of the stars to search for correlations. We detected RV variations up to ~5 m/s that we can attribute to activity cycle effects. However, only 36% of the stars with long-term activity variability appear to have their RV affected by magnetic cycles, on the typical timescale of ~6 years. Therefore, we suggest a careful analysis of activity data when searching for extrasolar planets using long-timespan RV data.

We report the results of an analysis of all Spitzer/MIPS 24 micron observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether we analyze 2.5 transits, 3 eclipses, and a 58-hour near-continuous observation designed to detect the planet’s thermal phase curve. The results of our analysis are: (1) A mean transit depth of 1.484% +/- 0.035%, consistent with previous measurements and showing no evidence of variability in transit depth at the 3% level. (2) A mean eclipse depth of 0.332% +/- 0.026%, somewhat higher than that previously reported for this system; this new value brings observations into better agreement with models. The dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the system occur 32 s +/- 129 s earlier than would be expected from a circular orbit, which constrains the orbital quantity (e cos omega) to be 0.00004 +/- 0.00033. This result is fully consistent with a circular orbit and sets an upper limit of 140 m/s (3 sigma) on any eccentricity-induced velocity offset during transit. The phase curve observations (including one of the transits) exhibit an anomalous trend similar to the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover the system parameters. The photometry which follows the ramp and transit exhibits a gradual, ~0.2% decrease in flux, similar to that seen in pre-launch calibration data. The large uncertainties associated with this poorly-understood, likely instrumental effect prevent us from usefully constraining the planet’s thermal phase curve. Our observations highlight the need for a thorough understanding of detector-related instrumental effects on long time scales when making the high-precision mid-infrared measurements planned for future missions such as EChO, SPICA, and JWST. [abridged]

Astrophysical reconnection takes place in a turbulent medium. The turbulence in most cases is pre-existing, not caused by the reconnection itself. The model of magnetic reconnection in Lazarian & Vishniac (1999) predicts that in the presence of turbulence the reconnection becomes fast, i.e. it is independent of resistivity, but dependent on the level of turbulence. Magnetic reconnection injects energy into plasmas through a turbulent outflow from the reconnection region and this outflow can enhance the level of turbulence creating bursts of reconnection. Magnetic reconnection in the presence of turbulence can accelerate energetic particles through the first order Fermi mechanism, as was discussed in Gouveia dal Pino & Lazarian (2005). We discuss new numerical results on particle acceleration in turbulent reconnection, compare the acceleration arising from turbulent reconnection to the acceleration of energetic particles in turbulent medium.

Colors of Trans Neptunian Objects (TNOs) are used to study the evolutionary processes of bodies in the outskirts of the Solar System, and to test theories regarding their origin. Here I describe a search for serendipitous Sloan Digital Sky Survey (SDSS) observations of known TNOs and Centaurs. I present a catalog of SDSS photometry, colors and astrometry of 388 measurements of 42 outer Solar-System objects. I find a weak evidence, at the ~2-sigma level (per trial), for a correlation between the g-r color and inclination of scattered disk objects and hot classical KBOs. I find a correlation between the g-r color and the angular momentum in the z direction of all the objects in this sample. Light curves as a function of phase angle are constructed for 13 objects. The steepness of the slopes of these light curves suggests that the coherent backscatter mechanism plays a major role in the reflectivity of outer Solar-System small objects at small phase angles. I find a weak evidence for an anti-correlation, significant at the 2-sigma confidence level (per trial), between the g-band phase angle slope parameter and the semi-major axis, as well as the aphelion distance, of these objects. I discuss the origin of this possible correlation and argue that if this correlation is real it probably indicates that “Sedna”-like objects have a different origin than other classes of TNOs. Finally, I identify several objects with large variability amplitudes (abridged).

Physical characterization of Trans-Neptunian objects, a primitive population of the outer solar system, may provide constraints on their formation and evolution. The goal of this work is to characterize a set of 15 scattered disk (SDOs) and detached objects, in terms of their size, albedo, and thermal properties. Thermal flux measurements obtained with the Herschel-PACS instrument at 70, 100 and 160 \mu m, and whenever applicable, with Spitzer-MIPS at 24 and 70 \mu m, are modeled with radiometric techniques, in order to derive the objects’ individual size, albedo and when possible beaming factor. Error bars are obtained from a Monte-Carlo approach. We look for correlations between these and other physical and orbital parameters. Diameters obtained for our sample range from 100 to 2400 km, and the geometric albedos (in V band) vary from 3.8 % to 84.5 %. The unweighted mean V geometric albedo for the whole sample is 11.2 % (excluding Eris); 6.9 % for the SDOs, and 17.0 % for the detached objects (excluding Eris). We obtain new bulk densities for three binary systems: Ceto/Phorcys, Typhon/Echidna and Eris/Dysnomia. Apart from correlations clearly due to observational bias, we find significant correlations between albedo and diameter (more reflective objects being bigger), and between albedo, diameter and perihelion distance (brighter and bigger objects having larger perihelia). We discuss possible explanations for these correlations.

It has long been thought that the presence of a giant planet is a pre-requisite for the development of life on potentially habitable planets. Without Jupiter, it was argued, the Earth would have been subject to a punishing impact regime, which would have significantly retarded or outright prevented the development of life on our planet. Although this idea is widely embraced, little research has previously been carried out to support it. Here, we present the results of several suites of dynamical integrations used to model the influence of Jupiter’s mass and orbit on the impact rate that would be experienced by the Earth. We find that, far from being a simple shield, Jupiter’s role in determining the terrestrial impact flux is significantly more complicated than previously thought. Far from being a simple friend, such giant planets are perhaps more likely to imperil the development of life on otherwise habitable planets.

Radio detection provides unique means to measure and study magnetic fields of the coolest brown dwarfs. Previous radio surveys have observed quiescent and flaring emission from brown dwarfs down to spectral type L3.5, but only upper limits have been established for even cooler objects. We report the detection of sporadic, circularly polarized flares from the T6.5 dwarf, 2MASS J1047+21, with the Arecibo radio telescope at 4.75 GHz. This is by far the coolest brown dwarf yet detected at radio frequencies. The fact that such an object is capable of generating observable, coherent radio emission, despite its very low, ~900 K temperature, demonstrates the feasibility of studies of brown dwarfs in the meagerly explored LTY spectral range, using radio detection as a tool.

Radio detection provides unique means to measure and study magnetic fields of the coolest brown dwarfs. Previous radio surveys have observed quiescent and flaring emission from brown dwarfs down to spectral type L3.5, but only upper limits have been established for even cooler objects. We report the detection of sporadic, circularly polarized flares from the T6.5 dwarf, 2MASS J1047+21, with the Arecibo radio telescope at 4.75 GHz. This is by far the coolest brown dwarf yet detected at radio frequencies. The fact that such an object is capable of generating observable, coherent radio emission, despite its very low, ~900 K temperature, demonstrates the feasibility of studies of brown dwarfs in the meagerly explored LTY spectral range, using radio detection as a tool.

This work deals with the comparison of different parameters (French/IMCCE and US/JPL ephemerides) used to calculate the extension of the umbral shadow and the location of the centre line in the total solar eclipse that took place on March, 29th 2006. We needed to know the exact points where the predicted duration of the solar eclipse was 0 seconds. This phenomenon allows, observed on this particular position, the measurement of the solar diameter. The work is composed by three parts. In the first one is described the evolution of the lunar theory from Hipparchus to the current theory (Bureau de Longitudes/IMCCE, Chapront). The second part is dedicated to geodesy (WGS84 datum) and there are examples of longitude’s calculations. In the third one there is the analysis of the parameters to understand the reasons about the different extensions of the umbral shadow and the different positions of the centre line in the solar eclipse calculated respectively with French and American ephemerides. The use of the center of the figure of the Moon with respect to the center of mass produces the differences in the ephemerides. This work prepared the observational mission in Egypt to observe the total solar eclipse from the shadow’s edges in Zawayet al Mahtallah, near Sidi Barrani, whose results were after published in Solar Physics vol. 259 p. 189 (2009).

The Alpha Centauri AB system is an attractive one for radial velocity observations to detect potential exoplanets. The high metallicity of both Alpha Centauri A and B suggest that they could have possessed circumstellar discs capable of forming planets. As the closest star system to the Sun, with well over a century of accurate astrometric measurements (and Alpha Centauri B exhibiting low chromospheric activity) high precision surveys of Alpha Centauri B’s potential exoplanetary system are possible with relatively cheap instrumentation. Authors studying habitability in this system typically adopt habitable zones (HZs) based on global radiative balance models that neglect the radiative perturbations of Alpha Centauri A. We investigate the habitability of planets around Alpha Centauri B using 1D latitudinal energy balance models (LEBMs), which fully incorporate the presence of Alpha Centauri A as a means of astronomically forcing terrestrial planet climates. We find that the extent of the HZ is relatively unchanged by the presence of Alpha Centauri A, but there are variations in fractional habitability for planets orbiting at the boundaries of the zone due to Alpha Centauri A, even in the case of zero eccentricity. Temperature oscillations of a few K can be observed at all planetary orbits, the strength of which varies with the planet’s ocean fraction and obliquity.

Recently the first Earth Trojan has been observed (Mainzer et al., ApJ 731) and found to be on an interesting orbit close to the Lagrange point L4 (Connors et al., Nature 475). In the present study we therefore perform a detailed investigation on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our results are derived using different approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the phase space geometry, including the indirect influence of the additional planets of our Solar system. b) We use precise numerical methods to integrate the orbit forward and backward in time in different dynamical models. Based on a set of 400 clone orbits we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive numerical investigation of the stability region of the Earth’s Lagrangian points. We present a detailed parameter study in the regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of stability.

Recently the first Earth Trojan has been observed (Connors et al., Nature 475) and found to be on an interesting orbit close to the Lagrange point L4. In the present study we therefore perform a detailed investigation on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our results are derived using different approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the phase space geometry, including the indirect influence of the additional planets of our Solar system. b) We use precise numerical methods to integrate the orbit forward and backward in time in different dynamical models. Based on a set of 400 clone orbits we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive numerical investigation of the stability region of the Earth’s Lagrangian points. We present a detailed parameter study in the regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of stability.

The electron-cyclotron maser instability is widespread in the Universe, producing, e.g., radio emission of the magnetized planets and cool substellar objects. Diagnosing the parameters of astrophysical radio sources requires comprehensive nonlinear simulations of the radiation process. We simulate the electron-cyclotron maser instability in a very low-beta plasma. The model used takes into account the radiation escape from the source region and the particle flow through this region. We developed a kinetic code to simulate the time evolution of an electron distribution in a radio emission source. The model includes the terms describing the particle injection to and escape from the emission source region. The spatial escape of the emission from the source is taken into account by using a finite amplification time. The unstable electron distribution of the horseshoe type is considered. A number of simulations were performed for different parameter sets typical of the magnetospheres of planets and ultracool dwarfs. The generated emission (corresponding to the fundamental extraordinary mode) has a frequency close to the electron cyclotron frequency and propagates across the magnetic field. Shortly after the onset of a simulation, the electron distribution reaches a quasi-stationary state. If the emission source region is relatively small, the resulting electron distribution is similar to that of the injected electrons; the emission intensity is low. In larger sources, the electron distribution may become nearly flat due to the wave-particle interaction, while the conversion efficiency of the particle energy flux into waves reaches 10-20%. We found good agreement of our model with the in situ observations in the source regions of auroral radio emissions of the Earth and Saturn. The expected characteristics of the electron distributions in the magnetospheres of ultracool dwarfs were obtained.

Most known extrasolar planets (exoplanets) have been discovered using the radial velocity$^{\bf 1,2}$ or transit$^{\bf 3}$ methods. Both are biased towards planets that are relatively close to their parent stars, and studies find that around 17–30% (refs 4, 5) of solar-like stars host a planet. Gravitational microlensing$^{\bf 6\rm{\bf -}\bf 9}$, on the other hand, probes planets that are further away from their stars. Recently, a population of planets that are unbound or very far from their stars was discovered by microlensing$^{\bf 10}$. These planets are at least as numerous as the stars in the Milky Way$^{\bf 10}$. Here we report a statistical analysis of microlensing data (gathered in 2002–07) that reveals the fraction of bound planets 0.5–10 AU (Sun–Earth distance) from their stars. We find that 17$_{\bf -9}^{\bf +6}$% of stars host Jupiter-mass planets (0.3–10 $\MJ$, where $\MJ {\bf = 318}$ $\Mearth$ and $\Mearth$ is Earth’s mass). Cool Neptunes (10–30 $\Mearth$) and super-Earths (5–10 $\Mearth$) are even more common: their respective abundances per star are 52$_{\bf -29}^{\bf +22}$% and 62$_{\bf -37}^{\bf +35}$%. We conclude that stars are orbited by planets as a rule, rather than the exception.

We present the first results of the Hubble Wide Field Camera 3 Test of Surfaces in the Outer Solar System (H/WTSOSS). The purpose of this survey was to measure the surface properties of a large number of Kuiper belt objects and attempt to infer compositional and dynamical correlations. We find that the Centaurs and the low-perihelion scattered disk and resonant objects exhibit virtually identical bifurcated optical colour distributions and make up two well defined groups of object. Both groups have highly correlated optical and NIR colours which are well described by a pair of two component mixture models that have different red components, but share a common neutral component. The small, $H_{606}\gtrsim5.6$ high-perihelion excited objects are entirely consistent with being drawn from the two branches of the mixing model suggesting that the colour bifurcation of the Centaurs is apparent in all small excited objects. On the other hand, objects larger than $H_{606}\sim5.6$ are not consistent with the mixing model, suggesting some evolutionary process avoided by the smaller objects. The existence of a bifurcation amongst all excited populations argues that the two separate classes of object existed in the primordial disk before the excited Kuiper belt was populated. The cold classical objects exhibit a different type of surface which has colours that are consistent with being drawn from the red branch of the mixing model, but with much higher albedos.

Precise astronomical spectroscopic analyses routinely assume that individual pixels in charge-coupled devices (CCDs) have uniform sensitivity to photons. Intra-pixel sensitivity (IPS) variations may already cause small systematic errors in, for example, studies of extra-solar planets via stellar radial velocities and cosmological variability in fundamental constants via quasar spectroscopy, but future experiments requiring velocity precisions approaching ~1 cm/s will be more strongly affected. Laser frequency combs have been shown to provide highly precise wavelength calibration for astronomical spectrographs, but here we show that they can also be used to measure IPS variations in astronomical CCDs in situ. We successfully tested a laser frequency comb system on the Ultra-High Resolution Facility spectrograph at the Anglo-Australian Telescope. By modelling the 2-dimensional comb signal recorded in a single CCD exposure, we find that the average IPS deviates by <8 per cent if it is assumed to vary symmetrically about the pixel centre. We also demonstrate that series of comb exposures with absolutely known offsets between them can yield tighter constraints on symmetric IPS variations from ~100 pixels. We discuss measurement of asymmetric IPS variations and absolute wavelength calibration of astronomical spectrographs and CCDs using frequency combs.

(Abridged) Planetesimals embedded in a protoplanetary disc are stirred by gravitational torques exerted by density fluctuations in the surrounding turbulence. In particular, planetesimals in a disc supporting fully developed magneto-rotational turbulence are readily excited to velocity dispersions above the threshold for catastrophic disruption, halting planet formation. We aim to examine the stirring of planetesimals lying instead in a magnetically-decoupled midplane dead zone, stirred only by spiral density waves propagating out of the disc’s magnetically-coupled turbulent surface layers. We extend previous studies to include a wider range of disc models, and explore the effects of varying the disc column density and external magnetic field strength. [...] The strength of the stirring is found to be independent of the gas surface density, which is contrary to the increase with disc mass expected from a simple linear wave picture. The discrepancy arises from the shearing out of density waves as they propagate into the dead zone, resulting in density structures near the midplane that exert weaker stochastic torques on average. We provide a simple analytic fit to our numerically obtained torque amplitudes that accounts for this effect. The stirring on the other hand depends sensitively on the net vertical magnetic flux, up to a saturation level above which magnetic forces dominate in the turbulent layers. For the majority of our models, the equilibrium planetesimal velocity dispersions lie between the thresholds for disrupting strong and weak aggregates, suggesting that collision outcomes will depend on material properties. However, discs with relatively weak magnetic fields yield reduced stirring, and their dead zones provide safe-havens even for the weakest planetesimals against collisional destruction.

A growing sample of white dwarfs (WDs) with metal-enriched atmospheres are accompanied by excess infrared emission, indicating that they are encircled by a compact dusty disk of solid debris. Such `WD debris disks’ are thought to originate from the tidal disruption of asteroids or other minor bodies, but the precise mechanism(s) responsible for transporting matter to the WD surface remains unclear, especially in those systems with the highest inferred metal accretion rates dM_Z/dt ~ 1e8-1e10 g/s. Here we present global time-dependent calculations of the coupled evolution of the gaseous and solid components of WD debris disks. Solids transported inwards (initially due to PR drag) sublimate at tens of WD radii, producing a source of gas that accretes onto the WD surface and viscously spreads outwards in radius, where it overlaps with the solid disk. If the aerodynamic coupling between the solids and gaseous disks is sufficiently strong (and/or the gas viscosity sufficiently weak), then gas builds up near the sublimation radius faster than it can viscously spread away. Since the rate of drag-induced solid accretion increases with gas density, this results in a runaway accretion process, during which the WD accretion rate reaches values orders of magnitude higher than can be achieved by PR drag alone. We explore the evolution of WD debris disks across a wide range of physical conditions and calculate the predicted distribution of observed accretion rates dM_Z/dt, finding reasonable agreement with the current sample. Although the conditions necessary for runaway accretion are at best marginally satisfied given the minimal level of aerodynamic drag between circular gaseous and solid disks, the presence of other stronger forms of solid-gas coupling—such as would result if the gaseous disk is only mildly eccentric—substantially increase the likelihood of runaway accretion.

We compare limb darkening laws derived from 3D hydrodynamical model atmospheres and 1D hydrostatic MARCS models for the host stars of the two transiting exoplanet systems HD 209458 and HD 189733. The surface brightness distribution of the stellar disks is calculated using 3D LTE spectrum formation and opacity sampling. We test our predictions using least-squares fits of model light curves to primary eclipses that were observed with the Hubble Space Telescope (HST). The limb darkening law derived from the 3D model of HD 209458 between 2900 A and 5700 A produces significantly better fits to the HST data, removing systematic residuals that were previously observed for model light curves based on 1D predictions. This difference arises mainly from the shallower mean temperature structure of the 3D model, which is a consequence of the explicit simulation of surface granulation. In the case of HD 189733, the model atmospheres produce practically equivalent limb darkening curves between 2900 A and 5700 A, partly due to obstruction by spectral lines, and the data are not sufficient to distinguish between the light curves. We also analyze HST observations between 5350 A and 10500 A for this star; the 3D model leads to a better fit compared to 1D predictions. The significant improvement of fit quality for the HD 209458 system demonstrates the higher degree of realism of 3D models and the importance of surface granulation for the formation of the atmospheric radiation field of late-type stars. This result agrees well with recent investigations of limb darkening in the solar continuum and other observational tests. The case of HD 189733 is no contradiction as the model light curves are less sensitive to the temperature stratification of the atmosphere and the observed data in the 2900 A – 5700 A region are not sufficient to distinguish more clearly between 3D and 1D limb darkening.

We re-analyze 4 years of HARPS spectra of the nearby M1.5 dwarf GJ 667C available through the ESO public archive. The new radial velocity (RV) measurements were obtained using a new data analysis technique that derives the Doppler measurement and other instrumental effects using a least-squares approach. Combining these new 143 measurements with 41 additional RVs from the Magellan/PFS and Keck/HIRES spectrometers, reveals 3 additional signals beyond the previously reported 7.2-day candidate, with periods of 28 days, 75 days, and a secular trend consistent with the presence of a gas giant (Period sim 10 years). The 28-day signal implies a planet candidate with a minimum mass of 4.5 Mearth orbiting well within the canonical definition of the star’s liquid water habitable zone, this is, the region around the star at which an Earth-like planet could sustain liquid water on its surface. Still, the ultimate water supporting capability of this candidate depends on properties that are unknown such as its albedo, atmospheric composition and interior dynamics. The 75-day signal is less certain, being significantly affected by aliasing interactions among a potential 91-day signal, and the likely rotation period of the star at 105 days detected in two activity indices. GJ 667C is the common proper motion companion to the GJ 667AB binary, which is metal poor compared to the Sun. The presence of a super-Earth in the habitable zone of a metal poor M dwarf in a triple star system, supports the evidence that such worlds should be ubiquitous in the Galaxy.

We present high-precision radial velocity (RV) observations of four solar-type (F7-G5) stars – HD 79498, HD 155358, HD 197037, and HD 220773 – taken as part of the McDonald Observatory Planet Search Program. For each of these stars, we see evidence of Keplerian motion caused by the presence of one or more gas giant planets in long-period orbits. We derive orbital parameters for each system, and note the properties (composition, activity, etc.) of the host stars. While we have previously announced the two-gas-giant HD 155358 system, we now report a shorter period for planet c. This new period is consistent with the planets being trapped in mutual 2:1 mean-motion resonance. We therefore perform an in-depth stability analysis, placing additional constraints on the orbital parameters of the planets. These results demonstrate the excellent long-term RV stability of the spectrometers on both the Harlan J. Smith 2.7 m telescope and the Hobby-Eberly telescope.

We quantified ten parent volatiles in comet C/2009 P1 (Garradd) before perihelion, through high-dispersion infrared spectra acquired with CRIRES at ESO’s VLT on UT 2011 August 07 (Rh = 2.4 AU) and September 17-21 (Rh = 2.0 AU). On August 07, water was searched but not detected at an upper limit (3{\sigma}) of 2.1 \times 10^28 s-1, while ethane was detected with a production rate of 6.1 \times 10^26 s-1 (apparent mixing ratio > 2.90%). On September 17-21, the mean production rate for water was 8.4 \times 10^28 s-1, and abundance ratios (relative to water) of detected trace species were: CO (12.51%), CH3OH (3.90%), CH4 (1.24%), C2H6 (1.01%) and HCN (0.36%). Upper limits (3{\sigma}) to abundances for four minor species were: NH3 (1.55%), C2H2 (0.13%), HDO (0.89%) and OCS (0.20%). Given the relatively large heliocentric distance, we explored the effect of water not being fully sublimated within our FOV and identified the ‘missing’ water fraction needed to reconcile the retrieved abundance ratios with the mean values found for “organics-normal”. The individual spatial profiles of parent volatiles and the continuum displayed rather asymmetric outgassing. Indications of H2O and CO gas being released in different directions suggest different active vents and/or the possible existence of polar and apolar ice aggregates in the nucleus. The high fractional abundance of CO identifies comet C/2009 P1 as a CO-rich comet.

The late-type dwarf GJ 436 is known to host a transiting Neptune-mass planet in a 2.6-day orbit. We present results of our interferometric measurements to directly determine the stellar diameter ($R_{\star} = 0.455 \pm 0.018 R_{\odot}$) and effective temperature ($T_{\rm EFF} = 3416 \pm 54$ K). We combine our stellar parameters with literature time-series data, which allows us to calculate physical and orbital system parameters, including GJ 436’s stellar mass ($M_{\star} = 0.472^{+ 0.0636}_{- 0.0566} M_{\odot}$), planetary radius ($R_{p} = 0.370^{+ 0.0149}_{- 0.0145} R_{Jupiter}$), planetary mass ($M_{p} = 0.075^{+ 0.0076}_{- 0.0072} M_{Jupiter}$), implying a mean planetary density of $\rho_{p} = 1.48^{+ 0.116}_{- 0.103} \rho_{Jupiter}$. These values are generally in good agreement with previous literature estimates based on assumed stellar mass and photometric light curve fitting. Finally, we examine the expected phase curves of the hot Neptune GJ 436b, based on various assumptions concerning the efficiency of energy redistribution in the planetary atmosphere, and find that it could be constrained with {\it Spitzer} monitoring observations.

We study the radius evolution of close-in extra-solar jupiters under Ohmic heating, a mechanism that was recently proposed to explain the large observed sizes of many of these planets. Planets are born with high entropy and they subsequently cool and contract. We focus on two cases: first, that ohmic heating commences when the planet is hot (high entropy); and second, that it commences after the planet has cooled. In the former case, we use analytical scalings and numerical experiments to confirm that Ohmic heating is capable of suspending the cooling as long as a few percent of the stellar irradiation is converted into Ohmic heating, and the planet has a surface wind that extends to pressures of ~10 bar or deeper. For these parameters, the radii at which cooling is stalled are consistent with (or larger than) the observed radii of most planets. The only two exceptions are WASP-17b and HAT-P-32b. In contrast to the high entropy case, we show that Ohmic heating cannot significantly re-inflate planets after they have already cooled. This leads us to suggest that the diversity of radii observed in hot jupiters may be partially explained by the different epochs at which they are migrated to their current locations.

Exoplanetary systems are found not only among single stars, but also binaries of widely varying parameters. Binaries with separations of 100–1000 au are prevalent in the Solar neighborhood; at these separations planet formation around a binary member may largely proceed as if around a single star. During the early dynamical evolution of a planetary system, planet–planet scattering can eject planets from a star’s grasp. In a binary, the motion of a planet ejected from one star has effectively entered a restricted three-body system consisting of itself and the two stars, and the equations of motion of the three body problem will apply as long as the ejected planet remains far from the remaining planets. Depending on its energy, escape from the binary as a whole may be impossible or delayed until the three-body approximation breaks down, and further close interactions with its planetary siblings boost its energy when it passes close to its parent star. Until then this planet may be able to transition from the space around one star to the other, and chaotically `bounce’ back and forth. In this paper we directly simulate scattering planetary systems that are around one member of a circular binary, and quantify the frequency of bouncing in scattered planets. We find that a great majority (70 to 85 per cent) of ejected planets will pass at least once through the space of it’s host’s binary companion, and depending on the binary parameters about 45 to 75 per cent will begin bouncing. The time spent bouncing is roughly log-normally distributed with a peak at about $10^4$ years, with only a small percentage bouncing for more than a Myr. This process may perturb and possibly incite instability among existing planets around the companion star. In rare cases, the presence of multiple planets orbiting both stars may cause post-bouncing capture or planetary swapping.

We present Hubble Space Telescope near-infrared transmission spectroscopy of the transiting exoplanet HD 189733b, using Wide Field Camera 3. This consists of time-series spectra of two transits, used to measure the wavelength dependence of the planetary radius. These observations aim to test whether the Rayleigh scattering haze detected at optical wavelengths extends into the near-infrared, or if it becomes transparent leaving molecular features to dominate the transmission spectrum. Due to saturation and non-linearity affecting the brightest (central) pixels of the spectrum, light curves were extracted from the blue and red ends of the spectra only, corresponding to wavelength ranges of 1.099-1.168 um and 1.521-1.693 um, respectively, for the first visit, and 1.082-1.128 um and 1.514-1.671 um for the second. The light curves were fitted using a Gaussian process model to account for instrumental systematics whilst simultaneously fitting for the transit parameters. This gives values of the planet-to-star radius ratio for the blue and red light curves of 0.15650\pm0.00048 and 0.15634\pm0.00032, respectively, for visit one and 0.15716\pm0.00078 and 0.15630\pm0.00037 for visit 2 (using a quadratic limb darkening law). The planet-to-star radius ratios measured in both visits are consistent, and we see no evidence for the drop in absorption expected if the haze that is observed in the optical becomes transparent in the infrared. This tentatively suggests that the haze dominates the transmission spectrum of HD 189733b into near-infrared wavelengths, although more robust observations are required to provide conclusive evidence.

We present new photometric data of the transiting planet HAT-P-12b observed in 2011. Our three transit curves are modelled using the JKTEBOP code and adopting the quadratic limb-darkening law. Including our measurements, 18 transit times spanning about 4.2 yr were used to determine the improved ephemeris with a transit epoch of 2,454,187.85560$\pm$0.00011 BJD and an orbital period of 3.21305961$\pm$0.00000035 d. The physical properties of the star-planet system are computed using empirical calibrations from eclipsing binary stars and stellar evolutionary models, combined with both our transit parameters and previously-known spectroscopic results. We found that the absolute dimensions of the host star are $M_{\rm A}$=0.73$\pm$0.02 M$_\odot$, $R_{\rm A}$=0.70$\pm$0.01 R$_\odot$, log $g_{\rm A}$=4.61$\pm$0.02, $\rho_{\rm A}$=2.10$\pm$0.09 $\rho_{\odot}$, and $L_{\rm A}$=0.21$\pm$0.01 L$_\odot$. The planetary companion has $M_{\rm b}$=0.21$\pm$0.01 M$_{\rm Jup}$, $R_{\rm b}$=0.94$\pm$0.01 R$_{\rm Jup}$, log $g_{\rm b}$=2.77$\pm$0.02, $\rho_{\rm b}$=0.24$\pm$0.01 $\rho_{\rm Jup}$, and $T_{\rm eq}$=960$\pm$14 K. Our results agree well with standard models of irradiated gas giants with a core mass of 11.3 M$_{\oplus}$.

The abundance of CO2 in comets has been derived using CO Cameron band (a3pi –> X1Sigma+) emission assuming that photodissociative excitation of CO2 is the main production process of CO(a3pi). On comet 1P/Halley the Cameron (1-0) band has been observed by International Ultraviolet Explorer (IUE) on several days in March 1986. A coupled chemistry-emission model is developed for comet 1P/Halley to assess the importance of various production and loss mechanisms of CO(a3pi) and to calculate the intensity of Cameron band emission on different days of IUE observation. Two different solar EUV flux models, EUVAC of Richards et al. (1994) and SOLAR2000 of Tobiska (2004), and different relative abundances of CO and CO2, are used to evaluate the role of photon and photoelectron in producing CO molecule in a3pi state in the cometary coma. It is found that in comet 1P/Halley 60–70% of the total intensity of the Cameron band emission is contributed by electron impact excitation of CO and CO2, while the contribution from photodissociative excitation of CO2 is small (20–30%). Thus, in the comets where CO and CO2 relative abundances are comparable, the Cameron band emission is largely governed by electron impact excitation of CO, and not by the photodissociative excitation of CO2 as assumed earlier. Model calculated Cameron band 1-0 emission intensity (40 R) is consistent with the observed IUE slit-averaged brightness (37 +/- 6 R) using EUVAC model solar flux on 13 March 1986, and also on other days of observations. Since electron impact excitation is the major production mechanism, the Cameron emission can be used to derive photoelectron density in the inner coma rather than the CO2 abundance.

We present a new model atmosphere analysis of the most metal contaminated white dwarf known, the DBZ SDSS J073842.56+183509.06. Using new high resolution spectroscopic observations taken with Keck and Magellan, we determine precise atmospheric parameters and measure abundances of 14 elements heavier than helium. We also report new Spitzer mid-infrared photometric data that are used to better constrain the properties of the debris disk orbiting this star. Our detailed analysis, which combines data taken from 7 different observational facilities (GALEX, Gemini, Keck, Magellan, MMT, SDSS and Spitzer) clearly demonstrate that J0738+1835 is accreting large amounts of rocky terrestrial-like material that has been tidally disrupted into a debris disk. We estimate that the body responsible for the photospheric metal contamination was at least as large Ceres, but was much drier, with less than 1% of the mass contained in the form of water ice, indicating that it formed interior to the snow line around its parent star. We also find a correlation between the abundances (relative to Mg and bulk Earth) and the condensation temperature; refractory species are clearly depleted while the more volatile elements are possibly enhanced. This could be the signature of a body that formed in a lower temperature environment than where Earth formed. Alternatively, we could be witnessing the remains of a differentiated body that lost a large part of its outer layers.

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the stars are among the brightest in the Kepler field. We present an analysis of three months of nearly uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 46 and 41 oscillation frequencies for the two components respectively, including a clear detection of octupole (l=3) modes in both stars. We derive the properties of each star independently using the Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by fitting each component individually yield a common age (t=6.8+/-0.4 Gyr) and initial composition (Z_i=0.024+/-0.002, Y_i=0.25+/-0.01) within the uncertainties, as expected for the components of a binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data sets that will ultimately become available will allow future studies of differential rotation, convection zone depths, and long-term changes due to stellar activity cycles.

The evolved solar-type stars 16 Cyg A & B have long been studied as solar analogs, yielding a glimpse into the future of our own Sun. The orbital period of the binary system is too long to provide meaningful dynamical constraints on the stellar properties, but asteroseismology can help because the stars are among the brightest in the Kepler field. We present an analysis of three months of nearly uninterrupted photometry of 16 Cyg A & B from the Kepler space telescope. We extract a total of 46 and 41 oscillation frequencies for the two components respectively, including a clear detection of octupole (l=3) modes in both stars. We derive the properties of each star independently using the Asteroseismic Modeling Portal, fitting the individual oscillation frequencies and other observational constraints simultaneously. We evaluate the systematic uncertainties from an ensemble of results generated by a variety of stellar evolution codes and fitting methods. The optimal models derived by fitting each component individually yield a common age (t=6.8+/-0.4 Gyr) and initial composition (Z_i=0.024+/-0.002, Y_i=0.25+/-0.01) within the uncertainties, as expected for the components of a binary system, bolstering our confidence in the reliability of asteroseismic techniques. The longer data sets that will ultimately become available will allow future studies of differential rotation, convection zone depths, and long-term changes due to stellar activity cycles.

Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material “rain out”, which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

We present a catalogue of rotation and chromospheric activity in a sample of 334 M dwarfs of spectral types M0–M4.5 populating the parameter space around the boundary to full convection. We obtained high-resolution optical spectra for 206 targets and determined projected rotational velocity, vsini, and Halpha emission. The data are combined with measurements of vsini in field stars of the same spectral type from the literature. Our sample adds 157 new rotation measurements to the existing literature and almost doubles the sample of available vsini. The final sample provides a statistically meaningful picture of rotation and activity at the transition to full convection in the solar neighborhood. We confirm the steep rise in the fraction of active stars at the transition to full convection known from earlier work. In addition, we see a clear rise in rotational velocity in the same stars. In very few stars, no chromospheric activity but a detection of rotational broadening was reported. We argue that all of them are probably spurious detections; we conclude that in our sample all significantly rotating stars are active, and all active stars are significantly rotating. The rotation-activity relation is valid in partially and in fully convective stars. Thus, we do not observe any evidence for a transition from a rotationally dominated dynamo in partially convective stars to a rotation-independent turbulent dynamo in fully convective stars; turbulent dynamos in fully convective stars of spectral types around M4 are still driven by rotation. Finally, we compare projected rotational velocities of 33 stars to rotational periods derived from photometry in the literature and determine inclinations for a few of them.

In early 2011, the discovery of two planets moving on surprisingly extreme orbits around the eclipsing polar cataclysmic variable system HU Aquraii was announced based on variations in the timing of mutual eclipses between the two central stars. We perform a detailed dynamical analysis of the stability of the exoplanet system as proposed in that work, revealing that it is simply dynamically unfeasible. We then apply the latest rigorous methods used by the Anglo-Australian Planet Search to analyse radial velocity data to re-examine the data used to make the initial claim. Using that data, we arrive at a significantly different orbital solution for the proposed planets, which we then show through dynamical analysis to be equally unfeasible. Finally, we discuss the need for caution in linking eclipse-timing data for cataclysmic variables to the presence of planets, and suggest a more likely explanation for the observed signal.

In early 2011, the discovery of two planets moving on surprisingly extreme orbits around the eclipsing polar cataclysmic variable system HU Aquraii was announced based on variations in the timing of mutual eclipses between the two central stars. We perform a detailed dynamical analysis of the stability of the exoplanet system as proposed in that work, revealing that it is simply dynamically unfeasible. We then apply the latest rigorous methods used by the Anglo-Australian Planet Search to analyse radial velocity data to re-examine the data used to make the initial claim. Using that data, we arrive at a significantly different orbital solution for the proposed planets, which we then show through dynamical analysis to be equally unfeasible. Finally, we discuss the need for caution in linking eclipse-timing data for cataclysmic variables to the presence of planets, and suggest a more likely explanation for the observed signal.

We present a statistical analysis that demonstrates that the overwhelming majority of Kepler candidate multiple transiting systems (multis) indeed represent true, physically-associated transiting planets. Binary stars provide the primary source of false positives among Kepler planet candidates, implying that false positives should be nearly randomly-distributed among Kepler targets. In contrast, true transiting planets would appear clustered around a smaller number of Kepler targets if detectable planets tend to come in systems and/or if the orbital planes of planets encircling the same star are correlated. There are more than one hundred times as many Kepler planet candidates in multi-candidate systems as would be predicted from a random distribution of candidates, implying that the vast majority are true planets. Most of these multis are multiple planet systems orbiting the Kepler target star, but there are likely cases where (a) the planetary system orbits a fainter star, and the planets are thus significantly larger than has been estimated, or (b) the planets orbit different stars within a binary/multiple star system. We use the low overall false positive rate among Kepler multis, together with analysis of Kepler spacecraft and ground-based data, to validate the closely-packed Kepler-33 planetary system, which orbits a star that has evolved somewhat off of the main sequence. Kepler-33 hosts five transiting planets with periods ranging from 5.67 to 41 days.

Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets’ masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems’ architectures, even in cases for which high-precision Doppler follow-up is impractical.

We present a method to confirm the planetary nature of objects in systems with multiple transiting exoplanet candidates. This method involves a Fourier-Domain analysis of the deviations in the transit times from a constant period that result from dynamical interactions within the system. The combination of observed anti-correlations in the transit times and mass constraints from dynamical stability allow us to claim the discovery of four planetary systems Kepler-25, Kepler-26, Kepler-27, and Kepler-28, containing eight planets and one additional planet candidate.

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.