Calibration of the period-luminosity relation (PLR) for Cepheids has always been one of the biggest goals of stellar astronomy. Among a considerable number of different approaches, the Baade-Becker-Wesselink (BBW) method stands in the foreground as one of the most universal and precise methods. We present a new realization of the BBW method which is considered to be the generalization of surface brightness (\citealt{Barnes+Evans+1976}) and \cite{Balona+1977} approaches first proposed by \cite{Rastorguev+Dambis+2010} and described in \cite{Paper1}. One of the main features of this method is using measured effective temperature variations to determine the main parameters of Cepheid, such as distance, radius, luminosity, colour excess, intrinsic colour. We apply this method to 45 Cepheids of Northern sky, for which multiphase temperature data are available. We take into account the effect of shock waves, whose presence in stellar atmosphere distorts the observational data and the calibrations used in this work. Within $0.0-0.87$ phase interval we derived PL relation $<M_V>_I=-(2.67 \pm 0.17)\cdot logP - (1.58 \pm 0.16)$. It was used to calculate the distances, rotation curve and kinematical parameters of the sample of 435 Cepheids with GAIA DR2 proper motions.

A new version of the Baade--Becker--Wesselink (hereafter BBW) method is proposed, based on direct spectroscopic measurements of effective temperatures of 45 northern Cepheids, made in different pulsating phases. By comparing the temperature estimates obtained from the calibration of effective temperature by normal color with real temperature measurements we were able not only to determine the color excess with an accuracy of 0.01 mag, but also to derive new color calibration of the effective temperature immediately for all available measurements with taking into account the differences in [Fe/H] and log g values: log Teff = 3.88 - 0.20 (B-V)o + 0.026 (B-V)o^2 + 0.009 log g - 0.010 (B-V)o log g - 0.051 [Fe/H] + 0.051 (B-V)o [Fe/H], which is accurate to about 1.1%. We also showed the complete identity of the two main versions of the BBW technique: surface brightness method proposed by Barnes and Evans (1976) and maximum likelyhood method of Balona (1977), refined later by Rastorguev and Dambis (2010).

V CVn is a semiregular variable star with a V-band amplitude of $\approx2$ mag. This star has an unusually high amplitude of polarimetric variability: up to 6 per cent. It also exhibits a prominent inverse correlation between the flux and the fraction of polarization and a substantial constancy of the angle of polarization. To figure out the nature of these features, we observed the object using the Differential Speckle Polarimetry at three bands centered on 550, 625 and 880 nm using the 2.5-m telescope of Sternberg Astronomical Institute. The observations were conducted on 20 dates distributed over three cycles of pulsation. We detected an asymmetric reflection nebula consisting of three regions and surrounding the star at the typical distance of 35 mas. The components of the nebula change their brightness with the same period as the star, but with significant and different phase shifts. We discuss several hypotheses that could explain this behavior.

V350 Sgr is a classical Cepheid suitable for mass determination. It has a hot companion which is prominent in the ultraviolet and which is not itself a binary. We have obtained two high resolution echelle spectra of the companion at orbital velocity maximum and minimum with the Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph (STIS) in the 1320 to 1510 \AA\/ region. By cross-correlating these spectra we obtained the orbital velocity amplitude of the companion with an uncertainty in the companion amplitude of 1.9 km sec$^{-1}$. This provides a mass ratio of the Cepheid to the companion of 2.1. The ultraviolet energy distribution of the companion provides the mass of the companion, yielding a Cepheid mass of 5.2 $\pm$ 0.3 M$_\odot$. This mass requires some combination of moderate main sequence core convective overshoot and rotation to match evolutionary tracks.

