The first known asteroid with the orbit inside that of Venus is 2020~AV$_{2}$. This may be the largest member of a new population of small bodies with the aphelion smaller than 0.718~au, called Vatiras. The surface of 2020~AV$_{2}$ is being constantly modified by the high temperature, by the strong solar wind irradiation that characterizes the innermost region of the Solar system, and by high-energy micrometeorite impacts. The study of its physical properties represents an extreme test-case for the science of near-Earth asteroids. Here, we report spectroscopic observations of 2020~AV$_{2}$ in the 0.5-1.5~$\mu m$ wavelength interval. These were performed with the Nordic Optical Telescope and the William Herschel Telescope. Based on the obtained spectra, we classify 2020~AV$_{2}$ as a Sa-type asteroid. We estimate the diameter of this Vatira to be $1.50_{-0.65}^{+1.10}$ km by considering the average albedo of A-type and S-complex asteroids ($p_V=0.23_{-0.08}^{+0.11}$), and the absolute magnitude (H=$16.40\pm0.78$ mag). The wide spectral band around 1~$\mu m$ shows the signature of an olivine rich composition. The estimated band centre $BIC = 1.08 \pm 0.02~\mu m$ corresponds to a ferroan olivine mineralogy similar to that of brachinite meteorites.

In this work, we present the results of an observational study of 2I/Borisov carried out with the 10.4-m Gran Telescopio Canarias (GTC) and the 3.6-m Telescopio Nazionale Galileo (TNG), both telescopes located at the Roque de Los Muchachos Observatory, in the island of La Palma (Spain). The study includes images in the visible and near-infrared, as well as visible spectra in the 3600 - 9200 A wavelength range. N-body simulations were also performed to explore its orbital evolution and Galactic kinematic context. The comet's dust continuum and near-infrared colours are compatible with those observed for Solar system comets. From its visible spectrum on the nights of 2019, September 24 and 26 we measured CN gas production rates Q(CN) = (2.3 +- 0.4) x 10^{24} mol/s and Q(CN) = (9.5 +- 0.2) x 10^{24} mol/s, respectively, in agreement with measurements reported by other authors on similar nights. We also obtained an upper limit for the C2 production rate of Q(C2) < (4.5 +- 0.1) x 10^{24} mol/s. Dust modelling results indicate a moderate dust production rate of about 50 kg/s at heliocentric distance r_h=2.6 au, with a differential power-law dust size distribution of index -3.4, within the range reported for many comet comae. Our simulations show that the Galactic velocity of 2I/Borisov matches well that of known stars in the solar neighbourhood and also those of more distant regions of the Galactic disc.

Any near-Earth object (NEO) following an Earth-like orbit may eventually be captured by Earth's gravity during low-velocity encounters. This theoretical possibility was first attested during the fly-by of 1991 VG in 1991-1992 with the confirmation of a brief capture episode ---for about a month in 1992 February. Further evidence was obtained when 2006 RH120 was temporarily captured into a geocentric orbit from 2006 July to 2007 July. Here, we perform a numerical assessment of the orbital evolution of 2020 CD3, a small NEO found recently that could be the third instance of a meteoroid temporarily captured by Earth's gravity. We confirm that 2020 CD3 is currently following a geocentric trajectory although it will escape into a heliocentric path by early May 2020. Our calculations indicate that it was captured by the Earth in 2016 (median). This episode is longer than that of 2006 RH120. Prior to its capture as a minimoon, 2020 CD3 was probably a NEO of the Aten type, but an Apollo type cannot be excluded; in both cases, the orbit was very Earth-like, with low eccentricity and low inclination, typical of an Arjuna-type meteoroid. A few clone orbits remained geocentric for nearly a century, opening the door to the existence of yet-to-be-detected minimoons that are relatively stable for time-scales comparable to those of unbound quasi-satellites such as (469219) Kamo`oalewa. In addition, nearly 10 per cent of the experiments led to brief moon-moon episodes in which the value of the selenocentric energy of 2020 CD3 became negative.

The innermost section of the Solar system has not been extensively studied because minor bodies moving inside Earth's orbit tend to spend most of their sidereal orbital periods at very low solar elongation, well away from the areas more frequently observed by programs searching for near-Earth objects. The survey carried out from the Zwicky Transient Facility (ZTF) is the first one that has been able to detect multiple asteroids well detached from the direct gravitational perturbation of the Earth-Moon system. ZTF discoveries include 2019 AQ3 and 2019 LF6, two Atiras with the shortest periods among known asteroids. Here, we perform an assessment of the orbital evolution of 2020 AV2, an Atira found by ZTF with a similarly short period but following a path contained entirely within the orbit of Venus. This property makes it the first known member of the elusive Vatira population. Genuine Vatiras, those long-term dynamically stable, are thought to be subjected to the so-called von Zeipel-Lidov-Kozai oscillation that protects them against close encounters with both Mercury and Venus. However, 2020 AV2 appears to be a former Atira that entered the Vatira orbital domain relatively recently. It displays an anticoupled oscillation of the values of eccentricity and inclination, but the value of the argument of perihelion may circulate. Simulations show that 2020 AV2 might reach a 3:2 resonant orbit with Venus in the future, activating the von Zeipel-Lidov-Kozai mechanism, which in turn opens the possibility to the existence of a long-term stable population of Vatiras trapped in this configuration.

The chance discovery of the first interstellar minor body, 1I/2017 U1 (`Oumuamua), indicates that we may have been visited by such objects in the past and that these events may repeat in the future. Unfortunately, minor bodies following nearly parabolic or hyperbolic paths tend to receive little attention: over 3/4 of those known have data-arcs shorter than 30 d and, consistently, rather uncertain orbit determinations. This fact suggests that we may have observed interstellar interlopers in the past, but failed to recognize them as such due to insufficient data. Early identification of promising candidates by using N-body simulations may help in improving this situation, triggering follow-up observations before they leave the Solar system. Here, we use this technique to investigate the pre- and post-perihelion dynamical evolution of the slightly hyperbolic comet C/2018 V1 (Machholz-Fujikawa-Iwamoto) to understand its origin and relevance within the context of known parabolic and hyperbolic minor bodies. Based on the available data, our calculations suggest that although C/2018 V1 may be a former member of the Oort Cloud, an origin beyond the Solar system cannot be excluded. If extrasolar, it might have entered the Solar system from interstellar space at low relative velocity with respect to the Sun. The practical feasibility of this alternative scenario has been assessed within the kinematic context of the stellar neighbourhood of the Sun, using data from Gaia second data release, and two robust solar sibling candidates have been identified. Our results suggest that comets coming from interstellar space at low heliocentric velocities may not be rare.

Context. Hypervelocity stars move fast enough to leave the gravitational field of their home galaxies and venture into intergalactic space. The most extreme examples known have estimated speeds in excess of 1000 km/s. These can be easily induced at the centres of galaxies via close encounters between binary stars and supermassive black holes; however, a number of other mechanisms operating elsewhere can produce them as well. Aims. Recent studies suggest that hypervelocity stars are ubiquitous in the local Universe. In the Milky Way, the known hypervelocity stars are anisotropically distributed, but it is unclear why. Here, we used Gaia Data Release 2 (DR2) data to perform a systematic exploration aimed at confirming or refuting these findings. Methods. We assume that the farther the candidate hypervelocity stars are, the more likely they are to be unbound from the Galaxy. We used the statistical analysis of both the spatial distribution and kinematics of these objects to achieve our goals. Results. Focussing on nominal Galactocentric distances greater than 30 kpc, which are the most distant candidates, we isolated a sample with speeds in excess of 500 km/s that exhibits a certain degree of anisotropy but remains compatible with possible systematic effects. We find that the effect of the Eddington-Trumpler-Weaver bias is important in our case: over 80% of our sources are probably located further away than implied by their parallaxes; therefore, most of our velocity estimates are lower limits. If this bias is as strong as suggested here, the contamination by disc stars may not affect our overall conclusions. Conclusions. The subsample with the lowest uncertainties shows stronger, but obviously systematic, anisotropies and includes a number of candidates of possible extragalactic origin and young age with speeds of up to 2000 km/s.

Orbiting the Sun at an average distance of 0.59 au and with the shortest aphelion of any known minor body, at 0.77 au, the Atira-class asteroid 2019 AQ3 may be an orbital outlier or perhaps an early indication of the presence of a new population of objects: those following orbits entirely encompassed within that of Venus, the so-called Vatiras. Here, we explore the orbital evolution of 2019 AQ3 within the context of the known Atiras to show that, like many of them, it displays a present-day conspicuous coupled oscillation of the values of eccentricity and inclination, but no libration of the value of the argument of perihelion with respect to the invariable plane of the Solar system. The observed dynamics is consistent with being the result of the combined action of two dominant perturbers, the Earth-Moon system and Jupiter, and a secondary one, Venus. Such a multiperturber-induced secular dynamics translates into a chaotic evolution that can eventually lead to a resonant behaviour of the Lidov-Kozai type. Asteroid 2019 AQ3 may have experienced brief stints as a Vatira in the relatively recent past and it may become a true Vatira in the future, outlining possible dynamical pathways that may transform Atiras into Vatiras and vice versa. Our results strongly suggest that 2019 AQ3 is only the tip of the iceberg: a likely numerous population of similar bodies may remain hidden in plain sight, permanently confined inside the Sun's glare.

The population of comets hosted by the Oort cloud is heterogeneous. Most studies in this area focused on highly active objects, those with small perihelion distances or examples of objects with peculiar physical properties and/or unusual chemical compositions. This may have produced a biased sample of Oort cloud comets in which the most common objects may be rare, particularly those with perihelia well beyond the orbit of the Earth. Within this context, the known Oort cloud comets may not be representative of the full sample. Here, we study the spectral properties in the visible region and the cometary activity of Comet C/2018 F4 (PANSTARRS). We also explore its orbital evolution with the aim of understanding its origin within the context of known minor bodies moving along nearly parabolic or hyperbolic paths. We present observations obtained with the 10.4 m Gran Telescopio Canarias (GTC), derive the spectral class and visible slope of C/2018 F4 and characterise its level of cometary activity. Direct N-body simulations are carried out to explore its orbital evolution. The absolute magnitude of C/2018 F4 is Hr=13.62+/-0.04. Assuming a pV=0.04 its diameter is D<10.4 km. The object presents a conspicuous coma, with a level of activity comparable to those of other comets observed at similar heliocentric distances. Comet C/2018 F4 has a visible spectrum consistent with that of an X-type asteroid, and has a spectral slope S'=4.0+/-1.0 %/1000\AA and no evidence of hydration. The spectrum matches those of well-studied primitive asteroids and comets. The analysis of its dynamical evolution prior to discovery suggests that C/2018 F4 is not of extrasolar origin. Although the present-day heliocentric orbit of C/2018 F4 is slightly hyperbolic, its observational properties and past orbital evolution are consistent with those of a dynamically old comet with an origin in the Oort cloud.

