Astronomical Distances
Download
1 / 92

Sternberg Astronomical Institute Moscow University - PowerPoint PPT Presentation


  • 105 Views
  • Uploaded on

Astronomical Distances or Measuring the Universe (Chapters 5 & 6) by Rastorguev Alexey, professor of the Moscow State University and Sternberg Astronomical Institute, Russia. Sternberg Astronomical Institute Moscow University. Content.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Sternberg Astronomical Institute Moscow University' - tommy


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

Astronomical Distancesor Measuring the Universe(Chapters 5 & 6)by Rastorguev Alexey,professor of the Moscow State University and Sternberg Astronomical Institute, Russia

Sternberg Astronomical

Institute

Moscow University


Content
Content

  • Chapter Five: Main-Sequence Fitting, or the distance scale of star clusters

  • Chapter Six:Statistical parallaxes


Chapter five
Chapter Five

Main-Sequence Fitting, or

the distance scale of star clusters

  • Open clusters

  • Globular clusters


  • Main idea: to use the advantages of measuring photometric parallax of a whole stellar sample, i.e. close group of stars of common nature: of the same

    • age,

    • chemical composition,

    • interstellar extinction,

      but of different initial masses


Advantages of using star clusters as the standard candles 1
Advantages of using star clusters as the “standard candles” - 1

  • (a) Large statistics (N~100-1000 stars) reduce random errors as ~N-1/2. All derived parameters are more accurate than for single star

  • (b) All stars are of the same age. Star clusters are the only objects that enable direct age estimate, study of the galactic evolution and the star-formation history

  • (c) All stars have nearly the same chemical composition, and the differences in the metallicity between the stars play no role


Advantages of using star clusters as the standard candles 2
Advantages of using star clusters as the “standard candles” - 2

  • (d) Simplify the identification of stellar populations seen on HRD

  • (e) Large statistics also enables reliable extinction measurements

  • (f) Can be distinguished and studied even at large (5-6 kpc, for open clusters) distances from the Sun

  • (g) Enable secondary luminosity calibration of some stars populated star clusters – Cepheids, Novae and other variables


  • DataBase on open clusters: candles” - 2W.Dias, J.Lepine, B.Alessi (Brasilia)

  • Latest Statistics - Version 2.9 (13/apr/2008):

  • Number of clusters: 1776

  • Size: 1774 (99.89%)

  • Distance: 1082 (60.92%)

  • Extinction: 1061 (59.74%)

  • Age: 949 (53.43%)

  • Distance, extinction and age: 936 (52.70%)

  • Proper motion (PM): 890 (50.11%)

  • Radial velocity (VR): 447 (25.17%)

  • Proper motion and radial velocity: 432 (24.32%)

  • Distance, age, PM and VR:379 (21.34%)

  • Chemical composition [Fe/H]:158 ( 8.90%)

  • “These incomplete results point out to the observers that a large effort is still needed to improve the data in the catalog” (W.Dias)


Astrophysical backgrounds of isochrone fitting technique
Astrophysical backgrounds of candles” - 2“isochrone fitting” technique:

  • (a) Distance measurements: photometric parallax, or magnitude difference (m-M)

  • (b) Extinction measurements: color change, or “reddening”

  • (c) Age measurements: different evolution rate for different masses, declared itself by the turn-off point color and luminosity

    -----------------------------------------------

  • Common solution can be found on the basis of stellar evolution theory, i.e. on the evolutional interpretation of the CMD


  • Difference with single-stars method: candles” - 2

  • Instead of luminosity calibrations of single stars, we have to make luminosity calibration of all Main Sequence as a whole: ZAMS (Zero-Age Main Sequence), and isochrones of different ages (loci of stars of different initial masses but of the same age and metallicity)


  • Important note: candles” - 2Theoretical evolutionary tracks and theoretical isochrones are calculated in lg Teff – Mbol variables

  • Prior to compare directly evolution calculations with observations of open clusters, we have to transform Teff to observed colors, (B-V) etc., and bolometric luminosities lg L/LSun and magnitudes Mbol to absolute magnitudes MV etc. in UBV… broad-band photometric system (or others)


  • Important and necessary step: candles” - 2 the empirical (or semi-empirical) calibration of “color-temperature” and “bolometric correction-temperature” relations from data of spectroscopically well-studied stars of

    • (a) different colors

    • (b) different chemical compositions

    • (c) different luminosities

      with accurately measured spectral energy distributions (SED),

      or calibration based on the principles of the “synthetic photometry”


