We propose a densification algorithm to improve the Line Of Variations (LOV) method for impact monitoring, which can fail when the information is too little, as it may happen in difficult cases. The LOV method uses a 1-dimensional sampling to explore the uncertainty region of an asteroid. The close approaches of the sample orbits are grouped by time and LOV index, to form the so-called returns, and each return is analysed to search for local minima of the distance from the Earth along the LOV. The strong non-linearity of the problem causes the occurrence of returns with so few points that a successful analysis can be prevented. Our densification algorithm tries to convert returns with length at most 3 in returns with 5 points, properly adding new points to the original return. Due to the complex evolution of the LOV, this operation is not necessarily achieved all at once: in this case the information about the LOV geometry derived from the first attempt is exploited for a further attempt. Finally, we present some examples showing that the application of our method can have remarkable consequences on impact monitoring results, in particular about the completeness of the virtual impactors search.

On May 30th, 2017 at about 21h 09m 17s UTC a green bright fireball crossed the sky of north-eastern Italy. The fireball path was observed from some all-sky cameras starting from a mean altitude of $81.1 \pm 0.2$ km (Lat. $44.369^{\circ} \pm 0.002^{\circ}$ N; Long. $11.859^{\circ} \pm 0.002^{\circ}$ E) and extinct at $23.3 \pm 0.2$ km (Lat. $45.246^{\circ} \pm 0.002^{\circ}$ N; Long. $12.046^{\circ} \pm 0.002^{\circ}$ E), between the Italian cities of Venice and Padua. In this paper, on the basis of simple physical models, we will compute the atmospheric trajectory, analize the meteoroid atmospheric dynamics, the dark flight phase (with the strewn field) and compute the best heliocentric orbit of the progenitor body. Search for meteorites on the ground has not produced any results so far.

We study the post-encounter evolution of fictitious small bodies belonging to the so-called Line of Variations (LoV) in the framework of the analytic theory of close encounters. We show the consequences of the encounter on the local minimum of the distance between the orbit of the planet and that of the small body, and get a global picture of the way in which the planetocentric velocity vector is affected by the encounter. The analytical results are compared with those of numerical integrations of the restricted 3-body problem.

Near-Earth asteroid (410777) 2009 FD is a potentially hazardous asteroid with potential impacts on Earth at the end of the 22nd century. The astrometry collected during the 2019 apparition provides information on the trajectory of (410777) by constraining the Yarkovsky effect, which is the main source of uncertainty for future predictions, and informing the impact hazard assessment. We included the Yarkovsky effect in the force model and estimated its magnitude from the fit to the (410777) optical and radar astrometric data. We performed the (410777) hazard assessment over 200 years by using two independent approaches: the NEODyS group adopted a generalisation of the Line Of Variations method in a 7-dimensional space, whereas the JPL team resorted to the Multilayer Clustered Sampling technique. We obtain a 4-sigma detection of the Yarkovsky effect acting on (410777), which corresponds to a semimajor axis drift of (3.6 +/- 0.9) x 10^{-3} au/Myr. As for the hazard results of both teams, the main impact possibility in 2185 is ruled out and the only remaining one is in 2190, but with a probability smaller than 10^{-8}.

The completeness limit is a key quantity to measure the reliability of an impact monitoring system. It is the impact probability threshold above which every virtual impactor has to be detected. A goal of this paper is to increase the completeness without increasing the computational load. We propose a new method to sample the Line Of Variations (LOV) with respect to the previously one used in NEODyS. The step-size of the sampling is not uniform in the LOV parameter, since the probability of each LOV segment between consecutive points is kept constant. Moreover, the sampling interval has been extended to the larger interval [-5,5] in the LOV parameter and a new decomposition scheme in sub-returns is provided to deal with the problem of duplicated points in the same return. The impact monitoring system CLOMON-2 has been upgraded with all these features, resulting in a decrease of the generic completeness limit IP* by a factor $\simeq 4$ and in an increase of the computational load by a factor $\simeq 2$. Moreover, we statistically investigate the completeness actually reached by the system with two different methods: a direct comparison with the results of the independent system Sentry at JPL and an empirical power-law to model the number of virtual impactors as a function of the impact probability. We found empirically that the number of detected virtual impactors with $IP>IP^*$ appears to grow according to a power-law, proportional to $IP^{-2/3}$. An analytical model explaining this power-law is currently an open problem, but we think it is related to the way the number of virtual impactors within a time $t_{rel}$ from the first observed close approach accumulates, for which we prove the existence of a power-law proportional to $t_{rel}^3$. The power-law allows us to detect a loss of efficiency in the virtual impactors search for impact probabilities near the generic completeness limit.

