Motions in wide pairs at various galactocentric distances

The motions of the components of wide binary stars in the regular Galactic gravitational field on time scales ~10¹⁰ yr at various Galactocentric distances R ₀ have been studied numerically. Near the Galactic center, the influence of the bar has been taken into account. The regions of restricted motions of the components in wide pairs have been found depending on the initial conditions: the magnitude of the relative velocity of the components, their mutual distance, and the inclination of the relative velocity vector to the Galactic plane. The shape and sizes of these regions are shown to depend significantly on R ₀: the sizes of the region of initial conditions corresponding to restricted motions increase with R ₀. Profound changes in the eccentricity of the binary orbit occur at inclinations close to 90?, which can lead to close approaches of the stars with a pericenter distance less than 1 AU. For retrograde motions (the binary rotates in a direction opposite to the Galactic rotation), there are elongated branches of the region of restricted motions in all cases, which extend at least to 10 pc in some cases.     

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.     

Kinematics of Tycho-2 red giant clump stars

Based on the Ogorodnikov-Milne model, we analyze the proper motions of 95 633 red giant clump (RGC) stars from the Tycho-2 Catalogue. The following Oort constants have been found: A = 15.9 ± 0.2 km s^?1 kpc^?1 and B = ?12.0±0.2 km s^?1 kpc^?1. Using 3632 RGC stars with known proper motions, radial velocities, and photometric distances, we show that, apart from the star centroid velocity components relative to the Sun, only the model parameters that describe the stellar motions in the XY plane differ significantly from zero. We have studied the contraction (a negative K effect) of the system of RGC stars as a function of their heliocentric distance and elevation above the Galactic plane. For a sample of distant (500–1000 pc) RGC stars located near the Galactic plane (|z| < 200 pc) with an average distance of d = 0.7 kpc, the contraction velocity is shown to be Kd = ?3.5 ±0.9 km s^?1; a noticeable vertex deviation, l _xy = 9 _· ^o 1 ± 0 _· ^o 5, is also observed for them. For stars located well above the Galactic plane (|z| ≥200 pc), these effects are less pronounced, Kd = ?1.7 ± 0.5 km s^?1 and l _xy = 4 _· ^o 9 ± 0 _· ^o 6. Using RGC stars, we have found a rotation around the Galactic X axis directed toward the Galactic center with an angular velocity of ?2.5 ± 0.3 km s^?1 kpc^?1, which we associate with the warp of the Galactic stellar-gaseous disk.     

Search for kinematic siblings of the sun based on data from the XHIP catalog

From the XHIP catalogue, we have selected 1872 F-G-K stars with relative parallax measurement errors <20% and absolute values of their space velocities relative to the Sun <15 km s^?1. For all these stars, we have constructed their Galactic orbits for 4.5 Gyr into the past using an axisymmetric Galactic potential model with allowance made for the perturbations from the spiral density wave. Parameters of the encounter with the solar orbit have been calculated for each orbit. We have detected three new stars whose Galactic orbits were close to the solar one during a long time interval in the past. These stars are HIP 43852, HIP 104047, and HIP 112158. The spectroscopic binary HIP 112158 is poorly suited for the role of a kinematic sibling of the Sun by its age and spectroscopic characteristics. For the single star HIP 43852 and the multiple system HIP 104047, this role is quite possible. We have also confirmed the status of our previously found candidates for close encounters, HIP 47399 and HIP 87382. The star HIP 87382 with a chemical composition very close to the solar one is currently the most likely candidate, because it persistently shows close encounters with the Sun on time scales of more than 3 Gyr when using various Galactic potential models both without and with allowance made for the influence of the spiral density wave.     

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.     

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.     

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.     

Kinematic control of the inertiality of the system of Tycho-2 and UCAC2 stellar proper motions