Context. The discovery and tracking of 2018 LA marks only the third instance in history that the parent body of a fireball has been identified before its eventual disintegration in our atmosphere. The subsequent recovery of meteorites from 2018 LA was only the second time materials from outer space that reached the ground could be linked with certitude to a particular minor body. However, meteoroids like 2018 LA and its forerunners, 2008 TC3 and 2014 AA, are perhaps fragments of larger members of the near-Earth object (NEO) population. As the processes leading to the production of such fragments are unlikely to spawn just one meteoroid per event, it is important to identify putative siblings and plausible candidates from which the observed meteoroids might have originated. Aims. Here, we study the pre-impact orbital evolution of 2018 LA to place this meteoroid within the dynamical context of other NEOs that follow similar trajectories. Methods. Our statistical analyses are based on the results of direct N-body calculations that use the latest orbit determinations and include perturbations by the eight major planets, the Moon, the barycentre of the Pluto-Charon system, and the three largest asteroids. A state-of-the-art NEO orbit model was used to interpret our findings and a randomization test was applied to estimate their statistical significance. Results. We find a statistically significant excess of NEOs in 2018 LA-like orbits; among these objects, we find one impactor, 2018 LA, and the fourth closest known passer-by, 2018 UA. A possible connection with the Chi-Scorpiids meteor shower is also discussed. The largest known NEO with an orbit similar to that of 2018 LA is the potentially hazardous asteroid (454100) 2013 BO73 and we speculate that they both originate from a common precursor via a collisional cascade.

YORP-induced fission events may form dynamically coherent pairs or even families of asteroids. The outcome of this process is well documented among members of the main asteroid belt, but not in the case of the near-Earth asteroid (NEA) population because their paths randomize very efficiently in a short time-scale. Mean-motion resonances (MMRs) may stabilize the orbits of small bodies by making them avoid close encounters with planets. In theory, YORP-induced fission of asteroids trapped in MMRs can preserve evidence of this process even in near-Earth space. Here, we show that two NEAs, 2017 SN16 and 2018 RY7, are currently following an orbital evolution in which their relative mean longitude does not exhibit any secular increase due to the stabilizing action of the 3:5 MMR with Venus. The mechanism that makes this configuration possible may be at work both in the Solar system and elsewhere. Our analysis suggests that the pair 2017 SN16-2018 RY7 may have had its origin in one out of two mechanisms: YORP-induced splitting or binary dissociation.

Stellar streams are ubiquitous in the Galactic halo and they can be used to improve our understanding of the formation and evolution of the Milky Way as a whole. The so-called Monoceros Ring might have been the result of satellite accretion. Guglielmo et al. have used N-body simulations to search for the progenitor of this structure. Their analysis shows that, if the Ring has a dwarf galaxy progenitor, it might be found in the background of one out of eight specific areas in the sky. Here, we use Gaia DR2 data to perform a systematic exploration aimed at confirming or rejecting this remarkable prediction. Focusing on the values of the radial velocity to uncover possible multimodal spreads, we identify a bimodal Gaussian distribution towards Galactic coordinates (l, b) = (271, +2) degrees in Vela, which is one of the locations of the progenitor proposed by Guglielmo et al. This prominent feature with central values 60+/-7 km/s and 97+/-10 km/s, may signal the presence of the long sought progenitor of the Monoceros Ring, but the data might also be compatible with the existence of an unrelated, previously unknown, kinematically coherent structure.

A number of intriguing patterns have been identified in the distributions of the orbital elements of the known extreme trans-Neptunian objects or ETNOs --those with semimajor axis greater than 150 au and perihelion distance greater than 30 au. The announcement of the discovery of 2015 BP519, also known as Caju, has been hailed by some as a consistent piece of robust evidence for the existence of a yet-to-be-discovered massive planet far beyond the trans-Neptunian belt. In this Note, we use the latest data available to provide evidence against 2015 BP519 having followed the same dynamical pathway that placed other ETNOs where they are now. As an statistical outlier within the 29 known ETNOs, 2015 BP519 cannot be used as a reference to further support the trends perhaps present for other ETNOs. Asteroid 2015 BP519 represents a case of extreme dynamical anticorrelation within the ETNO orbital realm and the value of the distance to its descending node, nearly 700 au from the Sun, is the largest of any known ETNO.

