Context. We report the exploitation of a sample of Solar System observations based on data from the third Gaia Data Release Gaia DR3 of nearly 157 000 asteroids. It extends the epoch astrometric solution over the time coverage planned for the Gaia DR4 which is not expected before the end of 2025. This data set covers more than one full orbital period for the vast majority of these asteroids. The orbital solutions are derived from the Gaia data alone over a relatively short arc compared to the observation history of many of these asteroids. BR  Aims The work aims to produce orbital elements for a large set of asteroids based on 66 months of accurate astrometry provided by Gaia and to assess the accuracy of these orbital solutions with a comparison to the best available orbits derived from independent observations. A second validation is performed with accurate occultation timings. BR  Methods We processed the raw astrometric measurements of Gaia to obtain astrometric positions of moving objects with 1D submas accuracy at the bright end. For each asteroid that we matched to the data an orbit fitting was attempted in the form of the best fit of the initial conditions at the median epoch. The force model included Newtonian and relativistic accelerations to derive the observation equations which were solved with a linear leastsquares fit. BR  Results Orbits are provided in the form of state vectors in the International Celestial Reference Frame for 156 764 asteroids including nearEarth objects mainbelt asteroids and Trojans. For the asteroids with the best observations the formal relative uncertainty SUBaSUBa is better than 10SUP10SUP. Results are compared to orbits available from the Jet Propulsion Laboratory and MPC. Their orbits are based on much longer data arcs but from positions of lower quality. The relative differences in semimajor axes have a mean of 5  10SUP10SUP and a scatter of 5  10SUP9SUP.
['astrometry', 'celestial mechanics', 'minor planets', 'asteroids: general', 'Astrophysics - Earth and Planetary Astrophysics', 'Astrophysics - Instrumentation and Methods for Astrophysics']
Gaia Focused Product Release Asteroid orbital solution. Properties and assessment
{'Properties and assessment': 'Gaia Collaboration, P. David 1 , F. Mignard 2 , D. Hestro ff er 1 , P. Tanga 2 , F. Spoto 3 , J. Berthier 1 , T. Pauwels 4 , W. Roux 5 , A. Barbier 5 , A. Cellino 6 , B. Carry 2 , M. Delbo 2 , A. Dell\'Oro 7 , C. Fouron 8 , L. Galluccio 2 , S.A. Klioner 9 , N. Mary 10 , K. Muinonen 11 , 12 , C. Ordenovic 2 , I. Oreshina-Slezak 2 , C. Panem 5 , J.-M. Petit 13 , J. Portell 14 , 15 , 16 , A.G.A. Brown 17 , W. Thuillot 1 , A. Vallenari 18 , T. Prusti 19 , J.H.J. de Bruijne 19 , F. Arenou 20 , C. Babusiaux 21 , M. Biermann 22 , O.L. Creevey 2 , C. Ducourant 23 , D.W. Evans 24 , L. Eyer 25 , R. Guerra 26 , A. Hutton 27 , C. Jordi 14 , 15 , 16 , U. Lammers 26 , L. Lindegren 28 , X. Luri 14 , 15 , 16 , S. Randich 7 , P. Sartoretti 20 , R. Smiljanic 29 , N.A. Walton 24 , C.A.L. Bailer-Jones 30 , U. Bastian 22 , M. Cropper 31 , R. Drimmel 6 , D. Katz 20 , C. Soubiran 23 , F. van Leeuwen 24 , M. Audard 25 , 32 , J. Bakker 26 , R. Blomme 4 , J. Castañeda 33 , 14 , 16 , F. De Angeli 24 , C. Fabricius 16 , 14 , 15 , M. Fouesneau 30 , Y. Frémat 4 , A. Guerrier 5 , E. Masana 16 , 14 , 15 , R. Messineo 34 , C. Nicolas 5 , K. Nienartowicz 35 , 32 , F. Pailler 5 , P. Panuzzo 20 , F. Riclet 5 , G.M. Seabroke 31 , R. Sordo 18 , F. Thévenin 2 , G. Gracia-Abril 36 , 22 , D. Teyssier 37 , M. Altmann 22 , 38 , K. Benson 31 , P.W. Burgess 24 , D. Busonero 6 , G. Busso 24 , H. Cánovas 37 , N. Cheek 39 , G. Clementini 40 , Y. Damerdji 41 , 42 , M. Davidson 43 , P. de Teodoro 26 , L. Delchambre 41 , E. Fraile Garcia 44 , D. Garabato 45 , P. García-Lario 26 , N. Garralda Torres 46 , P. Gavras 44 , R. Haigron 20 , N.C. Hambly 43 , D.L. Harrison 24 , 47 , D. Hatzidimitriou 48 , J. Hernández 26 , S.T. Hodgkin 24 , B. Holl 25 , 32 , S. Jamal 30 , S. Jordan 22 , A. Krone-Martins 49 , 50 , A.C. Lanzafame 51 , 52 , W. Lö ffl er 22 , A. Lorca 27 , O. Marchal 53 , P.M. Marrese 54 , 55 , A. Moitinho 50 , M. Nuñez Campos 27 , P. Osborne 24 , E. Pancino 7 , 55 , A. Recio-Blanco 2 , M. Riello 24 , L. Rimoldini 32 , A.C. Robin 13 , T. Roegiers 56 , L.M. Sarro 57 , M. Schultheis 2 , C. Siopis 58 , M. Smith 31 , A. Sozzetti 6 , E. Utrilla 27 , M. van Leeuwen 24 , K. Weingrill 59 , U. Abbas 6 , P. Ábrahám 60 , 61 , A. Abreu Aramburu 46 , C. Aerts 62 , 63 , 30 , G. Altavilla 54 , 55 , M.A. Álvarez 45 , J. Alves 64 , R.I. Anderson 65 , T. Antoja 14 , 15 , 16 , D. Baines 66 , S.G. Baker 31 , Z. Balog 22 , 30 , C. Barache 38 , D. Barbato 25 , 6 , M. Barros 67 , M.A. Barstow 68 , S. Bartolomé 16 , 14 , 15 , D. Bashi 69 , 70 , N. Bauchet 20 , N. Baudeau 8 , U. Becciani 51 , L.R. Bedin 18 , I. Bellas-Velidis 71 , M. Bellazzini 40 , W. Beordo 6 , 72 , A. Berihuete 73 , M. Bernet 14 , 15 , 16 , C. Bertolotto 34 , S. Bertone 6 , L. Bianchi 74 , A. Binnenfeld 75 , A. Blazere 76 , T. Boch 53 , A. Bombrun 77 , S. Bouquillon 38 , 78 , A. Bragaglia 40 , J. Braine 23 , L. Bramante 34 , E. Breedt 24 , A. Bressan 79 , N. Brouillet 23 , E. Brugaletta 51 , B. Bucciarelli 6 , 72 , A.G. Butkevich 6 , R. Buzzi 6 , E. Ca ff au 20 , R. Cancelliere 80 , S. Cannizzo 10 , R. Carballo 81 , T. Carlucci 38 , M.I. Carnerero 6 , J.M. Carrasco 16 , 14 , 15 , J. Carretero 82 , 83 , S. Carton 10 , L. Casamiquela 23 , 20 , M. Castellani 54 , A. Castro-Ginard 17 , V. Cesare 51 , P. Charlot 23 , L. Chemin 84 , V. Chiaramida 34 , A. Chiavassa 2 , N. Chornay 24 , 32 , R. Collins 43 , G. Contursi 2 , W.J. Cooper 85 , 6 , T. Cornez 10 , M. Crosta 6 , 86 , C. Crowley 77 , C. Dafonte 45 , P. de Laverny 2 , F. De Luise 87 , R. De March 34 , R. de Souza 88 , A. de Torres 77 , E.F. del Peloso 22 , A. Delgado 44 , T.E. Dharmawardena 30 , S. Diakite 89 , C. Diener 24 , E. Distefano 51 , C. Dolding 31 , K. Dsilva 58 , J. Durán 44 , H. Enke 59 , P. Esquej 44 , C. Fabre 76 , M. Fabrizio 54 , 55 , S. Faigler 69 , M. Fatovi\'c 90 , G. Fedorets 11 , 91 , J. Fernández-Hernández 44 , P. Fernique 53 , F. Figueras 14 , 15 , 16 , Y. Fournier 59 , M. Gai 6 , M. Galinier 2 , A. Garcia-Gutierrez 16 , 14 , 15 , M. García-Torres 92 , A. Garofalo 40 , E. Gerlach 9 , R. Geyer 9 , P. Giacobbe 6 , G. Gilmore 24 , 93 , S. Girona 94 , G. Giu ff rida 54 , R. Gomel 69 , A. Gomez 45 , J. González-Núñez 95 , I. González-Santamaría 45 , E. Gosset 41 , 96 , M. Granvik 11 , 97 , V. Gregori Barrera 16 , 14 , 15 , R. Gutiérrez-Sánchez 37 , M. Haywood 20 , A. Helmer 10 , A. Helmi 98 , K. Henares 66 , S.L. Hidalgo 99 , 100 , T. Hilger 9 , D. Hobbs 28 , C. Hottier 20 , H.E. Huckle 31 , M. Jabło\'nska 101 , 103 , F. Jansen 104 , Ó. Jiménez-Arranz 14 , 15 , 16 , J. Juaristi Campillo 22 , S. Khanna 6 , 98 , G. Kordopatis 2 , Á Kóspál 60 , 30 , 61 , Z. Kostrzewa-Rutkowska 17 , M. Kun 60 , S. Lambert 38 , A.F. Lanza 51 , J.-F. Le Campion 23 , Y. Lebreton 105 , 106 , T. Lebzelter 64 , S. Leccia 107 , I. Lecoeur-Taibi 32 , G. Lecoutre 13 , S. Liao 108 , 6 , 109 , L. Liberato 2 , 110 , E. Licata 6 , H.E.P. Lindstrøm 6 , 111 , 112 , T.A. Lister 113 , E. Livanou 48 , A. Lobel 4 , C. Loup 53 , L. Mahy 4 , R.G. Mann 43 , M. Manteiga 114 , J.M. Marchant 115 , M. Marconi 107 , D. Marín Pina 14 , 15 , 16 , S. Marinoni 54 , 55 , D.J. Marshall 116 , J. Martín Lozano 39 , J.M. Martín-Fleitas 27 , G. Marton 60 , A. Masip 16 , 14 , 15 , D. Massari 40 , A. Mastrobuono-Battisti 20 , T. Mazeh 69 , P.J. McMillan 28 , J. Meichsner 9 , S. Messina 51 , D. Michalik 19 , N.R. Millar 24 , A. Mints 59 , D. Molina 15 , 14 , 16 , R. Molinaro  \n107 , L. Molnár 60 , 117 , 61 , G. Monari 53 , M. Monguió 14 , 15 , 16 , P. Montegri ff o 40 , A. Montero 39 , R. Mor 118 , 15 , 16 , A. Mora 27 , R. Morbidelli 6 , T. Morel 41 , D. Morris 43 , N. Mowlavi 25 , D. Munoz 10 , T. Muraveva 40 , C.P. Murphy 26 , I. Musella 107 , Z. Nagy 60 , S. Nieto 44 , L. Noval 10 , A. Ogden 24 , C. Pagani 119 , I. Pagano 51 , L. Palaversa 90 , P.A. Palicio 2 , L. Pallas-Quintela 45 , A. Panahi 69 , S. Payne-Wardenaar 22 , L. Pegoraro 5 , A. Penttilä 11 , P. Pesciullesi 44 , A.M. Piersimoni 87 , M. Pinamonti 6 , F.-X. Pineau 53 , E. Plachy 60 , 117 , 61 , G. Plum 20 , E. Poggio 2 , 6 , D. Pourbaix † 58 , 96 , A. Prša 120 , L. Pulone 54 , E. Racero 39 , 121 , M. Rainer 7 , 122 , C.M. Raiteri 6 , P. Ramos 123 , 14 , 16 , M. Ramos-Lerate 37 , M. Ratajczak 101 , P. Re Fiorentin 6 , S. Regibo 62 , C. Reylé 13 , V. Ripepi 107 , A. Riva 6 , H.-W. Rix 30 , G. Rixon 24 , N. Robichon 20 , C. Robin 10 , M. Romero-Gómez 14 , 15 , 16 , N. Rowell 43 , F. Royer 20 , D. Ruz Mieres 24 , K.A. Rybicki 124 , G. Sadowski 58 , A. Sáez Núñez 16 , 14 , 15 , A. Sagristà Sellés 22 , J. Sahlmann 44 , V. Sanchez Gimenez 16 , 14 , 15 , N. Sanna 7 , R. Santoveña 45 , M. Sarasso 6 , C. Sarrate Riera 33 , 14 , 16 , E. Sciacca 51 , J.C. Segovia 39 , D. Ségransan 25 , S. Shahaf 124 , A. Siebert 53 , 125 , L. Siltala 11 , E. Slezak 2 , R.L. Smart 6 , 85 , O.N. Snaith 20 , 126 , E. Solano 127 , F. Solitro 34 , D. Souami 105 , 128 , J. Souchay 38 , L. Spina 18 , E. Spitoni 2 , 129 , L.A. Squillante 34 , I.A. Steele 115 , H. Steidelmüller 9 , J. Surdej 41 , L. Szabados 60 , F. Taris 38 , M.B. Taylor 130 , R. Teixeira 88 , K. Tisani\'c 90 , L. Tolomei 34 , F. Torra 33 , 14 , 16 , G. Torralba Elipe 45 , 131 , 132 , M. Trabucchi 133 , 25 , M. Tsantaki 7 , A. Ulla 134 , 135 , N. Unger 25 , O. Vanel 20 , A. Vecchiato 6 , D. Vicente 94 , S. Voutsinas 43 , M. Weiler 16 , 14 , 15 , Ł. Wyrzykowski 101 , H. Zhao 2 , 136 , J. Zorec 137 , T. Zwitter 138 , L. Balaguer-Núñez 16 , 14 , 15 , N. Leclerc 20 , S. Morgenthaler 139 , G. Robert 10 , and S. Zucker 75  \n- 1 IMCCE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Lille, 77 av. Denfert-Rochereau, 75014 Paris, France\n- 2 Université Côte d\'Azur, Observatoire de la Côte d\'Azur, CNRS, Laboratoire Lagrange, Bd de l\'Observatoire, CS 34229, 06304 Nice Cedex 4, France\n- 3 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 15, Cambridge, MA 02138, USA\n- 4 Royal Observatory of Belgium, Ringlaan 3, 1180 Brussels, Belgium\n- 5 CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France\n- 6 INAF - Osservatorio Astrofisico di Torino, via Osservatorio 20, 10025 Pino Torinese (TO), Italy\n- 7 INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy\n- 8 Telespazio for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France\n- 9 Lohrmann Observatory, Technische Universität Dresden, Mommsenstraße 13, 01062 Dresden, Germany\n- 10 Thales Services for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France\n- 11 Department of Physics, University of Helsinki, P.O. Box 64, 00014 Helsinki, Finland\n- 12 Finnish Geospatial Research Institute FGI, Vuorimiehentie 5, 02150 Espoo, Finland\n- 13 Institut UTINAM CNRS UMR6213, Université de Franche-Comté, OSU THETA Franche-Comté Bourgogne, Observatoire de Besançon, BP1615, 25010 Besançon Cedex, France\n- 14 Institut de Ciències del Cosmos (ICCUB), Universitat de Barcelona (UB), Martí i Franquès 1, 08028 Barcelona, Spain\n- 15 Departament de Física Quàntica i Astrofísica (FQA), Universitat de Barcelona (UB), c. Martí i Franquès 1, 08028 Barcelona, Spain\n- 16 Institut d\'Estudis Espacials de Catalunya (IEEC), c. Gran Capità, 2-4, 08034 Barcelona, Spain\n- 17 Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands\n- 18 INAF - Osservatorio astronomico di Padova, Vicolo Osservatorio 5, 35122 Padova, Italy\n- 19 European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands  \n20  \nGEPI, Observatoire de Paris, Université PSL, CNRS, 5 Place Jules Janssen, 92190 Meudon, France  \n- 21 Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France\n- 22 Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität Heidelberg, Mönchhofstr. 12-14, 69120 Heidelberg, Germany\n- 23 Laboratoire d\'astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geo ff roy Saint-Hilaire, 33615 Pessac, France\n- 24 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, United Kingdom\n- 25 Department of Astronomy, University of Geneva, Chemin Pegasi 51, 1290 Versoix, Switzerland\n- 26 European Space Agency (ESA), European Space Astronomy Centre (ESAC), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 27 Aurora Technology for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 28 Lund Observatory, Division of Astrophysics, Department of Physics, Lund University, Box 43, 22100 Lund, Sweden\n- 29 Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, ul. Bartycka 18, 00-716 Warsaw, Poland\n- 30 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany\n- 31 Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking, Surrey RH5 6NT, United Kingdom\n- 32 Department of Astronomy, University of Geneva, Chemin d\'Ecogia 16, 1290 Versoix, Switzerland\n- 33 DAPCOM Data Services, c. dels Vilabella, 5-7, 80500 Vic, Barcelona, Spain\n- 34 ALTEC S.p.a, Corso Marche, 79,10146 Torino, Italy\n- 35 Sednai Sàrl, Geneva, Switzerland\n- 36 Gaia DPAC Project O ffi ce, ESAC, Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain  \n- 37 Telespazio UK S.L. for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 38 SYRTE, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, LNE, 61 avenue de l\'Observatoire 75014 Paris, France\n- 39 Serco Gestión de Negocios for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 40 INAF - Osservatorio di Astrofisica e Scienza dello Spazio di Bologna, via Piero Gobetti 93 / 3, 40129 Bologna, Italy\n- 41 Institut d\'Astrophysique et de Géophysique, Université de Liège, 19c, Allée du 6 Août, B-4000 Liège, Belgium\n- 42 CRAAG - Centre de Recherche en Astronomie, Astrophysique et Géophysique, Route de l\'Observatoire Bp 63 Bouzareah 16340 Algiers, Algeria  \n43  \nInstitute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, United Kingdom  \n- 44 RHEA for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 45 CIGUS CITIC - Department of Computer Science and Information Technologies, University of A Coruña, Campus de Elviña s / n, A Coruña, 15071, Spain\n- 46 ATG Europe for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 47 Kavli Institute for Cosmology Cambridge, Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA\n- 48 Department of Astrophysics, Astronomy and Mechanics, National and Kapodistrian University of Athens, Panepistimiopolis, Zografos, 15783 Athens, Greece\n- 49 Donald Bren School of Information and Computer Sciences, University of California, Irvine, CA 92697, USA\n- 50 CENTRA, Faculdade de Ciências, Universidade de Lisboa, Edif. C8, Campo Grande, 1749-016 Lisboa, Portugal\n- 51 INAF - Osservatorio Astrofisico di Catania, via S. Sofia 78, 95123 Catania, Italy\n- 52 Dipartimento di Fisica e Astronomia ""Ettore Majorana"", Università di Catania, Via S. Sofia 64, 95123 Catania, Italy\n- 53 Université de Strasbourg, CNRS, Observatoire astronomique de Strasbourg, UMR 7550, 11 rue de l\'Université, 67000 Strasbourg, France\n- 54 INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 00078 Monte Porzio Catone (Roma), Italy\n- 55 Space Science Data Center - ASI, Via del Politecnico SNC, 00133 Roma, Italy\n- 56 HE Space Operations BV for European Space Agency (ESA), Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands\n- 57\n- Dpto. de Inteligencia Artificial, UNED, c / Juan del Rosal 16, 28040 Madrid, Spain\n- 58 Institut d\'Astronomie et d\'Astrophysique, Université Libre de Bruxelles CP 226, Boulevard du Triomphe, 1050 Brussels, Belgium\n- 59 Leibniz Institute for Astrophysics Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany\n- 60 Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Eötvös Loránd Research Network (ELKH), MTA Centre of Excellence, Konkoly Thege Miklós út 15-17, 1121 Budapest, Hungary\n- 61 ELTE Eötvös Loránd University, Institute of Physics, 1117, Pázmány Péter sétány 1A, Budapest, Hungary\n- 62 Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, 3001 Leuven, Belgium  \n63  \nDepartment of Astrophysics  \n/  \nIMAPP, Radboud University, P.O.Box 9010, 6500 GL Nijmegen, The Netherlands  \n- 64 University of Vienna, Department of Astrophysics, Türkenschanzstraße 17, A1180 Vienna, Austria\n- 65 Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny, 1290 Versoix, Switzerland\n- 66 Quasar Science Resources for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 67 LASIGE, Faculdade de Ciências, Universidade de Lisboa, Edif. C6, Campo Grande, 1749-016 Lisboa, Portugal\n- 68 School of Physics and Astronomy , University of Leicester, University Road, Leicester LE1 7RH, United Kingdom\n- 69 School of Physics and Astronomy, Tel Aviv University, Tel Aviv 6997801, Israel\n- 70 Cavendish Laboratory, JJ Thomson Avenue, Cambridge CB3 0HE, United Kingdom\n- 71 National Observatory of Athens, I. Metaxa and Vas. Pavlou, Palaia Penteli, 15236 Athens, Greece\n- 72 University of Turin, Department of Physics, Via Pietro Giuria 1, 10125 Torino, Italy\n- 73 Depto. Estadística e Investigación Operativa. Universidad de Cádiz, Avda. República Saharaui s / n, 11510 Puerto Real, Cádiz, Spain\n- 74 EURIX S.r.l., Corso Vittorio Emanuele II 61, 10128, Torino, Italy\n- 75 Porter School of the Environment and Earth Sciences, Tel Aviv University, Tel Aviv 6997801, Israel\n- 76 ATOS for CNES Centre Spatial de Toulouse, 18 avenue Edouard Belin, 31401 Toulouse Cedex 9, France\n- 77 HESpace Operations BV for European Space Agency (ESA), Camino bajo del Castillo, s / n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain\n- 78 LFCA / DAS,Universidad de Chile,CNRS,Casilla 36-D, Santiago, Chile\n- 79 SISSA - Scuola Internazionale Superiore di Studi Avanzati, via Bonomea 265, 34136 Trieste, Italy\n- 80 University of Turin, Department of Computer Sciences, Corso Svizzera 185, 10149 Torino, Italy\n- 81 Dpto. de Matemática Aplicada y Ciencias de la Computación, Univ. de Cantabria, ETS Ingenieros de Caminos, Canales y Puertos, Avda. de los Castros s / n, 39005 Santander, Spain\n- 82 Institut de Física d\'Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona), Spain\n- 83 Port d\'Informació Científica (PIC), Campus UAB, C. Albareda s / n, 08193 Bellaterra (Barcelona), Spain\n- 84 Instituto de Astrofísica, Universidad Andres Bello, Fernandez Concha 700, Las Condes, Santiago RM, Chile\n- 85 Centre for Astrophysics Research, University of Hertfordshire, College Lane, AL10 9AB, Hatfield, United Kingdom\n- 86 University of Turin, Mathematical Department ""G.Peano"", Via Carlo Alberto 10, 10123 Torino, Italy\n- 87 INAF - Osservatorio Astronomico d\'Abruzzo, Via Mentore Maggini, 64100 Teramo, Italy\n- 88 Instituto de Astronomia, Geofìsica e Ciências Atmosféricas, Universidade de São Paulo, Rua do Matão, 1226, Cidade Universitaria, 05508-900 São Paulo, SP, Brazil\n- 89 Mésocentre de calcul de Franche-Comté, Université de Franche-Comté, 16 route de Gray, 25030 Besançon Cedex, France\n- 90 Ru¯der Boškovi\'c Institute, Bijeniˇcka cesta 54, 10000 Zagreb, Croatia  \n- 91 Astrophysics Research Centre, School of Mathematics and Physics, Queen\'s University Belfast, Belfast BT7 1NN, UK\n- 92 Data Science and Big Data Lab, Pablo de Olavide University, 41013, Seville, Spain\n- 93 Institute of Astrophysics, FORTH, Crete, Greece\n- 94 Barcelona Supercomputing Center (BSC), Plaça Eusebi Güell 1-3, 08034-Barcelona, Spain\n- 95 ETSE Telecomunicación, Universidade de Vigo, Campus Lagoas-Marcosende, 36310 Vigo, Galicia, Spain\n- 96 F.R.S.-FNRS, Rue d\'Egmont 5, 1000 Brussels, Belgium\n- 97 Asteroid Engineering Laboratory, Luleå University of Technology, Box 848, S-981 28 Kiruna, Sweden\n- 98 Kapteyn Astronomical Institute, University of Groningen, Landleven 12, 9747 AD Groningen, The Netherlands\n- 99 IAC - Instituto de Astrofisica de Canarias, Via Láctea s / n, 38200 La Laguna S.C., Tenerife, Spain\n- 100 Department of Astrophysics, University of La Laguna, Via Láctea s / n, 38200 La Laguna S.C., Tenerife, Spain\n- 101 Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warszawa, Poland\n- 102 Research School of Astronomy\n- 103 Astrophysics, Australian National University, Cotter Road, Weston, ACT 2611, Australia\n- 104 European Space Agency (ESA, retired), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201AZ, Noordwijk, The Netherlands\n- 105 LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université de Paris, 5 Place Jules Janssen, 92190 Meudon, France\n- 106 Université Rennes, CNRS, IPR (Institut de Physique de Rennes) - UMR 6251, 35000 Rennes, France\n- 107 INAF - Osservatorio Astronomico di Capodimonte, Via Moiariello 16, 80131, Napoli, Italy\n- 108 Shanghai Astronomical Observatory, Chinese Academy of Sciences, 80 Nandan Road, Shanghai 200030, People\'s Republic of China\n- 109 University of Chinese Academy of Sciences, No.19(A) Yuquan Road, Shijingshan District, Beijing 100049, People\'s Republic of China\n- 110 São Paulo State University, Grupo de Dinâmica Orbital e Planetologia, CEP 12516-410, Guaratinguetá, SP, Brazil\n- 111 Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark\n- 112 DXC Technology, Retortvej 8, 2500 Valby, Denmark\n- 113 Las Cumbres Observatory, 6740 Cortona Drive Suite 102, Goleta, CA 93117, USA\n- 114 CIGUS CITIC, Department of Nautical Sciences and Marine Engineering, University of A Coruña, Paseo de Ronda 51, 15071, A Coruña, Spain\n- 115\n- Astrophysics Research Institute, Liverpool John Moores University, 146 Brownlow Hill, Liverpool L3 5RF, United Kingdom\n- 116 IRAP, Université de Toulouse, CNRS, UPS, CNES, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France\n- 117 MTA CSFK Lendület Near-Field Cosmology Research Group, Konkoly Observatory, MTA Research Centre for Astronomy and Earth Sciences, Konkoly Thege Miklós út 15-17, 1121 Budapest, Hungary\n- 118 Pervasive Technologies s.l., c. Saragossa 118, 08006 Barcelona, Spain\n- 119 School of Physics and Astronomy, University of Leicester, University Road, Leicester LE1 7RH, United Kingdom\n- 120 Villanova University, Department of Astrophysics and Planetary Science, 800 E Lancaster Avenue, Villanova PA 19085, USA\n- 121 Departmento de Física de la Tierra y Astrofísica, Universidad Complutense de Madrid, 28040 Madrid, Spain\n- 122 INAF - Osservatorio Astronomico di Brera, via E. Bianchi, 46, 23807 Merate (LC), Italy\n- 123 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan\n- 124 Department of Particle Physics and Astrophysics, Weizmann Institute of Science, Rehovot 7610001, Israel\n- 125 Centre de Données Astronomique de Strasbourg, Strasbourg, France\n- 126 University of Exeter, School of Physics and Astronomy, Stocker Road, Exeter, EX2 7SJ, United Kingdom\n- 127 Departamento de Astrofísica, Centro de Astrobiología (CSIC-INTA), ESA-ESAC. Camino Bajo del Castillo s / n. 28692 Villanueva de la Cañada, Madrid, Spain\n- 128 naXys, Department of Mathematics, University of Namur, Rue de Bruxelles 61, 5000 Namur, Belgium\n- 129 INAF. Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, 34131, Trieste, Italy\n- 130 H H Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, United Kingdom\n- 131 Escuela de Arquitectura y Politécnica - Universidad Europea de Valencia, Spain\n- 132 Escuela Superior de Ingeniería y Tecnología - Universidad Internacional de la Rioja, Spain\n- 133 Department of Physics and Astronomy G. Galilei, University of Padova, Vicolo dell\'Osservatorio 3, 35122, Padova, Italy\n- 134 Applied Physics Department, Universidade de Vigo, 36310 Vigo, Spain\n- 135 Instituto de F\'isica e Ciencias Aeroespaciais (IFCAE), Universidade de Vigo, Á Campus de As Lagoas, 32004 Ourense, Spain\n- 136 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210023, China\n- 137 Sorbonne Université, CNRS, UMR7095, Institut d\'Astrophysique de Paris, 98bis bd. Arago, 75014 Paris, France\n- 138 Faculty of Mathematics and Physics, University of Ljubljana, Jadranska ulica 19, 1000 Ljubljana, Slovenia\n- 139 Institute of Mathematics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Switzerland  \nReceived / Accepted', 'ABSTRACT': 'Context. We report the exploitation of a sample of Solar System observations based on data from the third Gaia Data Release (Gaia DR3) of nearly 157 000 asteroids that extends the epoch astrometric solution over the time coverage planned for the complete Gaia DR4, which is not expected before the end of 2025. This extended data set covers more than one full orbital period for the vast majority of these asteroids. The orbital solutions are derived from the Gaia data alone over a relatively short arc compared to the observation history of many of these asteroids.  \nAims. The work aims to produce orbital elements for a large set of asteroids based on 66 months of accurate astrometry provided by Gaia and to assess the accuracy of these orbital solutions with a comparison to the best available orbits estimated from independent observations. A second validation is performed with accurate occultation timings.  \nMethods. We processed the raw astrometric measurements of Gaia to obtain astrometric positions of moving objects with 1D sub-mas accuracy at the bright end. For each asteroid that we matched to the data, an orbit fitting was attempted in the form of the best fit of the initial conditions (state vector) at the median epoch. The force model included Newtonian and relativistic accelerations to derive the observation equations, which were solved with a linear least-squares fit.  \nResults. Orbits are provided in the form of state vectors in the International Celestial Reference Frame for 156 764 asteroids, including near-Earth objects, main-belt asteroids, and Trojans. For the asteroids with the best observations, the (formal) relative accuracy measured by σ a / a is better than 10 -10 . Results are compared to orbits available from the Jet Propulsion Laboratory and MPC. Their orbits are based on much longer data arcs, but the positions are of lower quality. The relative di ff erences in semi-major axes have a mean of 5 × 10 -10 and a scatter of 5 × 10 -9 .  \nKey words. Solar-System:minor-planets-orbits-astrometry; Gaiamission', '1. Introduction': 'In the course of its systematic scan of the sky, the ESA Gaia spacecraft has detected a wide variety of celestial sources. It will continue operating until the end of its operational life sometime in early 2025. The celestial sources are primarily stars from the Milky Way, which form the bulk of the massive Gaia catalogues that were released between 2016 and 2022 in data releases Gaia DR1, Gaia DR2, Gaia E-DR3, and Gaia DR3 proper (Gaia Collaboration et al. (2016a), Gaia Collaboration et al. (2016b), Gaia Collaboration et al. (2023)). However, sources closer to Earth are also caught in the Gaia net, with representatives of every category of small Solar System bodies (SSSBs), such as eear-Earth objects (NEOs), main-belt asteroids (MBA), Trojans, and a few more distant objects from the Kuiper belt and from the transNeptunian region. In addition, Gaia data include the largest planetary moons and a sample of periodic comets that were observed not too far from their perihelion passage.  \nIn this Gaia-coordinated Focus Product Release (Gaia FPR), the astrometry of 156,792 minor planets and 31 natural satellites is provided. However, this paper focuses only on the orbit computation of minor planets, and only this category of Solar System bodies is investigated. Likewise, photometric and spectral data are also beyond the scope of this work and have been largely presented with the Gaia DR3 results in Tanga et al. (2023).  \nGaia DR2 has shown the quality of Gaia observations and their huge value for asteroids. Based on the high-precision astrometry and because the time span covered by the observations was long enough (at least ≃ 1000 days), it was possible to derive accurate orbits for the observed asteroids, as was shown in Gaia DR3 (Tanga et al. 2023), who reported a summary of the observational technique and the epoch astrometric accuracy. More precisely, Gaia DR3 was built upon 34 months of Gaia operations extending from July 2014 to May 2017. A well-observed Solar System body came into the Gaia fields of view about 2530 times on average, resulting in up to nine accurate astrometric positions at every passage. 1 Tanga et al. (2023) ranked the orbits from excellent to very approximate, with a clear di ff erence in reliability as a function of the fraction of the orbital period covered over 34 months.  \nFor the majority of the MBAs, this time was too short to remove strong constraints between orbital elements. In the processing reported in this paper, which yielded the Focused Product Release, the data arc can be as long as 66 months. This typi-  \ncally is a little longer than the orbital period of an average asteroid that orbits in the main belt between Mars and Jupiter. Therefore, the situation considerably improved as soon as we had a few pairs of data points that lie one orbital period apart from each other. This is su ffi cient to remove the partial degeneracy that is observed when the arc was too short. Although we have almost twice as many observations as in Gaia DR3, the main source of improvement remains the complete coverage of the orbit in the data set, which is in sharp contrast with the previous solution in Gaia DR3, which was based on only 34 months of data. The length of arc that is used is the primary factor in the improved uncertainty of the orbit.  \nIn this paper, we aim to present and discuss the data that are released in the Gaia FPR, together with the results of the orbital fitting obtained in an extended time span, using a refined version of the fitting algorithm, with better weighting and an improved outlier rejection. Several shortcomings in the code were found and corrected in this new exploitation, and the orbits we provide now are the best we can achieve at the moment with the extremely high astrometric precision of Gaia. It is important to mention at this stage that the orbital solutions are derived from the Gaia data alone, meaning that the arc is still relatively short in comparison with the fits performed at MPC or JPL using all the available observations that so conveniently archived by the MPC, which were collected over decades or even centuries. As shown by the results, the relatively short time-span is somewhat balanced by the homogeneity of the data and the unparalleled single-observation accuracy and precision.  \nAs mentioned above, no new photometric data in G band or reflectance spectra are released in the FPR. To associate photometry with FPR data, the G fluxes and magnitudes published in DR3 can be adopted for the corresponding transits in common. New spectrophotometric processing will be available in DR4, which is currently scheduled for the last term of 2025. Likewise, new astrometry is published in the FPR for planetary satellites, but their orbits have not yet been adjusted on these new positions.  \nThe paper is broken down into a few sections. In Section 2 we discuss the observation material used in the paper, and then we very briefly recall the relevant properties of the astrometric solution in Section 3. In Section 4 we present the principle of the orbit fitting that produced the set of state vectors at a median epoch for each asteroid. The validation and discussion of the results are taken up in detail in Section 5, along with some recommendation for users.', '2. Properties of the data set': 'Gaia scans the sky throughout the year in a continuous way, without virtually any significant interruption beyond unavoidable but short dead times. The principles have been described in several papers from the Gaia Collaboration, and we refer to the  \nFig. 1. Distribution of the number of observations used in the orbital solution. Blue shows the first 100 000 numbered asteroids with a median of 330 individual observations (91 785 orbits), and red shows the numbered asteroids with numbers ≥ 100 000 and a median of 190 observations (64 979 orbits).  \n<!-- image -->  \nreliable information therein (Gaia Collaboration et al. 2016b). The relevant facts for the orbital solutions we use here are the repeated observations of the same source, with an interval of 6 to 10 weeks between two sequences. A typical sequence consists of a small set of successive passages in the Preceding Field of View (pfov) and Following Field of View ( ff ov). The most common sequence is just a simple sequence of two passages pfov to ff ov (PF) or ff ov to pfov (FP) within 6 hours, and then a few weeks before a new similar sequence. However, every asteroid during the lifetime of Gaia enjoyed one or two much longer sequences, such as PFPF... or FPFP... and so on. In some cases, there were more than 15 consecutive transits within 4 or 5 days in a row. For the well-observed asteroids, typically, 40 to 90 passages (field crossings) were made in 66 months, and the number of visibility periods (closely packed repeating sequences at distinct epochs) is 30 ± 10 on average. This has the strongest e ff ect for the orbit determination because it drives the sampling in orbital phases.  \nDuring each passage, a maximum of nine astrometric positions is obtained. These make up the observations at the CCD level. The number of observations that were used (after rejections in the iterative model fitting) is shown in Figure 1 for the asteroids that were observed best with numbers < 100 000 , and for the fainter and less frequently observed asteroids with numbers ≥ 100 000, typically with half as many transits ending with a successful Gaia detection. The transits were missed because the magnitude was too faint in some part of the orbits or because problems were found in the astrometric solution, which led to a rejection. The transition at 100 , 000 between the two groups is just a convenient number to obtain balanced distributions. It has no deeper significance.  \nBeyond the number of transits, a key parameter for the orbital solution is the time coverage, or in this context, more precisely, the fraction ( < 1 or > 1) of the orbit that is covered between the first and last observation. That is to say, the arc length measured in orbital periods. This is shown in blue in Figure 2 for the asteroids that were observed best with numbers < 100 , 000 , and the fainter and smaller asteroids with numbers ≥ 100 , 000 are superimposed in red. The typical coverage is about 1 . 2 orbital periods for the first group, and it is about 0 . 5 orbital periods for the Trojans. The typical coverage is shorter at about one orbital period for the fainter asteroids, with a marked tail of short arcs (less than one orbital period) that is caused not by the periods themselves, but by a systematic of the Gaia sensitivity thresh-  \nFig. 2. Observation coverage expressed in orbital periods. Asteroid numbers < 100 000 are plotted in blue, and the complementary set is shown in red.  \n<!-- image -->  \nd: The faint asteroids are not observed at all elongations that are theoretically allowed by the Gaia scan, but only when their distance to the Earth is small enough to have an apparent magnitude < 20 . 7 in the Gaia white band. This on-board brightness cut is imposed by the instrument sensitivity and the ground-link data rate. An orbit is better determined when the observation arc covers a complete orbit or more, and the solution is less constrained for short arcs. This qualitative insight is confirmed in Section 5, in which we analyse the solutions.', '3. Astrometric solution': "The astrometric processing of Solar System objects has been presented in Tanga et al. (2023, Sect. 3.4) and was described with more technical details in the Gaia DR3 documentation in Muinonen et al. (2022, Sect. 8.3.3). We refer to these papers for more detailed information. Here we just recall the key points needed for the readability of this paper and draw attention to the significant departures from the classical ground-based CCD or traditional photographic plate astrometry. As a mission destined to survey the sky several times, a scanning strategy was imposed early during the mission design. The Gaia time sampling is a direct consequence of the parameters selected to optimise this intricate scanning law. Because of the Gaia core science, the free parameters were not optimised with Solar System sources as the main targets, but for the best science return in stellar and fundamental physics. The Solar System objects are well sampled, however, just below the sky average, because their ecliptic latitude is predominantly low. This sky area is less frequently visited by Gaia than the mid latitudes.  \nEach passage of a source in the field of view nominally includes nine crossings of the astrometric CCDs either in the preceding or the following field of view. This provides as many 1D astrometric locations in the instrument frame (pixel number) that are perfectly time tagged. These repeated measurements form the basic astrometric data from Gaia. Based on the instrument calibration and the orientation parameters (attitude), these instrument-tied positions can be transformed into nine astrometric positions on the sky, with an excellent accuracy ( ≈ sub-milliarcsecon to milliarcsecond level) in the scan direction and much poorer accuracy (up to 0 . 5 '' ) in the transverse direction. When it is transformed into right ascension and declination, the uncertainty appears degraded in both coordinates, but the true statistical information is kept in the covariance matrix, where correlations are close to 1. This encodes the extreme accuracy in the along-scan direction. Therefore, a standard usage of the RA, DEC positions that does not take the correlation coef-  \nt into account would completely fail to bring the true Gaia accuracy.  \nThe positions used in this paper to fit the orbits are at the CCD level, meaning that the number of observations is always the number of validated CCD crossings. Each field transit has a maximum of nine crossings, but the actual average is between seven and eight crossings, so that there are about seven to eight times as many observations as field transits. These are also called observations in other contexts of the Gaia data processing. Consistency checks between the multiple 1D positions in the same transit were used to accept or reject observations for inclusion in the data release, taking the expected linear motion of the Solar System object during the ≈ 40 s of the field crossing into account. All these published observations are considered in the subsequent fits, and normal points for position and velocity were not computed.  \nWe implemented an updated version for the FPR reported in this paper for the astrometric calibration used in DR3. The global transformation that produces astrometric positions is also more accurate because it is now based on the full duration of the nominal Gaia mission of 66 months. Moreover, the FPR processing is mostly independent of the previous data releases. Even in the overlapping time between DR3 and the FPR, observations can be found in one data set but not in the other. This is a result of a di ff erent selection that was used by the internal validation procedures or residuals that are now below or above the acceptance threshold with the new orbit. These are rare occurrences, and they arise as a natural consequence of the reprocessing of the whole data set.  \nThe precision of the astrometric measurement in the alongscan direction has been much discussed in Tanga et al. (2023). The precision that is relevant here for the orbit fitting is the combination of a random and systematic part. Its variation with G magnitude is shown in Figure 3, where it is reproduced in the same way as in Sect. 3.4.2 of the Gaia DR3 paper The elementary measurement at the CCD level is better than 1 mas for G < 18 mag and grows to reach about 10 mas at G ≃ 20 mag. This makes Gaia data outstanding and unusual in comparison with the core of the historical observations: They are extremely accurate, rather dense, homogeneous, but over a limited time range that is not too short to preclude orbit fitting.  \nFig. 3. Error model in the along-scan direction for the SSO astrometry in Gaia DR3 as a function of the G magnitude. The total error is represented as given by the squared sum of the random and the systematic component. The colour represents the data density (yellow or light means a higher density). The thick line and the two thin lines on each side are the quantiles that correspond to the mean and the 1 σ level.  \n<!-- image -->  \nThe orbit can be fit either with observation equations expressed in the along-scan and across-scan directions, or in the more usual way in RA, DEC coordinates. The use of either coordinate set is rigorously equivalent theoretically in any case provided that the 2 × 2 variance or covariance matrix is well implemented in the global weight matrix. Tests have shown that identical results are obtained with both approaches, as expected, but great care must be exercised in using the statistical information in the covariance matrix in the right manner.", '4.1. Presentation': 'We only consider a subset of known asteroids here that have known (elliptical) orbits, as in Gaia DR3. We determined a leastsquares orbit solution for them that represents the best fit to Gaia data. Detection, follow-up, and initial orbit determination of newly discovered asteroids are made within another pipeline that was described in Carry et al. (2021).  \nThe overall principle of the orbit computation from the Gaia astrometry was summarised in Tanga et al. (2023). As a rule, the general procedure used in the Gaia DR3 release is also used here, but the code has been updated to correct minor bugs and improve the overall consistency between the fundamental constants in their form in astronomical or metric units.', '4.2. Units of time and length': 'Gaia time tagging is the Barycentric coordinate time (TCB) internally for the data collection and processing, and externally in the public releases. This means that, for example, the Solar System ephemerides have TCB as an independent variable and that all epoch data are given with a timing in TCB. For stellar astrometry, and even for epoch photometry, the di ff erence between TCB and TDB (or TT), which is about 20 s, is not an issue for the user. This is no longer true for Solar System objects moving at speed of several 10 mas per second, however. This has consequences for the orbit fitting. Observations are timed in TCB, and the equations of motions must be consistent with this timescale. More specifically, the solar mass parameter must be TCB compatible.  \nUnfortunately, a small glitch has remained in the software between the TCB and our using of an outdated Solar System GM ⊙ equal to the square of the Gauss constant k = 0 . 017 202 098 95. This is equivalent to our using a unit of length that is not strictly equal to the IAU astronomical unit. Because of the constraints of the complex integrated processing of the Gaia data and the huge validation pipeline, it was not possible to correct for the scale factor in the FPR results. It seems to be a benign correction, and it is, but nothing is simple in a complex system. Fundamentally, this is a mistake, and this issue has been fixed in the software for the planned Gaia DR4. The positive side of this unfortunate bug is that noticing a systematic e ff ect as small as 0 . 5 × 10 -8 in the results already means something good about the quality of the solution.  \nMixing the Gauss constant and the TCB timescale results in an auFPR = 149 597 871 473 . 216 m instead of the IAU fundamental constant for the scale length in SI equal to 149 597 870 700 . 0 m. The consequence is that the Gaia published state vectors of the asteroids must be scaled to scale them properly in TCB and au, as it should have been in the first place. Afterwards, the transformation to osculating elements provides  \na TCB-compatible semi-major axis. More explicitly, we have  \nρ = auFPR / au  \n= 149597871473 . 216 / 149597870700  \n- = 1 . 0000000051686297 ,  \nwhich with k 2 ≈ ( GM ⊙ )TDB , is essentially  \n( GM ⊙ )TCB ( GM ⊙ )TDB ! 1 / 3 .  \nNow the TCB compatible values in au and au / day can be computed from the published components X FPR , V FPR of the state vector with  \nX TCB = ρ × X FPR  \nV TCB = ρ × V FPR.  \nIt can also be added that when the FPR state vector in osculating elements is transformed before the scaling, only a single scale transformation has to be applied to the semi-major axis,  \na TCB = ρ × a FPR.  \nWe show in Section 5.3.1 how the Gaia TCB compatible orbits and the equivalent from JPL TDB compatible orbits are handled.', '4.3. Principles and algorithm': 'The Gaia-centric astrometric data in right ascension and declination at CCD-level are adjusted using a standard weighted linear least-squares (LLS) fit. The error model with non-diagonal covariance matrix as derived in Tanga et al. (2023) provides the weights to be applied for each observation equation. The same set of asteroids as in Gaia DR3 is considered for this new release, but over an observational time span that is roughly twice as long.  \nFor the LLS method, the theoretical expectation is required. To spell out this function, the equations of motion are integrated simultaneously with the variational equations. In the ICRF reference frame and on the TCB timescale, the heliocentric equation of motion is expressed as  \nr i = -( GM ⊙ + Gmi ) r i r 3 i + X p f p , (1)  \nwhere ri = | r i | and r i is the heliocentric position vector of asteroid i , GM ⊙ is the mass parameter of the Sun, and Gmi is the mass parameter of asteroid i , which can generally be neglected. The additional perturbing accelerations, added linearly, correspond to  \nX p f p = f p | Grav + f p | relat + · · · (2)  \nThey account as in the Gaia DR3 orbital fitting for the planetary perturbations and the relativistic terms, respectively. The part due to the gravitational perturbations (by a total of np planets, and potentially other perturbing asteroids) is given by  \nf p | Grav = -np X j , i , j ≥ 1 Gmj        r i -r j r 3 i j + r j r 3 j        . (3)  \nThe relativistic contribution to the equations of motion can be approximated by  \nf p | relat ≃ GM ⊙ c 2 r 3 i (" GM ⊙ ri ˙ r 2 i # r i + 4( r i · ˙ r i ) ˙ r i ) . (4)  \nArticle number, page 8 of 22  \nThe variational equations are generated by the derivation of Eq. (1) with respect to the parameter to be corrected. Beutler (2005) and Pontriaguine (1969) give a very thorough exposition of how the variational equations can be obtained from a general point of view.  \nThe heliocentric positions of the asteroids are computed from numerical integration 2 of Eq. (1). In Gaia DR3, the Solar System model is the planetary solution INPOP10e 3 . In the current solution, the Solar System model has been updated to INPOP19a, which presents some important di ff erences with INPOP10e. As with INPOP10e, the Earth-Moon barycentre with the seven other planets, dwarf planet Pluto, and a selection of 343 asteroids are included in the Solar System model (Fienga et al. 2013; Deram et al. 2022). For INPOP19a, the ten most massive trans-neptunian objects (TNO) have additionally been included as perturbing bodies. Furthermore, in INPOP19a, influences by the less massive TNOs are modelled by a ring. These added perturbations in INPOP19a cause a shift in the Solar System barycentre that a ff ects the barycentric position of the Gaia satellite (Fienga et al. 2019).  \nNo perturbation of asteroids is included yet in Eq. (3). This is still a limitation compared to other databases of asteroid orbits, which consider systematic perturbations by several massive asteroids (e.g. Ceres, Pallas, or Vesta) together with the major planets, or for consistency with INPOP19a. This approximation on the dynamical model introduces a small bias on the orbital solution (quantified in internal validations; see Sect. 5.3.1), but almost no e ff ect on the uncertainty of the orbits themselves. Thus we will compare and analyse the formal uncertainties of the orbital solution, not its accuracy. These systematic perturbations will be included in subsequent data releases (starting from DR4), together with other mutual perturbations, that should also allow the mass determination from a global inversion.  \nBecause mutual interactions between asteroids are neglected, equations (1) and (3) can be integrated simultaneously with the variational equations for each asteroid separately. The small corrections relative to the reference orbit are then computed at a given reference epoch by an LLS method that minimises the residuals on the astrometric position. The correction is determined for the six-dimension initial state-vector in Cartesian coordinates taken from an auxiliary file of orbital elements. In our case, it was derived from the astorb database 4 and was propagated to the reference epoch. The reference epoch for the orbital fit was taken to be the mid-point of the total time span of the observations available for the given asteroid as a proxy for the weighted mean time of the observational arc. This reference epoch is thus not common to all asteroids in the orbital data base.  \nTo be more concise, we must solve the following linear system of equations for each asteroid source:  \n√ W d λ = √ W                               ∆ α 0 cos( δ 0) ∆ δ 0 ∆ α 1 cos( δ 1) ∆ δ 1 . . . ∆ α N -1 cos( δ N -1) ∆ δ N -1                               = √ WA . d q , (5)  \nwhere the vector d λ ( t ) = ( O -C ) represents the value of the difference between the observed and computed values of the measured quantity, in this case, the right ascension and declination for each of the N observations for a given asteroid. d q is the differential correction to the parameters, here the components of the initial state-vector for the asteroid under consideration, and A is the matrix formed with the partial derivatives with respect to these parameters, projected to right ascension and declination. Furthermore, all equations are weighted by the matrix √ W obtained through the corresponding variance or covariance matrix for the right ascension and declination of the observations. √ W is computed as described in Levinger (1980). We accepted only real values and used only the positive square roots where these are needed in the formulae. Moreover, because the variance or covariance matrix for right ascension and declination of the observations is positive semi-definite and symmetric, √ W verifies √ √  \nW T W = W . We refer to Tanga et al. (2023) for a detailed explanation of the uncertainties involved in the data reduction. As explained therein, each observation of an asteroid at a given epoch, that is, a right ascension α and declination δ , has an associated variance or co-variance matrix of the form  \nγ = var( α ∗ ) covar( α ∗ , δ ) covar( α ∗ , δ ) var( δ ) ! = σα ∗ α ∗ σα ∗ δ σα ∗ δ σδδ ! , (6)  \nwith α ∗ = α cos ( δ ) in di ff erential quantities. Then the weights are obtained as W = γ -1 k for the kth asteroid. We recall that the o ff -diagonal elements play a prominent role in the uncertainty deduced for the sky coordinates, as already stressed in Section 3.  \nPlanetary aberration as well as the gravitational light bending are then added to the computed coordinates as so-called corrections to the observations to construct the (O -C) vector in Eq.(5). All corrections are computed in the general relativistic framework of Gaia data reduction, considering the source at finite distance (Klioner 2003). While the relativistic light deflection e ff ect, which can reach several milliarcseconds at moderate solar elongations, is not included in ground-based observations, it is mandatory in the treatment of Gaia data due to the submilliarcsecond astrometry in the position of the asteroids. We recall that the published right ascension and declination are corrected for the annual aberration, that is, the equivalent to the stellar aberration due only to the barycentric motion of the observer, without any contribution from the asteroid motion.  \nThe LLS is based on the optimisation of a target function, which in our case is defined by  \nQ = 1 Ni Ni -1 X j = 0 ( O -C ) T j W j ( O -C ) j , (7)  \nthat is, the sum of the squares of the weighted residuals (Milani & Gronchi 2010) of the superscript T meaning the transpose of the vector ( O -C ) j . As in the Gaia DR3 release, optimisation is performed for each asteroid individually, considering all astrometric observations N . There are no more general parameters that would couple observations from di ff erent  \nplanets in the fit. Thus, the LLS solution of the system in d ¯ q = ( A T WA ) -1 A T W d λ with corresponding variance / covariance matrix ( A T WA ) -1 where A T indicates the transpose of the matrix A .', '4.4. Iterations and convergence': 'The procedure of orbit correction is iterated until convergence is reached with a given tolerance p following Milani & Gronchi 2010, sec. 5.2; initially this parameter has a value of 10 -8 however, it can be increased in some cases. A preliminary outlier filtering is performed during the astrometric reduction (section 3), further outlier rejection at the observation level has been implemented and follows Carpino et al. (2003).  \nThe procedure is built around a post-fit χ 2 computed as χ 2 = ξ j γ -1 ξ j ξ T j ; with ξ the residuals of all observations evaluated with respect to the corresponding orbit and γ ξ j their expected covariance matrix. The threshold for ξ is a parameter of the fitting process. The tolerance test for stopping the computation is evaluated between iteration k and k + 1 as follows:  \nc 1 = | Qk + 1 -Qk | Qk + 1 < p , (8)  \nwhere p is a tolerance threshold, and Qk is given by Eq.(7) for the kth iteration.  \nIf the next iteration is unlikely to improve the corrections, the processing can also be stopped; this can be assessed by considering the size of the last correction using an appropriate norm (see Milani & Gronchi (2010)), for example,  \nc 2 = ∥ ∆ X ∥ = r ∆ X T M ∆ X N < p , (9)  \nwhere M = A T WA is the normal matrix, ∆ X is the correction vector, and N is the number of observations. In practice, this criterion is used far less frequently. In addition to the two criteria mentioned above, an absolute maximum number of iterations is also denoted nmax .  \nThe iteration starts with p = 10 -8 , and the maximum number of iterations is set to nmax = 15. The value of p = 10 -8 is somewhat arbitrary, and di ff erent values could be imposed (or even a di ff erent tolerance for each condition c 1 and c 2). However, according to our tests, it is unlikely that this would significantly a ff ect the overall statistical outcome of the orbit adjustment because typically, only a few iterations are required. When neither c 1 nor c 2 is lower than p after 15 iterations, the iteration is continued and p is replaced by 10 p . When another 15 iterations do not return a c 1 or c 2 that is lower than 10 p , this is repeated one more time for a last trial. In general, convergence is reached after a few iterations, and in practice, the value of p had to be changed very rarely.  \nThree situations are considered as a failure in the attempt to fit an orbit. First when the maximum number of iterations is reached, second if the orbital solution is not an ellipse and is rejected and third when all observations have finally been rejected by the rejection algorithm.  \nA total of 59 sources were rejected, but their astrometry still appears in the FPR data; these are identified with a NULL state vector in the tables. Interestingly, 31 of these 59 sources are moons of Mars, Jupiter, Saturn, Uranus, or Neptune. An attempt at orbit improvement was not even performed for these objects. Twenty-four of the remaining 28 asteroids are main-belt asteroids, three are near-Earth objects ((433) Eros, 2001 XV10, and  \nFig. 4. Ratio of the arc length of the FPR sources to that of the Gaia DR3 ( R ). The peak at ≃ 2 is marked, and the excursion about the mode is moderate.  \n<!-- image -->  \n2003 YV37), and the last source is a Mars crosser (2006 PG1). For two of the main-belt objects, 2005 UB160 and 2001 QG115, the high rejection rate that is probably due to a larger fraction of outliers makes convergence unstable or the orbit computation based on overly short arcs unreliable. No convergence was reached. A dedicated investigation of this subset of asteroids, which consists of only 0 . 2% of the orbital catalogue, is needed to fully understand why convergence has not been reached. This is beyond our concerns in this FPR, however, because (433) Eros is a di ffi cult case in Gaia DR3 and requires further investigations.', '4.5. Post-fit residuals': 'The full processing resulted in 156 762 orbits for the FPR. In contrast to Gaia DR3, there is no post-processing of the FPR orbit determination. In the previous catalogue, a further filtering of solutions was applied based on the relative uncertainty in the semi-major axis σ a / a and arc length. We refer to Tanga et al. 2023, Sect. 3.5 for an explanation.  \nThe time span covered by the FPR observations is roughly twice that of Gaia DR3. Although here the astrometric processing is more accurate (see Section 3), it was a reasonable guess to expect that the arc covered by any object would be at least as long as it was in the Gaia DR3 release. This is true for an overwhelming majority of the sources (99%) and even 76% having the ratio R of the total arc length of the FPR sample to that of the Gaia DR3 at least equal to 2. Figure 4 shows the normalised distribution of R where a sharp peak is seen around R = 2 followed by a steep decline. A substantial fraction of the sources have a ratio R > 2, while the global time coverage from Gaia has changed from 34 to 66 months. These happen to be sources with actual time coverage, between first and last observation, much below 34 months in Gaia DR3, thus increasing automatically the ratio in FPR. One can see also that there are very few sources with a ratio less than one, i.e. an unexpected longer arc in the Gaia DR3 data than in the FPR sample. This occurs when the orbit adjustment rejects a large part of the observations, resulting into a smaller orbital arc; only eight sources (over approx. 157 000) belong to this group. Figure 5 displays the along scan residuals against the G magnitude, Gmag . The 5 σ cut-o ff value for rejection is visible; it is the envelope of the accepted observations, i.e. the limit between the grey and coloured points. Figure 6 shows the corresponding histograms for normalised residuals ∆ AL /σ AL for sources with a Gmag less than or equal to 11 . 5, or greater than  \nFig. 5. Along-scan residuals against G magnitudes. The grey points represent all observations. The coloured points are the observations that were accepted by the orbit adjustment procedure after outlier rejection chosen at | O -C | > 5 σ . For Gmag < 11 . 5 , the distribution is not Gaussian. This is clearer in Figure 6. For Gmag ≥ 11 . 5 , the distribution is Gaussian. The orange lines in the plot are ± 1 σ , the blue lines are ± 2 σ , and the yellow lines are ± 3 σ . The number density is given by the scale on the right.  \n<!-- image -->  \nFig. 6. Histograms of the along-scan residuals normalised to the formal uncertainties for the accepted solutions with Gmag ≤ 11 . 5 shown in blue or Gmag > 11 . 5 shown in green. The histogram in magenta shows all sources (accepted or rejected) with Gmag ≤ 11 . 5. A Gaussian fit to the data is shown as well. We note that for sources with Gmag > 11 . 5 , the distribution is very well represented by a Gaussian ( µ = 0 . 018, σ = 1 . 08), but for Gmag ≤ 11 . 5 , the distribution is clearly not Gaussian.  \n<!-- image -->  \n11 . 5, respectively. For accepted observations with Gmag greater than 11 . 5, the normalised distribution is Gaussian with a mean of 0 . 018 and standard deviation of 1 . 08; for these sources the orbit adjustment is acceptable. Sources with Gmag less than 11 . 5 do not show a Gaussian distribution, even if the rejected solutions are included; in this magnitude range the model is not entirely satisfactory, but can be improved by including some other effects, as the photocentre o ff set for example. For Gmag > 11 . 5 the symmetry is striking.  \nFigure 7 displays the residuals in along and across scan. There is a small o ff set from zero, ( ∆ AL , ∆ AC ) ≈ (0 , -15) mas or ( ∆ AL /σ AL , ∆ AC /σ AC ) ≈ (0 . 05 , -0 . 028), this is also the true  \nFig. 7. Post-fit residuals of the across-scan direction vs. the along scan direction. Top: True angular values. Bottom: Normalised quantities. The orange lines show the median, and the yellow lines show the 5th (left) and 95th percentile (right). The density is given in log-scale, so that only the core of the plots is populated.  \n<!-- image -->  \nin Gaia DR3. In the ∆ AC /σ AC vs ∆ AL /σ AL plot the distribution is clearly asymmetric about ∆ AC /σ AC = 0 but appears to be much more symmetric about ∆ AL /σ AL = 0. Again, the same is seen in the Gaia DR3 sample and is probably due to the greater uncertainties in the across scan direction. For 90% of the sources ∆ AL is in the interval [ -8 , + 8] mas, ∆ AL /σ AL ∈ [ -1 . 8 , + 1 . 8].  \nArepresentation of the astrometry quality over a single transit can finally be obtained by considering the dispersion of the measurements that were recorded during each passage of a source in the focal plane in the along-scan direction (Fig. 8). This figure can be compared to Fig. 12 in Tanga et al. (2023) for Gaia DR3. The figures appear nearly equivalent, with a slight improvement around the G-magnitude interval of the best performance (Gmag ∼ 12-13).', '4.6. Orbital elements': 'The output of the orbit fitting for a particular asteroid consists of a 6D vector of the initial conditions at a particular reference epoch, the midpoint of the observation times, and is given in TCB. This vector combines a position and a velocity vector and is heliocentric. It is given on the ICRF axes. The transformation from the state vector to the osculating elements is in principle routine work in celestial mechanics, and fully tested software is widely available. However, here we work with sub-  \nFig. 8. Standard deviation of the values of the along-scan residuals for each CCD, computed over each transit. The colour corresponds to the point density. The lines represent the smoothed average, and the quantiles at the 1 σ level.  \n<!-- image -->  \nmilliarcsecond astrometry, or with orbits at the 10 -9 -10 -10 accuracy level ( as σ ( a ) / a ), and care must be taken everywhere to ensure that no loss of accuracy is associated with the transformation and that the underlying conventions are well understood and implemented. This is the object of this section, to explain all the conventions we used in our solution and describe where they may di ff er from other usages.  \nThe transformation from the state vectors to the osculating elliptical elements is done in two steps: first, the transformation from heliocentric ICRF frame heliocentric ecliptic frame for the Cartesian coordinates and then the transformation in the ecliptic frame from state vector to osculating elements. In addition, the TCB is used consistently for Gaia while all other sources of osculating elements used TDB in their publications and this impacts the scale length.  \nThe transformation of the state vector to the ecliptic frame depends on just one rotation matrix to connect the two frames. The main part is the rotation ϵ around the x -axis, where ϵ is the obliquity of the ecliptic. The devil lies in the details, as so often is the case. The intersection of the dynamical ecliptic with the ICRF fundamental plane has no reason to coincide with the ICRF x -axis, and this is not the case, either. In addition, the obliquity of the ecliptic relates the ecliptic (taken as a dynamical plane with a strict definition) to the celestial equator. The latter is again di ff erent from the ICRF reference plane, and the angle between the two planes of interest is not exactly the physical obliquity of the ecliptic. It is very similar, but not equal. As a consequence, and for lack of a well-agreed convention, several definitions are in use. From a single state vector, we can derive six sets of osculating elements that di ff er by several tens of milliarcseconds in the longitude of the ascending node or the orbital inclination. This motivated our choice to publish Gaia orbits in the form of a heliocentric state vector in the ICRF axes instead of the usual osculating elements. The former is unambiguous, and from this vector, the orbital elements can be computed so that they agree with the di ff erent choices for the ecliptic. They di ff er slightly with the choice of the ecliptic, although they represent the same orbit.  \nThe issue extends beyond the published osculating elements to the users of this information. They may start from osculating elements of an asteroid to propagate to another time and compute local coordinates using the true equator of the date. Then they will call a precession or nutation package, relating the J2000 celestial equator to another epoch, or the ICRF plane  \n<latexit sha1\\_base64="R+NicdZlo8K16nDvJaa4pChb/2Y=">AAAHmnicnVXbbtNAEJ22QEq5NAXe4MGQVOKhiRLUUh4reADES5GatlITItvZpFZ8k71uGyz/C6/wR/wB/AVnJ5uQewKOHM+enXN2dma8tkLXiWWl8nNtfePW7Tu5zbtb9+4/eLid33l0GgdJZIuaHbhBdG6ZsXAdX9SkI11xHkbC9CxXnFndd2r+7EpEsRP4J7IXioZndnyn7dimBNTMP6l3TM8zm3UpbmRa33tvOmbWzBcq5QpfxrRR1UaB9HUc7GwcUp1aFJBNCXkkyCcJ2yWTYvwuqEoVCoE1KAUWwXJ4XlBGW+Am8BLwMIF28d/B6EKjPsZKM2a2jVVc3BGYBu2CE8Avgq1WM3g+YWWFztNOWVPF1sPT0loeUEmXQJfxBp6r8yzEFWAsgVjwmL/vdKgnmeNq7925/tec557ezeKMSmrTG86kg8yGjKgc23ovCddG5c8Yya2KOgSm7BbmI9g2MwfVNpgTcwVUhU2e/8WeClVjW/sm9FvnygdyzTXzeNc+1FPgSsEGomLMGLnkDN7ALtIXjEN4OLx6ER7ZCmoxxzxLLR3TG7WzFdTbXJ9J1TLGz/nu60fgGOy5imIL9VqsWQfP4X2tojfI52xNDx2iautTifYxSv45w+1hjpdn4v8zHSOyWf1Qh6VOH7lSLySs0NezuTOTpawW8A6fOWqkujjmN0MM39V0rCYDPcVrM97v/JSzrc6n/gmo3m51YvqYqbL6NNcDNp+3OO4un1b9qPd4pJ6x7voSr1lcohHQFXeHxTXO+DRv6IgGfWMAK9MBx/iX4fL7KyZ8SyO+Bdjz5rK5K6gvVHXyezRtnL4qV1+XDz7vF44q+lu1SU/pBb2E1iEd0Qc6php2/JW+0Xf6kXuWe5v7mPvUd11f05zHNHblTv4Aek6CTg==</latexit>  \n<latexit sha1\\_base64="EUUBAwplnzGTJ/g4ElaQmnc+sWU=">AAAHnHicnVXbbtNAEJ22QEq5pSDxgoQMSSUemshBLeWxEkggIaQikbRSEyrb2aRWfJO9bhsZ/wyv8EP8AfwFZyebkKS5gSPHs2fnnJ2dGa/tyHMTaZo/19Y3bty8Vdi8vXXn7r37D4rbDxtJmMaOqDuhF8YntpUIzw1EXbrSEydRLCzf9sSx3Xuj5o8vRJy4YfBZ9iPR8q1u4HZcx5KAzoqPmyJKXA9mU4ormTV331mulZ8VS2bV5Mu4btS0USJ9HYXbGwfUpDaF5FBKPgkKSML2yKIEv1OqkUkRsBZlwGJYLs8LymkL3BReAh4W0B7+uxidajTAWGkmzHawioc7BtOgHXBC+MWw1WoGz6esrNB52hlrqtj6eNpaywcq6RzoMt7Qc3WejbhCjCUQGx7z952N9CRzPO29M9f/kvPc17tZnFFJHXrNmXSR2YgRlWNH7yXl2qj8GWO5VVFHwJTdxnwM22HmsNoGcxKugKqwxfO/2FOhauxo35R+61wFQC65Zj7vOoB6BlwpOEBUjDkj55zBK9hl+oJxBA+XVy/DI19BLeGYZ6llE3rjdr6CeofrM61axfgZ3wP9GByDPVdRbKNeizWb4Lm8r1X0hvmcremjQ1RtA6rQHkbpP2e4M8rx8kz8f6YTRDarH5qw1OkjV+qFlBUGeg53ZrqU1Qbe5TNHjVQXJ/xmiNG7mk3UZKineB3GB52fcbbV+TQ4AdXbrU7MADM1Vr/O9YHN5y2Ou8en1SDqXR6pZ6K7vsJrlpdohHTB3WFzjXM+zVs6omHfGMCqtM8x/mV4/P6KKd/KmG8J9ry5fO4K6gtVm/4eXTcaL6u1V9X9T3ulQ1N/qzbpCT2nF9A6oEN6T0dUx46/0jf6Tj8KTwtvCx8KHweu62ua84gmrkLjD8Qsg1k=</latexit>  \n<latexit sha1\\_base64="8BX89Oou3gKmdTnLX3hi9wKQoLI=">AAAHiHicnVXbbtNAEJ22gEu4tfDIiyGpxAOJHNQSeKvghccikbZSE5DtbBIrvsmXlsjyL/AKv8YfwF9wdrwJuTkJOHI8e3bO2dmZ8doKXSdODOPnzu7erdt3tP27lXv3Hzx8dHD4+DwO0sgWbTtwg+jSMmPhOr5oJ07iisswEqZnueLCGr2X8xfXIoqdwP+UjEPR9cyB7/Qd20wk1AmHzpeDqtEw+NKXjaYyqqSus+Bwr0Ud6lFANqXkkSCfEtgumRTjd0VNMigE1qUMWATL4XlBOVXATeEl4GECHeF/gNGVQn2MpWbMbBuruLgjMHU6AieAXwRbrqbzfMrKEi3TzlhTxjbG01JaHtCEhkA38Sae2/MsxBVgnACx4FG+72yqlzDHVd5Hpf43nOex2s36jCbUpzecSQeZDRmRObbVXlKujcyfPpNbGXUITNo9zEewbWZOqq0zJ+YKyAqbPP+LPSUqx7byTem3ypUP5IZr5vGufahnwKWCDUTGmDMy5Ax+hV2jzxiH8HB49Ro88i3UYo55lVo2pzdr51uo97k+i6oNjJ/xXehH4OjsuY1iD/Var9kBz+F9baM3yedqTQ8dImvrU52OMUr/OcP9aY43Z+L/Mx0jslX90IElT59kq15IWaHQs7kz042sHvABnzlyJLs45jdDTN/VbK4mEz3J6zNedH7G2ZbnU3ECyrdbnpg+Zpqsvsz1gJXz1sc94tOqiPolj+QzVl1f5zVrGzQCuubusLjGOZ/mXRXRpG90YA064Rj/Mlx+f8WCb33Gtwq7bC4vXSHHF6q5+D1aNs5fNZqvGycfj6unhvpW7dNTek4voNWiU/pAZ9TGjof0jb7TD62iGVpLe1u47u4ozhOau7R3fwAms3tJ</latexit>  \n<latexit sha1\\_base64="zWhuPxn5hCh2xxnXL4J8BkXHKxs=">AAAHmnicnVXbbtNAEJ22QEq5pcAbPBjSSjw0UYJaymMFQqLipUDTVmpCZDub1Ipvstdtg+V/4RX+iD+Av+DsZBNyT8CR49mzc87OzozXVug6sSyXf66srt24eSu3fnvjzt179x/kNx+exEES2aJqB24QnVlmLFzHF1XpSFechZEwPcsVp1bnrZo/vRRR7AT+seyGou6Zbd9pObYpATXyj2tt0/PMRk2Ka5nWdg7fffqcNfKFcqnMlzFpVLRRIH0dBZtr+1SjJgVkU0IeCfJJwnbJpBi/c6pQmUJgdUqBRbAcnheU0Qa4CbwEPEygHfy3MTrXqI+x0oyZbWMVF3cEpkHb4ATwi2Cr1QyeT1hZobO0U9ZUsXXxtLSWB1TSBdBFvL7n8jwLcQUYSyAWPGbvOx3oSea42nt7pv8V57mrdzM/o5Ja9Joz6SCzISMqx7beS8K1UfkzhnKrog6BKbuJ+Qi2zcx+tQ3mxFwBVWGT53+xp0LV2Na+Cf3WufKBXHHNPN61D/UUuFKwgagYM0YuOIPXsLfoC8YhPBxefQse2RJqMcc8TS0d0Ru2syXUW1yfcdUSxs/47ulH4BjsuYxiE/War1kDz+F9LaPXz+d0TQ8domrrU5F2MUr+OcOtQY4XZ+L/Mx0jsmn9UIOlTh+5VC8krNDTs7kzk4WsJvA2nzlqpLo45jdDDN7VdKQmfT3FazHe6/yUs63Op94JqN5udWL6mKmw+iTXAzabNz/uDp9Wvah3eKSese76Iq+5tUAjoEvuDotrnPFpXtcR9fvGAFaiPY7xL8Pl91eM+RaHfAuwZ81lM1dQX6jK+Pdo0jh5Waq8Ku193C0clPW3ap2e0HN6Aa19OqD3dERV7PgrfaPv9CP3NPcmd5j70HNdXdGcRzRy5Y7/ALiJgg8=</latexit>  \n<latexit sha1\\_base64="mk5/6ntOrk6znZYVj1A9Tc3ghWA=">AAAHiHicnVXbbtNAEJ22gEu4tfDIiyGpxAOJHNQSeKvghccikbZSE5DtbBIrvsmXlsjyL/AKv8YfwF9wdrwJuTkJOHI8e3bO2dmZ8doKXSdODOPnzu7erdt3tP27lXv3Hzx8dHD4+DwO0sgWbTtwg+jSMmPhOr5oJ07iisswEqZnueLCGr2X8xfXIoqdwP+UjEPR9cyB7/Qd20wk1Alj58tB1WgYfOnLRlMZVVLXWXC416IO9Sggm1LySJBPCWyXTIrxu6ImGRQC61IGLILl8LygnCrgpvAS8DCBjvA/wOhKoT7GUjNmto1VXNwRmDodgRPAL4ItV9N5PmVliZZpZ6wpYxvjaSktD2hCQ6CbeBPP7XkW4gowToBY8CjfdzbVS5jjKu+jUv8bzvNY7WZ9RhPq0xvOpIPMhozIHNtqLynXRuZPn8mtjDoEJu0e5iPYNjMn1daZE3MFZIVNnv/FnhKVY1v5pvRb5coHcsM183jXPtQz4FLBBiJjzBkZcga/wq7RZ4xDeDi8eg0e+RZqMce8Si2b05u18y3U+1yfRdUGxs/4LvQjcHT23Eaxh3qt1+yA5/C+ttGb5HO1pocOkbX1qU7HGKX/nOH+NMebM/H/mY4R2ap+6MCSp0+yVS+krFDo2dyZ6UZWD/iAzxw5kl0c85shpu9qNleTiZ7k9RkvOj/jbMvzqTgB5dstT0wfM01WX+Z6wMp56+Me8WlVRP2SR/IZq66v85q1DRoBXXN3WFzjnE/zropo0jc6sAadcIx/GS6/v2LBtz7jW4VdNpeXrpDjC9Vc/B4tG+evGs3XjZOPx9VTQ32r9ukpPacX0GrRKX2gM2pjx0P6Rt/ph1bRDK2lvS1cd3cU5wnNXdq7P3Upe1Q=</latexit>  \n<latexit sha1\\_base64="EQcwgZ6Qye2M1eBV5/B3bpPzf+k=">AAAHnHicnVXbbtNAEJ22QEq5NAWJFyRkSCrx0EQJaimPlQAJhJDKJWmlJlS2s0mt+CZ73TYy/hle4Yf4A/gLzk42Ic0dHDmePTvn7OzMeG2FrhPLSuXnyurates3cus3N27dvnN3M791rx4HSWSLmh24QXRsmbFwHV/UpCNdcRxGwvQsVxxZ3Zdq/uhcRLET+J9lLxRNz+z4TtuxTQnoNP+gIcLYcWE2pLiUaWPn7euPn7LTfKFSrvBlTBpVbRRIX4fB1to+NahFAdmUkEeCfJKwXTIpxu+EqlShEFiTUmARLIfnBWW0AW4CLwEPE2gX/x2MTjTqY6w0Y2bbWMXFHYFp0DY4Afwi2Go1g+cTVlboLO2UNVVsPTwtreUBlXQGdBFv4Lk8z0JcAcYSiAWP2ftOh3qSOa723p7pf8F57undzM+opDa94Ew6yGzIiMqxrfeScG1U/oyR3KqoQ2DKbmE+gm0zc1BtgzkxV0BV2OT5X+ypUDW2tW9Cv3WufCAXXDOPd+1DPQWuFGwgKsaMkTPO4CXsIn3BOISHw6sX4ZEtoRZzzNPU0it6o3a2hHqb6zOuWsb4Md99/Qgcgz2XUWyhXvM1G+A5vK9l9Ab5nK7poUNUbX0q0S5GyT9nuD3M8eJM/H+mY0Q2rR8asNTpI5fqhYQV+no2d2aykNUC3uEzR41UF8f8Zojhu5peqclAT/HajPc7P+Vsq/OpfwKqt1udmD5mqqw+yfWAzebNj7vLp1U/6h0eqWesu77EaxYXaAR0zt1hcY0zPs2bOqJB3xjAyrTHMf5luPz+ijHf0ohvAfasuWzmCuoLVR3/Hk0a9Wfl6vPy3ofdwkFFf6vW6SE9oafQ2qcDekOHVMOOv9I3+k4/co9yr3Lvcu/7rqsrmnOfrly5+h8CdoMa</latexit>  \n<latexit sha1\\_base64="VbBcPCd0SBaC5J0PAyOjKSN1Tio=">AAAHm3icnVXbbtNAEJ22QEq5pSCeEJIhqcRDEzlVS3ms1BdU8VAk0lZqQmU7m9SKb7LXbSPLH8MrfBF/AH/B2ckm5J6AI8ezZ+ecnZ0Zr+3IcxNpmj/X1jfu3X9Q2Hy49ejxk6fPitvPz5IwjR1Rd0IvjC9sKxGeG4i6dKUnLqJYWL7tiXO7e6zmz29EnLhh8EX2ItH0rU7gtl3HkoCuii8bHcv3rauGFHcya+ye7JmmmV8VS2bV5MuYNmraKJG+TsPtjUNqUItCciglnwQFJGF7ZFGC3yXVyKQIWJMyYDEsl+cF5bQFbgovAQ8LaBf/HYwuNRpgrDQTZjtYxcMdg2nQDjgh/GLYajWD51NWVug87Yw1VWw9PG2t5QOVdA10GW/guTrPRlwhxhKIDY/5+86GepI5nvbemet/y3nu6d0szqikNn3gTLrIbMSIyrGj95JybVT+jJHcqqgjYMpuYT6G7TBzUG2DOQlXQFXY4vlf7KlQNXa0b0q/da4CILdcM593HUA9A64UHCAqxpyRa87gHewyfcU4gofLq5fhka+glnDMs9SyMb1RO19Bvc31mVStYvyG775+DI7BnqsotlCvxZoN8Fze1yp6g3zO1vTRIaq2AVVoH6P0nzPcHuZ4eSb+P9MJIpvVDw1Y6vSRK/VCygp9PYc7M13KagHv8JmjRqqLE34zxPBdzcZqMtBTvDbj/c7PONvqfOqfgOrtVidmgJkaq09zfWDzeYvj7vJp1Y96l0fqmeiur/Ca5SUaId1wd9hc45xP86aOaNA3BrAqHXCMfxkev79iwrcy4luCPW8un7uC+kLVJr9H08bZXrX2vnrweb90ZOpv1Sa9orf0DlqHdEQf6ZTqnN1v9J1+FF4XjgsnhU991/U1zXlBY1eh/gfTi4Hy</latexit>  \n<latexit sha1\\_base64="jqZEUdY8rjmnLlEhga75apbzmto=">AAAHkXicnVXbbtNAEJ22QEq5tfSRF0NSiQcSOail8BZRCYF4oCDSVmpCZTub1Ipv2Ou2keXv4BU+iz+Av+DsZBOae8CR49mzc87OzozXduS5iTTNnyurazdu3iqs3964c/fe/QebWw+PkjCNHVF3Qi+MT2wrEZ4biLp0pSdOolhYvu2JY7t7oOaPL0ScuGHwWfYi0fStTuC2XceSgJofzhpSXMns3cGnN/nZZtGsmHwZk0ZVG0XS12G4tbZPDWpRSA6l5JOggCRsjyxK8DulKpkUAWtSBiyG5fK8oJw2wE3hJeBhAe3iv4PRqUYDjJVmwmwHq3i4YzAN2gEnhF8MW61m8HzKygqdpZ2xpoqth6ettXygks6BLuINPJfn2YgrxFgCseExe9/ZUE8yx9PeOzP9LznPPb2b+RmV1KaXnEkXmY0YUTl29F5Sro3Kn3EttyrqCJiyW5iPYTvMHFTbYE7CFVAVtnj+F3sqVI0d7ZvSb52rAMgl18znXQdQz4ArBQeIijFn5JwzeAW7RF8wjuDh8uoleORLqCUc8zS1bETvup0vod7m+oyrVjB+zHdfPwbHYM9lFFuo13zNBngu72sZvUE+p2v66BBV24DKtItR+s8Zbg9zvDgT/5/pBJFN64cGLHX6yKV6IWWFvp7DnZkuZLWAd/jMUSPVxQm/GWL4rmYjNRnoKV6b8X7nZ5xtdT71T0D1dqsTM8BMldUnuT6w2bz5cXf5tOpH/YxH6pnori/zmqUFGiFdcHfYXOOcT/OmjmjQNwawCu1xjH8ZHr+/Ysy3fM23CHvWXD5zBfWFqo5/jyaNo+eV6ovK3sfdYs3U36p1ekRP6Cm09qlGb+mQ6tjxV/pG3+lHYbvwqlArvO67rq5ozjaNXIX3fwBRvn7x</latexit>  \n<latexit sha1\\_base64="E9sOUlp1xXmW5vYMxcNT0xz6Tvc=">AAAHnXicnVXbbtNAEJ22QEq5pfDAAw8Ykko8NJFTtZTHSggJIYSKRJpKTahsZ5Na8U32um1k+Wt4hQ/iD+AvODvZhCTNDRw5nj075+zszHhtR56bSNP8uba+cev2ncLm3a179x88fFTcfnyShGnsiLoTemF8aluJ8NxA1KUrPXEaxcLybU807N5bNd+4FHHihsEX2Y9Ey7e6gdtxHUsCOi8+bYoocT2YTSmuZdbc/bBnmmZ+XiyZVZMv46ZR00aJ9HUcbm8cUpPaFJJDKfkkKCAJ2yOLEvzOqEYmRcBalAGLYbk8LyinLXBTeAl4WEB7+O9idKbRAGOlmTDbwSoe7hhMg3bACeEXw1arGTyfsrJC52lnrKli6+Npay0fqKQLoMt4Q8/VeTbiCjGWQGx4zN93NtKTzPG0985c/yvOc1/vZnFGJXXoDWfSRWYjRlSOHb2XlGuj8meM5VZFHQFTdhvzMWyHmcNqG8xJuAKqwhbP/2JPhaqxo31T+q1zFQC54pr5vOsA6hlwpeAAUTHmjFxwBq9hl+krxhE8XF69DI98BbWEY56llk3ojdv5Cuodrs+0ahXjF3wP9GNwDPZcRbGNei3WbILn8r5W0Rvmc7amjw5RtQ2oQvsYpf+c4c4ox8sz8f+ZThDZrH5owlKnj1ypF1JWGOg53JnpUlYbeJfPHDVSXZzwmyFG72o2UZOhnuJ1GB90fsbZVufT4ARUb7c6MQPM1Fj9JtcHNp+3OO4en1aDqHd5pJ6J7voKr1leohHSJXeHzTXO+TRv6YiGfWMAq9IBx/iX4fH7K6Z8K2O+Jdjz5vK5K6gvVG36e3TTONmr1l5XDz7vl45M/a3apGf0kl5B65CO6D0dU50z/o2+04/C88K7wsfCp4Hr+prmPKGJq9D4Ax5zgv0=</latexit>  \n<latexit sha1\\_base64="3iI+NoJoDxNeODpJXkqLEddOEaU=">AAAHmXicnVXbbtNAEJ22QEq5NKWPfTEklXhoowS1lMciXqqKhyBIW6kJke1sUiu+yV63jSx/C6/wSfwB/AVnJ5uQewKOHM+enXN2dma8tkLXiWW5/HNtfePBw0e5zcdbT54+e76d33lxEQdJZIuaHbhBdGWZsXAdX9SkI11xFUbC9CxXXFrdD2r+8lZEsRP4X2QvFA3P7PhO27FNCaiZ3613TM8zm3Up7mVaPzivfsya+UK5VObLmDYq2iiQvqrBzsYJ1alFAdmUkEeCfJKwXTIpxu+aKlSmEFiDUmARLIfnBWW0BW4CLwEPE2gX/x2MrjXqY6w0Y2bbWMXFHYFp0D44Afwi2Go1g+cTVlboPO2UNVVsPTwtreUBlXQDdBlv4Lk6z0JcAcYSiAWP+ftOh3qSOa723p/rf8d57undLM6opDa940w6yGzIiMqxrfeScG1U/oyR3KqoQ2DKbmE+gm0zc1BtgzkxV0BV2OT5X+ypUDW2tW9Cv3WufCB3XDOPd+1DPQWuFGwgKsaMkRvO4D3sIn3FOISHw6sX4ZGtoBZzzLPU0jG9UTtbQb3N9ZlULWH8ku++fgSOwZ6rKLZQr8WadfAc3tcqeoN8ztb00CGqtj4d0hFGyT9nuD3M8fJM/H+mY0Q2qx/qsNTpI1fqhYQV+no2d2aylNUC3uEzR41UF8f8Zojhu5qO1WSgp3htxvudn3K21fnUPwHV261OTB8zFVaf5nrA5vMWx93l06of9QGP1DPWXX/IaxaXaAR0y91hcY0zPs0bOqJB3xjASnTMMf5luPz+ignfwxHfAux5c9ncFdQXqjL5PZo2Lt6UKm9Lx5+OCqdl/a3apD16Ra+hdUKndEZVqmHHPfpG3+lHbi/3PneWO++7rq9pzi6NXbnPfwBhuIG4</latexit>  \n<latexit sha1\\_base64="sLNxx8+C8ZGdht3leYmno9B1bcY=">AAAHm3icnVXbbtNAEJ22gEu5pSCeEJIhqcRDEyWopTxW6guq+hAk0lZqQmU7m9SKb7LXbSPLH8MrfBF/AH/B2ckmJGlu4Mjx7Nk5Z2dnxms78txEVqs/19Y37t1/YGw+3Hr0+MnTZ4Xt56dJmMaOaDihF8bntpUIzw1EQ7rSE+dRLCzf9sSZ3TtS82fXIk7cMPgi+5Fo+VY3cDuuY0lAl4WXTRElrgezKcWtzJq7x/WT/LJQrFaqfJl3jZo2iqSveri9cUBNalNIDqXkk6CAJGyPLErwu6AaVSkC1qIMWAzL5XlBOW2Bm8JLwMMC2sN/F6MLjQYYK82E2Q5W8XDHYJq0A04Ivxi2Ws3k+ZSVFTpPO2NNFVsfT1tr+UAlXQFdxht6rs6zEVeIsQRiw2P+vrORnmSOp7135vrfcJ77ejeLMyqpQx85ky4yGzGicuzovaRcG5U/cyy3KuoImLLbmI9hO8wcVttkTsIVUBW2eP4XeypUjR3tm9JvnasAyA3XzOddB1DPgCsFB4iKMWfkijN4C7tEXzGO4OHy6iV45CuoJRzzLLVsQm/czldQ73B9plUrGL/he6Afg2Oy5yqKbdRrsWYTPJf3tYreMJ+zNX10iKptQGXawyj95wx3Rjlenon/z3SCyGb1QxOWOn3kSr2QssJAz+HOTJey2sC7fOaokerihN8MMXpXs4maDPUUr8P4oPMzzrY6nwYnoHq71YkZYKbG6ne5PrD5vMVx9/i0GkS9yyP1THTXl3nN0hKNkK65O2yucc6neUtHNOwbE1iF9jnGvwyP318x5Vse8y3CnjeXz11BfaFq09+ju8bp+0rtQ2X/817xsKq/VZv0it7SO2gd0CF9ojo1OLvf6Dv9MF4bR8axcTJwXV/TnBc0cRmNP6qbgsM=</latexit>  \nFig. 9. Relative positions of the ICRF and J2000 equators and of the ecliptic with the di ff erent conventions found for the ecliptic frame. The origin of the ICRF ( O ICRF) is o ff set from the inertial equinox γ J2000. The relevant numbers are given in Table 1.  \n<!-- image -->  \nto the celestial equator by using the frame bias above the precession. Therefore, the reference plane and the origin of longitudes for the osculating element are key ingredients for high-precision computation. In this context, high precision means better than 50 mas, which is the magnitude of the small angles in Table 1 or of their di ff erence between institutes.  \nWe have attempted to collate the various choices made for the reference plane and the origin of the longitude of node in this plane by the most relevant sources of osculating elements. Figure 9 shows the di ff erent equatorial planes and ecliptics, together with their associated obliquity. The angle ϕ is the value of γ xx O ICRF (xx: ICRF, Gaia, JPL) measured in the ICRF equator and taken positive, like right ascension from the local γ xx (the sign convention is opposite to Chapront et al. (2002) and similar to Hilton & Hohenkerk (2004)). The so-called IERS values are taken from the SOFA software and come from the precession formulae derived in Fukushima (2003). The angle ψ is the angular distance between the intersection of the ecliptic with the ICRF plane and the intersection with the J2000 equator, positive like an ecliptic longitude. The numerical values of the angles needed to make the transformation from the ICRF to the selected ecliptic are given in Table 1. The Gaia o ff set was taken from Chapront et al. (2002) before the adoption of the Fukushima precession by IERS and SOFA, and has not yet been changed for continuity reasons. The ecliptic plane and origin in the ICRF equator are extremely similar (within 3 mas) to the VLBI-derived ecliptic adopted in SOFA. However, the origin of the longitude of node is not on the J2000 equator, but in the ICRF equator.  \nTable 1. Relevant angles defined in Figure 9 that were used to relate the inertial ecliptic to the ICRF fundamental plane.  \nϵ = 84381 \'\' + ∆ ϵ  \nNothing is fundamentally wrong with these multiple realisations of the ecliptic, this is just inconvenient. This all stems from the fact that the ecliptic concept is ambiguous when high accuracy is required because the plane orthogonal to the mean angular momentum of the Earth-Moon system depends on what is included in this mean, over which time it is averaged, whether long-period wiggles from planetary perturbations are removed or not, and so on. A less physical definition can appear as well with a rigid link to the ICRF based on conventional angles, in a similar way as the Galactic coordinates are related to the equatorial frame at J2000, with no attempt to recreate a physical definition accurately.  \nIn principle, an unambiguous definition together with numerical constants have been provided in 2006 by the IAU Working Group on Precession and Ecliptic (Hilton et al. 2006) and the ensuing IAU Resolution B1 in 2006. However, a unique realisation is currently not agreed, and even within the Gaia Consortium, the ecliptic has been inherited from Hipparcos, and change in standards in a large data-processing system is always weighed against the loss of continuity and is decided only when unavoidable. It is important for the users to be aware of these di ff erences and to be on their guard when using and propagating orbits with a target accuracy lower than 100 mas. The di ff erences between osculating elements may be deceptive and may not reveal real di ff erences in the orbits. When properly transformed into state vectors in the ICRF, they may look much more similar than from their elliptic elements.  \nTable 1 clearly shows that all the major providers of orbits use the same definition, in which the ecliptic is very simply related to the ICRF equator by a single rotation about the ICRF Xaxis, with an angle equal to the obliquity ϵ = 84381 . 4480 \'\' . We have adopted this convention in the following sections to compute the orbital elements and compare them to other solutions. However, to let the users select their preferred option, and above all, to avoid incorrect assumptions about the ecliptic that is used for the orbital elements, we decided that primary Gaia results will be published as state vectors in ICRF instead of osculating elements. Nevertheless, the covariance matrix between osculating elements is provided because it is independent (to ≈ 10 -8 ) of the small di ff erences between the ecliptic frames of Table 1, and its computation from the state vector could be tricky without an easily available standard transformation routine.  \nIncidentally, the Gaia astrometric data for asteroids and Solar System objects themselves can provide the link between the kinematically non-rotating reference frame (ICRF) and a dynamically non-rotating reference frame (e.g. ECJ200, represented by this inertial ecliptic from the planetary solution, or an invariant plane).', '5.1. Overall statistics': 'The same set of asteroids as was selected for Gaia DR3 was used in this processing, but over a time span of 66 months instead of the 34 months for Gaia DR3. The number of asteroids with validated astrometry that were processed this pipeline is 156 825 , and we obtained converged solutions for 156 762 asteroids with an accepted state vector, or equivalently, a set of six orbital elements, at a particular reference epoch. The success rate may seem particularly high, and it is, but this results from the fact that the sources selected for this run were precisely those with a successfully computed orbit with the Gaia DR3 data, but with a shorter time range. We described the filters  \n<latexit sha1\\_base64="jPlLdmLO/S8BASHEhB7kbcINwiE=">AAAHinicnVXbbtNAEJ22gEsp0MIjL4akEg80cqqWgnipBA88FomklZqAbGeTmPgmX1oiy//AK/wZfwB/wdnxJuTmJODI8ezZOWdnZ8ZrK3SdODGMnxubW7du39G27+7c273/4OHe/qNmHKSRLRp24AbRpWXGwnV80UicxBWXYSRMz3LFhTV4K+cvrkUUO4H/MRmGou2ZPd/pOraZAGq2TDfsm5/3KkbN4EufN+rKqJC6zoP9rVNqUYcCsikljwT5lMB2yaQYvyuqk0EhsDZlwCJYDs8LymkH3BReAh4m0AH+exhdKdTHWGrGzLaxios7AlOnA3AC+EWw5Wo6z6esLNEy7Yw1ZWxDPC2l5QFNqA90FW/kuT7PQlwBxgkQCx7l+87GeglzXOV9UOp/w3keqt0sz2hCXXrFmXSQ2ZARmWNb7SXl2sj86RO5lVGHwKTdwXwE22bmqNo6c2KugKywyfO/2FOicmwr35R+q1z5QG64Zh7v2od6Blwq2EBkjDkjfc7gV9hV+oRxCA+HV6/CI19DLeaYF6llU3qTdr6GepfrM6taw/gp34V+BI7OnusodlCv5Zot8Bze1zp6o3wu1vTQIbK2Ph3SMUbpP2e4O87x6kz8f6ZjRLaoH1qw5OmTrNULKSsUejZ3ZrqS1QHe4zNHjmQXx/xmiPG7mk3VZKQneV3Gi87PONvyfCpOQPl2yxPTx0yd1ee5HrBy3vK4B3xaFVG/4JF8xqrrD3nN6gqNgK65Oyyucc6neVtFNOobHViNTjjGvwyX318x43s44VuBXTaXl66Q4wtVn/0ezRvNo1r9Ze3kw3HlzFDfqm16Qs/oObRO6Yze0zk1sOMv9I2+0w9tVzvSXmtvCtfNDcV5TFOX9u4P+/N8Ig==</latexit>  \nwe applied in the course of the processing to accept a solution in Section 4. Although in this sample, we have numbered asteroids  \nFigure 14 compares the relative uncertainty between the Gaia DR3 and this FPR solution. It clearly shows the gain resulting from a longer and better-sampled arc for the orbit computation. The uncertainty in Gaia DR3 was typically ≃ 50 times larger, and the distribution was bimodal, with a second population. The second smaller peak in the histogram stems from sources with the shortest arcs (around typically 100 days). This second peak has now disappeared in the FPR as a direct consequence of a fuller coverage of the trajectories. More precisely, the first peak in the Gaia DR3 computations was at ∼ 2 . 6 × 10 -8 and the second at ∼ 6 . 9 × 10 -6 , while the distribution peaks at ∼ 8 × 10 -10 in  \n<!-- image -->  \nFig. 10. Completeness level of successful orbital solutions as a function of the asteroid number per bin of 1000 (top) and 10,000 (bottom). Top: First 50 000 numbered asteroids, with 1000 successes achievable per bin at most. Bottom: Whole set of ≈ 157 000 solutions per group of 10,000.  \n<!-- image -->  \nup to = 400 000, many of the asteroids with large numbers are not observed by Gaia or are too poorly sampled for us to compute an orbit from the small number of available positions. This feature is shown in Figure 10, with the level of completeness as a function of the rank of the asteroid for bin of 1000 or 10,000 planets. The top plot includes up to number = 50 000 and is virtually 100% complete (99% of the first < 25 000 planets with an orbit, and 97 . 12% of the first 50,000). The bottom plot shows a similar histogram for the full data set, showing the sharp decrease in completeness level for smaller and fainter asteroids. It reaches below 50% (5000 solutions in a bin of 10,000 planets) above number 140 000 and reaches 10% from number 300 000. Not much more improvement is expected in the future releases (Gaia DR4 and Gaia DR5) because the missing asteroids are just not detected and will remain so. However, the selection boundary at number 400 000 will be abandoned, and asteroids with larger IDs will be included in the match algorithm, but they will surely be detected with a lower completeness level.', '5.2. Semi-major axis': 'Of the means contrived to quickly grasp the global quality of a new orbit solution of a large number of asteroids, the semi-major axis remains the most significant and favoured orbital parameter. It has the desirable feature over the angular parameters of being a true geometric quantity that is independent of the reference frame. This is particularly valuable when two solutions are compared for which small di ff erences in inclination or in longitude of node may be unconnected with the orbits, but show up just as consequences of tiny rotations between reference frames. The eccentricity arguably shares this geometric nature, but the third Kepler law, linking the orbit size to the period, endows the semi-  \nFig. 11. Number of elementary astrometric observations as function of the number of the asteroid. The black line shows the average values. Asteroids with large numbers discovered since the year 2000 are generally smaller and fainter. A large fraction cannot be seen by the Gaia telescopes and detectors, or is seen only during some limited and more favourable orbital phases.  \n<!-- image -->  \nmajor axis with a much deeper meaning and has been taken for years as the first-rate parameter to compare solutions or to assess a solution. Above all, however, a di ff erence in semi-major axis between two orbits has its counterpart in the mean motion and will soon show up as a secular drift in longitude, dominating any other di ff erence from the other orbital elements and driving the accuracy of the ephemeris a few years away from the epoch.  \nTherefore, we follow this practice here by emphasising the semi-major axis over the angular quantities, which are less significant and are too sensitive to small di ff erences between reference frames. The plots always take the dimensionless quantity σ a / a to express the relative precision, where σ a is the formal standard deviation of the semi-major axis computed by transforming the covariance matrix of the state vector to the matrix for the osculating elements. Similarly, di ff erences between two solutions are shown with the relative distance ∆ a / a or its unsigned value ≈ amplitude).  \nThe two plots in Figure 12 show the relative formal uncertainty as a function of the semi-major axis for the best solutions (top) and for the whole catalogue (bottom). The overall appearance is as expected based on the observational span of 66 months: the Trojan orbits cannot be retrieved with a quality comparable to that of the bright NEOs even after scaling by the semi-major axis. For the first group, the relative uncertainty ranges from 5 × 10 -11 to 2 × 10 -9 for the MBA, and the density peaks about 2 × 10 -10 .  \nThe actual distribution per bin of accuracy is shown in Figure 13 for the same selections. The uncertainties are given in log-scale. The centre of the distribution at ≈ 1 . 5 × 10 -10 for the bright population with numbers < 20 000 is better visible, as is the lower accuracy for the fainter group with a centre at ≈ 5 × 10 -10 in σ a / a . The formal uncertainty is directly derived from the fit, and the covariance matrix is scaled on the post-fit residuals.  \nFig. 12. Formal uncertainty of the orbits measured by σ a / a , where a stands for the semi-major axis of the osculating orbit, and σ a shows its standard deviation. Top: First 50 000 numbered asteroids. Bottom: Complete solution of 157 000 orbits. The vertical line at a = 3 . 15 au corresponds to a period of 66 months, which is the longest arc in the data. The vertical scale is adapted to each set.  \n<!-- image -->  \nFig. 13. Relative formal uncertainty of the Gaia orbits resulting from the fit to the observations over 66 months, divided between the bright end (numbers < 20 000) and the fainter asteroids. The bins add up to 100% for each group.  \n<!-- image -->  \nthis study. This means a gain between one to two orders of magnitude in the uncertainty, which is much larger than a ≃ 1 / √ 2 improvement that would be obtained just by doubling the number of data points used in the adjustment. It essentially results from the longer observational arc, in agreement with Desmars et al. (2013). A final plot for this section shown in Figure 15 gives the formal uncertainty of the semi-major axis as a function of the arc length measured in orbital periods. The improvement as soon as the observation arc covers a full orbit is clearly visible in the steep slope on the left in the diagram, compared to the near stationary or much slower decrease when the arc is ⪆ 1 . 2 orbital period. In this regime, the improvement will come mainly from the combination of the larger number of observations and the extended arc length, while in the small-arc regime,  \nFig. 14. Histogram of the relative uncertainty of the semi-major axis, normalised to maximum, for all asteroids in the middle main belt. Gaia DR3 is shown in green, and this work for the FPR is shown in red.  \n<!-- image -->  \nFig. 15. Relative formal uncertainty of the Gaia orbits as a function of the arc-length coverage expressed in orbital periods. The black dots show this solution, and the light blue dots show the Gaia DR3 solution over 34 months.  \n<!-- image -->  \nthe photon noise is not the main source of (in)accuracy. At mission completion, with 10 . 6 years of data, even the Trojans will have nearly a full orbit coverage, and their orbit uncertainty will decrease to the 1 × 10 -9 relative uncertainty. For the purpose of comparison, the light blue dots in Figure 15 are plotted from the Gaia DR3 orbits, when only 34 months of data were used in the orbital solution, and almost no orbit had a complete orbit coverage. The improvement by a factor about 50 is outstanding in this combined plot. The main-belt orbits are not even as good as the Trojan orbits in the 66-month solution.', '5.3. Comparisons to other solutions': 'As mentioned earlier, several other sources of orbital elements are built on di ff erent observational material, di ff erent time spans, and di ff erent dynamical models. Although we compared the solution from Gaia to the Astorb (Moskovitz et al. 2022), MPC, and JPL (Park et al. 2021) solutions, only the comparisons to JPL orbits are discussed in detail in this paper, and MPC is mentioned only in passing.', '5.3.1. Comparison to JPL orbits': 'The JPL service named Small-Body Data Base Query 5 allowed us to obtain orbital elements without a truncation in the output for a sample with a reasonable size. This was done for the first 50 000 numbered asteroids, which is enough to have a good sample covering a wide range of semi-major axes. The elements are provided at the epoch JD 2460000.5 TDB (25 February 2023) with all the digits, including a safety margin beyond the true precision. The epoch is o ff set by about 5 years from the Gaia mean epoch in our solution, however. A propagation was therefore needed of one of the files, or for both files midway. We chose to propagate the JPL reference data to the Gaia epoch, that is, backward 5 years. This propagation itself is the source of an additional di ff erence between the orbits becausee it is hard to certify that an accuracy better than ≈ 1 × 10 -10 is maintained over 5 years of numerical integration, and above all, that the Gaia and JPL dynamical models are close enough to reach this performance. This propagation for Gaia was made outside the main pipeline with an independent code. The dynamical model includes the complete form of the EIH equations and gravitational perturbations from 14 asteroids.  \nTo ascertain that no significant deviation arose during our numerical integration, we tested this propagation on a much smaller set of asteroids with JPL orbits requested at the Gaia epoch. We were able to reproduce JPL integration at a level of 8 × 10 -10 in the di ff erences ∆ a / a after 5 years. This is acceptable for the comparison, although we know that many of the Gaia orbits are better than this limit (at least in formal uncertainty). It is likely that these residual deviations mainly hold for the small di ff erences between the two dynamical models (not for the same perturbing asteroids, the same masses, or the same Solar System ephemeris, etc.) and not from the numerical integration proper.  \nFig. 16. Relative uncertainty of the reference orbits computed by the JPL before propagation to the mean Gaia epoch. The bins add up to 100% for each group.  \n<!-- image -->  \nWe now describe how Gaia used TCB for its time tagging, as explained in Section 4.2. For a meaningful comparison with JPL or MPC, the same units must be used. Starting with the Gaia TCB-compatible state vector, a standard TCB-to-TDB scaling is performed as  \na TDB = (1 -LB ) a TCB X TDB = (1 -LB ) X TCB V TDB = V TCB , (10)  \nas discussed in Klioner et al. (2010) or Klioner (2008). Here, LB = 1 . 550519768 × 10 -8 is a constant defined in the system of  \nastronomical units and is equal to the rate ( dTCB / dTDB -1). The scaling must not be applied again to the semi-major axis that was computed from the rescaled state vector. The e ff ect is already included in the spatial components of the state vector.  \nFig. 17. Di ff erences in semi-major axis (as ∆ a / a ) between the Gaia orbital solution and the JPL orbits transformed to the Gaia epoch for the first 20 000 asteroids. There are 19 822 Gaia solutions, 261 of which are outside the plot boundaries.  \n<!-- image -->  \nFig. 18. Di ff erences in semi-major axis (as ∆ a / a ) between the Gaia orbital solution and the JPL orbits for the first 50 000 asteroids. The bias increases slowly from 1 . 3 × 10 -10 to 9 . 5 × 10 -10 from the left to the right.  \n<!-- image -->  \nFig. 19. Normalised di ff erences in semi-major axis between the Gaia orbital solution and the JPL orbits transformed to the Gaia epoch for the first 20 000 asteroids. Compared to Figure 17 , 13 bodies lie outside the[ -10 , 10] interval.  \n<!-- image -->  \nWhen the two sets of osculating elements are brought to a common epoch, the di ff erences in semi-major axes as shown in Figure 17 are obtained in the form of the dimensionless quantity ∆ a / a , with a few summary statistics. They are shown in Figure 18 with a scatter plot for the first 50,000 numbered plan-  \nFig. 20. Absolute values of the di ff erences in semi-major axis between the Gaia orbital solution and the JPL orbits transformed to the Gaia epoch for the first 20 000 asteroids. Figure 10 shows 19 , 822 Gaia solutions with ID ≤ 20 000. The di ff erence of 19 , 805 comes from the largest outliers of Figure 21.  \n<!-- image -->  \ns. The agreement is outstanding, even though a slight bias of 5 . 5 × 10 -10 remains between the two solutions for the first 20,000 numbered planets. Expressed in angle, when ∆ a is taken as a positional uncertainty, which is close to the osculating epoch before the drift in mean anomaly takes over, this is equivalent to 0.1 mas. This is a very small deviation given the total independence in the construction of the two solutions. The di ff erences in the dynamical model, in the masses, and above all, in the set of observations (none in common) drive the fit over two very distinct time spans. For a first full Gaia solution, this exceeds the most optimistic expectations. The bias is even smaller at 1 . 3 × 10 -10 for the first 5000 planets and grows to 8 . 5 × 10 -10 when the full set of 50 , 000 is considered. This slight increase in asteroid number is hardly noticeable in the scatter in Figure 18. The scatter about the mean in this plot is otherwise remarkably uniform for the range of selected planets. Its value at ∼ 5 × 10 -9 agrees with the formal uncertainty of either catalogue, as shown in Figs. 13 -16 . It shows that the true accuracy of the solutions is not much di ff erent from the formal accuracy. We are probably very close to the limit of what can be learnt from this comparison, which is as interesting for JPL as it is for Gaia.  \nThe normalised di ff erences computed as ∆ a / q σ ( a ) 2 Gaia + σ ( a ) 2 JPL are shown in Figure 19 for the same set of asteroids. No change was applied to JPL-reported uncertainties to account for the propagation from the JPL to the Gaia epoch. The bias, albeit normalised, has the same origin as in Figure 17 and is statistically significant. The scatter at ∼ 1 . 45 is not a robust estimate and takes too much weight from the tails. The inter-quartile range (IQR) for a normal distribution (IQR = 1 . 35 σ ) gives a more robust estimate of the central width at 1 . 607 / 1 . 35 ≈ 1 . 2 and decreases to 1 . 13 with the first 50 , 000 numbered planets and even 1 . 08 with the numbered planets in [25 , 000 , 50 , 000]. This is a very satisfactory result and indicates that the JPL uncertainties are realistic and fully compatible with the scatter of ∆ a between the two solutions. Because the formal errors for Gaia are smaller than those in JPL, the normalised plot is much less sensitive to an inflation factor that should be applied to Gaia uncertainties, unless this factor is as large as 10, in which case, the random uncertainties would be fully incompatible with the post-fit residuals.', '5.3.2. Comparison to MPC orbits': 'The Minor Planet Center at Harvard (MPC) has recently provided access to its orbits with full accuracy, including the uncertainty on each element. A similar comparison was made for the MPC solution for 50 000 orbits, and it confirmed the absence of significant bias ( < 3 × 10 -9 in this case) and a scatter just twice as large as with JPL solution ( ≈ 1 . 1 × 10 -8 ) for ∆ a / a . These di ff erences between the comparison to JPL and MPC might be traced to the fact that JPL used more observations than MPC in their fit, and in particular, they used some of the Gaia DR2 observations. This good complement demonstrates the perfect consistency between three independent fits and establishes that the Gaia solution is currently better than MPC and at least at the level of JPL. Based strictly on formal uncertainties, the Gaia solution appears even better, but we remain cautious about this claim because systematics could be larger and our arc length remains limited. The real test would be to compare the Gaia orbits to accurate observations (e.g. occultations) well outside the time range of the fit. This is a work in progress. The first results presented in Section 5.6 show that this is a work in itself , and much more work remains to be done for this goal to be reached.', '5.4. Outliers from this comparison': 'The histogram in Figure 20 has a vertical linear scale that prevents us from seeing the extended tails when the relative di ff erences between Gaia and JPL orbits are greater than 1 × 10 -7 or 1 × 10 -6 . The orbital solution of Gaia has very few filters during the iterative process or in the final output. If the solution has converged, the orbit is output, regardless of the formal uncertainty. Given the possibility of very small e ff ective arcs for some asteroids, which are also associated with the small number of detections, incorrect orbits, if rare, may be obtained in the Gaia solution. These are more likely than in JPL, which was built upon a much longer time range for the first 50 000 numbered asteroids considered in this comparison. Setting the level of anomalous di ff erences at 1 × 10 -7 for the di ff erence | ∆ a | / a , deviations that exceed this threshold occur in only 110 orbits. This is indeed an astonishingly small number. It is close to 0 . 2% of the set and it is well accounted for in the residual analysis, as we show below. Many of these orbits may have been filtered out from the Gaia solution based on their anomalous formal uncertainty that is shwon in Figure 21.  \nThis plot shows two well separated groups. On the left, about 50 orbits with Gaia solutions that are acceptable and have a statistically significant di ff erence with JPL. This can be viewed as the normal extended tail of the distribution, and no conclusion can be drawn without considering all solutions individually. The right part is more interesting and shows a strict relation between the departure between the two orbits (Gaia and JPL) and the increasingly lower quality of the Gaia solutions. These ≈ 60 orbits could be removed for the published set without damage, but it is very satisfactory to see this clear feature for such a small number of orbits. All the others except for one agree better with the JPL solution than 1 . 7 × 10 -6 in relative di ff erence in semi-major axis, and they agree usually much better than this limit, as shown in Figure 20. Essentially, there are no unexplained anomalous orbits in the set we used for the comparison. Gaia and JPL can both be proud of this.  \nFig. 21. Anomalous di ff erences in semi-major axis between the Gaia orbital solution and the JPL orbits. The plot shows the relative di ff erence as a function of the Gaia formal uncertainty σ a / a for all solutions with | ∆ a / a | > 1 × 10 -7 . The left group ( σ a / a < 5 × 10 -8 ) is just the small statistical right tail of the distribution of Figure 20 and is determined by the selection threshold for the plot. The group on the right with σ a / a > 1 × 10 -7 comprises solutions with a poor fit performance and is not an issue: these orbits could have been filtered out in the pipeline and were excluded from the delivery.  \n<!-- image -->  \nFig. 22. Di ff erences in eccentricity (expressed in arcsec) and inclination between the Gaia and JPL solutions for the first 10 000 numbered asteroids.  \n<!-- image -->', '5.5. Other orbital elements': 'The comparison of the other orbital elements (eccentricity and the angular elements) is briefly mentioned now with additional plots. The main point is to draw attention to the details of the reference ecliptic used by either group. Any di ff erence in the definition of the ecliptic (both the plane and the origin of longitude) would result in spurious systematic deviations that are unrelated to a real di ff erence between the orbits. As said earlier, osculating elements of Gaia have been compared to the same ecliptic  \nFig. 23. Di ff erences in inclination between the Gaia and JPL solutions. The core is bounded with ± 10 mas, but there are extended non-Gaussian tails.  \n<!-- image -->  \nFig. 24. Di ff erences in longitude of node (as ∆Ω sin i in arcsec) between the Gaia and JPL solutions for the first 50,000 numbered minor planets.  \n<!-- image -->  \nplane as JPL, with origin of longitude at the node with the ICRF equator as defined in Figure 9.  \nWe show in Figure 22 histograms of the di ff erence in orbital eccentricity and inclination found between the Gaia and JPL solutions for the first 10 000 numbered asteroids. Angular units are used, with dimensionless ∆ e taken as radians and converted into arcseconds (as if a ∆ e were a position shift and ( a ∆ e ) / a = ∆ e were an angular shift). The two solutions are not biased, and the scatter about the mean is typically 5 to 10 mas or < 5 × 10 -8 , with a robust Gaussian scatter of 40 and 43 mas, respectively). This is larger than the equivalent angular deviation seen with the semi-major axis, but it probably gives a good upper bound for the true quality of the orbits. The prominent tails in the inclination are not Gaussian, however.  \nThe following comparisons in Figure 23, Figure 24, and Figure 25 show scatter plots of the di ff erences between Gaia and JPL for the three angular elements inclination, longitude of node, and argument of perihelion. The last two are scaled with the inclination and the eccentricity to show the equivalent of a displacement on the sky rather than the di ff erence in the less meaningful coordinates. The three plots are very similar. They have a core within 5 mas that is surrounded by an extended tail that increases with the rank of the asteroid. An additional analysis shows that the variation in inclination has a systematic sine wave of ∼ 10 mas amplitude with the longitude of node. This is a clear sign of a reference-frame e ff ect. However, this e ff ect is hardly visible for asteroid numbers < 10 , 000 and becomes more visible with larger numbers.  \nFig. 25. Di ff erences in argument of perihelion (as e ∆ ω in arcsec) between the Gaia and JPL solutions for the first 50,000 numbered minor planets.  \n<!-- image -->', '5.6. Comparison to stellar occultations': 'Afinal severe test of the quality of orbits published in the FPR is provided by stellar occultations. In principle, positions computed by using FPR orbits at the epoch of an observed occultation of a star present in the Gaia catalogue must match the astrometry derived from the event.  \nWe set up our test for the objects with an orbit in the FPR and that were successfully observed by occultation. The selected data set resulted in 978 asteroids, associated with 5,774 astrometric measurements by stellar occultations. Large and bright asteroids with better orbits in the MPC or JPL data bases dominate the sample because their occultations are easier to predict. In addition, because they are larger than the faint ones, their occultations last longer and are easier to catch.  \nWe propagated the FPR orbits to the epoch of each observed occultation and computed the corresponding position and its difference (O -C) in the directions of the equatorial frame with respect to the observed astrometry. We limit our study to occultations by main-belt asteroids, which are the largest majority of the sample for occultation astrometry and Gaia orbits. In the following, we discuss the total O -C, that is, the distance on the tangent sky plane between the observed and the computed position.  \nThe global distribution of the O -C values (Fig. 26) reveals a large spread that reaches several hundred milliarcseconds. Most observations lie below ≈ 50 mas. This corresponds to the order of magnitude of the expected uncertainty for the dominating fraction of occultations with a small number of observers. In this situation, the cross-track astrometric error is dominated by the apparent size of the asteroid (Ferreira et al. 2022).  \nHowever, when only occultation events contemporary to the Gaia mission (same figure) are selected, a clear subset of much smaller O -C appears. It peaks at about 10 mas and has a tail of high values that is substantially suppressed.  \nWe can then verify whether an equivalent improvement can be obtained by selecting only occultations that are related to objects with very high-quality orbits. We set a threshold on the maximum allowed σ a in such a way that the same number of occultation events as for the 2014-2020 period was selected (1819). The obtained curve shows that this selection is not su ffi cient to eliminate the high residuals, and the peak of the distribution does not move appreciably with respect to the overall distribution.  \nThe apparent asteroid size (mentioned above as the main source of uncertainty) and other factors contributing to errors of occultation astrometry were factored into the total error budget following the guidelines illustrated by Herald et al. (2020). A result that is more independent of the details of this error model  \n<!-- image -->  \nFig. 26. Histogram of the absolute total O -C (distance on the sky) between the predicted and observed positions derived from stellar occultations. The bottom panel presents the O -C values normalised by the uncertainty of the occultation astrometry for each event. The curve labelled "2014-2020" includes only occultations around those years, while "best orbit" considers only orbits with σ a < 1.6 × 10 -10 au.  \n<!-- image -->  \ncan be obtained by dividing the O -C values by the uncertainty associated with each position (Fig. 26, bottom panel). The obtained distribution is self-consistent and peaks around unity. The selection of the Gaia period clearly stands out as the best sample.  \nThis evidence highlights a trend with time that is revealed by Fig. 27. Uncertainties in the occultation astrometry in the 1990s may be a ff ected by worse statistics, but later observations are rather homogeneous (Fig. 2 in Ferreira et al. 2022). The minimum around the period of nominal operations of Gaia that was used to obtain the FPR orbits is thus not due to the properties of the errors in asteroid occultations, but is clearly related to the time interval covered by the Gaia data. Interestingly, this result allows us to quantify a degradation in the ephemeris of ∼ 10 mas / year. This clearly only holds for the asteroid sample that is observed via occultations, which is biased towards bright asteroids. Gaia astrometry is excellent but not fully optimal for these sources.  \nAt this stage, we also investigated this time trend for occultations involving objects with only the best orbits (same selection as in Fig. 26). We found no significant di ff erence. This probably  \nFig. 27. Absolute values of O -C for the occultation data set averaged on yearly bases. The error bars correspond to the standard deviation.  \n<!-- image -->  \nmeans that the degradation in our sample of orbits is more sensitive to the limited observing arc than to the internal di ff erences in quality.', '6. Conclusion': 'We have reported a new orbital solution for nearly 160 , 000 asteroids built on the Gaia astrometric measurements collected over 66 months of operation and showed its main properties in terms of accuracy and reliability. Comparisons with the best existing orbits published by the JPL or the MPC have shown the impressive consistency between these solutions, which arises from di ff erent observation sets and time coverage and dynamical modelling. We established the absence of systematic di ff erences larger than 5 × 10 -10 in relative di ff erence in the semi-major axes.  \nThis Gaia FPR orbital solution is the first within the Gaia releases whose time coverage is at least equal to the orbital period of main-belt asteroids. This feature is needed to ensure the accuracy of the orbit solution. These expectations were mostly shown to be correct in this paper, with an impressive overall improvement compared to Gaia DR3 solution by a factor almost 50 (see Figure 15), whereas the number of observations was just twice as large.  \nAnimportant issue that arose in the course of this work is that institutes that publish orbital elements should agree on the reference ecliptic (same plane and same origin of longitude), while di ff erent realisations are available that di ff er from each other by several tens milliarcseconds. With the current quality reached at JPL, MPC, NEODyS, and Gaia, di ff erences as small as this matter. Tiny deviations of few 1 × 10 -8 in dimensionless parameters or a few milliarcseconds in the angular elements on otherwise perfectly agreeing orbits are easily created with small di ff erences in reference frames. With sub-milliarcsecond astrometry and orbits approaching 1 × 10 -10 in relative accuracy from JPL or Gaia, the community needs to coordinate, possibly with IAU or IERS, in order to agree about the definition and realisation of the ecliptic, physical or conventional, and about its relation to ICRF and J2000. This goes much beyond the orbits of small bodies because it must be consistent with the precession-nutation model the user will apply to obtain the positions in the equator of date.  \nThe next step will be the future Gaia data releases, starting with Gaia DR4. They will complete the current set in two  \nways: (i) the photometric and spectrophotometric data, and data for comets and planetary satellites not analysed here; (ii) a much larger set of asteroids, about twice as large in number, because this has been limited in the FPR to the same set as selected for Gaia DR3, solely because the astrometric solution for the remainder was not yet available.  \nAfew years from today, the Gaia mission will finally provide data over approximately ten years, which will double the time interval of the FPR data we used here, and which will improve the quality, precision, and accuracy of the orbit determination of small Solar System bodies further. With this trend, Trojan orbits are expected to be as good as those of the main-belt asteroids are today. As of June 2023, about 85% of the expected data are already safely archived on the ground for the upcoming processing.  \nAcknowledgements. We are grateful to an anonymous referee of their kind remarks and constructive comments on several key-points of the paper. Their review has helped significantly improve the organisation and clarity of the paper. This work presents results from the European Space Agency (ESA) space mission Gaia. Gaia data are being processed by the Gaia Data Processing and Analysis Consortium (DPAC). Funding for the DPAC is provided by national institutions, in particular the institutions participating in the Gaia MultiLateral Agreement (MLA). The Gaia mission website is https://www.cosmos.esa.int/ gaia . The Gaia archive website is https://archives.esac.esa.int/gaia .', 'References': 'Astropy Collaboration, Price-Whelan, A. M., Sip"ocz, B. M., et al. 2018, AJ, 156, 123  \n- Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33\n- Beutler, G. 2005, Methods of Celestial Mechanics I: Physical, Mathematical, and\n- Numerical Principles (Germany, Berlin Heidelberg: Springer-Verlag) Carpino, M., Milani, A., & Chesley, S. R. 2003, Icarus, 166, 248\n- Carry, B., Thuillot, W., Spoto, F., et al. 2021, Astronomy & Astrophysics, 648, A96\n- Chapront, J., Chapront-Touzé, M., & Francou, G. 2002, A&A, 387, 700 Deram, P., Fienga, A., Verma, A. K., Gastineau, M., & Laskar, J. 2022, arXiv e-prints, arXiv:2203.01586\n- Desmars, J., Bancelin, D., Hestro ff er, D., & Thuillot, W. 2013, A&A, 554, A32 Ferreira, J. F., Tanga, P., Spoto, F., Machado, P., & Herald, D. 2022, A&A, 658, A73\n- Fienga, A., Deram, P., Viswanathan, V., et al. 2019, Notes scientifiques et techniques de l\'IMCCE\n- Fienga, A., Manche, H., Laskar, J., Gastineau, M., & Verma, A. 2013, arXiv e-prints, arXiv:1301.1510\n- Fukushima, T. 2003, AJ, 126, 494\n- Gaia Collaboration, Brown, A., Vallenari, A., et al. 2016a, A&A, 595, A2  \nGaia Collaboration, Prusti, T., de Bruijne, J., et al. 2016b, A&A, 595, A1  \n- Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1 Herald, D., Gault, D., Anderson, R., et al. 2020, Monthly Notices of the Royal Astronomical Society, 499, 4570\n- Hilton, J. L., Capitaine, N., Chapront, J., et al. 2006, Celestial Mechanics and Dynamical Astronomy, 94, 351\n- Hilton, J. L. & Hohenkerk, C. Y. 2004, A&A, 413, 765\n- Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90\n- Klioner, S. A. 2003, AJ, 125, 1580\n- Klioner, S. A. 2008, A&A, 478, 951\n- Klioner, S. A., Capitaine, N., Folkner, W. M., et al. 2010, Relativity in Fundamental Astronomy: Dynamics, Reference Frames, and Data Analysis, Proceedings, IAU Symposium, 261, 79\n- Levinger, B. W. 1980, Mathematics Magazine, 53, 222  \nMilani, A. & Gronchi, G. 2010, Theory of Orbit Determination (United Kingdom, Cambridge: Cambridge University Press)  \n- Moskovitz, N. A., Wasserman, L., Burt, B., et al. 2022, Astronomy and Computing, 41, 100661  \nMuinonen, K., Berthier, J., Cellino, A., et al. 2022, Gaia DR3 documentation Chapter 8: Solar System Objects, Gaia DR3 documentation, European Space Agency; Gaia Data Processing and Analysis Consortium. Online at \'https: // gea.esac.esa.int / archive / documentation / GDR3 / index.html\'  \nPark, R. S., Folkner, W. M., Williams, J. G., & Boggs, D. H. 2021, AJ, 161, 105 Pontriaguine, L. 1969, Equations Di ff érentielles Ordinaires (USSR, Moscow: Editions MIR)  \n- Tanga, P., Pauwels, T., Mignard, F., et al. 2022, A&A, in publication, doi:10.1051 / 0004-6361 / 202243796 https: // doi.org / 10.1051 / 00046361 / 202243796\n- Tanga, P., Pauwels, T., Mignard, F., et al. 2023, A&A, 674, A12\n- Taylor, M. B. 2005, Astronomical Data Analysis XIV proceedings, 347, 29', 'Access FPR asteroid data in the Gaia archive': "The data used in this article can be found at the Gaia mission archive of ESA 6 . By navigating to the '(Advanced) ADQL tab', the tree of the data base can be expanded. Asteroid data are in the focused\\_product\\_release branch, in tables named sso\\_source (orbits) and sso\\_observation (epoch astrometry).  \nThe meaning of the fields for the astrometry is extensively explained in the Gaia DR3 documentation and the related article (Tanga et al. 2022). The orbital elements have their equivalent in the gaiadr3.sso\\_orbits table, with all the relevant details provided in this article. The residuals used to derive data statistics can be computed starting from the published orbits.  \nAll the Solar System objects present in the FPR were also part of Gaia DR3. However, as the FPR processing is completely independent of the previous processing, the inverse is not strictly true and a small number of Gaia DR3 sources may be missing in the FPR. These di ff erences are very marginal, however.  \nWe provide below an example of an ADQL query that can be used to retrieve the state vector and semi-major axis uncertainty for all the valid orbital solutions, also excluding objects for which the heliocentric orbit cannot be computed (natural planetary satellites):  \n```\nSELECT number\\_mp, denomination, osc\\_epoch, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,1) AS x, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,2) AS y, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,3) AS z, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,4) AS vx, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,5) AS vy, ESDC\\_ARRAY\\_ELEMENT(h\\_state\\_vector,6) AS vz, SQRT(ESDC\\_ARRAY\\_ELEMENT (orbital\\_elements\\_var\\_covar\\_matrix, 1))\n```  \nAS sigma\\_a  \nFROM gaiafpr.sso\\_source WHERE ESDC\\_ARRAY\\_ELEMENT (orbital\\_elements\\_var\\_covar\\_matrix,1)>0. ORDER BY number\\_mp ASC  \nThe output is sorted by asteroid number and contains 156,762 orbital solutions. The content of the whole covariance matrix can be added by a similar method. The filter for the first element of the covariance matrix is sometimes negative when the orbit solution is poor.  \nAll the astrometric measurements for a given object can be retrieved based on their identifiers. For instance, for (704) Interamnia,  \nSELECT * FROM gaiafpr.sso\\_observation WHERE number\\_mp=704.  \nTo retrieve the astrometry for all asteroids associated with the best valid orbits with σ a < 1 × 10 -10 au, we can write  \n```\nSELECT astrom.* FROM gaiafpr.sso\\_observation AS astrom INNER JOIN gaiafpr.sso\\_source AS so USING(number\\_mp) WHERE ESDC\\_ARRAY\\_ELEMENT (so.orbital\\_elements\\_var\\_covar\\_matrix,1)>0. AND SQRT(ESDC\\_ARRAY\\_ELEMENT (so.orbital\\_elements\\_var\\_covar\\_matrix,1))<1.e-10 AND astrom.is\\_rejected = 'false'\n```  \nFinally, in the following example, we compute the time span covered by the observations of each asteroid,  \nSELECT number\\_mp, MAX(epoch)-MIN(epoch) AS t\\_delta\\_days FROM (SELECT astrom.number\\_mp, astrom.epoch FROM gaiafpr.sso\\_observation AS astrom INNER JOIN gaiafpr.sso\\_source AS so USING(number\\_mp)) AS subquery WHERE(number\\_mp>0) GROUP BY number\\_mp.", 'Appendix B:': 'Transformation of the state vector into ecliptic  \nThe state vector resulting from the orbit fitting is a 6D vector giving the initial condition at the epoch (one per body) as a combination of the position vector and the velocity vector into a single state vector. The origin is heliocentric, and the axes are aligned to the ICRF axes. The units are au and au / day, and the scaling of Equation 10 must be applied to obtain the TDB compatible values. By reference to Figure 9 and Table 1, the coordinates of the same vector in the heliocentric ecliptic frame are given by  \nr ecl = R z ( ψ ) R x ( ϵ ) R z ( -ϕ ) r ICRF (B.1) v ecl = R z ( ψ ) R x ( ϵ ) R z ( -ϕ ) v ICRF, (B.2)  \nwhere R x , y , z ( α ) denotes the passive rotation of angle α about the axes x , y , or z . In matrix form, this is  \nR x ( α ) =         1 0 0 0 cos α sin α -         \nR z ( α ) =         cos α sin α 0 -sin α cos α 0         .  \n0 sin α cos α 0 0 1  \nWith the state vector given in the ecliptic frame (more precisely, in one of the possible choices of the ecliptic), the transformation into osculating elements is common to all the cases and is normally available with reliable routines in any group handling orbits to perform the transformation in both directions for elliptical or hyperbolic orbits. However, the solar mass parameter must also be included in in this transformation to ensure that it agrees with the mass parameter that is used in the force model.  \nThe user who needs to compute a good precision (e.g. lower than 10 mas) position in the mean equator of the date must be aware that some precession routines may automatically include some of the rotations above in the form of a frame bias to relate the ICRF frame to the J2000 equatorial frame.', 'Appendix C:': "The Gaia mission and data processing have financially been supported by (in alphabetical order by country)  \n- -the Algerian Centre de Recherche en Astronomie, Astrophysique et Géophysique of Bouzareah Observatory;\n- -the Austrian Fonds zur Förderung der wissenschaftlichen Forschung (FWF) Hertha Firnberg Programme through grants T359, P20046, and P23737;\n- -the BELgian federal Science Policy O ffi ce (BELSPO) through various PROgramme de Développement d'Expériences scientifiques (PRODEX) grants and the Polish Academy of Sciences - Fonds Wetenschappelijk Onderzoek through grant VS.091.16N, and the Fonds de la  \nRecherche Scientifique (FNRS), and the Research Council of Katholieke Universiteit (KU) Leuven through grant C16 / 18 / 005 (Pushing AsteRoseismology to the next level with TESS, GaiA, and the Sloan DIgital Sky SurvEy PARADISE);  \n- -the Brazil-France exchange programmes Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) and Coordenação de Aperfeicoamento de Pessoal de Nível Superior (CAPES) - Comité Français d'Evaluation de la Coopération Universitaire et Scientifique avec le Brésil (COFECUB);\n- -the Chilean Agencia Nacional de Investigación y Desarrollo (ANID) through Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT) Regular Project 1210992 (L. Chemin);\n- -the National Natural Science Foundation of China (NSFC) through grants 11573054, 11703065, and 12173069, the China Scholarship Council through grant 201806040200, and the Natural Science Foundation of Shanghai through grant 21ZR1474100;\n- -the Tenure Track Pilot Programme of the Croatian Science Foundation and the École Polytechnique Fédérale de Lausanne and the project TTP-2018-07-1171 'Mining the Variable Sky', with the funds of the Croatian-Swiss Research Programme;\n- -the Czech-Republic Ministry of Education, Youth, and Sports through grant LG 15010 and INTER-EXCELLENCE grant LTAUSA18093, and the Czech Space O ffi ce through ESA PECS contract 98058;\n- -the Danish Ministry of Science;\n- -the Estonian Ministry of Education and Research through grant IUT40-1;\n- -the European Commission's Sixth Framework Programme through the European Leadership in Space Astrometry (ELSA) Marie Curie Research Training Network (MRTNCT-2006-033481), through Marie Curie project PIOFGA-2009-255267 (Space AsteroSeismology & RR Lyrae stars, SAS-RRL), and through a Marie Curie Transferof-Knowledge (ToK) fellowship (MTKD-CT-2004-014188); the European Commission's Seventh Framework Programme through grant FP7-606740 (FP7-SPACE-2013-1) for the Gaia European Network for Improved data User Services (GENIUS) and through grant 264895 for the Gaia Research for European Astronomy Training (GREAT-ITN) network;\n- -the European Cooperation in Science and Technology (COST) through COST Action CA18104 'Revealing the Milky Way with Gaia (MW-Gaia)';\n- -the European Research Council (ERC) through grants 320360, 647208, and 834148 and through the European Union's Horizon 2020 research and innovation and excellent science programmes through Marie Skłodowska-Curie grant 745617 (Our Galaxy at full HD - Gal-HD) and 895174 (The build-up and fate of self-gravitating systems in the Universe) as well as grants 687378 (Small Bodies: Near and Far), 682115 (Using the Magellanic Clouds to Understand the Interaction of Galaxies), 695099 (A sub-percent distance scale from binaries and Cepheids - CepBin), 716155 (Structured ACCREtion Disks - SACCRED), 951549 (Subpercent calibration of the extragalactic distance scale in the era of big surveys - UniverScale), and 101004214 (Innovative Scientific Data Exploration and Exploitation Applications for Space Sciences - EXPLORE);\n- -the European Science Foundation (ESF), in the framework of the Gaia Research for European Astronomy Training Research Network Programme (GREAT-ESF);\n- -the European Space Agency (ESA) in the framework of the Gaia project, through the Plan for European Cooperating States (PECS) programme through contracts C98090 and 4000106398 / 12 / NL / KML for Hungary, through contract 4000115263 / 15 / NL / IB for Germany, and through PROgramme de Développement d'Expériences scientifiques (PRODEX) grant 4000127986 for Slovenia;\n- -the Academy of Finland through grants 299543, 307157, 325805, 328654, 336546, and 345115 and the Magnus Ehrnrooth Foundation;\n- -the French Centre National d'Études Spatiales (CNES), the Agence Nationale de la Recherche (ANR) through grant ANR-10-IDEX-0001-02 for the 'Investissements d'avenir' programme, through grant ANR-15-CE31-0007 for project 'Modelling the Milky Way in the Gaia era' (MOD4Gaia), through grant ANR-14-CE33-0014-01 for project 'The Milky Way disc formation in the Gaia era' (ARCHEOGAL), through grant ANR-15-CE31-0012-01 for project 'Unlocking the potential of Cepheids as primary distance calibrators' (UnlockCepheids), through grant ANR-19-CE310017 for project 'Secular evolution of galaxies' (SEGAL), and through grant ANR-18-CE31-0006 for project 'Galactic Dark Matter' (GaDaMa), the Centre National de la Recherche Scientifique (CNRS) and its SNO Gaia of the Institut des Sciences de l'Univers (INSU), its Programmes Nationaux: Cosmologie et Galaxies (PNCG), Gravitation Références Astronomie Métrologie (PNGRAM), Planétologie (PNP), Physique et Chimie du Milieu Interstellaire (PCMI), and Physique Stellaire (PNPS), the 'Action Fédératrice Gaia' of the Observatoire de Paris, the Région de Franche-Comté, the Institut National Polytechnique (INP) and the Institut National de Physique nucléaire et de Physique des Particules (IN2P3) co-funded by CNES;\n- -the German Aerospace Agency (Deutsches Zentrum für Luft- und Raumfahrt e.V., DLR) through grants 50QG0501, 50QG0601, 50QG0602, 50QG0701, 50QG0901, 50QG1001, 50QG1101, 50QG1401, 50QG1402, 50QG1403, 50QG1404, 50QG1904, 50QG2101, 50QG2102, and 50QG2202, and the Centre for Information Services and High Performance Computing (ZIH) at the Technische Universität Dresden for generous allocations of computer time;\n- -the Hungarian Academy of Sciences through the Lendület Programme grants LP2014-17 and LP2018-7 and the Hungarian National Research, Development, and Innovation Office (NKFIH) through grant KKP-137523 ('SeismoLab');\n- -the Science Foundation Ireland (SFI) through a Royal Society - SFI University Research Fellowship (M. Fraser);\n- -the Israel Ministry of Science and Technology through grant 3-18143 and the Tel Aviv University Center for Artificial Intelligence and Data Science (TAD) through a grant;\n- -the Agenzia Spaziale Italiana (ASI) through contracts I / 037 / 08 / 0, I / 058 / 10 / 0, 2014-025-R.0, 2014-025-R.1.2015, and 2018-24-HH.0 to the Italian Istituto Nazionale di Astrofisica (INAF), contract 2014-049-R.0 / 1 / 2 to INAF for the Space Science Data Centre (SSDC, formerly known as the ASI Science Data Center, ASDC), contracts I / 008 / 10 / 0, 2013 / 030 / I.0, 2013-030-I.0.1-2015, and 2016-17-I.0 to the Aerospace Logistics Technology Engineering Company (ALTEC S.p.A.), INAF, and the Italian Ministry of Education, University, and Research (Ministero dell'Istruzione, dell'Università e della Ricerca) through the Premiale project 'MIning The Cosmos Big Data and Innovative Italian Technology for Frontier Astrophysics and Cosmology' (MITiC);  \n- -the Netherlands Organisation for Scientific Research (NWO) through grant NWO-M-614.061.414, through a VICI grant (A. Helmi), and through a Spinoza prize (A. Helmi), and the Netherlands Research School for Astronomy (NOVA);\n- -the Polish National Science Centre through HARMONIA grant 2018 / 30 / M / ST9 / 00311 and DAINA grant 2017 / 27 / L / ST9 / 03221 and the Ministry of Science and Higher Education (MNiSW) through grant DIR / WK / 2018 / 12;\n- -the Portuguese Fundação para a Ciência e a Tecnologia (FCT) through national funds, grants SFRH / -BD / 128840 / 2017 and PTDC / FIS-AST / 30389 / 2017, and work contract DL 57 / 2016 / CP1364 / CT0006, the Fundo Europeu de Desenvolvimento Regional (FEDER) through grant POCI-01-0145-FEDER-030389 and its Programa Operacional Competitividade e Internacionalização (COMPETE2020) through grants UIDB / 04434 / 2020 and UIDP / 04434 / 2020, and the Strategic Programme UIDB / -00099 / 2020 for the Centro de Astrofísica e Gravitação (CENTRA);\n- -the Slovenian Research Agency through grant P1-0188;\n- -the Spanish Ministry of Economy (MINECO / FEDER, UE), the Spanish Ministry of Science and Innovation (MICIN), the Spanish Ministry of Education, Culture, and Sports, and the Spanish Government through grants BES-2016-078499, BES-2017-083126, BES-C-20170085, ESP2016-80079-C2-1-R, ESP2016-80079-C2-2-R, FPU16 / 03827, PDC2021-121059-C22, RTI2018-095076-BC22, and TIN2015-65316-P ('Computación de Altas Prestaciones VII'), the Juan de la Cierva Incorporación Programme (FJCI-2015-2671 and IJC2019-04862-I for F. Anders), the Severo Ochoa Centre of Excellence Programme (SEV20150493), and MICIN / AEI / 10.13039 / 501100011033 (and the European Union through European Regional Development Fund 'A way of making Europe') through grant RTI2018095076-B-C21, the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia 'María de Maeztu') through grant CEX2019-000918-M, the University of Barcelona's o ffi cial doctoral programme for the development of an R + D + i project through an Ajuts de Personal Investigador en Formació (APIF) grant, the Spanish Virtual Observatory through project AyA201784089, the Galician Regional Government, Xunta de Galicia, through grants ED431B-2021 / 36, ED481A-2019 / 155, and ED481A-2021 / 296, the Centro de Investigación en Tecnologías de la Información y las Comunicaciones (CITIC), funded by the Xunta de Galicia and the European Union (European Regional Development Fund Galicia 2014-2020 Programme), through grant ED431G2019 / 01, the Red Española de Supercomputación (RES) computer resources at MareNostrum, the Barcelona Supercomputing Centre - Centro Nacional de Supercomputación (BSC-CNS) through activities AECT-2017-2-0002, AECT-2017-3-0006, AECT-2018-1-0017, AECT-2018-20013, AECT-2018-3-0011, AECT-2019-1-0010, AECT2019-2-0014, AECT-2019-3-0003, AECT-2020-1-0004, and DATA-2020-1-0010, the Departament d'Innovació, Universitats i Empresa de la Generalitat de Catalunya through grant 2014-SGR-1051 for project 'Models de Programació i Entorns d'Execució Parallels' (MPEXPAR), and Ramon y Cajal Fellowship RYC2018-025968-I funded by MICIN / AEI / 10.13039 / 501100011033 and the European Science Foundation ('Investing in your future');\n- -the Swedish National Space Agency (SNSA / Rymdstyrelsen);  \n- -the Swiss State Secretariat for Education, Research, and Innovation through the Swiss Activités Nationales Complémentaires and the Swiss National Science Foundation through an Eccellenza Professorial Fellowship (award PCEFP2\\_194638 for R. Anderson);\n- -the United Kingdom Particle Physics and Astronomy Research Council (PPARC), the United Kingdom Science and Technology Facilities Council (STFC), and the United Kingdom Space Agency (UKSA) through the following grants to the University of Bristol, the University of Cambridge, the University of Edinburgh, the University of Leicester, the Mullard Space Sciences Laboratory of University College London, and the United Kingdom Rutherford Appleton Laboratory (RAL): PP / D006511 / 1, PP / D006546 / 1, PP / D006570 / 1, ST / I000852 / 1, ST / J005045 / 1, ST / K00056X / 1, ST / K000209 / 1, ST / K000756 / 1, ST / L006561 / 1, ST / N000595 / 1, ST / N000641 / 1, ST / N000978 / 1, ST / N001117 / 1, ST / S000089 / 1, ST / S000976 / 1, ST / S000984 / 1, ST / S001123 / 1, ST / S001948 / 1, ST / -S001980 / 1, ST / S002103 / 1, ST / V000969 / 1, ST / W002469 / 1, ST / W002493 / 1, ST / W002671 / 1, ST / W002809 / 1, and EP / V520342 / 1.  \nThe GBOT programme uses observations collected at (i) the European Organisation for Astronomical Research in the Southern Hemisphere (ESO) with the VLT Survey Telescope (VST), under ESO programmes 092.B-0165, 093.B-0236, 094.B-0181, 095.B-0046, 096.B-0162, 097.B-0304, 098.B-0030, 099.B0034, 0100.B-0131, 0101.B-0156, 0102.B-0174, and 0103.B0165; and (ii) the Liverpool Telescope, which is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the United Kingdom Science and Technology Facilities Council, and (iii) telescopes of the Las Cumbres Observatory Global Telescope Network.  \nL. Liberato acknowledges support by the Coordena cão de Aperfei coamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001, also by CAPES-PRINT Process 88887.570251 / 2020-00.  \nThe authors want to acknowledge Valéry Lainey (IMCCE, Paris obseervatory) for providing the ephemerides of planetary satellites, and Josselin Desmars (IMCCE, Paris observatory) for providing extensive external and quality checks on the orbit computations from the NIMA software.  \nWe made use of the software products TOPCAT, (Taylor 2005); Matplotlib (Hunter 2007); Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013, 2018)."}