The GPS system is employed now as a very important element of the global information infrastructure e.g.  . Free and reliable access to this system led to the development of hundreds of applications that affects every aspect of modern life. GPS technology is now used almost everywhere, from precise clocks and cell phones, to bulldozers, sea vessels, military applications and bank computers. It increases productivity in various sectors of the economy, including agriculture, aviation, building, defense systems, all branches of Earth’s sciences, postal services, etc. The use of the GPS system saves lives, prevents traffic accidents, and promotes search and rescue operations, speeds up the delivery of emergency services and the elimination of the consequences of natural disasters. This system is vital for the next generation of the automatic transportation  , which should help to improve vehicle safety with increasing reliability and capacity. GPS is also of high importance for the majority of scientific, environmental, technical and medical tasks: weather (climate) forecasting e.g.  , earthquake monitoring, searching economic minerals, environmental protection, archaeological prospection, and safe driving  -  . The GPS system (space sensor net) is also crucial for any national security.
2. Problem Formulation
The theory of radiation, propagation and absorption of light originated in the beginning of the 21th century and is described in detail in the courses of common and theoretical physics e.g.  . In our problem, we will talk about the features of the propagation of photons radiated by satellite through a resonance quantum medium in the D and E layers of the upper atmosphere of the Earth. Suppose that the transmitter signal is a certain group N0 of photons sent in time t0 in a certain order. The energy of a photon (here and in the future, we will use the atomic system of units ) is defined as , where is the carrier frequency of the transmitter, which is assumed to be given. The total energy of the group will be , where N0 >> 1. For simplicity, we write the initial radiation power of the group (without taking into account the details of the emission of photons) as . Then the power of the radio emission flux has the form
where S is the area of the flux.
When a satellite signal propagates through the D and E layers of the atmosphere, photons (moving in its electromagnetic field) are scattered on the states of orbitally degenerate Rydberg complexes A**M  , which are split by inelastic interaction with the neutral medium. The symbol A** denotes the highly excited states of atmospheric molecules N2, O2 and NO, and M includes quadrupole diatomic molecules of nitrogen and oxygen. The total wave functions of these states is superposition over the angular momenta l of the weakly bound electron (for ) and are characterized by the angular momentum of the electron L with respect to the M molecule  . The values of the levels shift of the complex split off from the degenerate Coulomb levels due to the interaction of A** with the molecule of the medium M are noticeably greater (by an order of magnitude n2, where n is the principal quantum number) than the shifts of the isolated Rydberg levels (for )
where a is the electron scattering length on the molecule M.
Under these conditions, the influence of the neutral medium reduces to an irreversible process l of mixing with formation of the long-lived orbitally degenerate complexes A**M that are populated in a two-temperature plasma at the height of 80 - 110 km over the Earth  . At the lower boundary of this layer (domain), the Rydberg states are effectively quenched. The upper boundary of this layer corresponds to the restriction on the density of the medium for the mixing process l  . The characteristic frequency of transitions in the complex A**M (occurring without changing the principal quantum number n) in this region is of the order of 1 GHz (or 10−5 eV).
The strength of the satellite signal field is determined by its power and, according to Equation (1), is proportional to the number of photons . Since the wavelengths of the satellite GPS signal belong to the decimeter range and are much larger than the size of the complex, all the transitions can be considered using the dipole approximation. These transitions represent an electromagnetic interaction operator with the complex in the form of , where f is the quantity of the field strength operator determined by the source power, and D is the dipole moment operator of the complex A**M  .
The magnitude of the propagation time delay is described by a simple formula proposed in  :
where nis the radiofrequency emitted by the satellite (in Hz), is the ionospheric total electrons content (TEC) in the signal propagation path, along the column cross-section for one m2.
This model  is based on the assumption that the Earth’s ionosphere is a dispersive plasma medium that, strictly speaking, is valid only for the linear homogeneous media where harmonic fields independently exist. In essence, this is the simplest model of ionospheric plasma where the dispersion of the medium is associated with variations of the electron density at a fixed frequency. In real situations, we have more complex model of this phenomenon   .
The main goals of this investigation are following:
1) to establish the relationship between the observed features of the propagation of satellite signals through the atmosphere with the quantum dynamics of processes occurring within the orbitally degenerate Rydberg complexes,
2) to estimate the characteristic values of the shift and time delay of the signal in one act of interaction with the complex, which was not previously reported in the literature was considered.
It is important to note that Equation (3) is empirical and it turns out to be completely unsuitable for an increase in solar activity and occurrence of magnetic storms and disturbances   .
3. Types of Transition in Interaction of Photon with the Rydberg Complex
The description of the orbitally degenerate states of Rydberg complexes A**M requires, at a first stage, the introduction of diabatic vibronic potential energy surfaces (PES) , where the vector R is the distance between the diatomic molecules A** and M, u is their relative velocity, and the vector xM specifies the orientation of the molecule M in the coordinate system associated with the molecule A**.
Under vibronic we mean the electron-vibrational states of the Rydberg complex with the principal quantum number nv (depending on the vibrational quantum number v of the molecular ion A+) and the moment of the weakly bounded electron L relative to the molecule M. A regular method for construction of such calculation scheme is suggested in   . A rotation of the molecule M can be considered classically and we assume that the direction of its axis to be fixed, averaging over the orientations at the last stage of the calculations    .
In addition, we should average here the velocities, confining ourselves to the Maxwell distribution, since the main contribution is done by the domain of relative energy that is less or equal to the temperature of the neutral medium. Since the radiation processes in this formulation of the problem are described as quasi-static one, the frequency of radiation (or absorption) transitions is determined by the relation
Therefore, all subsequent averaging over orientations and relative velocities should be performed for a given transition with n frequency.
To understand the physics of the phenomenon, it suffices to confine ourselves to three types of transitions. The first concerns the intrinsic spontaneous microwave radiation of a neutral quantum medium, which is incoherent and propagates in an arbitrary direction. The dissociation of Rydberg complexes A**M is not taken into account in this case, since the distribution of populations of Rydberg complexes with respect to energy levels in a nonequilibrium plasma is assumed to be stationary and is maintained by a constant flux of free electrons  .
The second is a direct transition to the underlying state of the complex and corresponds to stimulated emission at the carrier frequency (Svanberg, 2001). It preserves the polarization of the photons of the satellite signal and can also be considered without taking into account a pre-dissociation. This leads to the creation of additional photons and the experimentally observed increase in the power of the received GPS signal. This effect was first discovered in the Afraimovich’s unique experiment  , which was held at the Cornell University (USA) and has not been explained until the present time. Direct measurements have shown that the power of the receiver’s received radiation at the frequency increases three-fold compared to the transmitter power on the satellite and is accompanied by signal failures at the receiver.
The third type of transition corresponds to the process of resonance scattering of a photon on a complex A**M, which proceeds through three successive stages. The first corresponds to stimulated absorption of a photon with transition to an intermediate overlying is the vibronic state of the complex. The next stage is connected with the quantum dynamics within the complex itself A**M, where along with the discrete states, the energy spectrum of the composite system contains autoionization states of the A** molecule. This is realized when the relative energy of the diatomic molecules A** and M in the complex is greater than the binding energy of the electron. Physically it corresponds to the case when , where Ta is the temperature of the neutral medium. A similar situation was considered earlier in  .
The process ends by the subsequent forced transition to the final state . We note that the general scheme of resonance scattering of photons on highly excited Rydberg molecules significantly differs from the traditional scheme of the processes described in  , which in principle cannot lead to a time delay. We also note that in the problem of resonance scattering of a Rydberg electron, the presence of predissociation channels should be taken into account when considering the dissociation of complexes for various fragments of a neutral medium independently  . Thus, the theory of propagation and distortion of satellite signals should include simultaneously the contributions of all three mentioned types of the radiation processes.
For a qualitative explanation of the features of GPS signal distortion during the propagation through a 60 - 110 km layer, it suffices to confine ourselves to the case of resonance (Rayleigh) elastic scattering (when the quantum state of the Rydberg complex A**M does not change) and analyze the resonant structure of the transition oscillator strength. The problem of the radiative transition to the ionization continuum of the Rydberg molecule (bound-free transitions) was considered in  . In this work, the total wave functions of the molecule are expressed in terms of the collision T-matrix elements describing the scattering of a slow electron on the ion core in the Coulomb potential and the force field of the M molecule  . Consequently, the dipole matrix element Dif becomes a superposition of the direct and two-stage (resonant) processes, and the strength of the transition oscillator acquires a typical form of the Fano contour.
Since in our case the initial and final electronic states of the Rydberg complex A**M are related, the frequency shift of the scattered photon and the delay time are formed in the intermediate second stage of the resonance process. This is due to the presence of a free term in the collision in the T-matrix, leading to a shift and broadening of the resonant lines  . They are determined by the multichannel electrostatic interaction of the Rydberg electron with the ion core and the configuration interaction with the dissociative continuum. Moreover, for the levels related to the vibrational excited states of the A+ ion (located near the corresponding levels of the series v=0), the connection with the dissociative continuum increases substantially  due to the decrease of l and the principal quantum number nv. Behavior of the mentioned factors is of high importance since the total wave function includes a superposition of all possible vibronic states with a given energy that also contribute to the shift of the photon frequency.
For a qualitative explanation of the propagation and distortion of the GPS signal, it suffices to limiting the case when the resonant (Rayleigh) elastic scattering (when the quantum state of the Rydberg complex A**M does not change) and analyze the resonant structure of the oscillator absorption strength of the photon .
Consequently, the dipole matrix element , becomes a superposition of the direct and two-stage (resonant) processes acquiring a typical form of the Fano contour. For pre-dissociation states, this situation arises with the participation of discrete vibronic levels of the vibrationally excited states of the A+ ion (located near the corresponding levels of the series ).It is caused by substantial increase of connection between the Rydberg configuration and the dissociative continuum due to the decrease of land the principal quantum number , since the total wave function includes the superposition all possible vibronic states with a given energy. As a result, the strength of the strength of the transition oscillator can be represented in the form
where is the strength of the direct transition oscillator, is the profile index, is the reduced shift of the carrier frequency, and is the width of the resonant level.
It is seen that the frequency shift of the scattered photon and the delay time in the intermediate state (due to the broadening of the autoionization of discrete levels) are completely determined by the strong nonadiabatic coupling of the electron and nuclear motions in the intermediate complex A**M and the configurational interaction with the dissociative continuum. The corresponding complex dipole transitions are formed by an entire set of principal quantum numbers n (Figure 1).
Curves shown in Figure 1 represent the frequency dependence of the radiation v on the principal quantum number calculated at the R = n distance. Lines begin from the position, which corresponds to a strong magnetic storm (with the electron density and temperature ). Over the
Figure 1. Lines of radiation (absorption) of quasimolecules A**N2 and A**O2.
frequency range up to 10 GHz, the distribution of emission lines depending on n and corresponding to contains four series of transition lines for each molecule. These series converge as L’ increases to the transition limit. ForN2 and O2 molecules, these lines are denoted by (marked in green) and (marked in blue), respectively. The frequency shift of and series limits with respect to each other result in a substantial nonuniformity of the resulting emission spectrum picture.
The selection rule here has the form of ( , where k is an integer), which is due to the lack of central symmetry for Rydberg molecules in the A** + M system  .
Preliminary estimates for single-resonant photon scattering give the following values for the shift of the time delay frequency Dt: ~ 10−2 GHz or ~ 10−10 s, respectively. Consequently, the formation of observable values of the delay time when a satellite signal passes through the atmospheric layer of 60 - 110 km should occur of no less than 106 scattering events of a group of photons. This agrees well with the previous calculations  , since the total number of Rydberg complexes in a column with an area per receiver antenna of 10 cm2, under normal conditions of the order of 107 - 108. Taking into account that the propagation velocity of a signal is much higher than the characteristic velocities of molecules inside the complexes, it sufficient to confine oneself to their quasi-static description. Then, taking the normal distribution for the Brownian motion, and taking into account that the probability of one event of a resonance scattering of a photon is , for 106 interactions we have Dt~ 10−7 s, i.e. about 100 nanoseconds.
Note that, up to now, the empirical Equation (3) is embedded in the GNSS positioning satellite systems for the delay time estimation. However, this expression was obtained as a result of statistical processing of the received signals under the normal atmospheric (physical-chemical) conditions. The situation changes significantly in the occurrence of magnetic storms (disturbances) when the concentration and flow of electrons increase noticeably. Attempts to describe such situations within the framework of classical electrodynamics in principle cannot lead to a delay, since it is well known and confirmed experimentally that the refractive index in the ionosphere (including plasma inhomogeneities) differs from “1.0” only in the fourth decimal place  .
We suggest that the resonant multiple interaction of the GPS signal with a quantum medium in atmosphere determines the delay process: this idea is firstly proposed in the scientific literature. The development of a general theory of signal propagation and signal delay should be based on an integral version of the theory of a multichannel quantum defect (MQD) and may be divided into three stages. At the first stage, it is assumed to construct a complete adiabatic wave function of the initial state of the Rydberg complex including the vibronic states and the exchange interaction (for the case where the Rydberg molecule and the neutral molecule of the medium coincide). To solve this problem, one can use the method developed in  for a simple system “Rydberg atom + atom with a filled electron shell”, and to modify it substantially.
At the second stage, it is necessary to develop a regular method for constructing the complete wave function of the complex taking into account the dynamics of its fragments, when the kinetic energy of their relative motion is included in the general analysis. This model also was not earlier discussed, and it proposed here for the first time. This will allow us to determine the corresponding dipole matrix elements of resonant scattering and to introduce into the theory at the third stage the power of the radiation obtained into receiver and the total delay time of the signal as a sum of the rescattering cascade along the propagation path of the signal. Simultaneously, it is necessary to develop a numerical program and calculate the total delay time of the satellite signal, depending on the main parameters of the nonequilibrium plasma (concentration, flux, and electron temperature) for various geomagnetic conditions in the ionosphere. We should note that a consistent quantum theory of the propagation of satellite signals for development of more precise model of the phenomenon must also take into account the Raman scattering of photons and the possibility of decay of A**M complexes, which have not been analyzed in the literature. This will allow calculating the reliable values of the carrier frequency shift and to compute the delay time of the incoming radio waves.
This work was supported by the Russian Foundation for Basic Research (Grant 16-05-00052).