It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$\alpha$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$\alpha$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$\alpha$ lines are delayed with a velocity gradient up to $\Delta v \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$\alpha$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere.

Astrobots are robotic artifacts whose swarms are used in astrophysical studies to generate the map of the observable universe. These swarms have to be coordinated with respect to various desired observations. Such coordination are so complicated that distributed swarm controllers cannot always coordinate enough astrobots to fulfill the minimum data desired to be obtained in the course of observations. Thus, a convergence verification is necessary to check the suitability of a coordination before its execution. However, a formal verification method does not exist for this purpose. In this paper, we instead use machine learning to predict the convergence of astrobots swarm. In particular, we propose a weighted $k$-NN-based algorithm which requires the initial status of a swarm as well as its observational targets to predict its convergence. Our algorithm learns to predict based on the coordination data obtained from previous coordination of the desired swarm. This method first generates a convergence probability for each astrobot based on a distance metric. Then, these probabilities are transformed to either a complete or an incomplete categorical result. The method is applied to two typical swarms including 116 and 487 astrobots. It turns out that the correct prediction of successful coordination may be up to 80% of overall predictions. Thus, these results witness the efficient accuracy of our predictive convergence analysis strategy.

Robotic fiber positioners play a vital role in the generation of massive spectroscopic surveys. The more complete a positioners set is coordinated, the more information its corresponding spectrograph receives during an observation. The complete coordination problem of positioners sets is studied in this paper. We first define the local and the global completeness problems and determine their relationship. We then propose a new artificial potential field according to which the convergences of a positioner and its neighboring positioners are cooperatively taken into account. We also discover the required condition for a complete coordination. We finally explain how the modifications of some of the parameters of a positioners set may resolve its incompleteness coordination scenarios. We verify our accomplishments using simulations.

Robotics have recently contributed to cosmological spectroscopy to automatically obtain the map of the observable universe using robotic fiber positioners. For this purpose, an assignment algorithm is required to assign each robotic fiber positioner to a target associated with a particular observation. The assignment process directly impacts on the coordination of robotic fiber positioners to reach their assigned targets. In this paper, we establish an optimal target assignment scheme which simultaneously provides the fastest coordination accompanied with the minimum of colliding scenarios between robotic fiber positioners. In particular, we propose a cost function by whose minimization both of the cited requirements are taken into account in the course of a target assignment process. The applied simulations manifest the improvement of convergence rates using our optimal approach. We show that our algorithm scales the solution in quadratic time in the case of full observations. Additionally, the convergence time and the percentage of the colliding scenarios are also decreased in both supervisory and hybrid coordination strategies.

Aims: From a coherent interpretation of the linear polarisation detected in the spectral lines of the Mira star $\chi$ Cyg, we derive information about the dynamics of the stellar photosphere, including pulsation. Methods: From spectropolarimetric observations of $\chi$ Cyg, we perform careful analysis of the polarisation signals observed in atomic and molecular lines, both in absorption and emission, using the radiative transfer in the polarisation context, through two mechanisms: intrinsic polarisation and continuum depolarisation. We also explain the observed line doubling phenomenon in terms of an expanding shell in spherical geometry, which allows us to pinpoint the coordinates over the stellar disk with enhanced polarisation. Results: We find that the polarised spectrum of $\chi$ Cyg is dominated by intrinsic polarisation, with a negligible continuum depolarisation. The observed polarised signals can only be explained by assuming that this polarisation is locally enhanced by velocity fields. During the pulsation, radial velocities are not homogeneous over the disk. We map these regions of enhanced velocities. Conclusions: We have set an algorithm to distinguish in any stellar spectra of linear polarisation the origin of this polarisation and the way to increase signal by coherently adding many lines with an appropriated weight. Applied to the Mira star $\chi$ Cyg, we reached the unexpected result that during the pulsation, velocities are radial but not homogeneous over the disk. The reason for these local velocity enhancements are probably related to the interplay between the atmospheric pulsation dynamics and the underlying stellar convection.

