How are accretion discs affected by their outflows? To address this question for white dwarfs accreting from cool giants, we performed optical, radio, X-ray, and ultraviolet observations of the outflow-driving symbiotic star MWC 560 (=V694 Mon) during its 2016 optical high state. We tracked multi-wavelength changes that signalled an abrupt increase in outflow power at the initiation of a months-long outflow fast state, just as the optical flux peaked: (1) an abrupt doubling of Balmer absorption velocities; (2) the onset of a $20$ $\mu$Jy/month increase in radio flux; and (3) an order-of-magnitude increase in soft X-ray flux. Juxtaposing to prior X-ray observations and their coeval optical spectra, we infer that both high-velocity and low-velocity optical outflow components must be simultaneously present to yield a large soft X-ray flux, which may originate in shocks where these fast and slow absorbers collide. Our optical and ultraviolet spectra indicate that the broad absorption-line gas was fast, stable, and dense ($\gtrsim10^{6.5}$ cm$^{-3}$) throughout the 2016 outflow fast state, steadily feeding a lower-density ($\lesssim10^{5.5}$ cm$^{-3}$) region of radio-emitting gas. Persistent optical and ultraviolet flickering indicate that the accretion disc remained intact. The stability of these properties in 2016 contrasts to their instability during MWC 560's 1990 outburst, even though the disc reached a similar accretion rate. We propose that the self-regulatory effect of a steady fast outflow from the disc in 2016 prevented a catastrophic ejection of the inner disc. This behaviour in a symbiotic binary resembles disc/outflow relationships governing accretion state changes in X-ray binaries.

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 continue our investigation of the bulk properties of asteroid dynamical families identified using only asteroid proper elements (Milani et al. 2014) to provide plausible collisional interpretations. We focus on cratering families consisting of a substantial parent body and many small fragments. We propose a quantitative definition of cratering families based on the fraction in volume of the fragments with respect to the parent body; fragmentation families are above this empirical boundary. We assess the compositional homogeneity of the families and their shape in proper element space by computing the differences of the proper elements of the fragments with respect to the ones of the major body, looking for anomalous asymmetries produced either by post-formation dynamical evolution, or by multiple collisional/cratering events, or by a failure of the Hierarchical Clustering Method (HCM) for family identification. We identified a total of 25 dynamical families with more than 100 members ranging from moderate to heavy cratering. For three families (4, 15 and 283) we confirm the occurrence of two separate cratering events. For family 20, we propose a double collision origin, not previously identified. In four cases (31, 480, 163 and 179) we performed a dedicated search for dynamical resonant transport mechanisms that could have substantially changed the shape of the family. By using a new synthetic method for computation of secular frequencies, we found possible solutions for families 31, 480, and 163, but not for family 179, for which we propose a new interpretation, based on a secular resonance contaminating this family: the family of 179 should be split into two separate clusters, one containing (179) itself and the other, family (9506) Telramund, of fragmentation type, for which we have computed an age.

We present an updated set of near-Earth asteroids with a Yarkovsky-related semimajor axis drift detected from the orbital fit to the astrometry. We find 87 reliable detections after filtering for the signal-to-noise ratio of the Yarkovsky drift estimate and making sure the estimate is compatible with the physical properties of the analyzed object. Furthermore, we find a list of 24 marginally significant detections, for which future astrometry could result in a Yarkovsky detection. A further outcome of the filtering procedure is a list of detections that we consider spurious because unrealistic or not explicable with the Yarkovsky effect. Among the smallest asteroids of our sample, we determined four detections of solar radiation pressure, in addition to the Yarkovsky effect. As the data volume increases in the near future, our goal is to develop methods to generate very long lists of asteroids with reliably detected Yarkovsky effect, with limited amounts of case by case specific adjustments. Furthermore, we discuss the improvements this work could bring to impact monitoring. In particular, we exhibit two asteroids for which the adoption of a non-gravitational model is needed to make reliable impact predictions.

