We report on the observation and measurement of astrometry, photometry, morphology, and activity of the interstellar object 3I/ATLAS, also designated C/2025 N1 (ATLAS), with the NSF-DOE Vera C. Rubin Observatory. The third interstellar object, comet 3I/ATLAS, was first discovered on UT 2025 July 1. Serendipitously, the Rubin Observatory collected imaging in the area of the sky inhabited by the object during regular commissioning activities. We successfully recovered object detections from Rubin visits spanning UT 2025 June 21 (10 days before discovery) to UT 2025 July 7. Facilitated by Rubin's high resolution and large aperture, we report on the detection of cometary activity as early as June 21st, and observe it throughout. We measure the location and magnitude of the object on 37 Rubin images in r, i, and z bands, with typical precision of about 20 mas (100 mas, systematic) and about 10 mmag, respectively. We use these to derive improved orbit solutions, and to show there is no detectable photometric variability on hourly timescales. We derive a V-band absolute magnitude of H_V = (13.7 +/- 0.2) mag, and an equivalent effective nucleus radius of around (5.6 +/- 0.7) km. These data represent the earliest observations of this object by a large (8-meter class) telescope reported to date, and illustrate the type of measurements (and discoveries) Rubin's Legacy Survey of Space and Time (LSST) will begin to provide once operational later this year.

We present here in full detail the evolution of the angular momentum deficit (AMD) during collisions as it was described in (Laskar, PRL,2000). Since then, the AMD has been revealed to be a key parameter for the understanding of the outcome of planetary formation models. We define here the AMD-stability criterion that can be easily verified on a newly discovered planetary system. We show how AMD-stability can be used to establish a classification of the multiplanet systems in order to exhibit the planetary systems that are long-term stable because they are AMD-stable, and those that are AMD-unstable which then require some additional dynamical studies to conclude on their stability. The AMD-stability classification is applied to the 131 multiplanet systems from The Extrasolar Planet Encyclopaedia database (this http URL) for which the orbital elements are sufficiently well known.

We present orbital fits and dynamical masses for HIP 113201AB and HIP 36985AB, two M1 + mid-M dwarf binary systems monitored as part of the SPHERE SHINE survey. To robustly determine ages via gyrochronology, we undertook a photometric monitoring campaign for HIP 113201 and for GJ 282AB, the two wide K star companions to HIP 36985, using the 40 cm Remote Observatory Atacama Desert (ROAD) telescope. We adopt ages of 1.2$\pm$0.1 Gyr for HIP 113201AB and 750$\pm$100 Myr for HIP 36985AB. To derive dynamical masses for all components of these systems, we used parallel-tempering Markov Chain Monte Carlo sampling to fit a combination of radial velocity, direct imaging, and Gaia and Hipparcos astrometry. Fitting the direct imaging and radial velocity data for HIP 113201 yields a primary mass of 0.54$\pm$0.03 M$_{\odot}$, fully consistent with its M1 spectral type, and a secondary mass of 0.145$\pm$ M$_{\odot}$. The secondary masses derived with and without including Hipparcos/Gaia data are more massive than the 0.1 M$_{\odot}$ estimated mass from the photometry of the companion. An undetected brown dwarf companion to HIP 113201B could be a natural explanation for this apparent discrepancy. At an age >1 Gyr, a 30 M$_{Jup}$ companion to HIP 113201B would make a negligible (<1$\%$) contribution to the system luminosity, but could have strong dynamical impacts. Fitting the direct imaging, radial velocity, and Hipparcos/Gaia proper motion anomaly for HIP 36985AB, we find a primary mass of 0.54$\pm$0.01 M$_{\odot}$ and a secondary mass of 0.185$\pm$0.001 M$_{\odot}$ which agree well with photometric estimates of component masses, the masses estimated from $M_{K}$-- mass relationships for M dwarf stars, and previous dynamical masses in the literature.

We present an example of low-carbon research activity carried out by astrophysicists and focused on ecology and environmental protection, with direct impacts on territories and society. This project serves as an illustration of action research for an astrophysics lab in the context of the current ecological crisis.

