GQ Lup B: Substellar Companion Insights

• GQ Lup B is a young substellar companion orbiting the T Tauri star GQ Lup A at roughly 100 AU, with mass estimates between 10 and 40 MJup that blur the line between brown dwarfs and planetary-mass objects.
• It is characterized by active, variable accretion, a low-gravity, late-M spectral atmosphere with effective temperatures around 2600-2700 K, and clear infrared evidence of a compact circumplanetary disk.
• Orbital analyses reveal a moderately eccentric and misaligned orbit with magnetospheric accretion signatures, offering key insights into alternative formation pathways beyond traditional core accretion.

GQ Lup B, also written GQ Lup b in much of the literature, is a young substellar companion to the classical T Tauri star GQ Lup A at a projected separation of about $0.7''$, corresponding to roughly $100$–$110$ AU in the cited studies. It has been described as a brown dwarf companion and as a planetary-mass companion, with reported mass estimates including $\sim 10$–$36~M_{\mathrm{Jup}}$, $\sim 20$–$40~M_{\mathrm{Jup}}$, and $\sim 10$–$30~M_J$ in different analyses (Horstman et al., 2024, Alcalá et al., 2020, Cugno et al., 2024). Because it is young, accreting, directly imaged, and associated with circumplanetary material, GQ Lup B has become a key system for studies of wide-orbit substellar formation, disk evolution, and the possible formation of satellites (Stolker et al., 2021, Horstman et al., 2024).

1. System context and basic properties

GQ Lup B is bound to GQ Lup A, a K7 T Tauri star of approximately $1\,M_\odot$, and the system age is quoted in the literature as $100$0–$100$1 Myr, $100$2–$100$3 Myr, and $100$4 Myr in different works (Alcalá et al., 2020, Demars et al., 2023, Schwarz et al., 2016). The companion itself has been reported at a projected separation of $100$5, $100$6 AU, $100$7 AU, and $100$8 AU depending on the assumed distance and epoch of analysis (Alcalá et al., 2020, Demars et al., 2023, Wu et al., 2017). These values place it in the class of wide-separation directly imaged substellar companions for which conventional close-in planet-formation analogies are difficult to apply straightforwardly.

The broader architecture is also important. A wide companion, 2MASS J15491331-3539118, was reported at $100$9, or about $110$0–$110$1 AU, with a mass of $110$2, making the system a probable hierarchical triple (Alcalá et al., 2020). That outer component is sometimes discussed informally as GQ Lup C. The resulting A+B+wide-companion configuration has been used as a test case for mixed formation scenarios in which stellar components arise from turbulent core fragmentation while the substellar companion forms in a circumprimary disk, although alternative interpretations remain in the literature (Alcalá et al., 2020).

Published observational characterizations consistently identify GQ Lup B as young and accreting. Signs of accretion were already associated with strong emission lines, and later work emphasized that the accretion is not steady (Alcalá et al., 2020). This combination of youth, wide separation, and ongoing accretion is central to the system’s significance: it allows simultaneous investigation of orbital dynamics, atmospheric composition, accretion physics, and circumplanetary material in a regime intermediate between giant planets and brown dwarfs.

2. Orbital architecture and dynamical measurements

Astrometric monitoring established GQ Lup B as a comoving companion and produced the first significant detection of orbital motion in the system. Using VLT/NACO together with earlier literature and archival data, one study measured a decline in separation of $110$3 mas/yr and an increase in position angle of $110$4yr, rejecting the background-star hypothesis at $110$5 in both coordinates (Ginski et al., 2014). A Least-Squares Monte-Carlo analysis then found best-fitting eccentric solutions with $110$6 between $110$7 and $110$8, while still allowing circular orbits for restricted high-inclination configurations (Ginski et al., 2014).

High-resolution spectroscopy added a line-of-sight dynamical constraint. CRIRES observations yielded a barycentric radial velocity of $110$9 km/s for the companion, which, when combined with then-available astrometric solutions, ruled out circular orbits and long-period low-eccentricity solutions with $\sim 10$0 AU and $\sim 10$1 (Schwarz et al., 2016). That same study identified three classes of allowed orbits, including one near $\sim 10$2 AU, $\sim 10$3, and $\sim 10$4, as well as broader families extending to higher eccentricity (Schwarz et al., 2016).