Short-arc orbit determination is crucial when an asteroid is first discovered. In these cases usually the observations are so few that the differential correction procedure may not converge. We have developed an initial orbit computation method, based on the systematic ranging, an orbit determination techniques which systematically explores a raster in the topocentric range and range-rate space region inside the admissible region. We obtain a fully rigorous computation of the probability for the asteroid that could impact the Earth within few days from the discovery, without any a priori assumption. We test our method on the two past impactors 2008 TC3 and 2014 AA, on some very well known cases, and on two particular objects observed by the ESA Gaia mission.

The JUICE mission will perform detailed measurements of the properties of the Galilean moons, with a nominal in-system science-mission duration of about 3.5 years. Using both the radio tracking data, and (Earth- and JUICE-based) optical astrometry, the dynamics of the Galilean moons will be measured to unprecedented accuracy. This will provide crucial input to the determination of the ephemerides and physical properties of the system, most notably the dissipation in Io and Jupiter. The data from Planetary Radio Interferometry and Doppler Experiment (PRIDE) will provide the lateral position of the spacecraft in the ICRF. In this article, we analyze the relative quantitative influence of the JUICE-PRIDE observables to the determination of the ephemerides of the Jovian system and the associated physical parameters. We perform a covariance analysis for a broad range of mission and system characteristics. We analyze the influence of VLBI data quality, observation planning, as well as the influence of JUICE orbit determination quality. This provides key input for the further development of the PRIDE observational planning and ground segment development. Our analysis indicates that the VLBI data are especially important for constraining the dynamics of Ganymede and Callisto perpendicular to their orbital planes. Also, the use of the VLBI data makes the uncertainty in the ephemerides less dependent on the error in the orbit determination of the JUICE spacecraft itself. Furthermore, we find that optical astrometry data of especially Io using the JANUS instrument will be crucial for stabilizing the solution of the normal equations. Knowledge of the dissipation in the Jupiter system cannot be improved using satellite dynamics obtained from JUICE data alone, the uncertainty in Io's dissipation obtained from our simulations is similar to the present level of uncertainty.

In this work we have estimated 10 collisional ages of 9 families for which for different reasons our previous attempts failed. In general, these are difficult cases that required dedicated effort, such as a new family classifications for asteroids in mean motion resonances, as well as a revision of the classification inside the $3/2$ resonance. Of the families locked in mean motion resonances, we succeeded in determining ages of the families of (1911) Schubart and of the "super-Hilda" family, assuming this is actually a severely eroded original family of (153) Hilda. In the Trojan region we found families with almost no Yarkovsky evolution, for which we could compute only physically implausible ages. Hence, we interpreted their modest dispersions of proper eccentricities and inclinations as implying that the Trojan asteroid families are fossil families, frozen at their proper elements determined by the original ejection velocity field. We have found a new family, among the Griquas locked in the 2/1 resonance with Jupiter: (11097) 1994 UD1. We have estimated the ages of 6 families affected by secular resonances: families of (5) Astraea, (25) Phocaea, (283) Emma, (363) Padua, (686) Gersuind, and (945) Barcelona. By using a numerical calibration method, we have shown that the secular resonances do not affect significanly the secular change of proper a. For the family of (145) Adeona we could estimate the age only after removal of a number of assumed interlopers. With the present paper we have concluded the series dedicated to the determination of asteroid ages with a uniform method. We computed the ages for a total of 57 families with $>100$ members. There remain families too small at present to provide reliable estimates, as well as some complex families (221, 135, 298) which may have more ages than we could currently estimate.

The modern optical telescopes produce a huge number of asteroid observations, that are grouped into very short arcs (VSAs), each containing a few observations of the same object in one single night. To decide whether two VSAs, collected in different nights, refer to the same observed object we can attempt to compute an orbit with the observations of both arcs: this is called the linkage problem. Since the number of orbit computations to be performed is very large, we need efficient methods of orbit determination. Using the first integrals of Kepler's motion we can write algebraic equations for the linkage problem, which can be put in polynomial form. The equations introduced in (Gronchi et al. 2015) can be reduced to a univariate polynomial of degree 9: the unknown is the topocentric distance $\rho$ of the observed body at the mean epoch of one of the VSAs. Using elimination theory we show an optimal property of this polynomial: it has the least degree among the univariate polynomials in the same variable that are consequence of the algebraic conservation laws and are obtained without squaring operations, that can be used to bring these algebraic equations in polynomial form. In this paper we also introduce a procedure to join three VSAs belonging to different nights: from the conservation of angular momentum at the three mean epochs of the VSAs, we obtain a univariate polynomial equation of degree 8 in the topocentric distance $\rho_2$ at the intermediate epoch. This algorithm has the same computational complexity as the classical method by Gauss, but uses more information, therefore we expect that it can produce more accurate results. For both methods, linking two and three VSAs, we also discuss how to select the solutions, making use of the full two-body dynamics, and show some numerical tests comparing the results with the ones obtained by Gauss' method.

