# The enigma of Gaia18cjb: A possible rare hybrid of FUor and EXor
## properties*,**

Eleonora Fiorellino<sup>1,2,3,●</sup>, Péter Ábrahám<sup>2,3,4,5,●</sup>, Ágnes Kóspál<sup>2,3,4,6</sup>, Mária Kun<sup>2,3</sup>, Juan M. Alcalá<sup>1,●</sup>,
Alessio Caratti o Garatti<sup>1</sup>, Fernando Cruz-Sáenz de Miera<sup>2,3,7</sup>, David García-Álvarez<sup>8,9</sup>, Teresa Giannini<sup>10</sup>,
Sunkyung Park<sup>2,3</sup>, Michał Siwak<sup>2,3</sup>, Máté Szilágyi<sup>2,3,4,●</sup>, Elvira Covino<sup>1,●</sup>, Gabor Marton<sup>2,3</sup>, Zsófia Nagy<sup>2,3,●</sup>,
Brunella Nisini<sup>10</sup>, Zsófia Marianna Szabó<sup>2,3,11,12,●</sup>, Zsófia Bora<sup>2,3,4</sup>, Borbála Cseh<sup>2,3,13</sup>, Csilla Kalup<sup>2,3,4</sup>,
Máté Krezinger<sup>2,3,4</sup>, Levente Kriskovics<sup>2,3</sup>, Waldemar Ogłozna<sup>14</sup>, András Pál<sup>2,3</sup>, Ádám Sódor<sup>2,3</sup>, Eda Sonbas<sup>15,16</sup>,
Róbert Szakáts<sup>2,3</sup>, Krisztián Vida<sup>2,3,4</sup>, József Vinkó<sup>2,3</sup>, Lukasz Wyrzykowski<sup>17</sup>, and Pawel Zielinski<sup>18</sup>

*(Affiliations can be found after the references)*

Received 22 August 2023 / Accepted 12 March 2024

## ABSTRACT

**Context.** Gaia18cjb is one of the Gaia-alerted eruptive young star candidates that has been experiencing a slow and strong brightening during the last 13 years, similarly to some FU Orionis-type objects.

**Aims.** The aim of this work is to derive the young stellar nature of Gaia18cjb and determine its physical and accretion properties to classify its variability.

**Methods.** We conducted monitoring observations using multi-filter optical and near-infrared (NIR) photometry, as well as NIR spectroscopy. We present an analysis of pre-outburst and outburst optical and IR light curves, color-magnitude diagrams in different bands, the detection of NIR spectral lines, and estimates of both stellar and accretion parameters during the burst.

**Results.** The optical light curve shows an unusually long (over 8 years) brightening event of 5 mag in the last 13 years, before reaching a plateau indicating that the burst is still ongoing, suggesting a FU Orionis-like (FUor-like) nature. The same outburst is less strong in the IR light curves. The NIR spectra, obtained during the outburst, exhibit emission lines typical of highly accreting low-intermediate mass young stars with typical EX Lupi-type (EXor) features. The spectral index of Gaia18cjb SED classifies it as a Class I in the pre-burst stage and a flat-spectrum young stellar object (YSO) during the burst.

**Conclusions.** Gaia18cjb is an eruptive YSO that exhibits FUor-like photometric features (in terms of brightening amplitude and length of the burst) as well as EXor-like spectroscopic features and accretion rate. Its nature appears similar to that of V350 Cep and V1647 Ori, which have been classified as objects in between FUors and EXors.

**Key words.** accretion, accretion disks – techniques: imaging spectroscopy – stars: formation – stars: low-mass – stars: pre-main sequence – stars: protostars

## 1. Introduction

Low-mass (<2 M<sub>☉</sub>) young stellar objects (YSOs) exhibit photometric variability in the optical and infrared (IR) bands on timescales spanning from minutes to centuries (Carpenter et al. 2001; Megeath et al. 2012; Cody et al. 2014; Siwak et al. 2018). Their photometric variability can be caused by changes in accretion rate, varying line-of-sight extinction, or rotating accretion hot or cold spots (see Fischer et al. 2022, for a review).

Among variable sources, the so-called eruptive young stars (EYSs) show the largest brightness increases. Historically, we classified EYSs on the amplitude of the brightness change. EX Lupi-type objects (EXors, Herbig 1989) present an EXor burst of about 1–2.5 mag on timescales ranging from several weeks up to one year in the framework of the magnetospheric accretion scenario. Their spectra are very similar to Class II pre-

main sequence stars (PMSs) or classical T Tauri stars (CTTSs). FU Orionis-type objects (FUors, Herbig 1977), on the other hand, present a stronger brightening, of 2.5–6 mag, and it takes from months to years to reach the peak. Spectroscopically, FUors show absorption-line profiles, formed because the disc atmosphere is heated from the mid-plane, and their spectral types depend on the observed wavelengths: F–G type in optical and K–M type in near-infrared (NIR; Hartmann & Kenyon 1996; Audard et al. 2014; Connelley & Reipurth 2018; Fischer et al. 2022). Most importantly, the nature of the FUor outburst cannot be explained by magnetospheric accretion and it is likely to be related to other accretion mechanisms; for example, in boundary layer accretion, where the mass accretion rate ($M_{acc}$) is so high that shrinks the stellar magnetic field down to the stellar radius. No matter the mechanism and the order of magnitude, both EXor and FUor outbursts fuel the build-up of the protostellar mass, dramatically increasing the mass accretion rate from the disk to the central forming star. Such strong bursts are thought to be a possible explanation to the protostellar luminosity spread (Kenyon et al. 1990; Evans et al. 2009; Fischer et al. 2022), consisting of the fact that YSO’s observed luminosity is larger by about an order of magnitude than what was expected by the standard steady-state collapse mode by Shu (1977).

Recent works show that some outbursting YSOs have peculiar spectroscopic properties and light-curves that do not fit any

<footnote>* Spectra from various sources (LT, LBT, GTC) are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/ viz-bin/cat/J/A+A/686/A160</footnote>

<footnote>** Based on observations collected at the Large Binocular Telescope under LBT programme IT-2019B-008, the European Southern Observatory under ESO/NTT programmes 105.203T.001 and 105.203T.003, and the Gran Telescopio Canarias under GTC program GTC29-22B.</footnote>
of the aforementioned two classical categories (see Fischer et al. 2022, and references therein). Detailed studies of individual objects are important to understand the physical mechanism, evolution, and properties of very different EYSs, as well as to understand how eruptive accretion impacts the protostellar mass building process and the interplay with the steady magnetospheric accretion process.

In the last decade, also thanks to the whole-sky monitoring of the Gaia space observatory and its Gaia Photometric Science Alerts Program<fnref n="1" /> (Hodgkin et al. 2021), many EYSs have been discovered. These new EYSs have been classified as FUors, such as Gaia17bpi (Hillenbrand et al. 2018), Gaia18dvy (Szegedi-Elek et al. 2020), Gaia21bty (Siwak et al. 2023), and Gaia21elv (Nagy et al. 2023); EXors, like Gaia18dvz (Hodapp et al. 2019), Gaia19fct (Park et al. 2022), and Gaia20eae (Cruz-Sáenz de Miera et al. 2022; Ghosh et al. 2022); and in between, such as Gaia19ajj (Hillenbrand et al. 2019), and Gaia19bey (Hodapp et al. 2020).

Our study of Gaia18cjb enters in this context. An alert was issued for this source by Gaia Photometric Science Alerts Program on 2018 August 21. After the alert, we started monitoring this source with photometric and spectroscopic follow-up observations.

Based on information from Gaia DR2, the Wide-field Infrared Survey Explorer (WISE), and Planck measurements, Marton et al. (2019) evaluated that the probability of Gaia18cjb being a YSO is 83%. It is located at (RA, Dec) = (06h39m07.54s, 00d08m54.49s). Gaia DR3 measured a parallax of $p = 0.83 \pm 0.69$ mas, corresponding to a distance range between 660 pc and 7140 pc. Bailer-Jones et al. (2021) proposed a prior-dominated distance estimate for sources with large uncertainties on their parallax, finding Gaia18cjb can be located at a distance of $d = 1.03_{-0.19}^{+0.20}$ kpc. However, this determination uses a prior for the absolute G-band which is dominated by main sequence (MS) stars. The issue of Gaia18cjb's distance is further discussed in Sect. 4.1.

In this work, we present the data we collected (Sect. 2) and our analysis of these data (Sect. 3). We discuss the distance, the YSO nature of Gaia18cjb, and its outburst on the basis of our results in Sect. 4. Our conclusions are highlighted in Sect. 5.

## 2. Observations and data reduction

To probe the eruptive nature of Gaia18cjb, we collected a plethora of new photometric and spectroscopic observations, along with archival data from a range of telescopes. In this section, we describe the dataset and its reduction.

### 2.1. Photometry

Technical details of the optical, IR, and archival photometric data are listed in Table 1. In the following subsections, we describe in details the observations and the data reduction of the collected photometry.

#### 2.1.1. Optical photometry

We have been monitoring Gaia18cjb since October 2020 with an approximately monthly cadence using the 80 cm Ritchey-Chrétien (RC80) telescope at the Piszkéstö mountain station of Konkoly Observatory (Hungary). The telescope is equipped with an FLI PL230 CCD camera and Johnson BV and Sloan $g'r'i'$ fil-