RRab stars are large amplitude pulsating stars in which the pulsation is associated with strong shock wake propagating in the atmosphere. The objective of this study is to provide a general overview of the dynamical structure of the atmosphere occurring over a typical pulsation cycle. We report new high-resolution observations with high time resolution of H$\alpha$ and sodium lines in the brightest RR Lyrae star of the sky: RR Lyr (HD 182989). A detailed analysis of line profile variations over the whole pulsation cycle is performed to understand the dynamical structure of the atmosphere. The main shock wave appears when it exits from the photosphere at $\varphi\simeq0.89$, i.e., when the main H$\alpha$ emission is observed. Whereas the acceleration phase of the shock is not observed, a significant deceleration of the shock front velocity is clearly present. The radiative stage of the shock wave is short: $4\%$ of the pulsation period ($0.892<\varphi<0.929$). A Mach number $M>10$ is required to get such a radiative shock. The sodium layer reaches its maximum expansion well before that of H$\alpha$ ($\Delta\varphi=0.135$). Thus, a rarefaction wave is induced between the H$\alpha$ and sodium layers. A strong atmospheric compression occurring around $\varphi=0.36$, which produces the third H$\alpha$ emission, takes place in the highest part of the atmosphere. The region located lower in the atmosphere where the sodium line is formed is not involved. The amplification of gas turbulence seems mainly due to strong shock waves propagating in the atmosphere rather than to the global compression of the atmosphere caused by the pulsation. It has not yet been clearly established whether the microturbulence velocity increases or decreases with height in the atmosphere. Furthermore, it seems very probable that an interstellar component is visible within the sodium profile.

The so-called H$\alpha$ third emission occurs around pulsation phase $\varphi$=0.30. It has been observed for the first time in 2011 in some RR Lyrae stars. The emission intensity is very weak, and its profile is a tiny persistent hump in the red side-line profile. We report the first observation of the H$\alpha$ third emission in RR Lyr itself (HD 182989), the brightest RR Lyrae star in the sky. New spectra were collected in 2013-2014 with the Aurelie}spectrograph (resolving power R=22$\,$700, T152, Observatoire de Haute-Provence, France) and in 2016-2017 with the eShel spectrograph (R=11$\,$000, T035, Observatoire de Chelles, France). In addition, observations obtained in 1997 with the Elodie spectrograph (R=42$\,$000, T193, Observatoire de Haute-Provence, France) were reanalyzed. The H$\alpha$ third emission is clearly detected in the pulsation phase interval $\varphi$=0.188-0.407, that is, during about 20% of the period. Its maximum flux with respect to the continuum is about 13%. The presence of this third emission and its strength both seem to depend only marginally on the Blazhko phase. The physical origin of the emission is probably due to the infalling motion of the highest atmospheric layers, which compresses and heats the gas that is located immediately above the rising shock wave. The infalling velocity of the hot compressed region is supersonic, almost 50 km$\cdot$s$^{-1}$, while the shock velocity may be much lower in these pulsation phases. When the H$\alpha$ third emission appears, the shock is certainly no longer radiative because its intensity is not sufficient to produce a blueshifted emission component within the H$\alpha$ profile. At phase $\varphi$=0.40, the shock wave is certainly close to its complete dissipation in the atmosphere.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

Massive spectroscopic survey are becoming trendy in astrophysics and cosmology, as they can address new fundamental knowledge such as Galactic Archaeology and probe the nature of the mysterious Dark Energy. To enable massive spectroscopic surveys, new technology are being developed to place thousands of optical fibers at a given position on a focal plane. These technology needs to be: 1) accurate, with micrometer positional accuracy; 2) fast to minimize overhead; 3) robust to minimize failure; and 4) low cost. In this paper we present the development of a new 8-mm in diameter fiber positionner robot using two 4mm DC-brushless gearmotors, developed in the context of the Dark Energy Spectroscopic Instrument. This development was conducted by a Spanish-Swiss (ES-CH) team led by the Instituto de F\'isica Te\'orica (UAM-CSIC) and the Laboratoire d'Astrophysique (EPFL), in collaboration with the AVS company in Spain and the Faulhaber group (MPS & FAULHABER-MINIMOTOR) in Switzerland.