We cross-match objects from several different astronomical catalogs to determine the absolute proper proper motions of stars within the 30-arcmin radius fields of 115 Milky-Way globular clusters with the accuracy of 1--2~mas/yr. The proper motions are based on positional data recovered from the USNO-B1, 2MASS, URAT1, ALLWISE, UCAC5, and GAIA DR1 surveys with up to 10 positions spanning an epoch difference of up to $\sim$~65~years, and reduced to GAIA DR1 TGAS frame using UCAC5 as the reference catalog. Cluster members are photometrically identified by selecting horizontal- and red-giant branch stars on color-magnitude diagrams, and the mean absolute proper motions of the clusters with a typical formal error of $\sim$~0.4~mas/yr are computed by averaging the proper motions of selected members. The inferred absolute proper motions of clusters are combined with available radial-velocity data and heliocentric distance estimates to compute the cluster orbits in terms of the Galactic potential models based on Miyamoto and Nagai disk, Hernquist spheroid, and modified isothermal dark-matter halo (axisymmetric model without a bar) and the same model + rotating Ferre's bar (non-axisymmetric). Five distant clusters have higher-than-escape velocities, most likely due to large errors of computed transversal velocities, whereas the computed orbits of all other clusters remain bound to the Galaxy. Unlike previously published results, we find the bar to affect substantially the orbits of most of the clusters, even those at large Galactocentric distances, bringing appreciable chaotization, especially in the portions of the orbits close to the Galactic center, and stretching out the orbits of some of the thick-disk clusters.

Kinetic processes in fractal stellar media are analysed in terms of the approach developed in our earlier paper (Chumak \& Rastorguev, 2015) involving a generalization of the nearest neighbour and random force distributions to fractal media. Diffusion is investigated in the approximation of scale-dependent conditional density based on an analysis of the solutions of the corresponding Langevin equations. It is shown that kinetic parameters (time scales, coefficients of dynamic friction, diffusion, etc.) for fractal stellar media can differ significantly both qualitatively and quantitatively from the corresponding parameters for a quasi-uniform random media with limited fluctuations. The most important difference is that in the fractal case kinetic parameters depend on spatial scale length and fractal dimension of the medium studied. A generalized kinetic equation for stellar media (fundamental equation of stellar dynamics) is derived in the Fokker-Planck approximation with the allowance for the fractal properties of the spatial stellar density distribution. Also derived are its limit forms that can be used to describe small departures of fractal gravitating medium from equilibrium.

We applied the currently most comprehensive version of the statistical-parallax technique to derive kinematical parameters of the maser sample with 136 sources. Our kinematic model comprises the overall rotation of the Galactic disk and the spiral density-wave effects. We take into account the variation of radial velocity dispersion with Galactocentric distance. The best description of the velocity field is provided by the model with constant radial and vertical velocity dispersions, $(\sigma U0, \sigma W0) \approx (9.4 \pm 0.9~, 5.9 \pm 0.8)~ km/s$. We compute flat Galactic rotation curve over the Galactocentric distance interval from 3 to 15 kpc and find the local circular rotation velocity to be $ V_0 \approx (235-238)$~ km/s $\pm 7$~ km/s. We also determine the parameters of the four-armed spiral pattern (pitch angle $i \approx (-10.4 \pm 0.3)^\circ$ and the phase of the Sun $\chi_0 \approx (125 \pm 10) ^\circ$). The radial and tangential spiral perturbations are about $f_R \approx (-6.9 \pm 1.4)$~km/s, $f_\Theta \approx (+2.8 \pm 1.0$) ~km/s. The kinematic data yield a solar Galactocentric distance of $R_0 \approx (8.24 \pm 0.12)~kpc$. Based on rotation curve parameters and the asymmetric drift we Infer the exponential disk scale $H_D \approx (2.7 \pm 0.2)$ ~kpc under assumption of marginal stability of the intermediate-age disk, and finally we estimate the minimum local surface disk density, $\Sigma (R_0) > (26 \pm 3) ~ M_\odot pc^{-2}$.

