Abstract
Ring laser gyroscopes are top sensitivity inertial sensors used in the measurement of angular rotation. It is well known that the response of such remarkable instruments can in principle access the very low frequency band, but the occurrence of nonlinear effects in the laser dynamics imposes severe limitations in terms of sensitivity and stability. We report here general relationships aimed at evaluating corrections able to effectively account for nonlinear laser dynamics. The so-derived corrections are applied to analyse thirty days of continuous operation of the large area ring laser gyroscope GINGERINO leading to duly reconstruct the Sagnac frequency \(\omega _s\). The analysis shows that the evaluated corrections affect the measurement of the Earth rotation rate \(\varOmega _{\oplus }\) at the level of 1 part in \(1.5\times 10^{3}\). The null shift term \(\omega _{ns}\) plays a non negligible role. It turns out proportional to the optical losses \(\mu \) of the ring cavity, which are changing in time at the level of \(10\%\) within the considered period of thirty days. The Allan deviation of estimated \(\varOmega _{\oplus }\) shows a remarkable long term stability, leading to a sensitivity better than \(10^{-10}\) rad/s with more than 10 s of integration time, and approaching \((8.5\pm 0.5)\times 10^{-12}\) rad/s with \(4.5\times 10^{5}\) s of integration time.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Ring laser gyroscopes (RLGs) are inertial sensors based on the Sagnac effect [1,2,3]. They are largely used for inertial navigation, and applications in geodesy, geophysics and even for General Relativity, where tests are foreseen [4]. Since 2011 we are studying the feasibility of the test of Lense–Thirring dragging of the rotating Earth at the level of \(1\%\) with an array of large frame RLGs [5,6,7]. For that purpose it is necessary to push the relative accuracy of the Earth rotation rate \(\varOmega _{\oplus }\) measurement in the range from 1 part in \(10^9\) up to 1 part in \(10^{12}\).
RLG consists of a laser with a cavity comprising of three or four mirrors, rigidly attached to a frame; large frame RLGs are utilised to measure the Earth rotation rate, being attached to the Earth crust. Because of the Sagnac effect, the two counter-propagating cavity beams have slightly different frequency, and the beat note of the two beams is proportional to the angular rotation rate of the ring cavity. Large frame RLGs are the most sensitive instruments for inertial angular rotation measurements. The Sagnac frequency of a RLG is in fact proportional to the component of the angular velocity \(\overrightarrow{\varOmega }\) felt by the instrument along the normal to the cavity plane:
where A is the area of the ring cavity, L is its perimeter, \(\lambda \) the wavelength of the light, and \(\theta \) is the angle between the area vector of the ring and \(\overrightarrow{\varOmega }\). For RLGs lying horizontally (area vector vertical) \(\theta \) is the co-latitude angle, while for RLGs aligned at the maximum Sagnac signal \(\theta = 0\). Equation 1 defines the scale factor S, which is a function of the geometry and of \(\lambda \), quantities than can be measured with a very high accuracy.
Further to sensitivity, other key points of such instruments rely on their broad bandwidth, which can span from
The laser dynamics is non linear and plays a role in determining the RLG signal. We have recently developed a model to reconstruct the Sagnac frequency \(\omega _s\) starting from the measured beat note \(\omega _m\) and using the mono-beam signals.Footnote 1
The Sagnac angular frequency can be expressed as the linear sum of several terms. As we have discussed in a recent paper [9], reconstruction of the Sagnac angular frequency is mostly affected by the correction \(\omega _{s0}\) associated with the so called back-scattering. In this paper we provide the necessary information to evaluate the other correction terms. Thirty days of continuous operation of GINGERINO are analysed taking into account the set of first order expansion terms, and results discussed.
Results demonstrate that the so called null shift term, the first order correction \(\omega _{ns1}\), plays a non negligible role in affecting both accuracy and sensitivity of the apparatus. The finding is in contrast with conventional treatments of RLG data, where null shift effects are typically neglected or considered simply proportional to the difference of the two mono-beam signals. On the contrary we found that including such correction leads to improve the measurement accuracy. Furthermore, since null shift effects can be related to the optical losses of the system, which inherently depend on operation time, a time dependent effect is found.
2 The analysis scheme to take into account laser dynamics
Our approach is based on the model of laser dynamics we recently developed in details in [9] based on the Aranowitz model of RLG [10, 11], where laser dynamics is described in terms of several dimensionless parameters (Lamb parameters). The polarization of the laser plasma is described up to the third order expansion in powers of the field [10, 11]. In the case of large frame RLGs, to avoid mode competition, a 50:50 mixture of two isotopes, \(^{20}Ne\) and \(^{22}Ne\), is utilised, and the laser operation is set close to threshold to guarantee single-mode operation. This particular choice allows for some simplifications: the Lamb parameters of cross-saturation (\(\theta _{12}\) and \(\theta _{21}\)) can be neglected, and we can assume the self-saturation terms equal to each other (i.e., \(\beta _1 = \beta _2 = \beta \)).
The general goal of the analysis is to evaluate the Sagnac frequency correcting systematics due to the laser nonlinear dynamics by using the available data: the measured beat note and the mono-beam signals [9, 12]. Figure 1 shows a scheme of the experimental apparatus, where photodiodes collecting the relevant signals are indicated.
The model provides two analytical expressions for the Sagnac angular frequency \(\omega _s = 2 \pi f_s\): a non-linear relationship (Eq. 8 of Ref. [9]). This analytical expression was expanded up to the second order in term of two parameters \(\delta _{ns}\) and K, and can be calculated from the Lamb parameters. Accordingly, \(\omega _s\) is expressed as the linear sum of a total of six terms:
As demonstrated in [9] \(\omega _{s0}\), accounting for back-scatter effects, is the dominant term in the expansion. The term \(\omega _{s0}\) does not depend on Lamb parameters and can be expressed starting from the knowledge of experimental parameters according to:
where \(\omega _m\) is the measured beat note (expressed as angular frequency), \(I_{1,2}\), \(I_{S1,S2}\), and \(\epsilon \) are the DC level of the mono-beam signals, their amplitude demodulated at the beat note frequency, and the dephasing between the two signals, respectively.
In order to better define \(\omega _{s0}\) and take into account experimental issues, such as dark currents in the photodiode signals, the term \(\omega _{s\xi }\), see Eq. 4, has been introduced, where \(\xi \) is a proportionality constant which could be evaluated by statistical means.
The other terms in Eq. 2 give smaller contributions. They can be evaluated as detailed in the Appendix, starting from the knowledge of the Lamb parameters [10], a careful estimation of all operating parameters of the system, and the evaluation of the total optical losses \(\mu \). The latter can be measured with a \(1\%\) accuracy via the ring down method [11], i.e., by suddenly switching off the laser discharge and measuring the time behaviour of the intracavity signals through photodiodes PH1 and PH2 in Fig. 1. The method cannot be applied continuously during the RLG operation, since it requires to switch off the laser. Moreover, optical losses can change during operation, therefore their evaluation should be repeated frequently.
Because of that, and also to avoid uncertainties and fluctuations related to estimation of operating parameters such as, laser beam size, pressure and polarization conditions of the discharge (see Appendix for their use in the relevant expressions), in the present paper we have adopted a statistical approach for evaluating \(\mu \). As shown in the next subsection, this parameter enters as a linear proportionality factor into \(\omega _{ns1}\). Therefore, we can express \(\omega _{ns1} = \mu \times {\tilde{\omega }}_{ns1}\), where \({\tilde{\omega }}_{ns1}\) is evaluated starting from known data and \(\mu \) determined through a linear regression procedure. The terms in Eq. 2 containing the subscript K require some further discussion. As reported in the appendix, the parameter K is proportional to oscillating terms at the frequency \(\omega _s\). The evaluation of the K terms is not straightforward, and several approaches are in principle feasible. In the current analysis, we assume that \(\omega _s\) is nearly constant in time, and as a consequence for frequencies much below \(\omega _s\) the contribution of these terms should average to zero. As reported in the appendix, assuming \(\langle \sin {\omega _s t} \rangle = \langle \cos {\omega _s} t \rangle = MK\) the equations is simplified, emerging a proportionality term MK, that is evaluated via the statistical approach.
In summary, the analysis procedure is as follows: in the first level, different terms \({{\tilde{\omega }}}_{ns1,2}\) are evaluated and stored on disk; in the second level, updated \(\xi \), \(\mu \) and MK coefficients are found with linear regression models, minimising the sum of Eq. 2. It is important to remark that \(\mu \) can be stable with time, or varying with time. In this second case it is necessary to add to the model the information necessary to follow the time behaviour. Section 2.3 presents the method used for GINGERINO. We remind that the available signals from the RLG are sampled at high frequency (typically 5
In the following, 30 days of data produced by GINGERINO will be analysed. In the analysis, corrections are included only up to the first order: \(\omega _{ns1}=\omega _{ns}\), \(\omega _{K1}=\omega _{K}\). Second order terms are in fact smaller and do not contribute to improve the results compared to the present state of accuracy. The Appendix contains however all information to evaluate the different terms in the general case.
2.1 Some details on the evaluation of first order expansion terms
As already mentioned, we start from the DC components, \(I_1\) and \(I_2\), of the mono-beam signals, their AC components at the beat frequency \(\omega _m /2 \pi \), \(I_{S1}\) and \(I_{S2}\), and the relative phase \(\epsilon \) between the two AC components.
Dissipative processes (like diffusion or absorption from the mirrors) can produce non reciprocal losses between the two counter-propagating beams, so that two distinct loss parameters \(\mu _{1} \ne \mu _{2}\) must be used. Without loss of generality, it is possible to take \(\mu _1=\mu \) and \(\mu _2=\mu +\delta \mu \), where \(\mu \) represents the reciprocal losses term and \(\delta \mu \) the non reciprocal ones. The plasma dispersion function depends on several parameters: gas pressure, line width of the excited isotopes, resonance frequency of the two counter propagating beams, area of the beam profile at the discharge, and temperature of the plasma. Assuming that \(\beta _1 = \beta _2 = \beta \) (a rough evaluation reported in [13] gives \(\beta _1 - \beta _2 \simeq 10^{-14}\)), combining the information of the mono-beam DC signals \(I_1\) and \(I_2\) with the plasma dispersion function, the gains of the two beams Gain1 and Gain2 (see Appendix) are evaluated. Since \(\beta _1 = \beta _2\) the non reciprocal loss term \(\delta \mu \) can be hence determined. Using other relationships reported in the literature [14, 15] it is possible to write \(\omega _{ns1,2}\), \(\omega _{K1,2}\) and \(\omega _{nsK}\) as a function of the known quantities \(\omega _m\), \(I_1\), \(I_2\), \(I_{S1}\), \(I_{S2}\), \(\epsilon \), and \(\mu \). As shown in the Appendix, losses appear as a multiplicative factor in \(\omega _{ns1,2}\).
Further to be dependent on losses, terms correcting the null shift depend also on \(\epsilon \), as well as on geometry and size of the cavity, in a rather complicated way. Therefore, minimization of \(\omega _{ns1,2}\) is not straightforward. To give an example using the parameters of GINGERINO, the minimum is found for \(\delta \mu = 0\); assuming typical values in physical units \(I_1=I_2=1\) V and \(I_{S1}=I_{S2}=0.01\) V (see Appendix for the calculation), the null shift correction is of the order of 0.4 \(\mu \)Hz, with a modulation induced by variations of \(\epsilon \) of the order of \(5\times 10^{-6}\) \(\upmu \)Hz. In general \(\delta \mu \ne 0 \), as can be inferred from the difference in mono-beam signals. While keeping all other parameters fixed as in the example above, a difference in the demodulated photodiode signals \(I_{S1}\) and \(I_{S2}\) of 1 mV leads to remarkable values of the null shift correction. In such conditions, the null shift contributes to around 500 \(\upmu \)Hz, with modulations induced by \(\epsilon \) of the order of 0.013 \(\upmu \)Hz.
The above examples suggest that, in general, variations with \(\epsilon \) are rather small and can be neglected. Therefore, in the ideal case of equal losses the effect of the null shift is a systematic term affecting the accuracy of the measurement at the level of \(1.5\times 10^{-9}\), assuming variations of the losses with time below \(10\%\), while when \(\delta \mu \) is not negligible the effect is orders of magnitude higher. In the considered case the requirement on the stability of \(\mu \) with time is of the order of 1 part in \(10^4\) in order to achieve a precision of 1 part in \(10^9\). As far as accuracy is concerned, the systematic effect is of the order of 1 part in \(10^6\), one thousand times the goal of GINGER (1 part in \(10^9\) at least).
While being less relevant than those related to null shift effects, also the terms with the K suffix deserve discussion. Equation 2 shows that such terms are linearly summed to \(\omega _s\). As a consequence, their contribution to low frequency signals, i.e., those stimulating the major interest in typical applications of RLGs, can be reduced by low pass filtering. Despite of the possibility to remove related contributions through low pass filtering, the evaluation of \(\omega _K\) requires some further assumptions. \(\omega _K\) oscillates at the Sagnac frequency, containing \(\sin (\omega _s t)\) and \(\cos (\omega _s t)\) terms. However, due to nonlinear dynamics, the long time average of such oscillating terms can deviate from zero. Therefore, we can express \(\omega _K\) as the product of a small parameter MK with a given function of known parameters. As already mentioned, the value of MK and its very slow variations can be found via linear regression.
A Mathematica notebook is reported in the Appendix to show details of the symbolic calculations leading to \(\omega _{ns}\) and \(\omega _K\) functions, to demonstrate how losses are a scale factor for the \(\omega _{ns}\) terms, and to provide numerical examples of the involved values calculated for the GINGERINO setup. We note that, according to our model, \(\omega _K\) can be written as the ratio of two polynomials where the parameter \(\mu \) enters in one term, preventing direct evaluation of such a parameter through a simple fit procedure. However, we have checked that, with typical \(\mu \) values of \(100-1\) ppm, variations in calculated results are not significant.
In the following, we focus onto calculations for a general RLG dynamics, and their implementation for the specific case of GINGERINO.
2.2 The selected playground and the first step of the analysis
Thirty days of continuous operation between June 16 and July 15 2018 have been selected for this test. After June 21, heavy operations took place in the underground laboratories where GINGERINO is placed for the building of another experiment. Disturbances induced in GINGERINO operation, consisting of mostly mechanical effects, are well visible in the second half of the data stream. We included also these data to show that, even in non ideal conditions, very high sensitivity and stability can be obtained. The data are acquired with our DAQ system at 5
2.3 The second step of the analysis
A second routine collects offline different days for longer analysis. The first operation is to identify the portions of data which are not at normal operation, typically less than 5% of the data are removed. To this aim, we make use of the fringe contrast of the interferogram acquired at 20 Hz rate. The choice of such a sampling rate minimizes the discarded sample. This operation associates the flag 1 to the good samples and the flag 0 to the bad ones; from this routine we get the mask to select data decimated at lower frequency.
In the present analysis, decimated at half an hour rate, above \(90\%\) of data are selected. Both the intensity of the mono-beam signal (\(I_2\), collected by the photodiode indicated as PH2 in Fig. 1), that is not used for the feedback control of the laser power, and of the Discharge Monitor (DM, i.e. a photodiode looking at the laser discharge) change with time. This is a clear indication that \(\mu \) changes with time; in particular, \(I_2\) indicates the presence of non reciprocal losses \(\delta \mu \), while DM is proportional to the total excited atoms, and therefore to the total losses. DM and \(I_2\) signals allow us to model the change of losses as a function of time through two form factors:
gain monitor form factor: \(FF_{DM} = \frac{DM- \langle DM \rangle }{ \langle DM \rangle }\)
mono-beam 2 form factor: \(FF_{PH2}=\frac{I_2- \langle I_2 \rangle }{\langle I_2 \rangle }\).
2.4 Reconstruction of the Sagnac frequency by a linear regression model
The vectors \(\omega _{ns}\), \(\omega _{ns}\times FF_{DM}\), \(\omega _{ns}\times FF_{PH2}\), \(\omega _{\xi }\), and \(\omega _{K}\) are collected, and the linear regression method is utilised to determine the unknown parameters \(\mu \), MK and \(\xi \), and evaluate \(\omega _s\) accordingly.
Let us remind that \(\omega _{s0}\) accounts for back-scatter noise, and an additional term \(\omega _\xi \) has been implemented to account for inaccuracies on the measured quantities \(I_{S1,2}\), \(I_{1,2}\) and \(\epsilon \) [9]. In summary, in the procedure \(\mu \) and MK are physical quantities which could be evaluated independently, while \(\xi \) compensates noises in the measured quantities.
Figure 4 shows data before (top panel) and after (bottom panel) the application of this procedure. The dispersion across the mean value is clearly suppressed using \(\omega _s\), meaning that the new regressive analysis, which includes the terms \(\omega _{ns}\) and \(\omega _K\), duly accounts for it. The Matlab function fitlm has been used to perform linear regression. Results indicate that the p-value, which tests the null hypothesis, is approximately zero for all vectors excluding the one related to \(\omega _K\), whose p-value is 0.038. On the other hand, the relative error in the estimation of the MK term involved in the \(\omega _K\) expression is around 50%, whereas it amounts to \(1-12\)% for all other vectors.Footnote 2
Figure 5 shows the contribution of the laser systematics. A dominant role is played by \(\omega _{ns}\), which produces a slowly variable level ranging, in units of frequency, between 180–200 mHz, with standard deviation 2.3 mHz. \(\omega _\xi \) gives a small correction with mean value \(-5.7\) mHz and standard deviation 2.3 mHz. These two terms affect residuals, whereas \(\omega _K\) is smaller by more than a factor 10; it has a mean value compatible with zero and standard deviation 0.05 mHz, see Fig. 6.
Figure 7 shows the evaluated losses. During the 30 day period considered, a variation of \(12\%\) is observed in \(\mu \), with a visible trend towards higher losses for increasing time. Figure 8 compares the distributions of the beat frequency \(\omega _m /2 \pi \), of the \(\omega _{s0} /2 \pi \) reconstructed in the first stage of our approach and that of the finally determined \(\omega _s / 2 \pi \). It should be noticed that the distribution of \(\omega _m\) is highly non-Gaussian while that of both reconstructed frequencies becomes close to Gaussian, indicating that the nonlinear terms of dynamics are correctly accounted for. Moreover, \(\omega _s\) is shifted towards lower frequencies as a consequence of the null shift term \(\omega _{ns}\).
Figure 9 reports the modified Allan deviation, expressed in angular rotation rate. We observe that, after the application of the procedure (red line), the variance is decreasing with the integration time up to more than 2 days. This suggests that with our procedure we are effectively correcting the long term laser dynamic effects. It is well known that RLGs are sensitive to Chandler Wobble effect, which is typically below 100 nrad. We have checked that the level reached in the long term stability is above the limit imposed by the Chandler Wobble effect. The so far obtained long term instability is about 10 times larger than the scale factor changes induced by temperature variations considering \(5\times 10^{-6}\)/oC the thermal expansion coefficient of granite, the material used for the RLG frame, and 0.02oC as typical temperature variation in the 30 day period. The obtained modified Allan deviation represents a fair improvement compared to the analysis in [12] (the improvement is around a factor 7) carried out by using available signals on a pure phenomenological basis, for instance the long time behaviour was cancelled using the DM signal.
3 Discussion and conclusions
The Sagnac frequency \(\omega _s\) of the GINGERINO RLG has been evaluated taking into account the laser dynamics. The model leads to a linear sum of several terms, which can be determined with the available signals of beat frequency and mono-beams. The so-called back-scatter noise, which has been so far considered the most severe limitation of RLG, is accounted for by terms \(\omega _{s0}\) and \(\omega _\xi \). It turns out that the null shift term \(\omega _{ns}\) is the second dominant contribution. It depends linearly on cavity losses \(\mu \); this term definitely affects accuracy of the apparatus, since in the best case, if the losses are constant, it represents a non-negligible DC contribution. Furthermore, it is likely that losses are not constant with time, hence \(\omega _{ns}\) affects stability. If a long term stability of 1 part in \(10^9\) is required for the \(\varOmega _{\oplus }\) measurement, the requirement for the stability of \(\mu \) with time, or equivalently for the relative accuracy in the evaluation of \(\mu \), is \(\frac{\omega _{ns}}{\omega _{s}}\times 10^9\). According to the model, it is straightforward to see that the null shift level is minimised when \(\delta \mu \simeq 0\).
Changes of \(\mu \) at the level of \(10\%\) are evident in GINGERINO. In order to follow those changes we have utilised the DC mono-beam signal of the uncontrolled beam, that is the signal which is not used in feedback loops controlling the laser operation (beam 2 in our case), and the discharge monitor DM, also called gain monitor. \(\omega _s\) is evaluated with a linear regression model.
Remarkably, during the measurements considered in this paper, in particular after the twentieth day of the 30 day record, heavy activity was present for the construction of a new experiment inside the underground laboratory hosting the RLG. Despite that, and the fact that GINGERINO is not equipped with control of the geometrical scale factor (accordingly its wavelength is not fixed), in 5.4 days of integration time it has reached the relative stability of \((1.7\pm 0.07)\times 10^{-7}\) [corresponding to a resolution in the angular velocity measurement of \((8.5 \pm 0.5)\times 10^{-12}\) rad/s].
The analysis poses the problem of the accuracy of angular velocity estimations; however, it is not straightforward to understand the accuracy of the final measurement, an issue to be addressed by future Monte Carlo studies to be carried out following the approach we have already employed in [13].
Main goal of the present analysis was to demonstrate that null shift effects can be evaluated in analytical terms and that they are not negligible at the present sensitivity level of our instrumentation. In order for the linear regression method to give reliable results, the model must be complete. Our data are affected by seismic contributions. They are obviously not considered in the model, but, due to their typical frequency, their role is suppressed thanks to low pass filtering used in our procedure. Furthermore, environmental conditions, in particular temperature, can affect RLG operation. We have verified that their role is negligible for the present data set. In order to check that the procedure does not artificially affect shape and properties of the input signals, we have added a spurious signal synthetically produced to the data, with amplitude of the order of fractions of nrad/s, and verified that it was not modified due to the procedure.
The analysis has also shown that losses are the main limitation of GINGERINO. The question is now which part of the apparatus is the main responsible for losses and for their behaviour as a function of time. The hypothesis is that the main contribution comes from the gain tube of the laser. In the standard high sensitivity RLG scheme, it acts both as gain tube and spatial filter. The laser discharge is in fact produced within a pyrex capillary, with inner diameter of 4 mm, which is a bit small compared to the waist of the beams.
In the near future, a new capillary tube with slightly larger inner diameter will be installed in GINGERINO to limit losses. Furthermore, additional data, such as current of the discharge tube and laser wavelength, will be continuously recorded in order to better characterise the instrument operation in the long term.
As far as the analysis is concerned, the effect of terms \(\theta _{1,2}\) will be eventually included in the model to account for possible deviations of the laser mixture behaviour from the expected one. Furthermore, the data analysis will be extended in order to evaluate the null shift second order terms \(\omega _{ns2}\) based on the recipes given in the present paper and to determine the associated correction in data from GINGERINO. Finally, we plan to devote efforts aimed to better understand the influence of each parameter, in particular of the change of wavelength, in the RLG operation, which for the present set of data matters, since GINGERINO is free running and the wavelength can change with even small temperature variations.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The data of this manuscript will not be deposited. The data can be made available from the corresponding author on reasonable request.]
Notes
For mono-beam signals we intend the signals observed by the two photodiodes that detect the laser intensities in the counter-propagating directions, see Fig. 1.
The coefficient of determination of the linear regression procedure and the standard deviation of the obtained residuals are \(R^2 = 0.952\) and \((770\pm 3)\) \(\mu \)Hz, respectively.
References
K.U. Schreiber, J.-P.R. Wells, Invited review article: large ring lasers for rotation sensing. Rev. Sci. Instrum. 84(4), 041101 (2013)
A. Gauguet (ed), The Sagnac effect: 100 years later. C. R. Phys. 15(10), iii – iv (2014)
A. Di Virgilio, M. Allegrini, A. Beghi, J. Belfi, N. Beverini, F. Bosi, B. Bouhadef, M. Calamai, G. Carelli, D. Cuccato, E. Maccioni, A. Ortolan, G. Passeggio, A. Porzio, M.L. Ruggiero, R. Santagata, A. Tartaglia, A ring lasers array for fundamental physics. C. R. Phys. 15(10), 866–874 (2014). The Sagnac effect: 100 years later / L’effet Sagnac : 100 ans après
M. Tercjak, A. Brzeziński, On the influence of known diurnal and subdiurnal signals in polar motion and ut1 on ring laser gyroscope observations. Pure Appl. Geophys. 174(7), 2719–2731 (2017)
F. Bosi, G. Cella, A. Di Virgilio, A. Ortolan, A. Porzio, S. Solimeno, M. Cerdonio, J.P. Zendri, M. Allegrini, J. Belfi, N. Beverini, B. Bouhadef, G. Carelli, I. Ferrante, E. Maccioni, R. Passaquieti, F. Stefani, M.L. Ruggiero, A. Tartaglia, K.U. Schreiber, A. Gebauer, J.-P.R. Wells, Measuring gravitomagnetic effects by a multi-ring-laser gyroscope. Phys. Rev. D 84, 122002 (2011)
A. Tartaglia, A. Di Virgilio, J. Belfi, N. Beverini, M. Luca Ruggiero, Testing general relativity by means of ring lasers. Eur. Phys. J. Plus 132(2), 73 (2017)
A.D.V. Di Virgilio, J. Belfi, W.-T. Ni, N. Beverini, G. Carelli, E. Maccioni, A. Porzio, GINGER: a feasibility study. Eur. Phys. J. Plus 132(4), 157 (2017)
A. Simonelli, H. Igel, J. Wassermann, J. Belfi, A. Di Virgilio, N. Beverini, G. De Luca, G. Saccorotti, Rotational motions from the 2016, central italy seismic sequence, as observed by an underground ring laser gyroscope. Geophys. J. Int. 214(1), 705–715 (2018)
A.D.V. Di Virgilio, N. Beverini, G. Carelli, D. Ciampini, F. Fuso, E. Maccioni, Analysis of ring laser gyroscopes including laser dynamics. Eur. Phys. J. C 79(7), 573 (2019)
F. Aronowitz, R.J. Collins, Lock-in and intensity-phase interaction in the ring laser. J. Appl. Phys. 41, 1 (1970)
D. Cuccato, A. Beghi, J. Belfi, N. Beverini, A. Ortolan, A. Di Virgilio, Controlling the nonlinear inter cavity dynamics of large he-ne laser gyroscopes. Metrologia 51, 97–107 (2014)
J. Belfi, N. Beverini, G. Carelli, A. Di Virgilio, U. Giacomelli, E. Maccioni, A. Simonelli, F. Stefani, G. Terreni, Analysis of 90 day operation of the GINGERINO gyroscope. Appl. Opt. 57(20), 5844–5851 (2018)
A. Beghi, J. Belfi, N. Beverini, B. Bouhadef, D. Cuccato, A. Di Virgilio, A. Ortolan, Compensation of the laser parameter fluctuations in large ring-laser gyros: a Kalman filter approach. Appl. Opt. 51(31), 7518–7528 (2012)
D. Cuccato, Modeling, estimation and control of ring laser gyroscopes for the accurate estimation of the Earth rotation. Ph.D. Thesis, University of Padova (2016)
Master thesis of Davide Cuccato. http://www.infn.it/thesis/thesis_dettaglio.php?tid=6109
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix A: Identification of He–Ne ring laser dynamics
Constants and parameters of He–Ne ring lasers
\(\text {Beam } - = \text {Cw}~, \text {Beam }+ = \text {Ccw} , \,\text {Isotope}\, 1 = \text {Ne20}, \text {Isotope} \,2 = \text {Ne22}\)
1.1 Electronic constants to express all in function of the monobeams voltage
1.2 Evaluation of output power and mean intensities from the above constants
1.3 Plasma dispersion function
1.4 Lamb parameters expressions for the single isotope case
1.5 Lamb parameters expressions for the double isotope case
20 and 22 are the isotopes; X and Y are the counter-propagating beams
1.6 Evaluating lamb parameters
evaluation of the corrective terms
Null Shift
1.7 K and \(\omega _K\) evaluations
mixed term
Approximated value of \(\omega _{\mathrm{k1}}\)
Assuming \(\langle \hbox {Cos}[\omega \hbox {s t}] \rangle = \langle \hbox {Sin}[\omega \hbox {s t}] \rangle = \hbox {MK}\), where MK is a proportionality constant estimated by statistics
1.8 Calculations of \(\alpha _{1,2}\) and \(r_{1,2}\) parameters
1.9 Operation parameters of GINGERINO
1.10 Calculation of the two functions used in the paper
In the following a few relations are explicitly reported using the parameters of GINGERINO
In this analysis \(\mu \) indicates the losses of mode 2, while \(\delta \mu \) is such that the losses of beam 1 is \(\mu +\delta \mu \)
\(\delta \mu \) is perfectly deducible: \(\delta \mu \) = DMU
\(\mu \) is a proportionality function for \(\omega _{\mathrm{ns1}}\)
\(\mu \) is a proportinality function for \(\delta _{\mathrm{ns}}\)
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Di Virgilio, A.D.V., Beverini, N., Carelli, G. et al. Identification and correction of Sagnac frequency variations: an implementation for the GINGERINO data analysis. Eur. Phys. J. C 80, 163 (2020). https://doi.org/10.1140/epjc/s10052-020-7659-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-7659-6