A new take on the low-mass brown dwarf companions on wide orbits in Upper-Scorpius

Context. The Upper-Scorpius association (5–11 Myr) contains a unique population of low-mass ( M ≤ 30


Introduction
The first brown dwarfs (BDs) were contemporaneously discovered with the first exoplanets at the end of the last millennium (Mayor & Queloz 1995;Nakajima et al. 1995;Rebolo et al. 1995).Since then, thousands of BDs have been detected, and Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO program 093.C-0769.studied in isolation (Cushing et al. 2005;Kirkpatrick 2005) in the field or as wide-or short-period companions to nearby stars (Sahlmann et al. 2011;Allers 2012).Bridging the gap between planets and stars, BDs are too light by definition to burn hydrogen but are massive enough to burn deuterium (Burrows et al. 1997).Despite two decades of intensive study of BDs, many fundamental questions remain unanswered regarding their formation and evolution processes, their physical and atmosphere properties, and their connection to stars and planets.Multiple stellar-like formation pathways have been proposed for A&A 633, A124 (2020) these objects: (i) turbulent fragmentation of molecular clouds (Padoan & Nordlund 2004), (ii) premature ejection of protostellar embryos (Bate et al. 2002), (iii) photo-erosion of prestellar cores (Hester et al. 1996), and (iv) disk instability (e.g., Stamatellos & Whitworth 2009).Studies of young star-forming regions are currently ongoing to identify differences in binary statistics linked to these different processes (e.g., Marks et al. 2017;Thies et al. 2015).Alternatively, we know from observations and theory that planetary formation mechanisms like core accretion are probably forming very massive giant planets, and populating the mass distribution of substellar companions up to M ≤ 35 M Jup (e.g., Mordasini et al. 2009Mordasini et al. , 2012)).There is therefore little doubt that both stellar and planetary formation mechanisms overlap in mass distribution, and an interesting question is to investigate whether stellar and planetary mechanisms might lead to different atmospheric properties that could be traced through observations.Core accretion might indeed lead to an overabundance of heavy elements in the atmosphere.Our ability to identify such chemical imprinting is very challenging, is beyond observational limitations, and is directly connected to our understanding and modeling of the physics of BD and exoplanet atmospheres.Confronting the latest predictions of substellar atmosphere models with high-quality optical and infrared spectra of young brown dwarfs is a key step toward this goal.
Evolutionary models predict how BDs contract and cool down with time, and how their fundamental parameters such as effective temperature, surface-gravity, radius, and luminosity evolve (Baraffe et al. 2003).The contraction leads to a decrease of the radius and an increase of the surface-gravity.This evolution impacts the bolometric luminosity and the spectrum with a modification of the pseudo-continuum and the appearance of atomic and molecular absorption lines at different evolutionary stages.Spectral morphology has been used for years to extend the old stellar spectral classification into the substellar one, from M-dwarfs to L-, T-, and Y-dwarfs reaching effective temperatures as cold as T eff ≤ 450 K (Kirkpatrick et al. 2012;Kopytova et al. 2014).The fine characterization of the BD atmospheres therefore improves our global understanding of their physical properties, and also of their formation and evolution.
There are currently two approaches for the spectral characterization of BDs.The first is empirical and based on comparison with libraries of known young and old BDs (e.g., Allers & Liu 2013;Bonnefoy et al. 2014;Bayo et al. 2017), and is tightly connected and limited by the size and the diversity of these libraries.A complementary alternative is the comparison of the observed spectra to the recent models of substellar atmospheres (e.g., Allard et al. 2001;Helling & Woitke 2006;Barman et al. 2011;Madhusudhan et al. 2011;Morley et al. 2012;Tremblin et al. 2017;Charnay et al. 2018) of BDs and giant planets.This offers the advantage that the physical parameters for a given model can be derived independently of other observations, and that the influence of new physical ingredients can be tested (e.g., new atomic and molecular line opacities, presence and properties of clouds, nonequilibrium chemistry, thermochemical instability, etc.).Nevertheless, this latter approach is hindered by our limited observational knowledge and constraints to disentangle the relative importance of the tested physical parameters, eventually leading to significant degeneracies and systematic errors when using different families of atmosphere models or simply different sets of physical parameters.Both approaches remain therefore highly complementary today.
The detection of forming companions is a new observational window on the initial conditions of planetary systems (formation zone, timescales, and modes; accretion physics; e.g., Mordasini et al. 2017).Thus far, H α (656.3 nm) and Pa β (1282.2nm) lines have been detected in the spectra of eight companions with masses below 30 M Jup and ages in the 1-10 Myr range (TWA 5 B, GQ Lup B, CT Cha B, USco CTIO 108 B, DH Tau B, GSC 06214-00210 B, SR 12 C, PDS 70 B, Neuhäuser et al. 2000;Seifahrt et al. 2007;Schmidt et al. 2008;Béjar et al. 2008;Bonnefoy et al. 2014;Bowler et al. 2011;Santamaría-Miranda et al. 2018;Wagner et al. 2018;Zhou et al. 2014;Wu et al. 2017b).These lines are known tracers of active accretion and of substellar chromospheric activity.Submillimeter observations of these objects have so far failed to reveal excess emission from the expected mass reservoir (circumplanetary disk) surrounding these objects (Dai et al. 2010;Bowler et al. 2015;MacGregor et al. 2017;Wolff et al. 2017;Wu et al. 2017a,b) and only one companion, GSC 06214-00210 B, shows clear excess emission at near-infrared (NIR) wavelengths (Bowler et al. 2011).This might indicate that these companions bear very compact and optically thick disks (Wu et al. 2017b).
The Upper-Scorpius subgroup (hereafter Upper-Sco) in the Scorpius-Centaurus OB association contains one of the nearest (d = 146 ± 3 pc; de Bruijne 1999;de Zeeuw et al. 1999;Galli et al. 2018) and richest populations of young stars and substellar objects (e.g.Ardila et al. 2000;Lodieu et al. 2007Lodieu et al. , 2018;;Luhman et al. 2018) down to the planetary-mass range.The extinction is low in this region (A V ≤ 2 mag; Walter et al. 1994;Lallement et al. 2019).At an estimated age of 5-11 Myr (Pecaut et al. 2012;David et al. 2019), stars harbor primordial, transitional, and debris disks (Luhman & Mamajek 2012;Esplin et al. 2018) suggesting planet or BD formation within disks at different completion levels.Upper-Sco also contains a large sample of low-mass BDs (M < 30 M Jup ) and planetary mass companions identified with deep-imaging and high-contrast imaging techniques.These companions are found over a wide range of projected separations (∼300-3400 au; Béjar et al. 2008;Lafrenière et al. 2008Lafrenière et al. , 2010Lafrenière et al. , 2011;;Ireland et al. 2011;Aller et al. 2013) from M 7-to B6-type stars.Some of the free-floating lowmass BDs and companions harbor disks (e.g., Bowler et al. 2011;Dawson et al. 2013) and are actively accreting (Herczeg et al. 2009;Bowler et al. 2011;Lodieu et al. 2018).This unique population of low-mass BDs and planetary mass objects with various configurations (as companions, binaries, free-floating) and likely diverse origins represents a unique test bed for planet and BD formation models.
In this paper, we present a new study of three young BD companions to the stars HIP 78530, HIP 77900, and USco 161031.9-16191305(hereafter USco 1610-1913), members of the Upper-Scorpius association.We obtain VLT/X-shooter 0.3-2.5 µm spectra of these objects, and use them to characterize their physical properties.The data notably extend previous analysis of these objects to the optical at medium resolving powers, thus enabling investigation of emission lines related to accretion, and testing of the atmospheric models of young BDs.We provide an up-todate description of our targets in Sect. 2. The observations and the data reduction are detailed in Sect.3. In Sect.4, we present our results using both empirical and synthetic model approaches.Our ForMoSA forward modeling tool is used to explore the different atmospheric models and determine the most probable physical properties of these three companions.In Sect.5, we focus our study on the emission-line properties observed for the three companions.In Sect.6, we finally summarize and discuss our results in the context of previous work, and present perspective of future studies.References.2013); (d)  Lachapelle et al. (2015).

