Method of ignoring the systematic variations of stellar parallaxes over the celestial sphere in a kinematic analysis of stellar proper motions

A method for determining the velocity field parameters free from the distortions due to the systematic variations of stellar parallaxes over the celestial sphere is proposed. The method is based on the approximation of parallaxes as a function of coordinates on the sphere using spherical harmonics and can be applied in those cases where the solar motion cannot be eliminated from the stellar proper motions. Numerical experiments have shown that our method is able to obtain accurate coordinates of the solar apex and to calculate the kinematic parameters of the Ogorodnikov-Milne model to within three coefficients of the decomposition of parallaxes into first-order spherical harmonics. Examples of applying the method to the stellar proper motions of the Hipparcos catalogue, which admits checking the results using trigonometric parallaxes, are provided. Such a check has been found to yield a positive result only for nearby stars at heliocentric distances that do not exceed 400 pc and for which the parallaxes were determined with a relative error of at least 30%. An interesting feature of this method is the possibility to construct the shape of the figure which is formed by the deviations of the parallaxes from the sphere corresponding to the average parallaxes of the stars under consideration. It should be specially emphasized that all of this is done in the complete absence of information about the stellar parallaxes. The “solar terms” of the stellar proper motions that are formed by the products of the parallaxes by the solar motion components relative to the centroid of stars are the main source of information about the parallaxes here.     

Kinematic analysis of stellar radial velocities by the spherical harmonics

A method for a kinematic analysis of stellar radial velocities using spherical harmonics is proposed. This approach does not depend on the specific kinematic model and allows both low-frequency and high-frequency kinematic radial velocity components to be analyzed. The possible systematic variations of distances with coordinates on the celestial sphere that, in turn, are modeled by a linear combination of spherical harmonics are taken into account. Theoretical relations showing how the coefficients of the decomposition of distances affect the coefficients of the decomposition of the radial velocities themselves have been derived. It is shown that the larger the mean distance to the sample of stars being analyzed, the greater the shift in the solar apex coordinates, while the shifts in the Oort parameter A are determined mainly by the ratio of the second zonal harmonic coefficient to the mean distance to the stars, i.e., by the degree of flattening of the spatial distribution of stars toward the Galactic plane. The distances to the stars for which radial velocity estimates are available in the CRVAD-2 catalog have been decomposed into spherical harmonics, and the existing variations of distances with coordinates are shown to exert no noticeable influence on both the solar motion components and the estimates of the Oort parameter A, because the stars from this catalog are comparatively close to the Sun (no farther than 500 pc). In addition, a kinematic component that has no explanation in terms of the three-dimensional Ogorodnikov-Milne model is shown to be detected in the stellar radial velocities, as in the case of stellar proper motions.     

Application of vector spherical harmonics for kinematic analysis of stars from zonal catalogues

We solve the problem on a kinematic analysis of the three-dimensional velocity field of stars from zonal catalogues, i.e., catalogues in which the stars are presented at all right ascensions in some declination zones. We have constructed a system of vector spherical harmonics with the properties of completeness and orthogonality for a chosen declination zone. We suggest a method that allows the Ogorodnikov-Milne model parameters in the Galactic coordinate system to be estimated by analyzing the proper motions and radial velocities of stars in the equatorial coordinate system. The vector spherical harmonics are shown to have the following advantages over the standard approach based on a direct leastsquares estimation of the parameters for a specific model. First, in contrast to the standard approach, the new method can reveal all systematic components of the velocity field irrespective of a particular model. Second, it allows one to get rid of the correlation between the sought-for parameters, which presents a serious problem for the conventional method in the case of zonal catalogues. Third, the method of vector spherical harmonics allows the kinematic parameters to be estimated at least by two techniques. Comparison of these two solutions makes it possible to test the standard kinematic model for compatibility with the observational data. The developed method has been tested on the basis of numerical experiments and applied for a kinematic analysis of the proper motions of Tycho-2 stars in the southern hemisphere for which the parallaxes can be estimated using data from the Tycho-2 Spectral Type Catalogue.     

The red giant branch in the Tycho-2 catalogue

