Substellar Research Articles

Exploring the planetary-mass population in the Upper Scorpius association

We aim at identifying very low-mass isolated planetary-mass member candidates in the nearest OB association to the Sun, Upper Scorpius (145 pc; 5-10 Myr), to constrain the form and shape of the luminosity function and mass spectrum in this regime. We conducted a deep multi-band ($Y$=21.2, $J$=20.5, $Z$=22.0 mag) photometric survey of six square degrees in the central region of Upper Scorpius. We extend the current sequence of astrometric and spectroscopic members by about two magnitudes in $Y$ and one magnitude in $J$, reaching potentially T-type free-floating members in the association with predicted masses below 5 Jupiter masses, well into the planetary-mass regime. We extracted a sample of 57 candidates in this area and present infrared spectroscopy confirming two of them as young L-type members with characteristic spectral features of 10 Myr-old brown dwarfs. Among the 57 candidates, we highlight 10 new candidates fainter than the coolest members previously confirmed spectroscopically. We do not see any obvious sign of decrease in the mass spectrum of the association, suggesting that star processes can form substellar objects with masses down to 4-5 Jupiter masses.

Large Adaptive optics Survey for Substellar Objects (LASSO) at wide separations around young, nearby, low-mass stars with Robo-AO

Large Adaptive optics Survey for Substellar Objects (LASSO) at wide separations around young, nearby, low-mass stars with Robo-AO

Speckle Imaging Characterization of Radial Velocity Exoplanet Systems

Abstract We conducted speckle imaging observations of 53 stellar systems that were members of long-term radial velocity (RV) monitoring campaigns and exhibited substantial accelerations indicative of planetary or stellar companions in wide orbits. Our observations were made with blue and red filters using the Differential Speckle Survey Instrument at Gemini-South and the NN-Explore Exoplanet Stellar Speckle Imager at the WIYN telescope. The speckle imaging identifies eight luminous companions within 2″ of the primary stars. In three of these systems—HD 1388, HD 87359, and HD 104304—the properties of the imaged companion are consistent with the RV measurements, suggesting that these companions may be associated with the primary and the cause of the RV variation. For all 53 stellar systems, we derive differential magnitude limits (i.e., contrast curves) from the imaging. We extend this analysis to include upper limits on companion mass in systems without imaging detections. In 25 systems, we rule out companions with masses greater than 0.2 M ⊙, suggesting that the observed RV signals are caused by late-M dwarfs or substellar (potentially planetary) objects. On the other hand, the joint RV and imaging analysis almost entirely rules out planetary explanations of the RV signal for HD 19522 and suggests that the companion must have an angular separation below a few tenths of an arcsecond. This work highlights the importance of combined RV and imaging observations for characterizing the outer regions of nearby planetary systems.

Exploring the planetary-mass population in the Upper Scorpius association

ABSTRACT We aim at identifying very low-mass isolated planetary-mass member candidates in the nearest OB association to the Sun, Upper Scorpius (USco) (145 pc; 5–10 Myr), to constrain the form and shape of the luminosity function and mass spectrum in this regime. We conducted a deep multiband (Y = 21.2, J = 20.5, Z = 22.0 mag) photometric survey of 6 deg2 in the central region of USco. We extend the current sequence of astrometric and spectroscopic members by about two magnitudes in Y and one magnitude in J, reaching potentially T-type free-floating members in the association with predicted masses below 5 Jupiter masses, well into the planetary-mass regime. We extracted a sample of 57 candidates in this area and present infrared spectroscopy confirming two of them as young L-type members with characteristic spectral features of 10-Myr-old brown dwarfs. Among the 57 candidates, we highlight 10 new candidates fainter than the coolest members previously confirmed spectroscopically. We do not see any obvious sign of decrease in the mass spectrum of the association, suggesting that star processes can form substellar objects with masses down to 4–5 Jupiter masses.

Large adaptive optics survey for substellar objects around Young, Nearby, Low-mass Stars with Robo-AO

Large adaptive optics survey for substellar objects around Young, Nearby, Low-mass Stars with Robo-AO

Search for associations containing young stars (SACY)