Subsequent orbit fitting with a longer astrometric baseline favored a relatively low-eccentricity but strongly misaligned orbit. Combining data from 2004–2019, one analysis reported $\sim 10$5 AU, $\sim 10$6, $\sim 10$7 deg, and a mutual inclination with the circumstellar disk of $\sim 10$8 deg (Stolker et al., 2021). A later study incorporating GRAVITY astrometry and companion radial velocities refined the adopted solution to $\sim 10$9 AU, $36~M_{\mathrm{Jup}}$0, $36~M_{\mathrm{Jup}}$1 deg, and a peri-center distance of $36~M_{\mathrm{Jup}}$2 AU, with a mutual inclination of $36~M_{\mathrm{Jup}}$3 deg relative to the circumstellar disk and $36~M_{\mathrm{Jup}}$4 deg relative to the host-star spin axis (Venkatesan et al., 24 Sep 2025). In that framework, the combination of radial velocity and astrometry was explicitly used to break the inclination–eccentricity degeneracy that affects long-period astrometry-only solutions (Venkatesan et al., 24 Sep 2025).

The mutual inclination is computed with the standard relation

$36~M_{\mathrm{Jup}}$5

where $36~M_{\mathrm{Jup}}$6 is the angle between two planes (Stolker et al., 2021, Venkatesan et al., 24 Sep 2025). The recurring result across orbital studies is that GQ Lup B is not simply coplanar with the circumstellar disk of GQ Lup A. This has been one of the principal arguments against uncomplicated in-disk assembly followed by quiescent evolution.

3. Atmosphere, chemistry, and rotation

Spectroscopic work has characterized GQ Lup B as a late-type, low-gravity object. Medium-resolution optical spectroscopy with VLT/MUSE found an M9 spectral type, weak and narrow Na I and K I alkali doublets, enhanced VO features, and an absence of FeH absorption, all consistent with a low-gravity atmosphere (Stolker et al., 2021). Joint optical/near-infrared modeling in that study gave $36~M_{\mathrm{Jup}}$7 K, $36~M_{\mathrm{Jup}}$8–$36~M_{\mathrm{Jup}}$9, and $\sim 20$0, together with a visual extinction of $\sim 20$1 mag (Stolker et al., 2021). J-band analysis in a separate variability study reported $\sim 20$2 K, $\sim 20$3, $\sim 20$4, $\sim 20$5 mag, and an atmospheric radius of $\sim 20$6 (Demars et al., 2023).

High-resolution spectroscopy has repeatedly shown slow rotation. CRIRES K-band measurements detected CO at $\sim 20$7 and H$\sim 20$8O at $\sim 20$9, yielding $40~M_{\mathrm{Jup}}$0 km/s and a barycentric radial velocity of $40~M_{\mathrm{Jup}}$1 km/s (Schwarz et al., 2016). A later CRIRES$40~M_{\mathrm{Jup}}$2 atmospheric retrieval found $40~M_{\mathrm{Jup}}$3 km/s and $40~M_{\mathrm{Jup}}$4 km/s, consistent with the earlier result and reinforcing the description of GQ Lup B as a slow rotator (Picos et al., 3 Jan 2025). The interpretation advanced in the earlier spin study was that the low projected spin is most likely a consequence of youth: the object is still accreting material and angular momentum and is expected to spin up as it contracts (Schwarz et al., 2016).

The most detailed chemical analysis so far used CRIRES$40~M_{\mathrm{Jup}}$5 K-band spectra fitted with petitRADTRANS and reported H$40~M_{\mathrm{Jup}}$6O, $40~M_{\mathrm{Jup}}$7CO, $40~M_{\mathrm{Jup}}$8CO, HF, Na, Ca, and Ti in the atmosphere of GQ Lup B (Picos et al., 3 Jan 2025). In that analysis, $40~M_{\mathrm{Jup}}$9CO was detected at $\sim 10$0, the carbon abundance was $\sim 10$1, the fitted surface gravity was $\sim 10$2, and the elemental ratio was $\sim 10$3, consistent with the solar value (Picos et al., 3 Jan 2025). The carbon isotope ratio was measured as $\sim 10$4 for GQ Lup B, while the host star gave $\sim 10$5 after explicit modeling of strong stellar veiling (Picos et al., 3 Jan 2025). The similarity of the stellar and companion isotope ratios was interpreted as evidence for a shared material reservoir rather than a chemically distinct core-accretion outcome (Picos et al., 3 Jan 2025).