The next generation of large-scale spectroscopic survey experiments such as DESI, will use thousands of fiber positioner robots packed on a focal plate. In order to maximize the observing time with this robotic system we need to move in parallel the fiber-ends of all positioners from the previous to the next target coordinates. Direct trajectories are not feasible due to collision risks that could undeniably damage the robots and impact the survey operation and performance. We have previously developed a motion planning method based on a novel decentralized navigation function for collision-free coordination of fiber positioners. The navigation function takes into account the configuration of positioners as well as their envelope constraints. The motion planning scheme has linear complexity and short motion duration (~2.5 seconds with the maximum speed of 30 rpm for the positioner), which is independent of the number of positioners. These two key advantages of the decentralization designate the method as a promising solution for the collision-free motion-planning problem in the next-generation of fiber-fed spectrographs. In a framework where a centralized computer communicates with the positioner robots, communication overhead can be reduced significantly by using velocity profiles consisting of a few bits only. We present here the discretization of velocity profiles to ensure the feasibility of a real-time coordination for a large number of positioners. The modified motion planning method that generates piecewise linearized position profiles guarantees collision-free trajectories for all the robots. The velocity profiles fit few bits at the expense of higher computational costs.

In the large-scale, Dark Energy Spectroscopic Instrument (DESI), thousands of fiber positioners will be used. Those are robotic positioners, with two axis, and having the size of a pen. They are tightly packed on the focal plane of the telescope. Dedicated micro-robots have been developed and they use 4mm brushless DC motors. To simplify the implementation and reduce the space occupancy, each actuator will integrate its own electronic control board. This board will be used to communicate with the central trajectory generator, manage low level control tasks and motor current feeding. In this context, we present a solution for a highly compact electronic. This electronic is composed of two layers. The first is the power stage that can drive simultaneously two brushless motors. The second one consists of a fast microcontroller and deals with different control tasks: communication, acquisition of the hall sensor signals, commutation of the motors phases, and performing position and current regulation. A set of diagnostic functions are also implemented to detect failure in the motors or the sensors, and to sense abnormal load change that may be the result of two robots colliding.

In order to complete the knowledge of the magnetic field and of its influence during the transition from Asymptotic Giant Branch to Planetary Nebulae stages, we have undertaken a search for magnetic fields at the surface of Mira stars. We used spectropolarimetric observations, collected with the Narval instrument at TBL, in order to detect - with Least Squares Deconvolution method - a Zeeman signature in the visible part of the spectrum. We present the first spectropolarimetric observations of the S-type Mira star chi Cyg, performed around its maximum light. We have detected a polarimetric signal in the Stokes V spectra and we have established its Zeeman origin. We claim that it is likely to be related to a weak magnetic field present at the photospheric level and in the lower part of the stellar atmosphere. We have estimated the strength of its longitudinal component to about 2-3 Gauss. This result favors a 1/r law for the variation of the magnetic field strength across the circumstellar envelope of chi Cyg. This is the first detection of a weak magnetic field at the stellar surface of a Mira star and we discuss its origin in the framework of shock waves periodically propagating throughout the atmosphere of these radially pulsating stars. At the date of our observations of chi Cyg, the shock wave reaches its maximum intensity, and it is likely that the shock amplifies a weak stellar magnetic field during its passage through the atmosphere. Without such an amplification by the shock, the magnetic field strength would have been too low to be detected. For the first time, we also report strong Stokes Q and U signatures (linear polarization) centered onto the zero velocity (i.e., at the shock front position). They seem to indicate that the radial direction would be favored by the shock during its propagation throughout the atmosphere.

Some of the next generation massive spectroscopic survey projects, such as DESI and PFS, plan to use thousands of fiber positioner robots packed at a focal plane to quickly move in parallel the fiber-ends from the previous to the next target points. The most direct trajectories are prone to collision that could damage the robots and impact the survey operation. We thus present here a motion planning method based on a novel decentralized navigation function for collision-free coordination of fiber positioners. The navigation function takes into account the configuration of positioners as well as the actuator constraints. We provide details for the proof of convergence and collision avoidance. Decentralization results in linear complexity for the motion planning as well as dependency of motion duration with respect to the number of positioners. Therefore the coordination method is scalable for large-scale spectrograph robots. The short in-motion duration of positioner robots (~2.5 seconds using typical actuator constraints), will thus allow the time dedicated for observation to be maximized.

