Archive for category Earth and Planetary

Today's Postings

On the unusual gas composition in the Beta Pictoris debris disk

The metallic gas associated with the Beta Pic debris disk is believed to not be primordial, but arise during the destruction of dust grains. Recent observations have shown that carbon and oxygen in this gas are exceptionally overabundant compared to other elements, by some 400 times. We study the origin of this enrichment under two opposing hypothesis, preferential production, where the gas is produced with the observed unusual abundance, and preferential depletion, where the gas evolves to the observed state from an original solar abundance under a number of dynamical processes. We include in our study the following processes: radiative blow-out of metallic elements, dynamical coupling between different species, and viscous accretion onto the star. We find that, if gas viscosity is sufficiently low (the conventional alpha parameter 1e-1, as expected for this largely ionized disk), gas is continuously produced and viscously accreted toward the star. This removal process does not discriminate between elements so the observed overabundance of C and O has to be explained by a preferential production that strongly favors C and O to other metallic elements. One such candidate is photo-desorption off the grains. We compare our calculation against all observed elements (~10) in the gas disk and find a mild preference for the second scenario, based on the abundance of Si alone. If true, Beta Pic should still be accreting at an observable rate, well after its primordial disk has disappeared.

Comet 22P/Kopff: Dust environment and grain ejection anisotropy from visible and infrared observations

We present optical observations and Monte Carlo models of the dust coma, tail, and trail structures of comet 22P/Kopff during the 2002 and 2009 apparitions. Dust loss rates, ejection velocities, and power-law size distribution functions are derived as functions of the heliocentric distance using pre- and post-perihelion imaging observations during both apparitions. The 2009 post-perihelion images can be accurately fitted by an isotropic ejection model. On the other hand, strong dust ejection anisotropies are required to fit the near-coma regions at large heliocentric distances (both inbound at $r_h$=2.5 AU and outbound at $r_h$=2.6 AU) for the 2002 apparition. These asymmetries are compatible with a scenario where dust ejection is mostly seasonally-driven, coming mainly from regions near subsolar latitudes at far heliocentric distances inbound and outbound. At intermediate to near-perihelion heliocentric distances, the outgassing would affect much more extended latitude regions, the emission becoming almost isotropic near perihelion. We derived a maximum dust production rate of 260 kg s$^{-1}$ at perihelion, and an averaged production rate over one orbit of 40 kg s$^{-1}$. An enhanced emission rate, accompanied also by a large ejection velocity, is predicted at $r_h>$2.5 pre-perihelion. The model has also been extended to the thermal infrared in order to be applied to available trail observations with IRAS and ISO spacecrafts of this comet. The modeled trail intensities are in good agreement with those observations, which is remarkable taking into account that those data are sensitive to dust ejection patterns corresponding to several orbits before the 2002 and 2009 apparitions.

Tidal synchronization of close-in satellites and exoplanets. A rheophysical approach [Replacement]

New theory of the dynamical tides of celestial bodies founded on a Newtonian creep instead of the classical delaying approach of the standard viscoelastic theories. The results of the theory derive mainly from the solution of a non-homogeneous ordinary differential equation. Lags appear in the solution, but as quantities determined from the solution of the equation and are not arbitrary external quantities plugged on an elastic model. The resulting lag of each tide component is an increasing function of its frequency (as in Darwin’s theory), and lags are not small quantities. The amplitudes of the tide components depend on the viscosity of the body and on their frequencies; they are not constants. The resulting stationary rotations (often called pseudo-synchronous) have an excess velocity roughly proportional to 6ne^2/(X^2+1/X^2) (X is the mean-motion in units of one relaxation factor inversely proportional to the viscosity) instead of the exact 6ne^2 of standard theories. The dissipation in the pseudo-synchronous solution is inversely proportional to (X+1/X); thus, in the inviscid limit it is roughly proportional to the frequency (as in standard theories), but that behavior is inverted when the viscosity is high and the relaxation factor much smaller than the tide frequency. For free rotating bodies, the dissipation is given by the same law, but now X is the frequency of the semidiurnal tide in units of the relaxation factor. This approach fails, however, to reproduce the actual tidal lags on Earth and on natural satellites. To reconcile theory and observations, in this case, we had to assume the coexistence of an elastic tide superposed to the creeping tide. The theory is applied to several Solar System and extrasolar bodies and values of the relaxation factor \gamma\ are derived for these bodies on the basis of currently available data.

Tidal synchronization of close-in satellites and exoplanets. A rheophysical approach

New theory of the dynamical tides of celestial bodies founded on a Newtonian creep instead of the classical delaying approach of the standard viscoelastic theories. The results of the theory derive mainly from the solution of a non-homogeneous ordinary differential equation. Lags appear in the solution, but as quantities determined from the solution of the equation and are not arbitrary external quantities plugged on an elastic model. The resulting lag of each tide component is an increasing function of its frequency (as in Darwin’s theory), and are not small quantities. The amplitudes of the tide components depend on the viscosity of the body and on their frequencies; they are not constants. The resulting stationary rotations (often called pseudo-synchronous) have an excess velocity roughly proportional to 6ne^2/(\chi^2+1/\chi^2) (\chi\ is the mean-motion in units of one relaxation factor inversely proportional to the viscosity) instead of the exact 6ne^2 of standard theories. The dissipation in the pseudo-synchronous solution is inversely proportional to (\chi+1/\chi); thus, in the inviscid limit it is roughly proportional to the frequency (as in standard theories), but that behavior is inverted when the viscosity is high and the relaxation factor much smaller than the tide frequency. For free rotating bodies, the dissipation is given by the same law, but now \chi\ is the frequency of the semidiurnal tide in units of the relaxation factor. This approach fails, however, to reproduce the actual tidal lags on Earth and on natural satellites. To reconcile theory and observations, in this case, we had to assume the coexistence of a small elastic tide superposed to the creeping tide. The theory is applied to several Solar System and extrasolar bodies and values of the relaxation factor \gamma\ are derived for these bodies on the basis of currently available data.

The Transiting Circumbinary Planets Kepler-34 and Kepler-35

Most Sun-like stars in the Galaxy reside in gravitationally-bound pairs of stars called “binary stars”. While long anticipated, the existence of a “circumbinary planet” orbiting such a pair of normal stars was not definitively established until the discovery of Kepler-16. Incontrovertible evidence was provided by the miniature eclipses (“transits”) of the stars by the planet. However, questions remain about the prevalence of circumbinary planets and their range of orbital and physical properties. Here we present two additional transiting circumbinary planets, Kepler-34 and Kepler-35. Each is a low-density gas giant planet on an orbit closely aligned with that of its parent stars. Kepler-34 orbits two Sun-like stars every 289 days, while Kepler-35 orbits a pair of smaller stars (89% and 81% of the Sun’s mass) every 131 days. Due to the orbital motion of the stars, the planets experience large multi-periodic variations in incident stellar radiation. The observed rate of circumbinary planets implies > ~1% of close binary stars have giant planets in nearly coplanar orbits, yielding a Galactic population of at least several million.

Mutual Events in the Cold Classical Transneptunian Binary System Sila and Nunam

Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season of mutual events, in which the two near-equal brightness bodies alternate in passing in front of one another as seen from Earth, is in progress right now, and on 2011 Feb. 1 UT, one such event was observed from two different telescopes. The mutual event season offers a rich opportunity to learn much more about this barely-resolvable binary system, potentially including component sizes, colors, shapes, and albedo patterns. The low eccentricity of the orbit and a photometric lightcurve that appears to coincide with the orbital period are consistent with a system that is tidally locked and synchronized, like the Pluto-Charon system. The orbital period and semimajor axis imply a system mass of (10.84 \pm 0.22) \times 10^18 kg, which can be combined with a size estimate based on Spitzer and Herschel thermal infrared observations to infer an average bulk density of 0.72 +0.37 -0.23 g cm^-3, comparable to the very low bulk densities estimated for small transneptunian binaries of other dynamical classes.

Stability Regions Around the Components of the Triple System 2001 SN263 [Replacement]

The NEAs (Near-Earth Asteroids) are good targets for spatial missions, since they periodically approach the orbit of the Earth. Recently, the NEA (153591) 2001 SN263 was chosen as the target of the ASTER MISSION- First Brazilian Deep Space Mission, planned to be launched in 2015. In February 2008, the radio astronomers from Arecibo-Puerto Rico concluded that (153591) 2001 SN263 is actually a triple system (Nolan et al., 2008). The announcement of the ASTER MISSION has motivated the development of the present work, whose goal is to characterize regions of stability and instability of the triple system (153591) 2001 SN263. The method adopted consisted in dividing the region around the system into four distinct regions. We have performed numerical integrations of systems composed by seven bodies: Sun, Earth, Mars, Jupiter and the three components of the system, and by thousands of particles randomly distributed within the demarcated regions, for the planar and inclined prograde cases. The results are diagrams of semi-major axis versus eccentricity, where it is shown the percentage of particles that survive for each set of initial conditions. The regions where 100% of the particles survive is defined as stable regions. We found that the stable regions are in the neighborhood of Alpha and Beta, and in the external region. It was identified resonant motion of the particles with Beta and Gamma in the internal regions, which lead to instability. For particles with I>45{\deg} in the internal region, where I is the inclination with respect to Alpha’s equator, there is no stable region, except for the particles placed really close to Alpha. The stability in the external region is not affected by the variation of inclination. We also present a discussion on the long-term stability in the internal region, for the planar and circular cases, with comparisons with the short-term stability.

Stability Regions Around the Components of the Triple System 2001 SN263

The NEAs (Near-Earth Asteroids) are good targets for spatial missions, since they periodically approach the orbit of the Earth. Recently, the NEA (153591) 2001 SN263 was chosen as the target of the ASTER MISSION- First Brazilian Deep Space Mission, planned to be launched in 2015. In February 2008, the radio astronomers from Arecibo-Puerto Rico concluded that (153591) 2001 SN263 is actually a triple system (Nolan et al., 2008). The announcement of the ASTER MISSION has motivated the development of the present work, whose goal is to characterize regions of stability and instability of the triple system (153591) 2001 SN263. The method adopted consisted in dividing the region around the system into four distinct regions. We have performed numerical integrations of systems composed by seven bodies: Sun, Earth, Mars, Jupiter and the three components of the system, and by thousands of particles randomly distributed within the demarcated regions, for the planar and inclined prograde cases. The results are diagrams of semi-major axis versus eccentricity, where it is shown the percentage of particles that survive for each set of initial conditions. The regions where 100% of the particles survive is defined as stable regions. We found that the stable regions are in the neighborhood of Alpha and Beta, and in the external region. It was identified resonant motion of the particles with Beta and Gamma in the internal regions, which lead to instability. For particles with I>45{\deg} in the internal region, where I is the inclination with respect to Alpha’s equator, there is no stable region, except for the particles placed really close to Alpha. The stability in the external region is not affected by the variation of inclination. We also present a discussion on the long-term stability in the internal region, for the planar and circular cases, with comparisons with the short-term stability.

Calibration of Equilibrium Tide Theory for Extrasolar Planet Systems II

We present a new empirical calibration of equilibrium tidal theory for extrasolar planet systems, extending a prior study by incorporating detailed physical models for the internal structure of planets and host stars. The resulting strength of the stellar tide produces a coupling that is strong enough to reorient the spins of some host stars without causing catastrophic orbital evolution, thereby potentially explaining the observed trend in alignment between stellar spin and planetary orbital angular momentum. By isolating the sample whose spins should not have been altered in this model, we also show evidence for two different processes that contribute to the population of planets with short orbital periods. We apply our results to estimate the remaining lifetimes for short period planets, examine the survival of planets around evolving stars, and determine the limits for circularisation of planets with highly eccentric orbits. Our analysis suggests that the survival of circularised planets is strongly affected by the amount of heat dissipated, which is often large enough to lead to runaway orbital inflation and Roche lobe overflow.

Nebular water depletion as the cause of Jupiter's low oxygen abundance

Motivated by recent spectroscopic observations suggesting that atmospheres of some extrasolar giant-planets are carbon-rich, i.e. carbon/oxygen ratio (C/O) $\ge$ 1, we find that the whole set of compositional data for Jupiter is consistent with the hypothesis that it be a carbon-rich giant planet. We show that the formation of Jupiter in the cold outer part of an oxygen-depleted disk (C/O $\sim$1) reproduces the measured Jovian elemental abundances at least as well as the hitherto canonical model of Jupiter formed in a disk of solar composition (C/O = 0.54). The resulting O abundance in Jupiter’s envelope is then moderately enriched by a factor of $\sim$2 $\times$ solar (instead of $\sim$7 $\times$ solar) and is found to be consistent with values predicted by thermochemical models of the atmosphere. That Jupiter formed in a disk with C/O $\sim$1 implies that water ice was heterogeneously distributed over several AU beyond the snow line in the primordial nebula and that the fraction of water contained in icy planetesimals was a strong function of their formation location and time. The Jovian oxygen abundance to be measured by NASA’s Juno mission en route to Jupiter will provide a direct and strict test of our predictions.

The stellar wind cycles and planetary radio emission of the Tau Boo system

Tau Boo is an intriguing planet-host star that is believed to undergo magnetic cycles similar to the Sun, but with a duration that is about one order of magnitude smaller than that of the solar cycle. With the use of observationally derived surface magnetic field maps, we simulate the magnetic stellar wind of Tau Boo by means of three-dimensional MHD numerical simulations. As the properties of the stellar wind depend on the particular characteristics of the stellar magnetic field, we show that the wind varies during the observed epochs of the cycle. Although the mass loss-rates we find (~2.7e-12 Msun/yr) vary less than 3 per cent during the observed epochs of the cycle, our derived angular momentum loss-rates vary from 1.1 to 2.2e32erg. The spin-down times associated to magnetic braking range between 39 and 78Gyr. We also compute the emission measure from the (quiescent) closed corona and show that it remains approximately constant through these epochs at a value of ~10^{50.6} cm^{-3}. This suggests that a magnetic cycle of Tau Boo may not be detected by X-ray observations. We further investigate the interaction between the stellar wind and the planet by estimating radio emission from the hot-Jupiter that orbits at 0.0462 au from Tau Boo. By adopting reasonable hypotheses, we show that, for a planet with a magnetic field similar to Jupiter (~14G at the pole), the radio flux is estimated to be about 0.5-1 mJy, occurring at a frequency of 34MHz. If the planet is less magnetised (field strengths roughly <4G), detection of radio emission from the ground is unfeasible due to the Earth's ionospheric cutoff. According to our estimates, if the planet is more magnetised than that and provided the emission beam crosses the observer line-of-sight, detection of radio emission from Tau Boo b is only possible by ground-based instruments with a noise level of < 1 mJy, operating at low frequencies.