Context. Nearby young associations offer one of the best opportunities for a detailed study of the properties of young stellar and substellar objects thanks to their proximity (<200 pc) and age (∼5−150 Myr). Previous works have identified spectroscopic (<5 au) binaries, close (5−1000 au) visual binaries, and wide or extremely wide (1000−100 000 au) binaries in the young associations. In most of the previous analyses, single-lined spectroscopic binaries (SB1) were identified based on radial velocities variations. However, this apparent variation may also be caused by mechanisms unrelated to multiplicity. Aims. We seek to update the spectroscopy binary fraction of the Search for Associations Containing Young stars (SACY) sample, taking into consideration all possible biases in our identification of binary candidates, such as activity and rotation. Methods. Using high-resolution spectroscopic observations, we produced ∼1300 cross-correlation functions (CCFs) to disentangle the previously mentioned sources of contamination. The radial velocity values we obtained were cross-matched with the literature and then used to revise and update the spectroscopic binary (SB) fraction in each object of the SACY association. In order to better describe the CCF profile, we calculated a set of high-order cross-correlation features to determine the origin of the variations in radial velocities. Results. We identified 68 SB candidates from our sample of 410 objects. Our results hint that at the possibility that the youngest associations have a higher SB fraction. Specifically, we found sensitivity-corrected SB fractions of 22−11+15% for ϵ Cha, 31−14+16% for TW Hya and 32−8+9% for β Pictoris, in contrast to the five oldest associations we have sampled (∼35−125 Myr) which are ∼10% or lower. This result seems independent of the methodology used to asses membership to the associations. Conclusions. The new CCF analysis, radial velocity estimates, and SB candidates are particularly relevant for membership revision of targets in young stellar associations. These targets would be ideal candidates for follow-up campaigns using high-resolution techniques to confirm binarity, resolve orbits, and, ideally, calculate dynamical masses. Additionally, if the results on the SB fraction in the youngest associations were confirmed, it could hint at a non-universal multiplicity among SACY associations.

TOI-519 b: A short-period substellar object around an M dwarf validated using multicolour photometry and phase curve analysis

Context. We report the discovery of TOI-519 b (TIC 218795833), a transiting substellar object (R = 1.07 RJup) orbiting a faint M dwarf (V = 17.35) on a 1.26 d orbit. Brown dwarfs and massive planets orbiting M dwarfs on short-period orbits are rare, but more have already been discovered than expected from planet formation models. TOI-519 is a valuable addition to this group of unlikely systems, and it adds towards our understanding of the boundaries of planet formation. Aims. We set out to determine the nature of the Transiting Exoplanet Survey Satellite (TESS) object of interest TOI-519 b. Methods. Our analysis uses a SPOC-pipeline TESS light curve from Sector 7, multicolour transit photometry observed with MuSCAT2 and MuSCAT, and transit photometry observed with the LCOGT telescopes. We estimated the radius of the transiting object using multicolour transit modelling, and we set upper limits for its mass, effective temperature, and Bond albedo using a phase curve model that includes Doppler boosting, ellipsoidal variations, thermal emission, and reflected light components. Results. TOI-519 b is a substellar object with a radius posterior median of 1.07 RJup and 5th and 95th percentiles of 0.66 and 1.20 RJup, respectively, where most of the uncertainty comes from the uncertainty in the stellar radius. The phase curve analysis sets an upper effective temperature limit of 1800 K, an upper Bond albedo limit of 0.49, and a companion mass upper limit of 14 MJup. The companion radius estimate combined with the Teff and mass limits suggests that the companion is more likely a planet than a brown dwarf, but a brown-dwarf scenario is a priori more likely given the lack of known massive planets in ≈ 1 day orbits around M dwarfs with Teff < 3800 K, and given the existence of some (but few) brown dwarfs.

Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B

Abstract We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The J, H, K1, and K2 spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 to 1800 K for effective temperature and from 3.0 to 5.0 for log(g), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of submicron dust particles, similar to dereddening approaches applied to the peculiar red field L dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from 4 yr of astrometric monitoring is consistent with a ∼30 yr period, an orbital inclination of 147°, and a semimajor axis of 10 au, well within the estimated disk inner radius of ∼50 au. As one of a very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion–disk interaction and provide a benchmark for substellar and planetary atmospheric study.

Evolution of spheroidal dust in electrically active sub-stellar atmospheres

