We present the results from four stellar occultations by (486958) Arrokoth, the flyby target of the New Horizons extended mission. Three of the four efforts led to positive detections of the body, and all constrained the presence of rings and other debris, finding none. Twenty-five mobile stations were deployed for 2017 June 3 and augmented by fixed telescopes. There were no positive detections from this effort. The event on 2017 July 10 was observed by SOFIA with one very short chord. Twenty-four deployed stations on 2017 July 17 resulted in five chords that clearly showed a complicated shape consistent with a contact binary with rough dimensions of 20 by 30 km for the overall outline. A visible albedo of 10% was derived from these data. Twenty-two systems were deployed for the fourth event on 2018 Aug 4 and resulted in two chords. The combination of the occultation data and the flyby results provides a significant refinement of the rotation period, now estimated to be 15.9380 $\pm$ 0.0005 hours. The occultation data also provided high-precision astrometric constraints on the position of the object that were crucial for supporting the navigation for the New Horizons flyby. This work demonstrates an effective method for obtaining detailed size and shape information and probing for rings and dust on distant Kuiper Belt objects as well as being an important source of positional data that can aid in spacecraft navigation that is particularly useful for small and distant bodies.

We present results from three world-wide campaigns that resulted in the detections of two single-chord and one multi-chord stellar occultations by the Plutino object (84922) 2003~VS$_2$. From the single-chord occultations in 2013 and 2014 we obtained accurate astrometric positions for the object, while from the multi-chord occultation on November 7th, 2014, we obtained the parameters of the best-fitting ellipse to the limb of the body at the time of occultation. We also obtained short-term photometry data for the body in order to derive its rotational phase during the occultation. The rotational light curve present a peak-to-peak amplitude of 0.141 $\pm$ 0.009 mag. This allows us to reconstruct the three-dimensional shape of the body, with principal semi-axes $a = 313.8 \pm 7.1$ km, $b = 265.5^{+8.8}_{-9.8}$ km, and $c = 247.3^{+26.6}_{-43.6}$ km, which is not consistent with a Jacobi triaxial equilibrium figure. The derived spherical volume equivalent diameter of $548.3 ^{+29.5}_{-44.6}$ km is about 5\% larger than the radiometric diameter of 2003~VS$_2$ derived from Herschel data of $523 \pm 35$ km, but still compatible with it within error bars. From those results we can also derive the geometric albedo ($0.123 ^{+0.015}_{-0.014}$) and, under the assumption that the object is a Maclaurin spheroid, the density $\rho = 1400^{+1000}_{-300}$ for the plutino. The disappearances and reappearances of the star during the occultations do not show any compelling evidence for a global atmosphere considering a pressure upper limit of about 1 microbar for a pure nitrogen atmosphere, nor secondary features (e.g. rings or satellite) around the main body.

The study of the visible colours of the trans-Neptunian objects opened a discussion almost 20 years ago which, in spite of the increase in the amount of available data, seems far from subside. Visible colours impose constraints to the current theories of the early dynamical evolution of the Solar System such as the environment of formation, initial surface composition, and how (if) they were scattered to regions closer to the inner planets. In this paper we present an updated version of our database of absolute colours and relative phase coefficients for 117 objects. We define the absolute colours as the difference of the absolute magnitudes $H_V-H_R$, and the relative phase coefficient as the difference of the slopes of the phase curves $\Delta\beta$. These were obtained joining our own observations plus data from the literature. The methodology has been introduced in previous works and here we expand in some interesting results, in particular the strong anti-correlation found between $H_V-H_R$ and $\Delta\beta$, which means that redder objects have steeper phase curves in the R filter, while bluer objects have steeper phase curves in the V filter. We analyse a series of results published in the literature in view of our database, which is free of phase effects, and show that their statistical meaning is not very strong. We point out that phase-colouring and observational errors play an important role in the understanding of these proposed relationships.