We present the main findings of a dynamical mapping performed in the Low Earth Orbit region. The results were obtained by propagating an extended grid of initial conditions, considering two different epochs and area-to-mass ratios, by means of a singly averaged numerical propagator. It turns out that dynamical resonances associated with high-degree geopotential harmonics, lunisolar perturbations and solar radiation pressure can open natural deorbiting highways. For area-to-mass ratios typical of the orbiting intact objects, these corridors can be exploited only in combination with the action exerted by the atmospheric drag. For satellites equipped with an area augmentation device, we show the boundary of application of the drag, and where the solar radiation pressure can be exploited.

The aim of this work is to highlight the crucial role that orbital resonances associated with solar radiation pressure can have in Low Earth Orbit. We review the corresponding literature, and provide an analytical tool to estimate the maximum eccentricity which can be achieved for well-defined initial conditions. We then compare the results obtained with the simplified model with the results obtained with a more comprehensive dynamical model. The analysis has important implications both from a theoretical point of view, because it shows that the role of some resonances was underestimated in the past, but also from a practical point of view in the perspective of passive deorbiting solutions for satellites at the end-of-life.

We use a secular model to describe the non-resonant dynamics of trans-Neptunian objects in the presence of an external ten-earth-mass perturber. The secular dynamics is analogous to an "eccentric Kozai mechanism" but with both an inner component (the four giant planets) and an outer one (the eccentric distant perturber). By the means of Poincar\'e sections, the cases of a non-inclined or inclined outer planet are successively studied, making the connection with previous works. In the inclined case, the problem is reduced to two degrees of freedom by assuming a non-precessing argument of perihelion for the perturbing body. The size of the perturbation is typically ruled by the semi-major axis of the small body: we show that the classic integrable picture is still valid below about 70 AU, but it is progressively destroyed when we get closer to the external perturber. In particular, for a>150 AU, large-amplitude orbital flips become possible, and for a>200 AU, the Kozai libration islands are totally submerged by the chaotic sea. Numerous resonance relations are highlighted. The most large and persistent ones are associated to apsidal alignments or anti-alignments with the orbit of the distant perturber.

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or 2012VP113, which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at semi-major axes in [100;300] AU, perihelion distances in [50;70] AU and inclinations [30;50]{\deg}. In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.

Two semi-analytical one-degree-of-freedom secular models are presented for the motion of small bodies beyond Neptune. A special attention is given to trajectories entirely exterior to the planetary orbits. The first one is the well-known non-resonant model of Kozai (1962) adapted to the transneptunian region. Contrary to previous papers, the dynamics is fully characterized with respect to the fixed parameters. A maximum perihelion excursion possible of 16.4 AU is determined. The second model handles the occurrence of a mean-motion resonance with one of the planets. In that case, the one-degree-of-freedom integrable approximation is obtained by postulating the adiabatic invariance, and is much more general and accurate than previous secular models found in the literature. It brings out in a plain way the possibility of perihelion oscillations with a very high amplitude. Such a model could thus be used in future studies to deeper explore that kind of motion. For complex resonant orbits (especially of type 1:k), a segmented secular description is necessary since the trajectories are only "integrable by parts". The two models are applied to the Solar System but the notations are kept general so that it could be used for any quasi-circular and coplanar planetary system.

The BEaTriX (Beam Expander Testing X-ray facility) project is an X-ray apparatus under construction at INAF/OAB to generate a broad (200 x 60 mm2), uniform and low-divergent X-ray beam within a small lab (6 x 15 m2). BEaTriX will consist of an X-ray source in the focus a grazing incidence paraboloidal mirror to obtain a parallel beam, followed by a crystal monochromation system and by an asymmetrically-cut diffracting crystal to perform the beam expansion to the desired size. Once completed, BEaTriX will be used to directly perform the quality control of focusing modules of large X-ray optics such as those for the ATHENA X-ray observatory, based on either Silicon Pore Optics (baseline) or Slumped Glass Optics (alternative), and will thereby enable a direct quality control of angular resolution and effective area on a number of mirror modules in a short time, in full X-ray illumination and without being affected by the finite distance of the X-ray source. However, since the individual mirror modules for ATHENA will have an optical quality of 3-4 arcsec HEW or better, BEaTriX is required to produce a broad beam with divergence below 1-2 arcsec, and sufficient flux to quickly characterize the PSF of the module without being significantly affected by statistical uncertainties. Therefore, the optical components of BEaTriX have to be selected and/or manufactured with excellent optical properties in order to guarantee the final performance of the system. In this paper we report the final design of the facility and a detailed performance simulation.

