Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.

In the framework of the Global Architecture of Planetary Systems (GAPS) project we collected more than 300 spectra with HARPS-N at the TNG for the bright G9V star HD164922. This target is known to host one gas giant planet in a wide orbit (Pb~1200 days, semi-major axis ~2 au) and a Neptune-mass planet with a period Pc ~76 days. We searched for additional low-mass companions in the inner region of the system. We compared the radial velocities (RV) and the activity indices derived from the HARPS-N time series to measure the rotation period of the star and used a Gaussian process regression to describe the behaviour of the stellar activity. We exploited this information in a combined model of planetary and stellar activity signals in an RV time-series composed of almost 700 high-precision RVs, both from HARPS-N and literature data. We performed a dynamical analysis to evaluate the stability of the system and the allowed regions for additional potential companions. Thanks to the high sensitivity of the HARPS-N dataset, we detect an additional inner super-Earth with an RV semi-amplitude of 1.3+/-0.2 m/s, a minimum mass of ~4+/-1 M_E and a period of 12.458+/-0.003 days. We disentangle the planetary signal from activity and measure a stellar rotation period of ~42 days. The dynamical analysis shows the long term stability of the orbits of the three-planet system and allows us to identify the permitted regions for additional planets in the semi-major axis ranges 0.18-0.21 au and 0.6-1.4 au. The latter partially includes the habitable zone of the system. We did not detect any planet in these regions, down to minimum detectable masses of 5 and 18 M_E, respectively. A larger region of allowed planets is expected beyond the orbit of planet b, where our sampling rules-out bodies with minimum mass > 50 M_E.

Binary systems are very common among field stars. While this relatively small number of planets in binaries is probably partly due to strong observational biases, there is, however, statistical evidence that planets are indeed less frequent in binaries with separations smaller than 100 au, strongly suggesting that the presence of a close in companion star has an adverse effect on planet formation. It is indeed possible for the gravitational pull of the second star to affect all the different stages of planet formation, from proto-planetary disk formation to dust accumulation into planetesimals, to the accretion of these planetesimals into large planetary embryos and, eventually, the final growth of these embryos into planets. For the crucial planetesimal accretion phase, the complex coupling between dynamical perturbations from the binary and friction due to gas in the protoplanetary disk suggests that planetesimal accretion might be hampered due to increased, accretion hostile impact velocities. Likewise, the interplay between the binary secular perturbations and mean motion resonances lead to unstable regions, where not only planet formation is inhibited, but where a massive body would be ejected from the system on a hyperbolic orbit. The amplitude of these two main effects is different for S and P type planets, so that a comparison between the two populations might outline the influence of the companion star on the planet formation process. Unfortunately, at present the two populations (circumstellar or circumbinary) are not known equally well and different biases and uncertainties prevent a quantitative comparison. We also highlight the long term dynamical evolution of both S and P type systems and focus on how these different evolutions influence the final architecture of planetary systems in binaries.

Pairs of planets in a system may end up close to their host star on eccentric orbits as a consequence of planet planet scattering, Kozai or secular migration. In this scenario, general relativity and secular perturbations have comparable timescales and may interfere with each other with relevant effects on the eccentricity and pericenter evolution of the two planets. We explore, both analytically and via numerical integration, how the secular evolution is changed by general relativity for a wide range of different initial conditions. We find that when the faster secular frequency approaches the general relativity precession rate, which tipically occurs when the outer planet moves away from the inner one, it relaxes to it and a significant damping of the proper eccentricity of the inner planet occurs. The proper eccentricity of the outer planet is reduced as well due to the changes in the secular interaction of the bodies. The lowering of the peak eccentricities of the two planets during their secular evolution has important implications on their stability. A significant number of two planet systems, otherwise chaotic because of the mutual secular perturbations, are found stable when general relativity is included.

The aim of this paper is to underline conditions necessary for the emergence and development of life. They are placed at local planetary scale, at Galactic scale and within the cosmological evolution, as pointed out by the Anthropic Cosmological Principle. We will consider the circumstellar habitable zone (CHZ) for planetary systems and a Galactic Habitable Zone (GHZ) including also a set of strong cosmological constraints to allow life (cosmological habitability (COSH)). Some requirements are specific of a single scale and its related physical phenomena, while others are due to the conspired effects occurring at more than one scale. The scenario emerging from this analysis is that all the habitability conditions here detailed must at least be met. Then life asks something which may appear as "a monstrous sequence of accidents" as Hoyle (1959) thought, or as providential collaborations leading to realize how finely tuned is the architecture within which Life, as precious stone, is embedded.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

The amount of dust present in circumstellar disks is expected to steadily decrease with age due to the growth from micron-sized particles to planetesimals and planets. Mature circumstellar disks, however, can be observed to contain significant amounts of dust and possess high dust-to-gas ratios. Using HD 163296 as our case study, we explore how the formation of giant planets in disks can create the conditions for collisionally rejuvenating the dust population, halting or reversing the expected trend. We combine N-body simulations with statistical methods and impact scaling laws to estimate the dynamical and collisional excitation of the planetesimals due to the formation of HD 163296's giant planets. We show that this process creates a violent collisional environment across the disk that can inject collisionally produced second-generation dust into it, significantly contributing to the observed dust-to-gas ratio. The spatial distribution of the dust production can explain the observed local enrichments in HD 163296's inner regions. The results obtained for HD 163296 can be extended to any disk with embedded forming giant planets and may indicate a common evolutionary stage in the life of such circumstellar disks. Furthermore, the dynamical excitation of the planetesimals could result in the release of transient, non-equilibrium gas species like H2O, CO2, NH3 and CO in the disk due to ice sublimation during impacts and, due to the excited planetesimals being supersonic with respect to the gas, could produce bow shocks in the latter that could heat it and cause a broadening of its emission lines.

Planet--Planet scattering is an efficient and robust dynamical mechanism for producing eccentric exoplanets. Coupled to tidal interactions with the central star, it can also explain close--in giant planets on circularized and potentially misaligned orbits. We explore scattering events occurring close to the star and test if they can reproduce the main features of the observed orbital distribution of giant exoplanets on tight orbits.In our modeling we exploit a numerical integration code based on the Hermite algorithm and including the effects of general relativity, dynamical tides and two--body collisions.We find that P--P scattering events occurring in systems with three giant planets initially moving on circular orbits close to their star produce a population of planets similar to the presently observed one, including eccentric and misaligned close--in planets. The contribution of tides and general relativity is relevant in determining the final outcome of the chaotic phase. Even if two--body collisions dominate the chaotic evolution of three planets in crossing orbits close to their star, the final distribution shows a significant number of planets on eccentric orbits. The highly misaligned close--in giant planets are instead produced by systems where the initial semi--major axis of the inner planet was around 0.2 au or beyond.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

We investigate via numerical modeling the effects of two planets locked in resonance, and migrating outward, on the dust distribution of the natal circumstellar disk. We aim to test whether the dust distribution exhibits peculiar features arising from the interplay among the gravitational perturbations of the planets in resonance, the evolution of the gas, and its influence on the dust grains' dynamics. We focus on the 3:2 and 2:1 resonance, where the trapping may be caused by the convergent migration of a Jupiter- and Saturn-mass planet, preceding the common gap formation and ensuing outward (or inward) migration. Models show that a common gap also forms in the dust component -- similarly to what a single, more massive planet would generate -- and that outward migration leads to a progressive widening of the dust gap and to a decoupling from the gas gap. As the system evolves, a significantly wider gap is observed in the dust distribution, which ceases to overlap with the gas gap in the inner disk regions. At the outer edge of the gas gap, outward migration of the planets produces an over-density of dust particles, which evolve differently in the 3:2 and 2:1 resonances. For the 3:2, the dust trap at the gap's outer edge is partly efficient and a significant fraction of the grains filters through the gap. For the 2:1 resonance, the trap is more efficient and very few grains cross the gap, while the vast majority accumulate at the outer edge of the gap.

The phase function of the dust coma of comet 67P has been determined from Rosetta/OSIRIS images \citep{Bertini17}. This function show a deep minimum at phase angles near 100$^\circ$, and a strong backscattering enhancement. These two properties cannot be reproduced by regular models of cometary dust, most of them based on wavelength-sized and randomly-oriented aggregate particles. We show, however, that an ensamble of oriented elongated particles of a wide variety of aspect ratios, with radii $r \gtrsim$10 $\mu$m, and whose long axes are perpendicular to the direction of the solar radiation, are capable of reproducing the observed phase function. These particles must be absorbing, with an imaginary part of the refractive index of about 0.1 to match the expected geometric albedo, and with porosity in the 60-70\% range.

The Rosetta space probe accompanied comet 67P/Churyumov-Gerasimenko for more than two years, obtaining an unprecedented amount of unique data of the comet nucleus and inner coma. This work focuses identifying the source regions of faint jets and outbursts and on studying the spectrophotometric properties of some outbursts. We use observations acquired with the OSIRIS/NAC camera during July-October 2015, that is, close to perihelion. More than 200 jets of different intensities were identified directly on the nucleus. Some of the more intense outbursts appear spectrally bluer than the comet dark terrain in the vivible-to-near-infrared region. We attribute this spectral behavior to icy grains mixed with the ejected dust. Some of the jets have an extremely short lifetime. They appear on the cometary surface during the color sequence observations, and vanish in less than some few minutes after reaching their peak. We also report a resolved dust plume observed in May 2016 at a resolution of 55 cm/pixel, which allowed us to estimate an optical depth of $\sim$0.65 and an ejected mass of $\sim$ 2200 kg. We present the results on the location, duration, and colors of active sources on the nucleus of 67P from the medium-resolution (i.e., 6-10 m/pixel) images acquired close to perihelion passage. The observed jets are mainly located close to boundaries between different morphological regions. Jets depart not only from cliffs, but also from smooth and dust-covered areas, from fractures, pits, or cavities that cast shadows and favor the recondensation of volatiles. This study shows that faint jets or outbursts continuously contribute to the cometary activity close to perihelion passage, and that these events are triggered by illumination conditions. Faint jets or outbursts are not associated with a particular terrain type or morphology.

The nearby Herbig Be star HD100546 is known to be a laboratory for the study of protoplanets and their relation with the circumstellar disk that is carved by at least 2 gaps. We observed the HD100546 environment with high contrast imaging exploiting several different observing modes of SPHERE, including datasets with/without coronagraphs, dual band imaging, integral field spectroscopy and polarimetry. The picture emerging from these different data sets is complex. Flux-conservative algorithms images clearly show the disk up to 200au. More aggressive algorithms reveal several rings and warped arms overlapping the main disk. The bright parts of this ring lie at considerable height over the disk mid-plane at about 30au. Our images demonstrate that the brightest wings close to the star in the near side of the disk are a unique structure, corresponding to the outer edge of the intermediate disk at ~40au. Modeling of the scattered light from the disk with a geometrical algorithm reveals that a moderately thin structure can well reproduce the light distribution in the flux-conservative images. We suggest that the gap between 44 and 113 au span between the 1:2 and 3:2 resonance orbits of a massive body located at ~70au that might coincide with the candidate planet HD100546b detected with previous thermal IR observations. In this picture, the two wings can be the near side of a ring formed by disk material brought out of the disk at the 1:2 resonance with the same massive object. While we find no clear evidence confirming detection of the planet candidate HD100546c in our data, we find a diffuse emission close to the expected position of HD100546b. This source can be described as an extremely reddened substellar object surrounded by a dust cloud or its circumplanetary disk. Its astrometry is broadly consistent with a circular orbital motion on the disk plane.