Orbit determination is possible for a chaotic orbit of a dynamical system, given a finite set of observations, provided the initial conditions are at the central time. In a simple discrete model, the standard map, we tackle the problem of chaotic orbit determination when observations extend beyond the predictability horizon. If the orbit is hyperbolic, a shadowing orbit is computed by the least squares orbit determination. We test both the convergence of the orbit determination iterative procedure and the behaviour of the uncertainties as a function of the maximum number $n$ of map iterations observed. When the initial conditions belong to a chaotic orbit, the orbit determination is made impossible by numerical instability beyond a computability horizon, which can be approximately predicted by a simple formula. Moreover, the uncertainty of the results is sharply increased if a dynamical parameter is added to the initial conditions as parameter to be estimated. The uncertainty of the dynamical parameter decreases like $n^a$ with $a<0$ but not large (of the order of unity). If only the initial conditions are estimated, their uncertainty decreases exponentially with $n$. If they belong to a non-chaotic orbit the computational horizon is much larger, if it exists at all, and the decrease of the uncertainty is polynomial in all parameters, like $n^a$ with $a\simeq 1/2$. The Shadowing Lemma does not dictate what the asymptotic behaviour of the uncertainties should be. These phenomena have significant implications, which remain to be studied, in practical problems of orbit determination involving chaos, such as the chaotic rotation state of a celestial body and a chaotic orbit of a planet-crossing asteroid undergoing many close approaches.

A new family classification, based on a catalog of proper elements with $\sim 384,000$ numbered asteroids and on new methods is available. For the $45$ dynamical families with $>250$ members identified in this classification, we present an attempt to obtain statistically significant ages: we succeeded in computing ages for $37$ collisional families. We used a rigorous method, including a least squares fit of the two sides of a V-shape plot in the proper semimajor axis, inverse diameter plane to determine the corresponding slopes, an advanced error model for the uncertainties of asteroid diameters, an iterative outlier rejection scheme and quality control. The best available Yarkovsky measurement was used to estimate a calibration of the Yarkovsky effect for each family. The results are presented separately for the families originated in fragmentation or cratering events, for the young, compact families and for the truncated, one-sided families. For all the computed ages the corresponding uncertainties are provided. We found 2 cases where two separate dynamical families form together a single V-shape with compatible slopes, thus indicating a single collisional event. We have also found 3 examples of dynamical families containing multiple collisional families, plus a dubious case. We have found 2 cases of families containing a conspicuous subfamily, such that it is possible to measure the slope of a distinct V-shape, thus the age of the secondary collision. We also provide data on the central gaps appearing in some families. The ages computed in this paper are obtained with a single and uniform methodology, thus the ages of different families can be compared, providing a first example of collisional chronology of the asteroid main belt.

Juno is a NASA mission launched in 2011 with the goal of studying Jupiter. The probe will arrive to the planet in 2016 and will be placed for one year in a polar high-eccentric orbit to study the composition of the planet, the gravity and the magnetic field. The Italian Space Agency (ASI) provided the radio science instrument KaT (Ka-Band Translator) used for the gravity experiment, which has the goal of studying the Jupiter's deep structure by mapping the planet's gravity: such instrument takes advantage of synergies with a similar tool in development for BepiColombo, the ESA cornerstone mission to Mercury. The Celestial Mechanics Group of the University of Pisa, being part of the Juno Italian team, is developing an orbit determination and parameters estimation software for processing the real data independently from NASA software ODP. This paper has a twofold goal: first, to tell about the development of this software highlighting the models used, second, to perform a sensitivity analysis on the parameters of interest to the mission.

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.