Context. Understanding the source of sub-stellar polarimetric observations in the optical and near-infrared is key to characterizing sub-stellar objects and developing potential diagnostics for determining properties of their atmospheres. Differential scattering from a population of aligned, non-spherical dust grains is a potential source of polarization that could be used to determine geometric properties of the dust clouds. Aims. This paper addresses the problem of the spheroidal growth of dust grains in electrically activated sub-stellar atmospheres. It presents the novel application of a mechanism whereby non-spherical, elongated dust grains can be grown via plasma deposition as a consequence of the surface electric field effects of charged dust grains. Methods. We numerically solve the differential equations governing the spheroidal growth of charged dust grains via plasma deposition as a result of surface electric field effects in order to determine how the dust eccentricity and the dust particle eccentricity distribution function evolve with time. From these results, we determine the effect of spheroidal dust on the observed linear polarization. Results. Numerical solutions show that e ≈ 0.94 defines a watershed eccentricity, where the eccentricity of grains with an initial eccentricity less than (greater than) this value decreases (increases) and spherical (spheroidal) growth occurs. This produces a characteristic bimodal eccentricity distribution function yielding a fractional change in the observed linear polarization of up to ≈0.1 corresponding to dust grains of maximal eccentricity at wavelengths of ≈1 μm, consistent with the near infrared observational window. Order of magnitude calculations indicate that a population of aligned, spheroidal dust grains can produce degrees of polarization P ≈ 𝒪(10−2 − 1%) consistent with observed polarization signatures. Conclusions. The results presented here are relevant to the growth of non-spherical, irregularly shaped dust grains of general geometry where non-uniform surface electric field effects of charged dust grains are significant. The model described in this paper may also be applicable to polarization from galactic dust and dust growth in magnetically confined plasmas.

Direct Radio Discovery of a Cold Brown Dwarf

Abstract Magnetospheric processes seen in gas giants such as aurorae and circularly polarized cyclotron maser radio emission have been detected from some brown dwarfs. However, previous radio observations targeted known brown dwarfs discovered via their infrared emission. Here we report the discovery of BDR J1750+3809, a circularly polarized radio source detected around 144 MHz with the Low-Frequency Array (LOFAR) telescope. Follow-up near-infrared photometry and spectroscopy show that BDR J1750+3809 is a cold methane dwarf of spectral type T6.5 ± 1 at a distance of . The quasi-quiescent radio spectral luminosity of BDR J1750+3809 is ≈5 × 1015 erg s−1 Hz−1, which is over two orders of magnitude larger than that of the known population of comparable spectral type. This could be due to a preferential geometric alignment or an electrodynamic interaction with a close companion. In addition, as the emission is expected to occur close to the electron gyrofrequency, the magnetic field strength at the emitter site in BDR J1750+3809 is B ≳ 25 G, which is comparable to planetary-scale magnetic fields. Our discovery suggests that low-frequency radio surveys can be employed to discover substellar objects that are too cold to be detected in infrared surveys.

A Novel Survey for Young Substellar Objects with the W-band Filter. I. Filter Design and New Discoveries in Ophiuchus and Perseus

We present the design and implementation of a medium-band near-IR filter tailored for detecting low-mass stars and brown dwarfs from the summit of Maunakea. The W-band filter is centered at 1.45 μm with a bandpass width of 6%, designed to measure the depth of the H2O water absorption prominent in objects with spectral types of M6 and later. When combined with standard J and H photometry, the W-band filter is designed to determine spectral types to ≈1.4 subtypes for late-M and L dwarfs, largely independent of surface gravity and reddening. This filter’s primary application is completing the census of young substellar objects in star-forming regions, using W-band selection to greatly reduce contamination by reddened background stars that impede broad-band imaging surveys. We deployed the filter on the UH 88 inch telescope to survey ∼3 degree2 of the NGC 1333, IC 348, and ρ Ophiuchus star-forming regions. Our spectroscopic followup of W-band selected candidates resulted in the confirmation of 48 ultracool dwarfs with a success rate of 89%, demonstrating the efficacy of this new filter and selection method.

The search for disks or planetary objects around directly imaged companions: a candidate around DH Tauri B

Context. In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. While the majority of the sample is populated by objects discovered using radial velocity and transit techniques, an increasing number have been directly imaged. These planets and brown dwarfs are extraordinary sources of information that help in rounding out our understanding of planetary systems. Aims. In this paper, we focus our attention on substellar companions detected with the latter technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery would shed light on many unresolved questions, particularly with regard to their possible formation mechanisms. Methods. To reveal bound features of directly imaged companions, whether for point-like or extended sources, we need to suppress the contribution from the source itself. Therefore, we developed a method based on the negative fake companion technique that first estimates the position in the field of view (FoV) and the flux of the imaged companion with high precision, then subtracts a rescaled model point spread function (PSF) from the imaged companion, using either an image of the central star or another PSF in the FoV. Next it performs techniques, such as angular differential imaging, to further remove quasi-static patterns of the star (i.e., speckle contaminants) that affect the residuals of close-in companions. Results. After testing our tools on simulated companions and disks and on systems that were chosen ad hoc, we applied the method to the sample of substellar objects observed with SPHERE during the SHINE GTO survey. Among the 27 planets and brown dwarfs we analyzed, most objects did not show remarkable features, which was as expected, with the possible exception of a point source close to DH Tau B. This candidate companion was detected in four different SPHERE observations, with an estimated mass of ~1MJup, and a mass ratio with respect to the brown dwarf of 1∕10. This binary system, if confirmed, would be the first of its kind, opening up interesting questions for the formation mechanism, evolution, and frequency of such pairs. In order to address the latter, the residuals and contrasts reached for 25 companions in the sample of substellar objects observed with SPHERE were derived. If the DH Tau Bb companion is real, the binary fraction obtained is ~7%, which is in good agreement with the results obtained for field brown dwarfs. Conclusions. While there may currently be many limitations affecting the exploration of bound features to directly imaged exoplanets and brown dwarfs, next-generation instruments from the ground and space (i.e., JWST, ELT, and LUVOIR) will be able to image fainter objects and, thus, drive the application of this technique in upcoming searches for exo-moons and circumplanetary disks.

