The New Horizons space probe led the first close flyby of one of the most primordial and distant objects left over from the formation of the solar system, the contact binary Kuiper Belt object (486958) Arrokoth, which is composed of two progenitors, the lobes nicknamed Ultima and Thule. In the current work, we investigated Arrokoth's surface in detail to identify the location of equilibrium points and also explore each lobe's individual dynamic features. We assume Arrokoth's irregular shape as a homogeneous polyhedra contact binary. We numerically explore its dynamic characteristics by computing its irregular binary geopotential to study its quantities, such as geometric height, oblateness, ellipticity, and zero-power curves. The stability of Arrokoth Hill was also explored through zero-velocity curves. Arrokoth's external equilibrium points have no radial symmetry due to its highly irregular shape. We identified even equilibrium points concerning its shape and spin rate: i.e., four unstable external equilibrium points and three inner equilibrium points, where two points are linearly stable, with an unstable central point that has a slight offset from its centroid. Moreover, the large and small lobes each have five equilibrium points with different topological structures from those found in Arrokoth. Our results also indicate that the equatorial region of Arrokoth's lobes is an unstable area due to the high rotation period, while its polar locations are stable resting sites for surface particles. Finally, the zero-power curves indicate the locations around Arrokoth where massless particles experience enhancing and receding orbital energy.

The (153591) 2001 SN263 asteroid system, target of the first Brazilian interplanetary space mission, is one of the known three triple systems within the population of NEAs. One of the mission objectives is to collect data about the formation of this system. The analysis of these data will help in the investigation of the physical and dynamical structures of the components (Alpha, Beta and Gamma) of this system, in order to find vestiges related to its origin. In this work, we assume the irregular shape of the 2001 SN263 system components as uniform density polyhedra and computationally investigate the gravitational field generated by these bodies. The goal is to explore the dynamical characteristics of the surface and environment around each component. Then, taking into account the rotational speed, we analyze their topographic features through the quantities geometric altitude, tilt, geopotential, slope, surface accelerations, among others. Additionally, the investigation of the environment around the bodies made it possible to construct zero-velocity curves, which delimit the location of equilibrium points. The Alpha component has a peculiar number of 12 equilibrium points, all of them located very close to its surface. In the cases of Beta and Gamma, we found four equilibrium points not so close to their surfaces. Then, performing numerical experiments around their equilibrium points, we identified the location and size of just one stable region, which is associated with an equilibrium point around Beta. Finally, we integrated a spherical cloud of particles around Alpha and identified the location on the surface of Alpha were the particles have fallen.

In this work is investigated the possibility of close-binary star systems having Earth-size planets within their habitable zones. First, we selected all known close-binary systems with confirmed planets (totaling 22 systems) to calculate the boundaries of their respective habitable zones (HZ). However, only eight systems had all the data necessary for the computation of the HZ. Then, we numerically explored the stability within the habitable zones for each one of the eight systems using test particles. From the results, we selected five systems that have stable regions inside the habitable zones (HZ), namely Kepler-34, 35, 38, 413 and 453. For these five cases of systems with stable regions in the HZ, we perform a series of numerical simulations for planet formation considering disks composed of planetary embryos and planetesimals, with two distinct density profiles, in addition to the stars and host planets of each system. We found that in the case of Kepler-34 and 453 systems no Earth-size planet is formed within the habitable zones. Although planets with Earth-like masses were formed in the Kepler-453, but they were outside the HZ. In contrast, for Kepler-35 and 38 systems, the results showed that potentially habitable planets are formed in all simulations. In the case of the Kepler-413 system, in just one simulation a terrestrial planet was formed within the habitable zone.

Radar observations show that (16) Psyche is one of the largest and most massive asteroids of the M-class located in the main belt, with a diameter of approximately 230 km. This fact makes Psyche a unique object since observations indicated an iron-nickel composition. It is believed that this body may be what was left of a metal core of an early planet that would have been fragmented over millions of years due to violent collisions. In this work we study a variety of dynamical aspects related to the surface, as well as, the environment around this asteroid. We use computational tools to explore the gravitational field generated by this body, assuming constant values for its density and rotation period. We then determine a set of physical and dynamical characteristics over its entire surface. The results include the geometric altitude, geopotential altitude, tilt, slope, among others. We also explore the neighborhood around the asteroid (16) Psyche, so that the location and linear stability of the equilibrium points were found. We found four external equilibrium points, two of them linearly stable. We confirmed the stability of these points by performing numerical simulations of massless particles around the asteroid, which also showed an asymmetry in the size of the stable regions. In addition, we integrate a cloud of particles in the vicinity of (16) Psyche in order to verify in which regions of its surface the particles are most likely to collide.