The LAMP (Lightweight Asymmetry and Magnetism Probe) X-ray telescope is a mission concept to measure the polarization of X-ray astronomical sources at 250 eV via imaging mirrors that reflect at incidence angles near the polarization angle, i.e., 45 deg. Hence, it will require the adoption of multilayer coatings with a few nanometers d-spacing in order to enhance the reflectivity. The nickel electroforming technology has already been successfully used to fabricate the high angular resolution imaging mirrors of the X-ray telescopes SAX, XMM-Newton, and Swift/XRT. We are investigating this consolidated technology as a possible technique to manufacture focusing mirrors for LAMP. Although the very good reflectivity performances of this kind of mirrors were already demonstrated in grazing incidence, the reflectivity and the scattering properties have not been tested directly at the unusually large angle of 45 deg. Other possible substrates are represented by thin glass foils or silicon wafers. In this paper we present the results of the X-ray reflectivity campaign performed at the BEAR beamline of Elettra - Sincrotrone Trieste on multilayer coatings of various composition (Cr/C, Co/C), deposited with different sputtering parameters on nickel, silicon, and glass substrates, using polarized X-rays in the spectral range 240 - 290 eV.

While one of the IAU's missions is to "serve as the internationally recognized authority for assigning designations to celestial bodies and surface features on them", the participation of the public in the naming of celestial objects has been a little-known, but decade-long tradition of the IAU. While reiterating its opposition to having the public pay to give a name to an exoplanet, the IAU Executive Committee nonetheless recognized the right of organizations to invite public, international exoplanet naming or voting campaigns. To this end, clear selection rules were to be defined by the IAU, inviting mutual collaboration, the goal being to sanction the campaign and officially approve the resulting names, for the sake of boosting the public's interest in astronomy and at the same time reaffirm the authority of the IAU. In no way were these names supposed to supersede the designations in use by professional astronomers. Since the field of exoplanet research was still "virgin" for naming but is now becoming scientifically mature (with hundreds of well-characterized, confirmed objects) and so appealing to the public's interest and imagination, IAU decided to organize a naming contest as the leader of a new project with a potentially high worldwide impact, the first since the International Year of Astronomy in 2009, called "NameExoWorlds".

Recent studies have shown that the medium-Earth orbit (MEO) region of the Global Navigation Satellite Systems is permeated by a devious network of lunisolar secular resonances, which can interact to produce chaotic and diffusive motions. The precarious state of the four navigation constellations, perched on the threshold of instability, makes it understandable why all past efforts to define stable graveyard orbits, especially in the case of Galileo, were bound to fail; the region is far too complex to allow of an adoption of the simple geosynchronous disposal strategy. We retrace one such recent attempt, funded by ESA's General Studies Programme in the frame of the GreenOPS initiative, that uses a systematic parametric approach and the straightforward maximum-eccentricity method to identify long-term stable regions, suitable for graveyards, as well as large-scale excursions in eccentricity, which can be used for post-mission deorbiting of constellation satellites. We then apply our new results on the stunningly rich dynamical structure of the MEO region toward the analysis of these disposal strategies for Galileo, and discuss the practical implications of resonances and chaos in this regime. We outline how the identification of the hyperbolic and elliptic fixed points of the resonances near Galileo can lead to explicit criteria for defining optimal disposal strategies.

