An ALMA Dynamical Mass Estimate of the Proposed Planetary-mass Companion FW Tau C
Ya-Lin Wu Patrick D. Sheehan





Steward Observatory, University of Arizona, Tucson, AZ 85721, USA
Accepted for Publication in ApJL



Dynamical mass estimates down to the planet-mass regime can help to understand planet formation. We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of FW Tau C, a proposed ∼10 M_ Jup planet-mass companion at ∼330 au from the host binary FW Tau AB. We spatially and spectrally resolve the accretion disk of FW Tau C in ^12CO (2–1). By modeling the Keplerian rotation of gas, we derive a dynamical mass of ∼0.1 M_. Therefore, FW Tau C is unlikely a planet, but rather a low-mass star with a highly inclined disk. This also suggests that FW Tau is a triple system consisting of three ∼0.1 M_ stars. 






§ 1. INTRODUCTION

In the past ∼10 years, many planet-mass companions at wide separations (few to tens of M_ Jup; tens to hundreds of astronomical units from host stars) have been discovered in direct imaging surveys. Their masses are usually determined by comparing observables like luminosity and effective temperature to predictions from evolutionary and atmospheric models (e.g., ), which vary widely depending on the formation pathway (e.g., core accretion versus gravitational instability). As a result, it would be valuable if masses could be dynamically measured. As some of these substellar companions are young (≲10 Myr) and have features associated with active accretion (e.g., ), their dynamical masses can be measured if their accretion disks can be spatially and spectrally resolved.

Among all of the known wide-separation companions, the proposed planet-mass object FW Tau C is the prime target to search for a Keplerian-rotating disk. FW Tau C is a tertiary companion located at ∼23 projected separation (∼330 au) to FW Tau AB, a close binary (∼008) of nearly equal-mass stars (∼0.1 M_) in the 2 Myr Taurus–Auriga star-forming region. FW Tau C was discovered by <cit.> in their survey of binary stars, and its common proper motion was recently confirmed <cit.>. Studies have shown that FW Tau C has a rather flat near-infrared continuum owing to accretion-induced veiling, as well as many emission lines indicative of outflow and accretion activities <cit.>. The accretion disk was previously detected by the Atacama Large Millimeter/submillimeter Array (ALMA) in 1.3 mm dust continuum <cit.> and ^12CO (2–1) <cit.>, and a 1–2 M_⊕ of dust was inferred. Despite these observational efforts, FW Tau C's nature remains enigmatic because strong veiling inhibits accurate mass estimates. <cit.> derived a planetary mass of 10±4 M_ Jup from its dereddened K' flux, while <cit.> suggested that FW Tau C could be a 0.03–0.15 M_ brown dwarf or low-mass star embedded in an edge-on disk in order to explain its flat K-band spectrum and faint optical and near-infrared brightness.

Accurate mass measurement is therefore key to distinguishing the planet-mass scenario from the stellar-mass scenario. Here, we present the new ALMA 1.3 mm data from Cycle 3. With a ∼02 beam and ∼0.3  velocity resolution, we spatially and spectrally resolve the gas disk and derive a dynamical mass of FW Tau C by modeling the Keplerian rotation. Parameters of the FW Tau system are summarized in Table <ref>. 

@lccl@

Properties of FW Tau 

Parameter    
FW Tau AB    
FW Tau C    
References 



Distance (pc)					    	2c∼140		    	1		
 

Age (Myr)						    	2c∼2		    	2		
 

Separation ()				    	2c∼2.3		    	3, 4		
 

PA ()						    	2c∼296		   	3, 4		
 

SpT							   	M6 ± 1					     	⋯ 	    	3	
 

A_V (mag)					   	∼0.4						   	⋯ 	   	3	
 

log(L/L_) 				   	∼-1.1a		    	⋯ 		   	5	
 