Context. Pluto's tenuous nitrogen (N2) atmosphere undergoes strong seasonal effects due to high obliquity and orbital eccentricity, and has been recently (July 2015) observed by the New Horizons spacecraft. Goals are (i) construct a well calibrated record of the seasonal evolution of surface pressure on Pluto and (ii) constrain the structure of the lower atmosphere using a central flash observed in 2015. Method: eleven stellar occultations by Pluto observed between 2002 and 2016 are used to retrieve atmospheric profiles (density, pressure, temperature) between $\sim$5 km and $\sim$380 km altitude levels (i.e. pressures from about 10 microbar to 10 nanobar). Results: (i) Pressure has suffered a monotonic increase from 1988 to 2016, that is compared to a seasonal volatile transport model, from which tight constraints on a combination of albedo and emissivity of N2 ice are derived; (ii) A central flash observed on 2015 June 29 is consistent with New Horizons REX profiles, provided that (a) large diurnal temperature variations (not expected by current models) occur over Sputnik Planitia and/or (b) hazes with tangential optical depth of about 0.3 are present at 4-7 km altitude levels and/or (c) the nominal REX density values are overestimated by an implausibly large factor of about 20% and/or (d) higher terrains block part of the flash in the Charon facing hemisphere.

Aims. The aim of this work is to study the dynamical effects of the Galaxy on binary star systems with physical and orbital charac- teristics similar to those of the population of known wide binary stars with exoplanets. As secondary goal we analyse the possible consequences on the stability of a hypothetical planetary system orbiting one of the stellar components. Methods. We numerically reproduced the temporal evolution of a sample of 3 * 10 5 binary star systems disturbed by the Galactic potential and passing stars in an environment similar to the solar neighbourhood. Results. Our results show that the dynamical evolution of the population of wide binary stars with exoplanets in the solar neigh- bourhood is modelled by the process of disruption of binary star systems induced by the Galaxy. We found that this process depends mainly on the separation between both stars, whereas it is almost independent of the initial orbital configuration. Moreover, our calcu- lations are in agreement with the results of previous works regarding the indirect influence of the Galaxy on the stability of planetary systems in wide binary stars. However, the effects on the planetary region show a dependence on the initial configuration of binary stars. Finally, we obtain an indirect test of the impulse approximation model for dynamical studies of binary star systems.

Context. Among the main effects that the Milky Way exerts in binary systems, the Galactic tide is the only one that is not probabilistic and can be deduced from a potential. Therefore, it is possible to perform an analysis of the global structure of the phase space of binary systems in the solar neighbourhood using the Galactic potential. Aims. The aim of this work is to obtain a simple model to study the collisionless dynamical evolution of generic wide binaries systems in the solar neighbourhood. Methods. Through an averaging process, we reduced the three-dimensional potential of the Galaxy to a secular one-degree of freedom model. The accuracy of this model was tested by comparing its predictions with numerical simulations of the exact equations of motion of a two-body problem disturbed by the Galaxy. Results. Using the one-degree of freedom model, we developed a detailed dynamical study, finding that the secular Galactic tide period changes as a function of the separation of the pair, which also gives a dynamical explanation for the arbitrary classification between "wide" and "tight" binaries. Moreover, the secular phase space for a generic gravitationally bound pair is similar to the dynamical structure of a Lidov-Kozai resonance, but surprisingly this structure is independent of the masses and semimajor axis of the binary system. Thus, the Galactic potential is able to excite the initially circular orbit of binary systems to high values of eccentricity, which has important implications for studies of binary star systems (with and without exoplanets), comets, and Oort cloud objects.