The realization of X-ray telescopes with imaging capabilities in the hard (> 10 keV) X-ray band requires the adoption of optics with shallow (< 0.25 deg) grazing angles to enhance the reflectivity of reflective coatings. On the other hand, to obtain large collecting area, large mirror diameters (< 350 mm) are necessary. This implies that mirrors with focal lengths >10 m shall be produced and tested. Full-illumination tests of such mirrors are usually performed with on- ground X-ray facilities, aimed at measuring their effective area and the angular resolution; however, they in general suffer from effects of the finite distance of the X-ray source, e.g. a loss of effective area for double reflection. These effects increase with the focal length of the mirror under test; hence a "partial" full-illumination measurement might not be fully representative of the in-flight performances. Indeed, a pencil beam test can be adopted to overcome this shortcoming, because a sector at a time is exposed to the X-ray flux, and the compensation of the beam divergence is achieved by tilting the optic. In this work we present the result of a hard X-ray test campaign performed at the BL20B2 beamline of the SPring-8 synchrotron radiation facility, aimed at characterizing the Point Spread Function (PSF) of a multilayer-coated Wolter-I mirror shell manufactured by Nickel electroforming. The mirror shell is a demonstrator for the NHXM hard X-ray imaging telescope (0.3 - 80 keV), with a predicted HEW (Half Energy Width) close to 20 arcsec. We show some reconstructed PSFs at monochromatic X-ray energies of 15 to 63 keV, and compare them with the PSFs computed from post-campaign metrology data, self-consistently treating profile and roughness data by means of a method based on the Fresnel diffraction theory. The modeling matches the measured PSFs accurately.

It has long been suspected that the Global Navigation Satellite Systems exist in a background of complex resonances and chaotic motion; yet, the precise dynamical character of these phenomena remains elusive. Recent studies have shown that the occurrence and nature of the resonances driving these dynamics depend chiefly on the frequencies of nodal and apsidal precession and the rate of regression of the Moon's nodes. Woven throughout the inclination and eccentricity phase space is an exceedingly complicated web-like structure of lunisolar secular resonances, which become particularly dense near the inclinations of the navigation satellite orbits. A clear picture of the physical significance of these resonances is of considerable practical interest for the design of disposal strategies for the four constellations. Here we present analytical and semi-analytical models that accurately reflect the true nature of the resonant interactions, and trace the topological organization of the manifolds on which the chaotic motions take place. We present an atlas of FLI stability maps, showing the extent of the chaotic regions of the phase space, computed through a hierarchy of more realistic, and more complicated, models, and compare the chaotic zones in these charts with the analytical estimation of the width of the chaotic layers from the heuristic Chirikov resonance-overlap criterion. As the semi-major axis of the satellite is receding, we observe a transition from stable Nekhoroshev-like structures at three Earth radii, where regular orbits dominate, to a Chirikov regime where resonances overlap at five Earth radii. From a numerical estimation of the Lyapunov times, we find that many of the inclined, nearly circular orbits of the navigation satellites are strongly chaotic and that their dynamics are unpredictable on decadal timescales.

Numerical simulations carried out over the past decade suggest that the orbits of the Global Navigation Satellite Systems are unstable, resulting in an apparent chaotic growth of the eccentricity. Here we show that the irregular and haphazard character of these orbits reflects a similar irregularity in the orbits of many celestial bodies in our Solar System. We find that secular resonances, involving linear combinations of the frequencies of nodal and apsidal precession and the rate of regression of lunar nodes, occur in profusion so that the phase space is threaded by a devious stochastic web. As in all cases in the Solar System, chaos ensues where resonances overlap. These results may be significant for the analysis of disposal strategies for the four constellations in this precarious region of space.

Asteroid 2009 FD could impact Earth between 2185 and 2196. The long term propagation to the possible impacts and the intervening planetary encounters make 2009 FD one of the most challenging asteroids in terms of hazard assessment. To compute accurate impact probabilities we model the Yarkovsky effect by using the available physical characterization of 2009 FD and general properties of the Near Earth Asteroid population. We perform the hazard assessment with two independent methods: the first method is a generalization of the standard impact monitoring algorithms in use by NEODyS and Sentry, while the second one is based on a Monte Carlo approach. Both methods generate orbital samples in a 7 dimensional space that includes orbital elements and the parameter characterizing the Yarkovsky effect. The highest impact probability is $2.7 \times 10^{-3}$ for an impact during the 2185 Earth encounter. Impacts after 2185 corresponding to resonant returns are possible, the most relevant being in 2190 with a probability of $3 \times 10^{-4}$. Both numerical methods can be used in the future to handle similar cases. The structure of resonant returns and the list of the possible keyholes on the Target Plane of the scattering encounter in 2185 can be predicted by an analytic theory.