Based on multicolor photometry from the 2MASS and Tycho-2 catalogues, we have produced a sample of 38 368 branch red giants that has less than 1% of admixtures and is complete within 500 pc of the Sun. The sample includes 30 671 K giants, 7544Mgiants, 49 C giants, and 104 suspected supergiants or S stars. The photometric distances have been calculated for K, M, and C stars with an accuracy of 40%. Tycho-2 proper motions and PCRV radial velocities are used to analyze the stellar kinematics. The decrease in the stellar distribution density with distance from the Galactic equator approximated by the barometric law, contrary to the Besanconmodel of the Galaxy, and the kinematic parameters calculated using the Ogorodnikov-Milne model characterize the overwhelming majority of the selected K and M giants as disk stars with ages of more than 3 Gyr. A small number of K and M giants are extremely young or, conversely, thick-disk ones. The latter show a nonuniform distribution in the phase space of coordinates and velocities, arguing against isothermality and full relaxation of the disk and for the theory of dynamical streams or superclusters. The spatial distribution and kinematics of the selected C stars force us to consider them as asymptotic branch giants with masses of more than 2M _⊙ and ages of less than 2 Gyr probably associated with the Gould Belt. The offset of the Sun above the Galactic equator has been found from the distribution of stars to be 13 ± 2 pc, which coincides with the previously obtained value for the clump red giants.     

Predicted star counts and mean parallaxes down to the 23rd magnitude

Star counts and mean parallaxes as a function of B, V, R magnitudes down to 23 are presented. The data were computed by the use of two fundamental equations of stellar statistics. The assumed model considers the Galaxy as a symmetrical system with respect to its rotation axis and to its equatorial plane and as composed of the thin disk (main sequence and red giants), the thick disk and spheroid populations. Numbers of stars and mean parallaxes were derived in bins of galactic longitude and latitude of 30° and 10°, respectively. For the computation of the mean parallaxes depending on Galactic coordinates and magnitudes, series of products of Hermite and Legendre polynomials and of Fourier terms were used. The results of this paper may help in the planning of future survey missions and in the design of new telescopes. In addition, mean parallaxes can be used to derive corrections to absolute parallaxes and proper motions for any position in the sky.     

Kinematics of Stars from the TGAS (Gaia DR1) Catalogue

Based on the stellar proper motions of the TGAS (Gaia DR1) catalogue, we have analyzed the velocity field of main-sequence stars and red giants from the TGAS catalogue with heliocentric distances up to 1.5 kpc. We have obtained four variants of kinematic parameters corresponding to different methods of calculating the distances from the parallaxes of stars measured with large relative errors. We have established that within the Ogorodnikov–Milne model changing the variant of distances affects significantly only the solar velocity components relative to the chosen centroid of stars, provided that the solution is obtained in narrow ranges of distances (0.1 kpc). The estimates of all the remaining kinematic parameters change little. This allows the Oort coefficients and related Galactic rotation parameters as well as all the remaining Ogorodnikov–Milne model parameters (except for the solar terms) to be reliably estimated irrespective of the parallax measurement accuracy. The main results obtained from main-sequence stars in the range of distances from 0.1 to 1.5 kpc are: A = 16.29 ± 0.06 km s^?1 kpc^?1, B = ?11.90 ± 0.05 km s^?1 kpc^?1, C = ?2.99 ± 0.06 km s^?1 kpc^?1, K = ?4.04 ± 0.16 km s^?1 kpc^?1, and the Galactic rotation period P = 217.41 ± 0.60 Myr. The analogous results obtained from red giants in the range from 0.2 to 1.6 kpc are: the Oort constants A = 13.32 ± 0.09 km s^?1 kpc^?1, B = ?12.71 ± 0.06 km s^?1 kpc^?1, C = ?2.04 ± 0.08 km s^?1 kpc^?1, K = ?2.72 ± 0.19 km s^?1 kpc^?1, and the Galactic rotation period P = 236.03 ± 0.98 Myr. The Galactic rotation velocity gradient along the radius vector (the slope of the Galactic rotation curve) is ?4.32 ± 0.08 km s^?1 kpc^?1 for main-sequence stars and ?0.61 ± 0.11 km s^?1 kpc^?1 for red giants. This suggests that the Galactic rotation velocity determined from main-sequence stars decreases with increasing distance from the Galactic center faster than it does for red giants.     

Analysis of the Systematic Differences between the Stellar Parallaxes of the TGAS and Hipparcos Catalogues Using Spherical Harmonics

The systematic differences between the trigonometric parallaxes of the Hipparcos and TGAS catalogues have been investigated using spherical harmonics. The most significant harmonics in the expansion have been determined. The distribution of the parallax difference dispersion in various regions of the celestial sphere has also been studied. The distribution of the rms deviation has the simplest form in the ecliptic coordinate system.     

Information from the Kinematics of F and G Stars in the Solar Neighborhood