By considering all asteroid linear polarization data available in the literature, it is possible to obtain updated phase - polarization curves for several tens of objects. In a separate paper (Cellino et al., 2015a, MNRAS, 451, 3473) we have produced new calibrations of different relations between the geometric albedo and several polarimetric parameters, based on an analysis of a limited sample of asteroids for which the albedo is known with sufficient accuracy. In this paper, we present the main polarization parameters and corresponding albedos for a larger dataset of asteroids which we did not use for calibration purposes. We find a good agreement between the albedo values computed using different polarization parameters. Conversely, in the case of the so-called Barbarian asteroids the situation is rather unclear. Moreover, we present an updated analysis of the distributions of different polarimetric parameters, including the so-called inversion angle and the solar phase angle corresponding to the extreme value of negative polarization, and study their mutual relations. We find that the above parameters can be used to clearly distinguish some unusual classes of asteroids. Polarimetric parameters are known to be related to physical properties of asteroid surfaces which are difficult to infer by means of other observing techniques. By using a much larger dataset, in our analysis we confirm and extend some results obtained in the past by other authors, and we explore more systematically some features that had been mostly unexplored before, mainly concerning the morphology of the negative polarization branch.

We analyze two multi-chord stellar occultations by Pluto observed on July 18th, 2012 and May 4th, 2013, and monitored respectively from five and six sites. They provide a total of fifteen light-curves, twelve of them being used for a simultaneous fit that uses a unique temperature profile, assuming a clear (no-haze) and pure N_2 atmosphere, but allowing for a possible pressure variation between the two dates. We find a solution that fits satisfactorily (i.e. within the noise level) all the twelve light-curves, providing atmospheric constraints between ~1,190 km (pressure ~ 11 \mubar) and ~ 1,450 km (pressure ~0.1 \mubar) from Pluto's center. Our main results are: (1) the best-fitting temperature profile shows a stratosphere with strong positive gradient between 1,190 km (at 36 K, 11 \mubar) and r = 1,215 km (6.0 \mubar), where a temperature maximum of 110 K is reached; above it is a mesosphere with negative thermal gradient of -0.2 K/km up to ~ 1,390 km (0.25 \mubar), where, the mesosphere connects itself to a more isothermal upper branch around 81 K; (2) the pressure shows a small (6 %) but significant increase (6-\sigma level) between the two dates; (3) without troposphere, Pluto's radius is found to be R_P = 1,190 +/- 5km. Allowing for a troposphere, R_P is constrained to lie between 1,168 and 1,195 km; (4) the currently measured CO abundance is too small to explain the mesospheric negative thermal gradient. Cooling by HCN is possible, but only if this species is largely saturated; Alternative explanations like zonal winds or vertical compositional variations of the atmosphere are unable to explain the observed mesospheric trend.

We present a new extensive analysis of the old problem of finding a satisfactory calibration of the relation between the geometric albedo and some measurable polarization properties of the asteroids. To achieve our goals, we use all polarimetric data at our disposal. For the purposes of calibration, we use a limited sample of objects for which we can be confident to know the albedo with good accuracy, according to previous investigations of other authors. We find a new set of updated calibration coefficients for the classical slope - albedo relation, but we generalize our analysis and we consider also alternative possibilities, including the use of other polarimetric parameters, one being proposed here for the first time, and the possibility to exclude from best-fit analyzes the asteroids having low albedos. We also consider a possible parabolic fit of the whole set of data.

Until now, rings have been detected in the Solar System exclusively around the four giant planets. Here we report the discovery of the first minor-body ring system around the Centaur object (10199) Chariklo, a body with equivalent radius 124$\pm$9 km. A multi-chord stellar occultation revealed the presence of two dense rings around Chariklo, with widths of about 7 km and 3 km, optical depths 0.4 and 0.06, and orbital radii 391 and 405 km, respectively. The present orientation of the ring is consistent with an edge-on geometry in 2008, thus providing a simple explanation for the dimming of Chariklo's system between 1997 and 2008, and for the gradual disappearance of ice and other absorption features in its spectrum over the same period. This implies that the rings are partially composed of water ice. These rings may be the remnants of a debris disk, which were possibly confined by embedded kilometre-sized satellites.

