10.02.2013 Views

Max Planck Institute for Astronomy - Annual Report 2005

Max Planck Institute for Astronomy - Annual Report 2005

Max Planck Institute for Astronomy - Annual Report 2005

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

64 III. Scientific Work<br />

z [AU]<br />

1<br />

0<br />

−1<br />

800<br />

600<br />

1000 1�10 –14<br />

800<br />

400<br />

0 0.5 1.0<br />

r [AU]<br />

1.5 2.0<br />

Fig. III.2.1: A dust torus is illuminated by a central source. The<br />

dotted lines indicate the density distribution of the dust. The<br />

colors illustrate the temperature distribution (bright � hot) determined<br />

with the »diffusion plus ray-tracing« algorithm. The<br />

corresponding lines of constant temperature are drawn as solid<br />

lines. For comparison, the dashed line shows the result from the<br />

more accurate continuum radiative transfer. The local error of<br />

the sound velocity in this test was always below 5 %.<br />

A Bit Closer to the Ideal Case: Inverse 3D Radiative<br />

Transfer<br />

When images of an object are available it would be<br />

desirable <strong>for</strong> the interpretation of the observational results<br />

to directly calculate a density and temperature distribution<br />

from the telescope images. However, the calculation of<br />

the radiation field <strong>for</strong> a given distribution is so complex<br />

that an inversion is possible only <strong>for</strong> very simply structured<br />

objects. For modeling telescope images one there<strong>for</strong>e<br />

has to rely on multiple iterations, varying the density<br />

and temperature distribution. A glance at the spatially<br />

600<br />

600<br />

resolved dust distributions around young stars suffices to<br />

realize that filamentary or ellipsoidal distributions have to<br />

be assumed in almost all cases. And <strong>for</strong> these, modeling<br />

is virtually out of question because it takes too long to<br />

calculate a configuration.<br />

Scientists at MPIA, together with colleagues from<br />

Bordeaux and Potsdam, succeeded in calculating the first<br />

3D-model of a condensation in the molecular cloud near<br />

rho Ophiuchi. This starless core is a good example of the<br />

initial stage of star <strong>for</strong>mation where a protostar is <strong>for</strong>ming<br />

by gravitational collapse of the core. For the core IsocaM<br />

maps at 7 and 15 µm wavelength were available, where<br />

it is visible in absorption against the infrared radiation of<br />

the background region as a complex structure divided in<br />

two main clumps. In a map obtained with the IraM-30-mtelescope<br />

in the long-wavelength mm-region, however,<br />

the core was detectable by emission of cold dust. The<br />

complex geometry of the distribution precluded a simple<br />

modeling with a 1D- or 2D-model. So 30 clumps with<br />

a three-dimensional Gaussian density distribution were<br />

used that were variable in their position, their axial ratio<br />

and their overall density. With the help of the optimizing<br />

algorithm »simulated annealing« the clumps were shifted<br />

and de<strong>for</strong>med in such a way that the measured absorption<br />

in the mid-infrared range at a given background radiation<br />

was exactly reproduced. However, the distribution of the<br />

clumps along the line of sight and their shapes could not<br />

be obtained by this alone. But since the dust emission in<br />

the mm-range is sensitive to the particular shape of the<br />

core it could be used to infer the missing 3D-in<strong>for</strong>mation.<br />

For the calculation of the mm-emission the scientists<br />

made use of the fact that the temperature witin the core<br />

mainly depends on the mean optical depth of the radiation<br />

entering from the outside. This dependence varied with<br />

the translation and distortion of the clumps by less than 1<br />

K in temperature and could be calculated be<strong>for</strong>e the actual<br />

modeling within the so-called T-tau method by using a<br />

3D-program <strong>for</strong> the transfer of continuum radiation.<br />

The density distribution determined this way shows<br />

many local maxima, with a compact southern absolute<br />

maximum. Comparison with hydrodynamic simulations,<br />

which also allow <strong>for</strong> turbulent motions, shows that such<br />

maxima can be of transient nature and may disappear on<br />

the dynamical timescale of the turbulence. But because<br />

of the high density and compactness of the southern<br />

maximum it can be expected not to dissolve and even to<br />

collapse.<br />

Complex structures as often found in star <strong>for</strong>mation<br />

regions have many free parameters which have to be<br />

determined by comparison with observations. The new<br />

inverse 3D-radiative-transport-method, however, is based<br />

on a simultaneous modeling of several maps containing<br />

several thousand pixels and is there<strong>for</strong>e well-defined. The<br />

development of three-dimensional models is also a natural<br />

consequence of the fact that many observations show<br />

a clearly resolved three-dimensional structure where conventional<br />

one-and two-dimensional models fail.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!