22.07.2013 Views

Annual Report 2011 Max Planck Institute for Astronomy

Annual Report 2011 Max Planck Institute for Astronomy

Annual Report 2011 Max Planck Institute for Astronomy

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

62 III. Selected Research Areas<br />

q<br />

0.9<br />

0.86<br />

0.82<br />

0.9<br />

0.86<br />

0.82<br />

a)<br />

c)<br />

–0.3<br />

–0.2 –0.1 0<br />

Fig. III.3.5: Posterior distribution <strong>for</strong> median projected axis ratio<br />

(shape) versus velocity anisotropy at the end of an MCMC<br />

run <strong>for</strong> Centauri with the points coloured according to their<br />

likelihood (red high, blue low). a): likelihood <strong>for</strong> x proper<br />

motions only, b) likelihood <strong>for</strong> y proper motions only, c) likelihood<br />

<strong>for</strong> line-of-sight velocities only, d) total likelihood. The<br />

<strong>for</strong>eground or background population) are removed via<br />

a series of cuts be<strong>for</strong>e proceeding with the modeling;<br />

this is a tricky process – make the cuts too conservative<br />

and true members will be excised as well as the<br />

contaminants, make the cuts too generous and contaminants<br />

will still remain. A better approach is to include<br />

all stars in a model, both likely members and suspected<br />

contaminants, and allow <strong>for</strong> the presence of a contaminant<br />

population in the likelihood calculations. As our<br />

models have many parameters a simple grid search is<br />

not feasible, so we turn to Markov-Chain Monte Carlo<br />

(MCMC) methods in order to sample the large parameter<br />

space efficiently.<br />

One prime example that highlights what we are currently<br />

able to achieve is the Galactic globular cluster<br />

w Centauri. Located only 5 kpc from the Sun, it is large<br />

and bright and has been observed many times with many<br />

different instruments, over a long time baseline. As a result,<br />

there are line-of-sight velocities and metal abundances<br />

available <strong>for</strong> thousands of stars, and proper motions<br />

measurements available <strong>for</strong> hundreds of thousands<br />

0.1<br />

ln L x<br />

ln L z<br />

b)<br />

d)<br />

ln L y<br />

ln L total<br />

0.2 –0.3<br />

b<br />

–0.2 –0.1 0 0.1 0.2<br />

degeneracy between shape and anisotropy is seen in the scatter<br />

of the points. It is clear from the first three panels that the<br />

individual velocities do not converge on the correct answer;<br />

however their combination (as shown in panel d) does break<br />

the degeneracy and converges nicely.<br />

of stars. w Centauri is an interesting object to study as<br />

it demonstrates many qualities in common with globular<br />

clusters, such as an apparent absence of dark matter,<br />

and also many qualities in common with dwarf spheroidal<br />

galaxies, such as a complex star <strong>for</strong>mation history.<br />

There is also an ongoing debate concerning the presence<br />

(or absence) of an intermediate-mass black hole (IMBH)<br />

at its center. High-quality data and sophisticated modeling<br />

techniques should enable us to finally answer some<br />

of the questions that linger over the nature of this curious<br />

object and how it has <strong>for</strong>med and evolved. The result<br />

of one MCMC model run <strong>for</strong> w Centauri is shown in<br />

Fig. III.3.5 <strong>for</strong> two of our parameters: the median deprojected<br />

axis ratio (q, a proxy <strong>for</strong> shape) and velocity anisotropy<br />

(b). These plots demonstrate the shape-anisotropy<br />

degeneracy that exists, and that proper motions and<br />

line-of-sight velocities must be used together in order to<br />

break the degeneracy.<br />

A second object to which we are applying our modeling<br />

tools is the prototypical core collapsed Galactic<br />

globular cluster M 15. A longstanding debate is wheth-<br />

Credit: Glenn van de Ven

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

Saved successfully!

Ooh no, something went wrong!