We use the tidal deformations of the Moon induced by the Earth and the Sun as a tool for studying the inner structure of our satellite. Based on measurements of the degree-two tidal Love numbers k2 and h2 and dissipation coefficients from the GRAIL mission, Lunar Laser Ranging and Laser Altimetry on board of the LRO spacecraft, we perform Monte Carlo samplings for 120,000 possible combinations of thicknesses and viscosities for two classes of the lunar models. The first one includes a uniform core, a low viscosity zone (LVZ) at the core-mantle boundary, a mantle and a crust. The second one has an additional inner core. All models are consistent with the lunar total mass as well as its moment of inertia. By comparing predicted and observed parameters for the tidal deformations we find that the existence of an inner core cannot be ruled out. Furthermore, by deducing temperature profiles for the LVZ and an Earth-like mantle, we obtain stringent constraints on the radius (500 +- 1) km, viscosity,21 (4.5 +- 0.8) x 10^16 Pa.s and the density (3400 +- 10) kg/m^3 of the LVZ. We also infer the first estimation for the outer core viscosity, (2.07 +- 1.03) x 10^17 Pa.s, for two different possible structures: a Moon with a 70 km thick outer core and a large inner core (290 km radius with a density of 6000 kg/m3), and a Moon with a thicker outer core (169 km thick) but a denser and smaller inner core (219 km radius for 8000 kg/m^3).

We report measured spectral properties for more than 1000 NEOs, representing>5% of the currently discovered population. Thermal flux detected below 2.5 {\mu}m allows us to make albedo estimates for nearly 50 objects, including two comets. Additional spectral data are reported for more than 350 Mars-crossing asteroids. Most of these measurements were achieved through a collaboration between researchers at the Massachusetts Institute of Technology and the University of Hawaii, with full cooperation of the NASA Infrared Telescope Facility (IRTF) on Mauna Kea. We call this project the MIT-Hawaii Near-Earth Object Spectroscopic Survey (MITHNEOS; myth-neos).

Determining habitable zones in binary star systems can be a challenging task due to the combination of perturbed planetary orbits and varying stellar irradiation conditions. The concept of "dynamically informed habitable zones" allows us, nevertheless, to make predictions on where to look for habitable worlds in such complex environments. Dynamically informed habitable zones have been used in the past to investigate the habitability of circumstellar planets in binary systems and Earth-like analogs in systems with giant planets. Here, we extend the concept to potentially habitable worlds on circumbinary orbits. We show that habitable zone borders can be found analytically even when another giant planet is present in the system. By applying this methodology to Kepler-16, Kepler-34, Kepler-35, Kepler-38, Kepler-64, Kepler-413, Kepler-453, Kepler-1647 and Kepler-1661 we demonstrate that the presence of the known giant planets in the majority of those systems does not preclude the existence of potentially habitable worlds. Among the investigated systems Kepler-35, Kepler-38 and Kepler-64 currently seem to offer the most benign environment. In contrast, Kepler-16 and Kepler-1647 are unlikely to host habitable worlds.

Earth-mass planets are expected to have atmospheres and experience thermal tides raised by the host star. These tides transfer energy to the planet that can counter the dissipation from bodily tides. Indeed, even a relatively thin atmosphere can drive the rotation of these planets away from the synchronous state. Here we revisit the dynamical evolution of planets undergoing thermal atmospheric tides. We use a novel approach based on a vectorial formalism, which is frame independent and valid for any configuration of the system, including any eccentricity and obliquity values. We provide the secular equations of motion after averaging over the mean anomaly and the argument of the pericenter, which are suitable to model the long-term spin and orbital evolution of the planet.

In this paper, we give a new proof of the classical KAM theorem on the persistence of an invariant quasi-periodic torus, whose frequency vector satisfies the Bruno-R\"ussmann condition, in real-analytic non-degenerate Hamiltonian systems close to integrable. The proof, which uses rational approximations instead of small divisors estimates, is an adaptation to the Hamiltonian setting of the method we introduced in a previous work for perturbations of constant vector fields on the torus.

Earth-like planets in the habitable zone of low-mass stars undergo strong tidal effects that modify their spin states. These planets are expected to host dense atmospheres that can also play an important role in the spin evolution. On one hand, gravitational tides tend to synchronise the rotation with the orbital mean motion, but on the other hand, thermal atmospheric tides push the rotation away and may lead to asynchronous equilibria. Here, we investigate the complete tidal evolution of Earth-like planets by taking into account the effect of obliquity and eccentric orbits. We adopted an Andrade rheology for the gravitational tides and benchmarked the unknown parameters with the present rotation of Venus. We then applied our model to Earth-like planets, and we show that asynchronous rotation can be expected for planets orbiting stars with masses between 0.4 and 0.9 $M_\odot$ and semi-major axes between 0.2 and 0.7 au. Interestingly, we find that Earth-like planets in the habitable zone of stars with masses $\sim 0.8$ $M_\odot$ may end up with an equilibrium rotation of 24 hours. We additionally find that these planets can also develop high obliquities, which may help sustain temperate environments.

