Accepted to ApJ: October 26, 2018 Preprint typeset using LATEX style AASTeX6 v. 1.0

ALMA DETECTION OF EXTENDED MILLIMETER HALOS IN THE HD 32297 AND HD 61005 DEBRIS DISKS

Meredith A. MacGregor1,2, Alycia J. Weinberger1, A. Meredith Hughes3, D. J. Wilner4, Thayne Currie5, John H. Debes6, Jessica K. Donaldson1, Seth Redfield3, Aki Roberge7, Glenn Schneider8

1Department of Terrestrial Magnetism, Carnegie Institution for Science, 5241 Broad Branch Road NW, Washington, DC 20015, USA 2NSF Astronomy and Astrophysics Postdoctoral Fellow 3Department of Astronomy, Wesleyan University, Van Vleck Observatory, 96 Foss Hill Drive, Middletown, CT 06457, USA 4Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA 5National Astronomical Observatory of Japan, Subaru Telescope, National Institutes of Natural Sciences, Hilo, HI 96720, USA 6Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD, 21218, USA 7Exoplanets and Stellar Astrophysics Lab, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 8Steward Observatory, The University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA

ABSTRACT We present ALMA 1.3 mm (230 GHz) observations of the HD 32297 and HD 61005 debris disks, two of the most iconic debris disks due to their dramatic swept-back wings seen in scattered light images. These observations achieve sensitivities of 14 and 13 µJy beam−1 for HD 32297 and HD 61005, respectively, and provide the highest resolution images of these two systems at millimeter wavelengths to date. By adopting a MCMC modeling approach, we determine that both disks are best described by a two-component model consisting of a broad (∆R/R > 0.4) planetesimal belt with a rising surface density gradient, and a steeply falling outer halo aligned with the scattered light disk. The inner and outer edges of the planetesimal belt are located at 78.5 ± 8.1 AU and 122 ± 3 AU for HD 32297, and 41.9 ± 0.9 AU and 67.0 ± 0.5 AU for HD 61005. The halos extend to 440 ± 32 AU and 188 ± 8 AU, respectively. We also detect 12CO J= 2 − 1 gas emission from HD 32297 co-located with the dust continuum. These new ALMA images provide observational evidence that larger, millimeter-sized grains may also populate the extended halos of these two disks previously thought to only be composed of small, micron-sized grains. We discuss the implications of these results for potential shaping and sculpting mechanisms of asymmetric debris disks. Keywords: circumstellar matter — : individual (HD 32297, HD 61005) — submillimeter: plane- tary systems

1. INTRODUCTION Debris disks are often seen as a sign of previous successful formation in a system (see recent reviews including Hughes et al. 2018; Wyatt 2018). Planetesimals, larger bodies leftover from the planet formation process, collide within these disks and grind down into smaller and smaller dust grains that are observable via scattered light or thermal emission. It is expected that any interior should interact with this dusty material, sculpting it through arXiv:1812.05610v1 [astro-ph.EP] 13 Dec 2018 gravitational perturbations. In a few systems, we can link the geometry and location of a dust belt back to known planets. For example, the presence of a giant planet was inferred from a warp in the β Pictoris and later confirmed through direct imaging (Mouillet et al. 1997; Lagrange et al. 2010). Given the known interactions, we could potentially use the observed structure of debris disks to place constraints on the orbits and masses of that would otherwise be undetectable with current methods. The HR 8799 system with four directly imaged planets between two dust belts offers a template case (e.g., Su et al. 2009); Wilner et al. (2018) use the location of the outer disk’s inner edge determined from millimeter observations to place a constraint on the mass of the outermost directly imaged planet. In both β Pictoris and HR 8799, the planet mass estimates derived from disk structure provide a separate check on limits derived from evolution, atmospheric modeling, and dynamical stability (Lagrange et al. 2009; Currie et al. 2011). Unfortunately, this seemingly simple picture is complicated since structure in debris disks can be created through alternative, non-planetary mechanisms including stellar flybys, interactions with the (ISM, e.g., Artymowicz & Clampin 1997), recent massive collisions (e.g., Mazoyer et al. 2014), or 2 instabilities from dust and gas interactions (e.g., Lyra & Kuchner 2013). Scattered light imaging has revealed a population of edge-on debris disks that exhibit dramatic asymmetric structures including bowed disk midplanes (e.g., ‘wings’ such as HD 32297 and HD 61005) and extensions (e.g., ‘needles’ such as HD 15115). To date, the dominant mechanism at work shaping these systems remains the subject of debate. The leading theories involve interactions between disk material and either the ISM or an interior giant planet. In the first scenario, ram pressure from interstellar gas can remove bound or unbound grains from the disk (Debes et al. 2009; Maness et al. 2009). On the other hand, an eccentric (e > 0.2) giant planet orbiting within the disk can secularly perturb a narrow ring of parent bodies producing small dust grains that are thrown outwards via stellar radiation pressure (Esposito et al. 2016; Lee & Chiang 2016). A planet could also directly perturb grains over many orbits, forcing them to have high eccentricities and inclinations (Maness et al. 2009). Since they both show wing-like features in scattered light observations, the HD 32297 (A5V, < 30 Myr, 133 ± 1 pc) +0.03 and HD 61005 (G8V, 40 Myr, 36.56−0.06 pc, nicknamed ‘The Moth’) debris disks are a useful pair to explore the mechanisms responsible for shaping asymmetric debris disks. We note that the recent Gaia Data Release 2 (Gaia Collaboration et al. 2016, 2018) gives an updated distance for the HD 32297 system of 133 pc compared with the Hipparcos result of 112 pc (Perryman et al. 1997). The difference between the Gaia and Hipparcos distances is smaller for HD 61005 with values of 36.6 and 34.5 pc, respectively. Both of these disks were first resolved in scattered light at 1.1 µm with HST NICMOS (Schneider et al. 2005; Hines et al. 2007), and have been subsequently imaged at millimeter wavelengths with low angular resolution (Maness et al. 2008; Ricarte et al. 2013; Steele et al. 2016). We know that the small grains that dominate scattered light observations are more easily affected by non-gravitational forces such as stellar radiation and winds. Larger, millimeter-sized grains are expected to be less affected by these interactions, and therefore more reliably trace the underlying planetesimal distribution in these disks. Indeed, interactions with the ISM will preferentially strip small grains from the disk, but are unable to affect larger grains. Thus, high resolution observations at millimeter wavelengths are especially critical to our understanding of the physical mechanisms shaping the structure of these asymmetric debris disks. Here, we present new observations from the Atacama Large Millimeter/submillimeter Array (ALMA) with the highest angular resolution to date of both the HD 32297 and HD 61005 debris disks. Our new ALMA images reveal that despite differences in spectral type, both systems are best described by a two-component structure with (1) a parent body belt, and (2) an outer halo coincident with the scattered light disk. In Section 2, we present the details of our observations. In Section 3, we discuss the structure of the millimeter continuum emission for both disks (3.1), the modeling approach we take to characterize the observed visibilities (3.2), and the detection of 12CO J= 2 − 1 emission from the HD 32297 system (3.3). In Section 4, we compare our best-fit models to previous work (4.1). We also explore the significance of our results in the context of debris disk shaping (4.2) and stirring (4.3) mechanisms, and the origin of gas in debris disks (4.4). Finally, in Section 5, we present our conclusions.