Mass (M_)				   	∼0.12a		    	∼0.1		    	5, 6
^a For each component.
(1) <cit.>, (2) <cit.>, (3) <cit.>, (4) <cit.>, (5) <cit.>, (6) This work.












§ 2. METHODOLOGY



 §.§ 2.1. Observations and Data Reduction

FW Tau was observed with ALMA Band 6 during Cycle 3, on 2016 September 14. During the observations thirty-six 12 m antennas were available, with baselines ranging from 15 to 3247 m. The Band 6 receiver was configured to have three basebands set up for dust continuum observations, centered at 233.0, 246.0, and 248.0 GHz and each with 2 GHz of bandwidth. The last baseband was configured with 3840 0.122 MHz channels centered at 230.538 GHz (0.32  velocity resolution; Hanning smoothed) in order to spatially and spectrally resolve ^12CO (2–1) emission from the disk. J0510+1800 was used as the bandpass and flux calibrator, and J0433+2905 was used as the the gain calibrator. The on-source time was ∼13 minutes. 

The data were reduced using the ALMA pipeline in the  package. We employed one iteration of phase-only self-calibration to the 2 GHz continuum basebands. The baseband containing CO emission was much narrower, so it had much lower signal-to-noise ratio in continuum emission and we were unable to obtain good self-calibration solutions. We then CLEANed the calibrated data using the multi-frequency synthesis mode and natural weighting to enhance sensitivity in the images. The continuum map (left panel of Figure <ref>) has a beam size of 029×015 with a position angle (PA) of 1608, and an rms of 35 μJy beam^-1. The CO channel maps (top panel of Figure <ref>) have a beam size of 030×016 with a PA of 1603, and a mean rms of 4.5 mJy beam^-1 in signal-free channels. The integrated moment-zero and moment-one maps of the CO emission are shown in Figure <ref>.


@lc@

FW Tau C Disk Properties

Parameter    
Value 



M_* (M_⊙)     0.098 ± 0.015 
 

M_disk,dust (M_)     1.15 ± 0.01  
 

M_disk,gas (M_)     0.58±0.13  
 

R_disk (au)	     141.9 ± 6.6  
 

γ      1.63 ± 0.14  
 

T_0 (K)      62.1 ± 34.6  
 

q      0.14 ± 0.10  
 

ξ ()      0.65 ± 0.08  
 

v_ sys ()      6.10 ± 0.02 
 

i () 	      62.8 ± 2.4  
 

PA ()		      40.6 ± 1.5 





 §.§ 2.2. Disk Modeling: Continuum


The continuum emission from FW Tau C is, at best, only marginally resolved by our observations, so we use a simple geometrical model to fit the data. We assume that the continuum emission traces a uniform brightness disk, with a 1.3 mm brightness F_ν, a radius R_ disk, and some inclination i and position angle PA. We also allow the centroid of the emission to vary in our fit. We fit the model directly to the continuum visibilities using the Markov Chain Monte Carlo (MCMC) package  <cit.>.

In order to measure the dust mass of the system, we assume that the continuum emission traces optically thin dust so that we can estimate the disk mass from

    M_ disk = F_ν D^2/κ_ν  B_ν(T)

<cit.>. We use standard assumptions of T = 20 K and κ_ν = 2.3 cm^2 g^-1 and a distance to FW Tau C of 140 pc. 



 §.§ 2.3. Disk Modeling: Keplerian Rotation


Unlike the 1.3 mm continuum emission, FW Tau C's gas disk is well resolved by our ^12CO (2–1) observations (see Figure <ref>), including a clear detection of spatially resolved Keplerian rotation. To model the data and determine disk and stellar parameters, we follow the modeling procedure described in <cit.> to fit our channel maps with synthetic channel maps produced from radiative transfer models. Such models can be used to measure disk parameters such as radius and inclination, as well as the stellar (or planetary) mass <cit.>. Although we follow the procedure outlined by <cit.>, we have developed our own codes to run and fit these models to our data set.