Bolometric magnitudes and bolometric corrections
Bolometric magnitudes and bolometric corrections candles” - 2

  • Bolometric Magnitude, Mbol, specifies total energy output of the star (to some constant):

  • Bolometric Correction,BCV, is defined as the correction to V magnitude:

>1

BCV≤ 0

By definition,Mbol = MV + BCV


Example: candles” - 2BCV vs lg Teff:unique relation for all luminosities

From P.Flower (ApJ

V.469, P.355, 1996)


  • Note: candles” - 2Maximum BCV ~0 at lgTeff~3.8-4.0 (for F3-F5 stars), when maximum of SED coincides with the maximum of V-band sensitivity curve

  • Obviously, the bolometric corrections can be calculated to the absolute magnitude defined in each band


  • For modern color-temperature and BC-temperature calibrations see papers by:

  • P. Flower (ApJ V.469, P.355, 1996):

    lgTeff - BCV – (B-V) from observations

  • T. Lejeune et al. (A&AS V.130, P.65, 1998):

    Multicolor synthetic-photometry approach;

    lgTeff–BCV–(U-B)-(B-V)-(V-I)-(V-K)-…-(K-L),

    for dwarf and giants with [Fe/H]=+1…-3

    (with step 0.5 in [Fe/H])


  • lgT see papers by:eff – (B-V)

  • for different luminosities; based on observations

  • (from P.Flower, ApJ V.469, P.355, 1996)

  • Shifted down by Δ lgTeff = 0.3 relative to next more luminous class for the sake of convenience



  • Before HIPPARCOS mission, astronomers used Hyades “convergent-point” distance as most reliable zero-point of the ZAMS calibration and the base of the distance scale of all open clusters

  • Recently, the situation has changed, but Hyades, along with other ~10 well-studied nearby open clusters, still play important role in the calibration of isochrones via their accurate distances



Pleiades problem: Leeuwen, 2007)

HST gives smaller

parallax (by ~8%)

ΔMHp≈ -0.17m

  • Combined MHp – (V-I) HRD for 8 nearby open clusters constructed by revised HIPPARCOS parallaxes of individual stars (from van Leeuwen, 2007) and corrected for small light extinction

  • Hyades MS shift (red squares) is due to

    • Larger [Fe/H]

    • Larger age ~650 Myr

  • Bottom envelope (----) can be treated as an observed ZAMS

MHp

(V-I)



Primary empirical calibration of the Hyades MS & isochrone for different colors, by HIPPARCOS parallaxes(M.Pinsonneault et al. ApJ V.600, P.946, 2004)

Solid line: theoretical isochrone with

Lejeune et al. (A&AS V.130, P.65, 1998)

color-temperature calibrations

MV


Zams and hyades isochrones sensitivity to the age for 650 100 myr from y lebreton 2001
ZAMS and Hyades isochrones: sensitivity to the age for 650±100 Myr (from Y.Lebreton, 2001)

  • Fitting color of the turn-off point

ZAMS


  • Best library of isochrones recommended to calculate cluster distances, ages and extinctions:

  • L.Girardi et al. “Theoretical isochrones in several photometric systems I. (A&A V.391, P.195, 2002)

  • Theoretical background:

    • (a) Evolution tracks calculations for different initial stellar masses (0.15-7MSun) and metallicities (-2.5…+0.5) (also including α-element enhanced models and overshooting)

    • (b) Synthetic spectra by Kurucz ATLAS9

    • (c) Synthetic photometry (bolometric corrections and color-temperature relations) calibrated by well-studied spectroscopic standards


  • L.Girardi et al. “Theoretical isochrones in several photometric systems I. (A&A V.391, P.195, 2002)

  • Distribution of spectra in Padova library on lg Teff – lg g plane for [Fe/H] from -2.5 to +0.5

  • Wide variety of stellar models, from giants to dwarfs and from hot to cool stars, to compare with observations in a set of popular photometric bands:

  • UBVRIJHK (Johnson-Cousins-Glass), WFPC2 (HST), …

Giants



Metallicity effects on isochrones modelling variables m bol t eff
Metallicity effects on isochrones: the diskmodelling variables, Mbol - Teff