Progress in astrometry and orbital modelling of planetary moons in the last decade enabled better determinations of their orbits. These studies need accurate positions spread over extended periods. We present the results of the 2014-2015 Brazilian campaign for 40 mutual events from 47 observed light curves by the Galilean satellites plus one eclipse of Amalthea by Ganymede. We also reanalysed and updated results for 25 mutual events observed in the 2009 campaign. All telescopes were equipped with narrow-band filters centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The albedos' ratio was determined using images before and after each event. We simulated images of moons, umbra, and penumbra in the sky plane, and integrated their fluxes to compute albedos, simulate light curves and fit them to the observed ones using a chi-square fitting procedure. For that, we used the complete version of the Oren-Nayer reflectance model. The relative satellite positions mean uncertainty was 11.2 mas ($\sim$35 km) and 10.1 mas ($\sim$31 km) for the 2014-2015 and 2009 campaigns respectively. The simulated and observed \textsc{ascii} light curve files are freely available in electronic form at the \textit{Natural Satellites DataBase} (NSDB). The 40/25 mutual events from our 2014-2015/2009 campaigns represent a significant contribution of 17%/15% in comparison with the PHEMU campaigns lead by the IMCCE. Besides that, our result for the eclipse of Amalthea is only the 4$^{th}$ such measurement ever published after the three ones observed by the 2014-2015 international PHEMU campaign. Our results are suitable for new orbital/ephemeris determinations for the Galilean moons and Amalthea.

Context. Bright stellar positions are now known with an uncertainty below 1 mas thanks to Gaia DR2. Between 2019-2020, the Galactic plane will be the background of Jupiter. The dense stellar background will lead to an increase in the number of occultations, while the Gaia DR2 catalogue will reduce the prediction uncertainties for the shadow path. Aims. We observed a stellar occultation by the Galilean moon Europa (J2) and propose a campaign for observing stellar occultations for all Galilean moons. Methods. During a predicted period of time, we measured the light flux of the occulted star and the object to determine the time when the flux dropped with respect to one or more reference stars, and the time that it rose again for each observational station. The chords obtained from these observations allowed us to determine apparent sizes, oblatness, and positions with kilometre accuracy. Results. We present results obtained from the first stellar occultation by the Galilean moon Europa observed on 2017 March 31. The apparent fitted ellipse presents an equivalent radius of 1561.2 $\pm$ 3.6 km and oblatenesses 0.0010 $\pm$ 0.0028. A very precise Europa position was determined with an uncertainty of 0.8 mas. We also present prospects for a campaign to observe the future events that will occur between 2019 and 2021 for all Galilean moons. Conclusions. Stellar occultation is a suitable technique for obtaining physical parameters and highly accurate positions of bright satellites close to their primary. A number of successful events can render the 3D shapes of the Galilean moons with high accuracy. We encourage the observational community (amateurs included) to observe the future predicted events.

PDS 110 is a young disk-hosting star in the Orion OB1A association. Two dimming events of similar depth and duration were seen in 2008 (WASP) and 2011 (KELT), consistent with an object in a closed periodic orbit. In this paper we present data from a ground-based observing campaign designed to measure the star both photometrically and spectroscopically during the time of predicted eclipse in September 2017. Despite high-quality photometry, the predicted eclipse did not occur, although coherent structure is present suggesting variable amounts of stellar flux or dust obscuration. We also searched for RV oscillations caused by any hypothetical companion and can rule out close binaries to 0.1 $M_\odot$. A search of Sonneberg plate archive data also enabled us to extend the photometric baseline of this star back more than 50 years, and similarly does not re-detect any deep eclipses. Taken together, they suggest that the eclipses seen in WASP and KELT photometry were due to aperiodic events. It would seem that PDS 110 undergoes stochastic dimmings that are shallower and shorter-duration than those of UX Ori variables, but may have a similar mechanism.