We assume that the ^12CO (2–1) emission comes from a flared accretion disk, with a density profile described by

    ρ(R,z) = Σ(R)/√(2π)  h(R) exp[-1/2(z/h(R))^2],

where R and z are defined in cylindrical coordinates, and Σ(R) and h(R) are the surface density and disk scale height, respectively. We assume that the disk has a power-law surface density profile: 

    Σ(R) = Σ_0  (R/R_0)^-γ,

and we calculate the CO column density from this surface density profile as

    N_ CO(R) = X_ CO Σ(R)/μ  m_ H.

Here, X_ CO = 1 × 10^-4 is the CO mass abundance fraction, and μ = 2.37 is the mean molecular weight. The disk is truncated at an inner radius of 0.1 au and an outer radius of R_ disk, beyond which the density drops to zero.

We also assume that the disk is in local thermodynamic equilibrium and is vertically isothermal, with a radial temperature profile of 

    T(R) = T_0  (R/1  au)^-q.

Under these assumptions, the scale height in the disk is set by the balance of thermal pressure and gravity such that

    h(R) = ( k_b   R^3   T(R)/G   M_*  μ  m_ H)^1/2,

where k_b is the Boltzmann constant, and G is the gravitational constant.

Finally, rotation in the disk is assumed to be Keplerian, with an azimuthal velocity of

    v_k = √(G   M_*/r).

We assume that the velocities in the radial and vertical directions are zero. We also include microturbulent line broadening, which we assume is uniform throughout the disk, with a value of ξ in units of . Finally, we allow the star to have a systemic velocity, v_ sys, which Doppler shifts the velocity center away from zero.

In all, the density, temperature, and velocity structure of the system are described by the following parameters: M_ disk, R_ disk, γ, T_0, q, ξ, M_*, and v_ sys. We also allow the viewing geometry of the system, the inclination, and position angle to vary. In our model, the position angle of the disk is defined as the angle east of north of the projection of the disk angular momentum vector on the sky and ranges from 0^∘ to 360^∘. As our detection of CO (2–1) is not high sensitivity, we allow the inclination of the fit to range from 0^∘ to 90^∘. With higher sensitivity observations, however, it may be possible to distinguish between i<90^∘ and i>90^∘ <cit.>. Initially, we center the CO data based on the results of our geometrical continuum modeling, but allow for a small deviation from that centering in our model. 

We fix the distance to the disk at 140 pc (see Table 1). The uncertainty in the distance, however, translates linearly into an uncertainty on the measured stellar mass <cit.>, so we add the uncertainty on the distance estimate in quadrature with the uncertainty derived from our fit. Here, we assume a conservative distance uncertainty of 20 pc following <cit.>.

We use the 3D radiative transfer modeling package  <cit.> to calculate the level populations in each cell and produce synthetic ^12CO (2–1) channel maps for a given set of model parameters. Those synthetic channel maps are Fourier transformed and fit directly to the visibilities using the MCMC fitting package  <cit.>, with uniform priors for all parameters.

Although these models are computationally expensive, they can be run on powerful supercomputers so that the computations are spread out over a large number of central processing units (CPUs). For this particular instance, we run the fit over 128 CPUs on the University of Arizona El Gato supercomputer, and the modeling took a few days to converge. The models were determined to be converged when the  walkers had reached a steady state with measured best-fit values changing minimally over a large number of steps.



§ 3. RESULTS


Figure <ref> shows the 1.3 mm dust continuum, the ^12CO (2–1) integrated intensity map (moment 0), and the intensity-weighted velocity map (moment 1) of the FW Tau C disk. Similar to <cit.> and <cit.>, no signal is found from the close binary FW Tau AB, and the 3σ upper limit for a unresolved source suggests a dust mass ≲0.07 M_⊕. It is known that close binaries can shorten disk lifetimes (e.g., ), so the disk around FW Tau AB may be depleted already. For the companion FW Tau C, its dust disk is compact and likely unresolved; in contrast, the gas disk is more extended and clearly shows a Keplerian rotation.

