Deformation and tidal evolution of close-in planets and satellites using a Maxwell viscoelastic rheology

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

Small particles in Pluto's environment: Effects of the solar radiation pressure

Impacts of micrometeoroids on the surfaces of the plutonian small satellites Nix and Hydra can generate dust particles. Even in this region so far from the Sun these tiny ejected particles are under the effects of the solar radiation pressure. In this work, we investigate the orbital evolution of the escaping ejecta from both the small satellites under the effects of the radiation pressure combined with the gravitational effects of Pluto, Charon, Nix and Hydra. The mass production rate of micron-sized dust particles generated by micrometeoroids hitting the satellites is obtained, and numerical simulations are performed to derive the lifetime of the ejecta. These pieces of information allow us to estimate the optical depth of a putative ring, which extends from the orbits of Nix to Hydra. The ejected particles, between the orbits of Nix and Hydra, form a wide ring of about 16 000 km. Collisions with the massive bodies and escape from the system are mainly determined by the effects of the solar radiation pressure. This is an important loss mechanism, removing 30 per cent of the initial set of 1 μm-sized particles in 1 yr. The surviving particles form a ring too faint to be detectable with the derived maximum optical depth of 4 × 10-11. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

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

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

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

Stable satellites around extrasolar giant planets

In this work, we study the stability of hypothetical satellites of extrasolar planets. Through numerical simulations of the restricted elliptic three-body problem we found the borders of the stable regions around the secondary body. From the empirical results, we derived analytical expressions of the critical semimajor axis beyond which the satellites would not remain stable. The expressions are given as a function of the eccentricities of the planet, e(P), and of the satellite, e(sat). In the case of prograde, satellites, the critical semimajor axis, in the units of Hill's radius, is given by a(E) approximate to 0.4895 (1.0000 - 1.0305e(P) - 0.2738e(sat)). In the case of retrograde satellites, it is given by a(E) approximate to 0.9309 (1.0000 - 1.0764e(P) - 0.9812e(sat)). We also computed the satellite stability region (a(E)) for a set of extrasolar planets. The results indicate that extrasolar planets in the habitable zone could harbour the Earth-like satellites.

Effects of the planetary migration on some primordial satellites of the outer planets

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

Nebular gas drag and planetary accretion with eccentric high-mass planets

Aims.We investigate the dynamics of pebbles immersed in a gas disk interacting with a planet on an eccentric orbit. The model has a prescribed gap in the disk around the location of the planetary orbit, as is expected for a giant planet with a mass in the range of 0.1-1 Jupiter masses. The pebbles with sizes in the range of 1 cm to 3 m are placed in a ring outside of the giant planet orbit at distances between 10 and 30 planetary Hill radii. The process of the accumulation of pebbles closer to the gap edge, its possible implication for the planetary accretion, and the importance of the mass and the eccentricity of the planet in this process are the motivations behind the present contribution. Methods. We used the Bulirsch-Stoer numerical algorithm, which is computationally consistent for close approaches, to integrate the Newtonian equations of the planar (2D), elliptical restricted three-body problem. The angular velocity of the gas disk was determined by the appropriate balance between the gravity, centrifugal, and pressure forces, such that it is sub-Keplerian in regions with a negative radial pressure gradient and super-Keplerian where the radial pressure gradient is positive. Results. The results show that there are no trappings in the 1:1 resonance around the L 4 and L5 Lagrangian points for very low planetary eccentricities (e2 < 0.07). The trappings in exterior resonances, in the majority of cases, are because the angular velocity of the disk is super-Keplerian in the gap disk outside of the planetary orbit and because the inward drift is stopped. Furthermore, the semi-major axis location of such trappings depends on the gas pressure profile of the gap (depth) and is a = 1.2 for a planet of 1 MJ. A planet on an eccentric orbit interacts with the pebble layer formed by these resonances. Collisions occur and become important for planetary eccentricity near the present value of Jupiter (e 2 = 0.05). The maximum rate of the collisions onto a planet of 0.1 MJ occurs when the pebble size is 37.5 cm ≤ s < 75 cm; for a planet with the mass of Jupiter, it is15 cm ≤ s < 30 cm. The accretion stops when the pebble size is less than 2 cm and the gas drag dominates the motion. © 2013 ESO.

Analysis of 25 mutual eclipses and occultations between the Galilean satellites observed from Brazil in 2009

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

The Gemini NICI planet-finding campaign: the frequency of giant planets around young B and A stars

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (≈1.5-2.5 M ☉) conducted to date and includes the planet hosts β Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58$^{+21}_{-20}$ M Jup and 55$^{+20}_{-19}$ M Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M ☉ stars can have giant planets greater than 4 M Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

The Mass of Kepler-93b and The Composition of Terrestrial Planets