The technique of mutual approximations accurately gives the central instant at the maximum apparent approximation of two moving natural satellites in the sky plane. This can be used in ephemeris fitting to infer the relative positions between satellites with high precision. Only the mutual phenomena -- occultations and eclipses -- may achieve better results. However, mutual phenomena only occur every six years in the case of Jupiter. Mutual approximations do not have this restriction and can be observed at any time along the year as long as the satellites are visible. In this work, we present 104 central instants determined from the observations of 66 mutual approximations between the Galilean moons carried out at different sites in Brazil and France during the period 2016--2018. For 28 events we have at least two independent observations. All telescopes were equipped with a narrow-band filter centred at 889 nm with a width of 15 nm to eliminate the scattered light from Jupiter. The telescope apertures ranged between 25--120 cm. For comparison, the precision of the positions obtained with classical CCD astrometry is about 100 mas, for mutual phenomena it can achieve 10 mas or less and the average internal precision obtained with mutual approximations was 11.3 mas. This new kind of simple, yet accurate observations can significantly improve the orbits and ephemeris of Galilean satellites and thus be very useful for the planning of future space missions aiming at the Jovian system.

In 2014, the discovery of two well-defined rings around the Centaur (10199) Chariklo were announced. This was the first time that such structures were found around a small body. In 2015, it was proposed that the Centaur (2060) Chiron may also have a ring. In a previous study, we analyzed how close encounters with giant planets would affect the rings of Chariklo. The most likely result is the survival of the rings. In the present work, we broaden our analysis to (2060) Chiron. In addition to Chariklo, Chiron is currently the only known Centaur with a presumed ring. By applying the same method as \cite{araujo2016}, we performed numerical integrations of a system composed of 729 clones of Chiron, the Sun, and the giant planets. The number of close encounters that disrupted the ring of Chiron during one half-life of the study period was computed. This number was then compared to the number of close encounters for Chariklo. We found that the probability of Chiron losing its ring due to close encounters with the giant planets is about six times higher than that for Chariklo. Our analysis showed that, unlike Chariklo, Chiron is more likely to remain in an orbit with a relatively low inclination and high eccentricity. Thus, we found that the bodies in Chiron-like orbits are less likely to retain rings than those in Chariklo-like orbits. Overall, for observational purposes, we conclude that the bigger bodies in orbits with high inclinations and low eccentricities should be prioritized.

The large increase in exoplanet discoveries in the last two decades showed a variety of systems whose stability is not clear. In this work we chose the $\upsilon$ Andromedae system as the basis of our studies in dynamical stability. This system has a range of possible masses, as a result of detection by radial velocity method, so we adopted a range of masses for the planets $c$ and $d$ and applied the secular theory. We also performed a numerical integration of the 3-body problem for the system over a time span of 30 thousand years. The results exposed similarities between the secular perturbation theory and the numerical integration, as well as the limits where the secular theory did not present good results. The analysis of the results provided hints for the maximum values of masses and eccentricities for stable planetary systems similar to $\upsilon$ Andromedae.

It is undeniable that life as we know it depends on liquid water. It is difficult to imagine any biochemical machinery that does not require water. On Earth, life adapts to the most diverse environments and, once established, it is very resilient. Considering that water is a common compound in the Universe, it seems possible (maybe even likely) that one day we will find life elsewhere in the universe. In this study, we review the main aspects of water as an essential compound for life: when it appeared since the Big Bang, and where it spread throughout the diverse cosmic sites. Then, we describe the strong relation between water and life, as we know it.

The Cassini spacecraft found a new and unique ring that shares the trajectory of Janus and Epimetheus, co-orbital satellites of Saturn. Performing image analysis, we found this to be a continuous ring. Its width is between 30% and 50% larger than previously announced. We also verified that the ring behaves like a firefly. It can only be seen from time to time, when Cassini, the ring and the Sun are arranged in a particular geometric configuration, in very high phase angles. Otherwise, it remains "in the dark", not visible to Cassini's cameras. Through numerical simulations, we found a very short lifetime for the ring particles, less than a couple of decades. Consequently, the ring needs to be constantly replenished. Using a model of particles production due to micrometeorites impacts on the surfaces of Janus and Epimetheus, we reproduce the ring, explaining its existence and the "firefly" behavior.