Water ice in the Kuiper belt

We examine a large collection of low resolution near-infrared spectra of Kuiper belt objects and centaurs in an attempt to understand the presence of water ice in the Kuiper belt. We find that water ice on the surface of these objects occurs in three separate manners: (1) Haumea family members uniquely show surfaces of nearly pure water ice, presumably a consequence of the fragmentation of the icy mantle of a larger differentiated proto-Haumea; (2) large objects with absolute magnitudes of $H<3$ (and a limited number to H=4.5) have surface coverings of water ice – perhaps mixed with ammonia – that appears to be related to possibly ancient cryovolcanism on these large objects; and (3) smaller KBOs and centaurs which are neither Haumea family members nor cold-classical KBOs appear to divide into two families (which we refer to as "neutral" and "red"), each of which is a mixture of a common nearly-neutral component and either a slightly red or very red component that also includes water ice. A model suggesting that the difference between neutral and red objects is due to formation in an early compact solar system either inside or outside, respectively, of the ~20 AU methanol evaporation line is supported by the observation that methanol is only detected on the reddest objects, which are those which would be expected to have the most of the methanol containing mixture.

Year 3 LUNAR Annual Report to the NASA Lunar Science Institute [Replacement]

The Lunar University Network for Astrophysics Research (LUNAR) is a team of researchers and students at leading universities, NASA centers, and federal research laboratories undertaking investigations aimed at using the Moon as a platform for space science. LUNAR research includes Lunar Interior Physics & Gravitation using Lunar Laser Ranging (LLR), Low Frequency Cosmology and Astrophysics (LFCA), Planetary Science and the Lunar Ionosphere, Radio Heliophysics, and Exploration Science. The LUNAR team is exploring technologies that are likely to have a dual purpose, serving both exploration and science. There is a certain degree of commonality in much of LUNAR’s research. Specifically, the technology development for a lunar radio telescope involves elements from LFCA, Heliophysics, Exploration Science, and Planetary Science; similarly the drilling technology developed for LLR applies broadly to both Exploration and Lunar Science.

Year 3 LUNAR Annual Report to the NASA Lunar Science Institute

The Lunar University Network for Astrophysics Research (LUNAR) is a team of researchers and students at leading universities, NASA centers, and federal research laboratories undertaking investigations aimed at using the Moon as a platform for space science. LUNAR research includes Lunar Interior Physics & Gravitation using Lunar Laser Ranging (LLR), Low Frequency Cosmology and Astrophysics (LFCA), Planetary Science and the Lunar Ionosphere, Radio Heliophysics, and Exploration Science. The LUNAR team is exploring technologies that are likely to have a dual purpose, serving both exploration and science. There is a certain degree of commonality in much of LUNAR’s research. Specifically, the technology development for a lunar radio telescope involves elements from LFCA, Heliophysics, Exploration Science, and Planetary Science; similarly the drilling technology developed for LLR applies broadly to both Exploration and Lunar Science.

Mantle Dynamics in Super-Earths: Post-Perovskite Rheology and Self-Regulation of Viscosity

Simple scalings suggest that super-Earths are more likely than an equivalent Earth-sized planet to be undergoing plate tectonics. Generally, viscosity and thermal conductivity increase with pressure while thermal expansivity decreases, resulting in lower convective vigor in the deep mantle. According to conventional thinking, this might result in no convection in a super-Earth’s deep mantle. Here we evaluate this. First, we here extend the density functional theory (DFT) calculations of post-perovskite activation enthalpy of to a pressure of 1 TPa. The activation volume for diffusion creep becomes very low at very high pressure, but nevertheless for the largest super-Earths the viscosity along an adiabat may approach 1030 Pa s in the deep mantle. Second, we use these calculated values in numerical simulations of mantle convection and lithosphere dynamics of planets with up to ten Earth masses. The models assume a compressible mantle including depth-dependence of material properties and plastic yielding induced plate tectonics. Results confirm the likelihood of plate tectonics and show a novel self-regulation of deep mantle temperature. The deep mantle is not adiabatic; instead internal heating raises the temperature until the viscosity is low enough to facilitate convective loss of the radiogenic heat, which results in a super-adiabatic temperature profile and a viscosity increase with depth of no more than ~3 orders of magnitude, regardless of the viscosity increase that is calculated for an adiabat. Convection in large super-Earths is characterised by large upwellings and small, time-dependent downwellings. If a super-Earth was extremely hot/molten after its formation, it is thus likely that even after billions of years its deep interior is still extremely hot and possibly substantially molten with a “super basal magma ocean” – a larger version of (Labrosse et al., 2007).

Mapping Earth-analogs from Photometric Variability: Spin-Orbit Tomography for Planets in Inclined Orbits

Aiming at obtaining detailed information of surface environment of Earth-analogs, Kawahara & Fujii 2011 proposed an inversion technique of annual scattered light curves named the spin-orbit tomography (SOT), which enables one to sketch a 2-dimensional albedo map from annual variation of the disk-integrated scattered light, and demonstrated the method with a planet in a face-on orbit. We extend it to be applicable to general geometric configurations, including low-obliquity planets like the Earth in inclined orbits. We simulate light curves of the Earth in an inclined orbit in three photometric bands (0.4-0.5um, 0.6-0.7um, and 0.8-0.9um) and show that the distribution of clouds, snow, and continents are retrieved with the aid of the SOT. We also demonstrate the SOT by applying it to an upright Earth, a tidally-locked Earth, and Earth-analogs with ancient continental configurations. The inversion is model-independent in the sense that we do not assume specific albedo models when mapping the surface, and hence applicable in principle to any kind of inhomogeneity. This method can potentially serve as a unique tool to investigate the exohabitats/exoclimes of Earth-analogs.

On the low-frequency boundary of Sun-generated MHD turbulence in the slow solar wind

New aspects of the slow solar wind turbulent heating and acceleration are investigated. A physical meaning of the lower boundary of the Alfv\’en wave turbulent spectra in the solar atmosphere and the solar wind is studied and the significance of this natural parameter is demonstrated. Via an analytical and quantitative treatment of the problem we show that a truncation of the wave spectra from the lower frequency side, which is a consequence of the solar magnetic field structure and its cyclic changes, results in a significant reduction of the heat production and acceleration rates. An appropriate analysis is presented regarding the link of the considered problem with existing observational data and slow solar wind initiation scenarios.

Non-detection of transits of the super-Earth HD 97658b with MOST photometry

The radial velocity-discovered exoplanet HD 97658b was recently announced to transit, with a derived planetary radius of 2.93 $\pm$ 0.28 R$_{\oplus}$. As a transiting super-Earth orbiting a bright star, this planet would make an attractive candidate for additional observations, including studies of its atmospheric properties. We present and analyze follow-up photometric observations of the HD 97658 system acquired with the MOST space telescope. Our results show no transit with the depth and ephemeris reported in the announcement paper. For the same ephemeris, we rule out transits for a planet with radius larger than 1.87 R$_{\oplus}$. We also report new radial velocity measurements which continue to support the existence of an exoplanet with a period of 9.5 days, and obtain improved orbital parameters.

Simultaneous Swift X-ray and UV views of comet C/2007 N3 (Lulin)

We present an analysis of simultaneous X-Ray and UV observations ofcomet C/2007 N3 (Lulin) taken on three days between January 2009 and March 2009 using the Swift observatory. For our X-ray observations, we used basic transforms to account for the movement of the comet to allow the combination of all available data to produce an exposure-corrected image. We fit a simple model to the extracted spectrum and measured an X-ray flux of 4.3+/-1.3 * 10^-13 ergs cm-2 s-1 in the 0.3 to 1.0 keV band. In the UV, we acquired large-aperture photometry and used a coma model to derive water production rates given assumptions regarding the distribution of water and its dissociation into OH molecules about the comet’s nucleus. We compare and discuss the X-ray and UV morphology of the comet. We show that the peak of the cometary X-ray emission is offset sunward of the UV peak emission, assumed to be the nucleus, by approximately 35,000 km. The offset observed, the shape of X-ray emission and the decrease of the X-ray emission comet-side of the peak, suggested that the comet was indeed collisionally thick to charge exchange, as expected from our measurements of the comet’s water production rate (6–8 10^28 mol. s-1). The X-ray spectrum is consistent with solar wind charge exchange emission, and the comet most likely interacted with a solar wind depleted of very highly ionised oxygen. We show that the measured X-ray lightcurve can be very well explained by variations in the comet’s gas production rates, the observing geometry and variations in the solar wind flux.

Dissipative Divergence of Resonant Orbits

A considerable fraction of multi-planet systems discovered by the observational surveys of extrasolar planets reside in mild proximity to first-order mean motion resonances. However, the relative remoteness of such systems from nominal resonant period ratios (e.g. 2:1, 3:2, 4:3) has been interpreted as evidence for lack of resonant interactions. Here we show that a slow divergence away from exact commensurability is a natural outcome of dissipative evolution and demonstrate that libration of critical angles can be maintained tens of percent away from nominal resonance. We construct an analytical theory for the long-term dynamical evolution of dissipated resonant planetary pairs and confirm our calculations numerically. Collectively, our results suggest that a significant fraction of the near-commensurate extrasolar planets are in fact resonant and have undergone significant dissipative evolution.

Chameleon effect and the Pioneer anomaly [Cross-Listing]

The possibility that the apparent anomalous acceleration of the Pioneer 10 and 11 spacecraft may be due, at least in part, to a chameleon field effect is examined. A small spacecraft, with no thin shell, can have a more pronounced anomalous acceleration than a large compact body, such as a planet, having a thin shell. The chameleon effect seems to present a natural way to explain the differences seen in deviations from pure Newtonian gravity for a spacecraft and for a planet, and appears to be compatible with the basic features of the Pioneer anomaly, including the appearance of a jerk term. However, estimates of the size of the chameleon effect indicate that its contribution to the anomalous acceleration is negligible. We conclude that any inverse-square component in the anomalous acceleration is more likely caused by an unmodelled reaction force from solar-radiation pressure, rather than a chameleon field effect.

The remarkable solar twin HIP 56948: a prime target in the quest for other Earths

We study HIP 56948, the best solar twin known to date, to determine with an unparalleled precision how similar is to the Sun in its physical properties, chemical composition and planet architecture. We explore whether the abundances anomalies may be due to pollution from stellar ejecta or to terrestrial planet formation. We perform a differential abundance analysis (both in LTE and NLTE) using high resolution (R = 100,000) high S/N (600) Keck HIRES spectra of the Sun and HIP 56948. We use precise radial velocity data from the McDonald and Keck observatories to search for planets around this star. We achieve a precision of sigma = 0.003 dex for several elements. Including errors in stellar parameters the total uncertainty is as low as sigma = 0.005 dex (1 %), which is unprecedented in elemental abundance studies. The similarities between HIP 56948 and the Sun are astonishing. HIP 56948 is only 17+/-7 K hotter than the Sun, and log g, [Fe/H] and microturbulence are only +0.02+/-0.02 dex, +0.02+/-0.01 dex and +0.01+/-0.01 km/s higher than solar, respectively. HIP 56948 has a mass of 1.02+/-0.02M_Sun and is 1 Gyr younger than the Sun. Both stars show a chemical abundance pattern that differs from most solar twins. The trend with T_cond in differential abundances (twins – HIP56948) can be reproduced very well by adding 3 M_Earth of a mix of Earth and meteoritic material, to the convection zone of HIP 56948. From our radial velocity monitoring we find no indications of giant planets interior to or within the habitable zone of HIP 56948. We conclude that HIP 56948 is an excellent candidate to host a planetary system like our own, including the possible presence of inner terrestrial planets. Its striking similarity to the Sun and its mature age makes HIP 56948 a prime target in the quest for other Earths and SETI endeavors.

Starbursts and High-Redshift Galaxies are Radioactive: High Abundances of ^{26}Al and Other Short Lived Radionuclides

Short lived radionuclides (SLRs) like ^{26}Al are synthesized by massive stars and are a byproduct of star formation. The abundances of SLRs in the gas of a star-forming galaxy is inversely proportional to its gas consumption time. The rapid evolution of specific star formation rate (SSFR) of normal galaxies implies they had mean SLR abundances ~10 times higher at z = 2. During the epoch of Solar System formation, the mean SLR abundances of the Galaxy were twice as high as at present, if SLR yields from massive stars do not depend on metallicity. If SLRs are well-mixed with the gas of galaxies, the high SSFRs of normal galaxies can partly explain the elevated abundance of SLRs like ^{60}Fe and ^{26}Al in the early Solar System. Starburst galaxies have much higher SSFRs still, and would have enormous mean abundances of ^{26}Al (^{26}Al/^{27}Al ~ 10^-3 for Solar metallicity gas). The main uncertainty is whether the SLRs are mixed with the molecular gas: they may decay before propagating from their origin sites, or be blown out by starburst winds. I show the enhanced ^{26}Al of starbursts can maintain moderate ionization rates (10^-18 – 10^-17 s^-1), possibly dominating ionization in dense clouds not penetrated by cosmic rays. Similar ionization rates would be maintained in protoplanetary disks of starbursts, and the radiogenic heating of planetesimals would likewise be much higher. In this way, galaxy evolution can affect the geological history of planetary systems.

Resonant Repulsion of Kepler Planet Pairs [Replacement]