Turn-off point




  • The corrections the diskΔM and ΔCI (CI –Color Index) vs Δ[Fe/H] or ΔZ to isochrones, taken for solar abundance, can be found either

    • from theoretical calculations,

    • or empirically, by comparing multicolor photometric data for clusters with different abundances and with very accurate trigonometric distances


  • Metallicity differences can be taken into account by the disk

    • (a) Adding the corrections to absolute magnitudes ΔM and to colors ΔCI to ZAMS and isochrone of solar composition. These corrections can follow both from observations and theory.

    • (b) Direct fitting of observed CMD by ZAMS and isochrone of the appropriate Z – now most common used technique

  • These methods are completely equivalent


  • Ideally, we should estimate [Fe/H] (or Z) the diskprior to fitting CMD by isochrones

  • If it is not the case, systematic errors in distances (again errors!) may result

  • Open question: differences in Helium content (Y). Theoretically, can play important role. As a rule, evolutionary tracks and isochrones of solar Helium abundance (Y=0.27-0.29) are used


  • L.Girardi the disk et al. (2002) database on isochrones and evolutionary tracks is of great value – it provides us with “ready-to-use” multicolor isochrones for a large variety of the parameters involved (age, [Fe/H], [α/Fe], convection,…)


  • Example: the disk Normalized transmission curves for two realizations of popular UBVRIJHK systems as compared to SED (spectral energy distributions) of some stars (from L.Girardi et al., 2002)

  • See next slides for ZAMS and some isochrones


0.1 the disk

1

  • Theoretical isochrones (color - MV magnitude diagrams) for solar composition (Z=0.019) and cluster ages 0.1 Gyr, 1 Gyr and 10 Gyr (L.Girardi et al., 2002, green solid lines)

10 Gyr


0.1 the disk

  • Theoretical isochrones (NIR color-magnitude diagrams) for solar composition (Z=0.019) and cluster ages 0.1 Gyr, 1 Gyr and 10 Gyr (L.Girardi et al., 2002, green solid lines)

1

What are fancy shapes !

1 Gyr


Girardi et al isochrones in modelling variables m bol lg t eff more detailed age grid
Girardi et al. isochrones in modelling variables the diskMbol – lg Teff (more detailed age grid)

ZAMS


Optics nir
Optics the disk NIR

  • The same but for “standard” multicolor system

ZAMS

ZAMS


How estimate age extinction and the distance 1 st variant
How estimate age, extinction and the distance? the disk1st variant

  • (a) Calculate color-excess CE for cluster stars on two-color diagram like (U-B) – (B-V). Statistically more accurate than for single star. Highly desirable to use a set of two-color diagrams as (U-B) – (B-V)and (B-V) – (V-I) etc., to reduce statistical and systematical errors


How estimate age extinction and the distance 1 st variant1
How estimate age, extinction and the distance? the disk1st variant

  • (b) If necessary, add corrections for [Fe/H] differences to ZAMS and isochrones family, constructed for solar abundance

  • (c) Shift observed CMD horizontally, the offset being equal to the color-excess found at (a) step, and then vertically, by ΔM, to fit proper ZAMS isochrone, i.e. cluster turn-off point. Calculate true distance modulus as (V-MV)0 = ΔV - RV∙E(B-V)

  • (for V–(B-V) CMD)


How estimate age extinction and the distance 2 nd variant
How estimate age, extinction and the distance? the disk2nd variant

  • (a) If necessary, add corrections for [Fe/H] differences to ZAMS and isochrones family, constructed for solar metallicity

  • (b) Match observed cluster CMD (color-magnitude diagram) to ZAMS and isochrone trying to fit cluster turn-off point

  • (c) Calculate horizontal and vertical offsets:

    H: Δ (color) = CE (color excess)

    V: (m-M) = (m-M)0 + R· CE

    (m-M)0 – true distance modulus


How estimate age extinction and the distance 2 nd variant1
How estimate age, extinction and the distance? the disk2nd variant

  • (d) Make the same procedure for all available observations in other photometric bands

  • (e) Compare all (m-M)0 and CE ratios. For MS fitting performed properly,

    • distances will be in general agreement,

    • CE ratios will be in agreement with accepted “standard” extinction law

You can start writing paper !


Ms fitting example pleiades good case
MS-fitting example: Pleiades, good case the disk

Magnitudes offset

gives

ΔV=(V-MV)0+RV∙E(B-V)

(m-M)0 = 5.60

E(B-V)=0.04

