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)

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.

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)

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.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS: THE MAGNETIC FIELD GEOMETRY AND THE INFLUENCE OF THE WIND ON GIANT PLANET MIGRATION

By means of numerical simulations, we investigate magnetized stellar winds of pre-main-sequence stars. In particular, we analyze under which circumstances these stars will present elongated magnetic features (e.g., helmet streamers, slingshot prominences, etc). We focus on weak-lined T Tauri stars, as the presence of the tenuous accretion disk is not expected to have strong influence on the structure of the stellar wind. We show that the plasma-beta parameter (the ratio of thermal to magnetic energy densities) is a decisive factor in defining the magnetic configuration of the stellar wind. Using initial parameters within the observed range for these stars, we show that the coronal magnetic field configuration can vary between a dipole-like configuration and a configuration with strong collimated polar lines and closed streamers at the equator (multicomponent configuration for the magnetic field). We show that elongated magnetic features will only be present if the plasma-beta parameter at the coronal base is beta(0) << 1. Using our self-consistent three-dimensional magnetohydrodynamics model, we estimate for these stellar winds the timescale of planet migration due to drag forces exerted by the stellar wind on a hot-Jupiter. In contrast to the findings of Lovelace et al., who estimated such timescales using the Weber and Davis model, our model suggests that the stellar wind of these multicomponent coronae are not expected to have significant influence on hot-Jupiters migration. Further simulations are necessary to investigate this result under more intense surface magnetic field strengths (similar to 2-3 kG) and higher coronal base densities, as well as in a tilted stellar magnetosphere.

On the orbits of the outer satellites of Jupiter

In this work we study the basic aspects concerning the stability of the outer satellites of Jupiter. Including the effects of the four giant planets and the Sun we study a large grid of initial conditions. Some important regions where satellites cannot survive are found. Basically these regions are due to Kozai and other resonances. We give an analytical explanation for the libration of the pericenters (ω) over bar - (ω) over bar (J). Another different center is also found. The period and amplitude of these librations are quite sensitive to initial conditions, so that precise observational data are needed for Pasiphae and Sinope. The effect of Jupiter's mass variation is briefly presented. This effect can be responsible for satellite capture and also for locking (ω) over bar - (ω) over bar (J) in temporary libration.

Time analysis for temporary gravitational capture: Satellites of Uranus

We study the problem of gravitational capture in the framework of the Sun-Uranus-particle system. Part of the space of initial conditions is systematically explored, and the duration of temporary gravitational capture is measured. The location and size of different capture-time regions are given in terms of diagrams of initial semimajor axis versus eccentricity. The other initial orbital elements - inclination (i), longitude of the node (Ω), argument of pericenter (ω), and time of pericenter passage (τ) - are first taken to be zero. Then we investigate the cases with ω = 90°, 180°, and 270°. We also present a sample of results for Ω = 90°, considering the cases i = 60°, 120°, 150°, and 180°. Special attention is given to the influence of the initial orbital inclination, taking orbits initially in opposition at pericenter. In this case, the initial inclination is varied from 0° to 180° in steps of 10°. The success of the final stage of the capture problem, which involves the transformation of temporary captures into permanent ones, is highly dependent on the initial conditions associated with the longest capture times. The largest regions of the initial-conditions space with the longest capture times occur at inclinations of 60°-70° and 160°. The regions of possible stability as a function of initial inclination are also delimited. These regions include not only a known set of retrograde orbits, but also a new sort of prograde orbit with inclinations greater than zero.

ORBITAL PERTURBATIONS OF THE GALILEAN SATELLITES DURING PLANETARY ENCOUNTERS

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

Finding exoplanets orbiting young active stars - I. Technique

Stellar activity, such as starspots, can induce radial velocity (RV) variations that can mask or even mimic the RV signature of orbiting exoplanets. For this reason RV exoplanet surveys have been unsuccessful when searching for planets around young, active stars and are therefore failing to explore an important regime which can help to reveal how planets form and migrate. This paper describes a new technique to remove spot signatures from the stellar line-profiles of moderately rotating, active stars (v sin i ranging from 10 to 50 km s(-1)). By doing so it allows planetary RV signals to be uncovered. We used simulated models of a G5V type star with differing dark spots on its surface along with archive data of the known active star HD 49933 to validate our method. The results showed that starspots could be effectively cleaned from the line-profiles so that the stellar RV jitter was reduced by more than 80 per cent. Applying this procedure to the same models and HD 49933 data, but with fake planets injected, enabled the effective removal of starspots so that Jupiter mass planets on short orbital periods were successfully recovered. These results show that this approach can be useful in the search for hot-Jupiter planets that orbit around young, active stars with a v sin i of similar to 10-50 km/s.

On the dynamics of some resonances of Phobos in the future