Binary stars constitute a large percentage of the stellar population, yet relatively little is known about the planetary systems orbiting them. Most constraints on circumbinary planets (CBPs) so far come from transit observations with the Kepler telescope, which is sensitive to close-in exoplanets but does not constrain planets on wider orbits. However, with continuous developments in high-contrast imaging techniques, this population can now be addressed through direct imaging. We present the full survey results of the Search for Planets Orbiting Two Stars (SPOTS) survey, which is the first direct imaging survey targeting CBPs. The SPOTS observational program comprises 62 tight binaries that are young and nearby, and thus suitable for direct imaging studies, with VLT/NaCo and VLT/SPHERE. Results from SPOTS include the resolved circumbinary disk around AK Sco, the discovery of a low-mass stellar companion in a triple packed system, the relative astrometry of up to 9 resolved binaries, and possible indications of non-background planetary-mass candidates around HIP 77911. We did not find any CBP within 300 AU, which implies a frequency upper limit on CBPs (1--15 $M_{\rm Jup}$) of 6--10 % between 30-300 AU. Coupling these observations with an archival dataset for a total of 163 stellar pairs, we find a best-fit CBP frequency of 1.9 % (2--15 $M_{\rm Jup}$) between 1--300 AU with a 10.5 % upper limit at a 95 % confidence level. This result is consistent with the distribution of companions around single stars.

The assessment of the frequency of planetary systems reproducing the Solar System's architecture is still an open problem. Detailed study of multiplicity and architecture is generally hampered by limitations in quality, temporal extension and observing strategy, causing difficulties in detecting low-mass inner planets in the presence of outer giant planetary bodies. We present the results of high-cadence and high-precision HARPS observations on 20 solar-type stars known to host a single long-period giant planet in order to search for additional inner companions and estimate the occurence rate $f_p$ of scaled Solar System analogs, i.e. systems featuring lower-mass inner planets in the presence of long-period giant planets. We carry out combined fits of our HARPS data with literature radial velocities using differential evolution MCMC to refine the literature orbital solutions and search for additional inner planets. We then derive the survey detection limits to provide preliminary estimates of $f_p$. We generally find better constrained orbital parameters for the known planets than those found in the literature. While no additional inner planet is detected, we find evidence for previously unreported long-period massive companions in systems HD 50499 and HD 73267. We finally estimate the frequency of inner low mass (10-30 M$_\oplus$) planets in the presence of outer giant planets as $f_p<9.84\%$ for P<150 days. Our preliminary estimate of $f_p$ is significantly lower than the values found in the literature; the lack of inner candidate planets found in our sample can also be seen as evidence corroborating the inward migration formation model for super-Earths and mini-Neptunes. Our results also underline the need for high-cadence and high-precision follow-up observations as the key to precisely determine the occurence of Solar System analogs.

We present 20 years of radial velocity (RV) measurements of the M1 dwarf Gl15A, combining 5 years of intensive RV monitoring with the HARPS-N spectrograph with 15 years of archival HIRES/Keck RV data. We carry out an MCMC-based analysis of the RV time series, inclusive of Gaussian Process (GP) approach to the description of stellar activity induced RV variations. Our analysis confirms the Keplerian nature and refines the orbital solution for the 11.44-day period super Earth, Gl15A\,b, reducing its amplitude to $1.68^{+0.17}_{-0.18}$ m s$^{-1}$ ($M \sin i = 3.03^{+0.46}_{-0.44}$ M$_\oplus$), and successfully models a long-term trend in the combined RV dataset in terms of a Keplerian orbit with a period around 7600 days and an amplitude of $2.5^{+1.3}_{-1.0}$ m s$^{-1}$, corresponding to a super-Neptune mass ($M \sin i = 36^{+25}_{-18}$ M$_\oplus$) planetary companion. We also discuss the present orbital configuration of Gl15A planetary system in terms of the possible outcomes of Lidov-Kozai interactions with the wide-separation companion Gl15B in a suite of detailed numerical simulations. In order to improve the results of the dynamical analysis, we derive a new orbital solution for the binary system, combining our RV measurements with astrometric data from the WDS catalogue. The eccentric Lidov-Kozai analysis shows the strong influence of Gl15B on the Gl15A planetary system, which can produce orbits compatible with the observed configuration for initial inclinations of the planetary system between $75^\circ$ and $90^\circ$, and can also enhance the eccentricity of the outer planet well above the observed value, even resulting in orbital instability, for inclinations around $0^\circ$ and $15^\circ - 30^\circ$. The Gl15A system is the multi-planet system closest to Earth, at $3.57$ pc, and hosts the longest period RV sub-jovian mass planet discovered so far.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

Context: In the early evolution of a planetary system, a pair of planets may be captured in a mean motion resonance while still embedded in their nesting circumstellar disk. Aims: The goal is to estimate the direction and amount of shift in the semimajor axis of the resonance location due to the disk gravity as a function of the gas density and mass of the planets. The stability of the resonance lock when the disk dissipates is also tested. Methods: The orbital evolution of a large number of systems is numerically integrated within a three-body problem in which the disk potential is computed as a series of expansion. This is a good approximation, at least over a limited amount of time. Results: Two different resonances are studied: the 2:1 and the 3:2. In both cases the shift is inwards, even if by a different amount, when the planets are massive and carve a gap in the disk. For super--Earths, the shift is instead outwards. Different disk densities, Sigma, are considered and the resonance shift depends almost linearly on Sigma. The gas dissipation leads to destabilization of a significant number of resonant systems, in particular if it is fast. Conclusions: The presence of a massive circumstellar disk may significantly affect the resonant behavior of a pair of planets by shifting the resonant location and by decreasing the size of the stability region. The disk dissipation may explain some systems found close to a resonance but not locked in it.

On 3 July 2016, several instruments on board ESA's Rosetta spacecraft detected signs of an outburst event on comet 67P, at a heliocentric distance of 3.32 AU from the sun, outbound from perihelion. We here report on the inferred properties of the ejected dust and the surface change at the site of the outburst. The activity coincided with the local sunrise and continued over a time interval of 14 - 68 minutes. It left a 10m-sized icy patch on the surface. The ejected material comprised refractory grains of several hundred microns in size, and sub-micron-sized water ice grains. The high dust mass production rate is incompatible with the free sublimation of crystalline water ice under solar illumination as the only acceleration process. Additional energy stored near the surface must have increased the gas density. We suggest a pressurized sub-surface gas reservoir, or the crystallization of amorphous water ice as possible causes.

A large number of systems harboring a debris disk show evidence for a double belt architecture. One hypothesis for explaining the gap between the belts is the presence of one or more planets dynamically carving it. This work aims to investigate this scenario in systems harboring two components debris disks. All the targets in the sample were observed with the SPHERE instrument which performs high-contrast direct imaging. Positions of the inner and outer belts were estimated by SED fitting of the infrared excesses or, when available, from resolved images of the disk. Very few planets have been observed so far in debris disks gaps and we intended to test if such non-detections depend on the observational limits of the present instruments. This aim is achieved by deriving theoretical predictions of masses, eccentricities and semi-major axes of planets able to open the observed gaps and comparing such parameters with detection limits obtained with SPHERE. The relation between the gap and the planet is due to the chaotic zone around the orbit of the planet. The radial extent of this zone depends on the mass ratio between the planet and the star, on the semi-major axis and on the eccentricity of the planet and it can be estimated analytically. We apply the formalism to the case of one planet on a circular or eccentric orbit. We then consider multi-planetary systems: 2 and 3 equal-mass planets on circular orbits and 2 equal-mass planets on eccentric orbits in a packed configuration. We then compare each couple of values (M,a), derived from the dynamical analysis of single and multiple planetary models, with the detection limits obtained with SPHERE. Our results show that the apparent lack of planets in gaps between double belts could be explained by the presence of a system of two or more planets possibly of low mass and on an eccentric orbits whose sizes are below the present detection limits.

We describe Castalia, a proposed mission to rendezvous with a Main Belt Comet (MBC), 133P/Elst-Pizarro. MBCs are a recently discovered population of apparently icy bodies within the main asteroid belt between Mars and Jupiter, which may represent the remnants of the population which supplied the early Earth with water. Castalia will perform the first exploration of this population by characterising 133P in detail, solving the puzzle of the MBC's activity, and making the first in situ measurements of water in the asteroid belt. In many ways a successor to ESA's highly successful Rosetta mission, Castalia will allow direct comparison between very different classes of comet, including measuring critical isotope ratios, plasma and dust properties. It will also feature the first radar system to visit a minor body, mapping the ice in the interior. Castalia was proposed, in slightly different versions, to the ESA M4 and M5 calls within the Cosmic Vision programme. We describe the science motivation for the mission, the measurements required to achieve the scientific goals, and the proposed instrument payload and spacecraft to achieve these.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The Southern hemisphere of the 67P/Churyumov-Gerasimenko comet has become visible from Rosetta only since March 2015. It was illuminated during the perihelion passage and therefore it contains the regions that experienced the strongest heating and erosion rate, thus exposing the subsurface most pristine material. In this work we investigate, thanks to the OSIRIS images, the geomorphology, the spectrophotometry and some transient events of two Southern hemisphere regions: Anhur and part of Bes. Bes is dominated by outcropping consolidated terrain covered with fine particle deposits, while Anhur appears strongly eroded with elongated canyon-like structures, scarp retreats, different kinds of deposits, and degraded sequences of strata indicating a pervasive layering. We discovered a new 140 m long and 10 m high scarp formed in the Anhur/Bes boundary during/after the perihelion passage, close to the area where exposed CO$_2$ and H$_2$O ices were previously detected. Several jets have been observed originating from these regions, including the strong perihelion outburst, an active pit, and a faint optically thick dust plume. We identify several areas with a relatively bluer slope (i.e. a lower spectral slope value) than their surroundings, indicating a surface composition enriched with some water ice. These spectrally bluer areas are observed especially in talus and gravitational accumulation deposits where freshly exposed material had fallen from nearby scarps and cliffs. The investigated regions become spectrally redder beyond 2 au outbound when the dust mantle became thicker, masking the underlying ice-rich layers.

We carried out a Bayesian homogeneous determination of the orbital parameters of 231 transiting giant planets (TGPs) that are alone or have distant companions; we employed DE-MCMC methods to analyse radial-velocity (RV) data from the literature and 782 new high-accuracy RVs obtained with the HARPS-N spectrograph for 45 systems over 3 years. Our work yields the largest sample of systems with a transiting giant exoplanet and coherently determined orbital, planetary, and stellar parameters. We found that the orbital parameters of TGPs in non-compact planetary systems are clearly shaped by tides raised by their host stars. Indeed, the most eccentric planets have relatively large orbital separations and/or high mass ratios, as expected from the equilibrium tide theory. This feature would be the outcome of high-eccentricity migration (HEM). The distribution of $\alpha=a/a_R$, where $a$ and $a_R$ are the semi-major axis and the Roche limit, for well-determined circular orbits peaks at 2.5; this also agrees with expectations from the HEM. The few planets of our sample with circular orbits and $\alpha >5$ values may have migrated through disc-planet interactions instead of HEM. By comparing circularisation times with stellar ages, we found that hot Jupiters with $a < 0.05$ au have modified tidal quality factors $10^{5} < Q'_p < 10^{9}$, and that stellar $Q'_s > 10^{6}-10^{7}$ are required to explain the presence of eccentric planets at the same orbital distance. As a by-product of our analysis, we detected a non-zero eccentricity for HAT-P-29; we determined that five planets that were previously regarded to have hints of non-zero eccentricity have circular orbits or undetermined eccentricities; we unveiled curvatures caused by distant companions in the RV time series of HAT-P-2, HAT-P-22, and HAT-P-29; and we revised the planetary parameters of CoRoT-1b.