Due to tidal interactions in the Earth-Moon system, the spin of the Earth slows down and the Moon drifts away. This recession of the Moon is now measured with great precision, but it has been realized, more than fifty years ago, that simple tidal models extrapolated back in time lead to an age of the Moon that is by far incompatible with the geochronological and geochemical evidence. In order to evade this problem, more elaborate models have been proposed, taking into account the oceanic tidal dissipation. However, these models did not fit both the estimated lunar age and the present rate of lunar recession simultaneously. Here we present a physical model that reconciles these two constraints and yields a unique solution of the tidal history. This solution fits well the available geological proxies for the history of the Earth-Moon system and consolidates the cyclostratigraphic method. The resulting evolution involves multiple crossings of resonances in the oceanic dissipation that are associated with significant and rapid variations in the lunar orbital distance, the Earth's length of the day, and the Earth's obliquity.

Lambert's theorem (1761) on the elapsed time along a Keplerian arc drew the attention of several prestigious mathematicians. In particular, they tried to give simple and transparent proofs of it (see our timeline \S 9). We give two new proofs. The first one (\S 4) goes along the lines of Hamilton's variational proof in his famous paper of 1834, but we shorten his computation in such a way that the hypothesis is now used without redundancy. The second (\S 6) is among the few which are close to Lambert's geometrical proof. It starts with the new remark that two Keplerian arcs related by the hypothesis of Lambert's theorem correspond to each other through an affine map. We also show (\S 7) that despite the singularities due to the occurrence of collisions, the classes of arcs related by Lambert's theorem all have the same topology. We give (\S 8) some simple related results about conic sections and affine transformations.

The population of large asteroids is thought to be primordial and they are the most direct witnesses of the early history of our Solar System. Those satellites allow study of the mass, and hence density and internal structure. We study here the properties of the triple asteroid (107) Camilla from lightcurves, stellar occultations, optical spectroscopy, and high-contrast and high-angular-resolution images and spectro-images. Using 80 positions over 15 years, we determine the orbit of its larger satellite to be circular, equatorial, and prograde, with RMS residuals of 7.8 mas. From 11 positions in three epochs only, in 2015 and 2016, we determine a preliminary orbit for the second satellite. We find the orbit to be somewhat eccentric and slightly inclined to the primary's equatorial plane, reminiscent of the inner satellites of other asteroid triple systems. Comparison of the near-infrared spectrum of the larger satellite reveals no significant difference with Camilla. Hence, these properties argue for a formation of the satellites by excavation from impact and re-accumulation of ejecta. We determine the spin and 3-D shape of Camilla. The model fits well each data set. We determine Camilla to be larger than reported from modeling of mid-infrared photometry, with a spherical-volume-equivalent diameter of 254 $\pm$ 36 km (3 $\sigma$ uncertainty), in agreement with recent results from shape modeling (Hanus2017+). Combining the mass of (1.12 $\pm$ 0.01) $\times$ 10$^{19}$ kg determined from the dynamics of the satellites and the volume from the 3-D shape model, we determine a density of 1,280 $\pm$ 130 SI. From this density, and considering Camilla's spectral similarities with (24) Themis and (65) Cybele (for which water ice coating on surface grains was reported), we infer a silicate-to-ice mass ratio of 1-6, with a 10-30% macroporosity.

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

Extrasolar planetary systems commonly exhibit planets on eccentric orbits, with many systems located near or within mean-motion resonances, showcasing a wide diversity of orbital architectures. Such complex systems challenge traditional secular theories, which are limited to first-order approximations in planetary masses or rely on expansions in orbital elements--eccentricities, inclinations, and semi-major axis ratios--that are subject to convergence issues, especially in highly eccentric, inclined, or tightly-packed systems. To overcome these limitations, we develop a numerical approach to second-order perturbation theory based on the Lie transform formalism. Our method avoids the need for expansions in orbital elements, ensuring broader applicability and more robust convergence. We first outline the Hamiltonian framework for the 3-body planetary problem, and apply a canonical transformation to eliminate fast angle dependencies, deriving the secular Hamiltonian up to second order in the mass ratio. We then use the fast Fourier transform algorithm to numerically simulate, in an accurate way, the long-term evolution of planetary systems near or away from mean-motion resonances. Finally, we validate our methods against well-known planetary configurations, such as the Sun-Jupiter-Saturn system, as well as to exoplanetary systems like WASP-148, TIC 279401253 and GJ 876, demonstrating the applicability of our models across a wide range of planetary configurations.