In the spectra of pulsating stars, especially Mira stars, the detection of intense hydrogen emission lines has been explained by a radiative shock wave, periodically propagating throughout the atmosphere. Previous observation of the Mira star omicron Ceti around a bright maximum of light led to the detection of a strong linear polarization associated to Balmer emissions, although the origin of this phenomenon is not fully explained yet. With the help of spectropolarimetry, we propose to investigate the nature of shock waves propagating throughout the stellar atmosphere and present, for omicron Ceti (the prototype of Mira stars), a full observational study of hydrogen emission lines formed in the radiative region of such a shock. Using the instrument NARVAL, we performed a spectropolarimetric monitoring of omicron Ceti during three consecutive pulsation cycles. All Stokes parameters were systematically collected, with a particular emphasis on the maxima of luminosity, i.e. when a radiative shock wave is supposed to emerge from the photosphere and propagate outward. On Balmer lines, over a large part of the luminosity cycle, we report detections in Stokes spectra which are evolving with time. These signatures appear to be strongly correlated to the presence of an intense shock wave responsible for the hydrogen emission lines. We establish that those lines are polarized by a process inherent to the mechanism responsible for the emission line formation: the shock wave itself. Two mechanisms are considered: a global one that implies a polarization induced by convective cells located around the photosphere and a local one that implies a charge separation due to the passage of the shock wave, inducing an electrical current. Combined with the existing turbulence, this may generate a magnetic field, hence polarization.

The projection factor is a key quantity for the interferometric Baade-Wesselink (hereafter IBW) and surface-brightness (hereafter SB) methods of determining the distance of Cepheids. Indeed, it allows a consistent combination of angular and linear diameters of the star. We aim to determine consistent projection factors that include the dynamical structure of the Cepheids' atmosphere. Methods. Hydrodynamical models of delta Cep and l Car have been used to validate a spectroscopic method of determining the projection factor. This method, based on the amplitude of the radial velocity curve, is applied to eight stars observed with the HARPS spectrometer. The projection factor is divided into three sub-concepts : (1) a geometrical effect, (2) the velocity gradient within the atmosphere, and (3) the relative motion of the "optical" pulsating photosphere compared to the corresponding mass elements (hereafter fo-g). Both, (1) and (3) are deduced from geometrical and hydrodynamical models, respectively, while (2) is derived directly from observations. The Fe I 4896.439 A line is found to be the best one to use in the context of IBW and SB methods. A coherent and consistent period-projection factor relation (hereafter Pp relation) is derived for this specific spectral line: p = [0.064 +- 0.020] log P + [1.376 +- 0.023]. This procedure is then extended to derive dynamic projection factors for any spectral line of any Cepheid. This Pp relation is an important tool for removing bias in the calibration of the period-luminosity relation of Cepheids. Moreover, it reveals a new physical quantity fo-g to investigate in the near future.

We considered the structure of steady-state plane-parallel radiative shock waves propagating through the partially ionized hydrogen gas of temperature T_1 = 3000K and density 1e-12 gm/cm^3 <= \rho_1 <= 1e-9 gm/cm^3. The upstream Mach numbers range within 6 <= M_1 <= 14. In frequency intervals of hydrogen lines the radiation field was treated using the transfer equation in the frame of the observer for the moving medium, whereas the continuum radiation was calculated for the static medium. Doppler shifts in Balmer emission lines of the radiation flux emerging from the upstream boundary of the shock wave model were found to be roughly one-third of the shock wave velocity. The gas emitting the Balmer line radiation is located at the rear of the shock wave in the hydrogen recombination zone where the gas flow velocity in the frame of the observer is approximately one-half of the shock wave velocity. The ratio of the Doppler shift to the gas flow velocity of 0.7 results both from the small optical thickness of the shock wave in line frequencies and the anisotropy of the radiation field typical for the slab geometry. In the ambient gas with density of \rho_1 >= 1e-11 gm/cm^3 the flux in the H-alpha frequency interval reveals the double structure of the profile. A weaker H-beta profile doubling was found for \rho_1 >= 1e-10 gm/cm^3 and U_1 <= 50 km/s. The unshifted redward component of the double profile is due to photodeexcitation accompanying the rapid growth of collisional ionization in the narrow layer in front of the discontinuous jump.

This paper reviews the results obtained so far with the newly developed tomographic technique, which probes the velocity field across the atmosphere of long-period variable (LPV) stars. The method cross-correlates the optical spectrum with numerical masks constructed from synthetic spectra and probing layers of increasing depths. This technique reveals that the line doubling often observed in LPV stars around maximum light is the signature of the shock wave propagating in the atmosphere of these pulsating stars. Applications to the Mira variables RT Cyg and RY Cep are presented, revealing the outward propagation of the shock.