The target asteroid of the OSIRIS-REx asteroid sample return mission, (101955) Bennu (formerly 1999 RQ$_{36}$), is a half-kilometer near-Earth asteroid with an extraordinarily well constrained orbit. An extensive data set of optical astrometry from 1999--2013 and high-quality radar delay measurements to Bennu in 1999, 2005, and 2011 reveal the action of the Yarkovsky effect, with a mean semimajor axis drift rate $da/dt = (-19.0 \pm 0.1)\times 10^{-4}$ au/Myr or $284\pm 1.5\;\rm{m/yr}$. The accuracy of this result depends critically on the fidelity of the observational and dynamical model. As an example, neglecting the relativistic perturbations of the Earth during close approaches affects the orbit with $3\sigma$ significance in $da/dt$. The orbital deviations from purely gravitational dynamics allow us to deduce the acceleration of the Yarkovsky effect, while the known physical characterization of Bennu allows us to independently model the force due to thermal emissions. The combination of these two analyses yields a bulk density of $\rho = 1260\pm70\,\rm{kg/m^3}$, which indicates a macroporosity in the range $40\pm10$% for the bulk densities of likely analog meteorites, suggesting a rubble-pile internal structure. The associated mass estimate is $(7.8\pm0.9)\times 10^{10}\, \rm{kg}$ and $GM = 5.2\pm0.6\,\rm{m^3/s^2}$. Bennu's Earth close approaches are deterministic over the interval 1654--2135, beyond which the predictions are statistical in nature. In particular, the 2135 close approach is likely within the lunar distance and leads to strong scattering and therefore numerous potential impacts in subsequent years, from 2175--2196. The highest individual impact probability is $9.5\times 10^{-5}$ in 2196, and the cumulative impact probability is $3.7\times 10^{-4}$, leading to a cumulative Palermo Scale of -1.70.

The number of asteroids with accurately determined orbits increases fast. The catalogs of asteroid physical observations have also increased, although the number of objects is still smaller than in the orbital catalogs. We developed a new approach to the asteroid family classification by combining the Hierarchical Clustering Method (HCM) with a method to add new members to existing families. This procedure makes use of the much larger amount of information contained in the proper elements catalogs, with respect to classifications using also physical observations for a smaller number of asteroids. Our work is based on the large catalog of the high accuracy synthetic proper elements (available from AstDyS). We first identify a number of core families; to these we attribute the next layer of smaller objects. Then, we remove all the family members from the catalog, and reapply the HCM to the rest. This gives both halo families which extend the core families and new independent families, consisting mainly of small asteroids. These two cases are discriminated by another step of attribution of new members and by merging intersecting families. By using information from absolute magnitudes, we take advantage of the larger size range in some families to analyze their shape in the proper semimajor axis vs. inverse diameter plane. This leads to a new method to estimate the family age (or ages). The results from the previous steps are then analyzed, using also auxiliary information on physical properties including WISE albedos and SDSS color indexes. This allows to solve some difficult cases of families overlapping in the proper elements space but generated by different collisional events. We analyze some examples of cratering families (Massalia, Vesta, Eunomia) which show internal structures, interpreted as multiple collisions. We also discuss why Ceres has no family.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

We assess the risk of an Earth impact for asteroid (99942) Apophis by means of a statistical analysis accounting for the uncertainty of both the orbital solution and the Yarkovsky effect. We select those observations with either rigorous uncertainty information provided by the observer or a high established accuracy. For the Yarkovsky effect we perform a Monte Carlo simulation that fully accounts for the uncertainty in the physical characterization, especially for the unknown spin orientation. By mapping the uncertainty information onto the 2029 b-plane and identifying the keyholes corresponding to subsequent impacts we assess the impact risk for future encounters. In particular, we find an impact probability greater than 10^-6 for an impact in 2068. We analyze the stability of the impact probability with respect to the assumptions on Apophis' physical characterization and consider the possible effect of the early 2013 radar apparition.