Images of the nucleus and the coma (gas and dust) of comet 67P/Churyumov- Gerasimenko have been acquired by the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras since March 2014 using both the Wide Angle Camera (WAC) and the Narrow Angle Camera (NAC). We use images from the NAC camera to study a bright outburst observed in the southern hemisphere on July 29, 2015. The high spatial resolution of the NAC is needed to localize the source point of the outburst on the surface of the nucleus. The heliocentric distance is 1.25 au and the spacecraft-comet distance is 186 km. Aiming to better understand the physics that led to the outgassing, we used the Direct Simulation Monte Carlo (DSMC) method to study the gas flow close to the nucleus and the dust trajectories. The goal is to understand the mechanisms producing the outburst. We reproduce the opening angle of the outburst in the model and constrain the outgassing ratio between the outburst source and the local region. The outburst is in fact a combination of both gas and dust, in which the active surface is approximately 10 times more active than the average rate found in the surrounding areas. We need a number of dust particles 7.83 $\times$ 10$^{11}$ - 6.90 $\times$ 10$^{15}$ (radius 1.97 - 185 {\mu}m), which corresponds to a mass of dust 220 - 21 $\times$ 10$^{3}$kg.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

Planetesimal belts are invoked to explain the prolonged existence of debris disks. Important parameters to model their collisional evolution and to compute the dust production rate are the intrinsic probability of collision $P_i$ and the mean impact velocity $U_c$. If a planet orbits close to the belt, the values of both these parameters are affected by its secular perturbations yielding a strong correlation between eccentricity $e$ and pericentre longitude $\varpi$. We adopt a new algorithm to compute both $P_i$ and $U_c$ in presence of various levels of secular correlation due to different ratios between proper and forced eccentricity. We tested this algorithm in a standard case with a Jupiter--sized planet orbiting inside a putative planetesimal belt finding that it is less collisionally active compared to a self--stirred belt because of the $e - \varpi$ coupling. The eccentricity of the planet is an important parameter in determining the amount of dust production since the erosion rate is 10 times faster when the planet eccentricity increases from 0.1 to 0.6. Also the initial conditions of the belt (either warm or cold) and its average inclination strongly affects $P_i$ and $U_c$ and then its long term collisional evolution in presence of the planet. We finally apply our method to the planetesimal belts supposedly refilling the dust disks around HD 38529 and $\epsilon$ Eridani. In the most collisionally active configurations, only a small fraction of bodies smaller than 100 km are expected to be fragmented over a time--span of 4 Gyr.

In order to understand the observed physical and orbital diversity of extrasolar planetary systems, a full investigation of these objects and of their host stars is necessary. Within this field, one of the main purposes of the GAPS observing project with HARPS-N@TNG is to provide a more detailed characterisation of already known systems. In this framework we monitored the star, hosting two giant planets, HD108874, with HARPS-N for three years in order to refine the orbits, to improve the dynamical study and to search for additional low-mass planets in close orbits. We subtracted the radial velocity (RV) signal due to the known outer planets, finding a clear modulation of 40.2 d period. We analysed the correlation between RV residuals and the activity indicators and modelled the magnetic activity with a dedicated code. Our analysis suggests that the 40.2 d periodicity is a signature of the rotation period of the star. A refined orbital solution is provided, revealing that the system is close to a mean motion resonance of about 9:2, in a stable configuration over 1 Gyr. Stable orbits for low-mass planets are limited to regions very close to the star or far from it. Our data exclude super-Earths with Msin i \gtrsim 5 M_Earth within 0.4 AU and objects with Msin i \gtrsim 2 M_Earth with orbital periods of a few days. Finally we put constraints on the habitable zone of the system, assuming the presence of an exomoon orbiting the inner giant planet.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

During its two years mission around comet 67P/Churyumov-Gerasimenko, ESA's Rosetta spacecraft had the unique opportunity to follow closely a comet in the most active part of its orbit. Many studies have presented the typical features associated to the activity of the nucleus, such as localized dust and gas jets. Here we report on series of more energetic transient events observed during the three months surrounding the comet's perihelion passage in August 2015. We detected and characterized 34 outbursts with the Rosetta cameras, one every 2.4 nucleus rotation. We identified 3 main dust plume morphologies associated to these events: a narrow jet, a broad fan, and more complex plumes featuring both previous types together. These plumes are comparable in scale and temporal variation to what has been observed on other comets. We present a map of the outbursts source locations, and discuss the associated topography. We find that the spatial distribution sources on the nucleus correlates well with morphological region boundaries, especially in areas marked by steep scarps or cliffs. Outbursts occur either in the early morning or shortly after the local noon, indicating two potential processes: Morning outbursts may be triggered by thermal stresses linked to the rapid change of temperature, afternoon events are most likely related to the diurnal or seasonal heat wave reaching volatiles buried under the first surface layer. In addition, we propose that some events can be the result of a completely different mechanism, in which most of the dust is released upon the collapse of a cliff.

When exploring the stability of multiplanet systems in binaries, two parameters are normally exploited: the critical semimajor axis ac computed by Holman and Wiegert (1999) within which planets are stable against the binary perturbations, and the Hill stability limit Delta determining the minimum separation beyond which two planets will avoid mutual close encounters. Our aim is to test whether these two parameters can be safely applied in multiplanet systems in binaries or if their predictions fail for particular binary orbital configurations. We have used the frequency map analysis (FMA) to measure the diffusion of orbits in the phase space as an indicator of chaotic behaviour. First we revisited the reliability of the empirical formula computing ac in the case of single planets in binaries and we find that, in some cases, it underestimates by 10-20% the real outer limit of stability. For two planet systems, the value of Delta is close to that computed for planets around single stars, but the level of chaoticity close to it substantially increases for smaller semimajor axes and higher eccentricities of the binary orbit. In these configurations ac also begins to be unreliable and non linear secular resonances with the stellar companion lead to chaotic behaviour well within ac, even for single planet systems. For two planet systems, the superposition of mean motion resonances, either mutual or with the binary companion, and non linear secular resonances may lead to chaotic behaviour in all cases. We have developed a parametric semiempirical formula determining the minimum value of the binary semimajor axis, for a given eccentricity of the binary orbit, below which stable two planet systems cannot exist.

Beginning in March 2014, the OSIRIS (Optical, Spectroscopic, and Infrared Remote Imaging System) cameras began capturing images of the nucleus and coma (gas and dust) of comet 67P/Churyumov-Gerasimenko using both the wide angle camera (WAC) and the narrow angle camera (NAC). The many observations taken since July of 2014 have been used to study the morphology, location, and temporal variation of the comet's dust jets. We analyzed the dust monitoring observations shortly after the southern vernal equinox on May 30 and 31, 2015 with the WAC at the heliocentric distance Rh = 1.53 AU, where it is possible to observe that the jet rotates with the nucleus. We found that the decline of brightness as a function of the distance of the jet is much steeper than the background coma, which is a first indication of sublimation. We adapted a model of sublimation of icy aggregates and studied the effect as a function of the physical properties of the aggregates (composition and size). The major finding of this article was that through the sublimation of the aggregates of dirty grains (radius a between 5 microm and 50 microm) we were able to completely reproduce the radial brightness profile of a jet beyond 4 km from the nucleus. To reproduce the data we needed to inject a number of aggregates between 8.5 x $10^{13}$ and 8.5 x $10^{10}$ for a = 5 microm and 50 microm respectively, or an initial mass of $H_2O$ ice around 22kg.

We present OSIRIS/NAC observations of decimetre-sized, likely ice-containing aggregates ejected from a confined region on the surface of comet 67P/Churyumov-Gerasimenko. The images were obtained in January 2016 when the comet was at 2 AU from the Sun out-bound from perihelion. We measure the acceleration of individual aggregates through a two-hour image series. Approximately 50% of the aggregates are accelerated away from the nucleus, and 50% towards it, and likewise towards either horizontal direction. The accelerations are up to one order of magnitude stronger than local gravity, and are most simply explained by the combined effect of gas drag accelerating all aggregates upwards, and the recoil force from asymmetric outgassing, either from rotating aggregates with randomly oriented spin axes and sufficient thermal inertia to shift the temperature maximum away from an aggregate's subsolar region, or from aggregates with variable ice content. At least 10% of the aggregates will escape the gravity field of the nucleus and feed the comet's debris trail, while others may fall back to the surface and contribute to the deposits covering parts of the northern hemisphere. The rocket force plays a crucial role in pushing these aggregates back towards the surface. Our observations show the future back fall material in the process of ejection, and provide the first direct measurement of the acceleration of aggregates in the innermost coma (<2km) of a comet, where gas drag is still significant.

The Wide Angle Camera of the OSIRIS instrument on board the Rosetta spacecraft is equipped with several narrowband filters that are centered on the emission lines and bands of various fragment species. These are used to determine the evolution of the production and spatial distribution of the gas in the inner coma of comet 67P with time and heliocentric distance, here between 2.6 - 1.3 AU pre-perihelion. Our observations indicate that the emission observed in the OH, OI, CN, NH, and NH2 filters is mostly produced by dissociative electron impact excitation of different parent species. We conclude that CO2 rather than H2O is a significant source of the [OI] 630 nm emission. A strong plume-like feature observed in the in CN and [OI] filters is present throughout our observations. This plume is not present in OH emission and indicates a local enhancement of the CO2/H2O ratio by as much as a factor of 3. We observed a sudden decrease in intensity levels after March 2015, which we attribute to decreased electron temperatures in the first kilometers above the nucleus surface.

Aims. In this work, we characterize the low mass companion of the A0 field star HR3549. Methods. We observed HR3549AB in imaging mode with the the NIR branch (IFS and IRDIS) of SPHERE@VLT, with IFS in YJ mode and IRDIS in the H band. We also acquired a medium resolution spectrum with the IRDIS long slit spectroscopy mode. The data were reduced using the dedicated SPHERE GTO pipeline, purposely designed for this instrument. We employed algorithms such as PCA and TLOCI to reduce the speckle noise. Results. The companion was clearly visible both with IRDIS and IFS.We obtained photometry in four different bands as well as the astrometric position for the companion. Based on our astrometry, we confirm that it is a bound object and put constraints on its orbit. Although several uncertainties are still present, we estimate an age of ~100-150 Myr for this system, yielding a most probable mass for the companion of 40-50MJup and T_eff ~300-2400 K. Comparing with template spectra points to a spectral type between M9 and L0 for the companion, commensurate with its position on the color-magnitude diagram.

From the results of a comprehensive asteroid population evolution model, we conclude that the YORP-induced rotational fission hypothesis can be consistent with the observed population statistics of small asteroids in the main belt including binaries and contact binaries. The foundation of this model is the asteroid rotation model of Marzari et al. (2011), which incorporates both the YORP effect and collisional evolution. This work adds to that model the rotational fission hypothesis and the binary evolution model of Jacobson & Scheeres (2011). The asteroid population evolution model is highly constrained by these and other previous works, and therefore it has only two significant free parameters: the ratio of low to high mass ratio binaries formed after rotational fission events and the mean strength of the binary YORP (BYORP) effect. We successfully reproduce characteristic statistics of the small asteroid population: the binary fraction, the fast binary fraction, steady-state mass ratio fraction and the contact binary fraction. We find that in order for the model to best match observations, rotational fission produces high mass ratio ($> 0.2$) binary components with four to eight times the frequency as low mass ratio ($< 0.2$) components, where the mass ratio is the mass of the secondary component divided by the mass of the primary component. This is consistent with post-rotational fission binary system mass ratio being drawn from either a flat or a positive and shallow distribution, since the high mass ratio bin is four times the size of the low mass ratio bin; this is in contrast to the observed steady-state binary mass ratio, which has a negative and steep distribution. This can be understood in the context of the BYORP-tidal equilibrium hypothesis, which predicts that low mass ratio binaries survive for a significantly longer period of time than high mass ratio systems.