We report on a pencil-beam survey of the Taurid Swarm, a possible concentration of bodies in the Taurid meteoroid stream associated with the 7:2 mean-motion resonance with Jupiter. Canada-France-Hawaii Telescope MegaCam observations reaching apparent magnitudes of 24.5 in the gri filter were taken over three nights. Rates of motion on the sky allowed for the quick elimination of main-belt objects from the over 1000 moving sources seen. Eight candidates with on-sky rates of motion consistent with Taurids were detected, but seven were subsequently shown to be non-Taurids (Hungarias, Mars-crossers, etc). One object might be a 60 m class Taurid but not enough data was collected and its orbit remains ambiguous. Our results are consistent with no Taurid Swarm members observed, and an upper limit of fewer than 3e3 - 3e4 objects down to H=25.6 +/- 0.3 (diameter of 34-76 m assuming a 2P/Encke-like albedo) at the 95% confidence level. While meteor observations have confirmed the Taurid Swarm's existence at meter and smaller sizes, our results indicate that the current mass budget of the swarm at 100 m sizes does not require an outsize parent to explain it.

Planetary systems orbiting close binaries are valuable testing grounds for planet formation and migration models. More detections with good mass measurements are needed. We present a new planet discovered during the BEBOP survey for circumbinary exoplanets using radial velocities. We use data taken with the SOPHIE spectrograph at the Observatoire de Haute-Provence, and perform a spectroscopic analysis to obtain high precision radial velocities. This planet is the first radial velocity detection of a previously unknown circumbinary system. The planet has a mass of $0.56$ $M_{Jup}$ and orbits its host binary in 550 days with an eccentricity of 0.25. Compared to most of the previously known circumbinary planets, BEBOP-3b has a long period (relative to the binary) and a high eccentricity. There also is a candidate outer planet with a $\sim1400$ day orbital period. We test the stability of potential further candidate signals inside the orbit of BEBOP-3b, and demonstrate that there are stable orbital solutions for planets near the instability region which is where the Kepler circumbinary planets are located. We also use our data to obtain independent dynamical masses for the two stellar components of the eclipsing binary using High Resolution Cross-Correlation Spectroscopy (HRCCS), and compare those results to a more traditional approach, finding them compatible with one another.

The Lunar Laser Ranging (LLR) experiment has accumulated 50 years of range data of improving accuracy from ground stations to the laser retroreflector arrays (LRAs) on the lunar surface. The upcoming decade offers several opportunities to break new ground in data precision through the deployment of the next generation of single corner-cube lunar retroreflectors and active laser transponders. This is likely to expand the LLR station network. Lunar dynamical models and analysis tools have the potential to improve and fully exploit the long temporal baseline and precision allowed by millimetric LLR data. Some of the model limitations are outlined for future efforts. Differential observation techniques will help mitigate some of the primary limiting factors and reach unprecedented accuracy. Such observations and techniques may enable the detection of several subtle signatures required to understand the dynamics of the Earth-Moon system and the deep lunar interior. LLR model improvements would impact multi-disciplinary fields that include lunar and planetary science, Earth science, fundamental physics, celestial mechanics and ephemerides.

Examples are given of solutions of the planar N-body problem which remain the same for at least two systems of masses with the same sum and same center of mass. The least value of N achieved up to now with this property is 474, a number which had been announced in the first author's thesis.

The planetary ephemerides play a crucial role for spacecraft navigation, mission planning, reduction and analysis of the most precise astronomical observations. The construction of such ephemerides is highly constrained by the tracking observations, in particular range, of the space probes collected by the tracking stations on the Earth. The present planetary ephemerides (DE, INPOP, EPM) are mainly based on such observations. However, the data used by the planetary ephemerides are not the direct raw tracking data, but measurements deduced after the analysis of raw data made by the space agencies and the access to such processed measurements remains difficult in terms of availability. The first part of the study deals with the analysis of the Mars Global Surveyor (MGS) tracking data. The tracking observations containing one-, two-, and three-way Doppler and two-way range are then used to reconstruct MGS orbit precisely. As a supplementary exploitation of MGS, we derived the solar corona model and estimated the average electron density along the line of sight separately for slow and fast wind regions. Fitting the planetary ephemerides, including additional data which were corrected for the solar corona perturbations, noticeably improves the extrapolation capability of the planetary ephemerides and the estimation of the asteroid masses (Verma et al., 2013). The second part of the thesis deals with the complete analysis of the MESSENGER tracking data. This analysis improved the Mercury ephemeris up to two order of magnitude compared to any latest ephemerides. Such high precision ephemeris, INPOP13a, is then used to perform general relativity tests of PPN-formalism. Our estimations of PPN parameters (? and ) are the most stringent than previous results (Verma et al., 2014).