Agekyan lambda-factor that accounts for the effect of multiple distant encounters with large impact factors is used for the first time to compute the diffusion coefficients in the velocity space of a stellar system. It is shown that in this case the cumulative effect - the total contribution of distant encounters to the change in the velocity of the test star - is finite, and the logarithmic divergence inherent to the classical description disappears, as also was earlier noted by Kandrup (1981). At the same time, the formulas for the diffusion coefficients, as before, contain the logarithm of the ratio of two independent scale factors that fully characterize the state of the stellar system: the average interparticle distance and the impact parameter of a close encounter. However, the physical meaning of this factor is no longer associated with the classical logarithmic divergence.

The nearest neighbor distribution (Chandrasekhar 1943) is generalized to fractal stellar systems.For such systems an asymptotic distribution of the magnitude of large random forces and a formula for the effective mean interparticle spacing are derived. It is shown that in the case of a power-law distribution of conditional density the derived asymptotic fully agrees with the results obtained in terms of a general approach. It is concluded that large random forces in a fractal stellar medium are due entirely to the nearest neighbors (clumps) located inside the sphere of the effective radius determined from the generalized Holtsmark distribution.

Stellar population and the interstellar gas-dust medium in the vicinity of the open star cluster vdB 130 are analysed using optical observations taken with the 6-m telescope of the SAO RAS and the 125-cm telescope of the SAI MSU along with the data of Spitzer and Herschel. Based on proper motions and BV and JHKs 2MASS photometric data, we select additional 36 stars as probable cluster members. Some stars in vdB 130 are classified as B stars. Our estimates of minimum colour excess, apparent distance modulus and the distance are consistent with young age (from 5 to 10 Myrs) of the cluster vdB 130. We suppose the large deviations from the conventional extinction law in the cluster direction, with $R_V$ ~ 4 - 5. The cluster vdB 130 appears to be physically related to the supershell around Cyg OB1, a cometary CO cloud, ionized gas, and regions of infrared emission. There are a few regions of bright mid-infrared emission in the vicinity of vdB 130. The largest of them is also visible on H-alpha and [SII] emission maps. We suggest that the infrared blobs that coincide in projection with the head of the molecular cloud are HII regions, excited by the cluster B-stars. Some signatures of a shock front are identified between these IR-bright regions.

We have examined high accuracy radial velocities of Cepheids to determine the binary frequency. The data are largely from the CORAVEL spectrophotometer and the Moscow version, with a typical uncertainty of $\leq1$~km~s$^{-1}$, and a time span from 1 to 20 years. A systemic velocity was obtained by removing the pulsation component using a high order Fourier series. From this data we have developed a list of stars showing no orbital velocity larger than $\pm1$~km~s$^{-1}$. The binary fraction was analyzed as a function of magnitude, and yields an apparent decrease in this fraction for fainter stars. We interpret this as incompleteness at fainter magnitudes, and derive the preferred binary fraction of $29\pm8$\% ( $20\pm6$\% per decade of orbital period) from the brightest 40 stars. Comparison of this fraction in this period range (1-20 years) implies a large fraction for the full period range. This is reasonable in that the high accuracy velocities are sensitive to the longer periods and smaller orbital velocity amplitudes in the period range sampled here. Thus the Cepheid velocity sample provides a sensitive detection in the period range between short period spectroscopic binaries and resolved companions. The recent identification of $\delta$ Cep as a binary with very low amplitude and high eccentricity underscores the fact that the binary fractions we derive are lower limits, to which other low amplitude systems will probably be added. The mass ratio (q) distribution derived from ultraviolet observations of the secondary is consistent with a flat distribution for the applicable period range (1 to 20 years).

We use the currently most complete collection of reliable Cepheid positions (565 stars) out to ~5 kpc based mostly on our photometric data to outline the spiral pattern of our Galaxy. We find the pitch-angle to be equal to 9--10 degrees with the most accurate estimate (i=9.5 +/-0.1 degrees) obtained assuming that the spiral pattern has a four-armed structure, and the solar phase angle in the spiral pattern to be chi_0 = 121+/-3 degrees. The pattern speed is found to be Omega_P=25.2+/-0.5km/s/kpc based on a comparison of the positions of the spiral arms delineated by Cepheids and maser sources and the age difference between these objects.

