Results of parametric modeling: (a) image of the dust continuum emission of TW Hya’s protoplanetary disk at a wavelength of 0.65 mm. The ellipse at the bottom left corner shows the beam (dimensions : 0.784″ × 0.612″; PA = 91.04°) used to convolve the image. The bar at the upper right corner shows the spatial scale which is equivalent to 10 au. (b) Unconvolved best-fit model of the 0.65 mm dust continuum emission fit with 3 gaussians using MCMC and galario (see Sect. 3.2.2). The substructures (D86 and B96) obtained from the intensity profile (see Fig. 2) are shown with dotted and solid curves respectively. (c) Best-fit model that has been convolved with the same beam as the image. (d) The residuals of the emission as multiples of the rms (σ). Residuals at the level of 3σ are shown with contours. It was obtained when the model was directly subtracted from the data. The beam used to generate this has the same dimensions as that used for the image.