Figure <ref> is the position–velocity diagram constructed along the major axis of the gas disk. On top of the PV diagram we also plot the Keplerian rotation curves for 10 M_ Jup and 0.1 M_ objects. It is clear that the velocity profile of the FW Tau C disk is incompatible with the 10 M_ Jup rotation curve, but more consistent with that of a 0.1 M_ star. Our disk modeling (see Figure <ref> and Table <ref>) also indicates that FW Tau C's mass is ∼0.1 M_, about 10 times higher than the 10 M_ Jup suggested by <cit.>. Hence, FW Tau C is not a planetary-mass object, but is rather a low-mass star with an inclined disk, as suggested by <cit.>. The high mass of FW Tau C is also in agreement with some features in its spectrum that closely resemble T Tauri stars <cit.>. Our observations therefore suggest that FW Tau is a young triple system in which each of the stars has a mass of ∼0.1 M_.

Searches for accretion disks in millimeter around directly imaged planet-mass companions have not yielded positive results (e.g., ), although the disk around the free-floating planet OTS 44 has recently been imaged <cit.>. 
Current flux upper limits imply that these wide companion disks might have <0.1 M_⊕ of dust, which in turn could imply a very short disk lifetime.
Future deep imaging down to a sublunar mass regime is perhaps needed to detect these disks, or place stringent constraints on their masses, and further elucidate the mass growth history of the planetary-mass companions at wide orbits. 

Finally, we summarize our disk modeling results as below and also in Table <ref>. Modeling the FW Tau C's dust disk as a circular Gaussian, we find that it has a 1.3 mm flux of ∼2.06 mJy, which equates to a dust mass of ∼1.15 M_⊕ and is consistent with previous measurements in <cit.> and <cit.>. For the gas disk, our radiative transfer modeling, as demonstrated in Figure <ref>, finds that the gas disk has a radius of ∼140 au, inclination of ∼63, PA of ∼41, systemic velocity of ∼6.1 , and gas mass of ∼0.58 M_⊕. Although recent studies have shown a low gas-to-dust ratio in protoplanetary disks (e.g., ), we caution that with a single optically thick CO line, we are very likely to underestimate the true gas mass by a factor of few or even orders of magnitude (e.g., ). As <cit.> pointed out, to accurately derive the gas mass, one should observe low-transition lines to ameliorate the optical-depth effects and should also use multiple CO isotopologues to characterize the chemical depletion and abundance variation across the disk.





§ 4. SUMMARY

FW Tau C is a wide-orbit companion at ∼330 au from the close binary FW Tau AB that has been suggested to have a planet-like mass of 10 M_ Jup. Here, we have used ALMA to detect FW Tau C's accretion disk in a 1.3 mm dust continuum and ^12CO (2–1) emission. We find that the gas motion is both spatially and spectrally resolved and clearly follows Keplerian rotation, enabling a dynamical mass estimate of the central object. We show that FW Tau C's mass is in fact ∼0.1 M_, so it is more likely a low-mass star embedded in an inclined accretion disk.  


We thank the referee for very helpful comments. We thank Laird Close, Johanna Teske, Yu-Cian Hong, and Jing-Hua Lin for discussions. Y.-L.W. is supported by the NASA Origins of Solar Systems award and the TRIF fellowship. This Letter makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00773.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.