We seek evidence of the Yarkovsky effect among Near Earth Asteroids (NEAs) by measuring the Yarkovsky-related orbital drift from the orbital fit. To prevent the occurrence of unreliable detections we employ a high precision dynamical model, including the Newtonian attraction of 16 massive asteroids and the planetary relativistic terms, and a suitable astrometric data treatment. We find 21 NEAs whose orbital fits show a measurable orbital drift with a signal to noise ratio (SNR) greater than 3. The best determination is for asteroid (101955) 1999 RQ36, resulting in the recovery of one radar apparition and an orbit improvement by two orders of magnitude. In addition, we find 16 cases with a lower SNR that, despite being less reliable, are good candidates for becoming stronger detections in the future. In some cases it is possible to constrain physical quantities otherwise unknown by means of the detected orbital drift. Furthermore, the distribution of the detected orbital drifts shows an excess of retrograde rotators that can be connected to the delivery mechanism from the most important NEA feeding resonances and allows us to infer the distribution for NEAs obliquity. We discuss the implications of the Yarkovsky effect for impact predictions. In particular, for asteroid (29075) 1950 DA our results favor a retrograde rotation that would rule out an impact in 2880.

The discovery of new objects in modern wide-field asteroid and comet surveys can be enhanced by first identifying observations belonging to known solar system objects. The assignation of new observations to a known object is an attribution problem that occurs when a least squares orbit already exists for the object but a separate fit is not possible to just the set of new observations. In this work we explore the strongly asymmetric attribution problem in which the existing least squares orbit is very well constrained and the new data are sparse. We describe an attribution algorithm that introduces new quality control metrics in the presence of strong biases in the astrometric residuals. The main biases arise from the stellar catalogs used in the reduction of asteroid observations and we show that a simple debiasing with measured regional catalog biases significantly improves the results. We tested the attribution algorithm using data from the PS1 survey that used the 2MASS star catalog for the astrometric reduction. We found small but statistically significant biases in the data of up to 0.1 arcsec that are relevant only when the observations reach the level of accuracy made possible by instruments like PS1. The false attribution rate was measured to be < 1/1,000 with a simple additional condition that can reduce it to zero while the attribution efficiency is consistent with 100%.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

We present the results of a large scale simulation, reproducing the behavior of a data center for the build-up and maintenance of a complete catalog of space debris in the upper part of the low Earth orbits region (LEO). The purpose is to determine the performances of a network of advanced optical sensors, through the use of the newest orbit determination algorithms developed by the Department of Mathematics of Pisa (DM). Such a network has been proposed to ESA in the Space Situational Awareness (SSA) framework by Carlo Gavazzi Space SpA (CGS), Istituto Nazionale di Astrofisica (INAF), DM, and Istituto di Scienza e Tecnologie dell'Informazione (ISTI-CNR). The conclusion is that it is possible to use a network of optical sensors to build up a catalog containing more than 98% of the objects with perigee height between 1100 and 2000 km, which would be observable by a reference radar system selected as comparison. It is also possible to maintain such a catalog within the accuracy requirements motivated by collision avoidance, and to detect catastrophic fragmentation events. However, such results depend upon specific assumptions on the sensor and on the software technologies.

We report the discovery of a large amplitude outburst from the young star HBC 722 (LkHA 188 G4) located in the region of NGC 7000/IC 5070. On the basis of photometric and spectroscopic observations, we argue that this outburst is of the FU Orionis type. We gathered photometric and spectroscopic observations of the object both in the pre-outburst state and during a phase of increase in its brightness. The photometric BVRI data (Johnson-Cousins system) that we present were collected from April 2009 to September 2010. To facilitate transformation from instrumental measurements to the standard system, fifteen comparison stars in the field of HBC 722 were calibrated in the BVRI bands. Optical spectra of HBC 722 were obtained with the 1.3-m telescope of Skinakas Observatory (Crete, Greece) and the 0.6-m telescope of Schiaparelli Observatory in Varese (Italy). The pre-outburst photometric and spectroscopic observations of HBC 722 show both low amplitude photometric variations and an emission-line spectrum typical of T Tau stars. The observed outburst started before May 2010 and reached its maximum brightness in September 2010, with a recorded Delta V~4.7 mag. amplitude. Simultaneously with the increase in brightness the color indices changed significantly and the star became appreciably bluer. The light curve of HBC 722 during the period of rise in brightness is similar to the light curves of the classical FUors - FU Ori and V1057 Cyg. The spectral observations during the time of increase in brightness showed significant changes in both the profiles and intensity of the spectral lines. Only H alpha remained in emission, while the H beta, Na I 5890/5896, Mg I triplet 5174, and Ba II 5854/6497 lines were in strong absorption.