(abridged) We present 2-dimensional hydrodynamic simulations using the Smoothed Particle Hydrodynamic (SPH) code, VINE, to model a self-gravitating binary system similar to the GG Tau A system. We simulate systems configured with semi-major axes of either $a=62$~AU (`wide') or $a=32$~AU (`close'), and with eccentricity of either $e=0$ or $e=0.3$. Strong spiral structures are generated with large material streams extending inwards. A small fraction accretes onto the circumstellar disks, with most returning to the torus. Structures also propagate outwards, generating net outwards mass flow and eventually losing coherence at large distances. The torus becomes significantly eccentric in shape. Accretion onto the stars occurs at a rate of a few $\times10^{-8}$\msun/yr implying disk lifetimes shorter than $\sim10^4$~yr, without replenishment. Only wide configurations retain disks by virtue of robust accretion. In eccentric configurations, accretion is episodic, occurs preferentially onto the secondary at wrates peaked near binary periapse. We conclude that the \ggtaua\ torus is strongly self gravitating and that a major contribution to its thermal energy is shock dissipation. We interpret its observed features as manifestations of spiral structures and the low density material surrounding it as an excretion disk created by outward mass flux. We interpret GG Tau A as a coplanar system with an eccentric torus, and account for its supposed mutual inclination as due to degeneracy between the interpretation of inclination and eccentricity. Although the disks persist for long enough to permit planet formation, the environment remains unfavorable due to high temperatures. We conclude that the GG Tau A system is in an eccentric, $a\sim62$~AU orbit.

Open cluster (OC) stars share the same age and metallicity, and, in general, their age and mass can be estimated with higher precision than for field stars. For this reason, OCs are considered an important laboratory to study the relation between the physical properties of the planets and those of their host stars, and the evolution of planetary systems. We started an observational campaign within the GAPS collaboration to search for and characterize planets in OCs We monitored the Praesepe member Pr0211 to improve the eccentricity of the Hot-Jupiter (HJ) already known to orbit this star and search for additional planets. An eccentric orbit for the HJ would support a planet-planet scattering process after its formation. From 2012 to 2015, we collected 70 radial velocity (RV) measurements with HARPS-N and 36 with TRES of Pr0211. Simultaneous photometric observations were carried out with the robotic STELLA telescope in order to characterize the stellar activity. We discovered a long-term trend in the RV residuals that we show to be due to the presence of a second, massive, outer planet. Orbital parameters for the two planets are derived by simultaneously fitting RVs and photometric light curves, with the activity signal modelled as a series of sinusoids at the rotational period of the star and its harmonics. We confirm that Pr0211b has a nearly circular orbit ($e = 0.02 \pm 0.01$), with an improvement of a factor two with respect to the previous determination of its eccentricity, and estimate that Pr0211c has a mass $M_p\sin i = 7.9 \pm 0.2 M_J$, a period $P>$3500 days and a very eccentric orbit ($e>$0.60). Such peculiar systems may be typical of open clusters if the planet-planet scattering phase leading to the formation of HJs is caused by stellar encounters rather than unstable primordial orbits. Pr0211 is the first multi-planet system discovered around an OC star. (abridged)

The ESA Rosetta spacecraft, currently orbiting around comet 67P, has already provided in situ measurements of the dust grain properties from several instruments, particularly OSIRIS and GIADA. We propose adding value to those measurements by combining them with ground-based observations of the dust tail to monitor the overall, time-dependent dust-production rate and size distribution. To constrain the dust grain properties, we take Rosetta OSIRIS and GIADA results into account, and combine OSIRIS data during the approach phase (from late April to early June 2014) with a large data set of ground-based images that were acquired with the ESO Very Large Telescope (VLT) from February to November 2014. A Monte Carlo dust tail code has been applied to retrieve the dust parameters. Key properties of the grains (density, velocity, and size distribution) were obtained from Rosetta observations: these parameters were used as input of the code to considerably reduce the number of free parameters. In this way, the overall dust mass-loss rate and its dependence on the heliocentric distance could be obtained accurately. The dust parameters derived from the inner coma measurements by OSIRIS and GIADA and from distant imaging using VLT data are consistent, except for the power index of the size-distribution function, which is $\alpha$=--3, instead of $\alpha$=--2, for grains smaller than 1 mm. This is possibly linked to the presence of fluffy aggregates in the coma. The onset of cometary activity occurs at approximately 4.3 au, with a dust production rate of 0.5 kg/s, increasing up to 15 kg/s at 2.9 au. This implies a dust-to-gas mass ratio varying between 3.8 and 6.5 for the best-fit model when combined with water-production rates from the MIRO experiment.

We analyze the physical properties and dynamical origin of a curved jet of comet 67P/Churyumov-Gerasimenko that was observed repeatedly in several nucleus rotations starting on May 30 and persisting until early August, 2015. We simulated the motion of dust grains ejected from the nucleus surface under the influence of the gravity and viscous drag effect of the expanding gas flow from the rotating nucleus. The formation of the curved jet is a combination of the size of the dust particles (~0.1-1 mm) and the location of the source region near the nucleus equator. This enhances the spiral feature of the collimated dust stream after the dust is accelerated to a terminal speed on the order of m/s.

The Search for Planets Orbiting Two Stars (SPOTS) survey aims to study the formation and distribution of planets in binary systems by detecting and characterizing circumbinary planets and their formation environments through direct imaging. With the SPHERE Extreme Adaptive Optics instrument, a good contrast can be achieved even at small (<300 mas) separations from bright stars, which enables studies of planets and disks in a separation range that was previously inaccessible. Here, we report the discovery of resolved scattered light emission from the circumbinary disk around the well-studied young double star AK Sco, at projected separations in the ~13--40 AU range. The sharp morphology of the imaged feature is surprising, given the smooth appearance of the disk in its spectral energy distribution. We show that the observed morphology can be represented either as a highly eccentric ring around AK Sco, or as two separate spiral arms in the disk, wound in opposite directions. The relative merits of these interpretations are discussed, as well as whether these features may have been caused by one or several circumbinary planets interacting with the disk.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.

Binary stars hosting exoplanets are a unique laboratory where chemical tagging can be performed to measure with high accuracy the elemental abundances of both stellar components, with the aim to investigate the formation of planets and their subsequent evolution. Here, we present a high-precision differential abundance analysis of the XO-2 wide stellar binary based on high resolution HARPS-N@TNG spectra. Both components are very similar K-dwarfs and host planets. Since they formed presumably within the same molecular cloud, we expect they should possess the same initial elemental abundances. We investigate if the presence of planets can cause some chemical imprints in the stellar atmospheric abundances. We measure abundances of 25 elements for both stars with a range of condensation temperature $T_{\rm C}=40-1741$ K, achieving typical precisions of $\sim 0.07$ dex. The North component shows abundances in all elements higher by $+0.067 \pm 0.032$ dex on average, with a mean difference of +0.078 dex for elements with $T_{\rm C} > 800$ K. The significance of the XO-2N abundance difference relative to XO-2S is at the $2\sigma$ level for almost all elements. We discuss the possibility that this result could be interpreted as the signature of the ingestion of material by XO-2N or depletion in XO-2S due to locking of heavy elements by the planetary companions. We estimate a mass of several tens of $M_{\oplus}$ in heavy elements. The difference in abundances between XO-2N and XO-2S shows a positive correlation with the condensation temperatures of the elements, with a slope of $(4.7 \pm 0.9) \times 10^{-5}$ dex K$^{-1}$, which could mean that both components have not formed terrestrial planets, but that first experienced the accretion of rocky core interior to the subsequent giant planets.

We explore the dynamical evolution of a planet embedded in a disk surrounding a star part of a binary system where the orbital plane of the binary is significantly tilted respect to the initial disk plane. Our aim is to test whether the planet remains within the disk and continues to migrate towards the star in a Type I/II mode in spite of the secular perturbations of the companion star. This would explain observed exoplanets with significant inclination respect to the equatorial plane of their host star. We have used two different SPH codes, vine and phantom, to model the evolution of a system star+disk+planet and companion star with time. After an initial coupled evolution, the inclination of the disk and that of the planet begin to differ significantly. The period of oscillation of the disk inclination, respect to the initial plane, is shorter than that of the planet which evolves independently after about 10^4 yr following a perturbed N-body behavior. However, the planet keeps migrating towards the star because during its orbital motion it crosses the disk plane and the friction with the gas causes angular momentum loss. Disk and planet in a significantly inclined binary system are not dynamically coupled for small binary separations but evolve almost independently. The planet abandons the disk and, due to the onset of a significant mutual inclination, it interacts with the gas only when its orbit intersects the disk plane. The drift of the planet towards the star is not due to type I/II with the planet embedded in the disk but to the friction with the gas during the disk crossing.

We determine the true and the projected obliquity of HAT-P-36 and WASP-11/HAT-P-10 systems, respectively, which are both composed of a relatively cool star and a hot-Jupiter planet. Thanks to the high-resolution spectrograph HARPS-N, we observed the Rossiter-McLaughlin effect for both the systems by acquiring precise radial-velocity measurements during planetary transit events. We also present photometric observations comprising six light curves covering five transit events, obtained using three medium-class telescopes and the telescope-defocussing technique. One transit of WASP-11/HAT-P-10 was followed simultaneously from two observatories. The three transit light curves of HAT-P-36b show anomalies that are attributable to starspot complexes on the surface of the parent star, in agreement with the analysis of its spectra that indicate a moderate activity. By analysing the complete HATNet data set of HAT-P-36, we estimated the stellar rotation period by detecting a periodic photometric modulation in the light curve caused by star spots, obtaining Prot=15.3 days, which implies that the inclination of the stellar rotational axis with respect to the line of sight is 65 degree. We used the new spectroscopic and photometric data to revise the main physical parameters and measure the sky-projected misalignment angle of the two systems. We found \lambda=-14 degree for HAT-P-36 and \lambda=7 degree for WASP-11/HAT-P-10, indicating in both cases a good spin-orbit alignment. In the case of HAT-P-36, we also measured its real obliquity, which turned out to be 25 degrees.

The Trojan asteroids provide a unique perspective on the history of Solar System. As a large population of small bodies, they record important gravitational interactions and dynamical evolution of the Solar System. In the past decade, significant advances have been made in understanding physical properties, and there has been a revolution in thinking about the origin of Trojans. The ice and organics generally presumed to be a significant part of Trojan compositions have yet to be detected directly, though low density of the binary system Patroclus (and possibly low density of the binary/moonlet system Hektor) is consistent with an interior ice component. By contrast, fine-grained silicates that appear to be similar to cometary silicates in composition have been detected, and a color bimodality may indicate distinct compositional groups among the Trojans. Whereas Trojans had traditionally been thought to have formed near 5 AU, a new paradigm has developed in which the Trojans formed in the proto-Kuiper Belt, and they were scattered inward and captured in the Trojan swarms as a result of resonant interactions of the giant planets. Whereas the orbital and population distributions of current Trojans are consistent with this origin scenario, there are significant differences between current physical properties of Trojans and those of Kuiper Belt objects. These differences may be indicative of surface modification due to the inward migration of objects that became the Trojans, but understanding of appropriate modification mechanisms is poor and would benefit from additional laboratory studies. Many open questions remain, and the future promises significant strides in our understanding of Trojans. The time is ripe for a spacecraft mission to the Trojans, to turn these objects into geologic worlds that can be studied in detail to unravel their complex history.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.

Context. Astrometric monitoring of directly-imaged exoplanets allows the study of their orbital parameters and system architectures. Because most directly-imaged planets have long orbital periods (>20 AU), accurate astrometry is challenging when based on data acquired on timescales of a few years and usually with different instruments. The LMIRCam camera on the LBT is being used for the LEECH survey to search for and characterize young and adolescent exoplanets in L' band, including their system architectures. Aims. We first aim to provide a good astrometric calibration of LMIRCam. Then, we derive new astrometry, test the predictions of the orbital model of 8:4:2:1 mean motion resonance proposed by Go\'zdziewski & Migaszewski, and perform new orbital fitting of the HR 8799 bcde planets. We also present deep limits on a putative fifth planet interior to the known planets. Methods. We use observations of HR 8799 and the Theta1 Ori C field obtained during the same run in October 2013. Results. We first characterize the distortion of LMIRCam. We determine a platescale and a true north orientation for the images of 10.707 +/- 0.012 mas/pix and -0.430 +/- 0.076 deg, respectively. The errors on the platescale and true north orientation translate into astrometric accuracies at a separation of 1 of 1.1 mas and 1.3 mas, respectively. The measurements for all planets are usually in agreement within 3 sigma with the ephemeris predicted by Go\'zdziewski & Migaszewski. The orbital fitting based on the new astrometric measurements favors an architecture for the planetary system based on 8:4:2:1 mean motion resonance. The detection limits allow us to exclude a fifth planet slightly brighter/more massive than HR 8799 b at the location of the 2:1 resonance with HR 8799 e (~9.5 AU) and about twice as bright as HR 8799 cde at the location of the 3:1 resonance with HR 8799 e (~7.5 AU).

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude $K=51\pm3$ m s$^{-1}$. The derived mass of TrES-4b is found to be $M_{\rm p} = 0.49\pm0.04 \rm M_{Jup}$, significantly lower than previously reported. Combined with the large radius ($R_{\rm p} = 1.84_{-0.09}^{+0.08} \rm R_{Jup}$) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets.

XO-2 is the first confirmed wide stellar binary system where the almost twin components XO-2N and XO-2S have planets. This stimulated a detailed characterization study of the stellar and planetary components based on new observations. We collected high-resolution spectra with the HARPS-N spectrograph and multi-band light curves. Spectral analysis led to an accurate determination of the stellar atmospheric parameters and characterization of the stellar activity. We collected 14 transit light curves of XO-2Nb used to improve the transit parameters. Photometry provided accurate magnitude differences between the stars and a measure of their rotation periods. The iron abundance of XO-2N was found to be +0.054 dex greater, within more than 3-sigma, than that of XO-2S. We confirm a long-term variation in the radial velocities of XO-2N, and we detected a turn-over with respect to previous measurements. We suggest the presence of a second massive companion in an outer orbit or the stellar activity cycle as possible causes of the observed acceleration. The latter explanation seems more plausible with the present dataset. We obtained an accurate value of the projected spin-orbit angle for the XO-2N system (lambda=7+/-11 degrees), and estimated the real 3-D spin-orbit angle (psi=27 +12/-27 degrees). We measured the XO-2 rotation periods, and found a value of P=41.6 days in the case of XO-2N, in excellent agreement with the predictions. The period of XO-2S appears shorter, with an ambiguity between 26 and 34.5 days that we cannot solve with the present dataset alone. XO-2N appears to be more active than the companion, and this could be due to the fact that we sampled different phases of their activity cycle, or to an interaction between XO-2N and its hot Jupiter that we could not confirm.

Aims. With the purpose of determining the orbital parameters of exoplanetary systems from observational data, we have developed a software, named TRADES (TRAnsits and Dynamics of Exoplanetary Systems), to simultaneously fit observed radial velocities and transit times data. Methods. We implemented a dynamical simulator for N-body systems, which also fits the available data during the orbital integration and determines the best combination of the orbital parameters using grid search, $\chi^2$ minimization, genetic algorithms, particle swarm optimization, and bootstrap analysis. Results. To validate TRADES, we tested the code on a synthetic three-body system and on two real systems discovered by the Kepler mission: Kepler-9 and Kepler-11. These systems are good benchmarks to test multiple exoplanet systems showing transit time variations (TTVs) due to the gravitational interaction among planets. We have found that orbital parameters of Kepler-11 planets agree well with the values proposed in the discovery paper and with a a recent work from the same authors. We analyzed the first three quarters of Kepler-9 system and found parameters in partial agreement with discovery paper. Analyzing transit times (T0s) covering 12 quarters of Kepler data, that we have found a new best-fit solution. This solution outputs masses that are about 55% of the values proposed in the discovery paper; this leads to a reduced semi-amplitude of the radial velocities of about 12.80 m/s.

Direct imaging surveys for exoplanets commonly exclude binary stars from their target lists, leaving a large part of the overall planet demography unexplored. To address this gap in our understanding of planet formation and evolution, we have launched the first direct imaging survey dedicated to circumbinary planets: SPOTS, the Search for Planets Orbiting Two Stars. In this paper, we discuss the theoretical context, scientific merit, and technical feasibility of such observations, describe the target sample and observational strategy of our survey, and report on the first results from our pilot survey of 26 targets with the VLT NaCo facility. While we have not found any confirmed substellar companions to date, a number of promising candidate companions remain to be tested for common proper motion in upcoming follow-up observations. We also report on the astrometry of the three resolved binaries in our target sample. This pilot survey constitutes a successful proof of concept for our survey strategy and paves the way for a second stage of exploratory observations with VLT SPHERE.

Planet--planet scattering is a major dynamical mechanism able to significantly alter the architecture of a planetary system. In addition to that, it may also affect the formation and retention of a debris disk by the system. A violent chaotic evolution of the planets can easily clear leftover planetesimal belts preventing the ignition of a substantial collisional cascade that can give origin to a debris disk. On the other end, a mild evolution with limited steps in eccentricity and semimajor axis can trigger the formation of a debris disk by stirring an initially quiet planetesimal belt. The variety of possible effects that planet--planet scattering can have on the formation of debris disks is analysed and the statistical probability of the different outcomes is evaluated. This leads to the prediction that systems which underwent an episode of chaotic evolution might have a lower probability of harboring a debris disk.

We performed an intensive radial velocity monitoring of XO-2S, the wide companion of the transiting planet-host XO-2N, using HARPS-N at TNG in the framework of the GAPS programme. The radial velocity measurements indicate the presence of a new planetary system formed by a planet that is slightly more massive than Jupiter at 0.48 au and a Saturn-mass planet at 0.13 au. Both planetary orbits are moderately eccentric and were found to be dynamically stable. There are also indications of a long-term trend in the radial velocities. This is the first confirmed case of a wide binary whose components both host planets, one of which is transiting, which makes the XO-2 system a unique laboratory for understanding the diversity of planetary systems.

Stellar activity can be a source of radial velocity (RV) noise and can reproduce periodic RV variations similar to those produced by an exoplanet. We present the vigorous activity cycle in the primary of the visual binary HD200466, a system made of two almost identical solar-type stars with an apparent separation of 4.6 arcsec at a distance of 44+/-2 pc. High precision RV over more than a decade, adaptive optics (AO) images, and abundances have been obtained for both components. A linear trend in the RV is found for the secondary. We assumed that it is due to the binary orbit and once coupled with the astrometric data, it strongly constrains the orbital solution of the binary at high eccentricities (e~0.85) and quite small periastron of ~21 AU. If this orbital motion is subtracted from the primary radial velocity curve, a highly significant (false alarm probability <0.1%) period of about 1300 d is obtained, suggesting in a first analysis the presence of a giant planet, but it turned out to be due to the stellar activity cycle. Since our spectra do not include the Ca~II resonance lines, we measured a chromospheric activity indicator based on the Halpha line to study the correlation between activity cycles and long-term activity variations. While the bisector analysis of the line profile does not show a clear indication of activity, the correlation between the Halpha line indicator and the RV measurements identify the presence of a strong activity cycle.

The dynamics of systems of two and three planets, initially placed on circular and nearly coplanar orbits, is explored in the proximity of their stability limit. The evolution of a large number of systems is numerically computed and their dynamical behaviour is investigated with the frequency map analysis as chaos indicator. Following the guidance of this analysis, it is found that for two-planet systems the dependence of the Hill limit on the planet mass, usually made explicit through the Hill's radius parametrization, does not appear to be fully adequate. In addition, frequent cases of stable chaos are found in the proximity of the Hill limit. For three-planet systems, the usual approach adopted in numerical explorations of their stability, where the planets are initially separated by multiples of the mutual Hill radius, appears too reducing. A detailed sampling of the parameter space reveals that systems with more packed inner planets are stable well within previous estimates of the stability limit. This suggests that a two-dimensional approach is needed to outline when three-planet systems are dynamically stable.

Stellar flybys in star clusters are suspected to affect the orbital architecture of planetary systems causing eccentricity excitation and orbital misalignment between the planet orbit and the equatorial plane of the star. We explore whether the impulsive changes in the orbital elements of planets, caused by an hyperbolic stellar flyby, can be fully damped by the circumstellar disk surrounding the star. The time required to disperse stellar clusters is in fact comparable to circumstellar disk's lifetime. We have modelled in 3D a system made of a solar type star surrounded by a low density disk with a giant planet embedded in it approached on a hyperbolic encounter trajectory by a second star, of similar mass and with its own disk. We focus on extreme configurations where a very deep stellar flyby perturbs a Jovian planet on an external orbit. This allows to test in full the ability of the disk to erase the effects of the stellar encounter. We find that the amount of mass lost by the disk during the stellar flyby is less than in 2D models where a single disk was considered due to the mass exchange between the two disks at the encounter. The damping in eccentricity is slightly faster than in 2D models and it occurs on timescales of the order of a few kyr. The only trace of the flyby left in the planet system, after about 10^4 yr, is a small misalignment, lower than 9 degrees, between the star equatorial plane and the planet orbit. In a realistic model based on 3D simulations of star--planet--disk interactions, we find that stellar flybys cannot excite significant eccentricities and inclinations of planets in stellar clusters. The circumstellar disks hosting the planets damp on a short timescale all the step changes in the two orbital parameters produced during any stellar encounter. All records of past encounters are erased.

The size distribution of small asteroids in the Main Belt is assumed to be determined by an equilibrium between the creation of new bodies out of the impact debris of larger asteroids and the destruction of small asteroids by collisions with smaller projectiles. However, for a diameter less than 6 km we find that YORP-induced rotational disruption significantly contributes to the erosion even exceeding the effects of collisional fragmentation. Including this additional grinding mechanism in a collision evolution model for the asteroid belt, we generate size-frequency distributions from either an accretional (Weidenschilling, 2011) or an "Asteroids were born big" (Morbidelli, 2009) initial size-frequency distribution that are consistent with observations reported in Gladman et al. (2009). Rotational disruption is a new mechanism that must be included in all future collisional evolution models of asteroids.

Circumstellar disks in binaries are perturbed by the companion gravity causing significant alterations of the disk morphology. Spiral waves due to the companion tidal force also develop in the vertical direction and affect the disk temperature profile. These effects may significantly influence the process of planet formation. We perform 3D numerical simulations of disks in binaries with different initial dynamical configurations and physical parameters. Our goal is to investigate their evolution and their propensity to grow planets. We use an improved version of the SPH code VINE modified to better account for momentum and energy conservation. The energy equation includes a flux--limited radiative transfer algorithm and the disk cooling is obtained via "boundary particles". We model a system made of star/disk + star/disk where the secondary star (and relative disk) is less massive than the primary. The numerical simulations performed for different values of binary separation and disk density show that the disk morphology is substantially affected by the companion perturbations. Trailing spiral shock waves develop when the stars approach their pericenter. Strong hydraulic jumps occur at the shock front creating breaking waves and a consistent mass stream between the two disks, significantly heating them. The high gas temperature may prevent the ice condensation by moving outward the "snow line". The hydraulic jumps may slow down or even halt the dust coagulation process. At apocenter these perturbations are reduced and the disks are cooled down and less eccentric. The strength of the hydraulic jumps, disk heating, and mass exchange depends on the binary separation, and for larger semi-major axes, the tidal spiral pattern is substantially reduced.

A long-term multi-purpose observational programme has started with HARPS-N@TNG aimed to characterise the global architectural properties of exoplanetary systems. In this first paper we fully characterise the transiting system Qatar-1. We exploit HARPS-N high-precision radial velocity measurements obtained during a transit to measure the Rossiter-McLaughlin effect in the Qatar-1 system, and out-of-transit measurements to redetermine the spectroscopic orbit. New photometric transit light-curves are analysed and a spectroscopic characterisation of the host star atmospheric parameters is performed based on various methods (line equivalent width ratios, spectral synthesis, spectral energy distribution). We achieved a significant improvement in the accuracy of the orbital parameters and derived the spin-orbit alignment of the system; this information, combined with the spectroscopic determination of the host star properties, allows us to derive the fundamental physical parameters for star and planet (masses and radii). The orbital solution for the Qatar-1 system is consistent with a circular orbit and the system presents a sky-projected obliquity of lambda = -8.4+-7.1 deg. The planet, with a mass of 1.33+-0.05 M_J, is found to be significantly more massive than previously reported. The host star is confirmed to be metal-rich ([Fe/H]= 0.20+-0.10) and slowly rotating (vsinI = 1.7+-0.3 km/s), though moderately active, as indicated by strong chromospheric emission in the Ca II H&K line cores (logR'_HK about -4.60). The system is well aligned and fits well within the general lambda vs Teff trend. We definitely rule out any significant orbital eccentricity. The evolutionary status of the system is inferred based on gyrochronology, and the present orbital configuration and timescale for orbital decay are discussed in terms of star-planet tidal interactions.

Recent observations from NASA's Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.

In the context of the program Global Architecture of Planetary Systems (GAPS), we have performed radial velocity monitoring of the metal-poor star HIP 11952 on 35 nights over about 150 days using the newly installed high resolution spectrograph HARPS-N at the TNG and HARPS at ESO 3.6m telescope. The radial velocities show a scatter of 7 m/s, compatible with the measurement errors for such a moderately warm metal-poor star (Teff = 6040+-120K; [Fe/H] =-1.9+-0.1). We then exclude the presence of the two giant planets with periods of 6.95+-0.01 d and 290.0+-16.2 d and radial velocity semi-amplitudes of 100.3+-19.4 m/s and 105.2+-14.7 m/s, respectively, which had recently been announced. This result is important considering that HIP 11952 was thought to be the most metal-poor star hosting a planetary system with giant planets, thus challenging some models of planet formation.

Most stars form in embedded clusters. Stellar flybys may affect the orbital architecture of the systems by exciting the eccentricity and causing dynamical instability. Since, incidentally, the timescale over which a cluster loses its gaseous component and begins to disperse is comparable to the circumstellar disk lifetime, we expect that closer, and more perturbing, stellar flybys occur when the planets are still embedded in their birth disk. We investigate the effects of the disk on the dynamics of planets after the stellar encounter to test whether it can damp the eccentricity and return the planetary system to a non-excited state. We use the hydrodynamical code FARGO to study the disk+planet(s) system during and after the stellar encounter in the context of evolved disk models whose superficial density is 10 times lower than that of the Minimum Mass Solar Nebula. The numerical simulations show that the planet eccentricity, excited during a close stellar flyby, is damped on a short timescale (~ 10 Kyr) in spite of the disk low initial density and subsequent tidal truncation. This damping is effective also for a system of 3 giant planets and the effects of the dynamical instability induced by the passing star are quickly absorbed. If the circumstellar disk is still present around the star during a stellar flyby, a planet (or a planetary system) is returned to a non-excited state on a short timescale. This does not mean that stellar encounters do not affect the evolution of planets, but they do it in a subtle way with a short period of agitated dynamical evolution. At the end of it, the system resumes a quiet evolution and the planetary orbits are circularized by the interaction with the disk.

The outward migration of a pair of resonant-orbit planets, driven by tidal interactions with a gas-dominated disk, is studied in the context of evolved solar nebula models. The planets' masses, M1 and M2, correspond to those of Jupiter and Saturn. Hydrodynamical calculations in two and three dimensions are used to quantify the migration rates and analyze the conditions under which the outward migration mechanism may operate. The planets are taken to be fully formed after 1e6 and before 3e6 years. The orbital evolution of the planets in an evolving disk is then calculated until the disk's gas is completely dissipated. Orbital locking in the 3:2 mean motion resonance may lead to outward migration under appropriate conditions of disk viscosity and temperature. However, resonance locking does not necessarily result in outward migration. This is the case, for example, if convergent migration leads to locking in the 2:1 mean motion resonance, as post-formation disk conditions seem to suggest. Accretion of gas on the planets may deactivate the outward migration mechanism by raising the mass ratio M2/M1 and/or by reducing the accretion rate toward the star, hence depleting the inner disk. For migrating planets locked in the 3:2 mean motion resonance, there are stalling radii that depend on disk viscosity and on stellar irradiation, when it determines the disk's thermal balance. Planets locked in the 3:2 orbital resonance that start moving outward from within 1-2 AU may reach beyond 5 AU only under favorable conditions. However, within the explored space of disk parameters, only a small fraction - less than a few per cent - of the models predict that the interior planet reaches beyond 4 AU.

We have performed H and Ks band observations of the planetary system around HR 8799 using the new AO system at the Large Binocular Telescope and the PISCES Camera. The excellent instrument performance (Strehl ratios up to 80% in H band) enabled detection the inner planet HR8799e in the H band for the first time. The H and Ks magnitudes of HR8799e are similar to those of planets c and d, with planet e slightly brighter. Therefore, HR8799e is likely slightly more massive than c and d. We also explored possible orbital configurations and their orbital stability. We confirm that the orbits of planets b, c and e are consistent with being circular and coplanar; planet d should have either an orbital eccentricity of about 0.1 or be non-coplanar with respect to b and c. Planet e can not be in circular and coplanar orbit in a 4:2:1 mean motion resonances with c and d, while coplanar and circular orbits are allowed for a 5:2 resonance. The analysis of dynamical stability shows that the system is highly unstable or chaotic when planetary masses of about 5 MJup for b and 7 MJup for the other planets are adopted. Significant regions of dynamical stability for timescales of tens of Myr are found when adopting planetary masses of about 3.5, 5, 5, and 5 Mjup for HR 8799 b, c, d, and e respectively. These masses are below the current estimates based on the stellar age (30 Myr) and theoretical models of substellar objects.

Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the disc's radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the disc's longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($\pi$). The disc's radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.

The European Space Agency's Rosetta spacecraft passed by the main belt asteroid (21) Lutetia the 10th July 2010. With its ~100km size, Lutetia is one of the largest asteroids ever imaged by a spacecraft. During the flyby, the on-board OSIRIS imaging system acquired spectacular images of Lutetia's northern hemisphere revealing a complex surface scarred by numerous impact craters, reaching the maximum dimension of about 55km. In this paper, we assess the cratering history of the asteroid. For this purpose, we apply current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models, coupled with appropriate crater scaling laws, allow us to interpret the observed crater size-frequency distribution (SFD) and constrain the cratering history. Thanks to this approach, we derive the crater retention age of several regions on Lutetia, namely the time lapsed since their formation or global surface reset. We also investigate the influence of various factors -like Lutetia's bulk structure and crater obliteration- on the observed crater SFDs and the estimated surface ages. From our analysis, it emerges that Lutetia underwent a complex collisional evolution, involving major local resurfacing events till recent times. The difference in crater density between the youngest and oldest recognized units implies a difference in age of more than a factor of 10. The youngest unit (Beatica) has an estimated age of tens to hundreds of Myr, while the oldest one (Achaia) formed during a period when the bombardment of asteroids was more intense than the current one, presumably around 3.6Gyr ago or older.

As part of our radial velocity planet-search survey performed with SARG at TNG, we monitored the components of HD 132563 for ten years. It is a binary system formed by two rather similar solar type stars with a projected separation of 4.1 arcsec, which corresponds to 400 AU at the distance of 96 pc. The two components are moderately metal-poor and the age of the system is about 5 Gyr. We detected RV variations of HD 132563B with period of 1544 days and semi-amplitude of 26 m/s. From the star characteristics and line profile measurements, we infer their Keplerian origin. Therefore HD 132563B turns out to host a planet with a projected mass msini=1.49 MJup at 2.6 AU with a moderately eccentric orbit (e=0.22). The planet around HD 132563B is one of the few that are known in triple stellar systems, as we found that the primary HD 132563A is itself a spectroscopic binary with a period longer than 15 years and an eccentricity higher than 0.65. The spectroscopic component was not detected in adaptive-optics images taken with AdOpt@TNG, since it expected at a projected separation that was smaller than 0.2 arcsec at the time of our observations. A small excess in K band difference between the components with respect to the difference in V band is compatible with a companion of about 0.55 Msun. A preliminary statistical analysis of the occurrence of planets in triple systems indicate a similar frequency of planets around the isolated component in a triple system, components of wide binaries and single stars. There is no significant iron abundance difference between the components. The lack of stars in binary systems and open clusters showing strong enhancements of iron abundance, comparable to the typical metallicity difference between stars with and without giant planets, agrees with the idea that accretion of planetary material producing iron abundance anomalies larger than 0.1 dex is rare.

The flux of ISM neutral atoms surrounding stars and their environment affects the motion of dust particles in debris disks, causing a significant dynamical evolution. Large values of eccentricity and inclination can be excited and strong correlations settle in among the orbital angles. This dynamical behaviour, in particular for bound dust grains, can potentially cause significant asymmetries in dusty disks around solar type stars which might be detected by observations. However, the amount of orbital changes due to this non--gravitational perturbation is strongly limited by the collisional lifetime of dust particles. We show that for large values of the disk's optical depth the influence of ISM flow on the disk shape is almost negligible because the grains are collisionally destroyed before they can accumulate enough orbital changes due to the ISM perturbations. On the other hand, for values smaller than $10^{-3}$, peculiar asymmetric patterns appear in the density profile of the disk when we consider 1-10 mum grains, just above the blow-out threshold. The extent and relevance of these asymmetries grow for lower values of the optical depth. An additional sink mechanism, which may prevent the formation of large clumps and warping in the disks is related to the fast inward migration due to the drag component of the forces. When a significant eccentricity is pumped up by the ISM perturbations, the drag forces (Poynting-Robertson and in particular ISM drag) drive the disk particles on fast migrating tracks leading them into the star on a short timescale. It is then expected that disks with small optical depth expand inside the parent body ring all the way towards the star while disks with large optical depth would not significantly extend inside.

The dynamical features of the irregular satellites of the giant planets argue against an in-situ formation and are strongly suggestive of a capture origin. Since the last detailed investigations of their dynamics, the total number of satellites have doubled, increasing from 50 to 109, and almost tripled in the case of Saturn system. We have performed a new dynamical exploration of Saturn system to test whether the larger sample of bodies could improve our understanding of which dynamical features are primordial and which are the outcome of the secular evolution of the system. We have performed detailed N--Body simulations using the best orbital data available and analysed the frequencies of motion to search for resonances and other possible perturbing effects. We took advantage of the Hierarchical Jacobian Symplectic algorithm to include in the dynamical model of the system also the gravitational effects of the two outermost massive satellites, Titan and Iapetus. Our results suggest that Saturn's irregular satellites have been significantly altered and shaped by the gravitational perturbations of Jupiter, Titan, Iapetus and the Sun and by the collisional sweeping effect of Phoebe. In particular, the effects on the dynamical evolution of the system of the two massive satellites appear to be non-negligible. Jupiter perturbs the satellites through its direct gravitational pull and, indirectly, via the effects of the Great Inequality, i.e. its almost resonance with Saturn. Finally, by using the Hierarchical Clustering Method we found hints to the existence of collisional families and compared them with the available observational data.

The origin of the irregular satellites of the giant planets has been long debated since their discovery. Their dynamical features argue against an in-situ formation suggesting they are captured bodies, yet there is no global consensus on the physical process at the basis of their capture. In this paper we explore the collisional capture scenario, where the actual satellites originated from impacts occurred within Saturn's influence sphere. By modeling the inverse capture problem, we estimated the families of orbits of the possible parent bodies and the specific impulse needed for their capture. The orbits of these putative parent bodies are compared to those of the minor bodies of the outer Solar System to outline their possible region of formation. Finally, we tested the collisional capture hypothesis on Phoebe by taking advantage of the data supplied by Cassini on its major crater, Jason. Our results presented a realistic range of solutions matching the observational and dynamical data.

Debris disc analysis and modelling provide crucial information about the structure and the processes at play in extrasolar planetary systems. In binary systems, this issue is more complex because the disc should in addition respond to the companion star's perturbations. We explore the dynamical evolution of a collisionally active debris disc for different initial parent body populations, diverse binary configurations and optical depths. We focus on the radial extent and size distribution of the disc at a stationary state. We numerically follow the evolution of $10^{5}$ massless small grains, initially produced from a circumprimary disc of parent bodies following a size distribution in $dN \propto s^{-3.5}$ds . Grains are submitted to both stars' gravity as well as radiation pressure. In addition, particles are assigned an empirically derived collisional lifetime. For all the binary configurations the disc extends far beyond the critical semimajor axis $a_crit$ for orbital stability. This is due to the steady production of small grains, placed on eccentric orbits reaching beyond $a_crit$ by radiation pressure. The amount of matter beyond acrit depends on the balance between collisional production and dynamical removal rates: it increases for more massive discs as well as for eccentric binaries. Another important effect is that, in the dynamically stable region, the disc is depleted from its smallest grains. Both results could lead to observable signatures. We have shown that a companion star can never fully truncate a collisionally active disc. For eccentric companions, grains in the unstable regions can significantly contribute to the thermal emission in the mid-IR. Discs with sharp outer edges, especially bright ones such as HR4796A, are probably shaped by other mechanisms.

Hydrodynamical simulations of two giant planets embedded in a gaseous disk have shown that in case of a smooth convergent migration they end up trapped into a mean motion resonance. These findings have led to the conviction that the onset of dynamical instability causing close encounters between the planets can occur only after the dissipation of the gas when the eccentricity damping is over. We show that a system of three giant planets may undergo planet-planet scattering when the gaseous disk, with density values comparable to that of the Minimum Mass Solar Nebula, is still interacting with the planets. The hydrodynamical code FARGO--2D--1D is used to model the evolution ofthe disk and planets, modified to properly handle close encounters between the massive bodies. Our simulations predict a variety of different outcomes of the scattering phase which includes orbital exchange, planet merging and scattering of a planet in a hyperbolic orbit. This implies thatthe final fate of a multiplanet system under the action of the disk torques is not necessarily a packed resonant configuration.

The Eurybates family is a compact core inside the Menelaus clan, located in the L4 swarm of Jupiter Trojans. Fornasier et al. (2007) found that this family exhibits a peculiar abundance of spectrally flat objects, similar to Chiron-like Centaurs and C-type main belt asteroids. On the basis of the visible spectra available in literature, Eurybates family's members seemed to be good candidates for having on their surfaces water/water ice or aqueous altered materials. To improve our knowledge of the surface composition of this peculiar family, we carried out an observational campaign at the Telescopio Nazionale Galileo (TNG), obtaining near-infrared spectra of 7 members. Our data show a surprisingly absence of any spectral feature referable to the presence of water, ices or aqueous altered materials on the surface of the observed objects. Models of the surface composition are attempted, evidencing that amorphous carbon seems to dominate the surface composition of the observed bodies and some amount of silicates (olivine) could be present.

The cratering history of main belt asteroid (2867) Steins has been investigated using OSIRIS imagery acquired during the Rosetta flyby that took place on the 5th of September 2008. For this purpose, we applied current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models coupled with appropriate crater scaling laws, allow the cratering history to be estimated. Hence, we derive Steins' cratering retention age, namely the time lapsed since its formation or global surface reset. We also investigate the influence of various factors -like bulk structure and crater erasing- on the estimated age, which spans from a few hundred Myrs to more than 1Gyr, depending on the adopted scaling law and asteroid physical parameters. Moreover, a marked lack of craters smaller than about 0.6km has been found and interpreted as a result of a peculiar evolution of Steins cratering record, possibly related either to the formation of the 2.1km wide impact crater near the south pole or to YORP reshaping.

The analysis of non-radiative sources of static or time-dependent gravitational fields in the Solar System is crucial to accurately estimate the free-fall orbits of the LISA space mission. In particular, we take into account the gravitational effects of Interplanetary Dust (ID) on the spacecraft trajectories. The perturbing gravitational field has been calculated for some ID density distributions that fit the observed zodiacal light. Then we integrated the Gauss planetary equations to get the deviations from the LISA keplerian orbits around the Sun. This analysis can be eventually extended to Local Dark Matter (LDM), as gravitational fields are expected to be similar for ID and LDM distributions. Under some strong assumptions on the displacement noise at very low frequency, the Doppler data collected during the whole LISA mission could provide upper limits on ID and LDM densities.

We investigate the dynamical evolution of a Jovian--mass planet injected into an orbit highly inclined with respect to its nesting gaseous disk. Planet--planet scattering induced by convergent planetary migration and mean motion resonances may push a planet into such an out of plane configuration with inclinations as large as $20^\circ-30^\circ$. In this scenario the tidal interaction of the planet with the disk is more complex and, in addition to the usual Lindblad and corotation resonances, it involves also inclination resonances responsible of bending waves. We have performed three--dimensional hydrodynamic simulations of the disk and of its interactions with the planet with a Smoothed Particle Hydrodynamics (SPH) code. A main result is that the initial large eccentricity and inclination of the planetary orbit are rapidly damped on a timescale of the order of $10^3$ yrs, almost independently of the initial semimajor axis and eccentricity of the planet. The disk is warped in response to the planet perturbations and it precesses. Inward migration occurs also when the planet is inclined and it has a drift rate which is intermediate between type I and type II migration. The planet is not able to open a gap until its inclination becomes lower than $\sim 10^\circ$ when it also begins to accrete a significant amount of mass from the disk.

We study the evolution of circumstellar massive disks around the primary star of a binary system focusing on the computation of disk eccentricity. In particular, we concentrate on its dependence on the binary eccentricity. Self-gravity is included in our numerical simulations. Our standard model assumes a semimajor axis for the binary of 30 AU, the most probable value according to the present binary statistics.

We explore planet formation in binary systems around the central star where the protoplanetary disk plane is highly inclined with respect to the companion star orbit. This might be the most frequent scenario for binary separations larger than 40 AU, according to Hale (1994). We focus on planetesimal accretion and compute average impact velocities in the habitable region and up to 6 AU from the primary.

Recent studies have shown that alpha Centauri B might be, from an observational point of view, an ideal candidate for the detection of an Earth-like planet in or near its habitable zone (0.5-0.9AU). We study here if such habitable planets can form, by numerically investigating the planet-formation stage which is probably the most sensitive to binarity effects: the mutual accretion of km-sized planetesimals. Using a state-of-the-art algorithm for computing the impact velocities within a test planetesimal population, we find that planetesimal growth is only possible, although marginally, in the innermost part of the HZ around 0.5AU. Beyond this point, the combination of secular perturbations by the binary companion and gas drag drive the mutual velocities beyond the erosion limit. Impact velocities might later decrease during the gas removal phase, but this probably happens too late for preventing most km-sized objects to be removed by inward drift, thus preventing accretion from starting anew. A more promising hypothesis is that the binary formed in a crowded cluster, where it might have been wider in its initial stages, when planetary formation was ongoing. We explore this scenario and find that a starting separation roughly 15 AU wider, or an eccentricity 2.5 times lower than the present ones are required to have an accretion-friendly environment in the whole HZ.

The five planets discovered around the main-sequence star 55 Cnc may represent a case of stable chaos. By using both the Frequency Map Analysis and MEGNO we find that about 15 % of the systems that can be build from the nominal orbital elements of the system are highly chaotic. However, in spite of the fast diffusion rate in the phase space, the planetary system is not destabilized over 400 Myr and close encounters between the planets are avoided.

We numerically explore planet formation around alpha Cen A by focusing on the crucial planetesimals-to-embryos phase. Our code computes the relative velocity distribution, and thus the accretion vs. fragmentation trend, of planetesimal populations having any given size distribution. This is a critical aspect of planet formation in binaries since the pericenter alignment of planetesimal orbits due to the gravitational perturbations of the companion star and to gas friction strongly depends on size. We find that, for the nominal case of a MMSN gas disc, the region beyond 0.5AU from the primary is hostile to planetesimal accretion. In this area, impact velocities between different-size bodies are increased, by the differential orbital phasing, to values too high to allow mutual accretion. For any realistic size distribution for the planetesimal population, this accretion-inhibiting effect is the dominant collision outcome and the accretion process is halted. Results are robust with respect to the profile and density of the gas disc: except for an unrealistic almost gas-free case, the inner accretion safe area never extends beyond 0.75AU. We conclude that planet formation is very difficult in the terrestrial region around alpha Cen A, unless it started from fast-formed very large (>30km) planetesimals. Notwithstanding these unlikely initial conditions, the only possible explanation for the presence of planets around 1 AU from the star would be the hypothetical outward migration of planets formed closer to the star or a different orbital configuration in the binary's early history. Our conclusions differ from those of several studies focusing on the later embryos-to-planets stage, confirming that the planetesimals-to-embryos phase is more affected by binary perturbations.

Numerical simulations of planetesimal accretion in circumprimary and circumbinary orbits are described. The secular perturbations by the com- panion star and gas drag are included in our models. We derive limits on the parameters of the binary system for which accretion and then planetary forma- tion are possible. In the circumbinary case we also outline the radial distance from the baricenter of the stars beyond which accumulation always occurs. Hy- drodynamical simulations are also presented to validate our N-body approach based on the axisymmetric approximation for the gas of the disk.

We study the dynamics of Jupiter Trojans in the early phase of the Solar system while the outer planets migrated due to their interaction with the planetesimal disk.

This paper explains why planets in binary star systems might have a lower frequency. A transient triple state of the binary causes the dispersal of planets.

We numerically investigate the possibility of planetesimal accretion in circumbinary disks, under the coupled influence of both stars' secular perturbations and friction due to the gaseous component of the protoplanetary disk. We focus on one crucial parameter: the distribution of encounter velocities between planetesimals in the 0.5 to 100km size range. An extended range of binary systems with differing orbital parameters is explored. The resulting encounter velocities are compared to the threshold velocities below which the net outcome of a collision is accumulation into a larger body instead of mass erosion. For each binary configuration, we derive the critical radial distance from the binary barycenter beyond which planetesimal accretion is possible. This critical radial distance is smallest for equal-mass binaries on almost circular orbits. It shifts to larger values for increasing eccentricities and decreasing mass ratio. The importance of the planetesimals' orbital alignments of planetesimals due to gas drag effects is discussed.

We study to which extent planetesimal accretion is affected by the perturbing presence of a compagnon star. We concentrate on one crucial parameter: the distribution of encounter velocities within the planetesimal swarm. We numerically explore the evolution of this parameter taking into account the secular perturbations of the binary and friction due to the very likely presence of gas in the disk. $<\Delta v>$ maps are derived, for planetesimals of different sizes, for a total of 120 binary configurations (eccentricity eb and separation ab). We identify for each case 3 different accreting behaviours. 1) In regions where no significant dV increase is observed, ``classical'' (i.e., single-star) accretion is possible. 2) In regions where $dV>v_{ero}$, the threshold velocity above which all impacts are eroding, no accretion is possible and planet growth is stopped. 3) In between these 2 limiting behaviours, a large fraction of binary configurations leads to significant dV increase, but still below the erosion threshold. In this intermediate case, planetesimal growth can occur, but proceeds slower than in the single-star case, possibly following the so-called type II runaway groth mode.