The Mars crossing region constitutes a path to deliver asteroids from the Inner Main Belt to the Earth crossing space. While both the Inner Main Belt and the population of Earth crossing asteroids contains a significant fraction of asteroids belonging to the V taxonomic class, only two of such V-type asteroids has been detected in the Mars crossing region up to now. In this work, we searched for asteroids belonging to the V class among the population of Mars crossing asteroids, in order to support alternative paths to the delivery of this bodies into the Earth crossing region. We selected 18 candidate V-type asteroids in the Mars crossing region using observations contained in the Sloan Digital Sky Survey Moving Objects Catalog. Then, we observed 4 of these candidates to take their visible spectra using the Southern Astrophysical Research Telescope (SOAR). We also performed the numerical simulation of the orbital evolution of the observed asteroids. We confirmed that 3 of the observed asteroids belong to the V class, and one of these may follow a path that drives it to an Earth collision in some tens of million years.

At present, approximately 1500 asteroids are known to evolve inside or sticked to the exterior 1:2 resonance with Mars at a = 2.418 AU, being (142) Polana the largest member of this group. The effect of the forced secular modes superposed to the resonance gives rise to a complex dynamical evolution. Chaotic diffusion, collisions, close encounters with massive asteroids and mainly orbital migration due to the Yarkovsky effect generate continuous captures to and losses from the resonance, with a fraction of asteroids remaining captured over long time scales and generating a concentration in the semimajor axis distribution that exceeds by 20% the population of background asteroids. The Yarkovsky effect induces different dynamics according to the asteroid size, producing an excess of small asteroids inside the resonance. The evolution in the resonance generates a signature on the orbits, mainly in eccentricity, that depends on the time the asteroid remains captured inside the resonance and on the magnitude of the Yarkovsky effect. The greater the asteroids, the larger the time they remain captured in the resonance, allowing greater diffusion in eccentricity and inclination. The resonance generates a discontinuity and mixing in the space of proper elements producing misidentification of dynamical family members, mainly for Vesta and Nysa-Polana families. The half-life of resonant asteroids large enough for not being affected by the Yarkovsky effect is about 1 Gyr. From the point of view of taxonomic classes, the resonant population does not differ from the background population and the excess of small asteroids is confirmed.

From CCD observations of a fixed and large star field that contained the binary TNO Orcus, we have been able to derive high-precision relative astrometry and photometry of the Orcus system with respect to background stars. The RA residuals of an orbital fit to the astrometric data revealed a periodicity of 9.7+-0.3 days, which is what one would expect to be induced by the known Orcus companion. The residuals are also correlated with the theoretical positions of the satellite with regard to the primary. We therefore have revealed the presence of Orcus' satellite in our astrometric measurements. The photocenter motion is much larger than the motion of Orcus around the barycenter, and we show here that detecting some binaries through a carefully devised astrometric technique might be feasible with telescopes of moderate size. We also analyzed the system's mid-term photometry to determine whether the rotation could be tidally locked to the satellite's orbital period. We found that a photometric variability of 9.7+-0.3 days is clear in our data, and is nearly coincident with the orbital period of the satellite. We believe this variability might be induced by the satellite's rotation. There is also a slight hint for an additional small variability in the 10 hr range that was already reported in the literature. This short-term variability would indicate that the primary is not tidally locked and therefore the system would not have reached a double synchronous state. Implications for the basic physical properties of the primary and its satellite are discussed. From angular momentum considerations we suspect that the Orcus satellite might have formed from a rotational fission. This requires that the mass of the satellite would be around 0.09 times that of the primary, close to the value that one derives by using an albedo of 0.12 for the satellite and assuming equal densities for both objects.