Considerations for atmospheric retrieval of high-precision brown dwarf spectra

ABSTRACT Isolated brown dwarfs provide remarkable laboratories for understanding atmospheric physics in the low-irradiation regime, and can be observed more precisely than exoplanets. As such, they provide a glimpse into the future of high-signal-to-noise ratio (SNR) observations of exoplanets. In this work, we investigate several new considerations that are important for atmospheric retrievals of high-quality thermal emission spectra of sub-stellar objects. We pursue this using an adaptation of the h y dra atmospheric retrieval code. We propose a parametric pressure–temperature (P–T) profile for brown dwarfs consisting of multiple atmospheric layers, parametrized by the temperature change across each layer. This model allows the steep temperature gradient of brown dwarf atmospheres to be accurately retrieved while avoiding commonly encountered numerical artefacts. The P–T model is especially flexible in the photosphere, which can reach a few tens of bar for T-dwarfs. We demonstrate an approach to include model uncertainties in the retrieval, focusing on uncertainties introduced by finite spectral and vertical resolution in the atmospheric model used for retrieval (∼8 per cent in the present case). We validate our retrieval framework by applying it to a simulated data set and then apply it to the HST/WFC3 (Hubble Space Telescope’s Wide-Field Camera 3) spectrum of the T-dwarf 2MASS J2339+1352. We retrieve sub-solar abundances of H2O and CH4 in the object at ∼0.1 dex precision. Additionally, we constrain the temperature structure to within ∼100 K in the photosphere. Our results demonstrate the promise of high-SNR spectra to provide high-precision abundance estimates of sub-stellar objects.

Techniques for Constraining the Abundance of Nonplanetary Substellar Objects

Abstract An interesting problem of relativistic astrophysics is the possible existence of non planetary Compact SubStellar Objects (CSSOs), such as black holes, but also non-standard objects such as boson stars. The challenge is to detect and distinguish them from ordinary planets. Here we present a new approach to this problem, which consists in detecting directly these CSSOs using exoplanet searches by transits, radial velocity measurements and astrometry and comparing the results to standard planet models.

Does gravitational radiation impact the stellar habitable zone?

Substellar objects of ultrashort periods (less than 50 min) orbiting central stars and many ultracompact star binary systems of even smaller periods (less than 10 min) have been optically detected in the last few years. Their gravitational luminosities can be of the same order or even much greater than their electromagnetic counterparts thereby also comprising a special class of targets for the Large Interferometer Space Antenna (LISA). In this paper, it is argued that the “gravitational illumination” of such systems may alter significantly the standard electromagnetic stellar habitable zone as long as a fraction of its gravitational energy is converted into heat. An analytical formula for the habitable zone is derived by assuming that the planet is simultaneously irradiated by electromagnetic and gravitational luminosities. The effect is quantified and some examples are discussed based on the available observations.

HST Survey of the Orion Nebula Cluster in the H2O 1.4 μm Absorption Band. I. A Census of Substellar and Planetary-mass Objects

Abstract In order to obtain a complete census of the stellar and substellar population, down to a few in the ∼1 Myr old Orion Nebula Cluster, we used the infrared channel of the Wide Field Camera 3 of the Hubble Space Telescope with the F139M and F130N filters. These bandpasses correspond to the 1.4 μm H2O absorption feature and an adjacent line-free continuum region. Out of 4504 detected sources, 3352 (about 75%) appear fainter than m 130 = 14 (Vega mag) in the F130N filter, a brightness corresponding to the hydrogen-burning limit mass ( ) at ∼1 Myr. Of these, however, only 742 sources have a negative F130M–F139N color index, indicative of the presence of H2O vapor in absorption, and can therefore be classified as bona fide M and L dwarfs, with effective temperatures T ≲ 2850 K at an assumed 1 Myr cluster age. On our color–magnitude diagram (CMD), this population of sources with H2O absorption appears clearly distinct from the larger background population of highly reddened stars and galaxies with positive F130M–F139N color index and can be traced down to the sensitivity limit of our survey, m 130 ≃ 21.5, corresponding to a 1 Myr old ≃3 planetary-mass object under about 2 mag of visual extinction. Theoretical models of the BT-Settl family predicting substellar isochrones of 1, 2, and 3 Myr down to ∼1 fail to reproduce the observed H2O color index at M ≲ 20 . We perform a Bayesian analysis to determine extinction, mass, and effective temperature of each substellar member of our sample, together with its membership probability.