Kepler-93b is a 1.478 ± 0.019 R ⊕ planet with a4.7 day period around a bright (V = 10.2), astroseismicallycharacterized host star with a mass of 0.911 ± 0.033 M⊙ and a radius of 0.919 ± 0.011 R⊙. Based on 86 radial velocity observations obtainedwith the HARPS-N spectrograph on the Telescopio Nazionale Galileo and 32archival Keck/HIRES observations, we present a precise mass estimate of4.02 ± 0.68 M ⊕. The corresponding high densityof 6.88 ± 1.18 g cm-3 is consistent with a rockycomposition of primarily iron and magnesium silicate. We compareKepler-93b to other dense planets with well-constrained parameters andfind that between 1 and 6 M ⊕, all dense planetsincluding the Earth and Venus are well-described by the same fixed ratioof iron to magnesium silicate. There are as of yet no examples of suchplanets with masses >6 M ⊕. All known planets inthis mass regime have lower densities requiring significant fractions ofvolatiles or H/He gas. We also constrain the mass and period of theouter companion in the Kepler-93 system from the long-term radialvelocity trend and archival adaptive optics images. As the sample ofdense planets with well-constrained masses and radii continues to grow,we will be able to test whether the fixed compositional model found forthe seven dense planets considered in this paper extends to the fullpopulation of 1-6 M ⊕ planets.Based on observations made with the Italian Telescopio Nazionale Galileo(TNG) operated on the island of La Palma by the Fundación GalileoGalilei of the INAF (Istituto Nazionale di Astrofisica) at the SpanishObservatorio del Roque de los Muchachos of the Instituto de Astrofisicade Canarias.

Ionization in atmospheres of brown dwarfs and extrasolar planets. V. Alfvén Ionization

Observations of continuous radio and sporadic X-ray emission from low-mass objects suggest they harbor localized plasmas in their atmospheric environments. For low-mass objects, the degree of thermal ionization is insufficient to qualify the ionized component as a plasma, posing the question: what ionization processes can efficiently produce the required plasma that is the source of the radiation? We propose Alfv´en ionization as a mechanism for producing localized pockets of ionized gas in the atmosphere, having sufficient degrees of ionization ( 10−7) that they constitute plasmas. We outline the criteria required for Alfv´en ionization and demonstrate its applicability in the atmospheres of low-mass objects such as giant gas planets, brown dwarfs, and M dwarfs with both solar and sub-solar metallicities. We find that Alfv´en ionization is most efficient at mid to low atmospheric pressures where a seed plasma is easier to magnetize and the pressure gradients needed to drive the required neutral flows are the smallest. For the model atmospheres considered, our results show that degrees of ionization of 10−6–1 can be obtained as a result of Alfv´en ionization. Observable consequences include continuum bremsstrahlung emission, superimposed with spectral lines from the plasma ion species (e.g., He, Mg, H2, or CO lines). Forbidden lines are also expected from the metastable population. The presence of an atmospheric plasma opens the door to a multitude of plasma and chemical processes not yet considered in current atmospheric models. The occurrence of Alfv´en ionization may also be applicable to other astrophysical environments such as protoplanetary disks.

Close-in planets around giant stars: lack of hot-Jupiters and prevalence of multiplanetary systems

Extrasolar planets abound in almost any possible configuration. However, until five years ago, there was a lack of planets orbiting closer than 0.5 au to giant or subgiant stars. Since then, recent detections have started to populated this regime by confirming 13 planetary systems. We discuss the properties of these systems in terms of their formation and evolution off the main sequence. Interestingly, we find that 70.0 ± 6.6% of the planets in this regime are inner components of multiplanetary systems. This value is 4.2σ higher than for main-sequence hosts, which we find to be 42.4 ± 0.1%. The properties of the known planets seem to indicate that the closest-in planets (a< 0.06 au) to main-sequence stars are massive (i.e., hot Jupiters) and isolated and that they are subsequently engulfed by their host as it evolves to the red giant branch, leaving only the predominant population of multiplanetary systems in orbits 0.06 and evolution of hot Jupiters.

Transiting exoplanets from the CoRoT space mission XI. CoRoT-8b: a hot and dense sub-Saturn around a K1 dwarf

Aims. We report the discovery of CoRoT-8b, a dense small Saturn-class exoplanet that orbits a K1 dwarf in 6.2 days, and we derive its orbital parameters, mass, and radius. Methods. We analyzed two complementary data sets: the photometric transit curve of CoRoT-8b as measured by CoRoT and the radial velocity curve of CoRoT-8 as measured by the HARPS spectrometer**. Results. We find that CoRoT-8b is on a circular orbit with a semi-major axis of 0.063 +/- 0.001 AU. It has a radius of 0.57 +/- 0.02 R(J), a mass of 0.22 +/- 0.03 M(J), and therefore a mean density of 1.6 +/- 0.1 g cm(-3). Conclusions. With 67% of the size of Saturn and 72% of its mass, CoRoT-8b has a density comparable to that of Neptune (1.76 g cm(-3)). We estimate its content in heavy elements to be 47-63 M(circle plus), and the mass of its hydrogen-helium envelope to be 7-23 M(circle plus). At 0.063 AU, the thermal loss of hydrogen of CoRoT-8b should be no more than similar to 0.1% over an assumed integrated lifetime of 3 Ga.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

Dynamics of two planets in co-orbital motion

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities; we also assume that both planets share the same orbital plane. Initially, we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analysed in more detail using a semi-analytical model. Apart from the well-known quasi-satellite orbits and the classical equilibrium Lagrangian points L(4) and L(5), we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at (delta lambda, delta pi) = (+/- 60 degrees, -/+ 120 degrees), where delta lambda is the difference in mean longitudes and delta pi is the difference in longitudes of pericentre. The position of these anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities and are found for eccentricities as high as similar to 0.7. Finally, we also applied a slow mass variation to one of the planets and analysed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.