The nova outburst experienced in 2010 by the symbiotic binary Mira V407 Cyg has been extensively studied at optical and infrared wavelengths with both photometric and spectroscopic observations. This outburst, reminiscent of similar events displayed by RS Oph, can be described as a very fast He/N nova erupting while being deeply embedded in the dense wind of its cool giant companion. The hard radiation from the initial thermonuclear flash ionizes and excites the wind of the Mira over great distances (recombination is observed on a time scale of 4 days). The nova ejecta is found to progressively decelerate with time as it expands into the Mira wind. This is deduced from line widths which change from a FWHM of 2760 km/s on day +2.3 to 200 km/s on day +196. The wind of the Mira is massive and extended enough for an outer neutral and unperturbed region to survive at all outburst phases.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The radioscience experiment is one of the on board experiment of the Mercury ESA mission BepiColombo that will be launched in 2014. The goals of the experiment are to determine the gravity field of Mercury and its rotation state, to determine the orbit of Mercury, to constrain the possible theories of gravitation (for example by determining the post-Newtonian (PN) parameters), to provide the spacecraft position for geodesy experiments and to contribute to planetary ephemerides improvement. This is possible thanks to a new technology which allows to reach great accuracies in the observables range and range rate; it is well known that a similar level of accuracy requires studying a suitable model taking into account numerous relativistic effects. In this paper we deal with the modelling of the space-time coordinate transformations needed for the light-time computations and the numerical methods adopted to avoid rounding-off errors in such computations.

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.

Given a set of astrometric observations of the same object, the problem of orbit determination is to compute the orbit and to assess its uncertainty and reliability. For the next generation surveys, with much larger number density of observed objects, new algorithms or substantial revisions of the classical ones are needed. The problem has three main steps, preliminary orbit, least squares orbit, and quality control. The classical theory of preliminary orbits was incomplete: the consequences of the topocentric correction had not been fully studied. We show that it is possible to account for this correction, possibly with an increase in the number of preliminary solutions, without impairing the overall orbit determination performance. We have developed modified least squares orbit determination algorithms that can be used to improve the reliability of the procedure. We have tested the complete procedure on two simulations with number densities comparable to that expected from the next generation surveys such as Pan-STARRS and LSST. To control the problem of false identifications we have introduced a quality control on the fit residuals based on an array of metrics and a procedure to remove duplications and contradictions in the output. The results confirm that large sets of discoveries can be obtained with good quality orbits and very high success rate losing only 0.6 to 1.3% of objects and a false identification rate in the range 0.02 to 0.06%.

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

The Panoramic Survey Telescope And Rapid Response System (Pan-STARRS) under development at the University of Hawaii's Institute for Astronomy is creating the first fully automated end-to-end Moving Object Processing System (MOPS) in the world. It will be capable of identifying detections of moving objects in our solar system and linking those detections within and between nights, attributing those detections to known objects, calculating initial and differentially-corrected orbits for linked detections, precovering detections when they exist, and orbit identification. Here we describe new kd-tree and variable-tree algorithms that allow fast, efficient, scalable linking of intra and inter-night detections. Using a pseudo-realistic simulation of the Pan-STARRS survey strategy incorporating weather, astrometric accuracy and false detections we have achieved nearly 100% efficiency and accuracy for intra-night linking and nearly 100% efficiency for inter-night linking within a lunation. At realistic sky-plane densities for both real and false detections the intra-night linking of detections into `tracks' currently has an accuracy of 0.3%. Successful tests of the MOPS on real source detections from the Spacewatch asteroid survey indicate that the MOPS is capable of identifying asteroids in real data.