Protoplanetary disks are a byproduct of the star formation process. In the dense mid-plane of these disks, planetesimals and planets are expected to form. The first step in planet formation is the growth of dust particles from submicrometer-sized grains to macroscopic mm-sized aggregates. The grain growth is accompanied by radial drift and vertical segregation of the particles within the disk. To understand this essential evolutionary step, spatially resolved multi-wavelength observations as well as photometric data are necessary which reflect the properties of both disk and dust. We present the first spatially resolved image obtained with NACO at the VLT in the L$_\text{p}$ band of the near edge-on protoplanetary disk FS Tau B. Based on this new image, a previously published Hubble image in H band and the spectral energy distribution from optical to millimeter wavelengths, we derive constraints on the spatial dust distribution and the progress of grain growth. For this purpose we perform a disk modeling using the radiative transfer code MC3D. Radial drift and vertical sedimentation of the dust are not considered. We find a best-fit model which features a disk extending from $2\,\text{AU}$ to several hundreds AU with a moderately decreasing surface density and $M_\text{disk}=2.8\,\times\,10^{-2}\,\text{M}_\odot$. The inclination amounts to $i=80^\circ$. Our findings indicate that substantial dust grain growth has taken place and that grains of a size equal to or larger than $1\,\text{mm}$ are present in the disk. In conclusion, the parameters describing the vertical density distribution are better constrained than those describing the radial disk structure.

We present a detailed comparison of two approaches, the use of a pre-calculated database and simulated annealing (SA), for fitting the continuum spectral energy distribution (SED) of astrophysical objects whose appearance is dominated by surrounding dust. While pre-calculated databases are commonly used to model SED data, only few studies to date employed SA due to its unclear accuracy and convergence time for this specific problem. From a methodological point of view, different approaches lead to different fitting quality, demand on computational resources and calculation time. We compare the fitting quality and computational costs of these two approaches for the task of SED fitting to provide a guide to the practitioner to find a compromise between desired accuracy and available resources. To reduce uncertainties inherent to real datasets, we introduce a reference model resembling a typical circumstellar system with 10 free parameters. We derive the SED of the reference model with our code MC3D at 78 logarithmically distributed wavelengths in the range [0.3um, 1.3mm] and use this setup to simulate SEDs for the database and SA. Our result shows directly the applicability of SA in the field of SED modeling, since the algorithm regularly finds better solutions to the optimization problem than a pre-calculated database. As both methods have advantages and shortcomings, a hybrid approach is preferable. While the database provides an approximate fit and overall probability distributions for all parameters deduced using Bayesian analysis, SA can be used to improve upon the results returned by the model grid.

Circumstellar disks are characteristic for star formation and vanish during the first few Myr of stellar evolution. During this time planets are believed to form in the dense midplane by growth, sedimentation and aggregation of dust. Indicators of disk evolution, such as holes and gaps, can be traced in the spectral energy distribution (SED) and spatially resolved images. We aim to construct a self-consistent model of HH 30 by fitting all available continuum observations simultaneously. New data sets not available in previous studies, such as high-resolution interferometric imaging with the Plateau de Bure Interferometer (PdBI) at lambda = 1.3 mm and SED measured with IRS on the Spitzer Space Telescope in the mid-infrared, put strong constraints on predictions and are likely to provide new insights into the evolutionary state of this object. A parameter study based on simulated annealing was performed to find unbiased best-fit models for independent observations made in the wavelength domain lambda ~ 1 micron ... 4 mm. The method essentially creates a Markov chain through parameter space by comparing predictions generated by our self-consistent continuum radiation transfer code MC3D with observations. We present models of the edge-on circumstellar disk of HH 30 based on observations from the near-infrared to mm-wavelengths that suggest the presence of an inner depletion zone with about 45 AU radius and a steep decline of mm opacity beyond 140 AU. Our modeling indicates that several modes of dust evolution such as growth, settling, and radial migration are taking place in this object. High-resolution observations of HH 30 at different wavelengths with next-generation observatories such as ALMA and JWST will enable the modeling of inhomogeneous dust properties and significantly expand our understanding of circumstellar disk evolution.