Modeling the vertical distribution of the Milky Way's flat subsystem objects


 This paper is an initial stage of consideration of the general problem of joint modeling of the vertical structure of a Galactic flat subsystem and the average surface of the disk of the Galaxy, taking into account the natural and measurement dispersions. We approximate the average surface of the Galactic disk in the region covered by the data with a general (polynomial) model and determine its parameters by minimizing the squared deviations of objects along the normal to the model surface. The smoothness of the model, i.e. its order~$n$, is optimized. An outlier elimination algorithm is applied. The developed method allows us to simultaneously identify significant details of the Galactic warping and estimate the offset~$z_{\sun}$ of the Sun relative to the average (in general, non-flat) surface of the Galactic disk and the vertical scale of the object system under consideration for an arbitrary area of the disk covered by data. The method is applied to data on classical Cepheids (\citeauthor{Berdnikov+2000}, \citeauthor{Mel'nik+2015}). Significant local extremes of the average disk surface model were found based on Cepheid data: the minimum in the first Galactic quadrant and the maximum in the second. A well-known warp (lowering of the disk surface) in the third quadrant has been confirmed. The optimal order of the model describing all these warping details was found to be $n_\text{o}=4$. The local (for a small neighborhood of the Sun, $n_\text{o}=0$) estimate of $z_{\sun} = 28.1 \pm \left.6.1\right|_{\text{stat.}}\left.{}^{+1.3}_{-1.3}\right|_{\text{cal.}}$\,pc is close to the non-local (taking into account warping, $n_\text{o}=4$) $z_{\sun} = 27.1 \pm \left.8.8\right|_{\text{stat.}}\left.{}^{+1.3}_{-1.2}\right|_{\text{cal.}}$\,pc (statistical and calibration uncertainties are indicated), which suggests that the proposed modeling method eliminates the influence of warping on the $z_{\sun}$ estimate. However, the non-local estimate of the vertical standard deviation of Cepheids $\sigma_{\rho} = 132.0 \pm \left.3.7\right|_{\text{stat.}}\left.{}^{+6.3}_{-5.9}\right|_{\text{cal.}}$\,pc differs significantly from the local $\sigma_{\rho} = \left.76.5 \pm 4.4\right|_{\text{stat.}}\left.{}^{+3.6}_{-3.4}\right|_{\text{cal.}}$\, pc, which means the need to introduce more complex models for the vertical distribution outside the Sun's vicinity.