Planetary systems discovered by the Kepler space telescope exhibit an intriguing feature. While the period ratios of adjacent low-mass planets appear largely random, there is a significant excess of pairs that lie just wide of resonances and a deficit on the near side. We demonstrate that this feature naturally arises when two near-resonant planets interact in the presence of weak dissipation that damps eccentricities. The two planets repel each other as orbital energy is lost to heat. This moves near-resonant pairs just beyond resonance, by a distance that reflects the integrated dissipation they experienced over their lifetimes. We find that the observed distances may be explained by tidal dissipation if tides are efficient (tidal quality factor ~10). Once the effect of resonant repulsion is accounted for, the initial orbits of these low mass planets show little preference for resonances. This is a strong constraint on their origin.

Resonant Repulsion of Kepler Planet Pairs

Planetary systems discovered by the Kepler space telescope exhibit an intriguing feature. While the period ratios of adjacent low-mass planets appear largely random, there is a significant excess of pairs that lie just wide of resonances and a deficit on the near side. We demonstrate that this feature naturally arises when two near-resonant planets interact in the presence of weak dissipation that damps eccentricities. The two planets repel each other as orbital energy is lost to heat. This moves near-resonant pairs just beyond resonance, by a distance that reflects the integrated dissipation they experienced over their lifetimes. We find that the observed distances may be explained by tidal dissipation if tides are efficient (tidal quality factor ~10). Once the effect of resonant repulsion is accounted for, the initial orbits of these low mass planets show little preference for resonances. This is a strong constraint on their origin.

An Analytic Method to determine Habitable Zones for S-Type Planetary Orbits in Binary Star Systems

With more and more extrasolar planets discovered in and around binary star systems, questions concerning the determination of the classical Habitable Zone arise. Do the radiative and gravitational perturbations of the second star influence the extent of the Habitable Zone significantly, or is it sufficient to consider the host-star only? In this article we investigate the implications of stellar companions with different spectral types on the insolation a terrestrial planet receives orbiting a Sun-like primary. We present time independent analytical estimates and compare these to insolation statistics gained via high precision numerical orbit calculations. Results suggest a strong dependence of permanent habitability on the binary’s eccentricity, as well as a possible extension of Habitable Zones towards the secondary in close binary systems.

Rotation-stimulated structures in the CN and C3 comae of comet 103P/Hartley 2 around the EPOXI encounter

In late 2010 a Jupiter Family comet 103P/Hartley 2 was a subject of an intensive world-wide investigation. On UT October 20.7 the comet approached the Earth within only 0.12 AU, and on UT November 4.6 it was visited by NASA’s EPOXI spacecraft. We joined this international effort and organized an observing campaign. The images of the comet were obtained through narrowband filters using the 2-m telescope of the Rozhen National Astronomical Observatory. They were taken during 4 nights around the moment of the EPOXI encounter. Image processing methods and periodicity analysis techniques were used to reveal transient coma structures and investigate their repeatability and kinematics. We observe shells, arc-, jet- and spiral-like patterns, very similar for the CN and C3 comae. The CN features expanded outwards with the sky-plane projected velocities between 0.1 to 0.3 km/s. A corkscrew structure, observed on November 6, evolved with a much higher velocity of 0.66 km/s. Photometry of the inner coma of CN shows variability with a period of 18.32+/-0.30 h (valid for the middle moment of our run, UT 2010 Nov. 5.0835), which we attribute to the nucleus rotation. This result is fully consistent with independent determinations around the same time by other teams. The pattern of repeatability is, however, not perfect, which is understendable given the suggested excitation of the rotation state, and the variability detected in CN correlates well with the cyclic changes in HCN, but only in the active phases. The revealed coma structures, along with the snapshot of the nucleus orientation obtained by EPOXI, let us estimate the spin axis orientation. We obtained RA=122 deg, Dec=+16 deg (epoch J2000.0), neglecting at this point the rotational excitation.

Magnitude and size distribution of long-period comets in Earth-crossing or approaching orbits

We analyse the population of near-Earth Long-Period Comets (LPCs) (perihelion distances q 10^3 yr). We have considered the sample of LPCs discovered during the period 1900-2009 and their estimated absolute total visual magnitudes H. For the period 1900-1970 we have relied upon historical estimates of absolute total magnitudes, while for the more recent period 1970-2009 we have made our own estimates of H based on Green’s photometric data base and IAU Circulars. We have also used historical records for the sample of brightest comets (H < 4.5) covering the period: 1500-1899, based mainly on Vsekhsvyatskii, Hasegawa and Kronk catalogues. We find that the cumulative distribution of H can be represented by a three-modal law of the form log_{10}N_{<H} = C + alpha times H, where the C's are constants for the different legs, and alpha \simeq 0.28 +/- 0.10 for H < 4.0, alpha \simeq 0.56 +/- 0.10 for 4.0 <= H < 5.8, and alpha \simeq 0.20 +/- 0.02 for 5.8 <= H <8.6. The large increase of the slope of the second leg of the H-distribution might be at least partially attributed to splitting of comet nuclei leading to the creation of two or more daughter comets. The cumulative H-distribution tends to flatten for comets fainter than H <= 8.6. LPCs fainter than H <= 12 (or diametres D \lesssim 0.5 km) are extremely rare, despite several sky surveys of near-Earth objects implemented during the last couple of decades, suggesting a minimum size for a LPC to remain active. We also find that about 30 % of all LPCs with q < 1.3 AU are new (original bound energies 0 < E_{or} < 10^{-4} AU^{-1}), and that among the new comets about half come from the outer Oort cloud (energies 0 \lesssim E_{or} \lesssim 0.3 times 10^{-4} AU^{-1}), and the other half from the inner Oort cloud (energies 0.3 times 10^{-4} \lesssim E_{or} \lesssim 10^{-4}AU^{-1}).

A Correlation Between the Eclipse Depths of Kepler Gas Giant Candidates and the Metallicities of their Parent Stars

Previous studies of the interior structure of transiting exoplanets have shown that the heavy element content of gas giants increases with host star metallicity. Since metal-poor planets are less dense and have larger radii than metal-rich planets of the same mass, one might expect that metal-poor stars host a higher proportion of gas giants with large radii than metal-rich stars. Here I present evidence for a negative correlation at the 2.3-sigma level between eclipse depth and stellar metallicity in the Kepler gas giant candidates. Based on Kendall’s tau statistics, the probability that eclipse depth depends on star metallicity is 0.981. The correlation is consistent with planets orbiting low-metallicity stars being, on average, larger in comparison with their host stars than planets orbiting metal-rich stars. Furthermore, since metal-rich stars have smaller radii than metal-poor stars of the same mass and age, a uniform population of planets should show a rise in median eclipse depth with [M/H]. The fact that I find the opposite trend indicates that substantial changes in gas giant interior structure must accompany increasing [M/H]. I investigate whether the known scarcity of giant planets orbiting low-mass stars could masquerade as an eclipse depth-metallicity correlation, given the degeneracy between metallicity and temperature for cool stars in the Kepler Input Catalog. While the eclise depth-metallicity correlation is not yet on firm statistical footing and will require spectroscopic [Fe/H] measurements for validation, it is an intriguing window into how the interior structure of planets and even the planet formation mechanism may be changing with Galactic chemical evolution.

Supercriticality to subcriticality in dynamo transitions [Cross-Listing]

In this paper, we present a three-mode dynamo model which describes both supercritical and subcritical dynamo transitions. The nature of dynamo transition changes from supercritical to subcritical as the magnetic Prandtl number is decreased, consistent with the numerical results in the spherical-shell and the Taylor-Green dynamo. We also perform a detailed analysis of the hysteresis zone of a subcritical dynamo using direct numerical simulations. Numerical simulation and model analysis show that the sets of initial conditions, called the basin of attraction, of the no-dynamo and the dynamo states are separated by an unstable manifold.

Star Hoppers: Planet Instability and Capture in Evolving Binary Systems

Many planets are observed in stellar binary systems, and their frequency may be comparable to that of planetary systems around single stars. Binary stellar evolution in such systems influences the dynamical evolution of the resident planets. Here we study the evolution of a single planet orbiting one star in an evolving binary system. We find that stellar evolution can trigger dynamical instabilities that drive planets into chaotic orbits. This instability leads to planet-star collisions, exchange of the planet between the binary stars (“star-hoppers”), and ejection of the planet from the system. The means by which planets can be recaptured is similar to the pull-down capture mechanism for irregular solar system satellites. Because planets often suffer close encounters with the primary on the asymptotic giant branch, captures during a collision with the stellar envelope are also possible. Such capture could populate the habitable zone around white dwarfs.

Time-resolved Ultraviolet Spectroscopy of the M-dwarf GJ 876 Exoplanetary System

Extrasolar planets orbiting M-stars may represent our best chance to discover habitable worlds in the coming decade. The ultraviolet spectrum incident upon both Earth-like and Jovian planets is critically important for proper modeling of their atmospheric heating and chemistry. In order to provide more realistic inputs for atmospheric models of planets orbiting low-mass stars, we present new near- and far-ultraviolet (NUV and FUV) spectroscopy of the M-dwarf exoplanet host GJ 876 (M4V). Using the COS and STIS spectrographs aboard the Hubble Space Telescope, we have measured the 1150-3140A spectrum of GJ 876. We have reconstructed the stellar HI LyA emission line profile, and find that the integrated LyA flux is roughly equal to the rest of the integrated flux (1150-1210A + 1220-3140A) in the entire ultraviolet bandpass (F(LyA)/F(FUV+NUV) ~0.7). This ratio is ~ 2500x greater than the solar value. We describe the ultraviolet line spectrum and report surprisingly strong fluorescent emission from hot H2 (T(H2) > 2000 K). We show the light-curve of a chromospheric + transition region flare observed in several far-UV emission lines, with flare/quiescent flux ratios >= 10. The strong FUV radiation field of an M-star (and specifically LyA) is important for determining the abundance of O2 — and the formation of biomarkers — in the lower atmospheres of Earth-like planets in the habitable zones of low-mass stars.

Mantle geoneutrinos in KamLAND and Borexino [Cross-Listing]

The KamLAND and Borexino experiments have observed, each at ~4 sigma level, signals of electron antineutrinos produced in the decay chains of thorium and uranium in the Earth’s crust and mantle (Th and U geoneutrinos). Various pieces of geochemical and geophysical information allow an estimation of the crustal geoneutrino flux components with relatively small uncertainties. The mantle component may then be inferred by subtracting the estimated crustal flux from the measured total flux. To this purpose, we analyze in detail the experimental Th and U geoneutrino event rates in KamLAND and Borexino, including neutrino oscillation effects. We estimate the crustal flux at the two detector sites, using state-of-the-art information about the Th and U distribution on global and local scales. We find that crust-subtracted signals show hints of a residual mantle component, emerging at ~2.4 sigma level by combining the KamLAND and Borexino data. The inferred mantle flux slightly favors scenarios with relatively high Th and U abundances, within +-1 sigma uncertainties comparable to the spread of predictions from recent mantle models.

The study of ground-level ozone in Kiev and its impact on public health [Replacement]

Ground-level ozone in Kiev for an episode of its high concentration in August 2000 was simulated with the model of the urban air pollution UAM-V (Urban Airshed Model). The study of total ozone over Kiev and its concentration changes with height in the troposphere is made on the basis of ground-based observations with the infrared Fourier spectrometer at the Main Astronomical Observatory of National Academy of Sciences of Ukraine as a part of the ESA-NIVR-KNMI no 2907. In 2008 the satellite Aura-OMI data OMO3PR on the atmosphere ozone profiles became available. Beginning in 2005, these data include the ozone concentration in the lower layer of the atmosphere and can be used for the evaluation of the ground-level ozone concentrations in all cities of Ukraine. Some statistical investigation of ozone air pollution in Kiev and medical statistics data on respiratory system was carried out with the application of the “Statistica” package. The regression analysis, prognostic regression simulation, and retrospective prognosis of the epidemiological situation with respect to respiratory system pathologies in Kiev during 2000-2007 were performed.

A study of ozone in the surface layer of Kiev and its impact on the human health

Ground-level ozone in Kiev for an episode of its high concentration in August 2000 was simulated with the model of the urban air pollution UAM-V (Urban Airshed Model). The study of total ozone over Kiev and its concentration changes with height in the troposphere is made on the basis of ground-based observations with the infrared Fourier spectrometer at the Main Astronomical Observatory of National Academy of Sciences of Ukraine as a part of the ESA-NIVR-KNMI no 2907. In 2008 the satellite Aura-OMI data OMO3PR on the atmosphere ozone profiles became available. Beginning in 2005, these data include the ozone concentration in the lower layer of the atmosphere and can be used for the evaluation of the ground-level ozone concentrations in all cities of Ukraine. Some statistical investigation of ozone air pollution in Kiev and medical statistics data on respiratory system was carried out with the application of the “Statistica” package. The regression analysis, prognostic regression simulation, and retrospective prognosis of the epidemiological situation with respect to respiratory system pathologies in Kiev during 2000-2007 were performed.

The phase equilibrium in a Lennard Jones fluid: possible applications in astrophysics [Cross-Listing]

Using standard thermodynamics and previous results of the author, this paper aims to discuss the conditions for phase equilibrium in a Lennard-Jones fluid. Possibilities of astrophysical applications of the results obtained here are discussed to some extent.

Transfer of Life-Bearing Meteorites from Earth to Other Planets

The probability is investigated that the meteorites originating on Earth are transferred to other planets in our Solar System and to extra solar planets. We take the collisional Chicxulub crater event, and material that was ejected as an example of Earth-origin meteors. If we assume the appropriate size of the meteorites as 1cm in diameter, the number of meteorites to reach the exoplanet system (further than 20 ly) would be much greater than one. We have followed the ejection and capture rates estimated by Melosh (2003) and the discussion by Wallis and Wickramasinghe (2004). If we consider the possibility that the fragmented ejecta (smaller than 1cm) are accreted to comets and other icy bodies, then buried fertile material could make the interstellar journey throughout Galaxy. If life forms inside remain viable, this would be evidence of life from Earth seeding other planets. We also estimate the transfer velocity of the micro-organisms in the interstellar space. In some assumptions, it could be estimated that, if life has originated $10^{10}$\ years ago anywhere in our Galaxy as theorized by Joseph and Schild (2010a, b), it will have since propagated throughout our Galaxy and could have arrived on Earth by 4.6 billion years ago. Organisms disperse.