4. Accretion phenomenology and line variability

Accretion is among the defining observational properties of GQ Lup B. Strong H$\sim 10$6 emission is seen in optical data, and Pa$\sim 10$7 was detected in archival observations (Stolker et al., 2021). MagAO imaging detected the companion in H$\sim 10$8, $\sim 10$9, $30~M_J$0, and $30~M_J$1, but not in the 643 nm continuum, indicating that the shortest-wavelength emission is dominated by H$30~M_J$2 (Wu et al., 2017). The optical spectrum also shows Ca II triplet emission, reinforcing the accretion interpretation (Stolker et al., 2021).

The accretion rate is method-dependent. From the measured H$30~M_J$3 flux and an adopted companion mass of $30~M_J$4, one study derived $30~M_J$5 using accretion-shock modeling (Stolker et al., 2021). An earlier MagAO-based estimate using the empirical relation

$30~M_J$6

yielded $30~M_J$7 to $30~M_J$8 for GQ Lup B, while also noting previous measurements as high as $30~M_J$9, implying variability and/or extinction effects (Wu et al., 2017). The literature therefore converges on active but non-steady accretion, while differing in its absolute calibration.

Time-domain spectroscopy has made the variability explicit. Monitoring of the HI Paschen $1\,M_\odot$0 line at $1\,M_\odot$1m found moderate short-timescale variability, $1\,M_\odot$2, and much larger changes on longer baselines, with line-flux variability reaching $1\,M_\odot$3 over a decade (Demars et al., 2023). The earliest epoch showed equivalent widths up to $1\,M_\odot$4 and integrated line fluxes of $1\,M_\odot$5–$1\,M_\odot$6 W m$1\,M_\odot$7, whereas 2017 epochs reached equivalent widths down to $1\,M_\odot$8 and line fluxes of $1\,M_\odot$9–$100$00 W m$100$01 (Demars et al., 2023). The line profiles are typically blue-shifted by $100$02 to $100$03 km/s and can reach FWHM values up to $100$04 km/s (Demars et al., 2023).

These line properties were compared with both magnetospheric-accretion and shock models. For GQ Lup B, the observed Paschen $100$05 blue-shifts were found to be compatible with magnetospheric accretion, whereas shock-only models could not reproduce the blue-shifted profiles (Demars et al., 2023). The study therefore favored magnetospheric accretion for the companion, while also emphasizing that higher spectral resolution, $100$06, is required to disentangle the line-formation mechanisms fully (Demars et al., 2023). This places GQ Lup B within an accretion phenomenology that resembles classical T Tauri stars, but at the low-mass companion boundary.

5. Circumplanetary material and the disk around GQ Lup B

The presence and nature of circumplanetary material around GQ Lup B are constrained differently at millimeter and infrared wavelengths. ALMA observations at $100$07m found no significant continuum emission at the location of the companion and set a $100$08 upper limit of $100$09 mJy (MacGregor et al., 2016). Under optically thin isothermal assumptions with $100$10 K, this implied $100$11; for $100$12 K, the limit became $100$13 (MacGregor et al., 2016). An optically thick interpretation constrained the emitting radius to $100$14 AU (MacGregor et al., 2016). A separate ALMA 1.3 mm analysis also reported non-detection, with a dust-mass upper limit of $100$15–$100$16 depending on assumed disk size, and emphasized that the deeper $100$17m limit was $100$18 (Wu et al., 2017).

Infrared studies, however, identify a compact warm disk. Optical-to-mid-infrared spectral energy distribution analysis showed that the $100$19 color is $100$20 mag redder than field dwarfs of similar spectral type and that the $100$21, NB4.05, and $100$22 bands exhibit significant excess emission over atmospheric predictions (Stolker et al., 2021). That excess was fit by a blackbody with $100$23 K and $100$24, interpreted as continuum emission from small grains in a protolunar disk (Stolker et al., 2021). The same work described the disk as transitional and suggested that the large inner cavity could be associated with satellite assembly or pebble depletion (Stolker et al., 2021).