Target description
The three companions have close physical properties and spectral types (M8-M9), yet they were selected because they come with different configurations (mass ratio with the host star, projected separation; see Table 1): HIP 78530 B. Was identified by Kouwenhoven et al. (2005) at a separation of 4.536 ± 0.006 from the ∼2.5 M B9V star (Houk & Smith-Moore 1988) HIP 78530 A. Lafrenière et al. (2011) confirmed the companion is co-moving with the primary star.The projected separation between the two objects now corresponds to 623 ± 8 au using the Gaia-DR2 distance (137.2 ± 1.5 pc; Gaia Collaboration 2018) The Banyan Σ tool (Gagné et al. 2018) and the DR2 astrometry confirms the star is a high-probability member (99.9%) of the Upper-Scorpius association.Lafrenière et al. (2011) has provided a medium-resolution (R λ ∼ 5300 to 6000) spectrum of the companion covering the 1.15−2.40µm wavelength range.Lachapelle et al. (2015) also presented a lower-resolution spectrum (R λ ∼ 1350) but extending down to 1 µm.The spectra confirm the companion is a young M 7 ± 0.5 dwarf.Lachapelle et al. (2015) estimate a T eff of 2700 ± 100 K and a luminosity of log (L/L ) = −2.53± 0.09 relying on the HIPPARCOS distance available at that time.These latter authors find a mass of 23 ± 2 M Jup assuming an age of 10 Myr for Upper-Scorpius.The primary star has no noticeable excess emission (Carpenter et al. 2009;Luhman & Mamajek 2012;McDonald et al. 2017) and low extinction (A v = 0.48 Neckel & Klare 1980).
HIP 77900 B. Was identified in UKIDSS and Pan-STARRS 1 images by Aller et al. (2013) from its red colors.It is located at a projected separation of 21.8 from the B6 star (Garrison 1967) HIP 77900, a high probability member of Upper-Sco (97.6% membership probability according to Banyan Σ).Unlike for the case of HIP 78530 B, the authors did not re-observe the target to check whether it shares the proper motion of the primary star.However, they obtained a low-resolution (R λ ∼ 100) spectrum of the source covering the 0.8−2.5 µm range whose features are indicative of a young M 9 (±0.5)object from Upper-Sco, and therefore argued that HIP 77900 B is bound to the star.We re-discuss the physical association of the two objects in Sect.6.1.
USco 1610-1913 B. Was identified by Kraus & Hillenbrand (2008) and confirmed to be co-moving with the K7 star GSC 06209-00691 (Kraus & Hillenbrand 2008, 2009;Kraus et al. 2014).The companion was last found at 5.837 ± 0.006 from the star, now corresponding to a separation of 779 ± 9 au at the Gaia-DR2 distance of the primary.Aller et al. (2013) presented a 0.8−2.5 µm low-resolution (R λ ∼ 100) spectrum of the object confirming its youth and substellar nature.These latter authors estimated a M9 ± 0.5 spectral type and found the companion to be four times overluminous with respect to HIP 77900 B which shared the same spectral properties at this spectral resolution.The spectral type was confirmed by Lachapelle et al. (2015) from a medium-resolution (R λ ∼ 5300 to 6000) NIR (1.15−2.40µm) spectrum of the companion.These latter author find T eff in the range 2300−2700 K for the object using DRIFT-PHOENIX and BT-SETTL atmospheric models with a large spread in T eff values depending on the wavelength interval considered for the fit.Their analysis also confirms the over-luminosity of the object for the inferred temperature.However, the star had no measured parallax at that time, and the average distance of Upper-Sco from de Zeeuw et al. (1999) was assumed.A second companion (hereafter USco 1610-1913 Ab) was discovered at a projected separation of 19.4 ± 0.3 au (Kraus et al. 2008) from USco 1610-1913A. USco 1610-1913 Ab has a mass of 103 ± 24 M Jup using the contrast reported in Kraus et al. (2008), the new Gaia-DR2 distance of USco 1610-1913 A (assuming it is unbiased; see Sect.6.1), and the Baraffe et al. (2015) tracks at an age of 5-11 Myr.
The six young and isolated free-floating objects (see Table 2) were selected from the sample of Lodieu et al. (2008) late-M/early-L Upper-Sco brown dwarfs.These latter authors also report low-resolution (R λ ∼ 1700) spectra over 1.15−2.50µm of the sources.Low-resolution (R λ ∼ 1350) multi-epoch redoptical (0.57−0.88 µm) spectra of these objects have also been collected by Lodieu et al. (2011).The spectra of USco J1610-2239 and USco J1607-2211 exhibit a H α line.No emission line was detected at that time in the spectra of the remaining objects.Dawson et al. (2013) and Luhman & Mamajek (2012) do not find noticeable infrared excess (up to 22 µm) for USco J1610-2239, USco J1607-2211 or USco J1607-2242.This is not the case for USco J1606-2335 and USco 1608-2315.The former is found to have Spitzer [4.5] and WISE W2 photometry indicative of a disk excess (Luhman & Mamajek 2012), and the latter is found to have a W1-W2 color suggestive of an excess (Dawson et al. 2013) but that excess is not confirmed by Luhman & Mamajek (2012) using similar data.Gaia-DR2 distances of 128.5 +15.4  −12.4 pc and 119.3 +20.8  −15.4 pc are available for USco J1610-2239 and USco J1607-2211.Using Banyan Σ (Gagné et al. 2018), we find that the two objects have a 99.6 and 94.9% chance, respectively, of belonging to Upper-Scorpius based on kinematics.USco J1608-2315 also has a Gaia-DR2 parallax value (π = 4.5073 ± 1.1874 mas).The object is found at a greater distance than the typical Upper-Sco members and therefore Banyan Σ estimates a lower probability of membership to the association (59.7%).The large error of DR2 for this object leads us to doubt this probability, noting that the target has spectral features clearly indicative of a membership to Upper-Sco (e.g., Lodieu et al. 2008;Bonnefoy et al. 2014).To conclude, adaptive-optics imaging and sparse aperture masking observations (Biller et al. 2011;Kraus & Hillenbrand 2012) of the sources exclude companions with similar masses down to ∼10 au.The optical and NIR spectra of the free-floating sources presented here have already been briefly introduced in Lodieu et al. (2018).These sources have increased spectral resolution or are extended at shorter wavelengths compared to previous spectroscopic data, and thus allow us to study the long-term variability of the H α line.We therefore re-analyze them in order to use them as empirical templates for the characterization of the three companions HIP 78530 B, HIP 77900 B, and USco 1610-1913 B.