ters. The aperture photometry for Gaia18cjb was obtained using 40 comparison stars located within a 12' × 12' box around our science target. We used an aperture radius of 5 pixels (2'75) and sky annulus between 20 and 40 pixels (11'' and 22''). The instrumental magnitudes were converted to the standard system using the APASS9 magnitudes (Henden et al. 2016) of the comparison stars and fitting a linear color term.

Gaia18cjb was also observed on three nights between 2021 January and February at the Mount Suhora Observatory of the Cracow Pedagogical University (Poland) with the 60 cm Carl-Zeiss telescope equipped with an Apogee Aspen-47 camera (Johnson BV and Sloan $g'r'i'$ filters). We obtained aperture photometry using the same aperture radius and sky annulus (in arcseconds) as for the RC80 data, and same 40 comparison stars for the photometric calibration.

Further observations were collected on six nights between 2020 November and 2021 February at Adiyaman University using ADYU60 Application and Research Center using ADYU60, a PlaneWave 60 cm f/6.5 corrected Dall-Kirkham Astrograph telescope. This telescope is equipped with an Andor iKon-M934 camera (Sloan $g'r'i'$ filters). Aperture photometry and calibration was done the same way as described above.

On 2021 May 9, we further observed Gaia18cjb with the 3.6m New Technology Telescope (NTT) located in La Silla Observatory (Chile). We used the ESO Faint Object Spectrograph and Camera (EFOSC, Arnaboldi et al. 2016, Bessell BVR filters). Aperture photometry was obtained as above. For the photometric calibration, we used five stars from the set above that were included in the smaller FoV of the detector without fitting a color term. We calculated their Bessell $R$ magnitudes by interpolating in their spectral energy distribution using their APASS9 BV$g'r'i'$ fluxes. In all the images that we acquired, the target appeared to be point-like and there were no issues with respect to the extended emission.

#### 2.1.2. IR photometry

On 2021 May 9, we observed Gaia18cjb with the ESO New Technology Telescope (NTT) IR spectrograph and imaging camera, namely, SofI (Arnaboldi et al. 2016). The data reduction of the photometry was done as prescribed by SOFI manual<fnref n="2" />. We first constructed a sky frame by taking the median of the five dither positions after scaling them to the same mean level. The flatfield frame was constructed by subtracting a lamp off image from a lamp on frame and correcting for the residual shade pattern. The resulting flat frames show vertical stripes on the left hand quadrants. To remove the stripes in the SW quadrant, we applied an empirical correction by taking the horizontal lines one-by-one from the affected upper-right quadrant. We subtracted the median value from each of the 512 pixel-long horizontal lines and added back the median of the remaining three quadrants of the whole image. Afterward, we divided the sky-subtracted and stripe-corrected frames by the flatfield. From the five dither positions, we constructed a mosaic per filter. Finally, we performed photometry on the mosaics. We used a 9 pixel radius aperture, with a sky ring between 7.8'' and 10.4''. For photometric calibration we used 30 2MASS sources with quality flag 'A'.

We also obtained JHK-band photometry of Gaia18cjb on 2022 October 10th and 26th using the 10.4m Gran Telescopio Canarias (GTC) installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias,

<footnote n="1">¹ https://gsaweb.ast.cam.ac.uk/alerts/home</footnote>

<footnote n="2">² https://www.eso.org/sci/facilities/lasilla/instruments/sofi/tools/reduction.html</footnote>
<table><caption>**Table 1.** Summary of photometric observations of Gaia18cjb.</caption><thead><tr><th>Telescope</th><th>Dates</th><th>Camera</th><th>Bands</th><th>Pixel scale</th><th>FoV</th></tr></thead><tbody><tr><td>RC80</td><td>2020–2023</td><td>FLI PL230 CCD</td><td>B,V,g',r',i'</td><td>0.55</td><td>18.8 × 18.8</td></tr><tr><td>Suhora</td><td>2021</td><td>Apogee Aspen-47</td><td>B,V,g',r',i'</td><td>1.116</td><td>19.0 × 19.0</td></tr><tr><td>ADYU60</td><td>2020–2021</td><td>Andor iKon-M934</td><td>g',r',i'</td><td>0.673</td><td>11.5 × 11.5</td></tr><tr><td>NTT</td><td>2021</td><td>EFOSC</td><td>B,V,R</td><td>0.24</td><td>4.1 × 4.1</td></tr><tr><td>NTT</td><td>2021</td><td>SOFI</td><td>J,H,K</td><td>0.273</td><td>4.5 × 4.5</td></tr><tr><td>GTC</td><td>2022</td><td>EMIR</td><td>Y,J,H,K</td><td>0.1945</td><td>6.67 × 6.67</td></tr><tr><td>ZTF</td><td>2017–2023</td><td>Cryogenic CCD</td><td>g',r'</td><td>1.0</td><td>47° × 47°</td></tr><tr><td>PanSTARRS</td><td>2010–2014</td><td>PS1 GPC1</td><td>g',r',i',z,Y</td><td>0.11</td><td>7° × 7°</td></tr><tr><td>2MASS</td><td>2000</td><td>3 NICMOS3</td><td>J,H,K<sub>s</sub></td><td>2.0</td><td>8.5 × 8.5</td></tr><tr><td>UKIDSS</td><td>2007–2010</td><td>WFCAM</td><td>J,H,K<sub>s</sub></td><td>0.4</td><td>0.21° × 0.21°</td></tr><tr><td>WISE</td><td>2010–2023</td><td>HAWAII 1RG</td><td>W1,W2</td><td>2.75"</td><td>47' × 47'</td></tr></tbody></table>

<table><caption>**Table 2.** Summary of spectroscopic observations of Gaia18cjb.</caption><thead><tr><th>Telescope</th><th>Dates</th><th>Spectrograph</th><th>Spectral range<br>µm</th><th>Resolution</th></tr></thead><tbody><tr><td>LT</td><td>2020 Nov</td><td>SPRAT</td><td>0.4020–0.7994</td><td>350</td></tr><tr><td>LBT</td><td>2021 Jan</td><td>LUCI</td><td>0.96–2.44</td><td>2000</td></tr><tr><td>NTT</td><td>2021 May</td><td>SOFI</td><td>0.95–2.52</td><td>1000</td></tr><tr><td>GTC</td><td>2022 Oct</td><td>EMIR</td><td>1.17–2.37</td><td>4000–5000</td></tr></tbody></table>

on the island of La Palma, with the wide field Espectrógrafo Multiobjeto Infra-Rojo (EMIR) Imager (Garzón et al. 2022). The images were reduced using dedicated python routines customised by the GTC staff for EMIR photometric data. The images were flat-fielded and the sky background was eliminated. The astrometric solution was also calculated. The final reduced image is the average of all the available images for each filter. We performed photometry on the reduced data using the same aperture selected for SOFI photometry (2.6"). For the photometric calibration, we used 30 2MASS sources with quality flag “A”.

### 2.1.3. Archival photometry

We complemented our monitoring data with public-domain photometry data. For the optical bands, we used the Gaia G band from the Gaia Photometric Science Alerts database, the Zwicky Transient Facility (ZTF, Masci et al. 2019) DR17 g'- and r'-band from the ZTF archive, and the Panoramic Survey Telescope & Rapid Response System (PanSTARRS).

We used the ZTF data with “catflags = 0”, that is, perfectly clean-extracted, to filter out bad-quality images. We downloaded the multi-epoch 2010–2014 Pan-STARRS images obtained in the g'r'i'zY bands<fnref n="3" />. We performed aperture photometry of Gaia18cjb using only the good quality images in 2.5" aperture. Then, to calibrate the results to absolute magnitudes, we used the magnitude zero points provided in the headers for each image and filter. While this approach may result in increased light losses during poor seeing conditions, we estimated that it does not exceed 0.2 mag in the g-band, and 0.05–0.1 mag for the other bands, which is sufficient for checking the light variability history before Gaia. As a consequence, we used a 0.2 dex of uncertainty for the g'-band photometry, and 0.1 dex for the other bands photometry. In addition, the smaller aperture

increased the number of useful images, as their certain fraction suffers from detector imperfections. We compared our results to APASS9 g'r'i' magnitudes of stars in the same field determined by Henden et al. (2016) finding that their measurements are reproduced with our method within 0.05 mag.

For IR wavelengths, we used JHK<sub>s</sub> photometric results from the Two Micron All Sky Survey (2MASS, Cutri et al. 2003) and downloaded processed JHK<sub>s</sub> images at two epochs from the UKIRT Infrared Deep Sky Survey (UKIDSS, Lucas et al. 2008). For UKIDSS images, to obtain results consistent with our SOFI and GTC IR photometry results (Sect. 2.1.2), we performed aperture photometry using the same aperture and sky radii (in arcseconds) as we had done for the SOFI and GTC data. Photometric calibration was performed via comparison with 15–25 2MASS sources, with the quality flag “A”. We also downloaded 3.35 µm (W1) and 4.60 µm (W2) photometry from the AllWISE Multiepoch Photometry Table and the NEOWISE-R Single Exposure (L1b) Source Table (Wright et al. 2010; Mainzer et al. 2011, 2014), available at the NASA/IPAC Infrared Science Archive<fnref n="4" />. After filtering out bad photometry, for instance, due to smeared images, images affected by the South Atlantic Anomaly or by scattered light from the moon, we averaged the magnitudes taken in each observing season (WISE scans the sky and takes a few frames of a certain area once every 6 months). Saturation correction was not necessary.

### 2.2. Spectroscopy

To probe the eruptive nature of Gaia18cjb, we observed it in the visible in one epoch and in the near-IR in three epochs, with different telescopes. Technical details on the spectroscopic data are listed in Table 2.

<footnote n="3">³ http://pslimages.stsci.edu/cgi-bin/pslcutouts</footnote>

<footnote n="4">⁴ https://irsa.ipac.caltech.edu/Missions/wise.html</footnote>
<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Hα line of Gaia18cjb from LT SPRAT spectrum.</figcaption></figure>

### 2.2.1. LT SPRAT

On 2020 November 11 we observed Gaia18cjb in the visible at low resolution (R = 350) with the SPectrograph for the Rapid Acquisition of Transients (SPRAT, Piascik et al. 2014), mounted on the Liverpool Telescope (LT)<fnref n="5" />. The spectrum was reduced and calibrated to absolute flux units by means of the dedicated SPRAT pipeline. The low resolution and high noise of the spectrum prevented us from using it for any quantitative analysis; however, it clearly shows the presence of the Hα line in emission (see Fig. 1).

### 2.2.2. LBT LUCI

We observed Gaia18cjb at the Large Binocular Telescope (LBT Observatory, Arizona) with LUCI1 and LUCI2 intruments in the zJ and HK filters on 2021 January 17. We reduced the data by using the SIPGI pipeline (Gargiulo et al. 2022), following the recipes indicated in the SIPGI documentation<fnref n="6" /> for the flat field correction, sky subtraction, wavelength calibration, and telluric correction. We then averaged the LUCI1 and LUCI2 spectra, to maximize the signal-to-noise ratio (S/N). Since we do not have contemporary zJHK photometry, we flux-calibrated the spectrum using SOFI continuum as both optical and NIR light curves suggest that the photometry does not vary much (e.g., Δ*m*<sub>r</sub> = 0.2 mag and Δ*m*<sub>w1</sub> = 0.1 mag) between January and May 2021 (see Fig. 2).

### 2.2.3. NTT SOFI

The blue part (0.95–1.64 µm) of the NTT/SOFI spectrum was observed on the 2021 May 5 and 8, while the red part (1.53–2.52 µm) only on the 8th. For the observed data, flat-fielding, bad pixel removal, sky subtraction, aperture tracing, and wavelength calibration were applied with the Image Reduction and Analysis Facility (IRAF, Tody 1993). For the wavelength calibration, a Xenon spectrum was used. The hydrogen lines in the telluric standard stars were removed by Gaussian fitting. The spectra were flux calibrated using the contemporary photometry (see Sect. 2.1). Technical details are reported in Table 2 and the spectra are presented in Sect. 3.

### 2.2.4. GTC EMIR

We obtained JHK-band spectroscopy of Gaia18cjb with the GTC/EMIR medium-resolution multi-object spectrograph (Espectrógrafo Multiobjeto Infra-Rojo, Garzón et al. 2022) in long slit mode on 2022 October 26. EMIR is equipped with a 2048 × 2048 Teledyne HAWAII-2 HgCdTe NIR-optimised chip with a pixel size of 0.2". The total exposure time on source was 1920 seconds per grism. One grism per each band (J, H and K) was used. A typical ABBA nodding pattern was applied. The seeing during the observations was FWHM ~ 1.2".

The spectra were reduced using several python routines customised by GTC staff for EMIR spectroscopic data. The sky background was first eliminated using consecutive A-B pairs. They were subsequently flat-fielded, calibrated in wavelength and combined to obtain the final spectrum. To correct for telluric absorption, we observed a telluric standard star with the same observing set up as the science target, right after the Gaia18cjb observations and at similar airmass. To apply the correction, we used a version of Xtellcor (Vacca et al. 2003) specifically modified to account for the atmospheric conditions of La Palma observatory (Ramos Almeida et al. 2009). We calibrated the flux scale of GTC spectra using the GTC contemporary photometry, see Sect. 2.1.2.

## 3. Results

Our photometric results, used for the optical and IR light curves (see Sect. 3.1), are shown in Appendix A. Uncertainties are computed as the quadratic sum of the uncertainty of the aperture photometry, the scatter of the individual exposures, and the uncertainty of the photometric calibration. Normalized and flux calibrated spectra are presented in Sect. 3.4 and Appendix B, respectively. They show emission lines typical of YSOs such as magnetospheric accretion features (HI), winds and outflows (H₂), and CO/NaI from the disk.

### 3.1. *Light curves*

Figure 2 shows the optical (top panel) and IR (bottom panel) light curve of Gaia18cjb.

#### 3.1.1. Optical light curve

The optical light curve shows the *Gaia* G-band brightening event of about 4.6 mag from 2015 February to 2020 May (black stars). Afterwards, the brightness still increased, albeit more slowly, by only 0.3 mag, up to the last epoch in March 2023. The overall brightening in the *Gaia* light curve is of about 4.9 mag. The outburst is still ongoing and there is no indication that the peak of the outburst has been reached. However, it looks as though Gaia18cjb reached a constant plateau in 2022-23. The photometric data, collected from other observatories after 2016, cover timescales shorter than *Gaia*. They confirm the brightening and the current “plateau” phase. Indeed, the r-band data (orange stars, filled circles, triangles, diamonds, and crosses, depending on the telescope) show a strong increase of brightness consistent with the *Gaia* data acquired with the G filter. The variability within 13 years is ~5.0 mag, a typical value for FUors variability (Fischer et al. 2022). However, the rise time (several years) is atypical for these kinds of objects, although there are FUors-like objects whose light curve has been increased for decades (e.g., V1515 Cyg, Hartmann & Kenyon 1996; Szabó et al. 2022). The ZTF g-band shows brightening of 2.8 mag within 4.75 years, and

<footnote n="5">⁵ ProgID: XOL20B01, PI: Pawel Zielinski.</footnote>

<footnote n="6">⁶ https://pandora.lambrate.inaf.it/docs/sipgi/</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Visible (top) and IR (bottom) light curves of Gaia18cjb. The red circle highlights the Gaia Alert trigger. The error bars smaller than the symbol size are not presented. Grey dashed lines correspond to the epochs for which we collected NIR spectroscopy.</figcaption></figure>

only of 0.5 mag after 2021 January (LUCI Epoch). NTT, RC80, Suhora, and ADYU60 light curves are in agreement within each other when the filter is the same and confirm the fact that since 2020 May the brightness of Gaia18cjb is still increasing but with a smaller amplitude than before (about 0.1 mag) in the *i*-, *r*-, *V*-, *g*-, and *B*- bands.

PanSTaRRS monitored Gaia18cjb prior to Gaia in the *g*, *r*, and *i* bands from 2010 to 2013. During these years the variation in the *r* band is only 0.65 mag, suggesting a quiescent phase before the ongoing outburst. The same behavior is also seen in other bands. We note that the last data points from the PanSTaRRS *r* band are brighter than the first data points of Gaia, suggesting the presence of some possible pre-burst variability. We also note that PanSTaRRS *r* and *g* variability is greater in a

specific epoch (i.e., May 2012) than in the three years overall. This can be a hint of ordinary small-amplitude accretion events as seen in CTTSs (see e.g., Fischer et al. 2022, and references therein).

### 3.1.2. IR light curve

We studied the variability also in the IR using WISE, UKIDSS, and PanSTaRRS archival data together with our acquired photometry (bottom panel of Fig. 2). The W1 channel at 3.4 µm (yellow filled dots) varies by 0.6 mag, while the variability is only 0.4 mag at 4.6 µm (W2 channel, blue dots) from 2010 March to 2022 October. The variability of both W1 and W2 channels reduced down to 0.1 since 2021 January. The W2 variability does
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Color-magnitude diagrams of Gaia18cjb based on data from Pan-STARRS data (yellow filled circles), ZTF, Konkoly, Mt Suhors, and WISE data before LUCI epoch (orange filled circles), and after LUCI epoch (blue filled circles). The arrows show the reddening vector, assuming extinction law from Cardelli et al. (1989) with $R_V = 3.1$ (black arrow) and $R_V = 5.5$ (red arrow), respectively.</figcaption></figure>

not seem to follow the variability we see in the optical bands and it hovers around a mean value between the pre-burst phase and the burst phase. We cannot confirm this trend in the W1 band since there were no W1 observations in the pre-burst phase. The IR light curve is similar to V1057 Cyg, one of the classical FUors (Szabó et al. 2021) post-1995, in displaying some kind of minor modulation.

The variability in the near-IR between UKIDSS (2007 April) and SOFI (2021 May) observations decreases with increasing wavelength: 2.70 mag, 2.14 mag, and 1.35 mag in J, H, and K bands, respectively. Interestingly, PanSTARRS light curves, observed from 2010 to 2014, invert the trend: in the Y band the variability is larger (0.89 mag) than in the z band (0.39 mag). This can be explained by the large variability in the first epoch (2010 February), which may be due to the presence of a small burst or to a sudden increase of the light curve, as is often the case before a strong FUor-like burst (Audard et al. 2014). However, this could also be the result of the different sampling for the PanSTARRS filters and, in turn, not necessarily associated with the intrinsic variability of the star.

We highlight that the photometric variability in the IR is moderate in general and it does not reflect the strong burst we observe in the visible. Indeed, the IR variability is as small as in typical non eruptive sources (see, e.g., Lorenzetti et al. 2012; Zsidi et al. 2022; Fiorellino et al. 2022). However, this result is puzzling since it suggests that the burst is restricted to the star and its immediate circumstellar environment and does not affect

the dusty disk at all. We discuss possible interpretations of this result in Sect. 4.

## 3.2. Color-magnitude diagrams

To investigate the photometric properties of Gaia18cjb in the visible and IR bands, we built color-magnitude diagrams (shown in Fig. 3). The top-left panel corresponds to the *g'* versus *g'* - *r'* diagram. It shows extinction related to color variations and a linear decreasing trend. This trend suggests that *g'*- and *r'*-band variations are correlated, but the *g'*-band variation is larger than *r'*-band variations (see light curves in Fig. 2). A similar trend is also present in the IR color-magnitude diagram in the bottom-right panel, where W1 and W1-W2 are strongly correlated following the direction of the reddening vectors. It appears that the source is redder when fainter, hence, the downward trend in the color-magnitude diagram. The top-right panel and the bottom-left panel present the V-band photometry versus *V* - *i'* and *V* - *r'*, respectively. In both cases, the distribution appears to be vertical, suggesting that the variability amplitude is similar in the V and *i'* and in the V and *r'* bands. In all the color-magnitude diagrams in Fig. 3 we note a temporal trend: the variation of the colors decreases with time as the source brightens. This produces a trend on the opposite direction of the reddening vectors. In general, the difference between the reddening law with $R_V = 3.1$ (black arrow) and $R_V = 5.5$ (red arrow) is small, and it is hard to distinguish which one of the two is followed better
<figure><img src="image_4.png" /></figure>

**Fig. 4.** NIR color-color diagram of Gaia18cjb. Purple and red circles are pre-outburst data from 2MASS and UKIDSS observations in 2000 and 2010, respectively, while the brown and pink dots correspond to SOFI and GTC photometry in the burst phase (2021 and 2022, respectively). The blue dashed line and light blue region represent the CTTS locus and its uncertainty (Meyer et al. 1997), and the tiny blue dotted line its extension toward larger color indices. The orange solid and green dotted-dashed lines correspond to the colors of ZAMS and GB stars, respectively. The black arrow shows the reddening vector for a source with 1 mag of extinction (Cardelli et al. 1989, $R_V = 3.1$).

from the data points. In the *g* versus *g* − *r* plot this difference is more evident and the black arrow better reproduces the variation of extinction from 2018 to 2020 (orange-blue filled dots), while the steeper slope of the red arrow better reproduces $A_V$ variations in the quiescent phase and in the latest epochs (from 2021 onward, yellow-blue filled dots). We discuss possible interpretations of this result in Sect. 4.

Figure 4 shows the *J* − *H* vs. *H* − *K* diagram, which can be used to classify a YSO and to estimate its extinction. We used UKIDSS and 2MASS archival data and the newly acquired SOFI and GTC data points. In the same figure, we plot the giant branch (GB) locus (green dot-dashed line), zero age main sequence (ZAMS, orange line), and classical T Tauri Stars (CTTS) locus (Meyer et al. 1997), where CTTS are not extincted ($A_V = 0$ mag, blue dashed line). By following the extinction vector by Cardelli et al. (1989) and $R_V = 3.1$ (black arrow), it is possible to compute the extinction of a CTTS. The 2MASS and UKIDSS data of Gaia18cjb lie within the CTTS locus, suggesting negligible extinction, while recent data points lie below the CTTS locus, making not possible to compute its extinction with this method and suggesting that Gaia18cjb is less evolved than a typical CTTS. However, in the case of GTC data, considering the error bars, it is possible that the Gaia18cjb colors lie on the CTTS locus, suggesting low extinction. This movement along the CTTS locus from redder to bluer is similar to what was found in Kóspál et al. (2011, see Fig. 7 therein) for HBC 722. In general, it is known that highly accreting stars become bluer in burst with respect to quiescence (see e.g., Lorenzetti et al. 2007, 2011; Hillenbrand et al. 2018; Szegedi-Elek et al. 2020), suggesting an accreting nature for Gaia18cjb.

### 3.3. Spectral energy distribution

From an observational point of view, YSOs are usually classified on the basis of their spectral energy distribution (SED, Lada & Wilking 1984; Lada et al. 1987; André 1995; Greene et al. 1994). Figure 5 shows the SED of Gaia18cjb, in the

<figure><img src="image_5.png" /></figure>

**Fig. 5.** Spectral energy distribution. Circles and stars are observations taken before and after 2020, respectively.

low brightening state (dots, before 2020), resembling the SED of an embedded object as a Class I YSOs, and in the outburst state (after 2020, star symbols), which is similar to a flat spectrum source. The data we used to build up the SEDs were observed at different periods and are listed in Appendix A.

We note that, based on the light curve, it is challenging to determine when the outburst begun. The Gaia alert was triggered in 2018, but the luminosity had already increased since 2015 (when first Gaia data are available), while the PanSTARRS monitoring between 2010 and 2014 does not provide definitive conclusions about the quiescent phase of this object, even if we can tentatively assume the average PanSTARRS magnitudes as values for the quiescent phase ($m_g = 21.8 \pm 0.4$, $m_r = 20.0 \pm 0.2$, $m_i = 19.1 \pm 0.1$, $m_z = 18.5 \pm 0.1$, and $m_Y = 18.2 \pm 0.2$). However, we distinguish here between the first part of the optical light curve, where the luminosity increases by 4.6 mag over about 10 years, and the second part, after the LUCI epoch, where the luminosity moderately increases by only 0.3 mag in three years (see the top panel of Fig. 2 and Sect. 3.1). We refer to the first one as the “brightening” or “low-state” phase and to the latter one as the “bursting” or “high-state” phase. We stress that this nomenclature is only useful to distinguish between the amplitude in the variability while in both the phases the burst is on, with the luminosity increasing.

We computed the spectral index $\alpha = \frac{d\log(F_{\lambda\alpha})}{d\log(\lambda)}$ between 2.16 µm and 22.1 µm using 2MASS and WISE data, respectively. For the pre-outburst phase, we obtain $\alpha = 0.55 \pm 0.28$, which classifies Gaia18cjb as a Class I protostar ($\alpha > 0.3$). We note that the spectral index is also compatible with a flat spectrum object ($-0.3 < \alpha < 0.3$) within the error. We conclude that Gaia18cjb is most likely a Class I YSOs, with a small chance that it is a flat spectrum. We also note that the classification we use computes the spectral index using $\lambda F_{\lambda}$ at 2 µm and 24 µm (not 22 µm as we do), but based on the SED shape, we do not expect a large difference between the flux at 22 µm and at 24 µm.

We also computed the spectral index using SOFI K-band and WISE W4 data to get an estimate of $\alpha$ for the eruptive phase. SOFI observations were taken during the eruptive phase, while
WISE observations are not, but since the low variability at long wavelengths and the lack of observations around 20 µm during the outburst, we used the same value as before (when considering the brightening SED). We obtained: $\alpha = 0.09 \pm 0.02$, which is compatible with a flat spectrum.

In both the low level and high level phases, the spectral index we compute does not describe Gaia18cjb as a Class II PMS young star. In the scenario where the empirical classification of YSOs based on the SED spectral index corresponds to an evolutionary classification, Gaia18cjb seems to be in between the protostellar phase, where the star is embedded in its envelope and the accretion is very high (Class 0 and I YSOs), and the pre-MS phase, where most of the stellar mass is accessed and the star is approaching the ZAMS (Class II and III YSOs). Estimates of Class 0 lifetimes are between 0.2 and 0.6 Myr, whereas based on *Spitzer* timescales, we see lifetimes of Class I and flat spectrum YSOs are 0.54 Myr and 0.40 Myr, respectively (Evans et al. 2009). Therefore, if the classification based on SED spectral index corresponds to evolutionary stages and our classification is correct, we can assume an age of about 1 Myr for Gaia18cjb.

### 3.4. Emission lines

The spectra of young solar-type stars show many features tracing the star-disk interaction as, for example, accretion and/or ejection. Our NIR spectra of Gaia18cjb covers several emission lines and, in particular, some accretion tracers, such as $B\gamma$ and $Pa\beta$. We describe here how we have measured the fluxes of these lines to be used for the analysis.

We measured the equivalent width ($W_{eq}$) and the observed line flux ($F_{obs}$) from the normalized and flux calibrated spectra, respectively, as follows. We fitted a linear curve to the local continuum in a wavelength range of $\Delta\lambda = 2$ nm, centered on the emission line wavelength $\lambda_0$. We slightly modified this range (when needed), taking the one most suitable for each line; for example, by avoiding other emission lines (if present) or telluric absorption lines. The line flux was determined by subtracting the local continuum from the spectra and integrating it over the line. We computed the noise of the line by multiplying the standard deviation of the local continuum (rms) for the wavelength element between two pixels, $\Delta\lambda$, and multiplying this by the square root of the number of pixels within the wavelength range ($N_{pix}$). We considered a line to be detected when its $S/N > 3$. For those lines that were detected in at least one epoch, we estimated the upper limits in the other epochs as three times the noise:

$$F_{\text{line}}^{\text{upp}} = 3 \times (\sqrt{N_{\text{pix}}} \times \text{rms} \times \Delta\lambda). \quad (1)$$

For LUCI spectra, the flux calibration with contemporary photometry was not possible. Therefore, we computed the observed flux by multiplying the equivalent width of each line by its local continuum in the SOFI spectrum. We stress that the acquisition of the LUCI spectrum was contemporaneous to that of the SOFI data. Equivalent widths and fluxes are reported in Tables 3 and 4, respectively.

**LBT spectrum: 2021 January.** Figure 6 reports the LBT normalized spectrum and the emission lines detected. This spectrum shows several accretion tracers such as $Pa\beta$, $Pay$, $Pa\delta$, $B\gamma$, $Br10$, $Br11$, $Br12$, $Br13$, $Br14$, and the NaI doublet. The $H_2$ line in emission suggests the presence of some shocked material in ejection, as winds. We also detected FeI in the emission.

<table><caption>**Table 3.** Equivalent width of Gaia18cjb emission lines for the three epochs we observed.</caption><thead><tr><th>Line</th><th>λ<br>[µm]</th><th>W<sub>eq</sub><sup>LUCI</sup></th><th>W<sub>eq</sub><sup>SOFI</sup></th><th>W<sub>eq</sub><sup>GTC</sup></th></tr></thead><tbody><tr><td>Paβ</td><td>1.281</td><td>18.41 ± 6.22</td><td>14.62 ± 1.66</td><td>21.99 ± 3.01</td></tr><tr><td>Pay</td><td>1.093</td><td>5.99 ± 0.96</td><td>7.20 ± 1.98</td><td>—</td></tr><tr><td>Paδ</td><td>1.004</td><td>4.58 ± 0.97</td><td>&lt;2.01</td><td>—</td></tr><tr><td>Bγ</td><td>2.166</td><td>11.90 ± 1.36</td><td>10.49 ± 3.08</td><td>10.51 ± 0.88</td></tr><tr><td>Br10</td><td>1.736</td><td>3.89 ± 0.66</td><td>&lt;3.93</td><td>3.12 ± 0.7</td></tr><tr><td>Br11</td><td>1.681</td><td>4.1 ± 1.12</td><td>&lt;6.99</td><td>2.06 ± 0.56</td></tr><tr><td>Br12</td><td>1.641</td><td>4.83 ± 1.33</td><td>&lt;4.77</td><td>1.61 ± 0.2</td></tr><tr><td>Br13</td><td>1.611</td><td>2.17 ± 0.56</td><td>&lt;8.94</td><td>1.14 ± 0.43</td></tr><tr><td>Br14</td><td>1.588</td><td>3.20 ± 1.12</td><td>&lt;1.60</td><td>2.77 ± 0.31</td></tr><tr><td>Br15</td><td>1.570</td><td>&lt;2.90</td><td>&lt;2.07</td><td>1.08 ± 0.12</td></tr><tr><td>Br16</td><td>1.556</td><td>&lt;3.81</td><td>1.16 ± 0.37</td><td>1.41 ± 0.16</td></tr><tr><td>Br17</td><td>1.544</td><td>&lt;2.28</td><td>&lt;2.31</td><td>1.04 ± 0.27</td></tr><tr><td>Br18</td><td>1.534</td><td>&lt;4.87</td><td>1.80 ± 0.42</td><td>&lt;2.73</td></tr><tr><td>Fe I</td><td>1.188</td><td>1.02 ± 0.15</td><td>&lt;5.39</td><td>2.3 ± 0.36</td></tr><tr><td>Fe I</td><td>1.197</td><td>2.84 ± 0.97</td><td>&lt;8.05</td><td>2.06 ± 0.65</td></tr><tr><td>Ca II</td><td>1.195</td><td>&lt;2.76</td><td>—</td><td>&lt;1.05</td></tr><tr><td>Si I</td><td>1.203</td><td>&lt;0.67</td><td>&lt;7.39</td><td>1.71 ± 0.64</td></tr><tr><td>Mg I</td><td>1.502</td><td>—</td><td>&lt;1.61</td><td>—</td></tr><tr><td>Mg I</td><td>1.504</td><td>—</td><td>3.45 ± 1.11</td><td>—</td></tr><tr><td>Mg I</td><td>1.576</td><td>&lt;2.72</td><td>&lt;10.94</td><td>0.96 ± 0.18</td></tr><tr><td>He I</td><td>2.033</td><td>—</td><td>&lt;0.85</td><td>1.73 ± 0.51</td></tr><tr><td>Ca I</td><td>2.101</td><td>&lt;1.35</td><td>4.20 ± 1.10</td><td>0.57 ± 0.22</td></tr><tr><td>Na I</td><td>2.206</td><td>2.08 ± 0.43</td><td>2.21 ± 0.49</td><td>&lt;1.35</td></tr><tr><td>Na I</td><td>2.208</td><td>2.18 ± 0.40</td><td>2.21 ± 0.49</td><td>0.92 ± 0.32</td></tr><tr><td>H<sub>2</sub></td><td>2.127</td><td>3.59 ± 0.71</td><td>2.71 ± 0.57</td><td>3.57 ± 0.45</td></tr></tbody></table>

Notes. The equivalent width unit is Å.

**SOFI spectrum: 2021 May.** Figure 7 shows the SOFI normalized spectrum. It confirms the presence of accretion ($Pa\beta$, $Pay$, $B\gamma$, $Br16$, $Br18$, NaI doublet) and ejection ($H_2$) tracers. Emission lines of MgI and CaI were also detected.

**GTC spectrum: 2022 October.** The GTC spectrum shows up to 9 lines of the Brackett series, but, due to the restricted wavelength range in the blue part of the spectrum, only the brightest line of the Paschen series ($Pa\beta$ see Fig. 8). Then, FeI, MgI, CaI, and NaI lines were detected, while $H_2$ emission was confirmed.

## 4. Discussion

In this section, we discuss evidences supporting that Gaia18cjb is a YSO and put its accretion properties into context with respect to other erupting YSOs.

### 4.1. On the distance of Gaia18cjb

Measuring the distance of an unknown object is a complex task. In particular, the case of Gaia18cjb presents its own set of challenges, given the large parallax uncertainty derived by Gaia. In this section, we discuss strengths and limitations of the photometric distance by Bailer-Jones et al. (2021).

The Gaia DR3 parallax for Gaia18cjb is $p = 0.83 \pm 0.69$ mas. With such a large uncertainty, the parallax method does not yield accurate results. Indeed, the main DR3 catalog does not provide a distance for sources whose parallax has a low S/N, including our target. For this reason, Bailer-Jones et al. (2021)
<table><caption>**Table 4.** Observed fluxes of the emission lines in the Gaia18cjb spectra.</caption><thead><tr><th>Line</th><th>λ<br>[µm]</th><th>F<sup>LUCI</sup><sub>obs</sub></th><th>F<sup>SOFI</sup><sub>obs</sub></th><th>F<sup>GTC</sup><sub>obs</sub></th></tr></thead><tbody><tr><td>Paβ</td><td>1.281</td><td>1.08 ± 0.07</td><td>1.96 ± 0.16</td><td>0.664 ± 0.046</td></tr><tr><td>Pay</td><td>1.093</td><td>0.73 ± 0.12</td><td>0.72 ± 0.10</td><td>—</td></tr><tr><td>Paδ</td><td>1.004</td><td>0.57 ± 0.12</td><td>&lt;0.25</td><td>—</td></tr><tr><td>Bry</td><td>2.166</td><td>0.88 ± 0.10</td><td>1.05 ± 0.11</td><td>0.680 ± 0.045</td></tr><tr><td>Br10</td><td>1.736</td><td>0.30 ± 0.07</td><td>&lt;0.25</td><td>0.078 ± 0.016</td></tr><tr><td>Br11</td><td>1.681</td><td>0.60 ± 0.08</td><td>&lt;0.29</td><td>0.051 ± 0.013</td></tr><tr><td>Br12</td><td>1.641</td><td>0.49 ± 0.18</td><td>&lt;0.29</td><td>0.041 ± 0.005</td></tr><tr><td>Br13</td><td>1.611</td><td>0.21 ± 0.05</td><td>&lt;0.43</td><td>0.029 ± 0.01</td></tr><tr><td>Br14</td><td>1.588</td><td>0.32 ± 0.11</td><td>&lt;0.40</td><td>0.072 ± 0.007</td></tr><tr><td>Br15</td><td>1.570</td><td>&lt;</td><td>—</td><td>0.029 ± 0.003</td></tr><tr><td>Br16</td><td>1.556</td><td>&lt;0.36</td><td>0.12 ± 0.04</td><td>0.037 ± 0.004</td></tr><tr><td>Br17</td><td>1.544</td><td>&lt;0.37</td><td>&lt;0.33</td><td>0.028 ± 0.007</td></tr><tr><td>Br18</td><td>1.534</td><td>&lt;0.37</td><td>0.20 ± 0.05</td><td>&lt;0.07</td></tr><tr><td>Fe I</td><td>1.188</td><td>0.08 ± 0.01</td><td>&lt;0.14</td><td>0.074 ± 0.01</td></tr><tr><td>Fe I</td><td>1.197</td><td>0.35 ± 0.12</td><td>&lt;0.84</td><td>0.068 ± 0.018</td></tr><tr><td>Ca II</td><td>1.195</td><td>&lt;0.12</td><td>&lt;0.80</td><td>0.032</td></tr><tr><td>Si I</td><td>1.203</td><td>&lt;0.25</td><td>&lt;0.80</td><td>0.055 ± 0.019</td></tr><tr><td>Mg I</td><td>1.502</td><td>—</td><td>&lt;0.56</td><td>—</td></tr><tr><td>Mg I</td><td>1.504</td><td>—</td><td>0.39 ± 0.13</td><td>—</td></tr><tr><td>Mg I</td><td>1.576</td><td>&lt;0.85</td><td>&lt;3.44</td><td>0.026 ± 0.005</td></tr><tr><td>He I</td><td>2.033</td><td>—</td><td>&lt;0.54</td><td>0.118 ± 0.032</td></tr><tr><td>Ca I</td><td>2.101</td><td>&lt;0.27</td><td>0.24 ± 0.05</td><td>0.037 ± 0.040</td></tr><tr><td>Na I</td><td>2.206</td><td>0.18 ± 0.04</td><td>0.19 ± 0.04</td><td>&lt;0.087</td></tr><tr><td>Na I</td><td>2.208</td><td>0.18 ± 0.03</td><td>0.19 ± 0.04</td><td>0.062 ± 0.021</td></tr><tr><td>H<sub>2</sub></td><td>2.127</td><td>0.25 ± 0.05</td><td>0.67 ± 0.18</td><td>0.232 ± 0.026</td></tr></tbody></table>

Notes. The observed flux unit is $10^{-14}$ erg s$^{-1}$ cm$^{-2}$. Observed fluxes are not corrected for the extinction. GTC and SOFI lines flux is computed from the flux calibrated spectra, while LUCI lines flux is computed from the $W_{eq}$ using the SOFI continuum for that line.

developed a probabilistic approach to estimate stellar distances. This method relies on the prior using a full range of populations in the mock catalogue, based on the Besancon Galaxy model. For each source, they computed two posterior probability distributions over distance: a geometric distance, based only on the prior and the parallax, and the photometric distance, based also on the G magnitude and the BP-RP color. Gaia18cjb's geometrical distance is $d_g = 3.20^{+0.66}_{-0.62}$ kpc, while the photometrical distance is $d_{pg} = 1.03^{+0.20}_{-0.19}$ kpc. The recommendation of Bailer-Jones et al. (2021) is to adopt the photometric distance for sources whose parallax has a low S/N. This is the case of Gaia18cjb.

The point of the distance inference in Bailer-Jones et al. (2021) is to provide a gradual transition from data-dominated to prior-dominated distance estimates. While the populations used by Bailer-Jones et al. (2021) contained YSOs, their results are not optimized for intrinsically red sources with high circumstellar extinction, such as YSOs. While a more appropriate photometric prior constructed solely for young stars would provide more accurate photometric distances, this is not the primary objective of our study. With the aim of remaining conservative, we discuss how our main results change, considering the spread in the distance suggested from Gaia parallax $d \sim 1/p = 660$ pc - 7140 pc. In the following, we discuss results assuming $d_{min} = 660$ pc, $d_{BJ} = 1030$ pc, and $d_{max} = 7140$ pc.

### 4.1.1. Interstellar extinction

To compute the reddening of Gaia18cjb due to the interstellar extinction, we used the 3D dust maps based on Gaia parallaxes and stellar photometry from Pan-STARRS 1 and 2MASS (Green et al. 2019). We ran the dustmaps python package<fnref n="7" />, using the most updated version Bayestar 2019, assuming $R_V = 3.1$, and the reddening law by Cardelli et al. (1989). We performed the same exercise using the SStructuring by Inversion the Local Interstellar Medium (STILISM, Capitanio et al. 2017; Lallement et al. 2018) with the online tool<fnref n="8" />. In both cases, we consider these values as lower limits for the extinction towards the source, because they are values averaged in a certain area and, most importantly, not sensitive to local circumstellar extinction. We repeat the calculation for $d_{min}$, $d_{BJ}$, and $d_{max}$. In Table 5 $A_V$ results are listed, ranging from 0.4 mag to 3.0 mag depending on the distance. We note that extinction estimates obtained with STILISM show slightly higher values, but compatible within the errors.

### 4.1.2. Bolometric luminosity

Using the fluxes reported in Appendix A, we computed the bolometric luminosity by integrating the SED fluxes ($F_\lambda$) with a dedicated python procedure, following the same method already used in literature (e.g., Antoniucci et al. 2008; Fiorellino et al. 2021, 2023). We integrated the SED from the shortest wavelength, interpolating with straight lines in the log $\lambda$-log $F_\lambda$ plane between the available SED points. We denote the value found with this method as $L_{bol}^{\text{pre}}$ to highlight that this is the bolometric luminosity corresponding to the pre-burst phase, obtained using the photometry observed before 2015 (filled circles in Fig. 5).

We added optical and NIR photometry acquired during the outburst (i.e., after 2015, filled stars). For the high-level phase, we computed the bolometric luminosity using Suhora, NTT/EFOSC, NTT/SOFI, and W1 photometry (blue, yellow, pink, and purple stars in the figure) instead of SDSS, PanSTARR, and 2MASS photometry (brown, green, and blue filled circles in the figure). As there is little or no variability in W2, we assumed that the source was constant within about 0.4 mag at 4.5 µm and we used the ALLWISE and AKARI points reward these wavelengths for the outburst SED, too. We name the value found with this method $L_{bol}^{\text{burst}}$, to highlight that this is the bolometric luminosity computed using the photometry observed during the outburst (see Fig. 5).

To build the SED and to measure the bolometric luminosity of Gaia18cjb we used the observed fluxes that were not corrected for extinction. We note that the longest wavelength data point we have is AKARI at 90 µm, and we assume a decrease in the emission as $1/\lambda^2$ after that. Even after applying this correction, we are aware that the bolometric luminosity derived so can be underestimated due to the lack of information at longer wavelengths. We also assume that beyond 3.6 µm the SED did not change during the burst.

It is desirable to correct the flux for extinction and compute the bolometric luminosity corrected for the extinction ($L_{bol}^{\text{corr}}$). Following Evans et al. (2009), since the spectral index suggest Gaia18cjb can be classified as a Class I/flat spectrum YSO, we used the interstellar $A_V$ to remove foreground extinction, but not local extinction from the surrounding envelope, which will be reradiated in the FIR. We used the mean value between the

<footnote n="7">⁷ http://argonaut.skymaps.info/usage</footnote>

<footnote n="8">⁸ https://stilism.obspm.fr/</footnote>
<figure><img src="image_6.png" /><figcaption>Fig. 6. Normalized LBT zJHK normalized spectrum. Emission lines reported in Table 3 are labeled.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 7. Normalized SOFI spectrum. Emission lines reported in Table 3 are labeled.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 8. Normalized GTC JHK spectrum. Emission lines reported in Table 3 are labeled.</figcaption></figure>
<table><caption>**Table 5.** Interstellar extinction and bolometric luminosity at different distances.</caption><thead><tr><th>Distance<br>pc</th><th>A<sup>B19</sup><sub>V</sub><br>mag</th><th>A<sup>STILIsM</sup><sub>V</sub><br>mag</th><th>A<sup>mean</sup><sub>V</sub><br>mag</th><th>L<sup>pre</sup><sub>bol</sub><br>L<sub>☉</sub></th><th>L<sup>pre_corr</sup><sub>bol</sub><br>L<sub>☉</sub></th><th>L<sup>burst</sup><sub>bol</sub><br>L<sub>☉</sub></th><th>L<sup>burst_corr</sup><sub>bol</sub><br>L<sub>☉</sub></th></tr></thead><tbody><tr><td>660</td><td>0.3 ± 0.1</td><td>0.5 ± 0.2</td><td>0.4 ± 0.2</td><td>1.4 ± 0.2</td><td>1.4 ± 0.2</td><td>3.0 ± 0.6</td><td>3.1 ± 0.6</td></tr><tr><td>1030</td><td>0.8 ± 0.2</td><td>1.0 ± 0.2</td><td>0.9 ± 0.3</td><td>3.4 ± 0.2</td><td>3.5 ± 0.2</td><td>8.0 ± 0.6</td><td>8.3 ± 0.6</td></tr><tr><td>7140</td><td>2.9 ± 0.4</td><td>3.1 ± 0.4</td><td>3.0 ± 0.6</td><td>162 ± 31</td><td>173 ± 31</td><td>353 ± 31</td><td>534 ± 31</td></tr></tbody></table>

values we obtain with the dustmaps and STILIsM softwares. Results are listed in Table 5.

The lowest possible distance and the BJ distance both results in pre-outburst luminosities (1.4–3.5 $L_{\odot}$) typical of a low-mass young protostar, and the outburst luminosities (3.1–8.3 $L_{\odot}$) are reasonable for modest accretion outbursts (cf. HBC 722 or some EXors). The luminosities obtained with the largest possible distance are in the several hundred $L_{\odot}$ range both in quiescence (173 $L_{\odot}$) and in outburst (534 $L_{\odot}$). While such a large outburst luminosity is possible in some of the most luminous and most highly accreting FUors (e.g., V1057 Cyg or FU Ori; Audard et al. 2014, and references therein), the quiescence luminosity is definitely too high for a low-mass YSO.

### 4.1.3. Stellar parameters and accretion rates

In this section, the stellar parameters (spectral type, extinction, luminosity, radius, and mass) and accretion properties (accretion luminosity and accretion rate) of Gaia18cjb are determined.

An estimate of the total extinction can be done by assuming that the accretion luminosity, $L_{\text{acc}}$, as derived from different accretion tracers, is the same. In order to derive $L_{\text{acc}}$, we used the Paβ and Bry lines and the empirical relationships between $L_{\text{acc}}$ and the line luminosity, $L_{\text{line}}$, by Alcalá et al. (2017). By minimizing the quantity $|L_{\text{acc}}(\text{Pa}\beta) - L_{\text{acc}}(\text{Bry})|$ it is possible to estimate the circumstellar extinction, $A_V$, as follows: the observed fluxes of the Paβ and Bry lines were dereddened with a series of $A_V$ values ranging from 0 to 200 mag in steps of 0.1 mag. The $A_V$ values minimizing the $|L_{\text{acc}}(\text{Pa}\beta) - L_{\text{acc}}(\text{Bry})|$ quantity are given in Table 6 and shown graphically in Fig. 9. For these, we determined values of $\log(L_{\text{acc}}/L_{\odot})$ between –0.07 and 2.63 dex, depending on the chosen distance (see Table 6). A robust value of $A_V$ for the eruptive phase can be then obtained by computing the average from the three individual estimates, the standard deviation being the error. This value ranges from $A_V = 10.5 \pm 2.7$ mag to $A_V = 15.0 \pm 0.3$ mag, depending on the distance (see Table 6).

Assuming that $L_{\text{bol}} = L_{\star} + L_{\text{acc}}$, we get a rough estimate of the stellar luminosity for Gaia18cjb by knowing the dereddened bolometric luminosity for the eruptive phase, $L_{\text{burst\_corr}}$, and the accretion luminosity, $L_{\text{acc}}$. We thus obtain $L_{\star} = 1.9 \pm 1.2 L_{\odot}$ for a $d_{\text{min}} = 660$ pc, $L_{\star} = 6.0 \pm 1.2 L_{\odot}$ for $d_{\text{BJ}} = 1030$ pc, and $L_{\star} = 107 \pm 62 L_{\odot}$ for $d_{\text{max}} = 7140$ pc.

Based on the spectral features, we discard the possibility for Gaia18cjb to be a high-mass star, since high excitation lines, such as He I are not detected. The stellar luminosity we obtain if assuming a distance of 7140 pc ($L_{\star} = 107 \pm 62 L_{\odot}$) is, therefore, not compatible with information contained in our spectra. Our spectra also discarded the possibility that the spectral type of Gaia18cjb could be later than K7. Indeed, M-type and later type spectra show photospheric absorption features such as the CO and VO bands. It is plausible to argue that the molecular bands could be veiled by the strong continuum emission; however, if that were the case, then the accretion luminosity should be extremely high – but the $W_{\text{eq}}$ value we obtained suggests this is not the case.

<table><caption>**Table 6.** Circumstellar extinction and accretion luminosity estimates for Gaia18cjb using empirical relations from Alcalá et al. (2017).</caption><thead><tr><th>Epoch</th><th>A<sub>V</sub><br>[mag]</th><th>log(L<sub>acc</sub>/L<sub>☉</sub>)<br>from Paβ</th><th>log(L<sub>acc</sub>/L<sub>☉</sub>)<br>from Bry</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center;">660 pc</td></tr><tr><td>LUCI</td><td>11.0 ± 0.3</td><td>-0.03 ± 0.21</td><td>-0.04 ± 0.24</td></tr><tr><td>SOFI</td><td>7.9 ± 0.3</td><td>-0.11 ± 0.17</td><td>-0.11 ± 0.24</td></tr><tr><td>GTC</td><td>12.5 ± 0.3</td><td>-0.08 ± 0.17</td><td>-0.08 ± 0.24</td></tr><tr><td>Mean</td><td>10.5 ± 2.7</td><td>-0.07 ± 0.11</td><td>-0.08 ± 0.14</td></tr><tr><td colspan="4" style="text-align: center;">1030 pc</td></tr><tr><td>LUCI</td><td>11.8 ± 0.3</td><td>0.39 ± 0.21</td><td>0.39 ± 0.31</td></tr><tr><td>SOFI</td><td>8.7 ± 0.3</td><td>0.34 ± 0.21</td><td>0.34 ± 0.32</td></tr><tr><td>GTC</td><td>13.3 ± 0.3</td><td>0.33 ± 0.12</td><td>0.33 ± 0.32</td></tr><tr><td>Mean</td><td>11.3 ± 0.3</td><td>0.35 ± 0.21</td><td>0.35 ± 0.21</td></tr><tr><td colspan="4" style="text-align: center;">7140 pc</td></tr><tr><td>LUCI</td><td>15.5 ± 0.3</td><td>2.68 ± 0.12</td><td>2.67 ± 0.18</td></tr><tr><td>SOFI</td><td>12.4 ± 0.3</td><td>2.59 ± 0.12</td><td>2.59 ± 0.18</td></tr><tr><td>GTC</td><td>17.0 ± 0.3</td><td>2.62 ± 0.12</td><td>2.62 ± 0.18</td></tr><tr><td>Mean</td><td>15.0 ± 0.3</td><td>2.63 ± 0.07</td><td>2.63 ± 0.10</td></tr></tbody></table>

Then, we used the evolutionary tracks by Siess et al. (2000), assuming an age of ~1 Myr (see Sect. 3.3) and considering only spectral type earlier than K7. We found that our results for $L_{\star}$ at 660 pc are compatible with a young star with stellar mass of $M_{\star} = (0.7-1.6) M_{\odot}$, a radius of $R_{\star} = (3.5-1.4) R_{\odot}$, and $T_{\text{eff}} = 4000-4600$ K, corresponding to SpT = K7–K3. The same exercise for $L_{\star}$ at 1030 pc leads to a young star with stellar mass of $M_{\star} = (0.7-2.2) M_{\odot}$, radius of $R_{\star} = (4.4-2.6) R_{\odot}$, and $T_{\text{eff}} = 4000-5080$ K, corresponding to SpT = K1–K7.

From the derived $L_{\text{acc}}$ and stellar parameters, we computed the mass accretion rate ($\dot{M}_{\text{acc}}$) in the bursting phase using the relation from Hartmann et al. (1998) in the framework of magnetospheric accretion scenario:

$$ \dot{M}_{\text{acc}} = \left(1 - \frac{R_{\star}}{R_{\text{in}}}\right)^{-1} \frac{L_{\text{acc}} R_{\star}}{GM_{\star}} \quad (2) $$

where $R_{\text{in}}$ is the inner-disk radius which we assume to be $R_{\text{in}} \sim 5R_{\star}$ as discussed by Hartmann et al. (2016, see their Eq. (3)) and for consistency with previous similar studies (e.g., Park et al. 2022; Cruz-Sáenz de Miera et al. 2022, 2023). We obtained $\dot{M}_{\text{acc}} = 1.8 \times 10^{-7} - 1.6 \times 10^{-7} M_{\odot} \text{yr}^{-1}$, for $d_{\text{min}}$; and $\dot{M}_{\text{acc}} = 1.0 - 5.3 \times 10^{-7} M_{\odot} \text{yr}^{-1}$, for $d_{\text{BJ}}$. Both the estimates are similar to the values found in highly accreting CTTSs and in EXors-like objects (e.g., Giannini et al. 2022).
## 4.2. On the YSO Nature of Gaia18cjb

Besides the photometric and spectroscopic features presented in the previous sections, another important aspect supporting the PMS nature of Gaia18cjb is the environment in which the object is located. We investigated the region nearby Gaia18cjb looking for its possible association with a star-forming region. Figure 10 shows the environment in the general direction of Gaia18cjb, with the nearby known star-forming regions highlighted. Gaia18cjb is located three degrees south-west of the Rosette star-forming region in the Monoceros constellation. Three degrees angular distance corresponds to ~50 pc at 1 kp. The distance to the Rosette complex is 1489 ± 32 pc (Mužić et al. 2022). Thus, it is unlikely that Gaia18cjb is connected with the Rosette star-forming region if its distance is $d \le d_{BJ}$. In contrast, if Gaia18cjb distance is similar to the Rosette star-forming region, then it might be part of the star-forming complex. The VDB 85 cluster is about ~1.5 degrees away from Gaia18cjb and its measured distance is about ~1.69 kpc (Cantat-Gaudin & Anders 2020), so, for the same reasoning applied to the Rosetta Complex, it is unlikely for Gaia18cjb to belong to this cluster. According to the 3D extinction map STILIsM, in correspondence of the TGU 1475 dark cloud position (l, b = 210.87°, -3.28°), there is one dark cloud at a distance of ~0.9–1.0 kpc, which we may suspect that corresponds to TGU 1475. The ellipse in Fig. 10 shows its nominal size from Dobashi et al. (2005). Although TGU 1475 and Gaia18cjb could have similar distances based on the assumption that $d \le d_{BJ}$, in looking at Fig. 10, we note that our target is outside the dark cloud. Therefore, we may conclude that Gaia18cjb is not a member of any of the known star-forming regions.

Recent investigations (Fratta et al. 2021; Prisinzano et al. 2022), based on $H\alpha$ surveys and on Gaia DR3, have shown evidence of populations of YSO candidates distributed in large volumes around the previously known star-forming regions. We consider whether Gaia18cjb might be part of such population in the general area of Monoceros. The sources showing $H\alpha$-excess are overplotted in Fig. 10, based on the areas selected by Fratta et al. (2021) and whose Gaia distance is between 0.95 kpc and 1.05 kpc, similar to the Gaia18cjb distance proposed by Bailer-Jones et al. (2021). The lack of $H\alpha$ emitting objects in the southern part is due to the latitude limit of the INT Photometric $H\alpha$ Survey of the north Galactic plane, IPHAS. We also distinguish in this plot objects having proper motions consistent, within $3\sigma$, to the proper motion of Gaia18cjb (orange dots) and those with inconsistent with proper motion (blue dots). These two samples are also plotted in the proper motion diagram in Fig. 11 (right panel). Gaia18cjb, shown as a red star symbol in the figure, is not associated with any source in this catalog.

Figure 11 (left panel) shows a color-magnitude diagram of the sample with $H\alpha$ excess, where the absolute magnitude computed from the Gaia G band as a function of the GBP-GRP color (blue and orange colors are as in Fig. 10). The black lines correspond to the ZAMS (solid) and to the 1 Myr isochrone (dashed, see Dell’Omodarme & Valle 2013; Randich et al. 2018; Tognelli et al. 2018, 2020). Models have been converted into the Gaia observed plane using the recent EDR3 filters passbands<fnref n="9" /> with the MARCS 2008 (Gustafsson et al. 2008) synthetic spectra for $T_{eff} \le 8000$ K and Castelli & Kurucz (2003) spectra for higher $T_{eff}$ (see e.g., Prisinzano et al. 2022). About 97% (91/94) of the $H\alpha$-excess sample lie above the ZAMS (blue dots), while only three lie below the ZAMS line. However, many of these objects, in particular, those piled-up close to the ZAMS, may

<figure><img src="image_9.png" /><figcaption>**Fig. 9.** Extinction estimates of Gaia18cjb computed from the empirical relations between the accretion luminosity and Bry and Paβ lines. Different colors represent $A_V$ estimates at different distances as described in the legend.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 10.** Background of the figure is the Planck 857 GHz map, 10 degrees in size, centred on Gaia18cjb, the red star symbol in the plot. Blue and orange dots are sources with $H\alpha$-excess from Fratta et al. (2021) whose distances are 0.95 pc < *d* < 1.05 pc. The orange dots are sources whose proper motions are within 3 mas yr⁻¹ from Gaia18cjb proper motions (see Fig. 11). The regions closest to Gaia18cjb are highlighted.</figcaption></figure>

be evolved $H\alpha$ emitting stars, such as novae and symbiotic stars (see Munari et al. 2022). We note, however, that symbiotic stars have giant companions and so are located near the giant sequence and quiescent novae resemble nova-like cataclysmic variables which lie between white dwarf and MS phases (see, e.g., Merc et al. 2020). Classical novae in outburst would probably be much brighter than the MS. We note that basically all the objects with proper motion compatible with those of Gaia18cjb (orange dots) lie well above the ZAMS hence, it is likely they are pre-MS (PMS) stars. However, wide-band spectroscopy of these objects is necessary to definitely establish their evolutionary status.

<footnote n="9">⁹ Gaia EDR3 passbands are available at https://www.cosmos.esa.int/web/gaia/edr3-passbands</footnote>
<figure><img src="image_11.png" /><figcaption>**Fig. 11.** Analysis on the young stellar nature of Gaia18cjb. **Left:** color-absolute magnitude diagram not corrected for the extinction. The red star symbol corresponds to Gaia18cjb. The black solid and dashed lines show the ZAMS and the 1Myr isochrone, respectively (Dell'Omodarme & Valle 2013). Dots are sources which show Hα-excess from Fratta et al. (2021) catalog whose distance is 0.95 kpc $< d < 1.05$ kpc. Orange dots are sources with proper motions consistent within $3\sigma$ the proper motion of Gaia18cjb (see right panel). The error bars smaller than the symbol size are not presented. The purple arrow represents an extinction vector of 1 mag. However, the $A_V$ vector depends on the color index and on the position in the diagram; see Prisinzano et al. (2022). **Right:** proper motions diagram of the Hα-excess sample. Symbols are the same as in the left panel.</figcaption></figure>

The red star symbol lying immediately above the ZAMS corresponds to the Gaia18cjb mean values of absolute magnitude and $G_{BP}-G_{RP}$ color provided by the Gaia archive: $G = 19.751 \pm 0.010$ mag, $G_{BP} = 20.65 \pm 0.77$ mag, and $G_{RP} = 18.16 \pm 0.77$ mag. These values, reported in Gaia EDR3 are averaged on the period 2014 July to 2017 May, when Gaia18cjb had already started brightening, but was still relatively faint; thus, it is more likely to display quiescent than outburst properties. We stress that the magnitudes of all these sources are not dereddened, implying that after correcting each source for its extinction, they will move up and left in the color-magnitude diagram. With respect to Gaia18cjb, we estimated an accretion luminosity in burst which is about 27% of the bolometric luminosity and $A_V \sim 8-12$ mag. With these values, the contribution of the extinction is larger than the contribution of the accretion, ending up moving Gaia18cjb up in the $M_{abs}$ versus $G_{BP}-G_{RP}$ plot. Therefore, we can speculate that these sources lie above the ZAMS, and are young stars. The study of the spatial distribution of YSO candidates within 1.5 kpc by Prisinzano et al. (2022), based on Gaia DR3, highlighted a population of distributed YSO candidates surrounding the previously known star-forming regions in the Monoceros region, but Gaia18cjb was not selected by their criterion, probably because it was too faint to be included in their catalog.

In summary, we selected 20 emission line stars (19 + Gaia18cjb) that: 1) lie above the ZAMS in the Gaia color-magnitude diagram; 2) are located at approximately the same distance ($0.95$ kpc $< d < 1.5$ kpc); and 3) have coherent proper motion. We therefore conclude that these likely YSOs, including Gaia18cjb might be part of the distributed population of the Monoceros complex. This results may indicate that the $d_{BJ}$ could be a good estimate of the real distance of Gaia18cjb. However, as discussed in Sect. 4.1, the photogeometric distance $d_{BJ}$ may be strongly affected by the prior for the absolute G-band mag, being dominated by the sample of MS stars.

Several pieces of evidence suggest that Gaia18cjb is a young star, based on: (i) the high probability of Gaia18cjb to be a YSO: 83%, computed by Marton et al. (2019); (ii) the similarity between Gaia18cjb SED and the SED of a young accreting object (Fig. 5); (iii) the spectral index which classifies Gaia18cjb as Class I YSO during the pre-burst phase and as a flat spectrum YSO during the burst (Sect. 3.3); (iv) the detection of CTTS typical features, as accretion and ejection tracers in the NIR spectra (see Sect. 3.4); and (v) the environment in which the object is projected on the sky. All these facts together lead us to conclude that Gaia18cjb is an accreting young star.

### 4.3. On the outburst of Gaia18cjb

The optical light curve (Fig. 2) presents a typical amplitude of a FUor-like object, on timescales that are compatible with EYOs, and longer than typical EXors and smaller than typical FUors. On the contrary, the WISE light curves in the MIR show no evidence of such FUor-like strong burst and the variability is smaller than expected for a FUor-like object, while it is compatible with the EXor-like YSOs or even CTTs variability.

If a burst is present, as the optical light curves suggest, it is not clear why it is less pronounced at NIR wavelengths and does not affect the SED at (and possibly beyond) the WISE band 2 (4.6 µm), since eruptive objects usually show similar variability at optical and IR wavelengths (see for example Gaia19fct, Park et al. 2022). Naively, we could conclude that the large optical and smaller amplitude in the IR for Gaia18cjb is due to variability in extinction. Looking at the color-magnitude diagrams in Fig. 3, we see that the extinction would only partially explain the variability of Gaia18cjb as well as its low-level and high-level colors.

For example, if the observed NIR variability ($\Delta W1 \sim 0.6$ mag, $\Delta m_K \sim 0.7$ mag, $\Delta m_H \sim 2.2$ mag, $\Delta m_J \sim 2.9$ mag) were caused only by extinction, this would correspond to a
visible amplitude of $\Delta m_V \sim 10-12$ mag, that is too large in comparison with the observed visual variability ($\Delta m_r \sim 5$ mag). This suggests that at least some of the observed blueing is intrinsic and not just linked to reddening. On the other hand, the accretion burst has to be relatively mild, as the continuum at 4.6 µm is barely or not affected at all. It is thus possible that the observed photometric variability of Gaia18cjb is a combination of enhanced accretion and reduced extinction, as inferred in Gaia19ajj (Hillenbrand et al. 2019). The reduced extinction might possibly rise from winds triggered by the enhanced accretion. These would partially clear the dust in the outflow cavities. This scenario is compatible with the H₂ line seen during the burst (see Sect. 3.4).

The burst weakness is also supported by the analysis of the spectra. The latter do not show typical FUors-like features such as the triangular shape in the *H* band or the undetected NIR HI and CO lines in absorption (Greene et al. 2008). On the contrary, accretion tracers are detected in emission in all the epochs we investigated, confirming this source is an accreting young star and strengthening the EXor nature of Gaia18cjb burst.

In Sect. 4.1.3, using empirical relations, we report that $L_{\text{acc}}$ ranges from 0.85 $L_{\odot}$ to 2.24 $L_{\odot}$, while $\dot{M}_{\text{acc}}$ ranges from $2.8 \times 10^{-8} M_{\odot} \text{yr}^{-1}$ to $1.0 \times 10^{-7} M_{\odot} \text{yr}^{-1}$ if $d_{\text{min}} \le d \le d_{\text{BJ}}$.

Calculating $L_{\text{acc}}$ using empirical relationships from Alcalá et al. (2017) is an indirect method that has shown success with CTTSs; however, some evidence has recently emerged that it may not be always applicable to eruptive sources, even if we expect magnetospheric accretion, as would be the case for EXors objects. For instance, a recent analysis of EX Lupi revealed that several emission lines indicated lower values of $L_{\text{acc}}$ compared to the slab model, the direct way to determine $L_{\text{acc}}$ in the framework of magnetospheric accretion (Cruz-Sáenz de Miera et al. 2023). They suggested that this discrepancy may arise from the challenge of identifying the continuum in a highly accreting objects. To investigate whether this holds true for Gaia18cjb, we would require wide-band spectroscopic data covering the Balmer continuum and Balmer jump to apply a slab model fit. Since we lack such spectroscopic data and, conversely, since we can only calculate $L_{\text{acc}}$ from the Paβ and Bry lines in the NIR, the estimate of $L_{\text{acc}}$ obtained from empirical relationships using these lines is currently the best estimate of $L_{\text{acc}}$. However, we emphasize that due to the likely overestimated value of Gaia18cjb's stellar luminosity and the uncertainty on the distance, there is a high probability that the accretion luminosity exceeds $L_{\text{acc}} = 2.24 L_{\odot}$. This needs to be tested through the direct measurement of $L_{\text{acc}}$ by studying the Balmer continuum excess emission and Balmer jump, using wide-band spectra.

Figure 12 shows the $L_{\text{acc}}$ vs. $L_{\star}$ (top panel) and $\dot{M}_{\text{acc}}$ vs. $M_{\star}$ (bottom panel) of Gaia18cjb at $d_{\text{min}} = 660$ pc and $d_{\text{BJ}} = 1.03$ kpc (red star symbols), compared to EX Lup (black stars, quiescent value from Cruz-Sáenz de Miera et al. 2023, and burst value from Fischer et al. 2022) and typical distributions for CTTSs in Lupus (grey, Alcalá et al. 2017), and Class I+flat spectrum protostars (pink, Fiorellino et al. 2023).

The top panel of Fig. 12 shows that Gaia18cjb is located in the Class I+flat spectrum YSO parameter space, right below the $L_{\text{acc}} = L_{\star}$ threshold (dashed blue line). Therefore, during the epochs we investigated with spectroscopy (i.e., during the burst), the ratio between the stellar and the accretion luminosity was more similar to a YSO more evolved than a protostar if $d_{\text{min}} \le d \le d_{\text{BJ}}$, and its accretion rate and stellar parameters are compatible with a Class I/flat spectrum object, suggesting Gaia18cjb's evolutionary stage lies between the protostellar and the PMS phases.

<figure><img src="image_12.png" /><figcaption>**Fig. 12.** Accretion parameters as a function of the stellar parameters. Top: $L_{\text{acc}} - L_{\star}$ diagram. The red and black star symbols correspond to Gaia18cjb and EX Lup, respectively. The blue dashed line represents the $L_{\text{acc}} = L_{\star}$ locus. The grey line is the best fit for the CTTSs of Lupus from Alcalá et al. (2017) and the grey region corresponds to the standard deviation of the t. Bottom: $\dot{M}_{\text{acc}} - \dot{M}_{\star}$ diagram. Star symbols are as in the top panel. The grey line is the best fit for the CTTSs of Lupus from Alcalá et al. (2017) and the grey region corresponds to the standard deviation of the t. The pink area corresponds to $\dot{M}_{\text{acc}}$ ranges for embedded Class I+flat spectrum YSOs from Fiorellino et al. (2023).</figcaption></figure>

This in agreement with the flat spectrum classification suggested by the spectral index (Sect. 3.3). We note, however, that the luminosity of Gaia18cjb is still compatible with that of typical CTTS, if the distance $d_{\text{min}} = 660$ pc is adopted.

Regarding the $\dot{M}_{\text{acc}}$ versus $M_{\star}$ plot, the bottom panel of Fig. 12 shows that Gaia18cjb accretes more than typical CTTSs of the same stellar mass. Gaia18cjb location in this plot is in agreement with bright Class I and flat spectrum sources and the mass accretion rate is similar to what derived for EX Lup during the burst. Indeed, while the optical light curve resemble a FUor, the spectroscopic features and the accretion rate suggest
an EXor-like nature of Gaia18cjb. Interestingly, there are at least
two young stars, V350 Cep and V1647 Ori, displaying photo-
metric characteristics of FUors in the visible bands and spec-
tral characteristics of EXors (Herbig 2008; Ibryamov et al. 2014;
Aspin & Reipurth 2009). In both cases, the conclusion was that
the sources were probably intermediate objects between FUors
and EXors.

The similarities between V350 Cep, V1647 Ori, and
Gaia18cjb, might suggest that also Gaia18cjb is an intermediate
object between FUors and EXors. However, there are a few
aspects to take into account. First, V350 Cep and V1647 Ori also
show PCyg Hα profile, and random fluctuations in brightness
with amplitudes of a few tenths of a magnitude in timescales
of several days that we do not see in Gaia18cjb. Second, the
current discussion about EYSs classification highlights that
the classification in very tight containers (e.g., FUor, EXor,
V1647 Or) is not sufficient. Gaia18cjb is another star suggesting
instead a spectrum of classifications in the EYSs. Last, our
results about accretion rates and stellar parameters should be
considered as rough estimates, because our analysis is based
only on NIR spectroscopy and because of the large uncertainty
on the distance. To provide accurate parameters and hence, a
robust characterization of Gaia18cjb, further observations and
analysis are needed:

- Continuing monitoring optical and IR light curves is crucial to follow the evolution of the Gaia18cjb burst. Also, *JHK* NIR IR monitoring is needed to study the light -curve variability in between the FUor-like visible light curve and the EXor-like IR light curves.

- By having a flux-calibrated intermediate and high-resolution spectrum from the UVB to the NIR (as ESO/VLT XShooter can provide), we will be able to study the Li I (λ6708 Å) absorption line to reinforce the YSOs nature of Gaia18cjb, to perform a direct spectral typing and determine in a more accurate way the stellar parameters, and hence, measure the accretion luminosity from the Balmer continuum excess emission. As a consequence, we will be able to provide robust estimates of the stellar and accretion parameters and of the mass accretion rate.

- Millimeter interferometry will provide estimates for disk mass and radius which are needed to prove the instability of the disk and believed to trigger the eruptive accretion.

- Accurate distance estimate, by using a pre-MS-dominated prior in the Bailer-Jones et al. (2021) prescription or a more accurate parallax estimate from GDR4.

**5. Conclusions**

The goal of this discovery paper is to determine the nature of the *Gaia*-alerted object Gaia18cjb. For this purpose, we collected all the available photometry in the archives and monitored Gaia18cjb in optical bands. NIR spectroscopy was also performed in three epochs between 2020 and 2022 and we studied the object’s variability, stellar and accretion parameters, and inspected the sky region nearby Gaia18cjb.

The large uncertainty on the *Gaia* parallax and the MS-
dominated prior of the Bailer-Jones et al. (2021) method, pre-
vented us from adopting a precise distance for this source.
Because of this, we cannot derive precise values, but only ranges
for the stellar and accretion parameters based on the compar-
ison between analytical results and spectral features. Compar-
ing the *Gaia* stellar parameters and accretion rates with those of
CTTSs, protostars, FUors, and other similar eruptive objects as
EX Lup, V350 Cep, V1647 Ori, and Gaia19ajj, our results sug-

gest Gaia18cjb is a young star experiencing eruptive accretion. It can be classified as in between FUors and EXors, namely: the length of the burst points to a FUor-like nature, whereas the strength of the burst and detected spectral features point to an EXor. However, further optical and NIR monitoring campaigns, based on at least one intermediate or high-resolution spectrum from UVB to NIR and millimeter interferometry, are needed to determine accurate stellar and accretion parameters and thereby provide robust conclusions on the Gaia18cjb variability classification.

Acknowledgements. We thank Dr. C. Bailer Jones for discussions on the parallax of the object and Dr. L. Prisinzano for providing the isochrones in the *Gaia* photometric system. We acknowledge Dr. András Ordasi for his contribution to this paper. This work has been supported by the project PRIN'INAF 2019 "Spectroscopically Tracing the Disk Dispersal Evolution (STRADE)" and by the INAF Large Grant 2022 "YSOs Outflow, Disks and Accretion (YODA)". This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme under grant agreement No 716155 (SACCRED). We acknowledge ESA Gaia, DPAC and the Photometric Science Alerts Team (https://gsaweb.ast.cam.ac.uk/alerts) This work is (partly) based on data obtained with the instrument EMIR, built by a Consortium led by the Instituto de Astrofísica de Canarias. EMIR was funded by GRANTECAN and the National Plan of Astronomy and Astrophysics of the Spanish Government. This work is (partly) based on data obtained at the Mount Suhora Observatory, Krakow Pedagogical University, Poland. This work has been supported by the Hungarian National Research, Development and Innovation Office grants OTKA K131508, OTKA K138962, the Élvonal KKP-143986 and KKP-137523 'SeismoLab' grants of the Hungarian Research, Development and Innovation Office (NKFIH). We acknowledge support from the ESA PRODEX contract nr. 4000132054. KV, LK, ZN, and GM are supported by the Bolyai János Research Scholarship of the Hungarian Academy of Sciences, KV is supported by the Bolyai+ grant ÚNKP-22-5-ELTE-1093. Authors acknowledge the financial support of the Austrian-Hungarian Action Foundation (112öu1). LK is supported by the Hungarian National Research, Development and Innovation Office grant PD-134784. G.M. acknowledges support from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 101004141. This work is supported by the Polish MNiSW grant DIR/WK/2018/12 and the European Union's Horizon 2020 research and innovation program under grant agreement No. 101004719 (OPTICON-RadioNet Pilot). Zs.M.Sz. acknowledges funding from a St Leonards scholarship from the University of St Andrews. For the purpose of open access, the author has applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising. BZs is supported by the ÚNKP-22-2 New National Excellence Program of the Ministry for Culture and Innovation from the source of the National Research, Development and Innovation Fund. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores University in the Spanish Observatory del Roque de los Muchachos of the Instituto de Astrofísica de Canarias with financial support from the UK Science and Technology Facilities Council. We acknowledge the Hungarian National Research, Development and Innovation Office grant OTKA FK 146023. FCSM received financial support from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (ERC Starting Grant "Chemtrip", grant agreement No 949278). This work was also supported by the NKFIH excellence grant TKP2021-NKTA-64.

**References**

Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20

André, P. 1995, Ap&SS, 224, 29

Antoniucci, S., Nisini, B., Giannini, T., & Lorenzetti, D. 2008, A&A, 479, 503

Arnaboldi, M., Delmotte, N., Geier, S., et al. 2016, in The Universe of Digital Sky Surveys, eds. N. R. Napolitano, G. Longo, M. Marconi, M. Paolillo, & E. Iodice, Astrophys. Space Sci. Proc., 42, 25

Aspin, C., & Reipurth, B. 2009, AJ, 138, 1137

Audard, M., Ábrahám, P., Dunham, M. M., et al. 2014, in Protostars and Planets VI, eds. H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning, 387

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Traedae, R. 2021, AJ, 161, 147

Cantat-Gaudin, T., & Anders, F. 2020, A&A, 633, A99

Capitanio, L., Lallement, R., Vergely, J. L., Elyajouri, M., & Monreal-Ibero, A. 2017, A&A, 606, A65

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

Carpenter, J. M., Hillenbrand, L. A., & Skrutskie, M. F. 2001, AJ, 121, 3160
Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, 210, A20

Cody, A. M., Stauffer, J., Baglin, A., et al. 2014, AJ, 147, 82

Connelley, M. S., & Reipurth, B. 2018, ApJ, 861, 145

Cruz-Sáenz de Miera, F., Kóspál, Á., Ábrahám, P., et al. 2022, ApJ, 927, 125

Cruz-Sáenz de Miera, F., Kóspál, Á., Ábrahám, P., et al. 2023, A&A, 678, A88

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, Vizier Online Data Catalog: II/246

Dell'Omodarme, M., & Valle, G. 2013, The R Journal, 5, 108

Dobashi, K., Uehara, H., Kandori, R., et al. 2005, Vizier Online Data Catalog: VII/244A

Evans, N. J., Dunham, M. M., Jorgensen, J. K., et al. 2009, Vizier Online Data Catalog: J/ApJS/181/321

Fiorellino, E., Manara, C. F., Nisini, B., et al. 2021, A&A, 650, A43

Fiorellino, E., Zsidi, G., Kóspál, Á., et al. 2022, ApJ, 938, 93

Fiorellino, E., Tychoniec, Ł., Cruz-Sáenz de Miera, F., et al. 2023, ApJ, 944, 135

Fischer, W. J., Hillenbrand, L. A., Herczeg, G. J., et al. 2022, ArXiv e-prints [arXiv:2203.11257]

Fratta, M., Scaringi, S., Drew, J. E., et al. 2021, Vizier Online Data Catalog: J/MNRAS/505/1135

Gargiulo, A., Fumana, M., Bisogni, S., et al. 2022, MNRAS, 514, 2902

Garzón, F., Balcells, M., Gallego, J., et al. 2022, A&A, 667, A107

Ghosh, A., Sharma, S., Ninan, J. P., et al. 2022, ApJ, 926, 68

Giannini, T., Giunta, A., Gangi, M., et al. 2022, ApJ, 929, 129

Green, G. M., Schlafly, E., Zucker, C., Speagle, J. S., & Finkbeiner, D. 2019, ApJ, 887, 93

Greene, T. P., Wilking, B. A., Andre, P., Young, E. T., & Lada, C. J. 1994, ApJ, 434, 614

Greene, T. P., Aspin, C., & Reipurth, B. 2008, AJ, 135, 1421

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hartmann, L., & Kenyon, S. J. 1996, ARA&A, 34, 207

Hartmann, L., Calvet, N., Gullbring, E., & D'Alessio, P. 1998, ApJ, 495, 385

Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135

Henden, A. A., Templeton, M., Terrell, D., et al. 2016, Vizier Online Data Catalog: II/246

Herbig, G. H. 1977, ApJ, 217, 693

Herbig, G. H. 1989, in European Southern Observatory Conference and Workshop Proceedings, 33, 233

Herbig, G. H. 2008, AJ, 135, 637

Hillenbrand, L. A., Contreras Peña, C., Morrell, S., et al. 2018, ApJ, 869, 146

Hillenbrand, L. A., Reipurth, B., Connelley, M., Cutri, R. M., & Isaacson, H. 2019, AJ, 158, 240

Hodapp, K. W., Reipurth, B., Pettersson, B., et al. 2019, AJ, 158, 241

Hodapp, K. W., Denneau, L., Tucker, M., et al. 2020, AJ, 160, 164

Hodgkin, S. T., Harrison, D. L., Breedt, E., et al. 2021, A&A, 652, A76

Ibryamov, S., Semkov, E., & Peneva, S. 2014, Res. Astron. Astrophys., 14, 1264

Kenyon, S. J., Hartmann, L. W., Strom, K. M., & Strom, S. E. 1990, AJ, 99, 869

Kóspál, Á., Ábrahám, P., Acosta-Pulido, J. A., et al. 2011, A&A, 527, A133

Lada, C. J. 1987, in Star Forming Regions, eds. M. Peimbert, & J. Jugaku, IAU Symp., 115, 1

Lada, C. J., & Wilking, B. A. 1984, ApJ, 287, 610

Lallement, R., Capitanio, L., Ruiz-Dern, L., et al. 2018, A&A, 616, A132

Lorenzetti, D., Giannini, T., Larionov, V. M., et al. 2007, ApJ, 665, 1182

Lorenzetti, D., Giannini, T., Larionov, V. M., et al. 2011, ApJ, 732, 69

Lorenzetti, D., Antoniucci, S., Giannini, T., et al. 2012, ApJ, 749, 188

Lucas, P. W., Hoare, M. G., Longmore, A., et al. 2008, MNRAS, 391, 136

Mainzer, A., Bauer, J., Grav, T., et al. 2011, ApJ, 731, 53

Mainzer, A., Bauer, J., Cutri, R. M., et al. 2014, ApJ, 792, 30

Marton, G., Ábrahám, P., Szegedi-Elek, E., et al. 2019, MNRAS, 487, 2522

Masci, F. J., Laher, R. R., Rusholme, B., et al. 2019, PASP, 131, 018003

Megeath, S. T., Gutermuth, R., Muzerolle, J., et al. 2012, AJ, 144, 192

Merc, J., Gális, R., Kára, J., Wolf, M., & Vrašták, M. 2020, MNRAS, 499, 2116

Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288

Munari, U., Alcalá, J. M., Frasca, A., et al. 2022, A&A, 661, A124

Mužić, K., Almendros-Abad, V., Bouy, H., et al. 2022, A&A, 668, A19

Nagy, Z., Park, S., Ábrahám, P., et al. 2023, MNRAS, 524, 3344

Park, S., Kóspál, Á., Ábrahám, P., et al. 2022, ApJ, 941, 165

Piascik, A. S., Steele, I. A., Bates, S. D., et al. 2014, in Ground-based and Airborne Instrumentation for Astronomy V, eds. S. K. Ramsay, I. S. McLean, & H. Takami, SPIE Conf. Ser., 9147, 91478H

Prisinzano, L., Damiani, F., Sciortino, S., et al. 2022, A&A, 664, A175

Ramos Almeida, C., Pérez García, A. M., & Acosta-Pulido, J. A. 2009, ApJ, 694, 1379

Randich, S., Tognelli, E., Jackson, R., et al. 2018, A&A, 612, A99

Shu, F. H. 1977, ApJ, 214, 488

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593

Siwak, M., Ogloza, W., Moffat, A. F. J., et al. 2018, MNRAS, 478, 758

Siwak, M., Hillenbrand, L. A., Kóspál, Á., et al. 2023, MNRAS, 524, 5548

Szabó, Z. M., Kóspál, Á., Ábrahám, P., et al. 2021, ApJ, 917, 80

Szabó, Z. M., Kóspál, Á., Ábrahám, P., et al. 2022, ApJ, 936, 64

Szegedi-Elek, E., Ábrahám, P., Wyrzykowski, Ł., et al. 2020, ApJ, 899, 130

Tody, D. 1993, in IRAF in the Nineties, eds. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, ASP Conf. Ser., 52, 173

Tognelli, E., Prada Moroni, P. G., & Degl'Innocenti, S. 2018, MNRAS, 476, 27

Tognelli, E., Prada Moroni, P. G., Degl'Innocenti, S., Salaris, M., & Cassisi, S. 2020, A&A, 638, A81

Vacca, W. D., Cushing, M. C., & Rayner, J. T. 2003, PASP, 115, 389

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

Zsidi, G., Fiorellino, E., Kóspál, Á., et al. 2022, ApJ, 941, 177

¹ INAF-Osservatorio Astronomico di Capodimonte, via Moariello 16, 80131 Napoli, Italy
e-mail: eleonora.fiorellino@inaf.it

² Konkoly Observatory, HUN-REN Research Centre for Astronomy and Earth Sciences, Konkoly-Thege Miklós út 15-17, 1121 Budapest, Hungary

³ CSFK, MTA Centre of Excellence, Konkoly-Thege Miklós út 15-17, 1121 Budapest, Hungary

⁴ ELTE Eötvös Loránd University, Institute of Physics and Astronomy, Pázmány Péter sétány 1/A, 1117 Budapest, Hungary

⁵ Institute for Astronomy (IfA), University of Vienna, Türkenschanzstrasse 17, 1180 Vienna, Austria

⁶ Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany

⁷ Institut de Recherche en Astrophysique et Planétologie, Université de Toulouse, UT3-PS, OMP, CNRS, 9 av. du Colonel Roche, 31028 Toulouse Cedex 4, France

⁸ Grantecan S.A., Centro de Astrofísica de La Palma, Cuesta de San José, 38712 Breña Baja, La Palma, Spain

⁹ Instituto de Astrofísica de Canarias, Avenida Vía Láctea, 38205 La Laguna, Tenerife, Spain

¹⁰ INAF-Osservatorio Astronomico di Roma, via di Frascati 33, 00078 Monte Porzio Catone, Italy

¹¹ Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany

¹² Scottish Universities Physics Alliance (SUPA), School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK

¹³ MTA-ELTE Lendület “Momentum” Milky Way Research Group, Budapest, Hungary

¹⁴ Mount Suhora Astronomical Observatory, Cracow Pedagogical University, ul. Podchorazych 2, 30-084 Kraków, Poland

¹⁵ Adiyaman University, Department of Physics, 02040 Adiyaman, Turkey

¹⁶ Astrophysics Application and Research Center, Adiyaman University, Adiyaman 02040, Turkey

¹⁷ Astronomical Observatory, University of Warsaw, Al. Ujazdowskie 4, 00-478 Warszawa, Poland

¹⁸ Institute of Astronomy, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland
## Appendix A: Photometry

Tables A.1 and A.2 collect the optical and NIR photometry of Gaia18cjb, respectively.

<table><caption>**Table A.1.** Optical photometry of Gaia18cjb we acquired and/or reduced. EFOSC *r* photometry was taken with Bessell *R* filter.</caption><thead><tr><th>Instr.</th><th>Date</th><th>HJD</th><th>B<br>mag</th><th>g<br>mag</th><th>V<br>mag</th><th>r<br>mag</th><th>i<br>mag</th></tr></thead><tbody><tr><td>PanSTaRRS</td><td>2011 Oct 25</td><td>2455860.13</td><td>—</td><td>—</td><td>—</td><td>20.2 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Oct 25</td><td>2455860.14</td><td>—</td><td>—</td><td>—</td><td>20.1 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Nov 01</td><td>2455566.95</td><td>—</td><td>22.1 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Nov 01</td><td>2455566.96</td><td>—</td><td>22.0 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Nov 11</td><td>2455876.06</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.2 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2011 Nov 11</td><td>2455876.07</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.4 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2012 Jan 04</td><td>2455930.94</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.1 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2012 Jan 04</td><td>2455930.96</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.3 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 14</td><td>2456215.08</td><td>—</td><td>21.8 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 14</td><td>2456215.08</td><td>—</td><td>21.7 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 14</td><td>2456215.10</td><td>—</td><td>22.4 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 14</td><td>2456215.10</td><td>—</td><td>21.7 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 19</td><td>2456220.09</td><td>—</td><td>—</td><td>—</td><td>20.1 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 19</td><td>2456220.10</td><td>—</td><td>—</td><td>—</td><td>20.0 ± 0.1</td><td>18.9 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 19</td><td>2456220.11</td><td>—</td><td>—</td><td>—</td><td>20.3 ± 0.1</td><td>19.0 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 08</td><td>2456240.09</td><td>—</td><td>—</td><td>—</td><td>19.7 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 08</td><td>2456240.09</td><td>—</td><td>—</td><td>—</td><td>20.1 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 08</td><td>2456240.10</td><td>—</td><td>—</td><td>—</td><td>19.8 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 08</td><td>2456240.10</td><td>—</td><td>—</td><td>—</td><td>20.1 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 09</td><td>2456241.05</td><td>—</td><td>—</td><td>—</td><td>20.1 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 09</td><td>2456241.06</td><td>—</td><td>—</td><td>—</td><td>20.0 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 11</td><td>2456243.04</td><td>—</td><td>21.6 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Nov 11</td><td>2456243.05</td><td>—</td><td>21.1 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Dec 23</td><td>2456284.96</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.0 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2013 Mar 05</td><td>2456356.78</td><td>—</td><td>21.9 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Mar 05</td><td>2456356.78</td><td>—</td><td>22.1 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Mar 05</td><td>2456356.79</td><td>—</td><td>22.6 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Mar 05</td><td>2456356.90</td><td>—</td><td>21.8 ± 0.2</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Dec 04</td><td>2456631.01</td><td>—</td><td>21.4 ± 0.2</td><td>—</td><td>19.9 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Dec 04</td><td>2456631.02</td><td>—</td><td>21.5 ± 0.2</td><td>—</td><td>19.8 ± 0.1</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Dec 23</td><td>2456649.99</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.0 ± 0.1</td></tr><tr><td>PanSTaRRS</td><td>2013 Dec 24</td><td>2456650.00</td><td>—</td><td>—</td><td>—</td><td>—</td><td>19.1 ± 0.1</td></tr><tr><td>Adiyaman</td><td>2020 Nov 13</td><td>2459167.49</td><td>—</td><td>17.688±0.052</td><td>—</td><td>16.227±0.026</td><td>15.564±0.060</td></tr><tr><td>Adiyaman</td><td>2020 Nov 27</td><td>2459180.51</td><td>—</td><td>17.617±0.065</td><td>—</td><td>16.149±0.025</td><td>15.495±0.054</td></tr><tr><td>Adiyaman</td><td>2020 Nov 29</td><td>2459183.39</td><td>—</td><td>17.425±0.082</td><td>—</td><td>16.162±0.031</td><td>15.491±0.045</td></tr><tr><td>Adiyaman</td><td>2021 Jan 10</td><td>2459225.28</td><td>—</td><td>17.498±0.073</td><td>—</td><td>16.141±0.036</td><td>15.504±0.047</td></tr><tr><td>Adiyaman</td><td>2021 Jan 24</td><td>2459239.45</td><td>—</td><td>17.575±0.086</td><td>—</td><td>16.213±0.033</td><td>15.503±0.052</td></tr><tr><td>Adiyaman</td><td>2021 Feb 10</td><td>2459256.40</td><td>—</td><td>17.613±0.058</td><td>—</td><td>16.174±0.025</td><td>15.503±0.056</td></tr><tr><td>EFOSC</td><td>2021 May 09</td><td>2459344.49</td><td>18.220±0.022</td><td>—</td><td>16.778±0.030</td><td>16.028±0.021</td><td>—</td></tr><tr><td>RC80</td><td>2020 Oct 22</td><td>2459145.56</td><td>18.32±0.18</td><td>17.597±0.035</td><td>16.903±0.032</td><td>16.242±0.014</td><td>15.516±0.012</td></tr><tr><td>RC80</td><td>2020 Nov 05</td><td>2459159.67</td><td>18.21±0.12</td><td>17.390±0.076</td><td>16.832±0.036</td><td>16.180±0.021</td><td>15.450±0.025</td></tr><tr><td>RC80</td><td>2020 Nov 22</td><td>2459175.65</td><td>18.50±0.15</td><td>17.562±0.038</td><td>16.880±0.041</td><td>16.213±0.012</td><td>15.470±0.011</td></tr><tr><td>RC80</td><td>2021 Feb 15</td><td>2459261.26</td><td>18.27±0.15</td><td>17.520±0.044</td><td>16.802±0.055</td><td>16.162±0.029</td><td>15.463±0.027</td></tr><tr><td>RC80</td><td>2021 Feb 18</td><td>2459264.27</td><td>18.17±0.11</td><td>17.511±0.079</td><td>16.830±0.051</td><td>16.167±0.031</td><td>15.461±0.030</td></tr><tr><td>RC80</td><td>2021 Feb 20</td><td>2459266.25</td><td>18.49±0.14</td><td>17.605±0.048</td><td>16.896±0.051</td><td>16.217±0.035</td><td>15.503±0.030</td></tr><tr><td>RC80</td><td>2021 Feb 23</td><td>2459269.25</td><td>17.80±0.25</td><td>17.29±0.15</td><td>16.54±0.14</td><td>16.20±0.13</td><td>15.296±0.088</td></tr><tr><td>RC80</td><td>2021 Feb 24</td><td>2459270.27</td><td>18.57±0.24</td><td>17.570±0.083</td><td>16.846±0.083</td><td>16.11±0.15</td><td>15.436±0.047</td></tr><tr><td>RC80</td><td>2021 Sep 08</td><td>2459465.61</td><td>18.30±0.13</td><td>17.530±0.053</td><td>16.785±0.032</td><td>16.137±0.016</td><td>15.410±0.018</td></tr><tr><td>RC80</td><td>2021 Oct 20</td><td>2459507.61</td><td>18.181±0.096</td><td>—</td><td>16.86±0.16</td><td>—</td><td>—</td></tr><tr><td>RC80</td><td>2021 Oct 23</td><td>2459510.52</td><td>—</td><td>17.548±0.052</td><td>—</td><td>16.128±0.021</td><td>15.415±0.016</td></tr><tr><td>RC80</td><td>2021 Nov 24</td><td>2459542.61</td><td>18.204±0.080</td><td>17.514±0.037</td><td>16.781±0.037</td><td>16.118±0.015</td><td>15.414±0.016</td></tr><tr><td>RC80</td><td>2022 Jan 07</td><td>2459586.59</td><td>18.19±0.10</td><td>17.464±0.047</td><td>16.688±0.036</td><td>16.013±0.017</td><td>15.319±0.017</td></tr><tr><td>RC80</td><td>2022 Feb 05</td><td>2459616.25</td><td>18.168±0.075</td><td>17.331±0.030</td><td>16.637±0.029</td><td>15.983±0.016</td><td>15.265±0.011</td></tr><tr><td>RC80</td><td>2022 Sep 24</td><td>2459846.63</td><td>18.16±0.11</td><td>17.313±0.034</td><td>16.636±0.030</td><td>15.982±0.026</td><td>15.290±0.016</td></tr><tr><td>RC80</td><td>2022 Oct 25</td><td>2459878.49</td><td>17.910±0.091</td><td>17.304±0.034</td><td>16.588±0.030</td><td>15.944±0.015</td><td>15.250±0.012</td></tr><tr><td>RC80</td><td>2022 Nov 25</td><td>2459908.62</td><td>18.172±0.085</td><td>17.478±0.032</td><td>16.731±0.029</td><td>16.090±0.015</td><td>15.405±0.016</td></tr><tr><td>RC80</td><td>2022 Dec 24</td><td>2459938.47</td><td>18.09±0.13</td><td>17.316±0.039</td><td>16.649±0.038</td><td>15.998±0.019</td><td>15.290±0.013</td></tr><tr><td>RC80</td><td>2023 Feb 06</td><td>2459982.34</td><td>18.06±0.19</td><td>17.363±0.094</td><td>16.599±0.056</td><td>15.933±0.023</td><td>15.220±0.017</td></tr><tr><td>RC80</td><td>2023 Mar 22</td><td>2460026.36</td><td>18.10±0.16</td><td>17.261±0.058</td><td>16.517±0.044</td><td>15.911±0.026</td><td>15.210±0.015</td></tr><tr><td>Suhora</td><td>2021 Dec 13</td><td>2459561.56</td><td>18.267±0.096</td><td>17.340±0.026</td><td>16.576±0.025</td><td>15.924±0.022</td><td>15.235±0.082</td></tr><tr><td>Suhora</td><td>2021 Jan 03</td><td>2459218.46</td><td>—</td><td>17.499±0.045</td><td>—</td><td>16.092±0.024</td><td>15.422±0.067</td></tr><tr><td>Suhora</td><td>2021 Jan 22</td><td>2459237.46</td><td>—</td><td>—</td><td>16.782±0.042</td><td>—</td><td>—</td></tr></tbody></table>
<table><caption>**Table A.2.** NIR photometry of Gaia18cjb we reduced.</caption><thead><tr><th>Instr.</th><th>Date</th><th>HJD</th><th>Y<br>mag</th><th>z<br>mag</th><th>J<br>mag</th><th>H<br>mag</th><th>K<br>mag</th></tr></thead><tbody><tr><td>UKIDSS</td><td>2007 Apr 15</td><td>2454205.5</td><td>—</td><td>—</td><td>—</td><td>—</td><td>12.620±0.040</td></tr><tr><td>UKIDSS</td><td>2010 Oct 24</td><td>2455493.5</td><td>—</td><td>—</td><td>16.210±0.040</td><td>14.860±0.050</td><td>13.260±0.040</td></tr><tr><td>PanSTaRRS</td><td>2010 Feb 15</td><td>2455242.8</td><td>18.7 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2010 Feb 15</td><td>2455242.8</td><td>18.3 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2010 Feb 15</td><td>2455242.8</td><td>18.6 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2010 Feb 15</td><td>2455242.8</td><td>18.1 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2010 Oct 22</td><td>2455492.1</td><td>18.2 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2010 Oct 22</td><td>2455492.1</td><td>18.1 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Sep 18</td><td>2455823.1</td><td>18.2 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Sep 18</td><td>2455823.1</td><td>17.9 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Oct 22</td><td>2455856.8</td><td>18.1 ± 0.1</td><td>18.6 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Oct 22</td><td>2455256.8</td><td>—</td><td>18.4 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Oct 27</td><td>2455862.1</td><td>—</td><td>18.5 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2011 Oct 27</td><td>2455862.1</td><td>—</td><td>18.6 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Feb 27</td><td>2455984.7</td><td>18.0 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Feb 27</td><td>2455984.7</td><td>18.1 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Feb 28</td><td>2455485.1</td><td>—</td><td>18.4 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Feb 28</td><td>2455485.1</td><td>—</td><td>18.5 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Mar 20</td><td>2456006.8</td><td>—</td><td>18.5 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Mar 20</td><td>2456006.8</td><td>—</td><td>18.6 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2012 Oct 22</td><td>2456223.1</td><td>—</td><td>18.3 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2013 Oct 20</td><td>2456586.1</td><td>18.0 ± 0.1</td><td>18.2 ± 0.1</td><td>—</td><td>—</td><td>—</td></tr><tr><td>PanSTaRRS</td><td>2014 Oct 08</td><td>2456939.1</td><td>18.0 ± 0.1</td><td>—</td><td>—</td><td>—</td><td>—</td></tr><tr><td>SOFI</td><td>2021 May 08</td><td>2459343.0</td><td>—</td><td>—</td><td>13.505±0.058</td><td>12.720±0.035</td><td>11.911±0.042</td></tr><tr><td>GTC</td><td>2022 Oct 10</td><td>2459862.5</td><td>—</td><td>—</td><td>13.354±0.080</td><td>12.672±0.073</td><td>11.952±0.106</td></tr><tr><td>GTC</td><td>2022 Oct 26</td><td>2459878.6</td><td>—</td><td>—</td><td>13.396±0.124</td><td>12.675±0.067</td><td>12.059±0.160</td></tr></tbody></table>

<table><caption>**Table A.3.** Photometry used for the SED of Gaia18cjb.</caption><thead><tr><th>Instr.</th><th>Year</th><th>λ<br>μm</th><th>Flux<br>mJy</th></tr></thead><tbody><tr><td>WISE</td><td>2013</td><td>11.55</td><td>93.1 ± 1.3</td></tr><tr><td>WISE</td><td>2013</td><td>22.09</td><td>145.0 ± 3.0</td></tr><tr><td>PanSTaRRS</td><td>2012</td><td>0.747</td><td>0.06580 ± 0.0014</td></tr><tr><td>PanSTaRRS</td><td>2012</td><td>0.865</td><td>0.1300 ± 0.0020</td></tr><tr><td>PanSTaRRS</td><td>2012</td><td>0.959</td><td>0.1620 ± 0.0040</td></tr><tr><td>WISE</td><td>2011</td><td>3.350</td><td>14.40 ± 0.30</td></tr><tr><td>WISE</td><td>2011</td><td>4.600</td><td>37.50 ± 0.70</td></tr><tr><td>AKARI</td><td>2010</td><td>8.610</td><td>114 ± 16</td></tr><tr><td>AKARI</td><td>2010</td><td>64.99</td><td>461 ± 76</td></tr><tr><td>AKARI</td><td>2010</td><td>90.00</td><td>506 ± 81</td></tr><tr><td>2MASS</td><td>2000</td><td>1.239</td><td>0.665 ± 0.052</td></tr><tr><td>2MASS</td><td>2000</td><td>1.649</td><td>1.430 ± 0.080</td></tr><tr><td>2MASS</td><td>2000</td><td>2.163</td><td>3.95 ± 0.15</td></tr><tr><td>SDSS</td><td>1998</td><td>0.477</td><td>0.00847 ± 0.00037</td></tr><tr><td>SDSS</td><td>1998</td><td>0.481</td><td>0.00849 ± 0.00037</td></tr><tr><td>SDSS</td><td>1998</td><td>0.612</td><td>0.03260 ± 0.00040</td></tr><tr><td>EFOSC</td><td>2021</td><td>0.44</td><td>1834 ± 38</td></tr><tr><td>Suhora</td><td>2021</td><td>0.52</td><td>0.373 ± 0.016</td></tr><tr><td>EFOSC</td><td>2021</td><td>0.5476</td><td>6331 ± 177</td></tr><tr><td>EFOSC</td><td>2021</td><td>0.6431</td><td>14823 ± 285</td></tr><tr><td>Suhora</td><td>2021</td><td>0.67</td><td>1.638 ± 0.036</td></tr><tr><td>Suhora</td><td>2021</td><td>0.79</td><td>3.226 ± 0.198</td></tr><tr><td>SOFI</td><td>2021</td><td>1.247</td><td>57.8 ± 3.1</td></tr><tr><td>SOFI</td><td>2021</td><td>1.653</td><td>75.7 ± 2.4</td></tr><tr><td>SOFI</td><td>2021</td><td>2.162</td><td>103.2 ± 4.0</td></tr></tbody></table>
**Appendix B: Spectra**

Figure B.1 shows the NIR flux calibrated spectra of Gaia18cjb.

<figure><img src="image_13.png" /><figcaption>Fig. B.1. Flux-calibrated NIR spectra of Gaia18cjb.</figcaption></figure>