A new set of atmosphere and evolution models for cool T–Y brown dwarfs and giant exoplanets

We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of 0.001–0.075 M⊙ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H–He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by ~1−2% at the stellar–substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the 3.5−5.5 μm flux window can be used to calibrate vertical mixing in cool T–Y spectral type objects.

Detection of Polarization due to Cloud Bands in the Nearby Luhman 16 Brown Dwarf Binary

Abstract Brown dwarfs exhibit patchy or spatially varying banded cloud structures that are inferred through photometric and spectroscopic variability modeling techniques. However, these methods are insensitive to rotationally invariant structures, such as the bands seen in Jupiter. Here, we present H-band Very Large Telescope/NaCo linear polarization measurements of the nearby Luhman 16 L/T transition binary, which suggest that Luhman 16A exhibits constant longitudinal cloud bands. The instrument was operated in pupil tracking mode, allowing us to unambiguously distinguish between a small astrophysical polarization and the ∼2% instrumental linear polarization. We measure the degree and angle of linear polarization of Luhman 16A and B to be p A = 0.031% ± 0.004% and ψ A = −32° ± 4°, and p B = 0.010% ± 0.004% and , respectively. Using known physical parameters of the system, we demonstrate that an oblate homogeneous atmosphere cannot account for the polarization measured in Luhman 16A, but could be responsible for that of the B component. Through a nonexhaustive search of banded cloud morphologies, we demonstrate a two-banded scenario that can achieve a degree of linear polarization of p = 0.03% and conclude that the measured polarization of the A component must be predominantly due to cloud banding. For Luhman 16B, either oblateness or cloud banding could be the dominant source of the measured polarization. The misaligned polarization angles of the two binary components tentatively suggest spin–orbit misalignment. These measurements provide new evidence for the prevalence of cloud banding in brown dwarfs while at the same time demonstrating a new method—complementary to photometric and spectroscopic variability methods—for characterizing the cloud morphologies of substellar objects without signs of variability.

Colour–colour and colour–magnitude diagrams for hot Jupiters

ABSTRACT We use ground-based and space-based eclipse measurements for the near-infrared (IR) bands (JHKs) and Spitzer 3.6- and 4.5-μm bands to construct colour–colour and colour–magnitude diagrams for hot Jupiters. We compare the results with previous observations of substellar objects and find that hot Jupiters, when corrected for their inflated radii, lie near the blackbody line and in the same region of the colour–magnitude diagrams as brown dwarfs, including low-gravity dwarfs that have been previously suggested as exoplanet analogues. We use theoretical emission spectra to investigate the effects of different metallicity, C/O ratios, and temperatures on the IR colours. In general, we find that while differences in C/O ratio and metallicity do correspond to different locations on these diagrams, the measurement errors are too large to use this method to put strong constraints on the composition of individual objects. However, as a class, hot Jupiters cluster around the location expected for solar metallicity and C/O ratio.

The TRENDS High-contrast Imaging Survey. VIII. Compendium of Benchmark Objects

Abstract The physical properties of faint stellar and substellar objects often rely on indirect, model-dependent estimates. For example, the masses of brown dwarfs are usually inferred using evolutionary models, which are age dependent and have yet to be properly calibrated. With the goal of identifying new benchmark objects to test low-mass stellar and substellar models, we have carried out a comprehensive adaptive optics survey as part of the TaRgetting bENchmark-objects with the Doppler Spectroscopy high-contrast imaging program. Using legacy radial velocity measurements from the High Resolution Echelle Spectrometer at Keck, we have identified several dozen stars that show long-term Doppler accelerations. We present follow-up high-contrast observations from the campaign and report the discovery of 31 comoving companions, as well as 11 strong candidate companions, to solar-type stars with well-determined parallax and metallicity values. Benchmark objects of this nature lend themselves to orbit determinations, dynamical mass estimates, and independent compositional assessment. This compendium of benchmark objects will serve as a convenient test group to substantiate theoretical evolutionary and atmospheric models near the hydrogen fusing limit.