We consider a network of telescopes capable of scanning all the observable sky each night and targeting Near-Earth objects (NEOs) in the size range of the Tunguska-like asteroids, from 160 m down to 10 m. We measure the performance of this telescope network in terms of the time needed to discover at least 50% of the impactors in the considered population with a warning time large enough to undertake proper mitigation actions. The warning times are described by a trimodal distribution and the telescope network has a 50% probability of discovering an impactor of the Tunguska class with at least one week of advance already in the first 10 yr of operations of the survey. These results suggest that the studied survey would be a significant addition to the current NEO discovery efforts.

This paper tackles important aspects of comets dynamics from a statistical point of view. Existing methodology uses numerical integration for computing planetary perturbations for simulating such dynamics. This operation is highly computational. It is reasonable to wonder whenever statistical simulation of the perturbations can be much more easy to handle. The first step for answering such a question is to provide a statistical study of these perturbations in order to catch their main features. The statistical tools used are order statistics and heavy tail distributions. The study carried out indicated a general pattern exhibited by the perturbations around the orbits of the important planet. These characteristics were validated through statistical testing and a theoretical study based on Opik theory.

The potentially hazardous asteroid (101955) 1999 RQ36 has the possibility of collision with the Earth in the latter half of the 22nd century, well beyond the traditional 100-year time horizon for routine impact monitoring. The probabilities accumulate to a total impact probability of approximately 10E-3, with a pair of closely related routes to impact in 2182 comprising more than half of the total. The analysis of impact possibilities so far in the future is strongly dependent on the action of the Yarkovsky effect, which raises new challenges in the careful assessment of longer term impact hazards. Even for asteroids with very precisely determined orbits, a future close approach to Earth can scatter the possible trajectories to the point that the problem becomes like that of a newly discovered asteroid with a weakly determined orbit. If the scattering takes place late enough so that the target plane uncertainty is dominated by Yarkovsky accelerations then the thermal properties of the asteroid,which are typically unknown, play a major role in the impact assessment. In contrast, if the strong planetary interaction takes place sooner, while the Yarkovsky dispersion is still relatively small compared to that derived from the measurements, then precise modeling of the nongravitational acceleration may be unnecessary.

Simbol-X will push grazing incidence imaging up to 80 keV, providing a strong improvement both in sensitivity and angular resolution compared to all instruments that have operated so far above 10 keV. The superb hard X-ray imaging capability will be guaranteed by a mirror module of 100 electroformed Nickel shells with a multilayer reflecting coating. Here we will describe the technogical development and solutions adopted for the fabrication of the mirror module, that must guarantee an Half Energy Width (HEW) better than 20 arcsec from 0.5 up to 30 keV and a goal of 40 arcsec at 60 keV. During the phase A, terminated at the end of 2008, we have developed three engineering models with two, two and three shells, respectively. The most critical aspects in the development of the Simbol-X mirrors are i) the production of the 100 mandrels with very good surface quality within the timeline of the mission; ii) the replication of shells that must be very thin (a factor of 2 thinner than those of XMM-Newton) and still have very good image quality up to 80 keV; iii) the development of an integration process that allows us to integrate these very thin mirrors maintaining their intrinsic good image quality. The Phase A study has shown that we can fabricate the mandrels with the needed quality and that we have developed a valid integration process. The shells that we have produced so far have a quite good image quality, e.g. HEW <~30 arcsec at 30 keV, and effective area. However, we still need to make some improvements to reach the requirements. We will briefly present these results and discuss the possible improvements that we will investigate during phase B.

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1 % for a < 20 000 AU to about 100 % for a >= 100 000 AU.

In this paper we present the observational campaign carried out at ESO NTT and VLT in April and May 2006 to investigate the nature and the structure of the Near Earth Object (144898) 2004 VD17. In spite of a great quantity of dynamical information, according to which it will have a close approach with the Earth in the next century, the physical properties of this asteroid are largely unknown. We performed visible and near--infrared photometry and spectroscopy, as well as polarimetric observations. Polarimetric and spectroscopic data allowed us to classify 2004 VD17 as an E-type asteroid. A good agreement was also found with the spectrum of the aubrite meteorite Mayo Belwa. On the basis of the polarimetric albedo (p_v=0.45) and of photometric data, we estimated a diameter of about 320 m and a rotational period of about 2 hours. The analysis of the results obtained by our complete survey have shown that (144898) 2004 VD17 is a peculiar NEO, since it is close to the breakup limits for fast rotator asteroids, as defined by Pravec and Harris (2000). These results suggest that a more robust structure must be expected, as a fractured monolith or a rubble pile in a "strength regime" (Holsapple 2002).