Observations and data reduction
We used the X-shooter seeing-limited medium-resolution spectrograph mounted at UT2 Cassegrain focus (Vernet et al. 2011).The wide wavelength coverage of the instrument (300-2480 nm) is ideally suited for the characterization of accreting BDs with emission line series.We chose the 1.6 , 1.5 , and 1.2 -wide slits for the UVB, VIS, and NIR arms corresponding to resolving powers R λ = λ/∆λ = 3300, 5400, and 4300, respectively.This setup was adopted for all our targets.The observing log is reported in Table 2.The slits were oriented perpendicular to the position angles of the companions in order to mitigate the flux contamination of the host stars.Each target was observed following an ABBA strategy to evaluate and remove the sky emission at the data processing step.Spectrophotometric standard stars were observed as part of the observatory calibration plan; they are not reported in the log as they were not used for our reduction or analysis.
We used the ESO reflex data-reduction environment (Freudling et al. 2013) to run the X-shooter pipeline version 2.9.3 on the raw data (Modigliani et al. 2010).The pipeline produces two-dimensional, curvature-corrected, and flux-calibrated spectra for each target and epoch of observation (trace).The spectra were extracted from the traces using a custom IDL script.The flux in each wavelength channel at the position of the source was averaged within 720 mas aperture in the UVB and VIS arms, and a 1120 mas aperture in the NIR arm.The script computed the noise at the position of the source into each spectral channel following the procedure described in Delorme et al. (2017a).The residual nonlinear pixels in the spectra were removed using the kappa-sigma clipping method.The telluric corrections were evaluated and removed using the molecfit package (Smette et al. 2015;Kausch et al. 2015).The spectra at each epoch were corrected from the barycentric velocity and renormalized using the epochs when the sky transmission was photometric as an anchor point.Our flux calibration was checked by computing the 2MASS or MKO synthetic photometry from the spectra and comparing the values to published ones (Lodieu et al. 2008).

Physical properties
4.1.Empirical analysis Lachapelle et al. (2015) derived the spectral types of USco 1610-1913 B and HIP 78530 B using NIR spectra of these companions.We noticed significant differences between them and the X-shooter spectrum of HIP 78530 B presented here (see Fig. 1).We noticed the same difference for the same target in comparing with the spectrum from Lafrenière et al. (2011).The X-shooter spectrum of USco 1610-1913 B is consistent with the lowresolution spectrum of Aller et al. (2013).We therefore conclude that the differences in the J-band spectrum of HIP 78530 B may stem from flux losses affecting long-slit observations with AO-fed spectrographs.In J-band, we identified several strong absorption lines, like the neutral sodium (Na I) doublet (1.138 and 1.141 µm), and the neutral potassium (K I) doublets (1.168, 1.177 µm and 1.243, 1.254 µm).Neutral iron (Fe I) lines (1.189 and 1.197 µm) are also present as well as weaker Na I (1.268 µm), magnesium (Mn I at 1.290 µm), and titanium (Ti I at 1.283 µm) lines.We notice the typical triangular shape of the H-band which is produced by wide H 2 O absorption bands and testifies to the young ages of our objects.In K-band we found the calcium (Ca I) triplet (at 1.98 µm) and a weak Na I doublet (2.206, 2.209 µm).We also detect CO bands (2.295, 2.322, 2.352 µm).

Spectral type and surface-gravity determination
We re-investigated the optical classification of our targets using a standard χ 2 comparison of our spectra to empirical templates from the Ultracool RIZzo Spectral Library1 .The RIZzo library is made of spectra of 265 M2-L5 brown dwarfs from Cruz et al. (2003), Cruz et al. (2007), and Reid et al. (2008).We restrained the fit to the 0.75-0.86µm range.The results of our spectral fitting are shown in Fig. 4 and are reported in Table 3.The spectral type errors are due to the subgroup increment defined in the library.For USco 1610-1913 B, HIP 77900 B, and HIP 78530 B, we find spectral types of M 9 ± 0.5, M 9 ± 0.5, and M 8 ± 0.5, respectively, with best fits 2MASS J11582484+1354456 (M 9) and 2MASS J07140394+3702459 (M 8) which are both freefloating objects (determined from BANYAN Σ Tool).We also classified our targets in the NIR using absorption lines respecting the Allers & Liu (2013) scheme.In NIR, we find later spectral types for USco 1610-1913 B and HIP 78530 B than the ones derived by Lachapelle et al. (2015), consistent with the redder slope of the X-shooter spectrum of HIP 78530 B, and because of the revised extinction A v values (see also Sect.4.3) considered for the two systems.The optical spectral types derived for the free-floating BDs are in agreement with those of Lodieu et al. (2018) within error bars.
Appendix B shows the systematic differences between each method.The spectral type derived from the H 2 O index seems to be one to two subtypes over the one from the visual method.Allers & Liu (2013) explain that the H 2 O index could be sensitive to gravity and so can be biased.We choose the spectral type from the visual comparison to avoid this bias.
In addition to the spectral-type determination, we also applied the surface-gravity classification formalized by Allers & Liu (2013).The results are reported in Table 3 and show that HIP 77900 B and HIP 78530 B are identified as young, intermediate-surface-gravity BDs.USco 1610-1913 B is confirmed as a very-low-surface-gravity BD.The gravity class of HIP 78530 B is consistent with that derived by Lachapelle et al. (2015).Our results are consistent with those of Lodieu et al. (2018) for the classification of the young free-floating objects.