Based on the Ogorodnikov-Milne model, we analyze the proper motions of Tycho-2 and UCAC2 stars. We have established that the model component that describes the rotation of all stars under consideration around the Galactic y axis differs significantly from zero at various magnitudes. We interpret this rotation found using the most distant stars as a residual rotation of the ICRS/Tycho-2 system relative to the inertial reference frame. For the most distant (d≈900 pc) Tycho-2 and UCAC2 stars, the mean rotation around the Galactic y axis has been found to be M ₁₃ ^? =?0.37±0.04 mas yr^?1. The proper motions of UCAC2 stars with magnitudes in the range 12–15^m are shown to be distorted appreciably by the magnitude equation in μ_α cos δ, which has the strongest effect for northern-sky stars with a coefficient of ?0.60±0.05 mas yr^?1 mag^?1. We have detected no significant effect of the magnitude equation in the proper motions of UCAC2 stars brighter than ≈11^m.     

Distortion of the stellar velocity field parameters due to systematic variations of parallaxes over the celestial sphere

We study the effect of systematic variations in stellar parallaxes over the celestial sphere on the results of a kinematic analysis of stellar proper motions. Our approach is based on the representation of stellar parallaxes by scalar spherical harmonics and on the decomposition of stellar proper motions into a system of vector spherical harmonics. We derive theoretical relations that relate the coefficients of the decomposition of stellar proper motions into toroidal and spheroidal harmonics to the coefficients of the decomposition of stellar parallaxes into scalar spherical harmonics. We have established that the systematic variations of parallaxes over the celestial sphere distort all parameters of the linear Ogorodnikov-Milne model and can be responsible for the appearance of beyond-the-model harmonics. We have performed a kinematic analysis of the proper motions of blue-white and red giants based on Hipparcos data. The parallaxes of blue-white giants show a strong dependence on Galactic latitude (with predominant contraction along the Galactic equator). In contrast, the deviations of the parallaxes from the mean for red giants are localized only in two regions of the celestial sphere. For these samples, the effect of parallax variations over the celestial sphere on kinematic parameters has turned out to be comparable to their rms errors. The global solutions performed using both samples have revealed strong beyond-the-model kinematic effects described by second-order toroidal harmonics and third-order spheroidal harmonics. Using the solutions performed separately in the northern and southern Galactic hemispheres, we have established that not the systematic variations of parallaxes over the celestial sphere but the retardation of Galactic rotation with increasing distance of stars from the principal Galactic plane is mainly responsible for the appearance of these harmonics. Based on these samples of stars, we have estimated the magnitude of the vertical Galactic rotation velocity gradient to be 18.0±2.9 and 22.7±2.2 km s^?1 kpc^?1, respectively.     

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.     

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.     

The formation of a bound star cluster: from the Orion nebula cluster to the Pleiades

Direct N -body calculations are presented of the formation of Galactic clusters using GasEx , which is a variant of the code Nbody6 . The calculations focus on the possible evolution of the Orion nebula cluster (ONC) by assuming that the embedded OB stars explosively drove out 2/3 of its mass in the form of gas about 0.4 Myr ago. A bound cluster forms readily and survives for 150 Myr despite additional mass loss from the large number of massive stars, and the Galactic tidal field. This is the very first time that cluster formation is obtained under such realistic conditions. The cluster contains about 1/3 of the initial 10⁴ stars, and resembles the Pleiades cluster to a remarkable degree, implying that an ONC-like cluster may have been a precursor of the Pleiades. This scenario predicts the present expansion velocity of the ONC, which will be measurable by upcoming astrometric space missions. These missions should also detect the original Pleiades members as an associated expanding young Galactic-field subpopulation. The results arrived at here suggest that Galactic clusters form as the nuclei of expanding OB associations.
The results have wide implications, also for the formation of globular clusters and the Galactic-field and halo stellar populations. In view of this, the distribution of binary orbital periods and the mass function within and outside the model ONC and Pleiades is quantified, finding consistency with observational constraints. Advanced mass segregation is evident in one of the ONC models. The calculations show that the primordial binary population of both clusters could have been much the same as is observed in the Taurus–Auriga star-forming region. The computations also demonstrate that the binary proportion of brown dwarfs is depleted significantly for all periods, whereas massive stars attain a high binary fraction.     

Red giant clump in the Tycho-2 catalogue