The kinematics and distribution of classical Cepheids within ~3 kpc from the Sun suggest the existence of the outer ring R1R2' in the Galaxy. The optimum value of the solar position angle with respect to the major axis of the bar, theta_b, providing the best agreement between the distribution of Cepheids and model particles is theta_b=37 +/- 13 degrees. The kinematical features obtained for Cepheids with negative Galactocentric radial velocity VR are consistent with the solar location near the descending segment of the outer ring R2. The sharp rise of extinction toward of the Galactic center can be explained by the presence of the outer ring R1 near the Sun.

The period - metallicity - WISE W1- and W2-band luminosity relations are derived for RR Lyrae stars based on WISE epoch photometry for 360 and 275 stars in 15 and 9 Galactic globular clusters, respectively. Our final relations have the form <MW1> = gamma(W1) - (2.381 +/- 0.097) log PF + (0.096 +/- 0.021)[Fe/H] and <MW2> = gamma(W2)-(2.269 +/- 0.127)log PF + (0.108 +/- 0.021)[Fe/H], where [Fe/H] values are on the scale of Carretta et al. (2009). We obtained two appreciably discrepant estimates for the zero points gamma(W1) and gamma(W2) of both relations: one based on a statistical-parallax analysis -- gamma(W1) = -0.829 +/- 0.093 and gamma(W2)=-0.776 +/- 0.093 and another, significantly brighter one, based on HST FGS trigonometric parallaxes -- gamma(W1, HST) =-1.150 +/- 0.077 and gamma(W2, HST) =-1.105 +/- 0.077. The period-metallicity-luminosity relations in the two bands yield highly consistent distance moduli for the calibrator clusters and the distance moduli computed using the W1- and W2-band relations with the HST zero points agree well with those computed by \citet{sollima} based on their derived period-metallicity-K-band luminosity relation whose zero point is tied to the HST trigonometric parallax of RR Lyrae itself (Delta DM0 = +0.04 and +0.06, respectively, with a scatter of only 0.06).

We use UCAC4 proper motions and WISE W1-band apparent magnitudes intensity-mean for almost 400 field RR Lyrae variables to determine the parameters of the velocity distribution of Galactic RR Lyrae population and constrain the zero points of the metallicity-<MV> relation and those of the period-metallicity-<MKs>-band and period-metallicity-<MW1>-band luminosity relations via statistical parallax. We find the mean velocities of the halo- and thick-disc RR Lyrae populations in the solar neighbourhood to be (U0(Halo), V0(Halo), W0(Halo)) = (-7 +/- 9, -214 +/- 10, -10 +/- 6) km/s and (U0(Disc), V0(Disc), W0(Disc)) =(-13 +/- 7, -37 +/- 6, -17 +/- 4) km/s, respectively, and the corresponding components of the velocity-dispersion ellipsoids, (sigma VR(Halo), sigma Vphi(Halo), sigma Vtheta(Halo)) = (153 +/- 9, 101 +/- 6, 96 +/- 5) km/s and (sigma VR(Disc), sigma Vphi(Disc), sigma Vtheta(Disc)) = (46 +/- 7, 37 +/- 5, 27 +/- 4) km/s, respectively. The fraction of thick-disc stars is estimated at 0.22 +/- 0.03. The corrected IR period-metallicity-luminosity relations are <MKs> = -0.769 +0.088 [Fe/H]- 2.33 mathoprm log PF and <MW1> = -0.825 + 0.088 [Fe/H] -2.33 mathoprm log PF, and the optical metallicity-luminosity relation, [Fe/H]-<MV>, is <MV> = +1.094 + 0.232 [Fe/H], with a standard error of +/- 0.089, implying an LMC distance modulus of 18.32 +/- 0.09, a solar Galactocentric distance of 7.73 +/- 0.36 kpc, and the M31 and M33 distance moduli of DM(M31) = 24.24 +/- 0.09 (D = 705 +/- 30 kpc) and DM(M33) = 24.36 +/- 0.09 (D = 745 +/- 31 kpc), respectively. Extragalactic distances calibrated with our RR Lyrae star luminosity scale imply a Hubble constant of ~80 km/s/Mpc. Our results suggest marginal prograde rotation for the population of halo RR Lyraes in the Milky Way.