In general, small bodies of the solar system, e.g., asteroids and comets, have a very irregular shape. This feature affects significantly the gravitational potential around these irregular bodies, which hinders dynamical studies. The Poincar\'e surface of sec- tion technique is often used to look for stable and chaotic regions in two-dimensional dynamic cases. In this work, we show that this tool can be useful for exploring the surroundings of irregular bodies such as the asteroid 4179 Toutatis. Considering a rotating system with a particle, under the effect of the gravitational field computed three-dimensionally, we define a plane in the phase space to build the Poincar\'e surface of sections. Despite the extra dimension, the sections created allow us to find trajec- tories and classify their stabilities. Thus, we have also been able to map stable and chaotic regions, as well as to find correlations between those regions and the contri- bution of the third dimension of the system to the trajectory dynamics as well. As examples, we show details of periodic(resonant or not) and quasi-periodic trajectories.

It is well accepted that knowing the composition and the orbital evolution of asteroids may help us to understand the process of formation of the Solar System. It is also known that asteroids can represent a threat to our planet. Such important role made space missions to asteroids a very popular topic in the current astrodynamics and astronomy studies. By taking into account the increasingly interest in space missions to asteroids, especially to multiple systems, we present a study aimed to characterize the stable and unstable regions around the triple system of asteroids (45) Eugenia. The goal is to characterize unstable and stable regions of this system and compare with the system 2001 SN263 - the target of the ASTER mission. Besides, Prado (2014) used a new concept for mapping orbits considering the disturbance received by the spacecraft from all the perturbing forces individually. This method was also applied to (45) Eugenia. We present the stable and unstable regions for particles with relative inclination between 0 and 180 degrees. We found that (45) Eugenia presents larger stable regions for both, prograde and retrograde cases. This is mainly because the satellites of this system are small when compared to the primary body, and because they are not so close to each other. We also present a comparison between those two triple systems, and a discussion on how these results may guide us in the planning of future missions.

The orbital structure of the asteroid belt holds a record of the Solar System's dynamical history. The current belt only contains ${\rm \sim 10^{-3}}$ Earth masses yet the asteroids' orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter's orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth's orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter and Saturn's early orbits were chaotic. Stochastic variations in the gas giants' orbits cause resonances to continually jump across the main belt and excite the asteroids' orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants' current orbits were achieved later, during an instability in the outer Solar System. Although it is well known that the present-day Solar System exhibits chaotic behavior, our results suggest that the early Solar System may also have been chaotic.

The Centaur population is composed by minor bodies wandering between the giant planets and that frequently perform close gravitational encounters with these planets, which leads to a chaotic orbital evolution. Recently, the discovery of two well-defined narrow rings was announced around the Centaur 10199 Chariklo. The rings are assumed to be in the equatorial plane of Chariklo and to have circular orbits. The existence a well-defined system of rings around a body in such perturbed orbital region poses an interesting new problem. Are the rings of Chariklo stable when perturbed by close gravitational encounters with the giant planets? Our approach to address this question consisted of forward and backward numerical simulations of 729 clones of Chariklo, with similar initial orbits, for a period of 100 Myrs. We found, on average, that each clone suffers along its lifetime more than 150 close encounters with the giant planets within one Hill radius of the planet in question. We identified some extreme close encounters able to significantly disrupt or to disturb the rings of Chariklo. About 3% of the clones lose the rings and about 4% of the clones have the ring significantly disturbed. Therefore, our results show that in most of the cases (more than 90%) the close encounters with the giant planets do not affect the stability of the rings in Chariklo-like systems. Thus, if there is an efficient mechanism that creates the rings, then these structures may be common among these kinds of Centaurs.

We present the results of an extensive study of the final stage of terrestrial planet formation in disks with different surface density profiles and for different orbits of Jupiter and Saturn. We carried out simulations for disk densities proportional to r^-0.5, r^-1, and r^-1.5, and also for partially depleted disks as in the recent model of Mars formation by Izidoro et al (2014). The purpose of our study is to determine how the final assembly of planets and their physical properties are affected by the total mass of the disk and its radial profile. Because of the important roles of secular resonances in orbits and properties of the final planets, we studied the effects of these resonances as well. We have divided this study into two parts. In Part 1, we are interested in examining the effects of secular resonances on the formation of Mars and orbital stability of terrestrial planets. In Part 2, our goal is to determine trends that may exist between the disk surface density profile and the final properties of terrestrial planets. In the context of the depleted disk model, results show that the nu_5 resonance does not have a significant effect on the final orbits of terrestrial planets. However, nu_6 and nu_16 resonances play important roles in clearing their affected areas ensuring that no additional mass will be scattered into the accretion zone of Mars so that it can maintain its mass and orbital stability. In Part 2, our results indicate that despite some small correlations, in general, no trend seems to exist between the disk surface density profile and the mean number of the final planets, their masses, time of formation, and distances to the central star. We present the results of our simulations and discuss their implications for the formation of Mars and other terrestrial planets, as well as the physical properties of these objects such as their masses and water contents.