The Tycho-2 proper motions and Tycho-2 and 2MASS photometry are used to select 97348 red giant clump (RGC) stars. The interstellar extinction and photometric distance are calculated for each of the stars. The selected stars are shown to form a selection-unbiased sample of RGC stars within about 350 pc of the Sun with the addition of more distant stars. The distribution of the selected stars in space and their motion are consistent with the assumption that the RGC contains Galactic disk stars with various ages and metallicities, including a significant fraction of stars younger than 1 Gyr with masses of more than 2M _⊙. These young stars show differences of their statistical characteristics from those of older RGC stars, including differences in the variations of their distribution density with distance from the Galactic plane and in the dispersion of their velocities found using radial velocities and proper motions. The Sun has been found to rise above the Galactic plane by 13 ± 1 pc. The distribution density of the stars under consideration in space is probably determined by the Local Spiral Arm and the distribution of absorbing matter in the plane of the Gould Belt.     

Velocity Space of Galactic O-B Stars

Based on the Hipparcos proper motions and available radial velocity data of O-B stars, we have re-examined the local kinematical structure of the young disk population of-1500 O-B stars not including the Gould-belt stars. A systematic warping motion of the stars about the direction to the Galactic center has been reconfirmed. A negative K-term implying a systematic contraction of stars in the solar vicinity has been detected. Two different distance scales are used in order to find out their impact on the kinematical parameters, and we conclude that the adopted distance scale plays an important role in characterizing the kinematical parameters at the present level of the measurement uncertainty.     

Proper motions,absolute magnitudes and spatial distribution of Zirconium stars

The computational algorithm to determine the the proper motions of Zirconium stars on the basis of catalogues “Carte du Ciel” and on the recent photographic observations carried out with the 70cm Abastumani meniscus telescope is presented. It allowed to determine the proper motions of 288 stars in the region around α Per with a rms error of ± 0,004 arcsec/yr. Applying the method proper motions of 74 Zirconium stars and 146 control stars have been obtained. The error of proper motions obtained for the North Zone (δ > −2°) 109 AGK₃ control stars is ± 0.006 arcsec/yr. On the basis of proper motins absolute magnitudes were separately calculated for the MCLPZS and LASZS. For the MCLPZS the average absolute visual magnitude at maximum, corresponding to the mean period of P = 350 days, equals −3^ϕm.9. For the LASZS the mean absolute visual magnitude, corresponding to the apparent median ones equals −1^ϕm.9. Low luminosity (M_v = −1^ϕmϕ9) Zirconium stars escape rather far (at a distance of up to 2 kpc) to the South from the Galactic plane into the region l ∼ 240 – 260°, where its assumed to be a connection with the Large Magellanic Cloud (LMC) begins to appear. Low luminosity Zirconium stars are weakly correlated with position of the Galaxy spiral arms. The MCLPZS show a somewhat other distribution.     

Proper motions and CCD photometry of stars in the region of the open cluster NGC 6866

We present the results of our comprehensive study of the Galactic open star cluster NGC 6866. The positions of stars in the investigated region have been obtained with the “Fantasy” automatic measuring machine from 10 plates of the normal astrograph at the Pulkovo Astronomical Observatory. The size of the investigated field is 40′ × 40′, the limiting magnitude is B ∼ 16_· ^m 6, and the maximum epoch difference is 79 yr. For 1202 field stars, we have determined the relative proper motions with an rms error of 2.5 mas yr⁻¹. Out of them, 423 stars may be considered cluster members with a probability P > 70% according to the astrometric criterion. Photometric diagrams have been used as an additional criterion. We have performed two-color BV CCD photometry of stars with the Pulkovo ZA-320M mirror astrograph. The U magnitudes from the literature have also been used to construct the two-color diagrams. A total of 267 stars have turned out to be members of NGC 6866 according to the two criteria. We present refined physical parameters of the cluster and its age estimate (5.6 × 10⁸ yr). The cluster membership of red and blue giants, variable, double, and multiple stars is considered. We have found an almost complete coincidence of the positions of one of the stars in the region (a cluster nonmember) and a soft X-ray source in the ROSAT catalog. The “Fantasy” automatic measuring machine is described in the Appendix.     

Kinematic parameters of the galactic spiral pattern from data on open star clusters and OB stars