Tidal evolution of exo-planetary systems: WASP-50, GJ 1214 and CoRoT-7

We perform numerical simulations to investigate tidal evolution of two single-planet systems, that is, WASP-50 and GJ 1214 and a two-planet system CoRoT-7. The results of orbital evolution show that tidal decay and circularization may play a significant role in shaping their final orbits, which is related to the initial orbital data in the simulations. For GJ 1214 system, different cases of initial eccentricity are also considered as only an upper limit of its eccentricity (0.27) is shown, and the outcome suggests a possible maximum initial eccentricity (0.4) in the adopted dynamical model. Moreover, additional runs with alternative values of dissipation factor $Q^\prime_1$ are carried out to explore tidal evolution for GJ 1214b, and these results further indicate that the real $Q^\prime_1$ of GJ 1214b may be much larger than its typical value, which may reasonably suggest that GJ 1214b bears a present-day larger eccentricity, undergoing tidal circularization at a slow rate. For the CoRoT-7 system, tidal forces make two planets migrating towards their host star as well as producing tidal circularization, and in this process tidal effects and mutual gravitational interactions are coupled with each other. Various scenarios of the initial eccentricity of the outer planet have also been done to investigate final planetary configuration. Tidal decay arising from stellar tides may still work for each system as the eccentricity decreases to zero, and this is in association with the remaining lifetime of each planet used to predict its future.

Theoretical Spectra of Terrestrial Exoplanet Surfaces

We investigate spectra of airless rocky exoplanets with a theoretical framework that self-consistently treats reflection and thermal emission. We find that a silicate surface on an exoplanet is spectroscopically detectable via prominent Si-O features in the thermal emission bands of 7 – 13 \mu m and 15 – 25 \mu m. The variation of brightness temperature due to the silicate features can be up to 20 K for an airless Earth analog, and the silicate features are wide enough to be distinguished from atmospheric features with relatively high-resolution spectra. The surface characterization thus provides a method to unambiguously identify a rocky exoplanet. Furthermore, identification of specific rocky surface types is possible with the planet’s reflectance spectrum in near-infrared broad bands. A key parameter to observe is the difference between K band and J band geometric albedos (A_g (K)-A_g (J)): A_g (K)-A_g (J) > 0.2 indicates that more than half of the planet’s surface has abundant mafic minerals, such as olivine and pyroxene, in other words primary crust from a magma ocean or high-temperature lavas; A_g (K)-A_g (J) < -0.09 indicates that more than half of the planet's surface is covered or partially covered by water ice or hydrated silicates, implying extant or past water on its surface. Also, surface water ice can be specifically distinguished by an H-band geometric albedo lower than the J-band geometric albedo. The surface features can be distinguished from possible atmospheric features with molecule identification of atmospheric species by transmission spectroscopy. We therefore propose that mid-infrared spectroscopy of exoplanets may detect rocky surfaces, and near-infrared spectrophotometry may identify ultramafic surfaces, hydrated surfaces and water ice.

Anelastic tidal dissipation in multi-layer planets

Earth-like planets have viscoelastic mantles, whereas giant planets may have viscoelastic cores. The tidal dissipation of such solid regions, gravitationally perturbed by a companion body, highly depends on their rheology and on the tidal frequency. Therefore, modelling tidal interactions presents a high interest to provide constraints on planets’ properties and to understand their history and their evolution, in our Solar System or in exoplanetary systems. We examine the equilibrium tide in the anelastic parts of a planet whatever the rheology, taking into account the presence of a fluid envelope of constant density. We show how to obtain the different Love numbers that describe its tidal deformation. Thus, we discuss how the tidal dissipation in solid parts depends on the planet’s internal structure and rheology. Finally, we show how the results may be implemented to describe the dynamical evolution of planetary systems. The first manifestation of the tide is to distort the shape of the planet adiabatically along the line of centers. Then, the response potential of the body to the tidal potential defines the complex Love numbers whose real part corresponds to the purely adiabatic elastic deformation, while its imaginary part accounts for dissipation. This dissipation is responsible for the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which we derive the characteristic evolution times. The rate at which the system evolves depends on the physical properties of tidal dissipation, and specifically on how the shear modulus varies with tidal frequency, on the radius and also the rheological properties of the solid core. The quantification of the tidal dissipation in solid cores of giant planets reveals a possible high dissipation which may compete with dissipation in fluid layers.

(1173) Anchises - Thermophysical and Dynamical Studies of a Dynamically Unstable Jovian Trojan

We have performed detailed thermophysical and dynamical modelling of Jovian Trojan (1173) Anchises. Our results reveal a most unusual object. By examining observational data taken by IRAS, Akari and WISE between 11.5 and 60 microns, along with variations in its optical lightcurve, we find Anchises is most likely an elongated body, with an axes-ratio of ~1.4. This yields calculated best-fit dimensions of 170×121x121km (an equivalent diameter of 136+18/-11km). We find the observations are best fit by Anchises having a retrograde sense of rotation, and an unusually high thermal inertia (25 to 100 Jm-2s-0.5K-1). The geometric albedo is found to be 0.027 (+0.006/-0.007). Anchises therefore has one of the highest published thermal inertias of any object larger than 100km in diameter, at such large heliocentric distances, and is one of the lowest albedo objects ever observed. More observations are needed to see if there is a link between the very shallow phase curve, with almost no opposition effect, and the derived thermal properties for this large Trojan asteroid. Our dynamical investigation of Anchises’ orbit has revealed it to be dynamically unstable on timescales of hundreds of Myr, similar to the unstable Neptunian Trojans 2001 QR322 and 2008 LC18. Unlike those objects, we find that Anchises’ dynamical stability is not a function of its initial orbital elements, the result of the exceptional precision with which its orbit is known. This is the first time that a Jovian Trojan has been shown to be dynamically unstable, and adds weight to the idea that planetary Trojans represent a significant ongoing contribution to the Centaur population, the parents of the short-period comets. The observed instability does not rule out a primordial origin for Anchises, but when taken in concert with the result of our thermophysical analysis, suggest that it would be a fascinating target for future study.

Circumbinary Planet Formation in the Kepler-16 system. I. N-body Simulations

The recently discovered circumbinary planets (Kepler-16 b, Kepler-34 b, Kepler-35 b) represent the first direct evidence of the viability of planet formation in circumbinary orbits. We report on the results of N-body simulations investigating planetesimal accretion in the Kepler-16 b system, focusing on the range of impact velocities under the influence of both stars’ gravitational perturbation and friction from a putative protoplanetary disk. Our results show that planet formation might be effectively inhibited for a large range in semi-major axis (1.75 < a_P < 4 AU), suggesting that the planetary core must have either migrated from outside 4 AU, or formed in situ very close to its current location.

Evidence for 9 planets in the HD 10180 system

We re-analyse the HARPS radial velocities of HD 10180 and calculate the probabilities of models with differing numbers of periodic signals in the data. We test the significance of the seven signals, corresponding to seven exoplanets orbiting the star, in the Bayesian framework and perform comparisons of models with up to nine periodicities. We use posterior samplings and Bayesian model probabilities in our analyses together with suitable prior probability densities and prior model probabilities to extract all the significant signals from the data and to receive reliable uncertainties for the orbital parameters of the six, possibly seven, known exoplanets in the system. According to our results, there is evidence for up to nine planets orbiting HD 10180, which would make this this star a record holder in having more planets in its orbits than there are in the Solar system. We revise the uncertainties of the previously reported six planets in the system, verify the existence of the seventh signal, and announce the detection of two additional statistically significant signals in the data. If of planetary origin, these two additional signals would correspond to planets with minimum masses of 5.1$^{+3.1}_{-3.2}$ and 1.9$^{+1.6}_{-1.8}$ M$_{\oplus}$ on orbits with 67.55$^{+0.68}_{-0.88}$ and 9.655$^{+0.022}_{-0.072}$ days periods (denoted using the 99% credibility intervals), respectively.

The Kilometer-Sized Main Belt Asteroid Population as Revealed by Spitzer

Multi-epoch Spitzer Space Telescope 24 micron data is utilized from the MIPSGAL and Taurus Legacy surveys to detect asteroids based on their relative motion. These infrared detections are matched to known asteroids and rotationally averaged diameters and albedos are derived using the Near Earth Asteroid Model (NEATM) in conjunction with Monte Carlo simulations for 1835 asteroids ranging in size from 0.2 to 143.6 km. A small subsample of these objects was also detected by IRAS or MSX and the single wavelength albedo and diameter fits derived from this data are within 5% of the IRAS and/or MSX derived albedos and diameters demonstrating the robustness of our technique. The mean geometric albedo of the small main belt asteroids in this sample is p_V = 0.138 with a sample standard deviation of 0.105. The albedo distribution of this sample is far more diverse than the IRAS or MSX samples. The cumulative size-frequency distribution of asteroids in the main belt at small diameters is directly derived. Completeness limits of the optical and infrared surveys are discussed.

Directly Imaging Tidally Powered Migrating Jupiters

We show that ongoing direct imaging experiments may detect a new class of long-period, highly luminous, tidally powered extrasolar gas giants. Even though they are hosted by Gyr-”old” main-sequence stars, they can be as “hot” as young Jupiters at ~100 Myr, the prime targets of direct imaging surveys. These planets, with years-long orbits, are presently migrating to “feed” the “hot Jupiters” in steady state. Their existence is expected from a class of “high-e” migration mechanisms, in which gas giants are excited to highly eccentric orbits and then shrink their semi-major axis by factor of ~ 10-100 due to tidal dissipation at successive close periastron passages. The dissipated orbital energy is converted to heat, and if it is deposited deep enough into the planet atmosphere, the planet likely radiates steadily at luminosity ~2-3 orders of magnitude larger than that of our Jupiter during a typical Gyr migration time scale. Their large orbital separations and expected high planet-to-star flux ratios in IR make them potentially accessible to high-contrast imaging instruments on 10m-class telescopes at present and in the near future. A dozen or so such planets are expected to exist around FGK dwarfs within ~50 pc. Long-period planets around nearby stars found by RV are viable candidates to follow up, and in particular, the highly eccentric planet HD 20782b at maximum angular separation ~ 0.08″ is the most promising candidate. Directly imaging these tidally powered Jupiters would enable a direct test of high-e migration mechanisms. Once detected, the luminosity would provide a direct measurement of the migration rate, and together with mass (and possibly radius) estimate, they would serve as a laboratory to study planetary spectral formation and tidal physics.

A New Disintegrative Capture Theory for the Origin of the Moon [Replacement]

The object that resulted in the creation of the Moon started in the same orbital path as Earth around the Sun, but at Earth’s L4. This proto-Moon (PM) was 4 times less massive than the usual Giant Impact (GI) object “Theia” and was captured into Earth orbit. It had a 32% Iron-Nickel-Sulfur core supporting a dynamo, which explains magnetized lunar rocks. Following capture, it was torn apart by tidal forces and its core of iron plastered itself, with some of its rock mantle, on the surface of Earth at a very flat angle (producing the “Late Veneer”). After tidal stripping, the remaining PM rock was driven away from Earth to about 3.8 times Earth’s radius and formed into what is now the Moon. The GI theory has several troubles: The violent collision melts the entire Earth, contrary to geological evidence. The Moon itself also has to condense out of the vapor cloud generated in the collision, but there is evidence that the Moon was not condensed out of vapor. In the new theory, the Moon as we know it may be only 3.8 – 3.9 billion years old, not 4.56 as usually assumed. That is the age of the PM. The minerals in the Moon would be about as old as the Earth, but would have been re-arranged in the capture and temporary disintegration process. If the Moon is as young as suggested, its origin would coincide with the beginning of life on Earth, which is unexplained in the GI theory. The manuscript asks, “Was the Moon Turned Inside-Out” and the answer is “Essentially, Yes.”

A New Disintegrative Capture Theory for the Origin of the Moon

The object that resulted in the creation of the Moon started in the same orbital path as Earth around the Sun, but at Earth’s L4. This proto-Moon (PM) was 4 times less massive than the usual Giant Impact (GI) object “Theia” and was captured into Earth orbit. It had a 32% Iron-Nickel-Sulfur core supporting a dynamo, which explains magnetized lunar rocks. Following capture, it was torn apart by tidal forces and its core of iron plastered itself, with some of its rock mantle, on the surface of Earth at a very flat angle (producing the “Late Veneer”). After tidal stripping, the remaining PM rock was driven away from Earth to about 3.8 times Earth’s radius and formed into what is now the Moon. The GI theory has several troubles: The violent collision melts the entire Earth, contrary to geological evidence. The Moon itself also has to condense out of the vapor cloud generated in the collision, but there is evidence that the Moon was not condensed out of vapor. In the new theory, the Moon as we know it may be only 3.8 – 3.9 billion years old, not 4.26 as usually assumed. That is the age of the PM. The minerals in the Moon would be about as old as the Earth, but would have been re-arranged in the capture and temporary disintegration process. If the Moon is as young as suggested, its origin would coincide with the beginning of life on Earth, which is unexplained in the GI theory. The manuscript asks, “Was the Moon Turned Inside-Out” and the answer is “Essentially, Yes.”

The roles of charge exchange and dissociation in spreading Saturn's neutral clouds [Cross-Listing]