JWST/MIRI spectroscopy extended this picture into the mid-infrared. Over $100$25–$100$26m, the disk spectrum showed no silicate emission feature at $100$27–$100$28m, a result interpreted as evidence for significant grain growth, with $100$29m, and potentially dust settling (Cugno et al., 2024). If the emission is dominated by an inner wall, the best-fit dust temperature was $100$30 K and the inferred cavity radius was $100$31, larger than the expected sublimation radius $100$32 (Cugno et al., 2024). By contrast, geometrically thin disk models gave a cavity of $100$33 and outer disk radii of $100$34–$100$35 (Cugno et al., 2024). The disk fractional luminosity, $100$36–$100$37, was described as consistent with a young or primordial disk rather than a debris disk (Cugno et al., 2024).

Taken together, these results indicate that the circumplanetary environment is warm and infrared-bright but faint in the millimeter. The cited studies explicitly connect that combination to compactness, grain growth, and dust settling, rather than to the absence of circumplanetary material per se (Cugno et al., 2024, MacGregor et al., 2016). This wavelength dependence is central to the current interpretation of GQ Lup B as an accreting substellar companion with a highly evolved or compact circumplanetary disk.

6. Satellite formation, exomoon searches, and formation scenarios

GQ Lup B is one of the few substellar companions with a detected circumplanetary disk, and the inferred cavity in that disk has motivated explicit exomoon searches (Horstman et al., 2024). Using the Keck Planet Imager and Characterizer, which feeds a high-resolution spectrograph covering $100$38–$100$39 microns at $100$40, one study presented the first dedicated radial-velocity observations around a high-contrast, directly imaged substellar companion in order to search for exo-satellites (Horstman et al., 2024). Across 11 epochs, the best and median RV errors were $100$41–$100$42 m/s, probably limited by systematic fringing from transmissive optics within KPIC (Horstman et al., 2024). At that precision, the search was sensitive to exomoons $100$43–$100$44 the mass of GQ Lup B, taking $100$45, for separations between the Roche limit and $100$46, the extent of the inferred cavity (Horstman et al., 2024).

The same study simulated the planned HISPEC instrument for Keck and estimated that future exomoon sensitivity could improve by more than an order of magnitude, extending sensitivity to less massive satellites potentially formed within the disk itself (Horstman et al., 2024). It also modeled the ability of satellites to carve cavities in a circumplanetary disk, finding that satellite-to-planet mass ratios $100$47 can create observable cavities and reporting a maximum cavity size of $100$48 carved by a satellite (Horstman et al., 2024). The summary of that work also gave the relation

$100$49

to describe the largest cavity associated with a satellite (Horstman et al., 2024).

These exomoon results intersect directly with the broader formation debate. Scattering scenarios have often been treated skeptically because several analyses of the circumprimary disk reported no gaps or inner cavity that would indicate a massive inner scatterer (Wu et al., 2017). One ALMA study of GQ Lup A resolved a compact dust disk with a radius of $100$50 AU and found no gaps or inner cavity, arguing that GQ Lup B might have formed in situ via disk fragmentation or by prestellar core collapse (Wu et al., 2017). Another high-resolution ALMA analysis of the circumprimary disk found a gap at $100$51 AU and a tentative gap at 32 AU, but did not attribute these annular features conclusively to the wide companion; instead, embedded planets within the primary disk were considered plausible explanations (Long et al., 2020). Meanwhile, the moderate eccentricity and strong misalignment of GQ Lup B’s orbit, particularly in the most recent 3D fits, have been interpreted as more consistent with cloud fragmentation or a star-like formation pathway than with core accretion or disk-driven migration (Venkatesan et al., 24 Sep 2025).

No single scenario is unambiguously established across all studies. Disc fragmentation, gravitational collapse, and mixed hierarchical formation remain viable in different combinations; scattering is generally disfavored where the circumprimary disk morphology shows no suitable massive perturber, but dynamical complexity is not excluded (Ginski et al., 2014, Wu et al., 2017). The chemically similar $100$52 ratios of GQ Lup A and B, the solar-like C/O ratio of the companion, the compact circumplanetary disk, and the large orbital misalignment collectively push current interpretation toward a formation history closer to fragmentation than to standard core accretion (Picos et al., 3 Jan 2025, Venkatesan et al., 24 Sep 2025). A plausible implication is that GQ Lup B occupies a transitional empirical regime: dynamically a wide substellar companion, spectroscopically a young low-gravity accretor, and circumplanetarily a system in which the observational signatures of disk evolution and possible satellite formation can be studied directly.