Over-luminosity of USco 1610-1913 B
Figure 5 shows the comparison between the visible and NIR calibrated flux of USco 1610-1913 B with the fluxes of young reference BDs of similar spectral types from our sample with Gaia parallaxes.We also include a comparison to the X-shooter spectrum of TWA 26 (Manara et al. 2013).All BDs have been scaled to the distance of USco 1610-1913A.All these objects reproduce the detailed absorptions and pseudo-continuum shape of USco 1610-1913 B spectrum provided that an extra scaling factor of two (TWA 26) and four (USco 1610-2239, HIP 77900 B) is considered.This over-luminosity of USco 1610-1913 B has already been noted by Aller et al. (2013) and Lachapelle et al. (2015).We confirm it over our extended wavelength range and resolution relying on an extended set of comparison objects from the association with now-published parallaxes.We discuss the possible origins of the over-luminosity in Sect.6.Table 3. Results of the empirical analysis.

Description of the atmospheric models
We used a forward-modeling approach to determine the atmospheric parameters of the young BDs observed in this study.Forward modeling codes enable comparison of the object spectrum to pre-computed grids of models which include our best knowledge of atmospheric physics.We used grids of synthetic spectra produced by the BT-SETTL15 model (Allard et al. 2012).This model handles the radiative transfer using the PHOENIX code (Hauschildt et al. 1997;Allard et al. 2001).It accounts for convection using the mixing-length theory, and works at hydrostatic and chemical equilibrium.The opacities are treated line by line (details on each element are given in Rajpurohit et al. 2018).The code models the condensation, coalescence, and mixing of 55 types of grains.The abundances of solids are determined comparing the timescales of these different processes at each layer.In this study, we considered the predictions of the BT-SETTL15 model with T eff ranging from 2100 to 3000 K (in steps of 100 K), and a range of log (g) from 3.5 to 5.5 dex (in steps of 0.5 dex).We assumed a solar metallicity M/H = 0.0, in broad agreement with the values reported in Sco-Cen (Bubar et al. 2011).These grids have been extensively used in previous studies of young BDs (Manjavacas et al. 2014;Bonnefoy et al. 2014;Delorme et al. 2017b;Bayo et al. 2017), but never as part of a Bayesian methodology, as developed here with the ForMoSA code.

The ForMoSA code
For this work, we chose to develop our own forward modeling code relying on the Nested Sampling procedure (Skilling 2006).1610-1913 B (black) to those of reference objects scaled to the distance of the companion (green).The companion is 2.0-3.9 times more luminous than objects with the same spectral type.
The method explores and recursively isolates different patches of likelihood levels in the parameter space.It offers several advantages over classical Markov-chain-Monte-Carlo algorithms also used for the forward modeling (e.g., Bonnefoy et al. 2018;Samland et al. 2017).The approach avoids missing local minimum within vast and degenerate parameter spaces while ensuring the convergence of the exploration.
The code ForMoSA (for FORward MOdeling for Spectral Analysis) takes as input an observed spectrum with associated error bars and any grid of synthetic spectra (here the BT-SETTL15 grid).The Nested Sampling is handled by the nestle Python module2 .ForMoSA re-samples the data and the models in order to make them comparable.The grid of model spectra are first interpolated onto the wavelength grid of the observation and degraded to the spectral resolution of the instrument that acquired the data.The optical spectra of our targets and the models were degraded to R λ = 3300, 5400 and 4300, similar to the UVB, VIS, and NIR parts of the spectrum, respectively.
In our case, we hypothesize that data points are independent in the observed spectra, and therefore the likelihood is derived from the χ 2 value.To compute χ 2 at each step, we need to generate a model spectrum for a set of free-parameters that does not necessarily exist in the original grid of spectra.Therefore, we generated a model spectrum on demand following a two-step process: A preliminary phase consists in reducing the grid meshes.To do so, ForMoSA interpolates and reduces each T eff increment to 10 K and log (g) increments to 0.01 dex.We considered the linear and bicubic spline interpolation approaches, and finally selected the bicubic spline interpolation which better accounts for the flux variation through the grid.This step needs to be taken only once and ensures a regular grid.In doing this, we increase the accuracy of the second interpolation phase.The second phase arises in the course of the nested sampling process when a new point in the parameter space is defined.The closest neighbors approach was found to provide the best tradeoff between the reliability of the interpolation process and the computation time needed to run the interpolation.
Each synthetic spectrum gives the flux at the top of the atmosphere.The comparison with the observed spectrum also requires multiplication of the model by a dilution factor , with R the object radius and d the distance.We adopt distances for the host stars from the Gaia DR2.We considered flat priors on T eff (2100-2900 K) and log (g) (3.5-5.5 dex).Here, T eff and radius are linked together by the luminosity (Boltzmann law) so we have chosen a flat prior on R (0.5-30.0 R Jup ) to be conservative and ensure we do not limit the exploration of T eff .We also provide the luminosity from the posterior distributions of R and T eff ForMoSA can also consider interstellar extinction (A v ) as a free parameter using the extinction law from Fitzpatrick & Massa (2007), the radial velocity (RV) by a Doppler shifting law 3 , and the projected rotational velocity v sin (i) according to the rotational broadening law from "The Observation and Analysis of Stellar Photospheres" by Gray 4 .Figure 6 illustrates a typical output of ForMoSA on HIP 78530 B. The method allows us to identify correlations between parameters such as T eff and R.