Neutrals sourced directly from Enceladus’s plumes are initially confined to a dense neutral torus in Enceladus’s orbit around Saturn. This neutral torus is redistributed by charge exchange, impact/photodissociation, and neutral-neutral collisions to produce Saturn’s neutral clouds. Here we consider the former processes in greater detail than in previous studies. In the case of dissociation, models have assumed that OH is produced with a single speed of 1 km/s, whereas laboratory measurements suggest a range of speeds between 1 and 1.6 km/s. We show that the high-speed case increases dissociation’s range of influence from 9 to 15 Rs. For charge exchange, we present a new modeling approach, where the ions are followed within a neutral background, whereas neutral cloud models are conventionally constructed from the neutrals’ point of view. This approach allows us to comment on the significance of the ions’ gyrophase at the moment charge exchange occurs. Accounting for gyrophase: (1) has no consequence on the H2O cloud; (2) doubles the local density of OH at the orbit of Enceladus; and (3) decreases the oxygen densities at Enceladus’s orbit by less than 10%. Finally, we consider velocity-dependent, as well as species-dependent cross sections and find that the oxygen cloud produced from charge exchange is spread out more than H2O, whereas the OH cloud is the most confined.

Using Satellites to Probe Extrasolar Planet Formation

Planetary satellites are an integral part of the heirarchy of planetary systems. Here we make two predictions concerning their formation. First, primordial satellites, which have an array of distinguishing characteristics, form only around giant planets. If true, the size and duration of a planetary system’s protostellar nebula, as well as the location of its snow line, can be constrained by knowing which of its planets possess primordial satellites and which do not. Second, all satellites around terrestrial planets form by impacts. If true, this greatly enhances the constraints that can be placed on the history of terrestrial planets by their satellites’ compositions, sizes, and dynamics.

Traditional formation scenarios fail to explain 4:3 mean motion resonances

At least two multi-planetary systems in a 4:3 mean motion resonance have been found by radial velocity surveys. These planets are gas giants and only stable when protected by a resonance. Additionally the Kepler mission has detected at least 4 strong candidate planetary systems with a period ratio close to 4:3. This paper investigates traditional dynamical scenarios for the formation of these systems. We systematically study migration scenarios with both N-body and hydro-dynamic simulations. We investigate scenarios involving the in-situ formation of two planets in resonance. We look at the results from finely tuned planet-planet scattering simulations with gas disk damping. Finally, we investigate a formation scenario involving isolation-mass embryos. Although the combined planet-planet scattering and damping scenario seems promising, none of the above scenarios is successful in forming enough systems in 4:3 resonance with planetary masses similar to the observed ones. This is a negative result but it has important implications for planet formation. Previous studies were successful in forming 2:1 and 3:2 resonances. This is generally believed to be evidence of planet migration. We highlight the main differences between those studies and our failure in forming a 4:3 resonance. We also speculate on more exotic and complicated ideas. These results will guide future investigators toward exploring the above scenarios and alternative mechanisms in a more general framework.

"TNOs are Cool": A survey of the trans-Neptunian region -- VII. Size and surface characteristics of (90377) Sedna and 2010 EK139

We present estimates of the basic physical properties (size and albedo) of (90377) Sedna, a prominent member of the detached trans-Neptunian object population and the recently discovered scattered disk object 2010 EK139, based on the recent observations acquired with the Herschel Space Observatory, within the “TNOs are Cool!” key programme. Our modeling of the thermal measurements shows that both objects have larger albedos and smaller sizes than the previous expectations, thus their surfaces might be covered by ices in a significantly larger fraction. The derived diameter of Sedna and 2010 EK139 are 995 +/- 80 km and 470 +35/-10 km, while the respective geometric albedos are pV 0.32 +/- 0.06 and 0.25 +0.02/-0.05. These estimates are based on thermophysical model techniques.

Strong tidal dissipation in Saturn and constraints on Enceladus' thermal state from astrometry

Tidal interactions between Saturn and its satellites play a crucial role in both the orbital migration of the satellites and the heating of their interiors. Therefore constraining the tidal dissipation of Saturn (here the ratio k2/Q) opens the door to the past evolution of the whole system. If Saturn’s tidal ratio can be determined at different frequencies, it may also be possible to constrain the giant planet’s interior structure, which is still uncertain. Here, we try to determine Saturn’s tidal ratio through its current effect on the orbits of the main moons, using astrometric data spanning more than a century. We find an intense tidal dissipation (k2/Q= (2.3 \pm 0.7) \times 10-4), which is about ten times higher than the usual value estimated from theoretical arguments. As a consequence, eccentricity equilibrium for Enceladus can now account for the huge heat emitted from Enceladus’ south pole. Moreover, the measured k2/Q is found to be poorly sensitive to the tidal frequency, on the short frequency interval considered. This suggests that Saturn’s dissipation may not be controlled by turbulent friction in the fluid envelope as commonly believed. If correct, the large tidal expansion of the moon orbits due to this strong Saturnian dissipation would be inconsistent with the moon formations 4.5 Byr ago above the synchronous orbit in the Saturnian subnebulae. But it would be compatible with a new model of satellite formation in which the Saturnian satellites formed possibly over longer time scale at the outer edge of the main rings. In an attempt to take into account for possible significant torques exerted by the rings on Mimas, we fitted a constant rate da/dt on Mimas semi-major axis, also. We obtained an unexpected large acceleration related to a negative value of da/dt= -(15.7 \pm 4.4) \times 10-15 au/day.

Optical and X-ray Transients from Planet-Star Mergers