We considered the structure of steady-state radiative shock waves propagating in the partially ionized hydrogen gas with density rho1 = 1e-10 gm/cm^3 and temperature 3000K <= T1 <= 8000K. The radiative shock wave models with electron thermal conduction in the vicinity of the viscous jump are compared with pure radiative models. The threshold shock wave velocity above of which effects of electron thermal conduction become perceptible is of U1=70 km/s and corresponds to the upstream Mach numbers from M1= 6 at T1=8000K to M1=11 at T1=3000K. In shocks with efficient electron heat conduction more than a half of hydrogen atoms are ionized in the radiative precursor, whereas behind the viscous jump the hydrogen gas undergoes the full ionization. The existence of the electron conductive precursor leads to the enhancement of the Lyman continuum flux trapped in the surroundings of the discontinuous jump. For upstream velocities ranged within 70 km/s <= U1 <= 85 km/s the partially ionized hydrogen gas of the radiative precursor undergoes the additional ionization (<= 5%), whereas the total radiave flux emerging from the shock wave increases by 10% <= delta(FRad) <= 25% .

A tomographic method is described that makes it possible to follow the propagation of shock waves across the photosphere of long-period variable stars. The method relies on the correlation of the observed spectrum with numerical masks probing layers of different atmospheric depths. The formation depth of spectral lines is derived from synthetic spectra of non-variable red giant stars. When applied to Mira stars around maximum light, the tomographic masks reveal that the deepest photospheric layers are generally characterized by blueshifted absorption lines (attesting their upward motion), whereas the uppermost layers generally exhibit redshifted absorption lines (indicating their infalling motion). Double absorption lines are found in intermediate layers, where the shock front is located. At later phases, the shock front is seen moving towards upper layers, until it leaves the photosphere.

This paper presents statistics of the line-doubling phenomenon in a sample of 81 long-period variable (LPV) stars of various periods, spectral types and brightness ranges. When correlated with a mask mimicking a K0III spectrum, 54% of the sample stars clearly showed a double-peaked cross-correlation profile around maximum light, reflecting double absorption lines. Several pieces of evidence are presented that point towards the double absorption lines as being caused by the propagation of a shock wave through the photosphere. The observation of the Balmer lines appearing in emission around maximum light in these stars corroborates the presence of a shock wave. The observed velocity discontinuities, ranging between 10 and 25 km/s, are not correlated with the brightness ranges. A comparison with the center-of-mass (COM) velocity obtained from submm CO lines originating in the circumstellar envelope reveals that the median velocity between the red and blue peaks is blueshifted with respect to the COM velocity, as expected if the shock moves upwards.

The grid of the models of radiative shock waves propagating through partially ionized hydrogen gas with temperature 3000K <= T_1 <= 8000K and density 10^{-12} gm/cm^3 <= \rho_1 <= 10^{-9}gm/cm^3 is computed for shock velocities 20 km/s <= U_1 <= 90 km/s. The fraction of the total energy of the shock wave irreversibly lost due to radiation flux ranges from 0.3 to 0.8 for 20 km/s <= U_1 <= 70 km/s. The postshock gas is compressed mostly due to radiative cooling in the hydrogen recombination zone and final compression ratios are within 1 <\rho_N/\rho_1 \lesssim 10^2, depending mostly on the shock velocity U_1. The preshock gas temperature affects the shock wave structure due to the equilibrium ionization of the unperturbed hydrogen gas, since the rates of postshock relaxation processes are very sensitive to the number density of hydrogen ions ahead the discontinuous jump. Both the increase of the preshock gas temperature and the decrease of the preshock gas density lead to lower postshock compression ratios. The width of the shock wave decreases with increasing upstream velocity while the postshock gas is still partially ionized and increases as soon as the hydrogen is fully ionized. All shock wave models exhibit stronger upstream radiation flux emerging from the preshock outer boundary in comparison with downstream radiation flux emerging in the opposite direction from the postshock outer boundary. The difference between these fluxes depends on the shock velocity and ranges from 1% to 16% for 20 km/s <= U_1 <= 60 km/s. The monochromatic radiation flux transported in hydrogen lines significantly exceeds the flux of the background continuum and all shock wave models demonstrate the hydrogen lines in emission.