We present the results of a visible spectroscopic and photometric survey of Jupiter Trojans belonging to different dynamical families carried out at the ESO-NTT telescope. We obtained data on 47 objects, 23 belonging to the L5 swarm and 24 to the L4 one. These data together with those already published by Fornasier et al. (2004a) and Dotto et al. (2006), constitute a total sample of visible spectra for 80 objects. The survey allows us to investigate six families (Aneas, Anchises, Misenus, Phereclos, Sarpedon, Panthoos) in the L5 cloud and four L4 families (Eurybates, Menelaus, 1986 WD and 1986 TS6). The sample that we measured is dominated by D--type asteroids, with the exception of the Eurybates family in the L4 swarm, where there is a dominance of C- and P-type asteroids. All the spectra that we obtained are featureless with the exception of some Eurybates members, where a drop--off of the reflectance is detected shortward of 5200 A. Similar features are seen in main belt C-type asteroids and commonly attributed to the intervalence charge transfer transition in oxidized iron. Our sample comprises fainter and smaller Trojans as compared to the literature's data and allows us to investigate the properties of objects with estimated diameter smaller than 40--50 km. The analysis of the spectral slopes and colors versus the estimated diameters shows that the blue and red objects have indistinguishable size distribution. We perform a statistical investigation of the Trojans's spectra property distributions as a function of their orbital and physical parameters, and in comparison with other classes of minor bodies in the outer Solar System. Trojans at lower inclination appear significantly bluer than those at higher inclination, but this effect is strongly driven by the Eurybates family.

