Issue |
A&A
Volume 697, May 2025
Euclid on Sky
|
|
---|---|---|
Article Number | A7 | |
Number of page(s) | 15 | |
Section | Galactic structure, stellar clusters and populations | |
DOI | https://doi.org/10.1051/0004-6361/202450793 | |
Published online | 30 April 2025 |
Euclid: Early Release Observations – A glance at free-floating newborn planets in the σ Orionis cluster★
1
Instituto de Astrofísica de Canarias,
Calle Vía Láctea s/n,
38204,
San Cristóbal de La Laguna, Tenerife,
Spain
2
Departamento de Astrofísica, Universidad de La Laguna,
38206
La Laguna, Tenerife,
Spain
3
Laboratoire d’Astrophysique de Bordeaux, CNRS and Université de Bordeaux, Allée Geoffroy St. Hilaire,
33165
Pessac,
France
4
Institut universitaire de France (IUF),
1 rue Descartes,
75231
Paris Cedex 05, France
5
Departamento de Física Fundamental. Universidad de Salamanca.
Plaza de la Merced s/n,
37008
Salamanca,
Spain
6
Centro de Astrobiología (CAB), CSIC-INTA, ESAC Campus,
Camino Bajo del Castillo s/n,
28692
Villanueva de la Cañada, Madrid,
Spain
7
Departamento de Inteligencia Artificial, Universidad Nacional de Educación a Distancia (UNED),
c/Juan del Rosal 16,
28040
Madrid, Spain
8
Departamento Física Aplicada, Universidad Politécnica de Cartagena,
Campus Muralla del Mar,
30202
Cartagena, Murcia,
Spain
9
European Space Agency/ESTEC,
Keplerlaan 1,
2201
AZ Noordwijk, The Netherlands
10
Aurora Technology for European Space Agency (ESA),
Camino bajo del Castillo, s/n, Urbanizacion Villafranca del Castillo,
Villanueva de la Cañada, 28692 Madrid,
Spain
11
Observatoire Astronomique de Strasbourg (ObAS), Université de Strasbourg – CNRS, UMR 7550,
Strasbourg,
France
12
Kapteyn Astronomical Institute, University of Groningen,
PO Box 800,
9700
AV Groningen, The Netherlands
13
Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM,
91191
Gif-sur-Yvette, France
14
Max-Planck-Institut für Astronomie,
Königstuhl 17,
69117
Heidelberg, Germany
15
Ohio University, Physics & Astronomy Department,1 Ohio University,
Athens,
OH
45701, USA
16
NSF’s NOIR,
Lab 950 N. Cherry Avenue,
Tucson,
AZ
85719, USA
17
Department of Physics, International University,
Ho Chi Minh City, Vietnam
18
Vietnam National University,
Ho Chi Minh City, Vietnam
19
International Space University,
1 rue Jean-Dominique Cassini,
67400
Illkirch-Graffenstaden, France
20
School of Physics and Astronomy, University of Leicester,
University Road,
Leicester
LE1 7RH, UK
21
Université de Franche-Comté, Institut UTINAM, CNRS UMR6213, OSU THETA Franche-Comté-Bourgogne, Observatoire de Besançon,
BP 1615,
25010
Besançon Cedex, France
22
INAF-Osservatorio Astrofisico di Torino,
Via Osservatorio 20,
10025
Pino Torinese (TO), Italy
23
School of Physics, Astronomy and Mathematics, University of Hertfordshire,
College Lane,
Hatfield
AL10 9AB, UK
24
Université Paris-Saclay, CNRS, Institut d’astrophysique spatiale,
91405,
Orsay,
France
25
ESAC/ESA, Camino Bajo del Castillo,
s/n, Urb. Villafranca del Castillo,
28692
Villanueva de la Cañada, Madrid,
Spain
26
INAF – Osservatorio Astronomico di Brera,
Via Brera 28,
20122
Milano, Italy
27
INAF – Osservatorio di Astrofisica e Scienza dello Spazio di Bologna,
Via Piero Gobetti 93/3,
40129
Bologna, Italy
28
Dipartimento di Fisica e Astronomia, Università di Bologna,
Via Gobetti 93/2,
40129
Bologna, Italy
29
INFN – Sezione di Bologna,
Viale Berti Pichat 6/2,
40127
Bologna, Italy
30
INAF – Osservatorio Astronomico di Padova,
Via dell’Osservatorio 5,
35122
Padova, Italy
31
Centre National d’Etudes Spatiales – Centre spatial de Toulouse,
18 avenue Edouard Belin,
31401
Toulouse Cedex 9, France
32
Max Planck Institute for Extraterrestrial Physics,
Giessenbachstr. 1,
85748
Garching, Germany
33
Universitäts-Sternwarte München, Fakultät für Physik, Ludwig-Maximilians-Universität München,
Scheinerstrasse 1,
81679
München, Germany
34
Dipartimento di Fisica, Università di Genova,
Via Dodecaneso 33,
16146
Genova, Italy
35
INFN – Sezione di Genova,
Via Dodecaneso 33,
16146
Genova, Italy
36
Department of Physics “E. Pancini", University Federico II,
Via Cinthia 6,
80126
Napoli, Italy
37
INAF – Osservatorio Astronomico di Capodimonte,
Via Moiariello 16,
80131
Napoli, Italy
38
INFN section of Naples,
Via Cinthia 6,
80126
Napoli, Italy
39
Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas,
4150-762
Porto, Portugal
40
Dipartimento di Fisica, Università degli Studi di Torino,
Via P. Giuria 1,
10125
Torino, Italy
41
INFN – Sezione di Torino,
Via P. Giuria 1,
10125
Torino, Italy
42
INAF – IASF Milano,
Via Alfonso Corti 12,
20133
Milano, Italy
43
Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT),
Avenida Complutense 40,
28040
Madrid, Spain
44
Port d’Informació Científica, Campus UAB,
C. Albareda s/n,
08193
Bellaterra (Barcelona), Spain
45
Institute for Theoretical Particle Physics and Cosmology (TTK), RWTH Aachen University,
52056
Aachen,
Germany
46
INAF-Osservatorio Astronomico di Roma,
Via Frascati 33,
00078
Monteporzio Catone, Italy
47
Dipartimento di Fisica e Astronomia “Augusto Righi” – Alma Mater Studiorum Università di Bologna,
Viale Berti Pichat 6/2,
40127
Bologna, Italy
48
Institute for Astronomy, University of Edinburgh, Royal Observatory,
Blackford Hill,
Edinburgh
EH9 3HJ, UK
49
Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, University of Manchester,
Oxford Road,
Manchester
M13 9PL, UK
50
European Space Agency/ESRIN,
Largo Galileo Galilei 1,
00044
Frascati, Roma,
Italy
51
Université Claude Bernard Lyon 1, CNRS/IN2P3, IP2I Lyon,
UMR 5822,
Villeurbanne,
69100,
France
52
Institute of Physics, Laboratory of Astrophysics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Observatoire de Sauverny,
1290
Versoix,
Switzerland
53
UCB Lyon 1, CNRS/IN2P3, IUF, IP2I Lyon,
4 rue Enrico Fermi,
69622
Villeurbanne,
France
54
Mullard Space Science Laboratory, University College London, Holmbury St Mary, Dorking,
Surrey
RH5 6NT, UK
55
Departamento de Física, Faculdade de Ciências, Universidade de Lisboa, Edifício C8, Campo Grande,
1749-016
Lisboa, Portugal
56
Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciên-cias, Universidade de Lisboa, Campo Grande,
1749-016
Lisboa, Portugal
57
Department of Astronomy, University of Geneva,
ch. d’Ecogia 16,
1290
Versoix, Switzerland
58
INAF-Istituto di Astrofisica e Planetologia Spaziali,
via del Fosso del Cavaliere, 100,
00100
Roma, Italy
59
INFN-Padova,
Via Marzolo 8,
35131
Padova, Italy
60
Institut d’Estudis Espacials de Catalunya (IEEC), Edifici RDIT, Campus UPC,
08860
Castelldefels, Barcelona,
Spain
61
Institut de Ciencies de l’Espai (IEEC-CSIC), Campus UAB, Carrer de Can Magrans,
s/n Cerdanyola del Vallés,
08193
Barcelona,
Spain
62
INAF-Osservatorio Astronomico di Trieste,
Via G. B. Tiepolo 11,
34143
Trieste, Italy
63
Aix-Marseille Université, CNRS/IN2P3, CPPM,
Marseille,
France
64
Istituto Nazionale di Fisica Nucleare, Sezione di Bologna,
Via Irnerio 46,
40126
Bologna, Italy
65
FRACTAL S.L.N.E.,
calle Tulipán 2, Portal 13 1A,
28231
Las Rozas de Madrid, Spain
66
Dipartimento di Fisica “Aldo Pontremoli”, Università degli Studi di Milano,
Via Celoria 16,
20133
Milano, Italy
67
Institute of Theoretical Astrophysics, University of Oslo,
P.O. Box 1029
Blindern, 0315 Oslo,
Norway
68
Leiden Observatory, Leiden University,
Einsteinweg 55,
2333
CC Leiden, The Netherlands
69
Jet Propulsion Laboratory, California Institute of Technology,
4800 Oak Grove Drive,
Pasadena,
CA,
91109,
USA
70
Department of Physics, Lancaster University,
Lancaster,
LA1 4YB,
UK
71
Felix Hormuth Engineering,
Goethestr. 17,
69181
Leimen, Germany
72
Technical University of Denmark,
Elektrovej 327,
2800
Kgs. Lyngby, Denmark
73
Cosmic Dawn Center (DAWN),
Denmark
74
Institut d’Astrophysique de Paris, UMR 7095, CNRS, and Sorbonne Université,
98 bis boulevard Arago,
75014
Paris,
France
75
NASA Goddard Space Flight Center,
Greenbelt,
MD
20771, USA
76
Department of Physics and Helsinki Institute of Physics,
Gustaf Hällströmin katu 2,
00014 University of Helsinki,
Finland
77
AIM, CEA, CNRS, Université Paris-Saclay, Université de Paris,
91191
Gif-sur-Yvette, France
78
Université de Genève, Département de Physique Théorique and Centre for Astroparticle Physics,
24 quai Ernest-Ansermet,
1211
Genève 4, Switzerland
79
Department of Physics,
PO Box 64,
00014 University of Helsinki,
Finland
80
Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University of Helsinki,
Helsinki,
Finland
81
Aix-Marseille Université, CNRS, CNES, LAM,
Marseille,
France
82
NOVA optical infrared instrumentation group at ASTRON,
Oude Hoogeveensedijk 4,
7991PD, Dwingeloo,
The Netherlands
83
INFN-Sezione di Milano,
Via Celoria 16,
20133
Milano, Italy
84
Universität Bonn, Argelander-Institut für Astronomie,
Auf dem Hügel 71,
53121
Bonn, Germany
85
Dipartimento di Fisica e Astronomia “Augusto Righi" – Alma Mater Studiorum Università di Bologna,
via Piero Gobetti 93/2,
40129
Bologna, Italy
86
Department of Physics, Centre for Extragalactic Astronomy, Durham University,
South Road,
DH1 3LE,
UK
87
Université Paris Cité, CNRS, Astroparticule et Cosmologie,
75013
Paris,
France
88
University of Applied Sciences and Arts of Northwestern Switzerland, School of Engineering,
5210
Windisch,
Switzerland
89
Institut d’Astrophysique de Paris,
98bis Boulevard Arago,
75014
Paris,
France
90
Institut de Física d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology,
Campus UAB,
08193
Bellaterra (Barcelona), Spain
91
Department of Physics and Astronomy, University of Aarhus,
Ny Munkegade 120,
8000
Aarhus C, Denmark
92
Waterloo Centre for Astrophysics, University of Waterloo,
Waterloo, Ontario
N2L 3G1, Canada
93
Department of Physics and Astronomy, University of Waterloo, Waterloo,
Ontario
N2L 3G1, Canada
94
Perimeter Institute for Theoretical Physics,
Waterloo, Ontario
N2L 2Y5, Canada
95
Space Science Data Center, Italian Space Agency,
via del Politecnico snc,
00133
Roma,
Italy
96
Institute of Space Science, Str. Atomistilor, nr. 409 Măgurele, Ilfov,
077125,
Romania
97
Dipartimento di Fisica e Astronomia “G. Galilei”, Università di Padova,
Via Marzolo 8,
35131
Padova, Italy
98
Departamento de Física, FCFM, Universidad de Chile,
Blanco Encalada 2008,
Santiago,
Chile
99
INFN-Sezione di Roma, Piazzale Aldo Moro 2,
c/o Dipartimento di Fisica, Edificio G. Marconi,
00185
Roma,
Italy
100
Satlantis, University Science Park,
Sede Bld 48940,
Leioa-Bilbao, Spain
101
Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans s/n,
08193
Barcelona,
Spain
102
Infrared Processing and Analysis Center, California Institute of Technology,
Pasadena,
CA
91125, USA
103
Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, Tapada da Ajuda,
1349-018
Lisboa, Portugal
104
Universidad Politécnica de Cartagena, Departamento de Electrónica y Tecnología de Computadoras,
Plaza del Hospital 1,
30202
Cartagena, Spain
105
Centre for Information Technology, University of Groningen,
PO Box 11044,
9700
CA Groningen, The Netherlands
106
Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, CNRS, UPS, CNES,
14 Av. Edouard Belin,
31400
Toulouse,
France
107
INFN – Bologna,
Via Irnerio 46,
40126
Bologna, Italy
108
IFPU, Institute for Fundamental Physics of the Universe,
via Beirut 2,
34151
Trieste, Italy
109
INFN, Sezione di Trieste,
Via Valerio 2,
34127
Trieste TS, Italy
110
SISSA, International School for Advanced Studies,
Via Bonomea 265,
34136
Trieste TS, Italy
111
Department of Physics and Astronomy, University of British Columbia,
Vancouver,
BC
V6T 1Z1, Canada
★★ Corresponding author; ege@iac.es
Received:
20
May
2024
Accepted:
25
July
2024
We provide an early assessment of the imaging capabilities of the Euclid space mission to deeply probe nearby star-forming regions and associated very young open clusters, and in particular, to determine to which extent it can shed light into the newborn free- floating planet population. This paper focusses on a low-reddening region observed in just one Euclid pointing. The dust and gas has been cleared out from the region by the hot σ Orionis star. One late-M and six known spectroscopically confirmed L-type ultracool members in the σ Orionis cluster were used as benchmarks to provide a high-purity procedure to select new candidate members with Euclid. The exquisite angular resolution and depth delivered by the Euclid instruments allowed us to focus on bona fide point sources. A cleaned sample of σ Orionis cluster substellar members was produced, and the initial mass function (IMF) was estimated by combining Euclid and Gaia data. Our σ Orionis substellar IMF is consistent with a power-law distribution without a significant steepening at the planetary-mass end. No evidence of a low-mass cutoff is found down to the detection limit of this study at 4 Jupiter masses in the very young σ Orionis open cluster.
Key words: catalogs / surveys / astrometry
© The Authors 2025
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1 Introduction
The nearest star-forming regions provide us with a natural laboratory to investigate the complex processes in detail that transform molecular clouds into stellar- and substellar-mass objects. In particular, one of the long-standing questions is whether there is a low-mass cutoff in the IMF, which was originally defined by Salpeter (1955) as a single power-law function over the mass range from 10 down to 0.4 solar masses (M⊙). While the early computations of spherical collapse including dust-grain opacities found a minimum mass of 0.1 M⊙ owing to opacity-limited fragmentation, that is, above the substellarmass limit (Silk 1977), recent calculations predict that the minimum fragment could reach down to 10−3 M⊙ (Mondal & Chattopadhyay 2019), which is well below the deuterium- burning mass limit. The thermonuclear fusion of deuterium, 2H(p,γ)3He, takes place at 106 K and can be strong in the early stages of evolution of objects with masses above 13 times the mass of Jupiter (1 MJ = 0.000955 M⊙) (see Bodenheimer 1966, Stahler 1988, and Chabrier et al. 2000).
Deep observations of stellar nurseries, very young open clusters, and young stellar associations have been made to search for the predicted low-mass cutoff of the IMF, and it was reported that the IMF extends smoothly into the realm of planetary masses, reaches down to the deuterium limit and overlaps with the masses of exoplanets. Various names have been used to refer to these unexpected substellar-mass objects, such as brown dwarfs (BDs) of planetary mass, sub-brown dwarfs, cluster planets, nomadic worlds, free-floating planets (FFPs), rogue planets, and planetary-mass objects (PMOs). Collectively, substellar-mass objects are ultracool dwarfs (UCDs) with very cool effective temperatures, late spectral types, small sizes, and faint luminosities that cause them to appear to be a tiny minority among the myriad stars and galaxies in deep astronomical surveys, even though their numbers can be significant.
Free-floating planets appear to be ubiquitous and numerous because they have been identified by direct imaging and spectroscopy in many different stellar cradles. Some examples of these targets are the Chamaeleon I star-forming region (Oasa et al. 1999; Luhman et al. 2004); the IC 348 and NGC 1333 clusters in Perseus (Esplin & Luhman 2017; Scholz et al. 2023), the Ophiuchus star-forming region (Chiang & Chen 2015; Bouy et al. 2022), the Orion nebula cluster (Lucas & Roche 2000; Lucas et al. 2001, 2006), the Lynds 1630 molecular clouds (Spezzi et al. 2015), the σ Orionis (Zapatero Osorio et al. 2000; Lodieu et al. 2009) and Collinder 69 (Bayo et al. 2011) young clusters in the Orion giant star formation complex, the Upper Sco OB association (Lodieu et al. 2018, 2021; Miret-Roig et al. 2022), and the Taurus dark clouds (Esplin & Luhman 2019). PMOs have also been found as wide companions to stars and BDs (Chauvin et al. 2005; Gauza et al. 2015), as members of young moving associations (Zhang et al. 2021), and as microlensing events towards the Galactic bulge (Mróz et al. 2018; Koshimoto et al. 2023; Sumi et al. 2023).
The existence of FFPs challenges models of star and planet formation. A variety of physical mechanisms have been proposed to explain the formation of substellar objects with masses well below the Jeans limit, the leading one being turbulent fragmentation (Padoan & Nordlund 2004; Hennebelle & Chabrier 2008), but others, including gravitational collapse in filaments, ejection from proto-planetary discs, and photo-erosion (Miret-Roig 2023), have not been discarded as potential players.
The cosmology-driven requirements of the Euclid mission (Laureijs et al. 2011) and the performance of its VISible instrument (VIS; Euclid Collaboration: Cropper et al. 2025) and Near-Infrared Spectrometer and Photometer (NISP; Euclid Collaboration: Jahnke et al. 2025) are expected to enable a major leap in sensitivity gain and area coverage that will foster the advance of many areas of legacy science in astrophysics (Euclid Collaboration: Mellier et al. 2025), including the detection of about one million UCDs over a large portion of the Milky Way (Solano et al. 2021; Martin et al. 2023), with spectroscopic reconnaissance spectra for thousands of them (Martín et al. 2021; Zhang et al. 2024). The Euclid reference observing sequence (ROS) is the main observation mode that is used for the wide and deep surveys. It is required to reach limiting AB magnitudes of 26.2 in the optical IE band and of 24.5 in the near-infrared NISP bands over a wide area (Euclid Collaboration: Scaramella et al. 2022).
The Euclid Early Release Observations (ERO) programme has been designed to be a showcase of the potential for legacy science across a wide range of sky regions. It demonstrates that Euclid brings a unique combination of unprecedented sensitivity, wide-area coverage, and high spatial resolution to the investigation of diverse science topics. The first ERO papers include studies of very high redshift objects (Weaver et al. 2025), clusters of galaxies (Atek et al. 2025; Kluge et al. 2025; Marleau et al. 2025; Saifollahi et al. 2025), nearby galaxies (Hunt et al. 2025), and galactic globular clusters (Massari et al. 2025).
This ERO paper investigates the power of Euclid to probe deep into very young regions over a wide area, reaching detection limits capable of revealing the FFP population and searching for the predicted low-mass cutoff of the IMF. The paper is structured as follows. In Sect. 2, the general Euclid ERO project 2 (ERO02; P.I. Martín) is presented. Five Euclid pointings were obtained, and this work focuses on about half of the area that is covered by one of them. In Sect. 3, the particular region that is the focus of this work is described, and in Sect. 4 we discuss previously known substellar-mass objects in the σ Orionis cluster and present the cuts we used to select new FFP candidates. Section 5 describes the revised IMF of the σ Orionis cluster in the area covered by the Euclid observations and compares it with the field IMF low-mass tail. Finally, Sect. 6 summarises our results and provides future prospects.
2 The Euclid Early Release Observations project of nearby star-forming regions
This Euclid ERO programme has targeted nearby (distance ≤400 pc) star-forming regions and very young open clusters (age <10 Myr) to explore their faint ultracool populations, search for FFPs, and determine whether there is an IMF low- mass cutoff. The total project consists of five Euclid pointings. The targets were the NGC 1333 cluster in Perseus (incomplete dataset), Barnard 30, Barnard 33 (the Horsehead nebula, which also includes the NGC 2023 embedded cluster and part of the σ Orionis open cluster), the Messier 78 dark clouds in the Orion star formation complex, and a field containing several dark clouds in the Taurus region. In this paper, we focus on one of these targets, called the Horsehead field, and in particular, we focus on about half of the area, hereafter called the ERO-SOri field. The other regions covered by the Euclid ERO pointings will be the subject of future studies.
The Euclid observation of the ERO-SOri field took place on 2 October 2023. A full ROS with good guiding was obtained. The centre coordinates of each of the four Euclid exposures that make up the ROS were 85°.150915, −2°.613342, 85°.167068, −2.582078, 85°.166265, −2°.551255, and 85°.182417, −2°.519991. The full field of view (FoV) of the ERO pointing presented here is displayed in Fig. 1 and covers an area of 0.58 square degrees. The FoV was chosen to avoid the blinding star σ Orionis and to include the Barnard 33 molecular cloud, the NGC 2023 cluster and reflection nebula, and the IC 434 H II region. The full Euclid ROS consisted of four dithered exposures in VIS and NISP using the nominal exposure times described in Euclid Collaboration: Cropper et al. (2025) and Euclid Collaboration: Jahnke et al. (2025), respectively. The dithering pattern was designed so that the gaps between the detectors can be covered when a stack of the four images is made. However, due to a failure in the implementation of the dithering during the science-verification phase, the pattern was not optimised, and there are some gaps in the mosaic of this ERO footprint. Furthermore, during the data processing, it was realised that about 5% of the FoV was covered by only one image and that cosmic rays could not be removed efficiently. After data reduction and image stacking, following the procedures described in Cuillandre et al. (2025), the data were validated and considered ready for scientific exploitation. In this work, the catalogues and images of the ERO public data release are used (Euclid Early Release Observations 2024). They do not include any spectroscopic data.
![]() |
Fig. 1 Multi-colour mosaic of the Euclid pointing studied in this work. The area covers 0.58 square degrees. The dark neck of the Horsehead points towards the bright σ Orionis star that is located just outside the field of view. The bright nebular emission crossing the image is the IC 434 H II region, and the bright concentration in the upper left corner is NGC 2023. This paper focuses on the low-reddening part of the field of view that was cleared out by the hot σ Orionis star. |
3 Region covered in the early release observation
The ERO pointing shown in Fig. 1 contains the complex region created by the interaction between the hot σ Orionis star and the Orion B giant molecular cloud (Lynds 1630). Extreme ultraviolet radiation from the O-type star σ Orionis creates a bright ionisation front that is known as the IC 434 H II region. A complicated pattern of bright and dark regions is clearly seen in fine detail in the Euclid mosaic (Fig. 1). The Horsehead nebula is projected in the foreground of the H II region at a distance of about 360 pc, and it points towards the ionising σ Orionis star that is in the background (Bally et al. 2018). Another ionisation source in the ERO-SOri FoV is the B-type star HD 37903, which illuminates a reflection nebula and is associated with the embedded open cluster NGC 2023, which contains very young low-mass stars (Depoy et al. 1990; Mookerjea et al. 2009; Kounkel et al. 2017). As a consequence of the complex previous and ongoing star-formation processes, there are patches with significant interstellar reddening. The AV extinction values for all the sources identified in this ERO pointing were calculated with two extinction maps from the literature: the generalised needlet internal linear combination map from Planck Collaboration Int. XLVIII (2016), which is a 2D extinction map, and the 3D extinction map Bayestar19 (Green 2018; Green et al. 2019), for which we assumed a mean distance of 400 pc to the Orion star-forming region. Both extinction maps were queried with the dustmap1 package. The cumulative distribution function of the AV extinction of our sources is shown in Fig. 2. The two extinction maps agree quite well and show that about 70% of the sources in the FoV have AV values between 0 and 1.7 mag. In the future, we plan to use more specific methods that take the extinction of the individual sources into account (e.g., Olivares et al. 2021).
![]() |
Fig. 2 Cumulative distribution of interstellar reddening in the Euclid ERO region shown in Fig. 1. Two different methods of estimating the reddening are compared. About 70% of the FoV has a modest reddening of less than 1.7 mag in the visible. This reddening value is used in the reddening vectors shown in the colour–magnitude and colour– colour diagrams of the highly reddened region. This work focuses on the σ Orionis cluster, which is located in the low-reddening region of the FoV. |
4 Euclid view of the σ Orionis substellar members
The Euclid footprint of the ERO pointing includes a portion of the well-known σ Orionis cluster, which has been a favourite hunting ground for very young substellar objects and FFPs for over two decades (Zapatero Osorio et al. 2000; Damian et al. 2023). A review of the σ Orionis cluster properties was provided by Walter et al. (2008). A recent assessment of cluster membership using the Gaia third data release (DR3, Gaia Collaboration 2023) has been carried out in a study of the young populations in the region (Žerjal et al. 2024). The ages of most σ Orionis cluster members are in the range from 1 Myr to 5 Myr (Zapatero Osorio et al. 2002). We adopted an age of (3 ± 2) Myr and a distance of 402.7±9.0 pc (Žerjal et al. 2024) for the σ Orionis cluster. The deepest survey carried out to date in the search for FFPs belonging to σ Orionis was reported by Peña Ramírez et al. (2012) using ground-based telescopes.
![]() |
Fig. 3 Mosaic of Euclid images centred on the benchmark object S Ori 52 in the four different photometric passbands. A clearly resolved visual companion is visible in the VIS image. |
4.1 Definition of benchmarks for Euclid based on confirmed σ Orionis substellar objects
We selected seven confirmed very cool members of the σ Orionis cluster with a ground-based low-resolution optical and near-infrared spectroscopic classification. Their names and coordinates are listed in Table A.1, together with the spectral types from the literature (Zapatero Osorio et al. 2000; Barrado y Navascués et al. 2001; Martín et al. 2001; Zapatero Osorio et al. 2017). The parameters of these seven benchmarks in the Euclid ERO catalogue are provided in Table A.2.
The values of the SPREAD-MODEL parameter for the benchmarks in all the Euclid passbands deviate very little from zero, as expected for bona fide point sources. This parameter was developed as a star and galaxy classifier by the data-management pipeline of the Dark Energy Survey (Mohr et al. 2012), and it was shown to be a good discriminant for point sources in nearby young clusters and stellar associations (Bouy et al. 2013). The parameter was adopted as one of the main selection criteria to separate point sources from galaxies.
We note that two benchmarks, namely S Ori 52 and 60, have slightly higher FWHM IMAGE I values in the ERO catalogue than the other four benchmarks, suggesting that they might be binaries with an angular separation close to the limit of the spatial resolution of the VIS data. Moreover, a resolved faint visual companion was spotted close to S Ori 52 at position J054009.36–022631.93 in the VIS image (see Fig. 3). SOri52 has an optical spectral type of L0.5 and a mass of about 15 MJ (Béjar et al. 2001). The candidate wide companion to S Ori 52 is 3.2 mag fainter in the IE passband than the primary, the angular separation is 0.96″ (387 au at 403 pc), and the position angle is 43.3°. The pair is barely resolved in the NISP images because they have a lower spatial resolution than the VIS image. Using PSF photometry, the difference in magnitude in the JE passband is 3.96. This difference is larger than in IE, indicating that the companion has a slightly bluer IE–JE colour than the primary, and casting some doubt on the physical association of these two objects. The possibility that Euclid may have found two substellar binaries close to the angular resolution of the VIS images in a sample of only seven benchmarks in the σ Orionis cluster is interesting and deserves further scrutiny. A Hubble Space Telescope (HST) imaging survey of wide binaries (projected semi-major axes between 100 and 1000 au) among pre-main-sequence (PMS) stars in the Orion star-forming complex found a binary frequency of (Kounkel et al. 2016), and recent work with the James Webb Space Telescope (JWST) suggested that substellar-mass binaries in the Trapezium cluster could be common (McCaughrean & Pearson 2023). On the other hand, no resolved binaries with separations >20 au were found in an imaging survey of 33 BDs in two young open clusters (ages in the range from 70 to 120 Myr) that was carried out with the HST (Martín et al. 2003).
4.2 Contamination estimates and definition of selection cuts for the Euclid data
To assess the likelihood that the substellar object candidates are contaminated by background extragalactic objects and by foreground ultracool dwarfs, all the objects that were not saturated in the Euclid images and listed by Peña Ramírez et al. (2012) in the ERO-SOri FoV were visually inspected in the VIS and NISP images and were cross-correlated with the ESO VISTA Hemisphere Survey (VHS) catalogue (McMahon et al. 2021) to check for proper motions. The total proper motion of true cluster members is expected to be ≤20 mas per year and thus should not be measurable when comparing VHS and NISP data. A summary of the results of this contamination assessment is provided in Table A.3. Sources that are spatially resolved as extended objects in any of the Euclid passbands are considered as non-members. They have values of the FWHM and SPREAD_MODEL parameters larger than those of the benchmarks. Sources that are detected to move by ≥100 mas from the VHS epoch to the Euclid epoch (baseline 14 years) are classified as non-members and are labelled as high proper motion. In particular, the Euclid VIS spatial resolution has been crucial to show that some of the very faint and red sources identified in Peña Ramírez et al. (2012) are likely to be galaxies. For the benchmarks, we confirmed that their coordinates match those from the VHS catalogue within 100 mas. As expected, the most frequent contamination comes from extended objects that are probably background galaxies (9/38 or 24%), particularly at the faint end of the sample. The ratio of extragalactic sources to ultracool dwarfs is expected to increase with increasing depth. It has recently been reported from JWST Near-Infrared Spectrograph spectroscopic follow-up of photometrically selected JWST Near-Infrared Camera compact sources that the ratio of extragalactic objects and ultracool dwarfs is 11/3 at depths fainter than those reached by the Euclid images (Langeroodi & Hjorth 2023).
The contamination by background extragalactic sources, the inhomogeneous interstellar extinction in star-forming regions, the possible presence of colour excesses owing to discs and accretion activity, and the low surface gravity and extreme youth of FFPs in Orion together make the selection of bona fide sub- stellar objects quite challenging. The Euclid passbands are not specifically designed to distinguish FFPs from other types of objects. They are broader than the passbands commonly used in ground-based surveys because they include spectral regions that are affected by saturated telluric water absorption. The calibrations available for this ERO study are scarce. Improved calibrations are expected in the future when Euclid photometry and spectroscopy of benchmark ultracool dwarfs become available. We limited the scope of this work to present a high-purity approach to select objects using the Euclid ERO catalogue and images that is anchored in the properties of the benchmarks described in the previous section. The selection cuts adopted in this work are presented in Table A.4. To arrive at these selection cuts, we calculated the 1 σ dispersion around the mean of the values for the benchmark sources and added it to both extremes of the distribution.
The ERO-SOri photometric catalogue was filtered using the cuts provided in Table A.4. The number of objects left after each step and the percentage with respect to the original sample are also given in the table. The percentage of sources detected in the JE band in the whole FoV is 76.83% of the total number of sources in the ERO catalogue. Sources not detected in the JE band were excluded from this work because young substel- lar objects are expected to be much brighter in the JE band than in the IE band. The CLASS-STAR classifier was found to be redundant with the SPREAD-MODEL parameter, and the latter was chosen because the values for the benchmarks are more stable. The FoV was divided into two regions separated by a constant RA value of 85.1875°. We call the low-reddening and dark-background part the σ Orionis region (RA < 85.1875°) and the high-reddening and bright-background part the Horsehead region. We compared the distribution of SPREAD-MODEL_J between these two parts of the FoV and found that it is narrower in the σ Orionis region than in the Horsehead region (Fig. B.1). This is likely due to the influence of a brighter background on the Horsehead side owing to light that is reflected in the nebulosity. Thus, we consider that the cuts defined in this work are valid only for regions with low interstellar background and negligible extinction. For the regions affected by high sky background, it will be important to obtain a new sample of substellar benchmarks using the Euclid NISP spectra and to consider the effects of variable background on the PSF of the sources.
4.3 Selection of new substellar member candidates in the σ Orionis cluster with Euclid data
After applying all the selection cuts defined above, only 2% of the sources in the initial catalogue remained. They are plotted in the IE versus IE–JE colour–magnitude diagram (CMD) following the approach of previous searches for substellar objects in the σ Orionis region (Peña Ramírez et al. 2012), and were compared with the 3-Myr isochrone provided by the ATMO models of Phillips et al. (2020), which were transformed into the Euclid photometric system for this work. These models were tested using the dynamic lithium-boundary method for brown dwarf binaries with dynamical masses and were found to fit the observational data better than other sets of models in the literature (Martín et al. 2022). The Euclid data and the CEQ (equilibrium chemistry) ATMO 3-Myr isochrone are shown in the CMD displayed in Fig. 4. The benchmarks clearly define the σ Orionis sequence, and other objects in the Euclid data that appear to follow this cluster sequence were identified previously as photometric candidate members (Peña Ramírez et al. 2012). Their Euclid coordinates and photometry are provided in Table A.5. Seven new objects were found to be located close to the cluster sequence and are well separated from the cloud of background sources, including two very faint objects that extend the sequence to fainter magnitudes than previous surveys. The three brightest objects in the NISP JE band were retrieved in the VHS catalogue, and their coordinates were found to agree within 100 mas. This means that they do not have a high proper motion. The coordinates and photometry of these seven new Euclid objects of interest identified in the σ Orionis cluster sequence are given in Table A.6.
Further examination of the cluster sequence, its degree of agreement with the ATMO isochrone, and the location of new candidate members was made in the colour–colour diagram shown in Fig. 5. The coolest benchmarks define a well-separated locus away from the cloud of contaminating sources. The behaviour of the benchmarks is qualitatively fairly well reproduced by the isochrone, although quantitatively, the fit could be improved because the isochrone does not reach as large a YE – HE colour as observed. The blueing of the isochrone in the YE – HE colour beyond IE – YE ≥ 3.5 is an effect of the appearance of methane in the transition from L- to T-type spectra. The new sources with codes D, E, and G fall within the L-type benchmark locus, making them strong FFP candidates. Source F is slightly bluer in YE–HE than the faintest benchmark, and it is also closer to the isochrone, suggesting that it might be the first L/T transition FFP identified in the σ Orionis cluster. Confirmation of these tentative assessments requires spectroscopy.
To determine the effects of reddening in the selection of substellar candidates, the same cuts as applied to the σ Orionis region were also applied to the Horsehead regions. The CMD and colour–colour diagrams are shown in Figs. 6 and 7, respectively. The separation between the cluster sequence defined by the benchmarks and the cloud of sources is no longer well defined in the CMD, and the locus of benchmarks in the colour–colour diagram is not well isolated. This example shows that it is difficult to select substellar candidates in regions with high interstellar reddening. Future work will address this issue.
To search for binaries, we show in Fig. B.2 the FWHM values (in pixels) versus aperture magnitudes measured in the VIS images for all the objects under study (benchmarks, confirmed candidates in the cluster sequence, and new discoveries). In addition to the two binary candidates among the benchmarks, one more candidate is found among the confirmed objects and another among the new objects that were found with Euclid. The object labelled G could be the first σ Orionis counterpart to the Jupiter-mass binary candidates reported in the Trapezium cluster (McCaughrean & Pearson 2023), but it needs confirmation with higher spatial resolution images that might be provided by HST optical or JWST near-infrared imaging observations.
![]() |
Fig. 4 IE vs. IE – JE colour–magnitude diagram for the σ Orionis part of the FoV. The black points are all the Euclid sources that remain after applying all the cuts. Benchmark objects are denoted with red circles. New objects near the cluster sequence are denoted with blue squares and labelled with capital letters. Known sources other than the benchmarks in the σ Orionis cluster are denoted with green circles. An ATMO CEQ isochrone (see Phillips et al. 2020) for an age of 3 Myr and a distance of 402.74 pc is shown. Theoretical masses are labelled on the isochrone. |
![]() |
Fig. 5 IE–YE vs. YE–HE Colour–colour diagram for the same region as in the previous figure (σ Orionis). All the symbols remain the same. |
![]() |
Fig. 6 IE vs. IE – JE colour–magnitude diagram for the Horsehead part of the FoV (not σ Orionis). The σ Orionis benchmarks are denoted with red circles. A reddening vector is shown as an arrow. |
![]() |
Fig. 7 IE – YE vs. YE – HE colour–colour diagram for the Horsehead part of the FoV (not σ Orionis). The σ Orionis benchmarks are denoted with red circles. A reddening vector is shown as an arrow. |
5 Euclid substellar initial mass function of the σ Orionis cluster
Our results are used here to revise the very low mass IMF of the σ Orionis cluster and try to extend it deeper into the planetary- mass regime. The Gaia sample covers the domain of very low mass stars and Euclid provides a continuation into the substellarmass regime, reaching down to about 4 MJ.
The mass-luminosity relation from the 3 Myr ATMO CEQ models was used for the substellar domain. The PAdova and TRieste Stellar Evolution Code (PARSEC) models (Bressan et al. 2012; Pastorelli et al. 2020) were used for the stellar domain. The IMF of the σ Orionis cluster using the Gaia-DR3 membership study by Žerjal et al. (2024), combined with the results of this work, is displayed in Fig. 8. Our results are consistent with a multi-power-law distribution for the IMF.
A comprehensive study of the IMF within a distance of 20 pc from the Sun has reported a change in the slope in the substellar domain (Kirkpatrick et al. 2024). The authors claimed that the solar vicinity IMF can be expressed as dN/dM = C M−α with four different values of the power-law exponent for different mass intervals. In particular, we are concerned with the low- mass tail of the IMF, where the slope estimated by Kirkpatrick et al. (2024) steepens from a value of α = 0.25 in the mass range 0.05 M⊙ < M < 0.22 M⊙ to α = 0.60 in the mass range 0.01 M⊙ < M < 0.05 M⊙.
We identified three different mass regimes here: the very low mass stellar domain from 0.15 to 0.1 M⊙ with α = 0.26±0.10, the brown dwarf domain from 0.1 to 0.011 M⊙ with α = 0.18±0.01, and the planetary-mass domain from 0.011 to 0.003 M⊙ with α = 0.12±0.02. These values were obtained with the linear fits that are shown in Fig. 8. We excluded from the fits the mass range between 0.1 M⊙ and 0.05 M⊙ because these objects are too faint to be complete for Gaia and too bright for Euclid. The error bars quoted for the IMF slopes were estimated by simulations of the effects of age, distance, and photometric uncertainties.
Our σ Orionis IMF results are consistent with the field in the very low mass stellar regime and extend deeper into the substel- lar regime than the field IMF. We do not confirm a steepening of the substellar IMF at the planetary-mass end. The slope of our IMF in the substellar regime is shallower than that reported in the previous study of the σ Orionis cluster by Peña Ramírez et al. (2012) and also shallower than the IMF study of the solar vicinity by Kirkpatrick et al. (2024), but it is unclear if the difference is significant because there are uncertainties that have not been included in this study. The error bars quoted in our IMF determination only take the internal uncertainty into account when the data were fit with power laws. The comparisons between different substellar IMFs are affected by low number statistics and other uncertainties (e.g. unresolved binaries and mass determination using evolutionary models), particularly at the planetary-mass end. The census of directly imaged FFPs must be increased significantly to investigate the possibility of substellar IMF variations in different environments that could be an indication of specific formation pathways in the planetarymass domain. These results demonstrate that Euclid can play a significant role in the detailed study of the low-mass shape of the IMF, and particularly in shedding light on the formation mechanisms of FFPs. Detailed theoretical models developed by different groups have indicated that the shape of the IMF is a useful indicator of the dominant mode of star formation in a given region (Adams & Fatuzzo 1996; Chabrier 2005; Thies et al. 2015), and that a multi-power-law IMF could arise from the interplay between the mass- and time-dependence of exponential growth in a distribution of accreting protostars (Essex et al. 2020).
The photometric completeness of the Euclid substellar IMF of the σ Orionis cluster was assessed using the source number count distribution shown in Fig. 9. Before any filtering of the catalogue, the depth reached with the VIS instrument indicated that the Euclid survey provides a complete detection of objects down to a VIS magnitude of about 27 and to NISP magnitudes of about 25. These values are consistent with the 5σ values provided by the ERO study of extragalactic fields (Atek et al. 2025). This photometric depth corresponds to 4 MJ for the age and distance of the σ Orionis cluster; but after all the filtering needed to select bona fide point sources, the number of objects retained in VIS started to decline at a VIS magnitude of 25, corresponding to masses of about 6 MJ . From these considerations, we estimate that the Euclid substellar IMF completeness limit of this work is located at 6 MJ , and the detection limit is at 4 MJ for the σ Orionis cluster.
The substellar IMF presented here does not take the likely presence of unresolved binaries in our sample into account. To assess the impact of unresolved binaries in our study, we generated random simulations of binaries with a frequency of 20% of equal-mass systems. In the Pleiades cluster, the substellar binary frequency with a semimajor axis in the range 7 to 12 au and mass ratios higher than 0.7 was estimated to be up to 20% in a survey carried out with the HST (Martín et al. 2003). A population of binaries in the a Orionis cluster like this would remain unresolved with Euclid. We recalculated the IMF for each random simulation of unresolved binaries. The average IMF resulting from ten simulations is shown in Fig. B.3. The main effect is that the slopes of the linear fits become steeper because there are more substellar objects in the sample and their masses tend to be lower than the unresolved binary system masses.
We note that both the photometric uncompleteness of the sample after filtering and the effects of unresolved binaries together contribute to our underestimating of the number of FFPs in our study, particularly at the low-mass end below 6 MJ . We therefore consider our conclusion to be solid that there is no evidence for an IMF cutoff down to the detection limit of this study at 4 MJ .
![]() |
Fig. 8 Combined very low mass Euclid–Gaia IMF of the σ Orionis cluster with linear fits in three different mass regions. |
![]() |
Fig. 9 Number counts of Euclid sources in the catalogue of the σ Orionis region before and after filtering. The correspondence between VIS magnitudes and masses has been made using the ATMO 3 Myr isochrone. |
6 Final remarks: Impact of Euclid on the study of free-floating planets in star-forming regions
This work is a showcase of the power of the Euclid mission to provide the area and depth required to explore the very low mass population, including FFPs of nearby star-forming regions and very young open clusters. In particular, for the well-known σ Orionis cluster, we showed that the sensitivity of the Euclid images is capable of probing down to FFPs that could have masses as low as 4 MJ according to theoretical models (for 3 Myr ages) and at a distance of 400 pc. This potential could be compromised by severe contamination from numerous background extragalactic sources if we were not to use stringent selection procedures. Using the Euclid data for seven benchmark objects in σ Orionis, we developed a high-purity method to filter out the contamination. This method is valid for regions of low reddening, but it needs additional work to be generalised to regions with any reddening. We note that the Euclid NISP spectra will likely play an important role in this effort. The ERO images, catalogues, and spectra processed with the official pipeline are planned to be released in March 2025 and they will enable a reassessment of the results presented here. Furthermore, deeper observations of another region in Orion will be made available to the community. Additionally multi-epoch observations with Euclid during the lifetime of the mission, possibly filling gaps in the cosmological surveys, can enable the study of proper motions and photometric variability that are useful probes for the study of the low-mass population in star-forming regions.
This paper is a showcase of the potential of Euclid to explore nearby star-forming regions. The observations presented here provide a first glimpse of the power of Euclid to shed light on the long-standing question of the putative low-mass cutoff of the IMF predicted a long time ago by the theory of opacitylimited fragmentation and collapse of molecular clouds. Our IMF for the σ Orionis open cluster extends previous studies to lower planetary masses and suggests that there is no indication that the predicted cutoff at the low-mass end is near the Euclid detection limit. Another open question that could be addressed in the future with Euclid is the degree of universality of the IMF slope in the planetary-mass regime when comparing different young open clusters, star-forming regions and the field. This work demonstrated the great potential of Euclid to determine the substellar IMF down to the FFP regime in nearby star-forming regions and very young open clusters, and it also paves the way to overcome the difficulties associated with the study of very young regions with deep and wide optical and near-infrared space-based imaging observations.
Acknowledgements
This work has made use of the Early Release Observations (ERO) data from the Euclid mission of the European Space Agency (ESA), 2024. We thank I. Baraffe and M. Phillips for making available a digitised version of their isochrones in the Euclid passbands and the referee (Joao Alves) for insightful comments that helped us to improve the manuscript. ELM., MŽ, AE, CDT, SMT, NS, ST, and NV are supported by the European Research Council Advanced grant SUBSTELLAR, project number 101054354. This research has made use of the Spanish Virtual Observatory (https://svo.cab.inta-csic.es) project funded by MCIN/AEI/10.13039/501100011033/ through grant PID2020-112949GB-I00 at Centro de Astrobiología (CSIC-INTA). DB and NH have been supported by PID2019-107061GB-C61 by the same agency. PMB is funded by Instituto Nacional de Técnica Aeroespacial through grant PRE-OVE. PC acknowledges financial support from the Spanish Virtual Observatory project (grant PID2020-112949GB-I00). NPB is funded by Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.99-2020.63. NL and VJSB acknowledge financial support from the Agencia Estatal de Investigación (AEI/10.13039/501100011033) of the Ministerio de Ciencia e Innovación and the ERDF ‘A way of making Europe’ through project PID2022-137241NB-C1. CdB acknowledges support from a Beatriz Galindo senior fellowship (BG22/00166) from the Spanish Ministry of Science, Innovation and Universities. AA acknowledges support from the Light Bridges Corporation. The Euclid Consortium acknowledges the European Space Agency and a number of agencies and institutes that have supported the development of Euclid, in particular the Agenzia Spaziale Italiana, the Austrian Forschungsförderungsgesellschaft funded through BMK, the Belgian Science Policy, the Canadian Euclid Consortium, the Deutsches Zentrum für Luft- und Raumfahrt, the DTU Space and the Niels Bohr Institute in Denmark, the French Centre National d’Etudes Spatiales, the Fundação para a Ciência e a Tecnologia, the Hungarian Academy of Sciences, the Ministerio de Ciencia, Innovación y Universidades, the National Aeronautics and Space Administration, the National Astronomical Observatory of Japan, the Netherlandse Onderzoekschool Voor Astronomie, the Norwegian Space Agency, the Research Council of Finland, the Romanian Space Agency, the State Secretariat for Education, Research, and Innovation (SERI) at the Swiss Space Office (SSO), and the United Kingdom Space Agency. A complete and detailed list is available on the Euclid web site (http://www.euclid-ec.org).
Appendix A Tables
Benchmark σ Orionis cluster members observed with Euclid.
Euclid parameters for σ Orionis benchmark objects. FWHM values are given in pixels.
Euclid-based assessment of previously identified candidate members of the σ Orionis cluster and very red sources in the region.
Euclid point-source selection criteria applied in this work. The cut in right ascension, RA, divides the area into two parts: the high- reddening region (RA < 85.1875); and the low reddening region (RA ≥ 85.1875).
Euclid photometry of previously known non-benchmark objects in the σ Orionis cluster sequence.
Euclid objects of interest in the σ Orionis region.
Appendix B Additional figures
![]() |
Fig. B.1 Cumulative distribution of spread model values for the Euclid JE band. The distribution of values in the low-reddening part corresponding to the σ Orionis cluster is sharper than in the high-reddening part. |
![]() |
Fig. B.2 Euclid VIS FWHM (in pixels) versus IE apparent magnitude for objects confirmed by Euclid to be in the σ Orionis cluster sequence. Both the FWHM and IE photometry values come from the ERO catalogue. Four objects were found to have FWHM values larger than the mean value of cluster members (more than 3 σ significance) and hence are considered as possible binaries that deserve further scrutiny. Two of them are benchmark objects, namely (S Ori 52 and S Ori 60), one is a known object (S Ori 38), and the faintest one is a new discovery (Euclid object of interest G). |
![]() |
Fig. B.3 Unresolved binary-corrected very low mass Euclid-Gaia IMF of the σ Orionis cluster with linear fits in three different mass regions. A 20% unresolved binary frequency of equal mass binaries has been assumed in the simulations. |
References
- Adams, F. C., & Fatuzzo, M. 1996, ApJ, 464, 256 [NASA ADS] [CrossRef] [Google Scholar]
- Atek, H., Gavazzi, R., Weaver, J. R., et al. 2025, A&A, 697, A15 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bally, J., Chambers, E., Guzman, V., et al. 2018, AJ, 155, 80 [Google Scholar]
- Barrado y Navascués, D., Zapatero Osorio, M. R., Béjar, V. J. S., et al. 2001, A&A, 377, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bayo, A., Barrado, D., Stauffer, J., et al. 2011, A&A, 536, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Béjar, V. J. S., Martín, E. L., Zapatero Osorio, M. R., et al. 2001, ApJ, 556, 830 [CrossRef] [Google Scholar]
- Bodenheimer, P. 1966, ApJ, 144, 103 [Google Scholar]
- Bouy, H., Bertin, E., Moraux, E., et al. 2013, A&A, 554, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bouy, H., Tamura, M., Barrado, D., et al. 2022, A&A, 664, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G. 2005, in Astrophysics and Space Science Library, 327, The Initial Mass Function 50 Years Later, eds. E. Corbelli, F. Palla, & H. Zinnecker, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 2000, ApJ, 542, L119 [CrossRef] [Google Scholar]
- Chauvin, G., Lagrange, A. M., Dumas, C., et al. 2005, A&A, 438, L25 [CrossRef] [EDP Sciences] [Google Scholar]
- Chiang, P., & Chen, W. P. 2015, ApJ, 811, L16 [Google Scholar]
- Cuillandre, J.-C., Bertin, E., Bolzonella, M., et al. 2025, A&A, 697, A6 (Euclid on Sky SI) [Google Scholar]
- Damian, B., Jose, J., Biller, B., et al. 2023, ApJ, 951, 139 [NASA ADS] [CrossRef] [Google Scholar]
- Depoy, D. L., Lada, E. A., Gatley, I., & Probst, R. 1990, ApJ, 356, L55 [NASA ADS] [CrossRef] [Google Scholar]
- Esplin, T. L., & Luhman, K. L. 2017, AJ, 154, 134 [Google Scholar]
- Esplin, T. L., & Luhman, K. L. 2019, AJ, 158, 54 [Google Scholar]
- Essex, C., Basu, S., Prehl, J., & Hoffmann, K. H. 2020, MNRAS, 494, 1579 [Google Scholar]
- Euclid Collaboration (Scaramella, R., et al.) 2022, A&A, 662, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Euclid Collaboration (Cropper, M. S., et al.) 2025, A&A, 697, A2 (Euclid on Sky SI) [Google Scholar]
- Euclid Collaboration (Jahnke, K., et al.) 2025, A&A, 697, A3 (Euclid on Sky SI) [Google Scholar]
- Euclid Collaboration (Mellier, Y., et al.) 2025, A&A, 697, A1 (Euclid on Sky SI) [Google Scholar]
- Euclid Early Release Observations 2024, https://doi.org/10.57780/esa-qmocze3 [Google Scholar]
- Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gauza, B., Béjar, V. J. S., Pérez-Garrido, A., et al. 2015, ApJ, 804, 96 [NASA ADS] [CrossRef] [Google Scholar]
- Green, G. M. 2018, J. Open Source Softw., 3, 695 [Google Scholar]
- Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93 [NASA ADS] [CrossRef] [Google Scholar]
- Hennebelle, P., & Chabrier, G. 2008, ApJ, 684, 395 [Google Scholar]
- Hunt, L. K., Annibali, F., Cuillandre, J.-C., et al. 2025, A&A, 697, A9 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kirkpatrick, J. D., Marocco, F., Gelino, C. R., et al. 2024, ApJS, 271, 55 [NASA ADS] [CrossRef] [Google Scholar]
- Kluge, M., Hatch, N. A., Montes, M., et al. 2025, A&A, 697, A13 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Koshimoto, N., Sumi, T., Bennett, D. P., et al. 2023, AJ, 166, 107 [NASA ADS] [CrossRef] [Google Scholar]
- Kounkel, M., Megeath, S. T., Poteet, C. A., Fischer, W. J., & Hartmann, L. 2016, ApJ, 821, 52 [NASA ADS] [CrossRef] [Google Scholar]
- Kounkel, M., Hartmann, L., Mateo, M., & Bailey, I., John I., 2017, ApJ, 844, 138 [NASA ADS] [CrossRef] [Google Scholar]
- Langeroodi, D., & Hjorth, J. 2023, ApJ, 957, L27 [NASA ADS] [CrossRef] [Google Scholar]
- Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, arXiv e-prints [arXiv:1110.3193] [Google Scholar]
- Lodieu, N., Zapatero Osorio, M. R., Rebolo, R., Martín, E. L., & Hambly, N. C. 2009, A&A, 505, 1115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lodieu, N., Zapatero Osorio, M. R., Béjar, V. J. S., & Peña Ramírez, K. 2018, MNRAS, 473, 2020 [Google Scholar]
- Lodieu, N., Hambly, N. C., & Cross, N. J. G. 2021, MNRAS, 503, 2265 [Google Scholar]
- Lucas, P. W., & Roche, P. F. 2000, MNRAS, 314, 858 [NASA ADS] [CrossRef] [Google Scholar]
- Lucas, P. W., Roche, P. F., Allard, F., & Hauschildt, P. H. 2001, MNRAS, 326, 695 [NASA ADS] [CrossRef] [Google Scholar]
- Lucas, P. W., Weights, D. J., Roche, P. F., & Riddick, F. C. 2006, MNRAS, 373, L60 [NASA ADS] [Google Scholar]
- Luhman, K. L., Peterson, D. E., & Megeath, S. T. 2004, ApJ, 617, 565 [Google Scholar]
- Marleau, F. R., Cuillandre, J.-C., Cantiello, M., et al. 2025, A&A, 697, A12 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martín, E. L., Zapatero Osorio, M. R., Barrado y Navascués, D., Béjar, V. J. S., & Rebolo, R. 2001, ApJ, 558, L117 [CrossRef] [Google Scholar]
- Martín, E. L., Barrado y Navascués, D., Baraffe, I., Bouy, H., & Dahm, S. 2003, ApJ, 594, 525 [CrossRef] [Google Scholar]
- Martín, E. L., Zhang, J. Y., Esparza, P., et al. 2021, A&A, 655, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martín, E. L., Lodieu, N., & del Burgo, C. 2022, MNRAS, 510, 2841 [Google Scholar]
- Martin, E. L., Bouy, H., Martin, D., et al. 2023, Proc., 30th Rencontres du Vietnam: Windows on the Universe [Google Scholar]
- Massari, D., Dalessandro, E., Erkal, D., et al. 2025, A&A, 697, A8 (Euclid on Sky SI) [Google Scholar]
- McCaughrean, M. J., & Pearson, S. G. 2023, A&A, submitted [arXiv:2310.03552] [Google Scholar]
- McMahon, R. G., Banerji, M., Gonzalez, E., et al. 2021, VizieR Online Data Catalog, II/367 [Google Scholar]
- Miret-Roig, N. 2023, Ap&SS, 368, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Miret-Roig, N., Bouy, H., Raymond, S. N., et al. 2022, Nat. Astron., 6, 89 [NASA ADS] [CrossRef] [Google Scholar]
- Mohr, J. J., Armstrong, R., Bertin, E., et al. 2012, SPIE Conf. Ser., 8451, 84510D [Google Scholar]
- Mondal, A., & Chattopadhyay, T. 2019, New A, 66, 45 [Google Scholar]
- Mookerjea, B., Sandell, G., Jarrett, T. H., & McMullin, J. P. 2009, A&A, 507, 1485 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mróz, P., Ryu, Y. H., Skowron, J., et al. 2018, AJ, 155, 121 [CrossRef] [Google Scholar]
- Oasa, Y., Tamura, M., & Sugitani, K. 1999, ApJ, 526, 336 [NASA ADS] [CrossRef] [Google Scholar]
- Olivares, J., Bouy, H., Sarro, L. M., et al. 2021, A&A, 649, A159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Padoan, P., & Nordlund, Å. 2004, ApJ, 617, 559 [NASA ADS] [CrossRef] [Google Scholar]
- Pastorelli, G., Marigo, P., Girardi, L., et al. 2020, MNRAS, 498, 3283 [Google Scholar]
- Peña Ramírez, K., Béjar, V. J. S., Zapatero Osorio, M. R., Petr-Gotzens, M. G., & Martín, E. L. 2012, ApJ, 754, 30 [Google Scholar]
- Phillips, M. W., Tremblin, P., Baraffe, I., et al. 2020, A&A, 637, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration Int. XLVIII. 2016, A&A, 596, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Saifollahi, T., Voggel, K., Lançon, A., et al. 2025, A&A, 697, A10 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
- Scholz, A., Muzic, K., Jayawardhana, R., Almendros-Abad, V., & Wilson, I. 2023, AJ, 165, 196 [Google Scholar]
- Silk, J. 1977, ApJ, 214, 152 [Google Scholar]
- Solano, E., Gálvez-Ortiz, M. C., Martín, E. L., et al. 2021, MNRAS, 501, 281 [Google Scholar]
- Spezzi, L., Petr-Gotzens, M. G., Alcalá, J. M., et al. 2015, A&A, 581, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Stahler, S. W. 1988, ApJ, 332, 804 [NASA ADS] [CrossRef] [Google Scholar]
- Sumi, T., Koshimoto, N., Bennett, D. P., et al. 2023, AJ, 166, 108 [NASA ADS] [CrossRef] [Google Scholar]
- Thies, I., Pflamm-Altenburg, J., Kroupa, P., & Marks, M. 2015, ApJ, 800, 72 [NASA ADS] [CrossRef] [Google Scholar]
- Walter, F. M., Sherry, W. H., Wolk, S. J., & Adams, N. R. 2008, in Handbook of Star Forming Regions, Volume I, 4, ed. B. Reipurth (ASP Conf. Ser.), 732 [Google Scholar]
- Weaver, J. R., Taamoli, S., McPartland, C. J. R., et al. 2025, A&A, 697, A16 (Euclid on Sky SI) [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zapatero Osorio, M. R., Béjar, V. J. S., Martín, E. L., et al. 2000, Science, 290, 103 [NASA ADS] [CrossRef] [Google Scholar]
- Zapatero Osorio, M. R., Béjar, V. J. S., Pavlenko, Y., et al. 2002, A&A, 384, 937 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zapatero Osorio, M. R., Béjar, V. J. S., & Peña Ramírez, K. 2017, ApJ, 842, 65 [CrossRef] [Google Scholar]
- Žerjal, M., Martín, E. L., & Pérez-Garrido, A. 2024, A&A, 686, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zhang, Z., Liu, M. C., Best, W. M. J., Dupuy, T. J., & Siverd, R. J. 2021, ApJ, 911, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Zhang, J. Y., Lodieu, N., & Martín, E. L. 2024, A&A, 686, A171 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Euclid parameters for σ Orionis benchmark objects. FWHM values are given in pixels.
Euclid-based assessment of previously identified candidate members of the σ Orionis cluster and very red sources in the region.
Euclid point-source selection criteria applied in this work. The cut in right ascension, RA, divides the area into two parts: the high- reddening region (RA < 85.1875); and the low reddening region (RA ≥ 85.1875).
Euclid photometry of previously known non-benchmark objects in the σ Orionis cluster sequence.
All Figures
![]() |
Fig. 1 Multi-colour mosaic of the Euclid pointing studied in this work. The area covers 0.58 square degrees. The dark neck of the Horsehead points towards the bright σ Orionis star that is located just outside the field of view. The bright nebular emission crossing the image is the IC 434 H II region, and the bright concentration in the upper left corner is NGC 2023. This paper focuses on the low-reddening part of the field of view that was cleared out by the hot σ Orionis star. |
In the text |
![]() |
Fig. 2 Cumulative distribution of interstellar reddening in the Euclid ERO region shown in Fig. 1. Two different methods of estimating the reddening are compared. About 70% of the FoV has a modest reddening of less than 1.7 mag in the visible. This reddening value is used in the reddening vectors shown in the colour–magnitude and colour– colour diagrams of the highly reddened region. This work focuses on the σ Orionis cluster, which is located in the low-reddening region of the FoV. |
In the text |
![]() |
Fig. 3 Mosaic of Euclid images centred on the benchmark object S Ori 52 in the four different photometric passbands. A clearly resolved visual companion is visible in the VIS image. |
In the text |
![]() |
Fig. 4 IE vs. IE – JE colour–magnitude diagram for the σ Orionis part of the FoV. The black points are all the Euclid sources that remain after applying all the cuts. Benchmark objects are denoted with red circles. New objects near the cluster sequence are denoted with blue squares and labelled with capital letters. Known sources other than the benchmarks in the σ Orionis cluster are denoted with green circles. An ATMO CEQ isochrone (see Phillips et al. 2020) for an age of 3 Myr and a distance of 402.74 pc is shown. Theoretical masses are labelled on the isochrone. |
In the text |
![]() |
Fig. 5 IE–YE vs. YE–HE Colour–colour diagram for the same region as in the previous figure (σ Orionis). All the symbols remain the same. |
In the text |
![]() |
Fig. 6 IE vs. IE – JE colour–magnitude diagram for the Horsehead part of the FoV (not σ Orionis). The σ Orionis benchmarks are denoted with red circles. A reddening vector is shown as an arrow. |
In the text |
![]() |
Fig. 7 IE – YE vs. YE – HE colour–colour diagram for the Horsehead part of the FoV (not σ Orionis). The σ Orionis benchmarks are denoted with red circles. A reddening vector is shown as an arrow. |
In the text |
![]() |
Fig. 8 Combined very low mass Euclid–Gaia IMF of the σ Orionis cluster with linear fits in three different mass regions. |
In the text |
![]() |
Fig. 9 Number counts of Euclid sources in the catalogue of the σ Orionis region before and after filtering. The correspondence between VIS magnitudes and masses has been made using the ATMO 3 Myr isochrone. |
In the text |
![]() |
Fig. B.1 Cumulative distribution of spread model values for the Euclid JE band. The distribution of values in the low-reddening part corresponding to the σ Orionis cluster is sharper than in the high-reddening part. |
In the text |
![]() |
Fig. B.2 Euclid VIS FWHM (in pixels) versus IE apparent magnitude for objects confirmed by Euclid to be in the σ Orionis cluster sequence. Both the FWHM and IE photometry values come from the ERO catalogue. Four objects were found to have FWHM values larger than the mean value of cluster members (more than 3 σ significance) and hence are considered as possible binaries that deserve further scrutiny. Two of them are benchmark objects, namely (S Ori 52 and S Ori 60), one is a known object (S Ori 38), and the faintest one is a new discovery (Euclid object of interest G). |
In the text |
![]() |
Fig. B.3 Unresolved binary-corrected very low mass Euclid-Gaia IMF of the σ Orionis cluster with linear fits in three different mass regions. A 20% unresolved binary frequency of equal mass binaries has been assumed in the simulations. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.