We have calculated the orbital parameters for 90 stars in Chen et al. and updated the kinematic data for stars in Edvardsson et al. by using the accurate Hipparcos parallaxes and proper motions, and recalculated the \\\\\\\\\\\\-element abundances in Edvardsson et al. in a way consistent with Chen et al. The two sets of data are combined in a study of stellar populations and characteristics of F & G stars in the solar neighborhood. We confirm the result of Chen et al. that a distinguishable group of stars may belong to the thick disk rather than the thin disk. The ages for the stars are determined using the theoretical isochrones of VandenBerg et al. The age-metallicity relation is investigated for different subgroups according to distance from the sun and galactic orbital parameters. It is found that a mixing of stars with different orbital parameters significantly affect the age-metallicity relation for the disk. Stars with orbits confined to the solar circle all have metallicities [Fe/H] > -0.3 irresp     

Investigation of the Galactic bar based on photometry and stellar proper motions

A new method for selecting stars in the Galactic bar based on 2MASS infrared photometry in combination with stellar proper motions from the Kharkiv XPM catalogue has been implemented. In accordance with this method, red clump and red giant branch stars are preselected on the color-magnitude diagram and their photometric distances are calculated. Since the stellar proper motions are indicators of a larger velocity dispersion toward the bar and the spiral arms compared to the stars with circular orbits, applying the constraints on the proper motions of the preselected stars that take into account the Galactic rotation has allowed the background stars to be eliminated. Based on a joint analysis of the velocities of the selected stars and their distribution on the Galactic plane, we have confidently identified the segment of the Galactic bar nearest to the Sun with an orientation of 20°–25° with respect to the Galactic center-Sun direction and a semimajor axis of no more than 3 kpc.     

Late-type members of young stellar kinematic groups – I. Single stars