lg (age) = 8.00

ZAMS

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


Young distant cluster good case
Young distant cluster, good case the disk

(m-M)0=12.55

E(B-V)=0.38

lg (age)=7.15

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


H per cluster
h Per cluster the disk

(m-M)0=13.65

E(B-V)=0.56

lg (age)=7.15

RSG

(Red

Super-

Giants)

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


RSG the disk

(m-M)0=12.10

E(B-V)=0.32

lg (age)=8.22

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


Older and older
Older and older… the disk

(m-M)0=7.88

E(B-V)=0.02

lg (age)=9.25

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


Very old open cluster m67
Very old open cluster, M67 the disk

(m-M)0=9.60

E(B-V)=0.03

lg (age)=9.60

G.Meynet et al.

(A&AS V.98, P.477,

1993)

Geneva isochrones


Optical data d an et al apj v 671 p 1640 2007 some open clusters populated with cepheid variables
Optical data: D.An et al. (ApJ V.671, P.1640, 2007) the disk(Some open clusters populated with Cepheid variables)




  • Important note: variables (from Open cluster field is often contaminated by large amount of foreground and background stars, nearby as well as more distant non-members

  • Prior to “MS-fitting” it is urgently recommended to “clean” CMD for field stars contribution, say, by selecting stars with similar proper motions on μx - μy vector-point diagram:

    (kinematic selection; reason –

    small velocity dispersion)

Field stars

Cluster stars


Ms fitting accuracy best case multicolor photometry d an et al apj v 655 p 233 2007
MS-fitting accuracy variables (from (best case, multicolor photometry)(D.An et al ApJ V.655, P.233, 2007)

  • Random error of MS-fitting

    • with spectroscopic [Fe/H]: δ(m-M)0 ≈ ±0.02m, i.e ~ 1% in the distance

  • Systematic errors due to uncertainties of calibrations, [Fe/H] and α-elements, field contamination and contribution of unresolved binaries

    • δ(m-M)0 ≈ ±0.04-0.06m, i.e. 2-4% in the distance

  • Uncertainties of Helium abundance may result in even larger systematic errors…


  • For distant clusters, with CMD contaminated by foreground/background stars, and uncertainties in [Fe/H], errors may increase to

    Δ(m-M)0≈±0.1m(random) ± 0.2m(systematic)

    Typical distance accuracy of remote open clusters is ~10-15%


  • Isochrones fitting is equally applicable to globular clusters, but this is not the only method of the distance estimates

  • Good idea to use additional horizontal branch luminosity

    indicators, including

    RR Lyrae variables

    (with nearly constant

    luminosity, see later)

RR Lyrae

BHB

(EHB)

TP


  • D.An clusters, but this is not the only method of the distance estimates et al. (arXiv:0808.0001v1)

  • Isochrones (MS + giant branch) for globular clusters of different [Fe/H] in (u g r i z)photometric bands (SDSS)

more metal-deficient

u 3551Å

g 4686Å

r 6165Å

i 7481Å

z 8931Å

Å


Isochrones fitting clusters, but this is not the only method of the distance estimates

example: M92

Age step 2 Gyr

Theoretical background of

this method is quite

straightforward

Galactic Globular Clusters

are distant objects and

very difficult to study,

even with HST

Reliable photometric data

exist mostly for brightest

stars: Horizontal Branch,

Red Giant Branch and

SubGiants



  • For CMDs of globular clusters, without pronounced Main Sequence, there are other methods of age estimates, based on

    • magnitude difference between Horizontal Branch and Turn-Off Point (“vertical method”)

    • color difference between Turn-Off Point and Giant Granch (“horizontal” method)



Gyr methods of age estimates of globular clusters

“Horizontal”

method

calibrations:

Color offset

vs [Fe/H]

for different

ages

Gyr

Gyr

Gyr


“Vertical” method calibrations: methods of age estimates of globular clustersmagnitude difference vs [Fe/H] for different ages

Gyr



ω because of Cen

Multiple populations ?

He abundance

differences ?

NGC 1806

(LMC)


Chapter six
Chapter Six because of

Statistical parallaxes


Astronomical background
Astronomical background because of

  • Statistical parallaxes provides very powerful tool used to refine luminosity calibrations of secondary “standard candles”, such as RR Lyrae variables, Cepheids, bright stars of constant luminosity, and isochrones applied for main-sequence fitting

  • Statistical parallax technique involves space velocities of uniform sample of objects – at first glance, it sounds as strange and unusual…


Main idea
Main idea because of

  • To match the tangential velocities (VT = k r μ, proportional to distance scale of the sample of studied stars) and radial velocities VR(independent on the distance scale), under three-dimensional normal (ellipsoidal) distribution of the residual velocities

VT=k r μ

VR

Sun

r


If all accepted distances are because of

systematically larger (shorter)

than true distances, then overestimated (underestimated)

tangential velocities will generally distort the ellipsoidal distribution of residual velocities, and the velocity ellipsoids will look like



  • One of the first attempts to calculate statistical parallax of stars has been made by E.Pavlovskaya in the paper entitled “Mean absolute magnitude and the kinematics of RR Lyrae stars” (Variable Stars V.9, P.349, 1953)

  • Her estimate <MV>RR ≈ +0.6m was widely used and kept before early 1980th and even recently, differ only slightly on modern value for metal-deficient RR Lyrae (~0.75m)


  • First rigorous formulation of modern statistical parallax technique have been done by:

  • S.Clube, J.Dawe in “Statistical Parallaxes and the Fundamental Distance Scale-I & II” (MNRAS V.190, P.575; P.591, 1980)

  • C.A.Murray in his book “Vectorial Astrometry” (Bristol: Adam Hilger, 1983)


  • Modern (3D) formulation of the statistical parallax technique enables

    • (a) To refine the accepted distance scale and absolute magnitude calibration used

    • (b) To take into account all observational errors

    • (c) To calculate full set of kinematical parameters of a given uniform stellar sample (space velocity of the Sun, rotation curve or other systemic velocity field, velocity dispersion etc.)

  • Advanced matrix algebra is required, so

    only brief description follows


  • Detailed description of the 3D statistical parallax technique can be found only in A.Rastorguev’s (2002) electronic textbook in

  • http://www.astronet.ru/db/msg/1172553

  • “The application of the maximum-likelyhood technique to the determination of the Milky Way rotation curve and the kinematical parameters and distance scale of the galactic populations”

  • (in russian)


Photometric distances are calculated by star’s apparent and absolute magnitudes. Absolute magnitudes are affected by random and systematic errors. The last can be treated as systematic offset of distance scale used, ΔM.

Statistical parallax technique distinguishes:

expected distancere, calculated by accepted mean absolute magnitude of the sample (after luminosity calibration);

refined distancer, calculated by refined mean absolute magnitude of the sample (after application of statistical parallax technique);

true distancert, appropriate to true absolute magnitude of the star (generally unknown).


Toy distribution and absolute magnitudes. Absolute magnitudes are affected by random and systematic errors. The last can be treated as systematic offset of distance scale used,

of accepted and

refined absolute

magnitudes

ΔMV

Refined

mean

Excpectedandrefined absolute magnitudes (distances) differ due to systematic offset of the absolute magnitude, ΔMV, just what we have to found

Trueandrefined absolute magnitudes differ due to random factors (chemistry, stellar rotation, extinction, age etc.). Random scatter can be described in terms of absolute magnitude standard (rms) variance,σM

Expected

mean

True MV


Kinematic model of the stellar sample
Kinematic model of the stellar sample: and absolute magnitudes. Absolute magnitudes are affected by random and systematic errors. The last can be treated as systematic offset of distance scale used,

Four components of 3D-velocity:

  • Local sample motionrelative the Sun,V0

  • Systematic motion, including differential rotation and noncircular motions, unified by the vectorVSYS

  • Ellipsoidal (3D-Gaussian) distribution of true residual velocities, manifested by star’s random velocity vectorη

  • Errors: in radial velocity and proper motions


L = <ΔV·ΔV T>, T –transposition sign

Vloc (re) is what we measure !


  • V predicted by the kinematical model is expected to have 3D-Gaussian distribution asloc(re) is defined in the local “astrocentric” coordinate system (see picture) via:

  • VR radial velocity, independent on distance re

  • Vl = kre μl velocity on the galactic longitude

  • Vb= kre μb velocity on the galactic latitude

Vb

VR

Vl

re

Galactic disk

Sun


After some advanced algebra

Covariance tensor predicted by the kinematical model is expected to have 3D-Gaussian distribution as

After some advanced algebra:

Observed errors

Ellipsoidal distribution

Systematic motion: (a) relative to the Sun and (b) rotation

where


Individual velocities of all stars are independent on each other; in this case full (N-body) distribution function is the product of N individual functions f,

where N is the number of stars, A is the “vector” of unknown parameters to be found. Maximum Likelihood principle states that observed set of velocity differences is most probable of all possible sets. The set of parameters,A, is calculated under assumption that F reaches its maximum (or minimum, for maximum-likelihood functionLF )


For 3d gaussian distributions functions f lf can be written as a function of a
For 3D-Gaussian distributions functions other; in this case full (N-body) distribution function is the f, LF can be written as a function of A

Here|L| is matrix determinant, L-1 is inverse matrix. By minimizing LF by A, we calculate all important parameters {A}, for example:


Robust statistical parallax method applied to local disk populations
Robust statistical parallax method: other; in this case full (N-body) distribution function is the applied to local disk populations

Astronomical background:

A, Oort constant, derived from proper motions alone, depends on the distance scale used, whereas A, derived only from radial velocities, do not depend on the distances

As a result, scale factor can be estimated by requirement that both A values are equal to each other


Local oort s approximation
Local Oort’s approximation other; in this case full (N-body) distribution function is the

Differential rotation contribution to space

velocity components in local approximation

r << R0(or|RP–R0 | << R0):

To first order by the ratio r/R0

in the expansion for the angular velocity:


Differential rotation effect to radial velocity vr
Differential rotation effect to radial velocity other; in this case full (N-body) distribution function is the Vr :

From 1st Bottlinger equation (for radial velocity)

calculate contribution of the differential rotation toVr :

Linearity onr,

“double wave” onl

AOort’s constant (definition)


Differential rotation effect to tangential velocity v l
Differential rotation effect to tangential velocity other; in this case full (N-body) distribution function is the Vl :

From 2nd Bottlinger equation (for velocity on l)

calculate contribution of the differential rotation toVl :

Linearity onr,

“double wave” onl


Oort constant a and the refinement of the distance scale
Oort constant other; in this case full (N-body) distribution function is the A and the refinement of the distance scale

A0Vrdepends on the distance scale: A0Vr~ p -1

(decreases with increasing distances)

A0μldo not depend on the distance scale,A0μl≈const

The requirementA0Vr(p) ≈A0μl– robust method of the

adjusting the scale factor p


Illustration of the robust technique
Illustration of the robust technique other; in this case full (N-body) distribution function is the

AVr

Optimal value of

the scale factor

Aμl


F a q how the corrections to absolute magnitudes are affected by the
F.A.Q. How the corrections to absolute magnitudes are affected by the:

  • (a) Shape of the velocity distribution (deviation from expected 3D-Gaussian form)

  • (b) Vertex deviation of the velocity ellipsoid (velocity-position correlations)

  • (c) Misestimates of the observation errors

  • (d) Non-uniform space distribution of stars

  • (e) Sample size

  • (f) Malmquist bias (excess of intrinsically bright stars in the magnitude-limited stellar sample)

  • (g) Interstellar extinction

  • (h) Misidentification of stellar populations


  • Possible factors of affected by the:systematic offsets have been analyzed by P.Popowski & A.Gould in the papers “Systematics of RR Lyrae statistical parallax. I-III” (ApJ V.506, P.259, P.271, 1998; ApJ V.508, P.844, 1998) (a) analytically and (b) by Monte-Carlo simulations, and applied to the sample of RR Lyrae variables


P popowski a gould 1998
P.Popowski & A.Gould (1998): affected by the:

  • “Statistical parallax method … is extremely robust and insensitive to several different categories of systematic effects”

  • “… statistical errors are dominated by the size of the stellar sample”

  • … sensitive to systematic errors in the observed data

  • … Malmquist bias should be taken into account prior to calculations


  • To eliminate the effects due to non-uniformity of the sample, bimodal versions of the statistical parallax method can be used (A.Rastorguev, A.Dambis & M.Zabolotskikh “The Three-Dimensional Universe with GAIA”,ESA SP-576, P.707, 2005)

  • Example: RR Lyrae sample of halo and thick disk stars


  • Statistical parallax technique is considered as sample, the absolute method of the distance scale calibration, though it exploits prior information on the adequate kinematic model of the sample studied

  • After HIPPARCOS, luminosities and distance scales of RR Lyrae stars, Cepheids and young open clusters have been analyzed by the statistical parallax technique


ad