The goal of the paper is to present new results on light echoes from young stellar objects. Broad band CCD images were obtained over three months at one-to-two week intervals for the field of NGC 6726, using the large field-of-view remotely-operated telescope on top of Cerro Burek. We detected scattered light echoes around two young, low-amplitude, irregular variable stars. Observations revealed not just one, but multiple light echoes from brightness pulses of the T Tauri star S CrA and the Herbig Ae/Be star R CrA. Analysis of S CrA's recurring echoes suggests that the star is located 138 +/- 16 pc from Earth, making these the closest echoes ever detected. The environment that scatters the stellar light from S CrA is compatible with an incomplete dust shell or an inclined torus some 10,000 AU in radius and containing $\sim$ $2 \times 10^{-3}$ $M_{\sun}$ of dust. The cause of such concentration at $\sim$ 10,000AU from the star is unknown. It could be the remnant of the envelope from which the star formed, but the distance of the cloud is remarkably similar to the nominal distance of the Oort cloud to the Sun, leading us to also speculate that the dust (or ice) seen around S CrA might have the same origin as the Solar System Oort cloud.

A wide observational campaign was carried out in 2004-2009 aimed to complete the ground-based investigation of Lutetia prior to the Rosetta fly-by in July 2010. We have obtained BVRI photometric and V-band polarimetric measurements over a wide range of phase angles, and visible and infrared spectra in the 0.4-2.4 micron range. We analyzed them together with previously published data to retrieve information on Lutetia's surface properties. Values of lightcurve amplitudes, absolute magnitude, opposition effect, phase coefficient and BVRI colors of Lutetia surface seen at near pole-on aspect have been determined. We defined more precisely parameters of polarization phase curve and showed their distinct deviation from any other moderate-albedo asteroid. An indication of possible variations both in polarization and spectral data across the asteroid surface was found. To explain features found by different techniques we propose that (i) Lutetia has a non-convex shape, probably due to the presence of a large crater, and heterogeneous surface properties probably related to surface morphology; (ii) at least part of the surface is covered by a fine-grained regolith with particle size less than 20 microns; (iii) the closest meteorite analogues of Lutetia's surface composition are particular types of carbonaceous chondrites or Lutetia has specific surface composition not representative among studied meteorites.

2003 EL61 is the largest member of a group of TNOs with similar orbits and 'unique' spectra (neutral slope in the visible and the deepest water ice absorption bands ever observed in the TNb). Studying the composition of the surface of 2003 EL61 provides useful constrains on the origin of this particular group of TNOs and on the outer Solar system's history. We present visible and near-infrared spectra of 2003 EL61 obtained with the 4.2m WHT and the 3.6m TNG at the Roque de los Muchachos Observatory (Canary Islands, Spain). Near infrared spectra were obtained at different rotational phases covering almost one complete rotational period. Spectra are fitted using Hapke scattering models and constraints on the surface composition are derived. No significant variations in the spectral slope and in the depth of the water ice absorption bands at different rotational phases are evident, suggesting that the surface of 2003 EL61 is homogeneous. The scattering models show that a 1:1 intimate mixture of crystalline and amorphous water ice is the most probable composition for the surface of this TNO, and constrain the presence of other minor constituents to a maximum of 8% The derived composition suggests that: a) cryovolcanism is unlikely to be the main resurfacing process responsible for the high presence of water ice on the surface of these bodies; b) the surface is older than 10^8 yr. Any catastrophic event, like the collision suggested to be the origin of this population, had to happen at least 10^8 yr ago; c) the surface of 2003 EL61 is depleted of carbon chains. According to the orbital parameters of this population, this makes it a possible source of carbon-depleted Jupiter family comets.