The Erigone family is a C-type group in the inner main belt. Its age has been estimated by several researchers to be less then 300 My, so it is a relatively young cluster. Yarko-YORP Monte Carlo methods to study the chronology of the Erigone family confirm results obtained by other groups. The Erigone family, however, is also characterized by its interaction with the $z_2$ secular resonance. While less than 15% of its members are currently in librating states of this resonance, the number of objects, members of the dynamical group, in resonant states is high enough to allow to use the study of dynamics inside the $z_2$ resonance to set constraints on the family age. Like the ${\nu}_{6}$ and $z_1$ secular resonances, the $z_2$ resonance is characterized by one stable equilibrium point at $\sigma = 180^{\circ}$ in the $z_2$ resonance plane $(\sigma, \frac{d\sigma}{dt})$, where $\sigma$ is the resonant angle of the $z_2$ resonance. Diffusion in this plane occurs on timescales of $\simeq 12$ My, which sets a lower limit on the Erigone family age. Finally, the minimum time needed to reach a steady-state population of $z_2$ librators is about 90 My, which allows to impose another, independent constraint on the group age.

Context: Co-orbital systems are bodies that share the same mean orbit. They can be divided into different families according to the relative mass of the co-orbital partners and the particularities of their movement. Janus and Epimetheus are unique in that they are the only known co-orbital pair of comparable masses and thus the only known system in mutual horseshoe orbit. Aims: We aim to establish whether the Janus-Epimetheus system might have formed by disruption of an object in the current orbit of Epimetheus. Methods: We assumed that four large main fragments were formed and neglected smaller fragments. We used numerical integration of the full N-body problem to study the evolution of different fragment arrangements. Collisions were assumed to result in perfectly inelastic merging of bodies. We statistically analysed the outcome of these simulations to infer whether co-orbital systems might have formed from the chosen initial conditions. Results: Depending on the range of initial conditions, up to 9% of the simulations evolve into co-orbital systems. Initial velocities around the escape velocity of Janus yield the highest formation probability. Analysis of the evolution shows that all co-orbital systems are produced via secondary collisions. The velocity of these collisions needs to be low enough that the fragments can merge and not be destroyed. Generally, collisions are found to be faster than an approximate cut-off velocity threshold. However, given a sufficiently low initial velocity, up to 15% of collisions is expected to result in merging. Hence, the results of this study show that the considered formation scenario is viable.

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disk between 1 and 3 AU. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 AU and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas-drag, as has been proposed to explain the formation of close-in super Earths. Here we show that the asteroid belt's orbital excitation provides a crucial constraint against this scenario for the Solar System. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from disks of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Disks with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogs significantly more massive than the real planet. In contrast, a disk with a surface density gradient proportional to $r^{-5.5}$ reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disk profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

The NEA 2001 SN263 is the target of the ASTER MISSION - First Brazilian Deep Space Mission. Araujo et al. (2012), characterized the stable regions around the components of the triple system for the planar and prograde cases. Knowing that the retrograde orbits are expected to be more stable, here we present a complementary study. We now considered particles orbiting the components of the system, in the internal and external regions, with relative inclinations between $90^{\circ}< I \leqslant180^{\circ}$, i.e., particles with retrograde orbits. Our goal is to characterize the stable regions of the system for retrograde orbits, and then detach a preferred region to place the space probe. For a space mission, the most interesting regions would be those that are unstable for the prograde cases, but stable for the retrograde cases. Such configuration provide a stable region to place the mission probe with a relative retrograde orbit, and, at the same time, guarantees a region free of debris since they are expected to have prograde orbits. We found that in fact the internal and external stable regions significantly increase when compared to the prograde case. For particles with $e=0$ and $I=180^{\circ}$, we found that nearly the whole region around Alpha and Beta remain stable. We then identified three internal regions and one external region that are very interesting to place the space probe. We present the stable regions found for the retrograde case and a discussion on those preferred regions. We also discuss the effects of resonances of the particles with Beta and Gamma, and the role of the Kozai mechanism in this scenario. These results help us understand and characterize the stability of the triple system 2001 SN263 when retrograde orbits are considered, and provide important parameters to the design of the ASTER mission.