Results
As a safety check, the ForMoSA code was applied to our three companions HIP 78530 B, HIP 77900 B, and USco 1610-1913 B and two well-known young BDs, TWA 26 and TWA 29, with similar spectral types and ages (∼8 Myr) and observed by Manara et al. (2013).We considered wavelengths from 1.0 to 2.5 µm in all our fits to avoid biases related to the residual contamination of HIP 78530 B and ran a homogeneous analysis on all objects.The results are shown in Figs.From these maps, we calculate the A 0 extinction at 550.0 nm and assume that the difference with the A v extinction is negligible (A 0 A v ).The fitting spectral range is indicated by the yellow background.
Panel e shows the best fit when using the complete J + H+ K-band spectral range, and adding the extinction value A v this time as a free parameter.We reddened the models with an extinction function in the code (Fitzpatrick & Massa 2007).
Panel f shows the best fit focused on the optical part (0.56-1.00 µm for all targets but 0.80-1.00µm for HIP 78530 B).We estimate the pseudo-continuum from the original X-shooter spectrum by degrading it to a very low spectral resolution (R λ = 100).We then subtract this pseudo-continuum from the original section to highlight line features.We add the v sin (i) as a free parameter with this fit.We determine the radius analytically with the relation from Cushing et al. (2008).
Panel g shows of zoomed view on a the gravity-and metallicity-sensitive lines of K I and Na I with the best-fitting solution.We subtracted the pseudo-continuum using the same 0.3 0.6 0.9 Av = 0.0 method as in panel f.The determination of the radius is also analytic.
Generally, the results show that the BT-SETTL15 models fail to consistently reproduce the pseudo-continuum of the X-shooter spectra of all objects at all wavelengths (visible and NIR) with the same range of physical parameters.The same problem arises when considering a χ 2 test on the original grid of BT-SETTL15 spectra.This is highlighted in panels a-d of all figures, where the best-fit solution varies considerably depending on which spectral range is fit.Similar discrepancies were evidenced by Manjavacas et al. (2014) using older releases of the models (Allard et al. 2011(Allard et al. , 2012(Allard et al. , 2013;;Lachapelle et al. 2015;Bayo et al. 2017).Therefore, the problem remains in the 2015 release of the models.
The posteriors of T eff are strongly tied to the pseudocontinuum shape.Consequently, they are affected by the choice of the spectral range used to estimate the best fit.For instance, the T eff is generally higher from a fit using the J-band compared to a fit using H-or K-band.The surface-gravity is also affected as it remains sensitive to the shape of the VO and H 2 O absorption bands.Notably, the radius, which is linked to the flux dilution factor, seems to be mostly consistent when using three different fitting bands.In addition to the best-fitting solution, the ForMoSA code provides the errors on the posterior solutions.The errors on the posterior solutions are tightly connected to the errors bars on the spectra themselves by the likelihood function.Given the high signal-to-noise ratio (S/N) of the X-shooter spectra for HIP 78530 B, HIP 77900 B, USco 1610-1913 B, TWA 26, and TWA 29, the resulting fitting error bars are about one to two orders of magnitude smaller than the systematic errors associated with the choice of the spectral fitting window and the non-perfect match of the BT-SETTL15 models with our spectra at all observed wavelengths (see Tables C.1-C.2).Our adopted solutions in T eff , log (g), and R are therefore derived from the average and dispersion of the solutions from the different spectral fitting windows with the extinction value given by Lallement et al. (2019).The bolometric luminosity is calculated with the Stefan-Boltzmann law.
In addition to the different spectral fitting window, the ForMoSA code was also applied with the extinction A v as a free parameter.The resulting posteriors can be directly compared with the ones obtained with the Gaia DR2 values at the location of our targets (Lallement et al. 2019).Surprisingly, with the simple use of an interstellar extinction law correction, we considerably improve the goodness of the fits of all targets (see Figs. objects are affected, it is very likely that the origin is not circumstellar, but linked to the current state-of-the-art of BT-SETTL15 models, as is for example the complexity of the dust formation process or the dust opacity (Pavlenko et al. 2007).Systematically testing how this deficiency evolves with lower effective temperature, age (surface-gravity), and gravity would be very interesting.
As an additional step to exploit the medium resolution of our X-shooter spectra, the radial and rotational velocity information (line profile and offset) has been incorporated into ForMoSA as additional fitting parameters.The RV can be adjusted on the full spectral window.For the adjustment of v sin (i), we restricted the fitting window to the optical part (0.56-1.00 µm) with the continuum subtracted to minimize the computation time.This choice is driven by the higher spectral resolution in optical (R = 5400) and the presence of various absorption lines (NaI, FeH or VO).In Tables C.1 and C.2, one can see that the RV results strongly depend on the choice of the fitting spectral window.We adopted the RV value from the optical fit for the same reasons as for the v sin (i).The two identified causes of this discrepancy are the slope issues in models and a systematic error in the wavelength calibration of the spectra.We therefore decided to independently compute the RVs using a cross-correlation approach between our observed spectra and a spectral template generated by the BT-SETTL15 code at T eff = 2400 K and a log (g) = 4.0 dex (see Fig. 12).For most cases, there is a good agreement between the cross-correlation results and the ForMoSA ones in the optical that we therefore adopt as final values.For the rotational velocity, the limited X-shooter spectral resolution does not allow the measurement of the low v sin (i) of our targets, but rather places an upper limit of roughly ≤50 km s −1 .This limit is consistent with the 55 km s −1 limit expected for a resolution of about R = 5400 in the visible arm.Considering the adopted values of T eff , log (g), R, RV, v sin (i) and L using the fixed extinction values of Lallement et al. (2019), we now detail the results of ForMoSA target by target.
TWA 26 and TWA 29.Both targets are well known and their X-shooter spectra have been analyzed as a sanity check of the ForMoSA code.They have been classified by Manara et al. (2015) as late-M dwarfs with spectral types M9γ and M9.5, respectively.Both are members of the TW Hydrae association, with an age estimate of 8 Myr, and are located at a distance of 49 pc and 83 pc, respectively.They have no extinction observed in their line of sight.Both were characterized by Filippazzo et al. (2015) using a complimentary approach to our work.Exploiting the spectral energy distribution from the visible to the NIR and the mid-infrared (MIR), Filippazzo et al. (2015) derived the bolometric luminosities of these two young brown dwarfs with known distances.Using evolutionary model   Filippazzo et al. (2015) derived the predicted radii and masses, and therefore the object effective temperatures using the Stefan-Boltzmann law as well as the surface gravities.The results are reported in Table 4, and roughly agree within the error bars with those obtained by the ForMoSA forward modeling of the X-shooter spectra of TWA 26 and TWA 29.As the three Upper-Sco companions HIP 78530 B, HIP 77900 B, and USco 1610-1913 B scan a similar range of temperatures and surface gravities as TWA 26 and TWA 29, this test appears to support the use of ForMoSA to explore their physical properties.These three companions have no known MIR photometry or spectroscopy, and therefore a direct determination of the bolometric luminosity as done by Filippazzo et al. (2015) is not possible here.Based on a fit using the optic wavelength range, we derive a v sin (i) = 66 ± 2 km s −1 for TWA 29.This is the only object for which we are not limited by the resolving power of the instrument.
USco 1610-1913 B .For their study of this companion, Lachapelle et al. (2015) used an incorrect value of A v = 1.1 given by Carpenter et al. (2009).Using the method described in Carpenter et al. (2008Carpenter et al. ( , 2009) ) derived this relatively high A v using the optical (Johnson B and V, Tycho B T and V T ) and NIR (2MASS J, H, and K s ) photometry and colors, but an incorrect spectral type (K7) for USco 1610-1913.As described in Appendix A, an updated analysis of the full spectral energy distribution (SED) indicates a later spectral type M0, leading to an extinction value of A v = 0.09 ± 0.01 mag, which is more consistent with the value from Lallement et al. (2019) : A v = 0.13.Using A v = 0.13 with ForMoSA, we derive an effective temperature of T eff = 2542 +68 −104 K and a surface-gravity of log (g) ≤ 4.17 dex, both compatible with a young M9 brown dwarf of very low-gravity for USco 1610-1913 B. However, we find clear nonphysical solutions for the radius and the luminosity (for a young late-M dwarf), which confirms the strong over-luminosity of this source, as already pointed out in Sect.4.1 and discussed further in Sect.6.We notice a difference between the RV found from ForMoSA and the one found from the cross-correlation algorithm.
HIP 78530 B. The extinction maps of Lallement et al. (2019) give a value of A v = 0.075 along the line of sight of HIP 78530, somewhat different from the value of A v = 0.5 used by Lachapelle et al. (2015) and from Carpenter et al. (2009).We adopted the former value in our analysis.ForMoSA finds T eff = 2679 +118 −119 K, in agreement with the M8 spectral type of the source (Luhman et al. 2003;Filippazzo et al. 2015, this work).The analysis of the continuum-subtracted spectrum yields a surface-gravity estimate compatible with the intermediategravity class of the object (Sect.4.1).The luminosity value of log (L/L ) = −2.87+0.15  −0.15 is also compatible with expected log (L/L ) = −2.5 to −3.0 found by Filippazzo et al. (2015) for young M8β and M8γ BDs.Finally, the effective temperature T eff and surface-gravity are consistent with the values from Lafrenière et al. (2011).
HIP 77900 B. With an extinction value of A v = 0.07 (Lallement et al. 2019), the resulting effective temperature of T eff = 2602 +117 −97 K given by ForMoSA is consistent with the later M9 ± 0.5 spectral type empirically derived for this companion.The estimated surface-gravity of log (g) ≤ 4.36 dex is consistent with the intermediate-gravity classification of Allers & Liu (2013), and the luminosity value of log (L/L ) = −2.89+0.15  −0.13 is compatible with the values found for young M9βγ brown dwarfs (Filippazzo et al. 2015).

Emission line properties
Our new high-S/N medium-resolution optical spectra allow us to identify the faint H α (656.28 nm) emission lines of the three companions and of four out of the six free-floating objects of our sample.Figure 14 shows the line profiles at each epoch.We checked for the emission line detection in the traces to ensure it was not produced by an uncorrected bad pixel (see Sect. 3).H α lines are commonly attributed to magnetospheric accretion (Manara et al. 2013(Manara et al. , 2015;;Natta et al. 2004) or chromospheric activity (e.g., White & Basri 2003).The three objects (USco 1610-1913B, USco 1606-2335, and USco 1608-2315) showing the strongest H α lines also display other Balmer lines, H β (486.1 nm) and H γ (434.0 nm), in their spectra.We tentatively identify Ca II-H (393.4 nm) and Ca II-K (396.9 nm) emission lines in the spectrum of USco 1610-1913 B as well.The Ca II lines can also be related to accretion or activity (Herczeg & Hillenbrand 2008).We report the apparent fluxes of these emission lines and 10% width of the H α line in Table C.3.The noisy spectral continua of our objects prevented us from obtaining a robust determination of the continuum and the equivalent widths.Both origins (accretion and activity) for our targets are investigated below.

Accretion rate determination
The accretion rates of each object are reported in Table C.5.They were derived following the relation of Alcalá et al. (2017): where L acc is the accretion luminosity, R in is the inner disk radius, R obj and M obj are the radius and mass of the objects, respectively.This is assuming R obj /R in = 0.2 (Alcalá et al. 2017) which holds in our mass range, irrespective of the object configuration (companion or isolated).The accretion luminosities were estimated beforehand from the line fluxes of the Ca II lines and the Balmer lines using the relationships given in Rigliaco et al. (2012) and Alcalá et al. (2017), respectively.Radii and masses were in addition inferred from the evolutionary model predictions as indicated in Table 5.Alternatively, we derived accretion rates based on 10% H α width following the relations of Natta et al. (2004): log ( Ṁacc ) = −12.89(± 0.3) + 9.7(± 0.7) × 10 −3 H α 10%. (2) The 10% H α accretion rates are also reported in Table C.5.They tend to be higher than the ones inferred from the Balmer and Ca II line fluxes.A similar divergence between the different accretion proxies is found for the young companion SR 12C (Santamaría-Miranda et al. 2018).The accretion rates of USco 1610-1913 B are compatible with those of the 1-5 Myr late-M sub-stellar companions FW Tau b (Bowler et al. 2014), DH Tau b (Zhou et al. 2014), and GQ Lup b (Wu et al. 2017b), which had the lowest recorded accretion rates among known accreting companions.The accretion rates of HIP 78530 B and HIP 77900 B are found to be an order of magnitude lower, indicating that these objects are nonaccretors.

Chromospheric activity
The weakness of the lines suggests the emission arises mostly from chromospheric activity.Manara et al. (2013) provides the L acc,noise /L bol values expected from chromospheric activity for young stars and BDs down to the M9 spectral type: log (L acc,noise /L bol ) = (6.17± 0.53) × log (T eff ) − (24.54 ± 1.88).
(3) All of our objects apart from USco 1608-2315 have log (L acc /L bol ) values below or comparable to the activity threshold log(L acc,noise /L bol ) corresponding to their T eff (estimated from ForMoSA), indicating that the lines are indeed likely related to chromospheric activity.Herczeg et al. (2009) also identify a H α line on USco 1607-2211 of similar intensity and came to the same conclusion.We note that this criterion applied to the H α line of GQ Lup b as observed in April 2016 (Wu et al. 2017b) indicates that the line mostly arises from chromospheric activity in spite of accretion at this epoch5 .The 10% width of the H α line may also discriminate between the accretion (large values) and chromospheric activity (small values).A threshold of 200 km s −1 is generally adopted for BDs (Jayawardhana et al. 2003).All of our objects but USco 1608-2315 have a 10% H α below the 200 km s −1 limit, supporting line production by chromospheric activity.Therefore, USco 1608-2315 is likely the only true accreting object in our sample.The source also has a tentative disk excess (Sect.2).Thus, we add USco 1608-2315 to the short list of free-floating accretors with masses below 30 M Jup (Herczeg et al. 2009;Bayo et al. 2012;Joergens et al. 2013;Alcalá et al. 2014;Boucher et al. 2016;Lodieu et al. 2018).
To further support these conclusions for a few specific objects, Fig. 13 (Manara et al. 2013).The H γ and H δ lines have similar strengths but the H β line appears more luminous in the companion spectrum.
Figure 14 shows that the two free-floating objects with the strongest lines (USco 1608-2315 and USco 1606-2335) display significant short-term variability.The lines are not detected in the optical spectra of these objects obtained in May 2007 (Lodieu et al. 2011).The H α lines of USco 1607-2211 and USco 1610-2239 do not show obvious variability in our data, but are not detected in May 2007 (Lodieu et al. 2011) either.It is possible that the nondetections from 2007 arise from a degraded sensitivity of the low-resolution data from Lodieu et al. (2011).Active M-dwarfs are known to display such variability from chromospheric activity (e.g., Gizis et al. 2002), which strengthens our conclusions on the line origins.The variability of USco 1608-2315 could also stem from variation of the accretion rate (Aguayo et al., in prep.; Bonnefoy et al., in prep.).

Revisiting the physical association of HIP 77900 B and
USco 1610-1913 B with Gaia HIP 77900 B and USco 1610-1913 B are sufficiently bright and distant from their primary for them to have reported Gaia parallaxes.This is to our knowledge the first case for which the physical association of young imaged BD companions with their primary stars can be investigated based on the individual five-parameters astrometric solutions of the system components.

The case of HIP 77900 B
The Gaia-DR2 parallax of HIP 77900 A confirms that this system would be extreme if bound (22.3 = 3375 au projected separation).Very wide systems such as HIP 77900 are expected to be rare (e.g., Baron et al. 2018).The companion has not been observed at multiple epochs so it could in principle be a background star.However, the close resemblance of the spectrum of the companion to those of free-floating analogs from our library of Upper-Sco objects identifies it as a likely member of this association.Therefore, HIP 77900 A and B are either coeval or aligned by chance within the association.Their respective distances (191 −30 +43 pc for A, and 151.4 +2.8 −2.7 pc for B) reveal a 1.4 σ difference.There is in addition a 1.8 σ difference between the proper motion in declination of the two objects.The p-value on the difference on the three astrometric parameters of A and B (0.06) does not allow us to firmly conclude whether or not the system is coeval.The Re-normalized Unit Weight Error (RUWE) values of A and B (0.97 and 1.21) indicate robust solutions.The galactic Cartesian coordinates of HIP 77900 A (X A = 137.8± 2.6 pc, Y A = −35.5 ± 0.7 pc, Z A = 51.8 ± 1.0 pc) computed following Gagné et al. (2014) are well within the range of expected values for Upper-Sco (Galli et al. 2018), contrary to those of B (X B = 173.8+39.1 −27.3 pc, Y B = −44.7 +7.1 −10.1 pc, Z B = 65.4 +14.8 −10.3 pc) which is at odds with the youth of the object.However, we find photometric distances of 151 +35 −24 pc and 150 −17 +22 pc for HIP 77900 B using the fluxcalibrated spectra of the spectrally analogous M8.5 dwarfs from our sample (USco 1607-2211 and USco 1610-2239, respectively) with Gaia parallaxes.These photometric distances are more consistent with the Gaia distance of A and would strengthen the case of a coeval system.However, we found a RV of 19.3 ± 1.2 km s −1 for the B component which is not consistent with the 1.3 ± 2.6 km s −1 found by Gontcharov (2006) for HIP 77900 A. This result raises severe doubts about the gravitational link between the two objects, and the next Gaia data releases of the relative astrometry of A and B will be needed to firmly conclude on their physical association.

The case of USco 1610-1913 B
The spectral properties and measured distances of USco 1610-1913 B clearly confirm its membership to Upper-Sco.USco 1610-1913 B and A have been shown to have a common proper motion from 2007 to 2012 and have been proposed to be bound (Kraus et al. 2014, and ref therein).However, the recent DR2 solutions of these two objects diverge in proper motion (2.0 and 2.4σ significance in µ α and µ δ .cos(δ),respectively) and in distance (1.4σ; B is 10.3 pc closer than A).A χ 2 test on the difference on the three astrometric parameters of each system component (parallaxes and proper motions) taking into account the correlations gives a p-value of 0.01 which is in favor of a physical bond between the two objects.The p-value accounts for the slight underestimation of the errors (Lindegren et al. 2018).The reliability of the Gaia astrometry may however be questioned for this particular system.Indeed, the ∼0.145M-dwarf companion to USco 1610-1913 A (Kraus et al. 2008) is unresolved in the DR2.The RUWE index is proposed as a more reliable and informative goodness-of-fit statistic than for instance the astrometric excess noise (Lindegren et al. 2018, Gaia technical note Gaia-C3-TN-LU-LL-124-01).The RUWE of USco 1610-1913 A (1.63) is above the threshold of 1.4 and indicates that the observations are inconsistent with a simple five-parameter astrometric model while B (RUWE = 1.19) shows a more reliable solution.The next releases of Gaia will here again solve this ambiguity.
If we keep the hypothesis that the system is bound (or even coeval) and located at the distance of the B component, USco 1610-1913 B remains 4.5 times more luminous than free-floating analogs from the Upper-Sco association.This is illustrated in Fig. 5 using USco 1607-2211 and USco 1610-2239 for comparison, which both have high-quality astrometric solutions.In addition, the probability of a chance alignment within the association is small.Aller et al. (2013)  Table 5.Comparison between parameters from ForMoSA and parameters in using evolutionary models (Chabrier et al. 2000) (BT-SETTL15/ForMoSA independent).

ForMoSA
Bolometric correction + Evolutionary models Object log (L/L ) R log (g) USco 1610-1913 B (1) −2.25 Notes.The luminosity has been determined using the bollometric corrections at J and K bands from Filippazzo et al. (2015).For the luminosity of our reference objects we directly use the Filippazzo et al. (2015) values.We also provide the mass values. (1)With the hypothesis that companions are at the same distance as the primary star.
and the cross-correlation approach (11.5 ± 0.4 and 0.9 km s −1 , respectively).Both values are in addition inconsistent with the RV derived by Dahm et al. (2012) (2013) proposed that differences in the accretion history of USco 1610-1913 B could play a role in such a discrepancy as proposed by Baraffe et al. (2012).For late K-and M-dwarfs, chromospheric activity as found in our spectra decreases the objects T eff and increases their radii (López-Morales 2007;Morales et al. 2008).Based on a sample of 669 M < 1 M nonaccreting stars from the Palomar/Michigan State University catalog (PMSU; Reid et al. 1995;Hawley et al. 1996), Stassun et al. (2012) proposed empirical relations to determine the bias induced by chromospheric activity: (5) Despite these variations of the radius and T eff of the atmosphere, the bolometric luminosity remains constant.This effect cannot explain the over-luminosity of USco 1610-1913 B. An alternative explanation would be that the companion is an unresolved quadruple system presumably with all components having similar spectral types in this case.Such systems are expected to be rare among A-and FG-type stars (Raghavan et al. 2010;De Rosa et al. 2014;Tokovinin 2014).The system would be even more exotic and rare considering that USco 1610-1913 A is a binary and that USco 1610-1913 B has to be composed of four identical objects to account for its overluminosity.Quadruples tend to be found as 2+2 tight binary systems.The discovery of a quadruple system made of two pairs of M5 eclipsing binaries in Upper-Sco (Wang et al. 2018) shows that such objects exist in isolation.To our knowledge, adaptive-optics images of the system have not displayed PSF-elongation or resolved the companion as a higher-order object but new NIR high-resolution (R 80 000 to 100 000) spectrographs (ESO/NIRPS, CRIRES+) could investigate the multiplicity of USco 1610-1913B in the near future.

Conclusion
Based on medium-resolution spectra obtained with the X-shooter spectrograph at the VLT, we carried out an indepth characterization of three low-mass BD companions on wide-orbits, namely USco 161031.9-16191305B, HIP 77900 B, and HIP 78530 B of the Upper-Scorpius association, together with six young isolated brown dwarfs of similar spectral types and ages.The X-shooter data yield the first medium-resolution optical spectra of the companion objects.We can summarize the main results as follows: 1.The re-investigation of the spectral classification of the three companions USco 1610-1913 B, HIP 77900 B and HIP 78530 B confirms that they have spectral types M9 ± 0.5, M9 ± 0.5, and M8 ± 0.5, respectively.HIP 77900 B and HIP 78530 B are identified as young, intermediatesurface-gravity brown dwarfs, whereas USco 1610-1913 B is confirmed as a very-low-surface-gravity brown dwarf.
2. The development and use of the ForMoSA forward modeling code relying on the Nested Sampling procedure enables us to infer posterior probability distributions of the physical properties (T eff , log (g), R, L and extinction) of USco 1610-1913 B, HIP 77900 B and HIP 78530 B using the BT-SETTL15 atmospheric models and the X-shooter spectra.We find that generally, the models fail to reproduce the pseudo-continuum of the X-shooter spectra over a broad range of wavelengths.Our solutions are mainly affected by the choice of the spectral range considered to estimate the best fit, and the fitting error bars remain about one to two orders of magnitude smaller given the high S/N of the X-shooter spectra.Finally, when the extinction is considered as a free parameter, the ForMoSA fitting solutions are considerably improved at all wavelengths with extinction values of 1.6-2.6 mag suggesting a clear deficiency in the dust grain modeling of the BT-SETTL15 atmospheric models.Notes. (a) Mean and standard deviation between J, H and K bands. (b) High value between J, H, K and JHK-cont bands. (c) Stefan−Boltzmann law in using T eff and radius ranges.

Fig. 2 .Fig. 3 .
Fig. 2. Comparison of the optical spectra normalized at 0.82 µm (left) and J-band normalized at 1.32 µm (right) of USco 1610-1913 B, HIP 77900 B, and HIP 78530 B (blue) to the free-floating object USco 1607-2211 from our original sample of spectra of Upper-Scorpius BDs and to the TW Hydrae association member TWA 26 (Manara et al. 2013) (black).All objects have a H α emission line (656.3nm; red label) in addition to molecular and atomic absorption lines typical of late-M dwarfs.

Fig. 4 .
Fig. 4. Chi-squared from the comparison between the X-shooter spectra of USco 1610-1913 B, HIP 77900 B, and HIP 78530 B to the optical spectra of the Ultracool RIZzo spectral library.We indicate the bestfitting object for each group in blue.The χ 2 minimum for all objects is represented in red.(a) : Cruz et al. (2003).(b) : Reid et al. (2008).

Fig. 5 .
Fig. 5. Comparison of the flux-calibrated 0.4-2.5 µm spectra of USco1610-1913 B (black)  to those of reference objects scaled to the distance of the companion (green).The companion is 2.0-3.9 times more luminous than objects with the same spectral type.

Fig. 6 .
Fig. 6.Posteriors of HIP 78530 B in using ForMoSA on the J + H + K and with the extinction parameter free.

Fig. 7 .
Fig. 7. Best fits of TWA 26 (in black) by BT-SETTL15 models (in blue) using different portions of the spectrum for the fit (shaded yellow rectangles).We define the fitting zones as follows: J = 1.0−1.32µm; H = 1.5−1.8µm; K = 2.0−2.48µm.We perform the following fits: (a) J + H + K; (b) J; (c) H; (d) K; (e) J + H + K with Av as a free parameter; (f ) the optical part (0.56−1.00 µm) of the flux-continuum and v sin (i) as a free parameter; and (g) examples of lines fitted on the flux-continuum.

Fig. 10 .
Fig. 10.Same as Fig. 7 but for HIP 78530 B. We avoid the flux inconsistency presented in Sect.4.1 by using the wavelength range 0.8-1.0µm in panel f.

Fig. 12 .
Fig.12.Correlation strength as a function of RV for each of our objects.We use the synthetic spectra at T eff = 2400 K and log (g) = 4.0 dex from BT-SETTL15.

Fig. 13 .
Fig. 13.Comparison of the UVB spectrum of USco 1610-1913 B (black) and TWA 26 (green) re-normalized to match the spectral continuum of the companion.Emission lines identified in the spectra are reported.
compares the UVB spectra of USco 1610-1913 B and TWA 26 scaled to the distance of USco 1610-1913 B. The flux from TWA 26 had to be multiplied by an additional factor of two to match the flux level of USco 1610-1913 B (see Sect. 6.2).The pseudo-continuum of USco 1610-1913 B is well reproduced by a template with emission lines due to chromospheric activity Fig. 14.H α line profiles for six of the targets with the strongest line emission.The dotted gray line corresponds to an estimate of the line spread function.
identified USco 1610-1913 B while searching for distant companions within 30 arcseconds of Upper-Sco stars and with optical (Pan-Starrs) and NIR (UKIDSS) colors compatible with young, cool objects.Using the now-available Pan-Starrs data (Chambers et al. 2016), we confirm that USco 1610-1913 B belongs to one of the three objects with the reddest i − z and i − y colors within 5 arcmin of USco 1610-1913 A. It is also the only object in this field with colors typical of M7-M9 Upper-Sco objects.Dahm et al. (2012) determined a RV of −6.91 ± 0.27 km s −1 for USco 1610-1913 A. USco 1610-1913 B is the unique target for which we find a strong discrepancy between the RV determined by ForMoSA A124, page 13 of 21 A&A 633, A124 (2020) (d) FromLallement et al. (2019)'s maps (derived fromGalli et al. 2018).

Table 4 .
Physical properties of TWA 26 and TWA 29.

Table C .
1. Results from fits with ForMoSA.

Table C .
2. Results from fits with ForMoSA.