A probable capture of Phobos into an interesting resonance was presented in our previous work. With a simple model, considering Mars in a Keplerian and circular orbit, it was shown that once captured in the resonance, the inclination of the satellite reaches very high values. Here, the integrations are extended to much longer times and escape situations are analyzed. These escapes are due to the interaction of new additional resonances, which appear as the inclination starts to increase reaching some specific values. Compared to classical capture in mean motion resonances, we see some interesting differences in this problem. We also include the effect of Mars' eccentricity in the process of the capture. The role played by this eccentricity becomes important, particularly when Phobos encounters a double resonance at a approximate to 2.619R(M). Planetary perturbations acting on Mars and variation of its equator are also included. In general, some possible scenarios of the future of Phobos are presented.

The effect of Jupiter's mass growth on satellite capture - Retrograde case

Gravitational capture can be used to explain the existence of the irregular satellites of giants planets. However, it is only the first step since the gravitational capture is temporary. Therefore, some kind of non-conservative effect is necessary to to turn the temporary capture into a permanent one. In the present work we study the effects of Jupiter mass growth for the permanent capture of retrograde satellites. An analysis of the zero velocity curves at the Lagrangian point L-1 indicates that mass accretion provides an increase of the confinement region ( delimited by the zero velocity curve, where particles cannot escape from the planet) favoring permanent captures. Adopting the restricted three-body problem, Sun-Jupiter-Particle, we performed numerical simulations backward in time considering the decrease of M-4. We considered initial conditions of the particles to be retrograde, at pericenter, in the region 100 R-4 less than or equal to a less than or equal to 400 R-4 and 0 less than or equal to e < 0.5. The results give Jupiter's mass at the moment when the particle escapes from the planet. Such values are an indication of the necessary conditions that could provide capture. An analysis of these results shows that retrograde satellites would be captured as soon as they get inside the Hills' radius and after that they keep migrating toward the planet while it is growing. For the region where the orbits of the four old retrograde satellites of Jupiter ( Ananke, Carme, Pasiphae and Sinope) are located we found that such satellites could have been permanently captured when Jupiter had between 62% and 93% of its present mass.

Time analysis for temporary gravitational capture - Stable orbits

In a previous work, Vieira Neto & Winter (2001) numerically explored the capture times of particles as temporary satellites of Uranus. The study was made in the framework of the spatial, circular, restricted three-body problem. Regions of the initial condition space whose trajectories are apparently stable were determined. The criterion adopted was that the trajectories do not escape from the planet during an integration of 10(5) years. These regions occur for a wide range of orbital initial inclinations (i). In the present work it is studied the reason for the existence of such stable regions. The stability of the planar retrograde trajectories is due to a family of simple periodic orbits and the associated quasi-periodic orbits that oscillate around them. These planar stable orbits had already been studied (Henon 1970; Huang & Innanen 1983). Their results are reviewed using Poincare surface of sections. The stable non-planar retrograde trajectories, 110 degrees less than or equal to i < 180, are found to be tridimensional quasi-periodic orbits around the same family of periodic orbits found for the planar case (i = 180 degrees). It was not found any periodic orbit out of the plane associated to such quasi-periodic orbits. The largest region of stable prograde trajectories occurs at i = 60 degrees. Trajectories in such region are found to behave as quasi-periodic orbits evolving similarly to the stable retrograde trajectories that occurs at i = 120 degrees.

Moonlets wandering on a leash-ring

Since the Voyager flybys, embedded moonlets have been proposed to explain some of the surprising structures observed in Saturn's narrow F ring. Experiments conducted with the Cassini spacecraft support this suggestion. Images of the F ring show bright compact spots, and seven occultations of stars by the F ring, monitored by ultraviolet and infrared experiments, revealed nine events of high optical depth. These results point to a large number of such objects, but it is not clear whether they are solid moonlets or rather loose particle aggregates. Subsequent images suggested an irregular motion of these objects so that a determination of their orbits consistent with the F ring failed. Some of these features seem to cross the whole ring. Here we show that these observations are explained by chaos in the F ring driven mainly by the 'shepherd' moons Prometheus and Pandora. It is characterized by a rather short Lyapunov time of about a few hundred orbital periods. Despite this chaotic diffusion, more than 93 per cent of the F-ring bodies remain confined within the F ring because of the shepherding, but also because of a weak radial mobility contrasted by an effective longitudinal diffusion. This chaotic stirring of all bodies involved prevents the formation of 'propellers' typical of moonlets, but their frequent ring crossings explain the multiple radial 'streaks' seen in the F ring. The related 'thermal' motion causes more frequent collisions between all bodies which steadily replenish F-ring dust and allow for ongoing fragmentation and re-accretion processes (ring recycling).