Models of terrestrial planet formation for our solar system have been successful in producing planets with masses and orbits similar to those of Venus and Earth. However, these models have generally failed to produce Mars-sized objects around 1.5 AU. The body that is usually formed around Mars' semimajor axis is, in general, much more massive than Mars. Only when Jupiter and Saturn are assumed to have initially very eccentric orbits (e $\sim$ 0.1), which seems fairly unlikely for the solar system, or alternately, if the protoplanetary disk is truncated at 1.0 AU, simulations have been able to produce Mars-like bodies in the correct location. In this paper, we examine an alternative scenario for the formation of Mars in which a local depletion in the density of the protosolar nebula results in a non-uniform formation of planetary embryos and ultimately the formation of Mars-sized planets around 1.5 AU. We have carried out extensive numerical simulations of the formation of terrestrial planets in such a disk for different scales of the local density depletion, and for different orbital configurations of the giant planets. Our simulations point to the possibility of the formation of Mars-sized bodies around 1.5 AU, specifically when the scale of the disk local mass-depletion is moderately high (50-75%) and Jupiter and Saturn are initially in their current orbits. In these systems, Mars-analogs are formed from the protoplanetary materials that originate in the regions of disk interior or exterior to the local mass-depletion. Results also indicate that Earth-sized planets can form around 1 AU with a substantial amount of water accreted via primitive water-rich planetesimals and planetary embryos. We present the results of our study and discuss their implications for the formation of terrestrial planets in our solar system.

One of the most important subjects of debate in the formation of the solar system is the origin of Earth's water. Comets have long been considered as the most likely source of the delivery of water to Earth. However, elemental and isotopic arguments suggest a very small contribution from these objects. Other sources have also been proposed, among which, local adsorption of water vapor onto dust grains in the primordial nebula and delivery through planetesimals and planetary embryos have become more prominent. However, no sole source of water provides a satisfactory explanation for Earth's water as a whole. In view of that, using numerical simulations, we have developed a compound model incorporating both the principal endogenous and exogenous theories, and investigating their implications for terrestrial planet formation and water-delivery. Comets are also considered in the final analysis, as it is likely that at least some of Earth's water has cometary origin. We analyze our results comparing two different water distribution models, and complement our study using D/H ratio, finding possible relative contributions from each source, focusing on planets formed in the habitable zone. We find that the compound model play an important role by showing more advantage in the amount and time of water-delivery in Earth-like planets.

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.

The relative distribution of abundances of refractory, intermediate, and volatile elements in stars with planets can be an important tool for investigating the internal migration of a giant planet. This migration can lead to the accretion of planetesimals and the selective enrichment of the star with these elements. We report on a spectroscopic determination of the atmospheric parameters and chemical abundances of the parent stars in transiting planets CoRoT-2b and CoRoT-4b. Adding data for CoRoT-3 and CoRoT-5 from the literature, we find a flat distribution of the relative abundances as a function of their condensation temperatures. For CoRoT-2, the relatively high lithium abundance and intensity of its Li I resonance line permit us to propose an age of 120 Myr, making this stars one of the youngest stars with planets to date. We introduce a new methodology to investigate a relation between the abundances of these stars and the internal migration of their planets. By simulating the internal migration of a planet in a disk formed only by planetesimals, we are able to separate the stellar fractions of refractory (R), intermediate (I), and volatile (V) rich planetesimals accreting onto the central star. Intermediate and volatile element fractions enriching the star are similar and much larger than those of pure refractory ones. We also show that these results are highly dependent on the model adopted for the disk distribution regions in terms of R, I, and V elements and other parameters considered. We note however, that this self-enrichment mechanism is only efficient during the first 20-30 Myr or later in the lifetime of the disk when the surface convection layers of the central star for the first time attain its minimum size configuration.