We evaluate the prompt observational signatures of the merger between a massive close-in planet (a `hot Jupiter’) and its host star, events with an estimated Galactic rate of ~0.1-1/yr. Depending on the ratio of the mean density of the planet rho_p to that of the star rho_star, a merger results in three possible outcomes. If rho_p/rho_star > 5, then the planet directly plunges below the stellar atmosphere before being disrupted by tidal forces. Dissipation of orbital energy creates a hot wake behind the planet, producing a EUV/soft X-ray transient as the planet sinks below the stellar surface. The peak luminosity L_X ~ 1e36 erg/s is achieved weeks-months prior to merger, after which the stellar surface is enshrouded by an outflow. The final inspiral is accompanied by an optical transient powered by the recombination of hydrogen in the outflow, which peaks at L~1e37-38 erg/s on a timescale ~days. If instead rho_planet/rho_star < 5, then Roche Lobe overflow occurs above the stellar surface. For rho_p/rho_star < 1, mass transfer is stable, resulting the planet being accreted on a relatively slow timescale. However, for 1 < rho_p/rho_star < 5, mass transfer may instead be unstable, resulting in the planet being dynamically disrupted into an accretion disk around the star. Super-Eddington outflows from the disk power an optical transient with L~1e37-38 erg/s and characteristic duration ~week-months. The disk itself becomes visible once the accretion rate become sub-Eddington, resulting in a bolometric brightening and spectral shift to the UV. Optical transients from planet merger events may resemble classical novae, but are distinguished by lower ejecta mass and velocity ~100s km/s, and by hard pre- and post-cursor emission, respectively. Promising search strategies include combined optical, UV, and X-ray surveys of nearby massive galaxies with cadences from days to months.

Sodium Atoms in the Lunar Exotail: Observed Velocity and Spatial Distributions

The lunar sodium tail extends long distances due to radiation pressure on sodium atoms in the lunar exosphere. Our earlier observations measured the average radial velocity of sodium atoms moving down the lunar tail beyond Earth (i.e., near the anti-lunar point) to be $\sim 12.5$ km/s. Here we use the Wisconsin H-alpha Mapper to obtain the first kinematically resolved maps of the intensity and velocity distribution of this emission over a $15 \times 15 \deg$ region on the sky near the anti-lunar point. We present both spatially and spectrally resolved observations obtained over four nights bracketing new Moon in October 2007. The spatial distribution of the sodium atoms is elongated along the ecliptic with the location of the peak intensity drifting $3 \deg$ east along the ecliptic per night. Preliminary modeling results suggest the spatial and velocity distributions in the sodium exotail are sensitive to the near surface lunar sodium velocity distribution. Future observations of this sort along with detailed modeling offer new opportunities to describe the time history of lunar surface sputtering over several days.

Constraints on Galileon-induced precessions from solar system orbital motions [Replacement]

We use latest data from solar system planetary orbital motions to put constraints on some Galileon-induced precessional effects. Due to the Vainshtein mechanism, the Galileon-type spherically symmetric field of a monopole induces a small, screened correction proprtional to \sqrt{r} to its usual r^-1 Newtonian potential which causes a secular precession of the pericenter of a test particle. In the case of our solar system, latest data from Mars allow to constrain the magnitude of such an interaction down to \alpha <= 0.3 level. Another Galileon-type effect which might impact solar system dynamics is due to an unscreened constant gradient induced by the peculiar motion of the Galaxy. The magnitude of such an effect, depending on the different gravitational binding energies of the Sun and the planets, is \xi <= 0.004 from the latest bounds on the supplementary perihelion precession of Saturn.

Constraints on Galileon-induced precessions from solar system orbital motions [Replacement]

We use latest data from solar system planetary orbital motions to put constraints on some Galileon-induced precessional effects. Due to the Vainshtein mechanism, the Galileon-type spherically symmetric field of a monopole induces a small, screened correction proprtional to \sqrt{r} to its usual r^-1 Newtonian potential which causes a secular precession of the pericenter of a test particle. In the case of our solar system, latest data from Mars allow to constrain the magnitude of such an interaction down to \alpha <= 0.3 level. Another Galileon-type effect which might impact solar system dynamics is due to an unscreened constant gradient induced by the peculiar motion of the Galaxy. The magnitude of such an effect, depending on the different gravitational binding energies of the Sun and the planets, is \xi <= 0.004 from the latest bounds on the supplementary perihelion precession of Saturn.

Constraints on Galileon-induced precessions from solar system orbital motions [Cross-Listing]

We use latest data from solar system planetary orbital motions to put constraints on some Galileon-induced precessional effects.

Chondrule Formation in Bow Shocks around Eccentric Planetary Embryos

Recent isotopic studies of Martian meteorites by Dauphas & Pourmond (2011) have established that large (~ 3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules – millimeter-sized igneous inclusions found in meteorites – were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We numerically model the hydrodynamics of the flow, and find that such large bodies retain an atmosphere, with Kelvin-Helmholtz instabilities allowing mixing of this atmosphere with the gas and particles flowing past the embryo. We calculate the trajectories of chondrules flowing past the body, and find that they are not accreted by the protoplanet, but may instead flow through volatiles outgassed from the planet’s magma ocean. In contrast, chondrules are accreted onto smaller planetesimals. We calculate the thermal histories of chondrules passing through the bow shock. We find that peak temperatures and cooling rates are consistent with the formation of the dominant, porphyritic texture of most chondrules, assuming a modest enhancement above the likely solar nebula average value of chondrule densities (by a factor of 10), attributable to settling of chondrule precursors to the midplane of the disk or turbulent concentration. We calculate the rate at which a planetary embryo’s eccentricity is damped and conclude that a single planetary embryo scattered into an eccentric orbit can, over ~ 10e5 years, produce ~ 10e24 g of chondrules. In principle, a small number (1-10) of eccentric planetary embryos can melt the observed mass of chondrules in a manner consistent with all known constraints.

"TNOs are Cool": A survey of the trans-Neptunian region VI. Herschel/PACS observations and thermal modeling of 19 classical Kuiper belt objects

Trans-Neptunian objects (TNO) represent the leftovers of the formation of the Solar System. Their physical properties provide constraints to the models of formation and evolution of the various dynamical classes of objects in the outer Solar System. Based on a sample of 19 classical TNOs we determine radiometric sizes, geometric albedos and beaming parameters. Our sample is composed of both dynamically hot and cold classicals. We study the correlations of diameter and albedo of these two subsamples with each other and with orbital parameters, spectral slopes and colors. We have done three-band photometric observations with Herschel/PACS and we use a consistent method for data reduction and aperture photometry of this sample to obtain monochromatic flux densities at 70.0, 100.0 and 160.0 \mu m. Additionally, we use Spitzer/MIPS flux densities at 23.68 and 71.42 \mu m when available, and we present new Spitzer flux densities of eight targets. We derive diameters and albedos with the near-Earth asteroid thermal model (NEATM). As auxiliary data we use reexamined absolute visual magnitudes from the literature and data bases, part of which have been obtained by ground based programs in support of our Herschel key program. We have determined for the first time radiometric sizes and albedos of eight classical TNOs, and refined previous size and albedo estimates or limits of 11 other classicals. The new size estimates of 2002 MS4 and 120347 Salacia indicate that they are among the 10 largest TNOs known. Our new results confirm the recent findings that there are very diverse albedos among the classical TNOs and that cold classicals possess a high average albedo (0.17 +/- 0.04). Diameters of classical TNOs strongly correlate with orbital inclination in our sample. We also determine the bulk densities of six binary TNOs.

On the Formation Age of the First Planetary System

Recently, it has been observed the extreme metal-poor stars in the Galactic halo, which must be formed just after Pop III objects. On the other hand, the first gas clouds of mass $\sim 10^6 M_{\odot}$ are supposed to be formed at $ z \sim $ 10, 20, and 30 for the $1\sigma$, $2\sigma $ and $3\sigma$, where the density perturbations are assumed of the standard $\Lambda$CDM cosmology. If we could apply this gaussian distribution to the extreme small probability, the gas clouds would be formed at $ z \sim $40, 60, and 80 for the $4\sigma$, $6\sigma$, and $8\sigma$. The first gas clouds within our galaxy must be formed around $z\sim 40$. Even if the gas cloud is metal poor, there is a lot of possibility to form the planets around such stars. The first planetary systems could be formed within $\sim 6\times 10^7$ years after the Big Bang in the universe. Even in our galaxies, it could be formed within $\sim 1.7\times 10^8$ years. It is interesting to wait the observations of planets around metal-poor stars. For the panspermia theory, the origin of life could be expected in such systems.

Lightcurve Survey of V-type Asteroids. I. Observations until Spring 2004

To examine the distribution of rotational rates for chips of asteroid 4 Vesta, lightcurve observation of seven V-type asteroids (2511 Patterson, 2640 Hallstorm, 2653 Principia, 2795 Lapage, 3307 Athabasca, 4147 Lennon, and 4977 Rauthgundis) were performed from fall 2003 to spring 2004. Distribution of spin rates of V-type main-belt asteroids from the past and our observations have three peaks. This result implies that age of catastrophic impact making Vesta family may be not as young as Karin and Iannini families but as old as Eos and Koronis families.

Presolar Grains from Novae: Evidence from Neon and Helium Isotopes in Comet Dust Collections [Replacement]

Presolar grains in meteorites and interplanetary dust particles (IDPs) carry non-solar isotopic signatures pointing to origins in supernovae, giant stars, and possibly other stellar sources. There have been suggestions that some of these grains condensed in the ejecta of classical nova outbursts, but the evidence is ambiguous. We report neon and helium compositions in particles captured on stratospheric collectors flown to sample materials from comets 26P/Grigg-Skjellerup and 55P/Tempel-Tuttle that point to condensation of their gas carriers in the ejecta of a neon (ONe) nova. The absence of detectable 3He in these particles indicates space exposure to solar wind (SW) irradiation of a few decades at most, consistent with origins in cometary dust streams. Measured 4He/20Ne, 20Ne/22Ne, 21Ne/22Ne and 20Ne/21Ne isotope ratios, and a low upper limit on 3He/4He, are in accord with calculations of nucleosynthesis in neon nova outbursts. Of these, the uniquely low 4He/20Ne and high 20Ne/22Ne ratios are the most diagnostic, reflecting the large predicted 20Ne abundances in the ejecta of such novae. The correspondence of measured Ne and He compositions in cometary matter with theoretical predictions is evidence for the presence of presolar grains from novae in the early solar system.

Presolar Grains from Novae: Evidence from Neon and Helium Isotopes in Comet Dust Collections

Presolar grains in meteorites and interplanetary dust particles (IDPs) carry non-solar isotopic signatures pointing to origins in supernovae, giant stars, and possibly other stellar sources. There have been suggestions that some of these grains condensed in the ejecta of classical nova outbursts, but the evidence is ambiguous. We report neon and helium compositions in particles captured on stratospheric collectors flown to sample materials from comets 26P/Grigg-Skjellerup and 55P/Tempel-Tuttle that point to condensation of their gas carriers in the ejecta of a neon (ONe) nova. The absence of detectable 3He in these particles indicates space exposure to solar wind (SW) irradiation of a few decades at most, consistent with origins in cometary dust streams. Measured 4He/20Ne, 20Ne/22Ne, 21Ne/22Ne and 20Ne/21Ne isotope ratios, and a low upper limit on 3He/4He, are in accord with calculations of nucleosynthesis in neon nova outbursts. Of these, the uniquely low 4He/20Ne and high 20Ne/22Ne ratios are the most diagnostic, reflecting the large predicted 20Ne abundances in the ejecta of such novae. The correspondence of measured Ne and He compositions in cometary matter with theoretical predictions is evidence for the presence of presolar grains from novae in the early solar system.

Cometary Charge Exchange Diagnostics in UV and X-ray

Since the initial discovery of cometary charge exchange emission, more than 20 comets have been observed with a variety of X-ray and UV observatories. This observational sample offers a broad variety of comets, solar wind environments and observational conditions. It clearly demonstrates that solar wind charge exchange emission provides a wealth of diagnostics, which are visible as spatial, temporal, and spectral emission features. We review the possibilities and limitations of each of those in this contribution.

Evolution of magnetic protection in potentially habitable terrestrial planets

We present a model for the evolution of the magnetic properties of habitable terrestrial planets and their effects on the protection of planetary atmosphere against the erosive action of stellar wind. Using up-to-date thermal evolution models and dynamo scaling laws we predict the evolution of the planetary dipole moment as a function of planetary mass and rotation rate. Combining these results with models for the evolution of the stellar wind, stellar XUV fluxes and planetary exosphere characteristics, we determine the properties of the magnetosphere and the exobase radius in order to estimate the level of atmospheric mass losses. We use this model to evaluate the magnetic protection of the potentially habitable super-Earths GJ 667Cc, Gl 581d and HD 85512b. We confirm that Earth-like planets, even under the highest attainable magnetic field strengths, will lose a significant fraction of their atmospheric volatiles if they are tidally locked in the habitable zone of dM stars, or even if having N/O-rich atmospheres they are in habitable zones closer than $\sim$ 0.8 AU. Similar mass-dependent inner limits have been found for super-Earths $M_p\gtrsim 3 M_\oplus$ that in any case seem to have better chances of preserving their atmospheres even if they are tidally locked. We predict that the atmosphere of GJ 667Cc has probably already been obliterated and it is presently uninhabitable. On the other hand, our model predicts that the atmospheres of Gl 581d and HD 85512b would be well protected by intrinsic magnetic fields, even under the worst expected conditions of stellar aggression. (abrigded abstract).

High-contrast spectroscopy of SCR J1845-6357 B

Spectral characterization of sub-stellar companions is essential to understand their composition and formation processes. However, the large contrast ratio of the brightness of each object to that of its parent star limits our ability to extract a clean spectrum, free from any significant contribution from the star. During the development of the long slit spectroscopy (LSS) mode of IRDIS, the dual-band imager and spectrograph of SPHERE, we proposed a data analysis method to estimate and remove the contributions of the stellar spectrum. This method has never been tested on real data because of the lack of instrumentation capable of combining adaptive optics (AO), coronagraphy, and LSS. Nonetheless, a similar attenuation of the star can be obtained using a particular observing configuration. Test data were acquired using the AO-assisted spectrograph VLT/NACO. We obtained new J- and H-band spectra of SCR J1845-6357 B, a T6 companion to a nearby (3.85\pm0.02 pc) M8 star. This system is a well-suited benchmark as it is relatively wide (~1.0″) with a modest contrast ratio (~4 mag), and a previously published JHK spectrum is available for reference. We demonstrate that (1) our method is efficient at estimating and removing the stellar contribution, (2) it allows to properly recover the spectral shape of the companion, and (3) it is essential to obtain an unbiased estimation of physical parameters. We also show that the slit configuration associated with this method allows us to use long exposure times with high throughput producing high signal-to-noise ratio data. However, the signal of the companion gets over-subtracted, particularly in our J-band data, compelling us to use a fake companion spectrum to estimate and compensate for the loss of flux. Finally, we report a new astrometric measurement of the position of the companion (sep = 0.817″, PA = 227.92 deg).

ExoMol: molecular line lists for exoplanet and other atmospheres

The discovery of extrasolar planets is one of the major scientific advances of the last two decades. Hundreds of planets have now been detected and astronomers are beginning to characterise their composition and physical characteristics. To do this requires a huge quantity of spectroscopic data most of which is not available from laboratory studies. The ExoMol project will offer a comprehensive solution to this problem by providing spectroscopic data on all the molecular transitions of importance in the atmospheres of exoplanets. These data will be widely applicable to other problems and will be used for studies on cool stars, brown dwarfs and circumstellar environments. This paper lays out the scientific foundations of this project and reviews previous work in this area. A mixture of first principles and empirically-tuned quantum mechanical methods will be used to compute comprehensive and very large rotation-vibration and rotation-vibration-electronic (rovibronic) line lists. Methodologies will be developed for treating larger molecules such as methane and nitric acid. ExoMol will rely on these developments and the use of state-of-the-art computing.

The Exozodiacal Dust Problem for Direct Observations of ExoEarths

Debris dust in the habitable zones of stars – otherwise known as exozodiacal dust – comes from extrasolar asteroids and comets and is thus an expected part of a planetary system. Background flux from the Solar System’s zodiacal dust and the exozodiacal dust in the target system is likely to be the largest source of astrophysical noise in direct observations of terrestrial planets in the habitable zones of nearby stars. Furthermore, dust structures like clumps, thought to be produced by dynamical interactions with exoplanets, are a possible source of confusion. In this paper, we qualitatively assess the primary impact of exozodical dust on high-contrast direct imaging at optical wavelengths, such as would be performed with a coronagraph. Then we present the sensitivity of previous, current, and near-term facilities to thermal emission from debris dust at all distances from nearby solar-type stars, as well as our current knowledge of dust levels from recent surveys. Finally, we address the other method of detecting debris dust, through high-contrast imaging in scattered light. This method is currently far less sensitive than thermal emission observations, but provides high spatial resolution for studying dust structures. This paper represents the first report of NASA’s Exoplanet Exploration Program Analysis Group (ExoPAG).

Constraining the Planetary System of Fomalhaut Using High-Resolution ALMA Observations

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring’s morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring’s morphology.

Ice Mineralogy Across and Into the Surfaces of Pluto, Triton, and Eris

We present three near-infrared spectra of Pluto taken with the IRTF and SpeX, an optical spectrum of Triton taken with the MMT and the Red Channel Spectrograph, and previously published spectra of Pluto, Triton, and Eris. We combine these observations with a two-phase Hapke model, and gain insight into the ice mineralogy on Pluto, Triton, and Eris. Specifically, we measure the methane-nitrogen mixing ratio across and into the surfaces of these icy dwarf planets. In addition, we present a laboratory experiment that demonstrates it is essential to model methane bands in spectra of icy dwarf planets with two methane phases – one highly-diluted by nitrogen and the other rich in methane. For Pluto, we find bulk, hemisphere-averaged, methane abundances of 9.1 \pm 0.5%, 7.1 \pm 0.4%, and 8.2 \pm 0.3% for sub-Earth longitudes of 10\degree, 125\degree, and 257\degree. Application of the Wilcoxon rank sum test to our measurements finds these small differences are statistically significant. For Triton, we find bulk, hemisphere-averaged, methane abundances of 5.0 \pm 0.1% and 5.3 \pm 0.4% for sub-Earth longitudes of 138\degree and 314\degree. Application of the Wilcoxon rank sum test to our measurements finds the differences are not statistically significant. For Eris, we find a bulk, hemisphere-averaged, methane abundance of 10 \pm 2%. Pluto, Triton, and Eris do not exhibit a trend in methane-nitrogen mixing ratio with depth into their surfaces over the few cm range probed by these observations. This result is contrary to the expectation (Grundy & Stansberry 2000) that since visible light penetrates deeper into a nitrogen-rich surface than the depths from which thermal emission emerges, net radiative heating at depth would drive preferential sublimation of nitrogen leading to an increase in the methane abundance with depth.

A peculiar class of debris disks from Herschel/DUNES - A steep fall off in the far infrared

Aims. We present photometric data of debris disks around HIP 103389 (HD 199260), HIP 107350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3 sigma sensitivity of a few mJy at 100 um and 160 um. In addition, we obtained Herschel/PACS photometric data at 70 um for HIP 103389. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70 um. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a very distinct range of grain sizes is implied to dominate the thermal emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented. (abridged)

Rethinking the Paleoproterozoic Great Oxidation Event: A Biological Perspective

Competing geophysical/geochemical hypotheses for how Earth’s surface became oxygenated – organic carbon burial, hydrogen escape to space, and changes in the redox state of volcanic gases – are examined and a more biologically-based hypothesis is offered in response. It is argued that organic carbon burial is of minor importance to the accumulation of oxygen in a mainly anoxic world where aerobic respiration is not globally significant. Thus, for the Great Oxidation Event (GOE) ~ 2.4 Gyr ago, an increasing flux of O2 due to its production by an expanding population of cyanobacteria is parameterized as the primary source of O2. Various factors would have constrained cyanobacterial proliferation and O2 production during most of the Archean and therefore a long delay between the appearance of cyanobacteria and oxygenation of the atmosphere is to be expected. Destruction of O2 via CH4 oxidation in the atmosphere was a major O2 sink during the Archean, and the GOE is explained to a significant extent by a large decline in the methanogen population and corresponding CH4 flux which, in turn, was caused primarily by partial oxygenation of the surface ocean. The partially oxygenated state of these waters also made it possible for an aerobic methanotroph population to become established. This further contributed to the large reduction in the CH4 flux to the atmosphere by increasing the consumption of CH4 diffusing upwards from the deeper anoxic depths of the water column as well as any CH4 still being produced in the upper layer. The reduction in the CH4 flux lowered the CH4 oxidation sink for O2 at about the same time the metamorphic and volcanic gas sinks for O2 also declined. As the O2 source increased from an expanding population of cyanobacteria – triggered by a burst of continent formation ~ 2.7-2.4 Gyr ago – the atmosphere flipped and became permanently oxygenated.

The dynamics of inner dead-zone boundaries in protoplanetary disks [Replacement]

In protoplanetary disks, the inner radial boundary between the MRI turbulent (`active’) and MRI quiescent (`dead’) zones plays an important role in models of the disk evolution and in some planet formation scenarios. In reality, this boundary is not well-defined: thermal heating from the star in a passive disk yields a transition radius close to the star (< 0.1 au), whereas if the disk is already MRI active, it can self-consistently maintain the requisite temperatures out to a transition radius of roughly 1 au. Moreover, the interface may not be static; it may be highly fluctuating or else unstable. In this paper, we study a reduced model of the dynamics of the active/dead zone interface that mimics several important aspects of a real disk system. We find that MRI-transition fronts propagate inward (a `dead front' suppressing the MRI) if they are initially at the larger transition radius, or propagate outward (an `active front' igniting the MRI) if starting from the smaller transition radius. In both cases, the front stalls at a well-defined intermediate radius, where it remains in a quasi-static equilibrium. We propose that it is this new, intermediate stalling radius that functions as the true boundary between the active and dead zones in protoplanetary disks. These dynamics are likely implicated in observations of variable accretion, such as FU Ori outbursts, as well as in those planet formation theories that require the accumulation of solid material at the dead/active interface.

The dynamics of inner dead-zone boundaries in protoplanetary disks