We present a comparative analysis of the spectral slope and color distributions of Jupiter Trojans, with particular attention to asteroid families. We use a sample of data from the Moving Object Catalogue of the Sloan Digital Sky Survey, together with spectra obtained from several surveys. A first sample of 349 observations, corresponding to 250 Trojan asteroids, were extracted from the Sloan Digital Sky Survey, and we also extracted from the literature a second sample of 91 spectra, corresponding to 71 Trojans. The spectral slopes were computed by means of a least-squares fit to a straight line of the fluxes obtained from the Sloan observations in the first sample, and of the rebinned spectra in the second sample. In both cases the reflectance fluxes/spectra were renormalized to 1 at 6230 $\textrm{\AA}$. We found that the distribution of spectral slopes among Trojan asteroids shows a bimodality. About 2/3 of the objects have reddish slopes compatible with D-type asteroids, while the remaining bodies show less reddish colors compatible with the P-type and C-type classifications. The members of asteroid families also show a bimodal distribution with a very slight predominance of D-type asteroids, but the background is clearly dominated by the D-types. The L4 and L5 swarms show different distributions of spectral slopes, and bimodality is only observed in L4. These differences can be attributed to the asteroid families since the backgraound asteroids show the same slope distribtuions in both swarms. The analysis of individual families indicates that the families in L5 are taxonomically homogeneous, but in L4 they show a mixture of taxonomic types. We discuss a few scenarios that might help to interpret these results.

The recent discovery of the first V-type asteroid in the middle belt, (21238) 1995WV7, located at ~2.54 AU, raises the question of whether it came from (4) Vesta or not. In this paper, we present spectroscopic observations indicating the existence of another V-type asteroid at ~2.53 AU, (40521) 1999RL95, and we investigate the possibility that these two asteroids evolved from the Vesta family to their present orbits by drifting in semi-major axis due to the Yarkovsky effect. The main problem with this scenario is that the asteroids need to cross the 3/1 mean motion resonance with Jupiter, which is highly unstable. Combining numerical simulations of the orbital evolution, that include the Yarkovsky effect, with Monte Carlo models, we compute the probability of an asteroid of given diameter D to evolve from the Vesta family and to cross over the 3/1 resonance, reaching a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 1995WV7 has a low probability of having evolved through this mechanism due to its large size (~5 km). However, the mechanism might explain the orbit of smaller bodies like (40521) 1999RL95 (~3 km), provided that we assume that the Vesta family formed > 3.5 Gy ago. We estimate that about 10% or more of the V-type bodies with D>1 km may come from the Vesta family by crossing over the 3/1 resonance. The remaining 90% must have a different origin.

Recent results suggest that there is a group of TNOs (2003 EL61 being the biggest member), with surfaces composed of almost pure water ice and with very similar orbital elements. We study the surface composition of another TNO that moves in a similar orbit, 2005 RR43, and compare it with the surface composition of the other members of this group. We report visible and near-infrared spectra, obtained with the 4.2m William Herschel Telescope and the 3.58m Telescopio Nazionale Galileo at the "Roque de los Muchachos" Observatory (La Palma, Spain). The spectrum of 2005 RR43 is neutral in color in the visible and dominated by very deep water ice absorption bands in the near infrared (D= 70.3 +/- 2.1 % and 82.8 +/- 4.9 % at 1.5 \mu and 2.0 \mu respectively). It is very similar to the spectrum of the group of TNOs already mentioned. All of them present much deeper water ice absorption bands (D>40 %) than any other TNO except Charon. Scattering models show that its surface is covered by water ice, a significant fraction in crytalline state with no trace (5 % upper limit) of complex organics. Possible scenarios to explain the existence of this population of TNOs are discussed: a giant collision, an originally carbon depleted composition, or a common process of continuous resurfacing. We conclude that TNO 2005 RR43 is member of a group, may be a population, of TNOs clustered in the space of orbital parameters that show abundant water ice and no signs of complex organics. The lack of complex organics in their surfaces suggests a significant smaller fraction of carbonaceous volatiles like CH4 in this population than in "normal" TNOs. A carbon depleted population of TNOs could be the origin of the population of carbon depleted Jupiter family comets already noticed by A'Hearn et al. (1995).