Source author record

F. Marzari

F. Marzari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

41works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

41 published item(s)

preprint2022arXiv

The GAPS Programme with HARPS-N at TNG. XXXV. Fundamental properties of transiting exoplanet host stars

Exoplanetary properties depend on stellar properties: to know the planet with accuracy and precision it is necessary to know the star as accurately and precisely as possible. Our immediate aim is to characterize in a homogeneous and accurate way a sample of 27 transiting planet-hosting stars observed within the GAPS program. We determined stellar parameters (effective temperature, surface gravity, rotational velocity) and abundances of 26 elements (Li,C,N,O,Na,Mg,Al,Si,S,Ca,Sc,Ti,V,Cr,Fe,Mn,Co,Ni,Cu,Zn,Y,Zr,Ba,La,Nd,Eu). Our study is based on high-resolution HARPS-N@TNG and FEROS@ESO spectra and uniform techniques. We derived kinematic properties from Gaia data and estimated for the first time in exoplanet host stars ages using elemental ratios as chemical clocks. Teff of our stars is of 4400-6700 K, while [Fe/H] is within -0.3 and 0.4 dex. Lithium is present in 7 stars. [X/H] and [X/Fe] abundances vs [Fe/H] are consistent with the Galactic Chemical Evolution. The dependence of [X/Fe] with the condensation temperature is critically analyzed with respect to stellar and kinematic properties. All targets with measured C and O abundances show C/O<0.8, compatible with Si present in rock-forming minerals. Most of targets show 1.0<Mg/Si<1.5, compatible with Mg distributed between olivine and pyroxene. HAT-P-26, the target hosting the lowest-mass planet, shows the highest Mg/Si ratio. From our chemo-dinamical analysis we find agreement between ages and position within the Galactic disk. We note a tendency for higher density planets to be around metal-rich stars and hints of higher stellar abundances of some volatiles for lower mass planets. We cannot exclude that part of our results could be also related to the location of the stars within the Galactic disk. We trace the planetary migration scenario from the composition of the planets related to the chemical composition of the hosting stars

preprint2020arXiv

The GAPS Programme at TNG -- XXIII. HD 164922 d: a close-in super-Earth discovered with HARPS-N in a system with a long-period Saturn mass companion

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.

preprint2020arXiv

The size, shape, density and ring of the dwarf planet Haumea from a stellar occultation

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-$σ$ limits), respectively, is detected.

preprint2016arXiv

Acceleration of Individual, Decimetre-sized Aggregates in the Lower Coma of Comet 67P/Churyumov-Gerasimenko

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.

preprint2016arXiv

Changes in the physical environment of the inner coma of 67P/Churyumov-Gerasimenko with decreasing heliocentric distance

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.

preprint2016arXiv

Characterizing HR3549B using SPHERE

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.

preprint2016arXiv

Decimetre-scaled spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

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°-33° 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.

preprint2016arXiv

Dynamics of Circumstellar Disks III: The case of GG Tau A

(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.

preprint2016arXiv

Stability of multiplanet systems in binaries

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.

preprint2016arXiv

Sublimation of icy aggregates in the coma of comet 67P/Churyumov-Gerasimenko detected with the OSIRIS cameras onboard Rosetta

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.

preprint2016arXiv

Summer fireworks on comet 67P

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.

preprint2016arXiv

The dust environment of comet 67P/Churyumov-Gerasimenko from Rosetta OSIRIS and VLT observations in the 4.5 to 2.9 au heliocentric distance range inbound

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 $α$=--3, instead of $α$=--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.

preprint2016arXiv

The GAPS programme with HARPS-N at TNG XI. Pr~0211 in M~44: the first multi-planet system in an open cluster

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)

preprint2015arXiv

Are fractured cliffs the source of cometary dust jets ? Insights from OSIRIS/Rosetta at 67P

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.

preprint2015arXiv

Comet 67P/Churyumov-Gerasimenko: Constraints on its origin from OSIRIS observations

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.

preprint2015arXiv

Gravitational slopes, geomorphology, and material strengths of the nucleus of comet 67P/Churyumov-Gerasimenko from OSIRIS observations

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.

preprint2015arXiv

Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft

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.

preprint2015arXiv

Temporal morphological changes in the Imhotep region of comet 67P/Churyumov-Gerasimenko

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.

preprint2015arXiv

The GAPS Programme with HARPS-N at TNG VIII: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting planetary systems HAT-P-36 and WASP-11/HAT-P-10

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 λ=-14 degree for HAT-P-36 and λ=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.

preprint2015arXiv

The GAPS Programme with HARPS-N at TNG. X. Differential abundances in the XO-2 planet hosting binary

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σ$ 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.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG V. A comprehensive analysis of the XO-2 stellar and planetary systems

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.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG VI: The Curious Case of TrES-4b

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.

preprint2015arXiv

The LEECH Exoplanet Imaging Survey. Further constraints on the planet architecture of the HR 8799 system

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ździewski & 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ździewski & 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).

preprint2014arXiv

A vigorous activity cycle mimicking a planetary system in HD200466

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.

preprint2014arXiv

Dynamical behaviour of multiplanet systems close to their stability limit

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.

preprint2014arXiv

The GAPS programme with HARPS-N@TNG IV: A planetary system around XO-2S

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.

preprint2013arXiv

The GAPS programme with HARPS-N at TNG. I: Observations of the Rossiter-McLaughlin effect and characterisation of the transiting system Qatar-1

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.

preprint2013arXiv

The GAPS Programme with HARPS-N@TNG II: No giant planets around the metal-poor star HIP 11952

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.

preprint2012arXiv

Eccentricity of radiative discs in close binary-star systems

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 ($π$). 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.

preprint2012arXiv

LBT observations of the HR 8799 planetary system: First detection of HR8799e in H band

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.

preprint2011arXiv

A giant planet in the triple system HD132563

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.

preprint2011arXiv

The Cratering History of Asteroid (21) Lutetia

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.

preprint2010arXiv

A new perspective on the irregular satellites of Saturn - I Dynamical and collisional history

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.

preprint2010arXiv

A new perspective on the irregular satellites of Saturn - II Dynamical and physical origin

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.

preprint2010arXiv

A Peculiar Family of Jupiter Trojans: the Eurybates

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.

preprint2010arXiv

Modulation of LISA free-fall orbits due to the Earth-Moon system

We calculate the effect of the Earth-Moon (EM) system on the free-fall motion of LISA test masses. We show that the periodic gravitational pulling of the EM system induces a resonance with fundamental frequency 1 yr^-1 and a series of periodic perturbations with frequencies equal to integer harmonics of the synodic month (9.92 10^-7 Hz). We then evaluate the effects of these perturbations (up to the 6th harmonics) on the relative motions between each test masses couple, finding that they range between 3mm and 10pm for the 2nd and 6th harmonic, respectively. If we take the LISA sensitivity curve, as extrapolated down to 10^-6 Hz, we obtain that a few harmonics of the EM system can be detected in the Doppler data collected by the LISA space mission. This suggests that the EM system gravitational near field could provide an absolute calibration for the LISA sensitivity at very low frequencies.

preprint2010arXiv

Planet--planet scattering in circumstellar gas disks

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.

preprint2010arXiv

The Cratering History of Asteroid (2867) Steins

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.

preprint2009arXiv

Interaction of a giant planet in an inclined orbit with a circum-stellar disk

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.

preprint2009arXiv

On the eccentricity of self-gravitating circumstellar disks in eccentric binary systems

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.

preprint2009arXiv

Planet formation in highly inclined binaries

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.