This is the first paper of a series aimed at studying the properties of late-type members of young stellar kinematic groups. We concentrate our study on classical young moving groups such as the Local Association (Pleiades moving group,     , IC 2391 supercluster (35 Myr), Ursa Major group (Sirius supercluster, 300 Myr), and Hyades supercluster (600 Myr), as well as on recently identified groups such as the Castor moving group (200 Myr). In this paper we compile a preliminary list of single late-type possible members of some of these young stellar kinematic groups. Stars are selected from previously established members of stellar kinematic groups based on photometric and kinematic properties as well as from candidates based on other criteria such as their level of chromospheric activity, rotation rate and lithium abundance. Precise measurements of proper motions and parallaxes taken from the Hipparcos Catalogue, as well as from the Tycho-2 Catalogue, and published radial velocity measurements are used to calculate the Galactic space motions ( U , V , W ) and to apply Eggen's kinematic criteria in order to determine the membership of the selected stars to the different groups. Additional criteria using age-dating methods for late-type stars will be applied in forthcoming papers of this series. A further study of the list of stars compiled here could lead to a better understanding of the chromospheric activity and their age evolution, as well as of the star formation history in the solar neighbourhood. In addition, these stars are also potential search targets for direct imaging detection of substellar companions.     

Trigonometric parallaxes of 29 stars with large proper motions

While implementing the first stage of the Pulkovo program of research on stars with large proper motions, we determined the trigonometric parallaxes of 29 stars (12 ^m < V < 16 ^m ) based on CCD observations with a 26-inch refractor. The mean standard error was 3.7 mas. Comparison of the Pulkovo parallaxes with those obtained at the Observatory of the Yale University and the US Naval Observatory (USNO) has shown that the parallax differences (Pulkovo-Yale/USNO) lie within the limits of their measurement errors in an overwhelming majority of cases. On average, they are −0.6 ± 1.0 mas. No systematic dependences on stellar distance, magnitude, and color in this set of differences have been found. Our comparisons show that the observing and data reduction techniques used in the Pulkovo program of research on fast stars allow highly accurate trigonometric parallaxes of these objects to be obtained. All program stars are within 50 pc of the Sun; most of them belong to the immediate solar neighborhood (D < 25 pc). For two stars (J0522+3814 and J1202+3636), the trigonometric parallaxes have been determined for the first time.     

Astrophysical supplements to the ASCC‐2.5. II. Membership probabilities in 520 Galactic open cluster sky areas

We present a catalogue (CSOCA) of stars residing in 520 Galactic open cluster sky areas which is the result of the kinematic (proper motion) and photometric member selection of stars listed in the homogeneous All‐sky Compiled Catalogue of 2.5Million Stars (ASCC‐2.5).We describe the structure and contents of the catalogue, the selection procedure applied, and the proper motion and photometric membership constraints adopted. In every cluster area the CSOCA contains the complete list of the ASCC‐2.5 stars regardless of their membership probability. For every star the CSOCA includes accurate J2000 equatorial coordinates, proper motions in the Hipparcos system, BV photometric data in the Johnson system, proper motion and photometric membership probabilities, as well as angular distances from the cluster centers for about 166 000 ASCC‐2.5 stars. If available, trigonometric parallaxes, spectral types, multiplicity and variability flags from the ASCC‐2.5, and radial velocities with their errors from the Catalogue of Radial Velocities of Galactic Stars with high precision Astrometric Data (CRVAD) are also given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of the three-dimensional stellar velocity field using vector spherical functions

We apply vector spherical functions to problems of stellar kinematics. Using these functions allows all of the systematic components in the stellar velocity field to be revealed without being attached to a specific physical model. Comparison of the theoretical decomposition coefficients of the equations for a particular kinematical model with observational data can provide precise information about whether the model is compatible with the observations and can reveal systematic components that are not described by this model. The formalism of vector spherical functions is particularly well suited for analyzing the present and future (e.g., GAIA) catalogs containing all three velocity vector components: the propermotions in both coordinates and the radial velocity. We show that there are systematic components in the proper motions of Hipparcos stars that cannot be interpreted in terms of the linear Ogorodnikov-Milne model. The same result is also confirmed by an analysis of the radial velocities for these stars.     

Kinematics and Dynamics of the Galactic Halo from RR Lyrae Variable Stars

Based on a sample of RR Lyrae variable stars including more than 9000 objects with proper motions and distances, we have investigated the kinematics of the Galactic halo from the two-dimensional velocity field. We have used both the proper motions deduced independently by us from the positional data taken from all-sky catalogues in a time interval up to 65 years and the proper motions taken from the Gaia DR2 catalogue. In addition, we have also studied the halo kinematics from the three-dimensional velocity field of ~850 RR Lyrae variables with distances, proper motions, and line-of-sight velocities. The kinematic parameters describing the velocity field have been estimated by the maximum-likelihood method; their change with Galactocentric distance has been investigated. The radial velocity dispersion in spherical coordinates σ_r ≈ 160?170 km s^?1 exceeds its values from previous papers approximately by 20 km s^?1, while the anisotropy parameter β ≈ 0.68?0.72 agrees satisfactorily with previous studies. When estimating the rotation velocity of the population of RR Lyrae stars, we identified the inner and outer halos with weak prograde and retrograde rotations, respectively.     

Astrophysical supplements to the ASCC‐2.5. I. Radial velocity data

We present a catalogue of radial velocities of Galactic stars with high precision astrometric data CRVAD which is the result of the cross‐identification of star lists from the General Catalog of Average Radial Velocities (GCRV) and from the homogeneous All‐sky Compiled Catalogue of 2.5Million Stars (ASCC‐2.5). The CRVAD includes accurate J2000 equatorial coordinates, proper motions and trigonometric parallaxes in the Hipparcos system, Johnson's BV photometric data, spectral types, multiplicity and variability flags from the ASCC‐2.5, and radial velocities, stellar magnitudes and spectral types from the GCRV for 34553 ASCC‐2.5 stars. The CRVAD was used for the construction of a sample of standard stars with accurate astrometric, photometric and radial velocity data for the RAVE project. A second application of the CRVAD , the radial velocity determination for 292 open clusters (including 97 with previously unknown radial velocities), using their newly defined members from proper motions and photometry in the ASCC‐2.5, is briefly described. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Investigation of the motions of fast stars based on observations with the Pulkovo normal astrograph

Astrometric CCD observations of 1123 stars with large proper motions (μ > 300 mas yr⁻¹) from the LSPM (I/298) catalog in the declination zone +30°–+70° have been carried out with the Pulkovo normal astrograph since 2006. The observational program includes mostly stars that previously have not entered into high-accuracy projects to determine the proper motions. Our studies are aimed at determining new proper motions of fast stars in the HCRF/UCAC3 system and searching for stars with invisible companions in the immediate Galactic neighborhoods of the Sun. Having analyzed about 10 000 CCD frames, we have obtained the equatorial coordinates of 414 program stars in the HCRF/UCAC3 system at an accuracy level of 10–50 mas and determined their new proper motions. To derive the proper motions, we have used the data from several star catalogs and surveys (M2000, CMC14, 2MASS, SDSS) as early epochs. The epoch differences range from 5 to 13 years (on average, about 10 years); the mean accuracy of the derived proper motions is 4–5 mas yr⁻¹. For 70 stars, we have revealed significant differences between the derived proper motions and those from the LSPM and I/306A catalogs (these proper motions characterize the mean motion of the photocenter in 50 years or more). Apart from systematic errors, these differences can result from the existence of invisible components of the program stars.     

On the nature of the purported common proper motion companions to the exoplanet host star 51 Peg

Greaves (2006) proposed that three red, high proper motion stars within 10° of 51 Peg (NLTT 54007, 54064, and 55547) are co‐moving companions to this famous exoplanet host star. While the stars clearly have proper motions similar to 51 Peg, the inferred kinematic parallaxes for these stars produce extremely inconsistent color‐magnitude positions 2 to 4 magnitudes below the main sequence. All three stars are likely to be background stars unrelated to 51 Peg (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Galactic orbits of Hipparcos stars: Classification of stars

The Galactic orbits of 27 440 stars of all classes with accurate coordinates and parallaxes of more than 3 mas from the Hipparcos catalogue, proper motions from the Tycho-2 catalogue, and radial velocities from the Pulkovo Compilation of Radial Velocities (PCRV) are analyzed. The sample obtained is much more representative than the Geneva-Copenhagen survey and other studies of Galactic orbits in the solar neighborhood. An estimation of the influence of systematic errors in the velocities on orbital parameters shows that the errors of the proper motions due to the duplicity of stars are tangible only in the statistics of orbital parameters for very small samples, while the errors of the radial velocities are noticeable in the statistics of orbital parameters for halo stars. Therefore, previous studies of halo orbits may be erroneous. The distribution of stars in selection-free regions of the multidimensional space of orbital parameters, dereddened colors, and absolute magnitudes is considered. Owing to the large number of stars and the high accuracy of PCRV radial velocities, nonuniformities of this distribution (apart from the well-known dynamical streams) have been found. Stars with their peri- and apogalacticons in the disk, perigalacticons in the bulge and apogalacticons in the disk, perigalacticons in the bulge and apogalacticons in the halo, and perigalacticons in the disk and apogalacticons in the halo have been identified. Thus, the bulge and the halo are inhomogeneous structures, each consisting of at least two populations. The radius of the bulge has been determined: 2 kpc.     

Programme of the Galactic meridian: Star counts and Galactic models 2. Distribution of absolute proper motions

A comparison of the observed distribution of absolute proper motions with a kinematical model of the Galaxy is presented. Proper motions with respect to galaxies were obtained for about 40 000 stars along the main Galactic meridian and in two fields near the North Galactic pole (programme MEGA). The Galaxy is considered as composed of the disk (main sequence and disk red giants), the thick disk and spheroid populations. For each subsystem, spatial velocity components and their dispersions were computed. The distribution of kinematical parameters were modelled for stars located in different directions of the Galaxy.     

Testing the Distance Scale of the Gaia TGAS Catalogue by the Kinematic Method

We have studied the simultaneous and separate solutions of the basic kinematic equations obtained using the stellar velocities calculated on the basis of data from the Gaia TGAS and RAVE5 catalogues. By comparing the values of Ω'₀ found by separately analyzing only the line-of-sight velocities of stars and only their proper motions, we have determined the distance scale correction factor p to be close to unity, 0.97 ± 0.04. Based on the proper motions of stars from the Gaia TGAS catalogue with relative trigonometric parallax errors less than 10% (they are at a mean distance of 226 pc), we have found the components of the group velocity vector for the sample stars relative to the Sun (U, V,W)_⊙ = (9.28, 20.35, 7.36) ± (0.05, 0.07, 0.05) km s^?1, the angular velocity of Galactic rotation Ω⁰ = 27.24 ± 0.30 km s?1 kpc?1, and its first derivative Ω'₀ = ?3.77 ± 0.06 km s^?1 kpc^?2; here, the circular rotation velocity of the Sun around the Galactic center is V₀ = 218 ± 6 km s^?1 kpc (for the adopted distance R₀ = 8.0 ± 0.2 kpc), while the Oort constants are A = 15.07 ± 0.25 km s^?1 kpc^?1 and B = ?12.17 ± 0.39 km s^?1 kpc^?1, p = 0.98 ± 0.08. The kinematics of Gaia TGAS stars with parallax errors more than 10% has been studied by invoking the distances from a paper by Astraatmadja and Bailer-Jones that were corrected for the Lutz–Kelker bias. We show that the second derivative of the angular velocity of Galactic rotation Ω'₀ = 0.864 ± 0.021 km s^?1 kpc^?3 is well determined from stars at a mean distance of 537 pc. On the whole, we have found that the distances of stars from the Gaia TGAS catalogue calculated using their trigonometric parallaxes do not require any additional correction factor.