Saturn is the only known planet to have coorbital satellite systems. In the present work we studied the process of mass accretion as a possible mechanism for coorbital satellites formation. The system considered is composed of Saturn, a proto-satellite and a cloud of planetesimals distributed in the coorbital region around a triangular Lagrangian point. The adopted relative mass for the proto-satellite was 10^-6 of Saturn's mass and for each planetesimal of the cloud three cases of relative mass were considered, 10^-14, 10^-13 and 10^-12 masses of Saturn. In the simulations each cloud of planetesimal was composed of 10^3, 5 x 10^3 or 10^4 planetesimals. The results of the simulations show the formation of coorbital satellites with relative masses of the same order of those found in the saturnian system (10^-13 - 10^-9). Most of them present horseshoe type orbits, but a significant part is in tadpole orbit around L_4 or L_5. Therefore, the results indicate that this is a plausible mechanism for the formation of coorbital satellites.

The origins of irregular satellites of the giant planets are an important piece of the giant "puzzle" that is the theory of Solar System formation. It is well established that they are not "in situ" formation objects, around the planet, as are believed to be the regular ones. Then, the most plausible hypothesis to explain their origins is that they formed elsewhere and were captured by the planet. However, captures under restricted three-body problem dynamics have temporary feature, which makes necessary the action of an auxiliary capture mechanism. Nevertheless, there not exist one well established capture mechanism. In this work, we tried to understand which aspects of a binary-asteroid capture mechanism could favor the permanent capture of one member of a binary asteroid. We performed more than eight thousand numerical simulations of capture trajectories considering the four-body dynamical system Sun, Jupiter, Binary-asteroid. We restricted the problem to the circular planar prograde case, and time of integration to 10^4 years. With respect to the binary features, we noted that 1) tighter binaries are much more susceptible to produce permanent captures than the large separation-ones. We also found that 2) the permanent capture probability of the minor member of the binary is much more expressive than the major body permanent capture probability. On the other hand, among the aspects of capture-disruption process, 4) a pseudo eastern-quadrature was noted to be a very likely capture angular configuration at the instant of binary disruptions. In addition, we also found that the 5) capture probability is higher for binary asteroids which disrupt in an inferior-conjunction with Jupiter. These results show that the Sun plays a very important role on the capture dynamic of binary asteroids.

he triple asteroidal system (87) Sylvia is composed of a 280-km primary and two small moonlets named Romulus and Remus (Marchis et al 2005). Sylvia is located in the main asteroid belt. The satellites are in nearly equatorial circular orbits around the primary. In the present work we study the stability of the satellites Romulus and Remus, in order to identify the effects and the contribution of each perturber. The results from the 3-body problem, Sylvia-Romulus-Remus, show no significant variation of their orbital elements. However, the inclinations of the satellites present a long period evolution, when the Sun is included in the system. Such amplitude is amplified when Jupiter is included. An analysis of these results show that Romulus and Remus are librating in a secular resonance and their longitude of the nodes are locked to each other. The satellites get caught in an evection resonance with Jupiter. However, the orbital evolutions of the satellites became completely stable when the oblateness of Sylvia is included in the simulations.

The photospheres of stars hosting planets have larger metallicity than stars lacking planets. In the present work we study the possibility of an earlier metal enrichment of the photospheres by means of impacting planetesimals during the first 20-30Myr. Here we explore this contamination process by simulating the interactions of an inward migrating planet with a disc of planetesimal interior to its orbit. The results show the percentage of planetesimals that fall on the star. We identified the dependence of the planet's eccentricity ($e_p$) and time scale of migration ($\tau$) on the rate of infalling planetesimals. For very fast migrations ($\tau=10^2$yr and $\tau=10^3$yr) there is no capture in mean motion resonances, independently of the value of $e_p$. Then, due to the planet's migration the planetesimals suffer close approaches with the planet and more than 80% of them are ejected from the system. For slow migrations ($\tau=10^5$yr and $\tau=10^6$yr) the percentage of collisions with the planet decrease with the increase of the planet's eccentricity. For $e_p=0$ and $e_p=0.1$ most of the planetesimals were captured in the 2:1 resonance and more than 65% of them collided with the star. Whereas migration of a Jupiter mass planet to very short pericentric distances requires unrealistic high disc masses, these requirements are much smaller for smaller migrating planets. Our simulations for a slowly migrating 0.1 $M_{\rm Jupiter}$ planet, even demanding a possible primitive disc three times more massive than a primitive solar nebula, produces maximum [Fe/H] enrichments of the order of 0.18 dex. These calculations open possibilities to explain hot Jupiters exoplanets metallicities.