2. OBSERVATIONS We observed both the HD 32297 and HD 61005 systems with ALMA in Band 6 (1.3 mm, 230 GHz, #2015.1.00633.S, PI: Alycia Weinberger). For HD 32297, ALMA executed two scheduling blocks (SBs) in Cycle 3 with 36 antennas in two configurations: (1) a 40-minute SB (22 minutes on-source) with baselines of 15 − 310 m on 2016 January 1, and (2) a 68-minute SB (44 minutes on-source) with baselines of 15 − 704 m on 2016 June 21. For HD 61005, compact observations (baseline lengths of 12 − 335 m) were taken previously in ALMA Cycle 1 (2012.1.00437.S, PI: David Rodriguez) on 2013 December 4 and 2014 March 20 with 27 and 31 antennas in the array, respectively. Analysis of these lower resolution observations was presented in Olofsson et al. (2016). ALMA executed an additional 68-minute SB (44 minutes on-source) in Cycle 3 with 36 antennas in the array and baselines extending to 650 m on 2016 June 18. Additional details (including precipitable water vapor and time on source) for all of these observations are listed in Table 1. The correlator for our new ALMA observations was set-up to maximize sensitivity to continuum dust emission. For HD 32297, we also covered the 12CO J= 2 − 1 line at 230.538 GHz with low spectral resolution (1.3 km s−1). We used four basebands centered at 215.5, 217.5, 230.538, and 232.538 GHz, and two polarizations (XX and YY). The baseband targeting the 12CO line had 3840 channels with a bandwith of 1.875 GHz. The three basebands targeting only continuum emission included 128 channels with a bandwidth of 2 GHz. The previous ALMA observations of HD 61005 included four basebands centered at 213, 215, 228, and 230.5 GHz. The first three basebands targeted continuum emission with 128 channels and 2 GHz bandwidth, while the last baseband targeted the 12CO J= 2 − 1 line using 3840 channels and a bandwidth of 1 GHz. Since no CO emission was detected in the earlier ALMA observations, we did not include the 12CO line in our more extended observations of HD 61005, and included only four continuum basebands 3 with a bandwidth of 2 GHz each centered at 224, 226, 240, and 242 GHz, in both the XX and YY polarizations. All of the raw data sets were calibrated with the ALMA pipeline using the CASA package (version 4.7.2). For HD 32297, observations of J0510+1800 were used for both bandpass and flux calibration. Observations of J0502+0609 were used to account for time-dependent gain variations due to instrumental and atmospheric effects. For HD 61005, observations of J0538-4405 were used for bandpass and flux calibration, while observations of J0747-3310 were used for gain calibration. The systematic uncertainty for absolute flux calibration is estimated to be < 10%. Visibility weights from ALMA are known to be inaccurate when reduced in versions of CASA prior to 4.4 (Ott & Kern 2017). This is not expected to be a problem with these data given that they were reduced using a more recent version of CASA. However, to be safe, we used the statwt task in CASA to recalculate the data weights according to their true scatter such that the weights are consistent with 1/σ2, where σ is the per-channel noise of the data, before proceeding with our analysis. The calibrated visibilities were averaged in 30 second intervals to reduce the total size of the measurement sets, resulting in a total number of visibility points of 591,384 and 613,038 for HD 32297 and HD 61005, respectively. To analyze the 12CO line emission from HD 32297, the continuum level was subtracted from the spectral window containing the line. Both continuum and line images were generated using the multi-frequency synthesis CLEAN algorithm in the CASA package.

3. RESULTS AND ANALYSIS 3.1. Continuum Emission Figure 1 shows our new ALMA 1.3 mm continuum images of HD 32297 (left) and HD 61005 (right) overlaid as contours on the previous Hubble Space Telescope (HST) STIS coronagraphic images of optical scattered light from Schneider et al. (2014). With robust = 0.5 weighting, the synthesized beam size is 000. 76 × 000. 51 (100 AU × 68 AU at 133 pc, position angle = −65◦) and the rms noise is 14 µJy beam−1 for HD 32297. For HD 61005, the synthesized beam size is 000. 49 × 000. 43 (18 AU × 16 AU at 36.6 pc, position angle = −71◦) and the rms is comparable, 13 µJy beam−1. For both sources, the millimeter continuum emission aligns well with the orientation of the disk major axis seen in the scattered light images. The HST images for these two disks are characterized by broad, wing-like structures extending to the northwest and southeast of the disk major axis for HD 32297 and HD 61005, respectively. While we do not detect such large wing structures for either source in our ALMA images, both do show indications of some extended emission (seen in the 3−12σ contours) beyond the disk ansae (the bright peaks seen on either side of the that are characteristic of limb brightening from an edge-on, optically thin disk), which we interpret as an outer halo. The presence of significant millimeter emission at large stellocentric distances is atypical for debris disks, most of which appear sharply truncated at an outer radius. For HD 61005, the extended halo emission matches well with the outer scattered light that is along the disk midplane. The structure of the halo is somewhat less apparent by-eye for HD 32297, given its much larger distance (113 pc compared to 36.6 pc for HD 61005) and thus lower resolution. Despite this, these new ALMA images indicate that both the HD 32297 and HD 61005 debris disks are best characterized by a two component model with (1) a planetesimal belt characteristic of most millimeter observations of debris disks and (2) an additional outer halo extending ∼ 200 − 400 radially. 3.2. Continuum Modeling Approach and Results We model the ALMA 1.3 mm continuum observations of both HD 32297 and HD 61005 with an optically thin, geometrically flat, axisymmetric belt model. In order to draw comparisons between the two systems and because their scattered light and millimeter images show similar structures, we adopt the same model for both disks. Model fits with only a single power law or Gaussian belt component (as we have used previously for other debris disks including AU Mic and  Eridani, see MacGregor et al. 2013, 2015) left significant residuals, most notably in the outer regions of the disk beyond the ansae. We interpret these residuals as evidence for an outer halo component not accounted for by simple models with symmetric inner and outer edges. Motivated by this apparent two-component structure of the ALMA continuum images (Figure 1) of HD 32297 and HD 61005, we adopt a parametric model with a radial temperature profile, T ∝ r−0.5, and a surface density profile described by a broken power law with Σ ∝ rx1 for x2 Rin

Figure 1. New ALMA images of HD 32297 (left) and HD 61005 (right) match closely to previous scattered light observations, and show evidence for an additional halo component not seen before in lower resolution millimeter images. The ALMA 1.3 mm continuum emission is overlaid as white contours on HST STIS images from Schneider et al. (2014). The extended halo emission is seen in the 2 − 6σ contours beyond the disk ansae (the bright peaks on either side of the star that are characteristic of limb brightening from an edge-on, optically thin disk). Contour levels in both panels are in steps of [−3, 3, 6, 12, 24, 36, 48, 60]× the rms noise of 14 and 13 µJy for HD 32297 and HD 61005, respectively. The dashed white ellipse in the lower left corner of both panels shows the beam size with robust= 0.5 weighting of 000. 76 × 000. 51 and 000. 49 × 000. 43, again for HD 32297 and HD 61005 respectively. discontinuity introduces an additional model parameter, and the continuous model fits demonstrate that this is not necessary. In addition, we fit for the belt geometry, namely the inclination, i, and position angle, PA. Although we explored more complicated models that allowed for different surface density profiles, gaps between the belt and halo components, or multiple rings, they did not yield substantially different results. In principle, other two-component parametric profiles, such as an asymmetric Gaussian, could be fitted to these data. We choose a broken power law as the simplest model justified by the sensitivity and resolution of our current dataset. To obtain best-fit values and uncertainties for all parameters, we adopt our previously described approach (e.g., MacGregor et al. 2013) and fit directly to the millimeter visibilities within a Markov Chain Monte Carlo (MCMC) framework making use of the emcee package (Foreman-Mackey et al. 2013). We compute synthetic visibilities for each model using the python package vis sample1, and then calculate a χ2 likelihood function, lnL = −χ2/2 to compare with the data. In order to fully explore the posterior probability distribution of the data conditioned on the model, we use 106 iterations (100 walkers and 10,000 steps each). To check for convergence, we compute the Gelman-Rubin statistic (Gelman & Rubin 1992) and ensure that R<ˆ 1.1 for all model parameters. We assume uniform priors, with limits that ensure each model is physical: Fbelt, Fhalo, Ftot > 0 and 0

1 vis sample is publicly available at https://github.com/AstroChem/vis_sample or in the Anaconda Cloud at https://anaconda.org/ rloomis/vis_sample.

6

4.0 HD 32297 6.0 HD 61005 3.0 4.0 2.0

2.0 1.0 real (mJy) real (mJy) real

0.0 0.0

-1.0 -2.0

0.5 0.5

0.0 0.0

-0.5 -0.5 imag.(mJy) imag.(mJy)

0 100 200 300 400 0 100 200 300 400

UV (kλ) UV (kλ)

Figure 3. The deprojected real and imaginary visibilities for HD 32297 (left) and HD 61005 (right). The 1.3 mm ALMA visibilities have been deprojected using the disk geometry determined by our model fits and averaged in bins of 10 kλ. In both plots, the solid blue line indicates the best-fit broken power law model, while the dashed red line shows a Gaussian fit to the data. The Gaussian model without a halo component fails to fit the data well for either disk, especially in the first null and second peak (50 − 250 kλ) of the real visibility function. For both sources, the imaginary visibilities show some low-amplitude structure indicative of asymmetric disk structure. been overlaid as contours on the intensity-weighted velocity (1st moment). The central velocity in the barycentric frame is 20.6±0.3kms−1, consistent with the measured stellar of ∼ 23 km s−1 (Torres et al. 2006). Converting to the LSRK frame, which accounts for the peculiar motion of the , yields a central velocity of 5.3 ± 0.1 km s−1. The channel maps show only the central nine channels in which emission is clearly detected above 3σ. The rms noise with robust = 0.5 weighting is 1.0 mJy beam−1. The peak and integrated intensity are 0.28 ± 0.003 Jy beam−1 and 1.02 ± 0.01 Jy km s−1 = 7.91 ± 0.08 × 10−21 W m−2, respectively, determined over the region above 3σ in the 0th moment map. No 12CO emission was detected from HD 61005 in the earlier Cycle 1 data, and Olofsson et al. (2016) −7 placed an upper limit on the CO gas mass assuming Local Thermodynamic Equilibrium (LTE) of 7.7 × 10 M⊕ (scaled to the new Gaia distance, see Section 1) for optically thin gas using these archival observations. In contrast, −6 Kral et al. (2017) predict a CO mass for HD 61005 in non-LTE of 1.8 × 10 M⊕. For HD 32297, we can calculate the total mass of CO probed by the J = 2 − 1 line transition (noted by the subscript 21 in the following equation) for optically thin line emission in LTE:

2 4π F21md MCO = (1) hν21 A21x2 −21 2 Here, F21 = 7.91 × 10 W m is the integrated flux density of the line, m is the mass of the CO molecule, d = 133 pc is the distance to the source, h is Planck’s constant, ν21 = 230.538 GHz is the rest frequency of the line, × −7 −1 N2 A21 = 6.91 10 s is the Einstein A coefficient (M¨uller et al. 2005), and x2 = Ntot is the fractional level population of the J = 2 rotational level (N2) relative to the total CO population (Ntot) evaluated at a characteristic temperature −4 of 30 K. This approach yields a total CO mass of MCO = 4.32 × 10 M⊕. If the CO were primordial and therefore −4 accompanied by H2, and assuming the canonical CO/H2 ratio of 10 typical for the ISM, the total mass of molecular hydrogen in HD 32297 would be ∼ 4.3 M⊕. We have made two significant assumptions in our analysis of the 12CO emission from the HD 32297 debris disk, which may need to be reevaluated: (1) the gas is optically thin and (2) the disk is in LTE. We can estimate the optical depth of the J = 2 − 1 line from the 3-D data cube following Matr`aet al. (2017) and again assuming LTE:

2 d F21 x1 τν = B12 − B21 (2) A21∆ν∆A x2  where d is the distance, F21 is the same as above, ∆ν is the line width, ∆A is the area of the emission on the sky 00 00 along the line of sight (equivalent to the area of the 0. 76 × 0. 51 synthesized beam), x1 and x2 are the fractional 7

12 Figure 4. We resolve the radial structure and rotation of the CO J= 2 − 1 emission from the HD 32297 system with ALMA. 12 (left) The CO J= 2 − 1 moment maps for HD 32297. Black contours show the 0th moment map (velocity-integrated intensity) 1 1 in steps of [3, 6, 12, 24, 36, 48, 60] × 1.0 mJy km s− beam− , the rms noise in the image, overlaid on the 1st moment (colors, intensity-weighted velocity). A scale bar for the first moment is shown in the upper right corner for reference. (right) The 12 1 corresponding channel maps for CO. Channels are 1.3 km s− wide with LSRK velocities labeled in the upper left corner of 1 1 each panel. Again, contours are in steps of [3, 6, 12, 24, 36, 48, 60] × 1.0 mJy km s− beam− . The leftmost panel in the bottom row shows the 000. 76 × 000. 51 synthesized beam size as a dashed ellipse. In both figures, the star symbol indicates the location of the central star, and the dashed gray line shows the best-fit position angle of the disk major axis determined from modeling the continuum visibilities, PA= 47◦. 9. populations of the lower and upper transition levels, and A21, B12, and B21 are the Einstein coefficients. To determine the fractional level populations, we assume a kinetic temperature of 30 K, appropriate for a belt at ∼ 80 AU. From this calculation, we obtain a peak optical depth of 0.3, reasonably consistent with optically thin emission. However, previous observations of gas in debris disks derive low excitation temperatures (e.g. K´osp´alet al. 2013; Flaherty et al. 2016; Hughes et al. 2017), suggesting that the assumption of LTE might not be appropriate. For non-LTE, a detailed radiative transfer model (e.g., the Line Modeling Engine or LIME, Brinch & Hogerheijde 2010) would be required to obtain accurate gas masses.

4. DISCUSSION We have made new high resolution maps at millimeter wavelengths of the HD 32297 and HD 61005 debris disks with ALMA at 1.3 mm (230 GHz). These two systems are notable for the dramatic swept-back wings evident in scattered light images. Surprisingly, our new ALMA images hint at a more extended disk structure in millimeter emission that aligns with the scattered light wings. Model fitting directly to the millimeter visibilities confirms that both disks are best described by a two-component structure with (1) a planetesimal belt, typical of debris disks observed at millimeter wavelengths, and (2) an extended outer halo, previously unseen in other millimeter images of these two systems with lower resolution. Single-component models with either a power law or Gaussian surface brightness profile leave significant residuals in the outer regions of the disk. Here, we compare our new observations to previous imaging studies, and discuss the implications for the origin of the halo component and stirring mechanisms. Given our detection of 12CO J= 2 − 1 emission from HD 32297, we also consider the properties of gas in debris disks.

4.1. Comparison to Previous Observations To place our new ALMA results into context, we compare our best-fit model structure and geometry for the HD 32297 and HD 61005 debris disks with previous imaging studies at shorter wavelengths. Optical scattered light observations trace small, micron-sized grains rather than the larger, millimeter-sized grains we expect to dominate emission at ALMA wavelengths. Comparing the disk geometry at these complementary wavelengths reveals the distribution of grain sizes within these disks. In all of our model fitting, we have adopted the new Gaia distances (discussed in Section 1). Thus, for the purposes of making a consistent comparison between our ALMA results and previous scattered light observations which assumed the Hipparcos distance, we have scaled all previous determinations of radii for both disks to the new Gaia distances. 8

4.1.1. HD 32297

Our model fits place strong constraints on the radial location of the HD 32297 debris disk. We note, however, that given the presence of residuals from the best-fit model (see Figure 2, it is likely that our quoted uncertainties are underestimates. The inner edge of the planetesimal belt is located at 78.5 ± 8.1 AU with an outer edge at 122 ± 3 AU. Esposito et al. (2014) detect a break in the surface brightness profile of the disk at ∼ 134 AU in 1.6 and 2.2 µm Keck NIRC2 images. Also using Keck NIRC2 coronagraphic imaging, Currie et al. (2012) see a similar break at ∼ 130 AU albeit with a more complicated surface brightness profile interior. Debes et al. (2009) again find a break in HST images, but with an asymmetric radial location of 109 and 134 AU on the NE and SW sides of the disk, respectively. In the ‘birth ring’ picture of debris disks, where small grains are produced through the collisional erosion of large planetesimals, we expect the outer edge of the planetesimal belt as traced by millimeter grains to line up with the break observed in the scattered light observations (e.g., Wilner et al. 2011). For HD 32297, there is good agreement between the outer edge of the planetesimal belt at 122 AU as seen by ALMA and the scattered light break observed by multiple telescopes at ∼ 134 AU. Fits to the spectral energy distribution (SED) in the far-infrared confirm this radial structure, with a best-fit model that includes a cold dust ring centered at 134 AU (Donaldson et al. 2013). Previous scattered light observations provide less information on the inner edge of the disk, mostly due to resolution challenges and large inner working angles. Moerchen et al. (2007) note evidence for a dust depletion interior to ∼ 85 AU in mid-infrared observations at 11.7 and 18.3 µm. Fitzgerald et al. (2007) find that the 11.2 µm emission exhibits a bi-lobed structure characteristic of a disk with an inner edge at ∼ 73 AU. Both of these inner edge locations are not well constrained with mid-infrared resolution of ∼ 30 − 50 AU at the distance of HD 32297, but the results are nonetheless similar to our best-fit model inner edge of ∼ 79 AU. Currie et al. (2012) note a plateau in the surface brightness profile of the disk in Keck/NIRC2 imaging, which they interpret as evidence for two dust populations: (1) located between ∼ 55 − 85 AU and (2) truncated at ∼ 133 AU. The radial locations of these two dust belts coincide with the inner edge of the planetesimal component and outer halo in our ALMA model fits. HD 32297 was also imaged in the early days of CARMA at 1.3 mm (Maness et al. 2008). Combining these low resolution millimeter images with Spitzer MIPS +20 photometry, Maness et al. (2008) determine an inner edge location of 85−10 AU, and also find evidence for two or three distinct grain populations within the disk. Thus, the geometry of the main planetesimal belt is consistent across scattered light and millimeter wavelength imaging. Beyond the outer edge of the planetesimal belt at 122 AU, the millimeter halo component extends to 440 ± 32 AU as determined from our ALMA observations. As discussed in Section 3.2, these fits could be influenced by the residual emission on the NW side of the disk since we do not account for vertical structure in our models. The uncertainty on this value is already quite large due to the low surface brightness in the outer regions of the disk in our model. The degeneracy between fitting for an emission gradient and the outer radius of a power law disk model is well known (e.g., Mundy et al. 1996). The first HST images of the HD 32297 disk showed scattered light emission from a nearly edge-on disk out to ∼ 480 AU (Schneider et al. 2005) and nebulosity that extends out to ∼ 2040 AU, deviating from the midplane disk axis by ∼ 15◦ on each side (Kalas 2005). Using new HST/STIS broadband optical coronagraphy, Schneider et al. (2014) find a similar extent of ∼ 1890 AU for the disk. It is clear that the scattered light emission from small grains has a much larger radial extent and deviation from the midplane axis than the millimeter emission probed by ALMA. The millimeter disk geometry is in relatively good agreement with previous measurements of the midplane disk geometry from scattered light observations. Our best-fit inclination and position angle are 83◦.6(+4◦.6, −0◦.4) and 47◦.9 ± 0◦.2, respectively. Schneider et al. (2014) find a position angle of 47◦.5 ± 1◦.0 from HST images, nearly identical to our result. Previous determinations of the disk inclination range from ∼ 80◦ (Kalas 2005) to ∼ 88◦ (Boccaletti et al. 2012), overlapping with our result. However, more recent scattered light images have preferred inclinations closer to ∼ 88◦ (e.g., Boccaletti et al. 2012; Currie et al. 2012), on the upper end of this range. We attribute this discrepancy to the 4σ residuals on the northwest side of the disk major axis seen in Figure 2. Our model fits favor lower inclinations in order to try and match this excess emission, and the asymmetric uncertainty on our reported best-fit inclination accounts for this discrepancy. From a simple χ2 analysis in the image plane, we estimate an additional 5% uncertainty towards higher inclinations, consistent with scattered light determinations. If we subtract a model with a higher inclination (i.e., 88◦), we are left with larger residuals on the northwest side of the disk in excess of 8σ, while no additional residuals are produced on the southeast side of the disk. This finding indicates that there is significant millimeter emission towards the northwest of the disk major axis coincident with the extended scattered light wings that our current symmetric model does not account for (first discussed in Section 3.2). 9

4.1.2. HD 61005 The best-fit model to the ALMA millimeter emission for HD 61005 has a planetesimal belt with an inner edge of 41.9 ± 0.9 AU and an outer edge of 67.0 ± 0.5 AU. Again, we note that the uncertainties quoted here are likely underestimates given the significant residuals from the best-fit model (see Figure 2). Gemini Planet Imager (GPI) total intensity observations at 1.6 µm probed down to projected separations of < 10 AU, revealing a sharp disk inner edge at ∼ 50 AU (Esposito et al. 2016), in agreement with the inner edge determined by ALMA of ∼ 42 AU. Using complementary Keck NIRC2 angular differential imaging of the disk at 1.2 − 2.3 µm, Esposito et al. (2016) note asymmetric breaks in the disk surface brightness profile beyond 55 and 65 AU on the NE and SW sides of the disk, respectively. As with HD 32297, the coincidence of the break observed in the scattered light surface brightness profile and the outer edge of the planetesimal belt determined by ALMA supports the conclusion that the observed small grains are produced through a collisional cascade occurring within a birth ring of planetesimals. In contrast to HD 32297, the HD 61005 system has been observed multiple times at millimeter wavelengths. Ricarte et al. (2013) first resolved the disk with the Submillimeter Array (SMA), and determined an inner edge for the disk ± +3.0 of 71 2 AU. A new analysis of the same SMA data by Steele et al. (2016) yielded a consistent result of 75.4−4.7 AU. Olofsson et al. (2016) fit models to the low resolution ALMA data included in this work (2012.1.00437.S, PI: David +6.1 Rodriguez), and obtain a similar value for the disk reference radius of 70.6−8.7 AU. Both of the SMA determinations place the inner disk edge farther out than our best-fit result, likely due to differences in model formulation and to the low angular resolution of these observations. The reference radius determined by Olofsson et al. (2016) indicates the location of a power law break in the disk model, and is thus most comparable to our planetesimal outer belt edge of ∼ 67 AU. None of these previous millimeter images show any indication of the extended wing structures seen in scattered light observations. The halo component, as newly detected by the higher resolution ALMA observations presented here, extends from the outer edge of the planetesimal belt at 67 AU to 188 ± 8 AU. Esposito et al. (2016) detect swept-back wings with a deflection angle of ∼ 22◦ in the KECK NIRC2 images, but they are outside of the GPI field of view. The wings extend from ∼ 62 to 135 AU on the NE side of the disk, and from ∼ 67 to 143 AU on the SW side. Buenzli et al. (2010) detect emission out to ∼ 150 AU with VLT/NACO, and HST images show extended emission at distances of & 200 AU from the star (Hines et al. 2007; Schneider et al. 2014). While the ALMA emission does not exhibit the significant deflection seen in scattered light observations, the extent of the halo component is similar at both wavelengths. The geometry of the planetesimal belt is also in good agreement. We determine a best-fit inclination of 85◦.6 ± 0◦.1 and position angle of 70◦.3 ± 0◦.1. Schneider et al. (2014) determine a nearly identical inclination and position angle for the disk major axis of 85◦.1 and 70◦.3, respectively.

4.2. On the Origin of the Millimeter Halo Component Since the wing-like structures of the HD 32297 and HD 61005 debris disks have previously only been seen in scattered light images, discussions of their origin have centered on mechanisms that affect predominantly small, micron-sized grains. The millimeter-sized grains that dominate emission at ALMA wavelengths are expected to closely trace the location of the parent bodies in the disk. Thus, we have typically assumed that these large grains should not populate extended halos. Few previous observations have challenged this assumption, with the exception of Marino et al. (2016) who detected low-level extended emission from HD 181327, a face-on debris disk, which they interpreted as either a halo or unresolved second ring. If the population of larger grains in HD 32297 and HD 61005 can be more extended, as our new ALMA observations suggest, the mechanisms invoked to produce the morphology of these disks must be re-evaluated. In this section, we present several potential mechanisms and discuss their applicability to millimeter grains. It is possible that the millimeter emission detected in the halo components of HD 32297 and HD 61005 could result from an abundance of micron-sized grains instead of millimeter-sized grains. To assess the likelihood of this scenario, we can estimate the mass of small grains needed in the halo to account for the observed millimeter flux. For optically 2 thin dust emission, the dust mass is Mdust = Fν D /(κν Bν (Tdust)), where Fν is the flux density, D is the distance (133 pc for HD 32297, 36.6 pc for HD 61005), κν is the dust opacity, and Bν is the Planck function at the dust temperature, Tdust. The main source of uncertainty in this calculation is the dust opacity. For millimeter-sized grains, 2 −1 a commonly used value in the literature is κν = 2.3 cm g at 1.3 mm (Beckwith et al. 1990). In comparison, Pollack et al. (1994) determine a dust opacity of 7.9 × 10−3 cm2 g−1 for their 300 µm radius composite (50%) grains and 4.4 × 10−3 cm2 g−1 for 3 µm grains with the same composition, several orders of magnitude smaller. In order to be able to compare to literature values for other debris disk systems, we adopt the typical value of 2.3 cm2 g−1 in our 10 calculations for millimeter-sized grains in the planetesimal belt, and 1.3 cm2 g−1 for micron-sized grains in the halo, determined using the same ratio between small and large grains as the Pollack et al. (1994) values. These values are similar to those shown in Andrews (2015) for grains that are a mixture of amorphous silicates, carbon, and vacuum with a filling factor of 0.5. We note, however, that if we instead used the Pollack et al. (1994) values, our calculated dust masses would be significantly higher. To obtain an estimate for the dust temperature, we assume radiative equilibrium and use the radial location from our best-fit models. Since small grains are likely hotter than the balckbody equilibrium temperature, we adopt a − − λ −β simplified emission efficiency, Qλ = 1 e ( 2πa ) , where a is the grain size and β is the power law index of grain emission efficiency at long wavelengths which we assume to be 0.7, comparable with measurements from other disks around main sequence stars (Williams et al. 2004). For HD 32297, the best fit inner edge of the planetesimal belt is 78.5 ± 8.1 AU and the start of the halo component is 122 ± 3 AU, which gives dust temperatures of ∼ 47 K and ∼ 38 K for each component, respectively. The total flux density of the HD 32297 planetesimal belt is 3.04 ± 0.21 mJy, which yields a total mass of millimeter-sized grains of 0.57 ± 0.07 M⊕. The halo component has a total flux density of 0.63 ± 0.26 mJy and a total mass of micron-sized grains of 0.29 ± 0.06 M⊕. For HD 61005, the planetesimal belt inner edge is at 41.9 ± 0.9 AU and the halo component starts at 67.0 ± 0.5 AU, implying dust temperatures of ∼ 37 K and ∼ 29 K, respectively. The total flux density of the planetesimal belt is 4.82 ± 0.29 mJy yielding a total mass of millimeter-sized grains of 0.088 ± 0.010 M⊕, and the total flux density of the halo is 1.54 ± 0.48 mJy for a mass of micron-sized grains of 0.070 ± 0.012 M⊕. We can then calculate the ratio of halo dust mass to planetesimal belt dust mass to be 0.51 and 0.79 for HD 32297 and HD 61005, respectively, a significant fraction for both systems. Recent work has shown that such high opacities may be unrealistic for micron-sized astro-silicate grains in the Rayleigh regime (Kataoka et al. 2014). If we assume a much lower opacity of 10−3 cm2 g−1, the implied mass fraction of the halo compared to the planetesimal belt is instead 800 and 500 for HD 32297 and HD 61005, respectively. Mass ratios of 0.5 − 800 (the full span of the opacity values we have considered here) are highly unphysical, and inconsistent with collisional evolution models. Thus, it is unlikely that a large enough mass of such grains will build up to explain the observed millimeter fluxes, disfavoring any proposed mechanisms that operate preferentially on small grains. A population of larger bodies with longer collisional lifetimes and higher opacities is needed in the halo of both disks. One mechanism that operates predominantly on small grains and that has been discussed extensively in the literature (e.g., Artymowicz & Clampin 1997; Debes et al. 2009; Maness et al. 2009; P´astor 2017) is interactions with the local interstellar medium (ISM), in which disk grains interact with gas in an ISM cloud, and are subsequently blown out due to gas drag or ram pressure stripping. Millimeter grains are not expected to feel a significant effect from interactions with the ISM. If the grains are highly porous as indicated for some debris disks (e.g., AU Mic, Graham et al. 2007), it is possible that ISM interactions could still play a role, but further investigation would be required to better evaluate this scenario. It has also been suggested that the observed halos could result from a density enhancement of small grains due to a recent massive collision, as has been invoked previously to explain asymmetric structure in the HD 15115 (Mazoyer et al. 2014) and β Pictoris (Dent et al. 2014) debris disks. Such an event would likely produce an asymmetric distribution of small grains and is expected to be transient (lasting for ∼ 1000 orbits or a few Myr at a radius of ∼ 50 AU, Jackson et al. 2014), neither of which is characteristic of the HD 32297 and HD 61005 systems. Interactions with massive planets are often invoked to explain observed disk structures, including asymmetries, warps, and offsets, and are able to sculpt both small and large grains. To date, there are no detected planets in either the HD 32297 or HD 61005 systems. However, scattered light imaging of both disks shows brightness asymmetries that could be indicative of eccentric structure. Currie et al. (2012) report that the southwest side of the HD 32297 debris disk is 50 − 100% brighter than the northeast side in Keck/NIRC2 coronagraphic images, which is mirrored in the early CARMA image at 1.3 mm (Maness et al. 2008). HST NICMOS and ACS observations of HD 61005 also show a brightness asymmetry, with the northeast side of the disk nearly twice as bright as the southwest (Hines et al. 2007). VLT NaCo (Buenzli et al. 2010) and SPHERE (Olofsson et al. 2016) observations of HD 61005 show an offset of the disk center to the west of the stellar position along the disk major axis giving an eccentricity for the disk of e = 0.045 ± 0.015 and e ∼ 0.1, respectively. The significant residuals that we see in our new ALMA observations (see Section 3.2 for discussion) further support the claim that the HD 61005 debris disk may be asymmetric. These ALMA observations do not detect the central star or any asymmetry between the two disk ansae for either HD 32297 and HD 61005, which would be required to place strong constraints on any disk eccentricity. Without detecting the star, we cannot determine whether there is an offset between the stellar position and the disk centroid. Without detecting a brightness asymmetry, we cannot make any claims of pericenter or apocenter glow. In an eccentric disk, the pericenter side of the disk is closer to the star and as a result glows more brightly at shorter wavelengths producing pericenter glow. At longer (millimeter and radio) wavelengths, the apocenter side of the disk appears brighter due 11 to a surface density enhancement (Pan et al. 2016; MacGregor et al. 2017). We did not attempt to model such a brightness asymmetry here, because it is not apparent in the residuals resulting from our simpler model. In principle, the argument of pericenter is unconstrained by these observations; the pericenter of the disk could be located along the observed disk major axis or be projected along the line of sight. To date, multiple models of secular perturbation by a planet have been proposed to explain the unusual morphology of the HD 32297 and HD 61005 debris disks. Esposito et al. (2016) consider a planet interior to the HD 61005 disk with an eccentricity of at least 0.2, are able to produce models that exhibit the same features as scattered light observations. A range of possible semi-major axes and masses are allowed, from an -mass planet farther from the star at 35 AU to a Jupiter-mass planet close-in at 5 AU. Lee & Chiang (2016) construct a more generalized debris disk model in which a possibly eccentric planet secularly perturbs a narrow ring of parent bodies producing small dust grains that are thrown outwards via stellar radiation pressure. The resulting scattered light disk morphology depends on several factors, including the planet eccentricity, disk inclination, and whether or not the dust grains are on apsidally aligned orbits. By varying these parameters, they are able to qualitatively reproduce the observed morphologies of both HD 32297 and HD 61005. To produce the ‘moth’-like shape of HD 61005, a highly eccentric planet (e ∼ 0.7) is needed along with grains on orbits that are apsidally aligned with the parent ring and with their apoastra directed towards the observer. The double wings of HD 32297 can be produced by a similar model, but with dust grains only launched at the parent body periastra. Both of these models treat only the small grains that dominate scattered light images, but their results could be extended in order to compare to millimeter images. Maness et al. (2009) discuss a slightly different process by which a planet could perturb the orbits of disk grains over time, slowly forcing them to have high orbital eccentricities and inclinations. It is also possible that an exterior planet could be responsible for sculpting these systems. In the models of Thebault et al. (2014), a planet exterior to an initially narrow ring produces a population of large grains outside of the planet’s orbit. Although these grains are eventually removed from the system, their lifetime is long enough to contribute to the optical depth of the disk at large radii. Our new ALMA observations indicate that millimeter-sized grains are present in the halos of the HD 32297 and HD 61005 debris disks. Since many of the previously proposed mechanisms operate predominantly on small grains, it is likely that interactions with planetary perturbers play a role in generating the wing structures seen at both optical and millimeter wavelengths. In the case of HD 32297, the high mass of gas within the disk may also contribute to the dynamics of the system. Ultimately, observations with higher angular resolution will reveal any asymmetric structure or disk centroid offsets that might be indicative of planetary sculpting in these disks. If an eccentric interior planet really is responsible for shaping these systems, we would expect to see additional substructure in the inner regions of the disk as a result of its gravitational interaction with disk material.

4.3. Constraints on Stirring Mechanisms from the Surface Density Gradient x Our model fits constrain the disk radial surface density profile, Σ ∝ r , with power law index x1 within the planetesimal belt and x2 for the halo. Both HD 32297 and HD 61005 exhibit a rising surface density profile for the planetesimal belt with x1 = 2.23 ± 0.61 and 2.87 ± 0.19, respectively. The surface density falls off steeply in the halo as x2 = −6.20 ± 0.30 for HD 32297 and −5.52 ± 0.14 for HD 61005. Resolved millimeter observations have placed constraints on the surface density gradients of a handful of other debris disks. Notably, a growing population exhibit rising surface density gradients (e.g., AU Mic and  Eridani, see MacGregor et al. 2015, 2013). Intriguingly, many of these disks appear radially broad like HD 32297 and HD 61005, which have fractional widths (∆R/R) of 0.43 ± 0.04 ± +0.13 +0.09 and 0.46 0.01, respectively. For AU Mic, ∆R/R = 1.28−0.01, and for  Eridani, ∆R/R = 0.31−0.13. In comparison, the fractional width of the confined Fomalhaut debris disk is only 0.10 ± 0.01 (MacGregor et al. 2017). Some of these systems may in fact exhibit more complicated radial structures. ALMA observations of the HD 107146 debris disk also indicate a rising surface density gradient, but with multiple concentric rings (Ricci et al. 2015). Measuring the surface density gradient of debris disks could provide a key observational constraint on the mechanisms operating in these systems to incite collisions. One possibility is that debris disks are ‘planet-stirred’ through the gravitational influence of giant planets interior to the disk (Mustill & Wyatt 2009). Kenyon & Bromley (2002, 2008) discuss an alternative ‘self-stirred’ model in which Pluto-sized bodies (radii ∼ 1000 km) are produced on the inner edge of the disk and trigger a collisional cascade. The growth time in a disk is proportional to P/Σ, where P ∝ r3/2 and Σ is the initial surface density. If the initial surface density decreases with radius, as is expected in a primordial disk, the growth time is shorter on the inner edge of the disk. For a disk with total dust mass comparable to HD 32297 and HD 61005, these models predict timescales of tens of Myr to form large bodies on the inner edge of a disk at tens of AU. On the outer edge of the disk (> 100 AU), the timescale is much longer, & 1 Gyr. Kennedy & Wyatt (2010) predict the resulting surface density profile in such a ‘self-stirred’ disk, finding Σ ∝ r+7/3. The ages of HD 32297 and 12

HD 61005 are estimated to be ∼ 30 − 40 Myr. At face value, this could allow enough time for larger bodies to form near the inner edge of the HD 61005 disk at ∼ 42 AU, but is not entirely consistent for HD 32297 with an inner edge nearly twice as far from the star at ∼ 79AU. There a few ways to solve this slight discrepancy in timescales. The age of the HD 32297 system is not well- constrained. Although stellar kinematics and the high disk luminosity suggest a young age (∼ 10 − 20 Myr), evolu- tionary tracks allow from ages up to ∼ 500 Myr (Rodigas et al. 2014), which may provide more time to form large bodies on the inner edge of the disk. In addition, recent work by Krivov & Booth (2018) considers the possibility that smaller (a few km to ∼ 200 km in size) planetesimals could effectively stir debris disks. Such bodies form much more rapidly in protoplanetary disks through pebble concentration than larger Pluto-sized bodies. If large bodies are in fact not needed to incite collisions in debris disks, both the HD 32297 and HD 61005 debris disks could be self-stirred. Of course, it is possible that an entirely different stirring mechanism may be at work in the gas-rich HD 32297 disk. There is some similarity between the surface density profile of this system and , a comparable ∼40 Myr-old gas-rich debris disk. Resolved ALMA observations of this disk are well-fit with a surface density profile that rises with a power law index of 1.44 in the inner part of the disk between 7 and 68 AU, and then falls off with an index of −2.5 out to 215 AU (Hughes et al. 2017). High resolution (sub)millimeter observations of a larger population of gas-rich debris disks are needed to better understand the stirring mechanisms in these systems.

4.4. Gas in the HD 32297 Debris Disk The HD 32297 system is one of three CO-rich debris disks around 30 − 40 Myr-old A-type stars, a group which also includes HD 21997 (A3V, 69.6 pc) and 49 Ceti (A1V, 57.1 pc). Our measured ALMA flux density makes HD 32297 one of the most CO-rich debris disks observed, comparable in integrated line flux to HD 21997 (K´osp´alet al. 2013) and 1.5× that of 49 Ceti (Hughes et al. 2017), after accounting for their different distances (see Mo´oret al. 2017). All three of these systems have a dust belt detected in addition to the observed gas emission. Although, for HD 21997, the gas disk appears to have an inner radius closer to the central star than the dust. In HD 32297, the dust and gas seem to be co-located, implying that both are likely created through collisions between the same planetesimals. To illustrate this, Figure 5 (left) shows the ALMA 12CO J= 2 − 1 emission overlaid as contours on the 1.3 mm continuum emission. Unlike the other two systems, HD 32297 is the only one to display millimeter gas and dust emission from an extended halo component. HD 32297 also has the latest spectral type of the group at A5V–A6V (often misquoted in the literature as A0V, see Debes et al. 2009, for further discussion) with a mass of ∼ 1.7 M and luminosity of ∼ 8 L , more comparable to the somewhat younger and more gas-poor β Pictoris system (A6V, 19.8 pc, 23 Myr). Figure 5 (right) shows the position-velocity diagram for the CO emission observed by ALMA. For reference, the expected Keplerian rotation curve is overlaid as a black dashed line, assuming a of 1.7 M , typical for an A5V–A6V star, and an inclination of 83◦.6, the best-fit from our continuum modeling. By eye, the rotation velocity of the gas at the peak locations appears to be ±2.5 km s−1 from the systemic velocity, below the ±4.4 km s−1 that is expected from Keplerian rotation for a ring at ∼ 80 AU (using a higher inclination of ∼ 88◦, as suggested by scattered light observations, makes the fit even worse). Although suggestive, this result is marginal given the low spectral resolution of our data. With a channel width of only 1.3 km s−1, the observed velocity difference is only ∼ 2σ. Intriguingly, sub-Keplerian rotation is also seen in the β Pictoris debris disk. By modeling CI and CII emission, (Cataldi et al. 2018) derive a dynamical mass of 0.8 M , significantly below the assumed stellar mass of 1.75 M (Crifo et al. 1997) but comparable to the dynamical mass derived from Na I emission (Olofsson et al. 2001). One explanation for sub-Keplerian rotation in debris disks could be radiation pressure, which acts opposite to a star’s gravity, making the star appear less massive. Mid to late A-type stars like HD 32297, however, do not emit significant UV radiation at the absorption wavelengths of CO, so this effect is likely have a negligible direct effect on neutral molecules. In the β Pictoris disk, CO is quickly converted to C and O through photodissociation on timescales of ∼ 50 . Fern´andez et al. (2006) suggest that these ionized species, which are directly affected by radiation pressure, are coupled into a single fluid through Coulomb interactions. As a result, neutral molecules like CO could also become coupled and experience indirect braking due to radiation pressure. Given their similar spectral types, it is possible a similar explanation could be invoked in the case of HD 32297. An alternative explanation was originally presented in Weidenschilling (1977), where a pressure gradient in a gas disk causes the gas to rotate more slowly than the free orbital velocity. The deviation of the gas velocity from the Keplerian orbital velocity, vk, is given by − ≈ − 2 − ∆v = vk vg (∆g/2g)vk, where g = vk/r is the central gravity and ∆g = nRT/µr is the residual gravity in a −n reference frame rotating with the gas. Here, n is the power law index of pressure, P = P0(r/r0) , which we take to be 3, R is the gas constant, T is the gas temperature, and µ is the mean molecular weight. Given the expected stellar mass (1.7 M ), gas temperature (30 K), and radius of the HD 32297 debris disk, we estimate the magnitude of this 13

13

80

9 60 ] ] 1 1 − −

5 40 [km s LSRK [mJybeam ν v M = 1.7 M

⊙ I 1 20

0 -3 −2 −1 0 1 2 Offset ["]

12 Figure 5. The dust and CO J= 2 − 1 gas emission originate from the same location in the HD 32297 disk at ∼ 80 AU, but the gas exhibits slower velocities than expected for Keplerian rotation. (left) The 1.3 mm continuum emission (color) with the 1 1 0th moment map overlaid as black contours in steps of [3, 6, 12, 24, 36, 48, 64] × 1.0 mJy km s− beam− , the rms noise in the 12 image, overlaid on the 1st moment. (right) The position-velocity diagram of the CO J= 2 − 1 emission from the HD 32297 1 1 debris disk. Contours are in steps of [3, 6, 12, 24, 36, 48, 64] × 1.0 mJy km s− beam− , the rms noise in the image. The expected ◦ Keplerian velocity profile for a disk inclination of 83. 6 (the best-fit to the dust continuum emission) and a stellar mass of 1.7 M (typical for an A5V–A6V star) is shown by the dashed black line. effect to be ∼ 0.1 − 0.5km s−1, not large enough to account for the ∼ 2km s−1 observed difference in orbital velocity. We note, however, that modifying the pressure gradient and mean molecular weight to account for a primordial or secondary origin of the gas could change this calculation significantly. Although we do not discuss them at length here, two other possible explanations for the observed sub-Keplerian rotation of the gas could include non-uniform self-absorption or a stellar dynamical mass that disagrees with pre-main sequence models. Circumstellar gas has been detected previously from the HD 32297 system. Redfield (2007) observed absorption from the Na I D doublet (5895.9242 and 5889.9510 A)˚ with a column density of logNNaI ∼ 11.4, the strongest Na I absorption measured toward any nearby debris disk. A 3.7σ detection of [CII] emission at 158 µm was also made with the Herschel PACS Spectrometer (Donaldson et al. 2013). Most recently, Greaves et al. (2016) detected 12CO J= 2 − 1 and tentatively J=3 − 2 emission with JCMT SCUBA-2 (taken several years earlier in 2013). Their measured 12 −21 −2 −3 integrated intensity for CO J= 2−1 is 3.7±0.4×10 Wm , giving a total CO mass of 1.3×10 M⊕. This result is 5× (7σ) higher than our measured ALMA integrated intensity for the same transition of 7.91 ± 0.08 × 10−21 Wm−2 −4 with corresponding CO mass of 4.32 × 10 M⊕. The difference between these two measurements could easily stem from the lower signal-to-noise and resolution of the earlier JCMT observations. However, Greaves et al. (2016) also found a higher velocity for the gas, which implies a 2× difference in the inner radius of the gas. Given this, the CO production rate and location in the HD 32297 disk could be variable on the ∼ 3 timescale between these two observations, with a higher production rate in 2013, although this would be challenging to explain given that the time between observations is much shorter than the orbital timescale for the system. Given the relative wealth of gas detections from the HD 32297 debris disk, it is clear that this system offers us a unique opportunity to explore the properties of gas in main sequence circumstellar disks. By combining our determinations −3 of the total gas mass (4.3 × 10 M⊕ for CO and 4.3 M⊕ for H2, calculated by assuming an ISM-like CO/H2 ratio of −4 10 expected only for a primordial system, see Section 3.3) and the total dust mass (0.72 ± 0.13 M⊕, Section 4.2), we can calculate the gas-to-dust ratio for the HD 32297 debris disk to be 6 × 10−4 if CO dominates the gas mass (possibly higher if the gas is non-LTE, but still < 1) or 6 if H2 is present in the disk. As expected for an evolved system, the result assuming a primordial disk is significantly below the canonical value for the ISM of ∼ 100 (Bohlin et al. 1978), but is comparable to recent values determined from surveys of protoplanetary disks in much younger star forming regions (e.g., Ansdell et al. 2016; Miotello et al. 2017). Both results are also comparable to 49 Ceti, which has a gas to dust-ratio of either 0.008 or 6, depending on whether CO or H2 makes the largest contribution to the gas mass, respectively (Hughes et al. 2017). K´osp´alet al. (2013) calculate a much higher gas-to-dust ratio for HD 21997 14 of 300 − 700. This difference points to a hybrid scenario for HD 21997, where some primordial gas remains in the system, but dust is generated through secondary collisions between planetesimals. For 49 Ceti and HD 32297, where the dust and gas are co-located, it seems more likely that both are created through secondary collisions between icy planetesimals. Indeed, Kral et al. (2017) show that the mass of CO gas detected around HD 32297 is completely consistent with a secondary origin when considering non-LTE. If this is the case, it may be illustrative to compare the CO-to-dust ratio in HD 32297 with Solar System comets; recent results from the Philae lander suggest that Comet 67P/Churyumov-Gerasimenko is more dust than ice with a volumetric ice-to-dust ratio of 0.4 − 2.6 (Kofman et al. 2015). Additional observations with higher spectral resolution will be critical to confirm the origin of the gas in the HD 32297 system, as well as to examine the somewhat confusing gas dynamics and whether or not the CO production rate is indeed variable.

5. CONCLUSIONS We present new ALMA observations of the HD 32297 and HD 61005 debris disks at 1.3 mm (230 GHz). These images give us the highest angular resolution view to date of these two edge-on systems at millimeter wavelengths. To analyze these data, we adopt an MCMC approach and fit directly to the millimeter visibilities. The main results of this analysis are as follows.

1. The 1.3 mm continuum emission from HD 32297 and HD 61005 is best described by a two-component power law model that includes both a planetesimal belt and an additional outer halo. Simpler models with either a single power law or Gaussian leave significant residuals in the outer regions of the disk. For HD 32297, the inner edge of the planetesimal belt is located at 78.5 ± 8.1 AU with an outer edge at 122 ± 3 AU, and a halo out to 440 ± 32 AU. For HD 61005, the planetesimal belt extends between 41.9 ± 0.9 AU and 67.0 ± 0.5 AU, with a halo out to 188 ± 8 AU. Although the best-fit models provide reasonably good fits to the data for both disks, we note some significant residuals at the 4 − 8σ level that could indicate additional asymmetric substructure not resolved by our current observations and likely make our quoted uncertainties underestimates.

2. Our ALMA observations detect millimeter emission from an extended halo component in a debris disk. Pre- viously, the halos of the HD 32297 and HD 61005 debris disks were only seen in scattered light observations, implying that they were populated by small, micron-sized grains easily removed by interactions with the ISM or an interior planet. If, instead, millimeter grains are also present in these halos, we must reevaluate how such structures are created. Since interactions with the ISM are only expected to affect small grains, these new observations favor sculpting via planetary perturbations. It is also likely that the presence of gas in a disk, as is the case for HD 32297, will have an effect on the system’s dynamics, since outflowing gas can entrain grains.

3. For both HD 32297 and HD 61005, the planetesimal belt component has a rising surface density gradient with power law indices of x1 = 2.23 ± 0.61 and 2.87 ± 0.19, respectively. Similar to other systems that exhibit rising surface density gradients (e.g., AU Mic and  Eridani), both belts are broad with fractional widths, ∆R/R > 0.4. These results contribute to a growing population of broad debris disks that exhibit rising surface density gradients, consistent with predictions of ‘self-stirred’ models for debris disks where ongoing formation of small bodies at the inner edge of the disk triggers a collisional cascade.

4. We detect 12CO J= 2 − 1 emission from the HD 32297 debris disk co-located with the dust continuum emission. The integrated intensity is 7.91 ± 0.08 × 10−21 W m−2. Assuming that the gas is optically thin and in LTE, this −4 flux implies a total CO mass of 4.32 × 10 M⊕, comparable to or more massive than other CO-rich debris disks. Intriguingly, the measured ALMA integrated intensity is a factor of 5× lower than previous measurements for the same system (Greaves et al. 2016). Although this difference could stem from the lower signal-to-noise and resolution of the earlier observations, it could also suggest variability in the CO production rate in this system. In addition, the CO gas appears to exhibit sub-Keplerian rotation, although the significance of this result is unclear given the low spectral resolution of our data (channel width of 1.3 km s−1).

The ALMA observations presented here provide new insights into the structure and dynamics of the HD 32297 and HD 61005 debris disks. Ultimately, millimeter observations with higher angular and spectral resolution are required to confirm the dominant mechanism responsible for sculpting these systems. A detection of additional substructure in either disk could reveal the presence and influence of a giant interior or exterior planet. 15

The authors thank David Rodriguez for leading the Cycle 1 ALMA proposal to observe HD 32297, Ian Czekala for helpful discussions regarding the HD 32297 gas detection, Gaspard Duchˆene for helpful discussions regarding the inclination of HD 32297, and the anonymous referee for comments that improved the clarity of the manuscript. M.A.M. acknowledges support from a National Science Foundation Astronomy and Astrophysics Postdoctoral Fellowship under Award No. AST-1701406. A.M.H. acknowledges support from NSF grant AST-1412647. This paper makes use of the following ALMA data: ADS/JAO.ALMA #2015.1.00633.S and #2012.1.00437.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and 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. This work has also made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/ dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

REFERENCES

Andrews, S. M. 2015, PASP, 127, 961 Hines, D. C., Schneider, G., Hollenbach, D., et al. 2007, ApJL, Ansdell, M., Williams, J. P., van der Marel, N., et al. 2016, ApJ, 671, L165 828, 46 Hughes, A. M., Duchene, G., & Matthews, B. 2018, ArXiv Artymowicz, P., & Clampin, M. 1997, ApJ, 490, 863 e-prints, arXiv:1802.04313 Ballering, N. P., Rieke, G. H., Su, K. Y. L., & Montiel, E. 2013, Hughes, A. M., Lieman-Sifry, J., Flaherty, K. M., et al. 2017, ApJ, 775, 55 ApJ, 839, 86 Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. Jackson, A. P., Wyatt, M. C., Bonsor, A., & Veras, D. 2014, 1990, AJ, 99, 924 MNRAS, 440, 3757 Boccaletti, A., Augereau, J.-C., Lagrange, A.-M., et al. 2012, Kalas, P. 2005, ApJL, 635, L169 A&A, 544, A85 Kataoka, A., Okuzumi, S., Tanaka, H., & Nomura, H. 2014, Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 A&A, 568, A42 Brinch, C., & Hogerheijde, M. R. 2010, A&A, 523, A25 Kennedy, G. M., & Wyatt, M. C. 2010, MNRAS, 405, 1253 Buenzli, E., Thalmann, C., Vigan, A., et al. 2010, A&A, 524, L1 Kenyon, S. J., & Bromley, B. C. 2002, ApJL, 577, L35 Cataldi, G., Brandeker, A., Wu, Y., et al. 2018, ArXiv e-prints, —. 2008, ApJS, 179, 451 arXiv:1804.06479 Kofman, W., Herique, A., Barbin, Y., et al. 2015, Science, 349, Chen, C. H., Mittal, T., Kuchner, M., et al. 2014, ApJS, 211, 25 http://science.sciencemag.org/content/349/6247/aab0639.full.pdf Crifo, F., Vidal-Madjar, A., Lallement, R., Ferlet, R., & K´osp´al, A.,´ Mo´or, A., Juh´asz, A., et al. 2013, ApJ, 776, 77 Gerbaldi, M. 1997, A&A, 320, L29 Kral, Q., Matr`a, L., Wyatt, M. C., & Kennedy, G. M. 2017, Currie, T., Burrows, A., Itoh, Y., et al. 2011, ApJ, 729, 128 MNRAS, 469, 521 Currie, T., Rodigas, T. J., Debes, J., et al. 2012, ApJ, 757, 28 Krivov, A. V., & Booth, M. 2018, MNRAS, arXiv:1806.05431 Debes, J. H., Weinberger, A. J., & Kuchner, M. J. 2009, ApJ, Lagrange, A.-M., Gratadour, D., Chauvin, G., et al. 2009, A&A, 702, 318 493, L21 Dent, W. R. F., Wyatt, M. C., Roberge, A., et al. 2014, Science, Lagrange, A.-M., Bonnefoy, M., Chauvin, G., et al. 2010, 343, 1490 Science, 329, 57 Donaldson, J. K., Lebreton, J., Roberge, A., Augereau, J.-C., & Lee, E. J., & Chiang, E. 2016, ApJ, 827, 125 Krivov, A. V. 2013, ApJ, 772, 17 Lyra, W., & Kuchner, M. 2013, Nature, 499, 184 Esposito, T. M., Fitzgerald, M. P., Graham, J. R., & Kalas, P. MacGregor, M. A., Wilner, D. J., Andrews, S. M., Lestrade, 2014, ApJ, 780, 25 Esposito, T. M., Fitzgerald, M. P., Graham, J. R., et al. 2016, J.-F., & Maddison, S. 2015, ApJ, 809, 47 AJ, 152, 85 MacGregor, M. A., Wilner, D. J., Rosenfeld, K. A., et al. 2013, Fern´andez, R., Brandeker, A., & Wu, Y. 2006, ApJ, 643, 509 ApJL, 762, L21 Fitzgerald, M. P., Kalas, P. G., & Graham, J. R. 2007, ApJ, 670, MacGregor, M. A., Matr`a, L., Kalas, P., et al. 2017, ApJ, 842, 8 557 Maness, H. L., Fitzgerald, M. P., Paladini, R., et al. 2008, ApJL, Flaherty, K. M., Hughes, A. M., Andrews, S. M., et al. 2016, 686, L25 ApJ, 818, 97 Maness, H. L., Kalas, P., Peek, K. M. G., et al. 2009, ApJ, 707, Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 1098 2013, PASP, 125, 306 Marino, S., Matr`a, L., Stark, C., et al. 2016, MNRAS, 460, 2933 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, Matr`a, L., Dent, W. R. F., Wyatt, M. C., et al. 2017, MNRAS, ArXiv e-prints, arXiv:1804.09365 464, 1415 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, Mazoyer, J., Boccaletti, A., Augereau, J.-C., et al. 2014, A&A, A&A, 595, A1 569, A29 Gelman, A., & Rubin, D. B. 1992, Statistical Science, 7, 457 Miotello, A., van Dishoeck, E. F., Williams, J. P., et al. 2017, Graham, J. R., Kalas, P. G., & Matthews, B. C. 2007, ApJ, 654, A&A, 599, A113 595 Moerchen, M. M., Telesco, C. M., De Buizer, J. M., Packham, Greaves, J. S., Holland, W. S., Matthews, B. C., et al. 2016, C., & Radomski, J. T. 2007, ApJL, 666, L109 MNRAS, 461, 3910 Mo´or, A., Cur´e, M., K´osp´al, A.,´ et al. 2017, ApJ, 849, 123 16

Morales, F. Y., Rieke, G. H., Werner, M. W., et al. 2011, ApJL, Ricci, L., Carpenter, J. M., Fu, B., et al. 2015, ApJ, 798, 124 730, L29 Rodigas, T. J., Debes, J. H., Hinz, P. M., et al. 2014, ApJ, 783, Mouillet, D., Larwood, J. D., Papaloizou, J. C. B., & Lagrange, 21 A. M. 1997, MNRAS, 292, 896 Schneider, G., Silverstone, M. D., & Hines, D. C. 2005, ApJL, M¨uller, H. S., Schl¨oder, F., Stutzki, J., & Winnewisser, G. 2005, 629, L117 Journal of Molecular Structure, 742, 215 Schneider, G., Grady, C. A., Hines, D. C., et al. 2014, AJ, 148, 59 Mundy, L. G., Looney, L. W., Erickson, W., et al. 1996, ApJL, Steele, A., Hughes, A. M., Carpenter, J., et al. 2016, ApJ, 816, 27 464, L169 Su, K. Y. L., Rieke, G. H., Stapelfeldt, K. R., et al. 2009, ApJ, Mustill, A. J., & Wyatt, M. C. 2009, MNRAS, 399, 1403 705, 314 Olofsson, G., Liseau, R., & Brandeker, A. 2001, ApJL, 563, L77 Thebault, P., Kral, Q., & Augereau, J.-C. 2014, A&A, 561, A16 Olofsson, J., Samland, M., Avenhaus, H., et al. 2016, A&A, 591, Torres, C. A. O., Quast, G. R., da Silva, L., et al. 2006, A&A, A108 460, 695 Ott, J., & Kern, J. 2017, CASA User Reference & Cookbook, Weidenschilling, S. J. 1977, MNRAS, 180, 57 NRAO Williams, J. P., Najita, J., Liu, M. C., et al. 2004, ApJ, 604, 414 Pan, M., Nesvold, E. R., & Kuchner, M. J. 2016, ApJ, 832, 81 Wilner, D. J., Andrews, S. M., & Hughes, A. M. 2011, ApJL, P´astor, P. 2017, A&A, 604, A61 727, L42 Perryman, M. A. C., Lindegren, L., Kovalevsky, J., et al. 1997, Wilner, D. J., MacGregor, M. A., Andrews, S. M., et al. 2018, A&A, 323, L49 ApJ, 855, 56 Pollack, J. B., Hollenbach, D., Beckwith, S., et al. 1994, ApJ, Wyatt, M. C. 2018, ArXiv e-prints, arXiv:1804.08636 421, 615 Redfield, S. 2007, ApJL, 656, L97 Ricarte, A., Moldvai, N., Hughes, A. M., et al. 2013, ApJ, 774, 80 17

Table 1. ALMA Observations

Target Observation #of Projected PrecipitableWater Timeon Date Antennas Baselines(m) Vapor(mm) Target(min) HD32297 2016Jan1 36 15 − 310 1.26 22.2 2016Jun21 36 15 − 704 1.15 44.4 a HD 61005 2013 Dec 4 27 12 − 330 2.99 49.0 a 2014 Mar 20 31 12 − 335 0.79 26.6 2016Jun18 36 15 − 650 0.73 44.4

a Note— Archival observations from ALMA project 2012.1.00437.S (PI: David Rodriguez)

Table 2. Best-fit Model Parameters

Parameter Description HD32297 HD61005 Best-fit Best-fit

Fbelt Planetesimalbeltfluxdensity[mJy] 3.04 ± 0.21 4.82 ± 0.29 Fhalo Halo flux density [mJy] 0.63 ± 0.26 1.54 ± 0.48 Rin Planetesimalbeltinneredge[AU] 78.5 ± 8.1 41.9 ± 0.9 Rhalo Planetesimal belt outer edge and disk halo inner edge [AU] 122 ± 3 67.0 ± 0.5 Rout Halo outer edge [AU] 440 ± 32 188 ± 8 x1 Planetesimalbeltpowerlawgradient 2.23 ± 0.61 2.87 ± 0.19 x2 Halo power law gradient −6.20 ± 0.30 −5.52 ± 0.14 a i Disk inclination [◦] 83.6(+4.6, −0.4) 85.6 ± 0.1 PA Disk position angle [◦] 47.9 ± 0.2 70.3 ± 0.1

a 2 Note— We estimate an additional 5% uncertainty towards higher inclinations from a simple χ analysis in the image plane (see Section 4.1.1 for discussion)