The present dynamical configuration of planets in binary star systems may not reflect their formation process since the binary orbit may have changed in the past after the planet formation process was completed. An observed binary system may have been part of a former hierarchical triple that became unstable after the planets completed their growth around the primary star. Alternatively, in a dense stellar environment even a single stellar encounter between the star pair and a singleton may singificantly alter the binary orbit. In both cases the planets we observe at present would have formed when the dynamical environment was different from the presently observed one. We have numerically integrated the trajectories of the stars (binary plus singleton) and of test planets to investigate the abovementioned mechanisms. Our simulations show that the circumstellar environment during planetary formation around the primary was gravitationally less perturbed when the binary was part of a hierarchical triple because the binary was necessarely wider and, possibly, less eccentric. This circumstance has consequences for the planetary system in terms of orbital spacing, eccentricity, and mass of the individual planets. Even in the case of a single stellar encounter the present appearance of a planetary system in a binary may significantly differ from what it had while planet formation was ongoing. However, while in the case of instability of a triple the trend is always towards a tighter and more eccentric binary system, when a single stellar encounter affects the system the orbit of the binary can become wider and be circularized.

We investigate classical planetesimal accretion in a binary star system of separation ab<50AU by numerical simulations, with particular focus on the region at a distance of 1 AU from the primary. The planetesimals orbit the primary, are perturbed by the companion and are in addition subjected to a gas drag force. We concentrate on the problem of relative velocities dv among planetesimals of different sizes. For various stellar mass ratios and binary orbital parameters we determine regions where dv exceed planetesimal escape velocities v_esc (thus preventing runaway accretion) or even the threshold velocity v_ero for which erosion dominates accretion. Gaseous friction has two crucial effects on the velocity distribution: it damps secular perturbations by forcing periastron alignment of orbits, but at the same time the size--dependence of this orbital alignment induces a significant dv increase between bodies of different sizes. This differential phasing effect proves very efficient and almost always increases dv to values preventing runaway accretion, except in a narrow domain of almost circular companion orbits. The erosion threshold dv>v_ero is reached in a wide (ab,eb) space for small (<10km) planetesimals, but in a much more limited region for bigger ~50km objects. In the intermediate v_esc<dv < v_ero domain, a possible growth mode would be the type II runaway growth identified by Kortenkmap et al.(2001)