Apollo meteoroid 2018 LA has become only the third natural object ever to be discovered prior to causing a meteor airburst and just the second one to have its meteorites recovered (at Botswana's Central Kalahari Game Reserve). Here, we use the latest orbit determination of 2018 LA (solution date 18-July-2018) to search for minor bodies moving in paths comparable to that of 2018 LA using the D-criteria, which are metrics to study orbit similarity, and N-body simulations. Our results further confirm the existence of a dynamical grouping of asteroids that might be related to 2018 LA and show that the impactor could be a recent fragment spawned by a larger object, the 550-m wide, potentially hazardous asteroid (454100) 2013 BO73. Spectroscopic observations of 454100 during its next flyby with our planet (brightest at an apparent visual magnitude of 18.4 on 2018 mid-November) may confirm or deny a putative similar chemical composition to that of the recovered meteorites of 2018 LA.

On 2018 June 2, meteoroid 2018 LA became the third natural body ever to be observed before entering our atmosphere ---small asteroids 2014 AA and 2008 TC3 had stricken the Earth on 2014 January 2 and 2008 October 7, respectively. Here, we explore the pre-impact orbital evolution of 2018 LA and investigate the possible presence of known NEOs moving in similar orbits using N-body simulations and the D-criteria. We identify several objects moving in orbits similar to that of 2018 LA and focus on three of them: (454100) 2013 BO73, which is the largest of the group and a PHA, 2016 LR, and 2018 BA5, which follows a path very close to that of 2018 LA in terms of semimajor axis, eccentricity, and inclination. All these objects could be part of a dynamical grouping and their orbital evolution is rather chaotic, experiencing close encounters with Venus, the Earth-Moon system, and Mars. NEO encounters take place at the node and, on the short-term, the relative positions of our planet in its orbit around the Sun repeat every year. Besides the bright fireball observed over Botswana on 2018 June 2, three other bolides were observed early in June in recent years: Crete on 2002 June 6, Washington State on 2004 June 3, and Reisadalen on 2007 June 7.

Atira-class asteroids have aphelia greater than 0.718 au and smaller than 0.983 au. It has been found that the orbits of many of them remain relatively stable for at least 1 Myr, but it is unclear what can make them so stable when they cross the orbit of Venus and some may reach, in heliocentric terms, as far as Earth's perihelion and as close as inside Mercury's perihelion. The analysis of their orbits shows that many of them have values of the argument of perihelion close to 0 or 180 degrees. Such an orbital arrangement means that the nodes are located at perihelion and at aphelion, i.e. for many known Atiras, away from the path of Venus. Here, we show that the orbits of some Atiras are strongly affected by the Kozai-Lidov resonance, which translates into enhanced long-term stability.

2MASS J06562998+3002455 or PSS 309-6 is a high proper-motion star that was discovered during a survey with the 2.1 m telescope at Kitt Peak National Observatory. Here, we reevaluate the status of this interesting star using Gaia DR2. Our results strongly suggest that PSS 309-6 could be a Population II star as the value of its V component is close to -220 km/s, which is typical for halo stars in the immediate solar neighborhood. Kapteyn's star is the nearest known halo star and PSS 309-6 exhibits similar kinematic and photometric signatures. Its properties also resemble those of 2MASS J15484023-3544254, which was once thought to be the nearest cool white dwarf but was later reclassified as K-type subdwarf. Although it is virtually certain that PSS 309-6 is not a nearby white dwarf but a more distant Population II subdwarf, further spectroscopic information, including radial velocity measurements, is necessary to fully characterize this probable member of the Galactic halo.

LP 543-25 or PSS 544-7 is a high proper-motion star located 458 pc from the Sun in the constellation of Canis Minor; it has been argued that it could be a candidate cannonball star ejected by a star cluster. Here, we revisit the issue of the kinematics of this interesting star using Gaia DR2. The heliocentric Galactic velocity components are (U, V, W) = (206, -289, 30) km/s; the corresponding Galactocentric Galactic velocity components show that LP 543-25 is moving in the Galactic plane and away from the Galactic Center at a rate of nearly 200 km/s, which is compatible with an origin in one of the multiple star clusters that inhabit the inner regions of the Milky Way. LP 543-25 appears to be a member of an elusive class of stars, the low-mass runaway stars. It is perhaps one of the closest and less massive runaway stars identified so far.

Gliese 710 is a K7V star located 19 pc from the Sun in the constellation of Serpens Cauda, which is headed straight for the solar system. Berski & Dybczynski (2016) used data from Gaia DR1 to show that this star will be 13366 AU from the Sun in 1.35 Myr from now. Here, we present an independent confirmation of this remarkable result using Gaia DR2. Our approach is first validated using as test case that of the closest known stellar flyby, by the binary WISE J072003.20-084651.2 or Scholz's star. Our results confirm, within errors, those in Berski & Dybczynski (2016), but suggest a somewhat closer, both in terms of distance and time, flyby of Gliese 710 to the solar system. Such an interaction might not significantly affect the region inside 40 au as the gravitational coupling among the known planets against external perturbation can absorb efficiently such a perturbation, but it may trigger a major comet shower that will affect the inner solar system.

At discovery time, C/2017 K2 (PANSTARRS) was the second most distant inbound active comet ever observed. It has been argued that this object is in the process of crossing the inner Solar System for the first time, but other authors have concluded that it is dynamically old. We have performed full N-body simulations for 3 Myr into the past using the latest public orbit determination for this object and most of them, 67%, are consistent with a bound and dynamically old Oort cloud comet, but about 29% of the studied orbits are compatible with an interstellar origin. Our independent calculations strongly suggest that C/2017 K2 is not a dynamically new Oort cloud comet.

Observed hyperbolic minor bodies might have an interstellar origin, but they can be natives of the Solar system as well. Fly-bys with the known planets or the Sun may result in the hyperbolic ejection of an originally bound minor body; in addition, members of the Oort cloud could be forced to follow inbound hyperbolic paths as a result of secular perturbations induced by the Galactic disc or, less frequently, due to impulsive interactions with passing stars. These four processes must leave distinctive signatures in the distribution of radiants of observed hyperbolic objects, both in terms of coordinates and velocity. Here, we perform a systematic numerical exploration of the past orbital evolution of known hyperbolic minor bodies using a full N-body approach and statistical analyses to study their radiants. Our results confirm the theoretical expectations that strong anisotropies are present in the data. We also identify a statistically significant overdensity of high-speed radiants towards the constellation of Gemini that could be due to the closest and most recent known fly-by of a star to the Solar system, that of the so-called Scholz's star. In addition to and besides 1I/2017 U1 (`Oumuamua), we single out eight candidate interstellar comets based on their radiants' velocities.

The present-day orbit of 2017 UV43 resembles the average one of the fragments of comet Shoemaker-Levy 9. Here, we present numerical evidence suggesting that 2017 UV43 might be related to Shoemaker-Levy 9. With our exploratory results we intend to encourage a search for precovery images of 2017 UV43 and perhaps even follow-up observations that may help in improving its poorly determined orbit so its true dynamical nature can be understood.

The organization of the orbits of most minor bodies in the Solar system seems to follow random patterns, the result of billions of years of chaotic dynamical evolution. Much as heterogeneous orbital behaviour is ubiquitous, dynamically coherent pairs and groups of objects are also present everywhere. Although first studied among the populations of asteroids and comets that inhabit or traverse the inner Solar system, where they are very numerous, at least one asteroid family has been confirmed to exist in the outer Solar system and two other candidates have been proposed in the literature. Here, we perform a systematic search for statistically significant pairs and groups of dynamically correlated objects through those with semimajor axis greater than 25 au, applying a novel technique that uses the angular separations of orbital poles and perihelia together with the differences in time of perihelion passage to single out pairs of relevant objects. Our analysis recovers well-known, dynamically coherent pairs and groups of comets and trans-Neptunian objects and uncovers a number of new ones, prime candidates for further spectroscopic study.

The discovery of 1991 VG on 1991 November 6 attracted an unprecedented amount of attention as it was the first near-Earth object (NEO) ever found on an Earth-like orbit. At that time, it was considered by some as the first representative of a new dynamical class of asteroids, while others argued that an artificial (terrestrial or extraterrestrial) origin was more likely. Over a quarter of a century later, this peculiar NEO has been recently recovered and the new data may help in confirming or ruling out early theories about its origin. Here, we use the latest data to perform an independent assessment of its current dynamical status and short-term orbital evolution. Extensive N-body simulations show that its orbit is chaotic on time-scales longer than a few decades. We confirm that 1991 VG was briefly captured by Earth's gravity as a minimoon during its previous fly-by in 1991-1992; although it has been a recurrent transient co-orbital of the horseshoe type in the past and it will return as such in the future, it is not a present-day co-orbital companion of the Earth. A realistic NEO orbital model predicts that objects like 1991 VG must exist and, consistently, we have found three other NEOs --2001 GP2, 2008 UA202 and 2014 WA366-- which are dynamically similar to 1991 VG. All this evidence confirms that there is no compelling reason to believe that 1991 VG is not natural.

Here, we explore some peculiar orbital features of the recently discovered asteroid A/2017 U1, which is a clear outlier when considering the average value of the eccentricity of known hyperbolic comets. As for the orientation of its orbit in space, the orbital plane of A/2017 U1 seems to be away from any obvious clusters present for this population. The orbital nodes of A/2017 U1 are well away from the paths of the planets of the Solar System and the Sun. All these orbital properties appear to confirm A/2017 U1 as the first known interstellar asteroid.

Venus has no known satellites, but has four known co-orbitals: (322756) 2001 CK32, 2002 VE68, 2012 XE133, and 2013 ND15. Here, we present numerical evidence suggesting that 2015 WZ12 is a possible Venus co-orbital; it might have been until recently a transient Trojan. Follow-up observations of this target in the near future will be difficult, though.

The first direct detection of the asteroidal YORP effect, a phenomenon that changes the spin states of small bodies due to thermal reemission of sunlight from their surfaces, was obtained for (54509) YORP 2000 PH5. Such an alteration can slowly increase the rotation rate of asteroids, driving them to reach their fission limit and causing their disruption. This process can produce binaries and unbound asteroid pairs. Secondary fission opens the door to the eventual formation of transient but genetically-related groupings. Here, we show that the small near-Earth asteroid (NEA) 2017 FZ2 was a co-orbital of our planet of the quasi-satellite type prior to their close encounter on 2017 March 23. Because of this flyby with the Earth, 2017 FZ2 has become a non-resonant NEA. Our N-body simulations indicate that this object may have experienced quasi-satellite engagements with our planet in the past and it may return as a co-orbital in the future. We identify a number of NEAs that follow similar paths, the largest named being YORP, which is also an Earth's co-orbital. An apparent excess of NEAs moving in these peculiar orbits is studied within the framework of two orbit population models. A possibility that emerges from this analysis is that such an excess, if real, could be the result of mass shedding from YORP itself or a putative larger object that produced YORP. Future spectroscopic observations of 2017 FZ2 during its next visit in 2018 (and of related objects when feasible) may be able to confirm or reject this interpretation.

Asteroids that follow similar orbits may have a dynamical connection as their current paths could be the result of a past interaction with a massive perturber. The pair of extreme trans-Neptunian objects or ETNOs (474640) 2004 VN112-2013 RF98 exhibits peculiar relative orbital properties, including a difference in longitude of the ascending node of just 1.61 degrees and 3.99 degrees in inclination. In addition, their reflectance spectra are similar in the visible portion of the spectrum. The origin of these similarities remains unclear. Neglecting observational bias, viable scenarios that could explain this level of coincidence include fragmentation and binary dissociation. Here, we present results of extensive direct N-body simulations of close encounters between wide binary ETNOs and one trans-Plutonian planet. We find that wide binary ETNOs can dissociate during such interactions and the relative orbital properties of the resulting unbound couples match reasonably well those of several pairs of known ETNOs, including 474640-2013 RF98. The possible presence of former binaries among the known ETNOs has strong implications for the interpretation of the observed anisotropies in the distributions of the directions of their orbital poles and perihelia.

It is a well-known fact that the presence of a massive perturber interacting with a population of minor bodies following very eccentric orbits can strongly affect the distribution of their nodal distances. The details of this process have been explored numerically and its outcome confirmed observationally in the case of Jupiter, where a bimodal distribution of nodal distances of comets has been found. Here, we show evidence for a possible bimodal distribution of the nodal distances of the extreme trans-Neptunian objects (ETNOs) in the form of a previously unnoticed correlation between nodal distance and orbital inclination. This proposed correlation is unlikely to be the result of observational bias as data for both large semimajor axis Centaurs and comets fit well into the pattern found for the ETNOs, and all these populations are subjected to similar background perturbations when moving well away from the influence of the giant planets. The correlation found is better understood if these objects tend to avoid a putative planet with semimajor axis in the range 300-400 au.

In the outer Solar system, primordial Trojan asteroids may have remained dynamically stable for billions of years. Several thousands of them accompany Jupiter in its journey around the Sun and a similarly large population may be hosted by Neptune. In addition, recently captured or transient Jovian and Neptunian Trojans are not uncommon. In contrast, no Trojans of Saturn have been found yet and just one Uranian Trojan is known, 2011 QF99. Here, we discuss the identification of a second Trojan of Uranus: 2014 YX49. Like 2011 QF99, 2014 YX49 is a transient L4 Trojan although it orbits at higher inclination (25.55 degrees versus 10.83 degrees), is larger (absolute magnitude of 8.5 versus 9.7) and its libration period is slightly shorter (5.1 versus 5.9 kyr); contrary to 2011 QF99, its discovery was not the result of a targeted survey. It is less stable than 2011 QF99; our extensive N-body simulations show that 2014 YX49 may have been following a tadpole trajectory ahead of Uranus for about 60 kyr and it can continue doing so for another 80 kyr. Our analysis suggests that it may remain as co-orbital for nearly 1 Myr. As in the case of 2011 QF99, the long-term stability of 2014 YX49 is controlled by Jupiter and Neptune, but it is currently trapped in the 7:20 mean motion resonance with Saturn. Consistently, the dynamical mechanism leading to the capture into and the ejection from the Trojan state involves ephemeral multibody mean motion resonances.

The existence of significant anisotropies in the distributions of the directions of perihelia and orbital poles of the known extreme trans-Neptunian objects (ETNOs) has been used to claim that trans-Plutonian planets may exist. Among the known ETNOs, the pair (474640) 2004 VN112-2013 RF98 stands out. Their orbital poles and the directions of their perihelia and their velocities at perihelion/aphelion are separated by a few degrees, but orbital similarity does not necessarily imply common physical origin. In an attempt to unravel their physical nature, visible spectroscopy of both targets was obtained using the OSIRIS camera-spectrograph at the 10.4 m Gran Telescopio Canarias (GTC). From the spectral analysis, we find that 474640-2013 RF98 have similar spectral slopes (12 vs. 15 %/0.1um), very different from Sedna's but compatible with those of (148209) 2000 CR105 and 2012 VP113. These five ETNOs belong to the group of seven linked to the Planet Nine hypothesis. A dynamical pathway consistent with these findings is dissociation of a binary asteroid during a close encounter with a planet and we confirm its plausibility using N-body simulations. We thus conclude that both the dynamical and spectroscopic properties of 474640-2013 RF98 favour a genetic link and their current orbits suggest that the pair was kicked by a perturber near aphelion.

On 2008 October 7, small asteroid 2008 TC3 turned itself into the parent body of the first meteor ever to be predicted before entering the Earth's atmosphere. Over five years later, the 2014 AA event became the second instance of such an occurrence. The uncertainties associated with the pre-impact orbit of 2008 TC3 are relatively small because thousands of observations were made during the hours preceding the actual meteor airburst. In sharp contrast, 2014 AA was only observed seven times before impact and consequently its trajectory is somewhat uncertain. Here, we present a recalculation of the impact parameters -location and timing- of this meteor based on infrasound recordings. The new values and their uncertainties together with Monte Carlo and N-body techniques, are applied to obtain an independent determination of the pre-impact orbit of 2014 AA. Our orbital solution is used to investigate the possible presence of known near-Earth objects (NEOs) moving in similar orbits. Prior to impact, 2014 AA was subjected to a web of overlapping secular resonances and it followed a path similar to those of 2011 GJ3, 2011 JV10, 2012 DJ54, and 2013 NJ4. NEOs in this transient group have their orbits controlled by close encounters with the Earth-Moon system at perihelion and Mars at aphelion, perhaps constituting a dynamical family. Extensive comparison with other studies is also presented.

A number of Earth co-orbital asteroids experience repeated transitions between the quasi-satellite and horseshoe dynamical states. Asteroids 2001 GO2, 2002 AA29, 2003 YN107 and 2015 SO2 are well-documented cases of such a dynamical behaviour. These transitions depend on the gravitational influence of other planets, owing to the overlapping of a multiplicity of secular resonances. Here, we show that the recently discovered asteroid (469219) 2016 HO3 is a quasi-satellite of our planet -the fifth one, joining the ranks of (164207) 2004 GU9, (277810) 2006 FV35, 2013 LX28 and 2014 OL339. This new Earth co-orbital also switches repeatedly between the quasi-satellite and horseshoe configurations. Its current quasi-satellite episode started nearly 100 yr ago and it will end in about 300 yr from now. The orbital solution currently available for this object is very robust and our full N-body calculations show that it may be a long-term companion (time-scale of Myr) to our planet. Among the known Earth quasi-satellites, it is the closest to our planet and as such, a potentially accessible target for future in situ study. Due to its presumably lengthy dynamical relationship with the Earth and given the fact that at present and for many decades this transient object remains well positioned with respect to our planet, the results of spectroscopic studies of this small body, 26-115 m, may be particularly useful to improve our understanding of the origins -local or captured- of Earth's co-orbital asteroid population. The non-negligible effect of the uncertainty in the value of the mass of Jupiter on the stability of this type of co-orbitals is also briefly explored.

In the Solar system, a quasi-satellite is an object that follows a heliocentric path with an orbital period that matches almost exactly with that of a host body (planetary or not). The trajectory is of such nature that, without being gravitationally attached, the value of the angular separation between host and quasi-satellite as seen from the Sun remains confined within relatively narrow limits for time-spans that exceed the length of the host's sidereal orbital period. Here, we show that under these conditions, a quasi-satellite traces an analemma in the sky as observed from the host in a manner similar to that found for geosynchronous orbits. The analemmatic curve (figure-eight-, teardrop-, ellipse-shaped) results from the interplay between the tilt of the rotational axis of the host and the properties of the orbit of the quasi-satellite. The analemma criterion can be applied to identify true quasi-satellite dynamical behaviour using observational or synthetic astrometry and it is tested for several well-documented quasi-satellites. For the particular case of 15810 (1994 JR1), a putative accidental quasi-satellite of dwarf planet Pluto, we show explicitly that this object describes a complex analemmatic curve for several Plutonian sidereal periods, confirming its transient quasi-satellite status.

The distribution of the orbital elements of the known extreme trans-Neptunian objects or ETNOs has been found to be statistically incompatible with that of an unperturbed asteroid population following heliocentric or, better, barycentric orbits. Such trends, if confirmed by future discoveries of ETNOs, strongly suggest that one or more massive perturbers could be located well beyond Pluto. Within the trans-Plutonian planets paradigm, the Planet Nine hypothesis has received much attention as a robust scenario to explain the observed clustering in physical space of the perihelia of seven ETNOs which also exhibit clustering in orbital pole position. Here, we revisit the subject of clustering in perihelia and poles of the known ETNOs using barycentric orbits, and study the visibility of the latest incarnation of the orbit of Planet Nine applying Monte Carlo techniques and focusing on the effects of the apsidal anti-alignment constraint. We provide visibility maps indicating the most likely location of this putative planet if it is near aphelion. We also show that the available data suggest that at least two massive perturbers are present beyond Pluto.

Both Centaurs and trans-Neptunian objects (TNOs) are minor bodies found in the outer Solar System. Centaurs are a transient population that moves between the orbits of Jupiter and Neptune, and they probably diffused out of the TNOs. TNOs move mainly beyond Neptune. Some of these objects display episodic cometary behaviour; a few percent of them are known to host binary companions. Here, we study the light-curves of two Centaurs -2060 Chiron (1977 UB) and 10199 Chariklo (1997 CU26)- and three TNOs -38628 Huya (2000 EB173), 28978 Ixion (2001 KX76), and 90482 Orcus (2004 DW)- and the colours of the Centaurs and Huya. Precise, ~1%, R-band absolute CCD photometry of these minor bodies acquired between 2006 and 2011 is presented; the new data are used to investigate the rotation rate of these objects. The colours of the Centaurs and Huya are determined using BVRI photometry. The point spread function of the five minor bodies is analysed, searching for signs of a coma or close companions. Astrometry is also discussed. A periodogram analysis of the light-curves of these objects gives the following rotational periods: 5.5+-0.4 h for Chiron, 7.0+-0.6 h for Chariklo, 4.45+-0.07 h for Huya, 12.4+-0.3 h for Ixion, and 11.9+-0.5 h for Orcus. The colour indices of Chiron are found to be B-V=0.53+-0.05, V-R=0.37+-0.08, and R-I=0.36+-0.15. The values computed for Chariklo are V-R=0.62+-0.07 and R-I=0.61+-0.07. For Huya, we find V-R=0.58+-0.09 and R-I=0.64+-0.20. We find very low levels of cometary activity (if any) and no sign of close or wide binary companions for these minor bodies.

The Planet Nine hypothesis has now enough constraints to deserve further attention in the form of detailed numerical experiments. The results of such studies can help us improve our understanding of the dynamical effects of such a hypothetical object on the extreme trans-Neptunian objects or ETNOs and perhaps provide additional constraints on the orbit of Planet Nine itself. Here, we present the results of direct N-body calculations including the latest data available on the Planet Nine conjecture. The present-day orbits of the six ETNOs originally linked to the hypothesis are evolved backwards in time and into the future under some plausible incarnations of the hypothesis to investigate if the values of several orbital elements, including the argument of perihelion, remain confined to relatively narrow ranges. We find that a nominal Planet Nine can keep the orbits of (90377) Sedna and 2012 VP113 relatively well confined in orbital parameter space for hundreds of Myr, but it may make the orbits of 2004 VN112, 2007 TG422 and 2013 RF98 very unstable on time-scales of dozens of Myr, turning them retrograde and eventually triggering their ejection from the Solar system. Far more stable orbital evolution is found with slightly modified orbits for Planet Nine.

Many asteroids in the main and trans-Neptunian belts are trapped in mean motion resonances with Jupiter and Neptune, respectively. As a side effect, they experience accidental commensurabilities among themselves. These commensurabilities define characteristic patterns that can be used to trace the source of the observed resonant behaviour. Here, we explore systematically the existence of commensurabilities between the known ETNOs using their heliocentric and barycentric semimajor axes, their uncertainties, and Monte Carlo techniques. We find that the commensurability patterns present in the known ETNO population resemble those found in the main and trans-Neptunian belts. Although based on small number statistics, such patterns can only be properly explained if most, if not all, of the known ETNOs are subjected to the resonant gravitational perturbations of yet undetected trans-Plutonian planets. We show explicitly that some of the statistically significant commensurabilities are compatible with the Planet Nine hypothesis; in particular, a number of objects may be trapped in the 5:3 and 3:1 mean motion resonances with a putative Planet Nine with semimajor axis of about 700 au.

Planet Nine is a hypothetical planet located well beyond Pluto that has been proposed in an attempt to explain the observed clustering in physical space of the perihelia of six extreme trans-Neptunian objects or ETNOs. The predicted approximate values of its orbital elements include a semimajor axis of 700 au, an eccentricity of 0.6, an inclination of 30 degrees, and an argument of perihelion of 150 degrees. Searching for this putative planet is already under way. Here, we use a Monte Carlo approach to create a synthetic population of Planet Nine orbits and study its visibility statistically in terms of various parameters and focusing on the aphelion configuration. Our analysis shows that, if Planet Nine exists and is at aphelion, it might be found projected against one out of four specific areas in the sky. Each area is linked to a particular value of the longitude of the ascending node and two of them are compatible with an apsidal anti-alignment scenario. In addition and after studying the current statistics of ETNOs, a cautionary note on the robustness of the perihelia clustering is presented.

It is widely accepted that a quasi-steady-state flux of minor bodies moving in and out of the co-orbital state with the Earth may exist. Some of these objects are very good candidates for future in situ study due to their favourable dynamical properties. In this paper, we show that the recently discovered near-Earth asteroids 2015 XX169, 2015 YA and 2015 YQ1 are small transient Earth co-orbitals. These new findings increase the tally of known Earth co-orbitals to 17. The three of them currently exhibit asymmetric horseshoe behaviour subjected to a Kozai resonance and their short-term orbital evolution is rather unstable. Both 2015 YA and 2015 YQ1 may leave Earth's co-orbital zone in the near future as they experience close encounters with Venus, the Earth-Moon system and Mars. Asteroid 2015 XX169 may have remained in the vicinity of, or trapped inside, the 1:1 mean motion resonance with our planet for many thousands of years and may continue in that region for a significant amount of time into the future.

Among the near-Earth object (NEO) population there are comets and active asteroids which are sources of fragments that initially move together; in addition, some NEOs follow orbits temporarily trapped in a web of secular resonances. These facts contribute to increasing the risk of meteoroid strikes on Earth, making its proper quantification difficult. The identification and subsequent study of groups of small NEOs that appear to move in similar trajectories are necessary steps in improving our understanding of the impact risk associated with meteoroids. Here, we present results of a search for statistically significant dynamical groupings among the NEO population. Our Monte Carlo-based methodology recovers well-documented groupings like the Taurid Complex or the one resulting from the split comet 73P/Schwassmann-Wachmann 3, and new ones that may have been the source of past impacts. Among the most conspicuous are the Mjolnir and Ptah groups, perhaps the source of recent impact events like Almahata Sitta and Chelyabinsk, respectively. Meteoroid 2014 AA, that hit the Earth on 2014 January 2, could have its origin in a marginally significant grouping associated with Bennu. We find that most of the substructure present within the orbital domain of the NEOs is of resonant nature, probably induced by secular resonances and the Kozai mechanism that confine these objects into specific paths with well-defined perihelia.

Earth co-orbitals of the horseshoe type are interesting objects to study for practical reasons. They are relatively easy to access from our planet and that makes them attractive targets for sample return missions. Here, we show that near-Earth asteroid (NEA) 2015 SO2 is a transient co-orbital to the Earth that experiences a rather peculiar orbital evolution characterised by recurrent, alternating horseshoe and quasi-satellite episodes. It is currently following a horseshoe trajectory, the ninth asteroid known to do so. Besides moving inside the 1:1 mean motion resonance with the Earth, it is subjected to a Kozai resonance with the value of the argument of perihelion librating around 270 degrees. Contrary to other NEAs, asteroid 2015 SO2 may have remained in the vicinity of Earth's co-orbital region for a few hundreds of thousands of years.

On 2013 February 15 a small asteroid rammed against the atmosphere above the region of Chelyabinsk in Russia, producing the most powerful superbolide since the Tunguska event in 1908. Lacking proper astrometric observations, the pre-impact orbit of this object has been determined using videos, satellite images, and pure geometry. Unfortunately, more than two years after the event, the published estimates vary so much that there is no clear orbital solution that could be used to investigate the origin of the impactor and the existence of dynamically, or perhaps even genetically, related asteroids. Here, we revisit this topic using a full N-body approach. A robust statistical test is applied to published solutions to discard those unable to produce a virtual impact at the observed time (03:20:20.8 s UTC). The same N-body methodology and the latest ephemerides are used to compute a new orbital solution: a=1.6247 AU, e=0.5318, i=3.9750 degrees, Omega=326.4607 degrees and omega=109.7012 degrees. This new solution --which has an impact probability > 0.99999 and uncertainties in time and space of 0.2 s and 6 km, respectively-- is utilized to explore the past orbital evolution of the impactor as well as the presence of near-Earth objects moving in similar paths. A dynamical link between asteroid 2011 EO40 and the Chelyabinsk impactor is confirmed. Alternative orbital solutions are extensively explored.

The detection of high-energy astrophysical neutrinos of extraterrestrial origin by the IceCube neutrino observatory in Antarctica has opened a unique window to the cosmos that may help to probe both the distant Universe and our cosmic backyard. The arrival directions of these high-energy events have been interpreted as uniformly distributed on the celestial sphere. Here, we revisit the topic of the putative isotropic angular distribution of these events applying Monte Carlo techniques to investigate a possible anisotropy. A modest evidence for anisotropy is found. An excess of events appears projected towards a section of the Local Void, where the density of galaxies with radial velocities below 3000 km/s is rather low, suggesting that this particular group of somewhat clustered sources are located either very close to the Milky Way or perhaps beyond 40 Mpc. The results of further analyses of the subsample of southern hemisphere events favour an origin at cosmological distances with the arrival directions of the events organized in a fractal-like structure. Although a small fraction of closer sources is possible, remote hierarchical structures appear to be the main source of these very energetic neutrinos. Some of the events may have their origin at the IBEX ribbon.

Minor bodies trapped in 1:1 co-orbital resonances with a host planet could be relevant to explain the origin of captured satellites. Among the giant planets, Uranus has one of the smallest known populations of co-orbitals, three objects, and all of them are short-lived. Asteroid 2015 DB216 has an orbital period that matches well that of Uranus, and here we investigate its dynamical state. Direct N-body calculations are used to assess the current status of this object, reconstruct its immediate dynamical past, and explore its future orbital evolution. A covariance matrix-based Monte Carlo scheme is presented and applied to study its short-term stability. We find that 2015 DB216 is trapped in a temporary co-orbital resonance with Uranus, the fourth known minor body to do so. A detailed analysis of its dynamical evolution shows that it is an unstable but recurring co-orbital companion to Uranus. It currently follows an asymmetric horseshoe trajectory that will last for at least 10 kyr, but it may remain inside Uranus' co-orbital zone for millions of years. As in the case of other transient Uranian co-orbitals, complex multibody ephemeral mean motion resonances trigger the switching between the various resonant co-orbital states. The new Uranian co-orbital exhibits a secular behaviour markedly different from that of the other known Uranian co-orbitals because of its higher inclination, nearly 38 degrees. Given its rather unusual discovery circumstances, the presence of 2015 DB216 hints at the existence of a relatively large population of objects moving in similar orbits.

Context. In the Milky Way, most globular clusters are highly conspicuous objects that were found centuries ago. However, a few dozen of them are faint, sparsely populated systems that were identified largely during the second half of the past century. One of the faintest is ESO 37-1 (E 3) and as such it remains poorly studied, with no spectroscopic observations published so far, although it was discovered in 1976. Aims. We investigate the globular cluster E 3 in an attempt to better constrain its fundamental parameters. Spectroscopy of stars in the field of E 3 is shown here for the first time. Methods. Deep, precise VI CCD photometry of E 3 down to V=26 mag is presented and analysed. Low-resolution, medium signal-to-noise ratio spectra of nine candidate members are studied to derive radial velocity and metallicity. Proper motions from the UCAC4 catalogue are used to explore the kinematics of the bright members of E 3. Results. Isochrone fitting indicates that E 3 is probably very old, with an age of about 13 Gyr; its distance from the Sun is nearly 10 kpc. It is also somewhat metal rich with [Fe/H]=-0.7. Regarding its kinematics, our tentative estimate for the proper motions is (-7.0+/-0.8, 3.5+/-0.3) mas/yr (or a tangential velocity of 382+/-79 km/s) and for the radial velocity is 45+/-5 km/s, in the solar rest frame. Conclusions. E 3 is one of the most intriguing globular clusters in the Galaxy. Having an old age and being metal rich is clearly a peculiar combination, only seen in a handful of objects like the far more conspicuous NGC 104 (47 Tucanae). In addition, its low luminosity and sparse population make it a unique template for the study of the final evolutionary phases in the life of a star cluster. Unfortunately, E 3 is among the most elusive and challenging known globular clusters because field contamination severely hampers spectroscopic studies.

Context. Asteroids with semi-major axes very close to that of a host planet can avoid node crossings when their nodal points are at perihelion and at aphelion. This layout protects the asteroids from close encounters, and eventual collisions, with the host planet. Aims. Here, we study the short-term dynamical evolution of four recently discovered near-Earth asteroids (NEAs) --2012 FC71, 2014 EK24, 2014 QD364, and 2014 UR-- that follow very Earth-like orbits. Methods. Our analysis is based on results of direct N-body calculations that use the most updated ephemerides and include perturbations from the eight major planets, the Moon, the barycentre of the Pluto-Charon system, and the three largest asteroids. Results. These four NEAs exhibit an orbital evolution unlike any other known near-Earth object (NEO). Beyond horseshoe, tadpole, or quasi-satellite trajectories, they follow co-orbital passing orbits relative to the Earth within the Kozai domain. Our calculations show that secular interactions induce librations of their relative argument of perihelion with respect to our planet but also to Venus, Mars, and Jupiter. Secular chaos is also present. The size of this transient population is probably large. Conclusions. Although some of these NEAs can remain orbitally stable for many thousands of years, their secular dynamics are substantially more complicated than commonly thought and cannot be properly described within the framework of the three-body problem alone owing to the overlapping of multiple secular resonances. Objects in this group are amongst the most atypical NEOs regarding favourable visibility windows because these are separated in time by many decades or even several centuries.

Arjuna-type orbits are characterized by being Earth-like, having both low-eccentricity and low-inclination. Objects following these trajectories experience repeated trappings in the 1:1 commensurability with the Earth and can become temporary Trojans, horseshoe librators, quasi-satellites, and even transient natural satellites. Here, we review what we know about this peculiar dynamical group and use a Monte Carlo simulation to characterize geometrically the Arjuna orbital domain, studying its visibility both from the ground and with the European Space Agency Gaia spacecraft. The visibility analysis from the ground together with the discovery circumstances of known objects are used as proxies to estimate the current size of this population. The impact cross-section of the Earth for minor bodies in this resonant group is also investigated. We find that, for ground-based observations, the solar elongation at perigee of nearly half of these objects is less than 90 degrees. They are best observed by space-borne telescopes, but Gaia is not going to improve significantly the current discovery rate for members of this class. Our results suggest that the size of this population may have been underestimated by current models. On the other hand, their intrinsically low encounter velocities with the Earth induce a 10-1000-fold increase in the impact cross-section with respect to what is typical for objects in the Apollo or Aten asteroid populations. We estimate that their probability of capture as transient natural satellites of our planet is about 8%.

Nearly all known extreme trans-Neptunian objects (ETNOs) have argument of perihelion close to 0 degrees. An existing observational bias strongly favours the detection of ETNOs with arguments of perihelion close to 0 degrees and 180 degrees yet no objects have been found at 180 degrees. No plausible explanation has been offered so far to account for this unusual pattern. Here, we study the dynamical evolution of comet 96P/Machholz 1, a bizarre near-Earth object (NEO) that may provide the key to explain the puzzling clustering of orbits around argument of perihelion close to 0 degrees recently found for the population of ETNOs. Comet 96P/Machholz 1 is currently locked in a Kozai resonance with Jupiter such that the value of its argument of perihelion is always close to 0 degrees at its shortest possible perihelion (highest eccentricity and lowest inclination) and about 180 degrees near its shortest aphelion (longest perihelion distance, lowest eccentricity and highest inclination). If this object is a dynamical analogue (albeit limited) of the known ETNOs, this implies that massive perturbers must keep them confined in orbital parameter space. Besides, its future dynamical evolution displays orbital flips when its eccentricity is excited to a high value and its orbit turns over by nearly 180 degrees, rolling over its major axis. This unusual behaviour, that is preserved when post-Newtonian terms are included in the numerical integrations, may also help understand the production of NEOs on retrograde orbits.

It is customarily assumed that Earth-striking meteoroids are completely random, and that all the impacts must be interpreted as uncorrelated events distributed according to Poisson statistics. If this is correct, their impact dates must be uniformly spread throughout the year and their impact coordinates must be evenly scattered on the surface of our planet. Here, we use a time- and yield-limited sample of Earth-impacting superbolides detected since 2000 to explore statistically this critical though frequently overlooked topic. We show that the cadence of these multi-kiloton impact events is incompatible with a random fall pattern at the 0.05 significance level or better. This result is statistically robust and consistent with the observed distribution of the longitudes of the ascending nodes of near-Earth objects (NEOs). This lack of randomness is induced by planetary perturbations, in particular Jupiter's, and suggests that some of the recent, most powerful Earth impacts may be associated with resonant groups of NEOs and/or very young meteoroid streams. An intriguing consequence of this scenario is that the impact hazard of Chelyabinsk-like objects should peak at certain times in the year.

Our planet has one permanently bound satellite -the Moon-, a likely large number of mini-moons or transient irregular natural satellites, and three temporary natural retrograde satellites or quasi-satellites. These quasi-moons -(164207) 2004 GU9, (277810) 2006 FV35 and 2013 LX28- are unbound companions to the Earth. The orbital evolution of quasi-satellites may transform them into temporarily bound satellites of our planet. Here, we study the dynamical evolution of the recently discovered Aten asteroid 2014 OL339 to show that it is currently following a quasi-satellite orbit with respect to the Earth. This episode started at least about 775 yr ago and it will end 165 yr from now. The orbit of this object is quite chaotic and together with 164207 are the most unstable of the known Earth quasi-satellites. This group of minor bodies is, dynamically speaking, very heterogeneous but three of them exhibit Kozai-like dynamics: the argument of perihelion of 164207 oscillates around -90 degrees, the one of 277810 librates around 180 degrees and that of 2013 LX28 remains around 0 degrees. Asteroid 2014 OL339 is not currently engaged in any Kozai-like dynamics.

The existence of an outer planet beyond Pluto has been a matter of debate for decades and the recent discovery of 2012 VP113 has just revived the interest for this controversial topic. This Sedna-like object has the most distant perihelion of any known minor planet and the value of its argument of perihelion is close to 0 degrees. This property appears to be shared by almost all known asteroids with semimajor axis greater than 150 au and perihelion greater than 30 au (the extreme trans-Neptunian objects or ETNOs), and this fact has been interpreted as evidence for the existence of a super-Earth at 250 au. In this scenario, a population of stable asteroids may be shepherded by a distant, undiscovered planet larger than the Earth that keeps the value of their argument of perihelion librating around 0 degrees as a result of the Kozai mechanism. Here, we study the visibility of these ETNOs and confirm that the observed excess of objects reaching perihelion near the ascending node cannot be explained in terms of any observational biases. This excess must be a true feature of this population and its possible origin is explored in the framework of the Kozai effect. The analysis of several possible scenarios strongly suggest that at least two trans-Plutonian planets must exist.

The Chelyabinsk superbolide was the largest known natural object to enter the Earth's atmosphere since the Tunguska event in 1908 and it has become a template to understand, manage and mitigate future impacts. Although the event has been documented in great detail, the actual pre-impact orbit of the parent body is still controversial. Here, we revisit this topic using an improved Monte Carlo approach that includes the coordinates of the impact point to compute the most probable solution for the pre-impact orbit (a = 1.62 au, e = 0.53, i = 3.97 degrees, \Omega; = 326.45 and \omega; = 109.71). We also check all the published solutions using a simple yet robust statistical test to show that many of them have problems to cause an impact at the right time. We use the improved orbit and N-body simulations to revisit the dynamical status of a putative Chelyabinsk asteroid family and confirm that it could be linked to resonant asteroids 2007 BD7 and 2011 EO40. In addition, and as the classification of Chelyabinsk meteorites is well established, a search for meteorite falls of the same chondrite group and petrologic type gives some evidence for the existence of an associated LL5 chondrite cluster.

Among all the asteroid dynamical groups, Centaurs have the highest fraction of objects moving in retrograde orbits. The distribution in absolute magnitude, H, of known retrograde Centaurs with semi-major axes in the range 6-34 AU exhibits a remarkable trend: 10% have H < 10 mag, the rest have H > 12 mag. The largest objects, namely (342842) 2008 YB3, 2011 MM4 and 2013 LU28, move in almost polar, very eccentric paths; their nodal points are currently located near perihelion and aphelion. In the group of retrograde Centaurs, they are obvious outliers both in terms of dynamics and size. Here, we show that these objects are also trapped in retrograde resonances that make them unstable. Asteroid 2013 LU28, the largest, is a candidate transient co-orbital to Uranus and it may be a recent visitor from the trans-Neptunian region. Asteroids 342842 and 2011 MM4 are temporarily submitted to various high-order retrograde resonances with the Jovian planets but 342842 may be ejected towards the trans-Neptunian region within the next few hundred kyr. Asteroid 2011 MM4 is far more stable. Our analysis shows that the large retrograde Centaurs form an heterogeneous group that may include objects from various sources. Asteroid 2011 MM4 could be a visitor from the Oort cloud but an origin in a relatively stable closer reservoir cannot be ruled out. Minor bodies like 2011 MM4 may represent the remnants of the primordial planetesimals and signal the size threshold for catastrophic collisions in the early Solar System.

Uranus has three known co-orbitals: 83982 Crantor (2002 GO9), 2010 EU65 and 2011 QF99. Here, we perform a comparative analysis of the orbital evolution of these transient co-orbitals to understand better how they got captured in the first place and what makes them dynamically unstable. We also look for additional temporary Uranian co-orbital candidates among known objects. Our N-body simulations show that the long-term stability of 2011 QF99 is controlled by Jupiter and Neptune; it briefly enters the 1:7 mean motion resonance with Jupiter and the 2:1 with Neptune before becoming a Trojan and prior to leaving its tadpole orbit. During these ephemeral two-body mean motion resonance episodes, apsidal corotation resonances are also observed. For known co-orbitals, Saturn is the current source of the main destabilizing force but this is not enough to eject a minor body from the 1:1 commensurability with Uranus. These objects must enter mean motion resonances with Jupiter and Neptune in order to be captured or become passing Centaurs. Asteroid 2010 EU65, a probable visitor from the Oort cloud, may have been stable for several Myr due to its comparatively low eccentricity. In addition, we show that the orbital properties and discovery circumstances of known objects can be used to outline a practical strategy by which additional Uranus' co-orbitals may be found.

In the last decade direct detection Dark Matter (DM) experiments have increased enormously their sensitivity and ton-scale setups have been proposed, especially using germanium and xenon targets with double readout and background discrimination capabilities. In light of this situation, we study the prospects for determining the parameters of Weakly Interacting Massive Particle (WIMP) DM (mass, spin-dependent (SD) and spin-independent (SI) cross section off nucleons) by combining the results of such experiments in the case of a hypothetical detection. In general, the degeneracy between the SD and SI components of the scattering cross section can only be removed using targets with different sensitivities to these components. Scintillating bolometers, with particle discrimination capability, very good energy resolution and threshold and a wide choice of target materials, are an excellent tool for a multitarget complementary DM search. We investigate how the simultaneous use of scintillating targets with different SD-SI sensitivities and/or light isotopes (as the case of CaF2 and NaI) significantly improves the determination of the WIMP parameters. In order to make the analysis more realistic we include the effect of uncertainties in the halo model and in the spin-dependent nuclear structure functions, as well as the effect of a thermal quenching different from 1.

Venus has three known co-orbitals: (322756) 2001 CK32, 2002 VE68 and 2012 XE133. The first two have absolute magnitudes 18 < H < 21. The third one, significantly smaller at H = 23.4 mag, is a recent discovery that signals the probable presence of many other similar objects: small transient companions to Venus that are also potentially hazardous asteroids (PHAs). Here, we study the dynamical evolution of the recently discovered asteroid 2013 ND15. At H = 24.1 mag, this minor body is yet another small Venus co-orbital and PHA, currently close to the Lagrangian point L4 and following the most eccentric path found so far for objects in this group. Due to its high eccentricity (0.6), its dynamics is different from that of the other three known Venus co-orbitals even if they all are near-Earth objects (NEOs). A Monte Carlo simulation that uses the orbital data and discovery circumstances of the four objects as proxies to estimate the current size of this population, indicates that the number of high-eccentricity, low-inclination Venus co-orbital NEOs may have been greatly underestimated by current models. Three out of four known objects were discovered with solar elongation at perigee greater than 135 degrees even if visibility estimates show that less than 4 per cent of these objects are expected to reach perigee at such large elongations. Nearly 70 per cent of the objects discussed here have elongation at perigee < 90 degrees and 65 per cent are prospective PHAs.

Historically, it has been assumed that globular and open clusters never interact. However, recent evidence suggests that: globular clusters passing through the disk may be able to perturb giant molecular clouds (GMCs) triggering formation of open clusters and some old open clusters may be linked to accreted globulars. Here, we further explore the existence of possible dynamical connections between globular and open clusters, and realize that the most obvious link must be in the form of gravitational interactions. If open clusters are born out of GMCs, they have to move in similar orbits. If we accept that globulars can interact with GMCs, triggering star formation, it follows that globular and open clusters must also interact. Consistently, theoretical arguments as well as observational evidence, show that globular and open clusters certainly are interacting populations and their interactions are far more common than usually thought, especially for objects part of the bulge/disk. Monte Carlo calculations confirm that conclusion. Globular clusters seem capable of not only inducing formation of open clusters but, more often, their demise. Relatively frequent high speed cluster encounters or cluster harassment may also cause, on the long-term, slow erosion and tidal truncation on the globulars involved. The disputed object FSR 1767 (2MASS-GC04) may be, statistically speaking, the best example of an ongoing interaction.

Bright fireballs or bolides are caused by meteoroids entering the Earth's atmosphere at high speed. On 2013 February 15, a superbolide was observed in the skies near Chelyabinsk, Russia. Such a meteor could be the result of the decay of an asteroid and here we explore this possibility applying a multistep approach. First, we use available data and Monte Carlo optimization (validated using 2008 TC3 as template) to obtain a robust solution for the pre-impact orbit of the Chelyabinsk impactor (semimajor axis = 1.62 au, eccentricity = 0.53, inclination = 3.82 deg, longitude of the ascending node = 326.41 deg and argument of perihelion = 109.44 deg). Then, we use this most probable orbit and numerical analysis to single out candidates for membership in, what we call, the Chelyabinsk asteroid family. Finally, we perform N-body simulations to either confirm or reject any dynamical connection between candidates and impactor. We find reliable statistical evidence on the existence of the Chelyabinsk cluster. It appears to include multiple small asteroids and two relatively large members: 2007 BD7 and 2011 EO40. The most probable parent body for the Chelyabinsk superbolide is 2011 EO40. The orbits of these objects are quite perturbed as they experience close encounters not only with the Earth-Moon system but also with Venus, Mars and Ceres. Under such conditions, the cluster cannot be older than about 20-40 kyr.

The object G2 will pass within nearly 100 au from Sgr A* in 2014. Due to its very short periapse, the study of the dynamical evolution of this object in the short-term future may offer some insight into the region surrounding the supermassive black hole at the centre of the Galaxy. With this scenario in mind, it has recently been proposed by Bartos et al. (arXiv:1302.3220) that, prior to its perinigricon, G2 will likely experience multiple encounters with members of the black hole and neutron-star populations believed to orbit near the Galactic Centre. Here, we further explore this possibility and study the general case for collisions with the G2 object using the latest orbital solutions provided by Phifer et al. (arXiv:1304.5280) and Gillessen et al., (arXiv:1306.1374) and a Monte Carlo approach to estimate the minimum orbit intersection distance (MOID) with G2 as a function of the orbital parameters of the incoming body. Our results indicate that encounters at distances closer than 100 au started to become statistically significant only during the last few years or so. MOIDs under 100 au are statistically more probable for the most dynamically cold orbits. If there is a population of objects moving in low-inclination, low-eccentricity orbits around the central black hole, the highest probability for a close encounter with G2 is found to be in the period 2014 January-March but enhanced activity due to encounters may start as early as 2013 July-August.

If stars form in clusters but most stars belong to the field, understanding the details of the transition from the former to the latter is imperative to explain the observational properties of the field. Aging open clusters are one of the sources of field stars. The disruption rate of open clusters slows down with age but, as an object gets older, the distinction between the remaining cluster or open cluster remnant (OCR) and the surrounding field becomes less and less obvious. As a result, finding good OCR candidates or confirming the OCR nature of some of the best candidates still remains elusive. One of these objects is NGC 1252, a scattered group of about 20 stars in Horologium. Here we use new wide-field photometry in the UBVI pass-bands, proper motions from the Yale/San Juan SPM 4.0 catalogue and high resolution spectroscopy concurrently with results from N-body simulations to decipher NGC 1252's enigmatic character. Spectroscopy shows that most of the brightest stars in the studied area are chemically, kinematically and spatially unrelated to each other. However, after analysing proper motions, we find one relevant kinematic group. This sparse object is relatively close (about 1 kpc), metal poor and is probably not only one of the oldest clusters (3 Gyr) within 1.5 kpc from the Sun but also one of the clusters located farthest from the disc, at an altitude of nearly -900 pc. That makes NGC 1252 the first open cluster that can be truly considered a high Galactic altitude OCR: an unusual object that may hint at a star formation event induced on a high Galactic altitude gas cloud. TW Hor and HD 20286 are unlikely to be part of NGC 1252. NGC 1252 17 is identified as an unrelated, Population II cannonball star moving at about 400 km/s.

Near-Earth objects (NEOs) moving in resonant, Earth-like orbits are potentially important. On the positive side, they are the ideal targets for robotic and human low-cost sample return missions and a much cheaper alternative to using the Moon as an astronomical observatory. On the negative side and even if small in size (2-50 m), they have an enhanced probability of colliding with the Earth causing local but still significant property damage and loss of life. Here, we show that the recently discovered asteroid 2013 BS45 is an Earth co-orbital, the sixth horseshoe librator to our planet. In contrast with other Earth's co-orbitals, its orbit is strikingly similar to that of the Earth yet at an absolute magnitude of 25.8, an artificial origin seems implausible. The study of the dynamics of 2013 BS45 coupled with the analysis of NEO data show that it is one of the largest and most stable members of a previously undiscussed dynamically cold group of small NEOs experiencing repeated trappings in the 1:1 commensurability with the Earth. This new resonant family is well constrained in orbital parameter space and it includes at least 10 other transient members: 2003 YN107, 2006 JY26, 2009 SH2 and 2012 FC71 among them. 2012 FC71 represents the best of both worlds as it is locked in a Kozai resonance and is unlikely to impact the Earth. Objects in this group could be responsible for the production of Earth's transient irregular natural satellites.

Apart from Mercury that has no known co-orbital companions, Venus remains as the inner planet that hosts the smallest number of known co-orbitals (two): (322756) 2001 CK32 and 2002 VE68. Both objects have absolute magnitudes 18 < H < 21 and were identified as Venus co-orbitals in 2004. Here, we analyse the orbit of the recently discovered asteroid 2012 XE133 with H = 23.5 mag to conclude that it is a new Venus co-orbital currently following a transitional trajectory between Venus' Lagrangian points L5 and L3. The object could have been a 1:1 librator for several thousand years and it may leave the resonance with Venus within the next few hundred years, after a close encounter with the Earth. Our calculations show that its dynamical status as co-orbital, as well as that of the two previously known Venus co-orbitals, is controlled by the Earth-Moon system with Mercury playing a secondary role. The three temporary co-orbitals follow rather chaotic but similar trajectories with e-folding times of the order of 100 yr. Out of the three co-orbitals, 2012 XE133 currently follows the most perturbed path. An actual collision with the Earth during the next 10000 yr cannot be discarded. Extrapolation of the number distribution of Venus co-orbitals as a function of the absolute magnitude suggests that dozens of objects similar to 2012 XE133 could be transient companions to Venus. Some additional objects that were or will be transient co-orbitals to Venus are also briefly discussed.

Mars was second to Jupiter in being recognized as the host of a population of Trojan minor bodies. Since 1990, five asteroids - 5261 Eureka, (101429) 1998 VF31, (121514) 1999 UJ7, 2001 DH47 and (311999) 2007 NS2 - have been identified as Mars Trojans, one L4 and four L5. Dynamical and spectroscopic evidence suggests that some Mars Trojans may be remnants of the original planetesimal population that formed in the terrestrial planets region. Here we revisit the long-term dynamical evolution of the previously known Mars Trojans and show that 2011 SC191, 2011 SL25 and 2011 UN63 are also trailing (L5) Mars Trojans. They appear to be as stable as Eureka and may have been Trojans over the age of the Solar system. The fact that five Trojans move in similar orbits and one of them is a binary may point to the disruption of a larger body early in the history of the Solar system. Such a catastrophic event may also explain the apparently strong asymmetry in terms of number of objects (one versus seven) between the L4 and L5 regions. Future spectroscopic observations should be able to reject or confirm a putative common chemical signature that may lend further support to collisional scenario.

Stable co-orbital motion with Uranus is vulnerable to planetary migration but temporary co-orbitals may exist today. So far only two candidates have been suggested, both moving on horseshoe orbits: 83982 Crantor (2002 GO9) and 2000 SN331. (83982) Crantor is currently classified in the group of the Centaurs by the MPC although the value of its orbital period is close to that of Uranus. Here we revisit the topic of the possible 1:1 commensurability of (83982) Crantor with Uranus and also explore its dynamical past and look into its medium-term stability and future orbital evolution. (83982) Crantor currently moves inside Uranus' co-orbital region on a complex horseshoe orbit. The motion of this object is primarily driven by the influence of the Sun and Uranus, although Saturn plays a significant role in destabilizing its orbit. The precession of the nodes of (83982) Crantor, which is accelerated by Saturn, controls its evolution and short-term stability. Although this object follows a temporary horseshoe orbit, more stable trajectories are possible and we present 2010 EU65 as a long-term horseshoe librator candidate in urgent need of follow-up observations. Available data indicate that the candidate 2000 SN331 is not a Uranus' co-orbital.

Numerical simulations suggest that Neptune primordial co-orbitals may outnumber the equivalent population hosted by Jupiter, yet the objects remain elusive. Since the first discovery in 2001 just 10 minor planets have been identified as Neptune co-orbitals. On the other hand, some simulations predict that a negligible fraction of passing bodies are captured into the 1:1 commensurability with Neptune today. Hundreds of objects have been discovered in the outer solar system during the various wide-field surveys carried out during the past decade, and many of them have been classified using cuts in the pericentre and other orbital elements. This leads to possible misclassifications of resonant objects. Here, we explore this possibility to uncover neglected Neptune co-orbitals. We confirm that 4 objects previously classified as Centaurs by the MPC currently are temporary Neptune co-orbitals. (148975) 2001 XA255 is the most dynamically unstable of the four. It appears to be a relatively recent (50 kyr) visitor from the scattered disk in its way to the inner solar system. (310071) 2010 KR59 is following a complex horseshoe orbit, (316179) 2010 EN65 is in the process of switching from leading to trailing Trojan, and 2012 GX17 is a promising trailing Trojan candidate in urgent need of follow-up. The four objects move in highly inclined orbits and have significant eccentricities. These dynamically hot objects are not primordial 1:1 librators, but are captured and likely originated from beyond Neptune.

In the solar system, quasi-satellites move in a 1:1 mean motion resonance going around their host body like a retrograde satellite but their mutual separation is well beyond the Hill radius and the trajectory is not closed as they orbit the Sun not the host body. So far, minor bodies temporarily trapped in the quasi-satellite dynamical state have been identified around Venus, Earth, the dwarf planet (1) Ceres, the large asteroid (4) Vesta, Jupiter, Saturn and Neptune. Using computer simulations, Tiscareno and Malhotra have predicted the existence of a small but significant population of minor bodies moving in a 1:1 mean motion resonance with Pluto. Here we show using N-body calculations that the Plutino 15810 (1994 JR1) is currently an accidental quasi-satellite of Pluto and it will remain as such for nearly 350,000 yr. By accidental we mean that the quasi-satellite phase is triggered (or terminated) not by a direct gravitational influence in the form of a discrete close encounter but as a result of a resonance. The relative mean longitude of the Plutino circulates with a superimposed libration resulting from the oscillation of the orbital period induced by the 2:3 mean motion resonance with Neptune. These quasi-satellite episodes are recurrent with a periodicity of nearly 2 Myr. This makes 15810 the first minor body moving in a 1:1 mean motion resonance with Pluto. It also makes Pluto the second dwarf planet, besides Ceres, to host a quasi-satellite. Our finding confirms that the quasi-satellite resonant phase is not restricted to small bodies orbiting major planets but is possible for dwarf planets/asteroids too. Moreover, 15810 could be considered as a possible secondary target for NASA's Pluto-Kuiper Belt Mission New Horizons after the main Pluto flyby in 2015.

Upon discovery, asteroid (309239) 2007 RW10 was considered a Neptune Trojan candidate. The object is currently listed by the Minor Planet Center as a Centaur but it is classified as a Scattered Disk or Trans-Neptunian Object by others. Now that its arc-length is 8,154 days and has been observed for more than 20 years, a more robust classification should be possible. Here we explore the orbital behaviour of this object in order to reveal its current dynamical status; we perform N-body simulations in both directions of time to investigate the evolution of its orbital elements. In particular, we study the librational properties of the mean longitude that currently librates around the value of the mean longitude of Neptune with an amplitude of nearly 50 degrees and a period of about 7,500 years. Our calculations show that it has been in its present dynamical state for about 12,500 years and it will stay there for another 12,500 years. Therefore, its current state is relatively short-lived. Due to its chaotic behaviour, the object may have remained in the 1:1 mean motion resonance with Neptune for several 100,000 years at most, undergoing transitions between the various resonant states. (309239) 2007 RW10 is currently a quasi-satellite, the first object of this dynamical class to be discovered around Neptune. With a diameter of about 250 km, it is the largest known co-orbital in the Solar System. Its significant eccentricity (0.30) and orbital inclination (36 degrees), strongly suggest that it did not form in situ but was captured, likely from beyond Neptune. With an apparent magnitude of 21.1 at opposition (October), it is well suited for spectroscopic observations that may provide information on its surface composition and hence eventually its origin.

Minor planet 2002 VE68 was identified as a quasi-satellite of Venus shortly after its discovery. At that time its data-arc span was only 24 days, now it is 2,947 days. Here we revisit the topic of the dynamical status of this remarkable object as well as look into its dynamical past and explore its future orbital evolution which is driven by close encounters with both the Earth-Moon system and Mercury. In our calculations we use a Hermite integration scheme, the most updated ephemerides and include the perturbations by the eight major planets, the Moon and the three largest asteroids. We confirm that 2002 VE68 currently is a quasi-satellite of Venus and it has remained as such for at least 7,000 yr after a close fly-by with the Earth. Prior to that encounter the object may have already been co-orbital with Venus or moving in a classical, non-resonant Near-Earth Object (NEO) orbit. The object drifted into the quasi-satellite phase from an L4 Trojan state. We also confirm that, at aphelion, dangerously close encounters with the Earth (under 0.002 AU, well inside the Hill sphere) are possible. We find that 2002 VE68 will remain as a quasi-satellite of Venus for about 500 yr more and its dynamical evolution is controlled not only by the Earth, with a non-negligible contribution from the Moon, but by Mercury as well. 2002 VE68 exhibits resonant (or near resonant) behavior with Mercury, Venus and the Earth. Our calculations indicate that an actual collision with the Earth during the next 10,000 yr is highly unlikely but encounters as close as 0.04 AU occur with a periodicity of 8 years.

The basic properties of the candidate binary cluster population in the Magellanic Clouds and Galaxy are similar. The fraction of candidate binary systems is $\sim$10% and the pair separation histogram exhibits a bimodal distribution commonly attributed to their transient nature. However, if primordial pairs cannot survive for long as recognizable bound systems, how are they ending up? Here, we use simulations to confirm that merging, extreme tidal distortion and ionization are possible depending on the initial orbital elements and mass ratio of the pair. The nature of the dominant evolutionary path largely depends on the strength of the local tidal field. Merging is observed for initially close primordial binary clusters but also for wider pairs in nearly parabolic orbits. Its characteristic timescale depends on the initial orbital semi-major axis, eccentricity, and cluster pair mass ratio, becoming shorter for closer, more eccentric equal mass pairs. Shredding or extreme tidal distortion of the less massive cluster and subsequent separation is observed in all pairs with appreciably different masses. Wide pairs steadily evolve into the separated twins state characterized by the presence of tidal bridges and separations of 200-500 pc after one Galactic orbit. In the Galaxy, the vast majority of observed binary candidates appear to be following this evolutionary path which translates into the dominant peak (25-30 pc) in the pair separation distribution. The secondary peak at smaller separations (10-15 pc) can be explained as due to close pairs in almost circular orbits and/or undergoing merging. Merged clusters exhibit both peculiar radial density and velocity dispersion profiles shaped by synchronization and gravogyro instabilities. Both simulations and observations show that, for the range of parameters studied here, long term binary cluster stability in the Galactic disk is unlikely.

Shortly after birth, open clusters start dissolving; gradually losing stars into the surrounding star field. The time scale for complete disintegration depends both on their initial membership and location within the Galaxy. Open clusters undergoing the terminal phase of cluster disruption or open cluster remnants (OCRs) are notoriously difficult to identify. From an observational point, a combination of low number statistics and minimal contrast against the general stellar field conspire to turn them into very challenging objects. To make the situation even worst, random samples of field stars often display features that may induce to classify them erroneously as extremely evolved open clusters. In this paper, we provide a detailed study of the stellar content and kinematics of NGC 6863, a compact group of a few stars located in Aquila and described by the POSS as a non existent cluster. Nonetheless, this object has been recently classified as OCR. The aim of the present work is to either confirm or disprove its OCR status by a detailed star-by-star analysis. The analysis is performed using wide-field photometry in the UBVI pass-band, proper motions from the UCAC3 catalogue, and high resolution spectroscopy as well as results from extensive $N$-body calculations. Our results show that the four brightest stars commonly associated to NGC 6863 form an asterism, a group of non-physically associated stars projected together, leading to the conclusion that NGC 6863 is not a real open cluster.

Hierarchical star formation leads to a progressive decrease in the clustering of star clusters both in terms of spatial scale and age. Consistently, the statistical analysis of positions and ages of clusters in the Milky Way disk strongly suggests that a correlation between the duration of star formation in a region and its size does exist. The average age difference between pairs of open clusters increases with their separation as the ~0.16 power. In contrast and for the Large Magellanic Cloud, Efremov & Elmegreen (1998) found that the age difference scales with the ~0.35 power of the region size. This discrepancy may be tentatively interpreted as an argument in support of intrinsically shorter (faster) star formation time-scales in smaller galaxies. However, if both the effects of cluster dissolution and incompleteness are taken into consideration, the average age difference between cluster pairs in the Galaxy increases with their separation as the ~0.4 power. This result implies that the characteristic time-scale for coherent, clustered-mode star formation is nearly 1 Myr. Therefore, the overall consequence of ignoring the effect of cluster dissolution is to overestimate the star formation time-scale. On the other hand, in the Galactic disk and for young clusters separated by less than three times the characteristic cluster tidal radius (10 pc), the average age difference is 16 Myr, which suggests common origin. (Abridged)

Observations indicate that the fraction of potential binary star clusters in the Magellanic Clouds is about 10%. In contrast, it is widely accepted that the binary cluster frequency in the Galaxy disk is much lower. Here we investigate the multiplicity of clusters in the Milky Way disk to either confirm or disprove this dearth of binaries. We quantify the open cluster multiplicity using complete, volume-limited samples from WEBDA and NCOVOCC. At the Solar Circle, at least 12% of all open clusters appear to be experiencing some type of interaction with another cluster; i.e., are possible binaries. As in the Magellanic Clouds, the pair separation histogram hints of a bimodal distribution. Nearly 40% of identified pairs are probably primordial. Most of the remaining pairs could be undergoing some type of close encounter, perhaps as a result of orbital resonances. Confirming early theoretical predictions, the characteristic time scale for destruction of bound pairs in the disk is 200 Myr or one galactic orbit. Our results show that the fraction of possible binary clusters in the Galactic disk is comparable to that in the Magellanic Clouds.

Observations indicate that present-day star formation in the Milky Way disk takes place in stellar ensembles or clusters rather than in isolation. Bound, long lived stellar groups are known as open clusters. They gradually lose stars and in their final evolutionary stages they are severely disrupted leaving an open cluster remnant made of a few stars. In this paper, we study in detail the stellar content and kinematics of the poorly populated star cluster NGC1901. This object appears projected against the Large Magellanic Cloud. The aim of the present work is to derive the current evolutionary status, binary fraction, age and mass of this stellar group. These are fundamental quantities to compare with those from N-body models in order to study the most general topic of star cluster evolution and dissolution.The analysis is performed using wide-field photometry in the UBVI pass-band, proper motions from the UCAC.2 catalog, and 3 epochs of high resolution spectroscopy, as well as results from extensive N-body calculations.The star group NGC1901 is found to be an ensemble of solar metallicity stars, 400+/-100 Myr old, with a core radius of 0.23 pc, a tidal radius of 1.0 pc, and located at 400+/-50 pc from the Sun. Out of 13 confirmed members, only 5 single stars have been found. Its estimated present-day binary fraction is at least 62%. The calculated heliocentric space motion of the cluster is not compatible with possible membership in the Hyades stream.Our results show that NGC1901 is a clear prototype of open cluster remnant characterized by a large value of the binary fraction and a significant depletion of low-mass stars. In the light of numerical simulations, this is compatible with NGC1901 being what remains of a larger system initially made of 500-750 stars.

Shaviv (2003) has shown evidence for a correlation between variations in the Galactic cosmic ray flux reaching Earth and the glaciation period record on Earth during the last 2 Gyr. If the flux of cosmic rays is mainly the result of Type II supernovae, an additional correlation between the star formation history of the Solar Neighbourhood and the timing of past ice ages is expected. Higher star formation rate implies increased cosmic ray flux and this may translate into colder climate through a rise in the average low altitude cloud cover. Here we reanalyze the correlation between this star formation history and the glaciation period record on Earth using a volume limited open cluster sample. Numerical modeling and recent observational data indicate that the correlation is rather strong but only if open clusters within 1.5 kpc from the Sun are considered.

We have derived the star formation history of the Milky Way disk over the last 2 Gyr from the age distribution diagram of a large sample of open clusters comprising more than 580 objects. By interpreting the age distribution diagram using numerical results from an extensive library of N-body calculations carried out during the last ten years, we reconstruct the recent star formation history of the Milky Way disk. Our analysis suggests that superimposed on a relatively small level of constant star formation activity mainly in small-N star clusters, the star formation rate has experienced at least 5 episodes of enhanced star formation lasting about 0.2 Gyr with production of larger clusters. This cyclic behavior seems to show a period of 0.4+/-0.1 Gyr.

We use detailed numerical simulations and theoretical estimates to show that, if confirmed, the unusually brief microlensing events observed by Sahu et al. (2001) in the field of the globular cluster M22 might be explained as a result of microlensing by a population of clustered MACHOs, a dark cluster or RAMBO, not associated with the globular cluster. If real, this dark cluster would be located between M22 and the Galactic bulge and could include at least $10^6$ substellar members with a typical size of 1-3 pc. Bound planets in wide or/and eccentric orbits are also able to reproduce the observed microlensing behaviour, but only if multiplanet systems (including large Kuiper-belt-like objects) are abundant, although, our calculations argue against the latter scenario as the ionization rate in M22 is very high. Dynamically ejected or lone planets are, in principle, incompatible with the observational findings as they either escape their parent cluster in a relatively short time-scale after ejection or segregate toward the outskirts of the cluster. We discuss additional implications of the dark cluster scenario, including the existence of a population of RAMBOs toward the Galactic bulge.