Data on the positions, radial velocities, and proper motions of open star clusters and OB stars are used to obtain the rotation curve of the Galaxy fitted by a polynomial in inverse powers of the distances from the Galactic rotation axis. We determine the locations of the corotation region and the inner and outer Lindblad resonances using a previously estimated pattern speed. Based on data for objects of the Carina-Sagittarius and Orion arms, we have determined the distortion amplitudes of the velocity field of the Galactic disk, ?_R = ?3.97±4.79 km s^?1 and f_θ=+13.27±2.57 km s^?1.     

Modeling the spatial distribution of neutron stars in the Galaxy

In this paper we investigate the space and velocity distributions of old neutron stars (aged 10⁹ to 10¹⁰?yr) in our Galaxy. Galactic old Neutron Stars (NSs) population fills a torus-like area extending to a few tens kiloparsecs above the galactic plane. The initial velocity distribution of NSs is not well known, in this work we adopt a three component initial distribution, as given by the contribution of kick velocities, circular velocities and Maxwellian velocities. For the spatial initial distribution we use a Γ function. We then use Monte Carlo simulations to follow the evolution of the NSs under the influence of the Paczyński Galactic gravitational potential. Our calculations show that NS orbits have a very large Galactic radial expansion and that their radial distribution peak is quite close to their progenitors’ one. We also study the NS vertical distribution and find that it can well be described by a double exponential low. Finally, we investigate the correlation of the vertical and radial distribution and study the radial dependence of scale-heights.     

Light scattering by moving dust grains in the immediate vicinity of young stars

We consider the problem of the distortion of the photospheric spectrum for a young star as its light is scattered in the inner accretion disk in the dust grain evaporation region. In T Tauri stars, this region is at a distance of the order of several stellar radii and is involved in the large-scale motions of matter with velocities of ～100 km s^?1 or higher. The light scattering in such a medium causes the frequency of the scattered radiation to be shifted due to the Doppler effect. We analyze the influence of this effect on the absorption line profiles in the spectra of T Tauri stars using classical results of the theory of radiative transfer. We consider two models of a scattering medium: (i) a homogeneous cylindrical surface and (ii) a cylindrical surface with an azimuth-dependent height (such conditions take place during the accretion of matter onto a star with an oblique magnetic dipole). We show that in the first case, the scattering of the photospheric radiation causes the absorption lines to broaden. If the motion of the circumstellar matter in the dust evaporation region is characterized by two velocity components, then the line profile of the scattered radiation is asymmetric, with the pattern of the asymmetry depending on the direction of the radial velocity. In the second case, the scattered radiation can cause periodic shifts of the absorption line centroid, which can be perceived by an observer as periodic radial-velocity variations in the star. We suggest that precisely this effect is responsible for the low-amplitude radial-velocity variations with periods close to the stellar rotation periods that have recently been found in some of the T Tauri stars.     

Galactic kinematics from a sample of young massive stars

Based on published sources, we have created a kinematic database on 220 massive (> 10 M _⊙) young Galactic star systems located within ≤3 kpc of the Sun. Out of them, ≈100 objects are spectroscopic binary and multiple star systems whose components are massive OB stars; the remaining objects are massive Hipparcos B stars with parallax errors of no more than 10%. Based on the entire sample, we have constructed the Galactic rotation curve, determined the circular rotation velocity of the solar neighborhood around the Galactic center at R ₀ = 8kpc, V ₀ = 259±16 km s^?1, and obtained the following spiral density wave parameters: the amplitudes of the radial and azimuthal velocity perturbations f _R = ?10.8 ± 1.2 km s^?1 and f _θ = 7.9 ± 1.3 km s^?1, respectively; the pitch angle for a two-armed spiral pattern i = ?6.0° ± 0.4°, with the wavelength of the spiral density wave near the Sun being λ = 2.6 ± 0.2 kpc; and the radial phase of the Sun in χ _⊙ = ?120° ± 4°. We show that such peculiarities of the Gould Belt as the local expansion of the system, the velocity ellipsoid vertex deviation, and the significant additional rotation can be explained in terms of the density wave theory. All these effects decrease noticeably once the influence of the spiral density wave on the velocities of nearby stars has been taken into account. The influence of Gould Belt stars on the Galactic parameter estimates has also been revealed. Eliminating them from the kinematic equations has led to the following new values of the spiral density wave parameters: f _θ = 2.9 ± 2.1 km s^?1 and χ _⊙ = ?104° ± 6°.