[Andrews et al.(2013)]A13Andrews, S. M., Rosenfeld, K. A., Kraus, A. L., & Wilner, D. J. 2013, , 771, 129
[Ansdell et al.(2016)]A16Ansdell, M., Williams, J. P., van der Marel, N., et al. 2016, , 828, 46
[Baraffe et al.(2015)]Baraffe15Baraffe, I., Homeier, D., Allard, F., & Chabrier, G. 2015, , 577, A42
[Bayo et al.(2017)]B17Bayo, A., Joergens, V., Liu, Y., et al. 2017, , 841, L11
[Beckwith et al.(1990)]B90Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, , 99, 924
[Bowler et al.(2015)]B15Bowler, B. P., Andrews, S. M., Kraus, A. L., et al. 2015, , 805, L17
[Bowler et al.(2014)]B14Bowler, B. P., Liu, M. C., Kraus, A. L., & Mann, A. W. 2014, , 784, 65
[Bowler et al.(2011)]Bowler11Bowler, B. P., Liu, M. C., Kraus, A. L., Mann, A. W., & Ireland, M. J. 2011, , 743, 148
[Caceres et al.(2015)]C15Caceres, C., Hardy, A., Schreiber, M. R., et al. 2015, , 806, L22
[Chabrier et al.(2000)]C00Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 2000, , 542, 464
[Cieza et al.(2009)]Cieza09Cieza, L. A., Padgett, D. L., Allen, L. E., et al. 2009, , 696, L84
[Czekala et al.(2015)]Cz15Czekala, I., Andrews, S. M., Jensen, E. L. N., et al. 2015, , 806, 154
[Czekala et al.(2016)]Cz16Czekala, I., Andrews, S. M., Torres, G., et al. 2016, , 818, 156
[Dullemond(2012)]Dullemond12Dullemond, C. P. 2012, RADMC-3D: A Multi-purpose Radiative Transfer Tool, Astrophysics Source Code Library, ascl:1202.015
[Foreman-Mackey et al.(2013)]FM13Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, , 125, 306
[Isella et al.(2014)]I14Isella, A., Chandler, C. J., Carpenter, J. M., Pérez, L. M., & Ricci, L. 2014, , 788, 129
[Kenyon et al.(1994)]K94Kenyon, S. J., Dobrzycka, D., & Hartmann, L. W. 1994, , 108, 1872
[Kraus et al.(2015)]K15Kraus, A. L., Andrews, S. M., Bowler, B. P., et al. 2015, , 798, L23
[Kraus & Hillenbrand(2009)]K09Kraus, A. L., & Hillenbrand, L. A. 2009, , 704, 531
[Kraus et al.(2014)]K14Kraus, A. L., Ireland, M. J., Cieza, L. A., et al. 2014, , 781, 20
[Lynden-Bell & Pringle(1974)]LBP74Lynden-Bell, D., & Pringle, J. E. 1974, , 168, 603
[MacGregor et al.(2017)]M17MacGregor, M. A., Wilner, D. J., Czekala, I., et al. 2017, , 835, 17
[Marley et al.(2007)]M07Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2007, , 655, 541
[Ricci et al.(2017)]R17Ricci, L., Cazzoletti, P., Czekala, I., et al. 2017, , 154, 24
[Rosenfeld et al. (2013)]R13Rosenfeld, K. A., Andrews, S. M., Hughes, A. M., et al. 2013, , 774, 16
[Spiegel & Burrows(2012)]SB12Spiegel, D. S., & Burrows, A. 2012, , 745, 174
[White & Ghez(2001)]WG01White, R. J., & Ghez, A. M. 2001, , 556, 265
[Williams & Best(2014)]WB14Williams, J. P., & Best, W. M. J. 2014, , 788, 59
[Wolff et al.(2017)]Wolff17Wolff, S. G., Ménard, F., Caceres, C., et al. 2017, , 154, 26
[Wu et al.(2017)]W17Wu, Y.-L., Sheehan, P. D., Males, J. R., et al. 2017, , 836, 223
[Yu et al.(2017)]Y17Yu, M., Evans, N. J., Dodson-Robinson, S. E., Willacy, K., & Turner, N. J. 2017, , 841, 39
[Zhou et al.(2014)]Z14Zhou, Y., Herczeg, G. J., Kraus, A. L., Metchev, S., & Cruz, K. L. 2014, , 783, L17