This paper is the first one in a series devoted to the study of the dynamics of the atmospheres of long-period variable stars. Results from a two-month-long monitoring of the Mira variables RT Cyg and X Oph around maximum light with the ELODIE spectrograph at the Haute-Provence Observatory are presented. The monitoring covers phases 0.80 to 1.16 for RT Cyg and phases 0.83 to 1.04 for X Oph. The cross-correlation profile of the spectrum of RT Cyg with a K0 III mask confirms that the absorption lines of RT Cyg in the optical domain appear double around maximum light. No line doubling was found in the optical spectrum of X Oph around maximum light, indicating that this feature is not common to all long-period variables. This paper also presents the application to RT Cyg of a new tomographic technique deriving the velocity field across the atmosphere by cross-correlating the optical spectrum with numerical masks constructed from synthetic spectra and probing layers of increasing depths. This technique reveals that both the temporal evolution of the line doubling, and its variation with depth in the atmosphere of RT Cyg, are consistent with the ``Schwarzschild scenario''. This scenario relates the temporal evolution of the red and blue peaks of the double absorption lines to the progression of a shock wave in the atmosphere. The temporal evolution of the Balmer Halpha, H beta, Hgamma and Hdelta emission lines around maximum light is also presented for RT Cyg and X Oph. The velocity variations of Halpha and of the absorption lines are discussed in the framework of two competing models for the formation of Balmer emission lines in long-period variable stars.

Models of steady-state plane-parallel shock waves propagating through the unperturbed hydrogen gas of temperature T=6000K and density rho = 1e-10 gm/cm^3 are computed for upstream velocities from 15 km/s to 70 km/s. The shock wave structure is considered in terms of the self-consistent solution of the radiation transfer, fluid dynamics and rate equations for 2 <= L <= 4 atomic bound levels with a continuum. The radiative flux Frad emergent from the shock wave was found to be independent of the lower limit nu_L of the frequency range provided that nu_L < nu_2, where nu_2 is the Balmer continuum head frequency, At the same time the decrease of nu_L is accompanied by decrease of the Lyman continuum flux and leads to smaller heating and weaker ionization of the hydrogen gas in the radiative precursor. For all models the size of the radiative precursor is of 10^4 cm and corresponds to several mean free paths of photons at the frequency of the Lyman continuum edge nu_1. The compression ratio at the discontinuous jump gradually increases with increasing upstream velocity U_1, reaches the maximum of \rho^+/\rho^- = 3.62 at U_1 = 55 km/s and slowly decreases for larger U_1 due to the strong rize of the preshock gas temperature. The radiative flux from the shock wave was determined as a function of the upstream velocity U_1 and its ratio to the total energy flux in the shock wave C_2 was found to range within 0.18 < Frad/C_2 < 0.92 for 15 km/s <= U_1 <= 65 km/s. Thus, at upstream velocities U_1 > 60 km/s the shock wave losses more than 90% of its total energy due to radiation. For all shock wave models the role of collisional processes in both bound-bound and bound-free atomic transitions was found to be negligible in comparison with corresponding radiative processes.

The structure of steady plane-parallel radiative shock waves propagating through the hydrogen gas undergoing partial ionization and excitation of bound atomic states is investigated in terms of the self-consistent solution of the equations of fluid dynamics, radiation transfer and atomic kinetics. The shock wave model is represented by a flat finite slab with no incoming radiation from external sources at both its boundaries. The self-consistent solution is obtained using the global iteration procedure each step of which involves (1) integration of the fluid dynamics and rate equations for the preshock and postshock regions, consecutively, both solutions being fitted by the Rankine-Hugoniot relations at the discontinuous jump; (2) solution of the radiation transfer equation for the whole slab. The global iteration procedure is shown to converge to the stable solution which allows for the strong coupling of the gas flow and the radiation field produced by this flow. Application of the method is demonstrated for the shock waves with upstream velocities from 15 km/s to 60 km/s (that is with upstream Mach numbers from 2.3 to 9.3) and the hydrogen gas of unperturbed temperature T=3000K and density rho = 1e-10 gm/cm^3.