In protoplanetary disks, the inner radial boundary between the MRI turbulent (`active’) and MRI quiescent (`dead’) zones plays an important role in models of the disk evolution and in some planet formation scenarios. In reality, this boundary is not well-defined: thermal heating from the star in a passive disk yields a transition radius close to the star ($<0.1$ au), whereas if the disk is already MRI active, it can self-consistently maintain the requisite temperatures out to a transition radius of roughly 1 au. Moreover, the interface may not be static; it may be highly fluctuating or else unstable. In this paper, we study a reduced model of the dynamics of the active/dead zone interface that mimics several important aspects of a real disk system. We find that MRI-transition fronts propagate inward (a `dead front' suppressing the MRI) if they are initially at the larger transition radius, or propagate outward (an `active front' igniting the MRI) if starting from the smaller transition radius. In both cases, the front stalls at a well-defined intermediate radius, where it remains in a quasi-static equilibrium. We propose that it is this new, intermediate stalling radius that functions as the true boundary between the active and dead zones in protoplanetary disks. These dynamics are likely implicated in observations of variable accretion, such as FU Ori outbursts, as well as in those planet formation theories that require the accumulation of solid material at the dead/active interface.

State-resolved valence shell photoionization of Be-like ions: experiment and theory

High-resolution photoionization experiments were carried out using beams of Be-like C$^{2+}$, N$^{3+}$, and O$^{4+}$ ions with roughly equal populations of the $^1$S ground-state and the $^3$P$^o$ manifold of metastable components. The energy scales of the experiments are calibrated with uncertainties of 1 to 10 meV depending on photon energy. Resolving powers beyond 20,000 were reached allowing for the separation of contributions from the individual metastable $^3$P$^o_0$, $^3$P$^o_1$, and $^3$P$^o_2$ states. The measured data compare favourably with semi-relativistic Breit-Pauli R-matrix

Observations of the 18-cm OH lines of comet 103P/Hartley 2 at Nan\c{c}ay in support to the EPOXI and Herschel missions

The 18-cm radio lines of the OH radical were observed in comet 103P/Hartley 2 with the Nan\c{c}ay radio telescope in support to its flyby by the EPOXI mission and to observations with the Herschel Space Observatory. The OH lines were detected from 24 September to 15 December 2010. These observations are used to estimate the gas expansion velocity within the coma to 0.83 \pm 0.08 km/s in October 2010. The water production increased steeply but progressively before perihelion, and reached 1.9 \pm 0.3 X 10E28 s-1 just before the EPOXI flyby.

The large-scale energetic ion layer in the high latitude Jovian magnetosphere as revealed by Ulysses/HI-SCALE cross-field intensity-gradient measurements

Ulysses investigated the high latitude Jovian magnetosphere for a second time after Pioneer 11 mission and gave us the opportunity to search the structure and the dynamics of this giant magnetosphere above the magnetodisc. Kivelson(1976) and Kennel & Coroniti(1979) reported that Pioneer 11 observed energetic particle intensities at high latitudes at the same level with those measured in the plasma sheet and inferred that they were not consistent with the magnetodisc model. Ulysses observations supported the idea about a large-scale layer of energetic ions and electrons in the outer high latitude Jovian magnetosphere (Cowley et al.1996; Anagnostopoulos et al. 2001). This study perform a number of further tests for the existence of the large scale layer of energetic ions in the outer high latitude Jovian magnetosphere by studying appropriate cross-B field anisotropies in order to monitor the ion northward/southward intensity gradients. In particular, we examined Ulysses/HI-SCALE observations of energetic ions with large gyro-radius (0.5-1.6MeV protons and >2.5MeV heavy(Z>5) ions) in order to compare instant intensity changes with remote sensing intensity gradients. Our analysis confirms the existence of an energetic particle layer in the north hemisphere, during the inbound trajectory of Ulysses traveling at moderate latitudes, and in the south high-latitude duskside magnetosphere, during the outbound segment of the spacecraft trajectory. Our Ulysses/HI-SCALE data analysis also provides evidence for the detection of an energetic proton magnetopause boundary layer during the outbound trajectory of the spacecraft. During Ulysses flyby of Jupiter the almost permanent appearance of alternative northward and southward intensity gradients suggests that the high latitude layer appeared to be a third major area of energetic particles, which coexisted with the radiation belts and the magnetodisc.

A Resolved Census of Millimeter Emission from Taurus Multiple Star Systems

We present a high angular resolution millimeter-wave dust continuum imaging survey of circumstellar material associated with individual components of 23 multiple star systems in the Taurus-Auriga young cluster. Combined with previous measurements, these new data permit a comprehensive look at how millimeter luminosity (a tracer of disk mass) relates to the separation and mass of a stellar companion. Approximately one third (28-37%) of individual stars in multiples have detectable millimeter emission, a rate half that for single stars (~62%). There is a strong correlation between the luminosity and projected separation (a_p) of a stellar pair. Wide pairs (a_p > 300 AU) have a similar luminosity distribution as single stars, medium pairs (a_p ~ 30-300 AU) are a factor of 5 fainter, and close pairs (a_p < 30 AU) are ~ 5 times fainter yet (aside from a small population of bright circumbinary disks). In most cases, the emission is dominated by a disk around the primary (or a wide tertiary in triples), but there is no clear relationship between luminosity and stellar mass ratio. A direct comparison of resolved disk sizes with predictions from tidal truncation models yields mixed results; some disks are larger than expected given their companion separations. We suggest that the presence of a stellar companion impacts disk properties at a level comparable to the internal evolution mechanisms operating in isolated systems, with both the multiple star formation process itself and star-disk tidal interactions likely playing important roles in the evolution of disk material. From the perspective of the mass content of the disk, we expect that (giant) planet formation is inhibited around the components of close pairs or secondaries, but should be as likely as for single stars around the primaries (or wide tertiaries in hierarchical triples) in more widely-separated multiple star systems.

Circumstellar disks and planets. Science cases for next-generation optical/infrared long-baseline interferometers

We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.

Resolving HD 100546 disc in the mid-infrared: Small and asymmetric inner disc inside a bright symmetric edge of the outer disc

A region of roughly half of the Solar system scale around the star HD 100546 is largely cleared of gas and dust, in contrast to the outer disc extending to about 400 AU. However, some material is observed in the immediate vicinity of the star, called the inner disc. Studying the structure of the inner and the outer disc is a first step to establish the origin of the gap between them and possibly link it to presence of planets. We answer how the dust is distributed within and outside the gap, and constrain the disc geometry. We observe the disc with VLTI interferometer N-band instrument MIDI. At these wavelengths disc completely dominates over the stellar emission. With baseline lengths of 40m our long baseline observations (8.2m telescopes) are most sensitive to the inner few AU from the star, and we combine them with observations at shorter, 15m baselines (1.8m telescopes), to probe emission beyond the gap at up to 20AU from the star. We derive an upper limit of 0.7AU for the mid-infrared size of the inner disc, from our longest baseline data. The chromatic phases show that the N-band brightness of the inner disc is not point-symmetric. Our short baseline data place a bright symmetric ring of emission at 11AU. This is consistent with prior observations of the transition region between the gap and the outer disc, known as the disc wall. The ring inclination and position angles are constrained by our data to i=53+-8deg and PA=145+-5deg. These values are close to known estimates of the rim and disc geometry and suggest co-planarity. Micron-sized dust is distributed asymmetrically in the region from the dust sublimation radius to less than 0.7AU from HD100546 in observations from 2004 Jun to 2005 Dec. This small dusty disc is separated from the symmetric edge of the outer disc by a large, ~10AU wide gap cleared of micron-sized dust but possibly populated by planetesimals and/or planets.

Dispelling the myth of robotic efficiency: why human space exploration will tell us more about the Solar System than will robotic exploration alone

There is a widely held view in the astronomical community that unmanned robotic space vehicles are, and will always be, more efficient explorers of planetary surfaces than astronauts (e.g. Coates, 2001; Clements 2009; Rees 2011). Partly this is due to a common assumption that robotic exploration is cheaper than human exploration (although, as we shall see, this isn’t necessarily true if like is compared with like), and partly from the expectation that continued developments in technology will relentlessly increase the capability, and reduce the size and cost, of robotic missions to the point that human exploration will not be able to compete. I will argue below that the experience of human exploration during the Apollo missions, more recent field analogue studies, and trends in robotic space exploration actually all point to exactly the opposite conclusion.

Saturn's icy satellites and rings investigated by Cassini - VIMS. III. Radial compositional variability

In the last few years Cassini-VIMS, the Visible and Infared Mapping Spectrometer, returned to us a comprehensive view of the Saturn’s icy satellites and rings. After having analyzed the satellites’ spectral properties (Filacchione et al. (2007a)) and their distribution across the satellites’ hemispheres (Filacchione et al. (2010)), we proceed in this paper to investigate the radial variability of icy satellites (principal and minor) and main rings average spectral properties. This analysis is done by using 2,264 disk-integrated observations of the satellites and a 12×700 pixels-wide rings radial mosaic acquired with a spatial resolution of about 125 km/pixel. The comparative analysis of these data allows us to retrieve the amount of both water ice and red contaminant materials distributed across Saturn’s system and the typical surface regolith grain sizes. These measurements highlight very striking differences in the population here analyzed, which vary from the almost uncontaminated and water ice-rich surfaces of Enceladus and Calypso to the metal/organic-rich and red surfaces of Iapetus’ leading hemisphere and Phoebe. Rings spectra appear more red than the icy satellites in the visible range but show more intense 1.5-2.0 micron band depths. The correlations among spectral slopes, band depths, visual albedo and phase permit us to cluster the saturnian population in different spectral classes which are detected not only among the principal satellites and rings but among co-orbital minor moons as well. Finally, we have applied Hapke’s theory to retrieve the best spectral fits to Saturn’s inner regular satellites using the same methodology applied previously for Rhea data discussed in Ciarniello et al. (2011).

Accretion Outbursts in Circumplanetary Disks

We describe a model for the long term evolution of a circumplanetary disk that is fed mass from a circumstellar disk and contains regions of low turbulence (dead zones). We show that such disks can be subject to accretion driven outbursts, analogous to outbursts previously modeled in the context of circumstellar disks to explain FU Ori phenomena. Circumplanetary disks around a proto-Jupiter can undergo outbursts for infall accretion rates onto the disks in the range ~10^{-9} to 10^{-7} M_sun/yr, typical of accretion rates in the T Tauri phase. During outbursts, the accretion rate and disk luminosity increases by several orders of magnitude. Most of the planet mass growth during planetary gas accretion may occur via disk outbursts involving gas that is considerably hotter than predicted by steady state models. For low infall accretion rates less than ~10^{-10} M_sun/yr that occur in late stages of disk accretion, disk outbursts are unlikely to occur, even if dead zones are present. Such conditions are favorable for the formation of icy satellites.

Stellar wobble in triple star systems

The radial velocity method for detecting extra-solar planets relies on measuring the star’s wobble around the system’s center of mass. Since this is an indirect method, we may ask if there are other dynamical effects that can mimic such wobble. In recent articles\cite{1,2,3}, we modeled the effect of a nearby binary system on a star’s radial velocity. We showed that, if we are unaware of this nearby binary, for instance because one component is unresolved or both components are faint stars, the binary’s effect may mimic a planet. Here, we review this work, explaining in which circumstances the binary’s effect may mimic a planet and we discuss what can be done in practice in order to distinguish between these two scenarios (planet or nearby binary).

Sterile and Fertile Planetary Systems - Statistical Analysis of Multi-Planet Systems in Kepler's data

The Kepler mission has discovered a large number of planetary systems. We analyze the implications of the discovered single/multi-exoplanet systems from Kepler’s data. We test a simple model in which the intrinsic occurrence of plant is an independent process, and with equal probability around all planet producing stars. This leads to a poisson distribution for the intrinsic number of planets around each host. However, the possibility of zero/low mutual inclination is taken into account, creating a correlation between detecting different planets in a given solar system, leading to a non poisson distribution for the number of transiting planets per system. Comparing the model’s predictions with the observations made by Kepler, we find that the correlation produced by planarity is insufficient and a higher correlation is needed. This implies that either the formation of one planet in the system enhances the likelihood of other planets to form, or that some stars are considerably more fertile than others. Followup observations on Kepler planet’s hosts can help pinpoint the physical nature of this correlation.

The EVIL-MC Model for Ellipsoidal Variations of Planet-Hosting Stars and Applications to the HAT-P-7 System [Replacement]

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet’s dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

The EVIL-MC Model for Ellipsoidal Variations of Planet-Hosting Stars and Applications to the HAT-P-7 System

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet’s dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

Thermal emission from WASP-24b at 3.6 and 4.5 {\mu}m

Aims. We observe occultations of WASP-24b to measure brightness temperatures and to determine whether or not its atmosphere exhibits a thermal inversion (stratosphere). Methods. We observed occultations of WASP-24b at 3.6 and 4.5 {\mu}m using the Spitzer Space Telescope. It has been suggested that there is a correlation between stellar activity and the presence of inversions, so we analysed existing HARPS spectra in order to calculate log R’HK for WASP-24 and thus determine whether or not the star is chromospherically active. We also observed a transit of WASP-24b in the Str\”{o}mgren u and y bands, with the CAHA 2.2-m telescope. Results. We measure occultation depths of 0.159 \pm 0.013 per cent at 3.6 {\mu}m and 0.202 \pm 0.018 per cent at 4.5 {\mu}m. The corresponding planetary brightness temperatures are 1974 \pm 71 K and 1944 \pm 85 K respectively. Atmosphere models with and without a thermal inversion fit the data equally well; we are unable to constrain the presence of an inversion without additional occultation measurements in the near-IR. We find log R’HK = -4.98 \pm 0.12, indicating that WASP-24 is not a chromospherically active star. Our global analysis of new and previously-published data has refined the system parameters, and we find no evidence that the orbit of WASP-24b is non-circular. Conclusions. These results emphasise the importance of complementing Spitzer measurements with observations at shorter wavelengths to gain a full understanding of hot Jupiter atmospheres.

Disk-satellite interaction in disks with density gaps [Replacement]

Gravitational coupling between a gaseous disk and an orbiting perturber leads to angular momentum exchange between them which can result in gap opening by planets in protoplanetary disks and clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks. Understanding the co-evolution of the disk and the orbit of the perturber in these circumstances requires knowledge of the spatial distribution of the torque exerted by the latter on a highly nonuniform disk. Here we explore disk-satellite interaction in disks with gaps in linear approximation both in Fourier and in physical space, explicitly incorporating the disk non-uniformity in the fluid equations. Density gradients strongly displace the positions of Lindblad resonances in the disk (which often occur at multiple locations), and the waveforms of modes excited close to the gap edge get modified compared to the uniform disk case. The spatial distribution of the excitation torque density is found to be quite different from the existing prescriptions: most of the torque is exerted in a rather narrow region near the gap edge where Lindblad resonances accumulate, followed by an exponential fall-off with the distance from the perturber. Despite these differences, for a given gap profile the full integrated torque exerted on the disk agrees with the conventional uniform disk theory prediction at the level of ~10%. The nonlinearity of the density wave excited by the perturber is shown to decrease as the wave travels out of the gap, slowing down its nonlinear evolution and damping. Our results suggest that gap opening in protoplanetary disks and gas clearing around SMBH binaries can be more efficient than the existing theories predict. They pave the way for self-consistent calculations of the gap structure and the orbital evolution of the perturber using accurate prescription for the torque density behavior.

Disk-satellite interaction in disks with density gaps

Gravitational coupling between a gaseous disk and an orbiting perturber leads to angular momentum exchange between them which can result in gap opening by planets in protoplanetary disks and clearing of gas by binary supermassive black holes (SMBHs) embedded in accretion disks. Understanding the co-evolution of the disk and the orbit of the perturber in these circumstances requires knowledge of the spatial distribution of the torque exerted by the latter on a highly nonuniform disk. Here we explore disk-satellite interaction in disks with gaps in linear approximation both in Fourier and in physical space, explicitly incorporating the disk non-uniformity in the fluid equations. Density gradients strongly displace the positions of Lindblad resonances in the disk (which often occur at multiple locations), and the waveforms of modes excited close to the gap edge get modified compared to the uniform disk case. The spatial distribution of the excitation torque density is found to be quite different from the existing prescriptions: most of the torque is exerted in a rather narrow region near the gap edge where Lindblad resonances accumulate, followed by an exponential fall-off with the distance from the perturber. Despite these differences, for a given gap profile the full integrated torque exerted on the disk agrees with the conventional uniform disk theory prediction at the level of ~10%. The nonlinearity of the density wave excited by the perturber is shown to decrease as the wave travels out of the gap, slowing down its nonlinear evolution and damping. Our results suggest that gap opening in protoplanetary disks and gas clearing around SMBH binaries can be more efficient than the existing theories predict. They pave the way for self-consistent calculations of the gap structure and the orbital evolution of the perturber using accurate prescription for the torque density behavior.

Thermo-magneto coupling in a dipole plasma [Cross-Listing]

On a dipole plasma, we observe the generation of magnetic moment, as the movement of the levitating magnet-plasma compound, in response to electron-cyclotron heating and the increase of $\beta$ (magnetically-confined thermal energy). We formulate a thermodynamic model with interpreting heating as injection of microscopic magnetic moment; the corresponding chemical potential is the ambient magnetic field.

The Photoeccentric Effect and Proto-Hot-Jupiters I. Measuring photometric eccentricities of individual transiting planets

Exoplanet orbital eccentricities offer valuable clues about the origins and orbital evolution of planetary systems. Eccentric, Jupiter-sized planets are particularly interesting: they may link the “cold” Jupiters beyond the ice line to hot Jupiters at a fraction of an AU, where they are unlikely to have formed in situ. To date, all eccentricities of individual planets come from radial velocity measurements. Kepler has discovered hundreds of transiting Jupiters spanning a range of periods, but the faintness of the host stars precludes radial velocity follow-up of most. Here we demonstrate a Bayesian method of measuring an individual planet’s eccentricity solely from its transit light curve using prior knowledge of its host star’s density. We show that eccentric Jupiters are readily identified by their short ingress/egress/total transit durations — the “photoeccentric effect” — even with long-cadence Kepler photometry and loosely-constrained stellar parameters. A Markov Chain Monte Carlo exploration of parameter posteriors naturally marginalizes over the periapse angle and automatically accounts for the transit probability. As a demonstration, we use four published transit light curves of HD 17156 b to measure an eccentricity of e = 0.72 +0.14/-0.09, in good agreement with the discovery value e = 0.67 +/- 0.08 based on 33 radial-velocity measurements. We present two additional tests using actual Kepler data. In each case the technique proves to be a viable method of measuring exoplanet eccentricities and their confidence intervals. Finally, we argue that this method is the most efficient, effective means of identifying the extremely eccentric, proto-hot-Jupiters predicted by Socrates and collaborators.

K-shell photoionization of ground-state Li-like carbon ions [C$^{3+}$]: experiment, theory and comparison with time-reversed photorecombination

Absolute cross sections for the K-shell photoionization of ground-state Li-like carbon [C$^{3+}$(1s$^2$2s $^2$S)] ions were measured by employing the ion-photon merged-beams technique at the Advanced Light Source. The energy ranges 299.8–300.15 eV, 303.29–303.58 eV and 335.61–337.57 eV of the [1s(2s2p)$^3$P]$^2$P, [1s(2s2p)$^1$P]$^2$P and [(1s2s)$^3$S 3p]$^2$P resonances, respectively, were investigated using resolving powers of up to 6000. The autoionization linewidth of the [1s(2s2p)$^1$P]$^2$P resonance was measured to be $27 \pm 5$ meV and compares favourably with a theoretical result of 26 meV obtained from the intermediate coupling R-Matrix method. The present photoionization cross section results are compared with the outcome from photorecombination measurements by employing the principle of detailed balance.

K-shell photoionization of ground-state Li-like boron ions [B$^{2+}$]: Experiment and Theory

Absolute cross sections for the K-shell photoionization of ground-state Li-like boron [B$^{2+}$(1s$^2$2s $^2$S)] ions were measured by employing the ion-photon merged-beams technique at the Advanced Light Source synchrotron radiation facility. The energy ranges 197.5–200.5 eV, 201.9–202.1 eV of the [1s(2s\,2p)$^3$P]$^2$P${\rm ^o}$ and [1s(2s\,2p)$^1$P] $^2$P${\rm ^o}$ resonances, respectively, were investigated using resolving powers of up to 17\,600. The energy range of the experiments was extended to about 238.2 eV yielding energies of the most prominent [1s(2$\ell$\,n$\ell^{\prime}$)]$^2$P$^o$ resonances with an absolute accuracy of the order of 130 ppm. The natural linewidths of the [1s(2s\,2p)$^3$P] $^2$P${\rm ^o}$ and [1s(2s\,2p)$^1$P] $^2$P${\rm ^o}$ resonances were measured to be $4.8 \pm 0.6$ meV and $29.7 \pm 2.5$ meV, respectively, which compare favourably with theoretical results of 4.40 meV and 30.53 meV determined using an intermediate coupling R-matrix method.

Numerical Modeling of the 2009 Impact Event on Jupiter

We have investigated the 2009 July impact event on Jupiter using the ZEUS-MP 2 three-dimensional hydrodynamics code. We studied the impact itself and the following plume development. Eight impactors were considered: 0.5 km and 1 km porous (\rho = 1.760 g cm^{-3}) and non-porous (\rho = 2.700 g cm^{-3}) basalt impactors, and 0.5 km and 1 km porous (\rho = 0.600 g cm^{-3}) and non-porous \rho = 0.917 g cm^{-3}) ice impactors. The simulations consisted of these bolides colliding with Jupiter at an incident angle of \theta = 69 degrees from the vertical and with an impact velocity of v = 61.4 km s^{-1}. Our simulations show the development of relatively larger, faster plumes created after impacts involving 1 km diameter bodies. Comparing simulations of the 2009 event with simulations of the Shoemaker-Levy 9 events reveals a difference in plume development, with the higher incident angle of the 2009 impact leading to a shallower terminal depth and a smaller and slower plume. We also studied the amount of dynamical chaos present in the simulations conducted at the 2009 incident angle. Compared to the chaos of the SL9 simulations, where \theta is approximately 45 degrees, we find no significant difference in chaos at the higher 2009 incident angle.