The BBW method remains one of most demanded tool to derive full set of Cepheid astrophysical parameters. Surface brightness version of the BBW technique was preferentially used during last decades to calculate Cepheid radii and to improve PLC relations. Its implementation requires a priory knowledge of Cepheid reddening value. We propose a new version of the Baade--Becker--Wesselink technique, which allows one to independently determine the colour excess and the intrinsic colour of a radially pulsating star, in addition to its radius, luminosity, and distance. It is considered to be a generalization of the Balona light curve modelling approach. The method also allows the function F(CI_0) = BC + 10 log Teff for the class of pulsating stars considered to be calibrated. We apply this technique to a number of classical Cepheids with very accurate light and radial-velocity curves. The new technique can also be applied to other pulsating variables, e.g. RR Lyraes. We discuss also possible dependence of the projection factor on the pulsation phase.

The BVRcIc CCD photometry in the fields of six open clusters toward the Perseus spiral arm is presented. These data along with JHKs magnitudes taken from 2MASS catalog, have been used to determine cluster's ages, distances and color excesses. In addition, the gaps in mass function of Be 97, King 12 and NGC 7788 clusters have been revealed in mass intervals of 1.3-1.5, 1.4-1.6 and 1.5-1.7 solar masses, respectively.

We propose a new version of the Baade--Becker--Wesselink technique, which allows one to independently determine the colour excess and the intrinsic colour of a radially pulsating star, in addition to its radius, luminosity, and distance. It is considered to be a generalization of the Balona approach. The method also allows the function F(CI) = BC + 10 log (Teff) for the class of pulsating stars considered to be calibrated. We apply this technique to a number of classical Cepheids with very accurate light and radial-velocity curves and with bona fide membership in open clusters (SZ Tau, CF Cas, U Sgr, DL Cas, GY Sge), and find the results to agree well with the reddening estimates of the host open clusters. The new technique can also be applied to other pulsating variables, e.g. RR Lyrae and RV Tauri.

Current state of the variable stars studies, including halo and globular-cluster variables is briefly discussed, as well as the halo kinematics and modern problems of the Universal distance scale. I propose new approach to the projection factor (PF) calculation which takes into account its variation with limb darkening, photosphere velocity and the intrinsic line width of the spectrograph used for radial velocity measurements. Simple analytic approximation for PF is given, useful for estimates of the pulsating variables radii via Baade-Wesselink-Balona technique.

An analysis of the residual-velocity field of OB associations within 3 kpc of the Sun has revealed periodic variations in the radial residual velocities along the Galactic radius vector with a typical scale length of lambda=2.0(+/-0.2) kpc and a mean amplitude of fR=7(+/-1) km/s. The fact that the radial residual velocities of almost all OB-associations in rich stellar-gas complexes are directed toward the Galactic center suggests that the solar neighborhood under consideration is within the corotation radius. The azimuthal-velocity field exhibits a distinct periodic pattern in the region 0<l<180 degrees, where the mean azimuthal-velocity amplitude is ft=6(+/-2) km/s. There is no periodic pattern of the azimuthal-velocity field in the region 180<l<360 degrees. The locations of the Cygnus arm, as well as the Perseus arm, inferred from an analysis of the radial- and azimuthal-velocity fields coincide. The periodic patterns of the residual-velocity fields of Cepheids and OB associations share many common features.