We numerically investigate under which conditions the planet detected at 2.1 AU of Gamma-Cephei could form through the core-accretion scenario despite the perturbing presence of the highly eccentric companion star. We first show that the initial stage of runaway accretion of kilometer-sized planetesimals is possible within 2.5 AU from the central star only if large amounts of gas are present. In this case, gaseous friction induces periastron alignment of the orbits which reduces the otherwise high mutual impact velocities due to the companion's secular perturbations. The following stage of mutual accretion of large embryos is also modeled. According to our simulations, the giant impacts among the embryos always lead to a core of 10 Mearth within 10 Myr, the average lifetime of gaseous discs. However, the core always ends up within 1.5 AU from the central star. Either the core grows more quickly in the inner region of the disc, or it migrates inside by scattering the residual embryos

The star HD 219542 B has been reported by us (Desidera et al. 2003) to show low-amplitude radial velocity variations that could be due to the presence of a Saturn-mass planetary companion or to stellar activity phenomena. In this letter we present the results of the continuation of the radial velocity monitoring as well as a discussion of literature determinations of the chromospheric activity of the star (Wright et al. 2004). These new data indicate that the observed radial velocity variations are likely related to stellar activity. In particular, there are indications that HD 219542 B underwent a phase of enhanced stellar activity in 2002 while the activity level has been lower in both 2001 and 2003. Our 2003 radial velocity measurements now deviate from our preliminary orbital solution and the peak in the power spectrum at the proposed planet period is severely reduced by the inclusion of the new data. We therefore dismiss the planet hypothesis as the cause of the radial velocity variations.