A Possible Detection of Occultation by a Proto-planetary Clump in GM Cephei

GM Cep in the young (~4 Myr) open cluster Trumpler 37 has been known to be an abrupt variable and to have a circumstellar disk with very active accretion. Our monitoring observations in 2009-2011 revealed the star to show sporadic flare events, each with brightening of < 0.5 mag lasting for days. These brightening events, associated with a color change toward the blue, should originate from an increased accretion activity. Moreover, the star also underwent a brightness drop of ~1 mag lasting for about a month, during which the star became bluer when fainter. Such brightness drops seem to have a recurrence time scale of a year, as evidenced in our data and the photometric behavior of GM Cep over a century. Between consecutive drops, the star brightened gradually by about 1 mag and became blue at peak luminosity. We propose that the drop is caused by obscuration of the central star by an orbiting dust concentration. The UX Orionis type of activity in GM Cep therefore exemplifies the disk inhomogeneity process in transition between grain coagulation and planetesimal formation in a young circumstellar disk.

Dynamical analysis and constraints for the HD 196885 system

The HD\,196885 system is composed of a binary star and a planet orbiting the primary. The orbit of the binary is fully constrained by astrometry, but for the planet the inclination with respect to the plane of the sky and the longitude of the node are unknown. Here we perform a full analysis of the HD\,196885 system by exploring the two free parameters of the planet and choosing different sets of angular variables. We find that the most likely configurations for the planet is either nearly coplanar orbits (prograde and retrograde), or highly inclined orbits near the Lidov-Kozai equilibrium points, i = 44^{\circ} or i = 137^{\circ} . Among coplanar orbits, the retrograde ones appear to be less chaotic, while for the orbits near the Lidov-Kozai equilibria, those around \omega= 270^{\circ} are more reliable, where \omega_k is the argument of pericenter of the planet’s orbit with respect to the binary’s orbit. From the observer’s point of view (plane of the sky) stable areas are restricted to (I1, \Omega_1) \sim (65^{\circ}, 80^{\circ}), (65^{\circ},260^{\circ}), (115^{\circ},80^{\circ}), and (115^{\circ},260^{\circ}), where I1 is the inclination of the planet and \Omega_1 is the longitude of ascending node.

Tidal Venuses: Triggering a Climate Catastrophe via Tidal Heating

Traditionally stellar radiation has been the only heat source considered capable of determining global climate on long timescales. Here we show that terrestrial exoplanets orbiting low-mass stars may be tidally heated at high enough levels to induce a runaway greenhouse for a long enough duration for all the hydrogen to escape. Without hydrogen, the planet no longer has water and cannot support life. We call these planets “Tidal Venuses,” and the phenomenon a “tidal greenhouse.” Tidal effects also circularize the orbit, which decreases tidal heating. Hence, some planets may form with large eccentricity, with its accompanying large tidal heating, and lose their water, but eventually settle into nearly circular orbits in the habitable zone (HZ). However, these planets are not habitable as past tidal heating desiccated them, and hence should not be ranked highly for detailed follow-up observations aimed at detecting biosignatures. We simulate the evolution of hypothetical planetary systems in a quasi-continuous parameter distribution and find that we can constrain the history of the system by statistical arguments. Planets orbiting stars with masses <0.3 solar masses may be in danger of desiccation via tidal heating. We apply these concepts to Gl 667C c, a ~4.5 earth-mass planet orbiting a 0.3 solar mass star at 0.12 AU. We find that it probably did not lose its water via tidal heating as orbital stability is unlikely for the high eccentricities required for the tidal greenhouse. As the inner edge of the HZ is defined by the onset of a runaway or moist greenhouse powered by radiation, our results represent a fundamental revision to the HZ for non-circular orbits. In the appendices we review a) the moist and runaway greenhouses, b) stellar mass-radius and mass-luminosity relations, c) terrestrial planet mass-radius relations, and d) linear tidal theories.

Spectroscopic identification of habitable extra-solar planets

An Earth-like extra-solar planet emits light which is many orders of magnitude fainter than that of the parent star. We propose a method of identifying bio-signature spectral lines in light from known extra-solar planets based on Fourier spectroscopy in the infra-red, using an off-center part of a Fourier interferogram only. This results in superior sensitivity to narrower molecular-type spectral bands, which are expected in the planet spectrum but are absent in the parent star. We support this idea by numerical simulations which include photon and thermal noise, and show it to be feasible at a luminosity ratio of 10^-6 for a Sun-like parent star in the infra-red. We also carried out a laboratory experiment to illustrate the method. The results suggest that this method should be applicable to real planet searches.

Rossiter-McLaughlin Effect Measurements for WASP-16, WASP-25 and WASP-31

We present new measurements of the Rossiter-McLaughlin (RM) effect for three WASP planetary systems, WASP-16, WASP-25 and WASP-31, from a combined analysis of their complete sets of photometric and spectroscopic data. We find a low amplitude RM effect for WASP-16 (Teff = 5700 \pm 150K), suggesting that the star is a slow rotator and thus of an advanced age, and obtain a projected alignment angle of lambda = -4.2 degrees +11.0 -13.9. For WASP-25 (Teff = 5750\pm100K) we detect a projected spin-orbit angle of lambda = 14.6 degrees \pm6.7. WASP-31 (Teff = 6300\pm100K) is found to be well-aligned, with a projected spin-orbit angle of lambda = 2.8degrees \pm3.1. A circular orbit is consistent with the data for all three systems, in agreement with their respective discovery papers. We consider the results for these systems in the context of the ensemble of RM measurements made to date. We find that whilst WASP-16 fits the hypothesis of Winn et al. (2010) that ‘cool’ stars (Teff < 6250K) are preferentially aligned, WASP-31 has little impact on the proposed trend. We bring the total distribution of the true spin-orbit alignment angle, psi, up to date, noting that recent results have improved the agreement with the theory of Fabrycky & Tremaine (2007) at mid-range angles. We also suggest a new test for judging misalignment using the Bayesian Information Criterion, according to which WASP-25 b's orbit should be considered to be aligned.

 

You need to log in to vote

The blog owner requires users to be logged in to be able to vote for this post.

Alternatively, if you do not have an account yet you can create one here.

Powered by Vote It Up

^ Return to the top of page ^