In this paper we present iron abundance analysis for 23 wide binaries with main sequence components in the temperture range 4900-6300 K, taken from the sample of the pairs currently included in the radial velocity planet search on going at the TNG. The use of a line-by-line differential analysis technique between the components of each pair allows to reach errors of about 0.02 dex in the iron content difference. Most of the pairs have abundance differences lower than 0.02 dex and there are no pairs with differences larger than 0.07 dex. The four cases of differences larger than 0.02 dex may be spouriuos because of the larger error bars affecting pairs with large temperature difference, cold stars and rotating stars. The pair HD 219542, previously reported by us to have different composition, results instead normal. For non-rotating stars warmer than 5500 K, we are able to exclude in most cases the ingestion of more than 1 Earth Mass of iron (about 5 Earth masses of meteoritic material) during the main sequence lifetime of the stars, placing more stringent limits (about 0.4 Earth masses of iron) in five cases of warm stars. This latter limit is similar to the estimates of rocky material accreted by the Sun during its main sequence lifetime. Combining the results of the present analysis with those for the Hyades and Pleiades, we conclude that the hypothesis that pollution by planetary material is the only responsible for the highest metallicity of the stars with planets may be rejected at more than 99% level of confidence if the incidence of planets in these samples is as high as 8% and similar to the field stars included in current radial velocity surveys. However, the significance of this result drops considerably if the incidence of planets around stars in binary systems and clusters is less than a half than around normal field stars.

The components of the wide binary HD 219542 were recently found to differ in metallicity by about 0.1 dex (Gratton et al. A&A 377, 123). In this paper, we present the results of 2 years of high precision radial velocity monitoring of these stars performed at the Telecopio Nazionale Galileo (TNG) using the high resolution spectrograph SARG. No indication for radial velocity variations above the measurement errors (5 m/s) was found for the metal richer component A. This allows us to place upper mass-limits for planets around this star. HD 219542B instead shows a low amplitude variation with a 112 day period at a confidence level of about 96-97%. This might suggest the presence of a Saturn-mass planet, although it is still possible that these variations are due to moderate activity of the star. Tests based on variations of bisectors, stellar magnitude and line equivalent widths were inconclusive so far.

The Planck mission, originally devised for cosmological studies, offers the opportunity to observe Solar System objects at millimetric and submillimetric wavelengths. We concentrate in this paper on the asteroids of the Main Belt. We intend to estimate the number of asteroids that can can be detected during the mission and to evaluate the strength of their signal. We have rescaled the instrument sensitivities, calculated by the LFI and HFI teams for sources fixed in the sky, introducing some degradation factors to properly account for moving objects. In this way a detection threshold is derived for asteroidal detection that is related to the diameter of the asteroid and its geocentric distance. We have developed a numerical code that models the detection of asteroids in the LFI and HFI channels during the mission. This code perfoprms a detailed integration of the orbits of the asteroids in the timespan of the mission and identifies those bodies that fall in the beams of Planck and their signal stenght. According to our simulations, a total of 397 objects will be observed by Planck and an asteroidal body will be detected in some beam in 30% of the total sky scan--circles. A significant fraction (in the range from ~50 to 100 objects) of the 397 asteroids will be observed with a high S/N ratio. Flux measurements of a large sample of asteroids in the submillimeter and millimeter range are relevant since they allow to analyze the thermal emission and its relation to the surface and regolith properties. Furthermore, it will be possible to check on a wider base the two standard thermal models, based on a nonrotating or rapidly rotating sphere. Our method can also be used to separate Solar System sources from cosmological sources in the survey. This work is based on Planck LFI activities.

At present the possible existence of planets around the stars of a close binary system is still matter of debate. Can planetary bodies form in spite of the strong gravitational perturbations of the companion star? We study in this paper via numerical simulation the last stage of planetary formation, from embryos to terrestrial planets in the Alpha Cen system, the prototype of close binary systems. We find that Earth class planets can grow around Alpha Cen A on a time-scale of 50 Myr. In some of our numerical models the planets form directly in the habitable zone of the star in low eccentric orbits. In one simulation two of the final planets are in a 2:1 mean motion resonance that, however, becomes unstable after 200 Myr. During the formation process some planetary embryos fall into the stars possibly altering their metallicity.