Method for determining uncertainty and error in the process of ophthalmic lens calibration

Método para la determinación de la incertidumbre y el error en el proceso de calibración de lentes oftálmicas

DOI 10.22430/22565337.1910 Table Figure

Received: March 9, 2021
Accepted: Augst 25, 2021
Available: Novemver 30, 2021

 

How to cite / Cómo citar
A. Salgar-Marín; J. A. Vargas; A. F. Ramirez-Barrera, “Method for determining uncertainty and error in the process of ophthalmic lens calibration”, TecnoLógicas, vol. 24, nro. 52, e1910, 2021. https://doi.org/10.22430/22565337.1910

 

Abstract

In the present investigation, a scientific procedure was developed, and a mathematical model was proposed, with the objective of determining, under standard conditions, the uncertainty, and the measurement of dioptric power in ophthalmic lenses. The methodology of the scientific procedure is based on the fundamentals of geometric optics, this process guarantees and establishes a standardized uncertainty measure in repeatable and reproducible processes. The methodology is complemented with a proposed mathematical model based on the guide for the expression of uncertainty in measurement - GUM. This model can be applied to lenses used for calibrating eye care equipment (such as lensometers, which are used to diagnose myopia and farsightedness) by evaluating the lenses without having direct contact with patients. When the proposed mathematical model was applied, its experimental result was a maximum expanded uncertainty of ± 0.0079 diopters in a 0.5-diopter lens. This is optimal compared to the result of other authors this article, who reported a maximum expanded uncertainty of ± 0.0086 diopters. In conclusion, the application of this scientific procedure provides manufacturers and users of this type of lenses with a reliable measurement thanks to a calibration process based on geometrical optics and centered on patient safety.


Highlights

  • The uncertainty estimation about the measurement of the dioptrical power impacts the confiability of the lens manufacturing

  • The calibration of the ophthalmic lenses impacts in the results of the diagnosis in optometry

  • The most important contribution of measurement uncertainty in the calibration of ophthalmic lenses is given by the deviation of the measurement.

  • The precision of the optical alignment contributes to the assurance of the measurement

  • Keywords: Optical metrology, Calibration function, Lens power, Focal length, Measurement uncertainty.

    Resumen

    En la presente investigación se desarrolló un procedimiento científico, y se propuso un modelo matemático, con el objetivo de determinar, bajo condiciones estándar, la incertidumbre y la medida de potencia dióptrica en lentes oftalmológicos. La metodología del procedimiento científico está basada en los fundamentos de la óptica geométrica, este proceso garantiza y establece una medida de incertidumbre estandarizada en procesos repetibles y reproducibles. La metodología se complementa con una propuesta de modelo matemático basado en la guía para la expresión de la incertidumbre en la medida - GUM. Este modelo se puede aplicar a los lentes que se utilizan para la calibración de equipos de salud visual, como los lensómetros, los cuales se emplean para el diagnóstico de la miopía e hipermetropía por medio de la evaluación de los lentes sin tener contacto directo con los pacientes. Al aplicar el modelo matemático propuesto, y de acuerdo con los datos experimentales, se obtuvieron resultados óptimos en su incertidumbre máxima expandida de aproximadamente 0,0079 dioptrías en una lente de 0,5 dioptrías, comparados con el reporte realizado por los autores, dado que su trabajo reporta una incertidumbre máxima expandida cercana 0,0086 dioptrías, obteniendo como conclusión que la aplicación de este procedimiento científico permite a los fabricantes, y a los usuarios de este tipo de lentes, una confiabilidad en sus mediciones por medio de un proceso de calibración basado en la óptica geométrica en torno a la seguridad del paciente.

    Palabras clave: Metrología óptica, Función de Calibración, Potencia de la lente, Distancia focal, incertidumbre de medición.

    1. INTRODUCTION

    Since 2019, the Colombian Ministry of Health and Social Protection has established several regulations aimed at surveilling medical devices for human use, according to the provisions set out in Resolution 3100 of 2019 [1]. This resolution states that the performance of biomedical equipment, such as ophthalmic devices, should be evaluated regularly. Said regulations seek to guarantee the reliability of the measurements in the diagnoses made using these devices and eliminate possible errors in diagnoses associated with any malfunction of the equipment or deviation from its original function. Only one accredited laboratory in Colombia has complied with these new regulations and demonstrated the reliability of the results in this type of equipment [2]. Such laboratory calibrates ophthalmic equipment by directly measuring the curvatures of the lenses according to their manufacturer’s specifications. However, it lacks a process that guarantees the traceability of the lenses used to calibrate this equipment. Therefore, it carries out validation processes because, currently, there is no standard for this purpose at the local or international level.

    Most of the gaps in the literature in this area are due to the fact that there is no record of documents that present a standardized method to calibrate ophthalmic lenses, which are, in turn, used to calibrate ophthalmic equipment such as lensometers or keratometers. Similarly, the existing protocols for different biomedical equipment exhibit gaps, and, for that reason, a management model for legal metrological control and conformity assessment of biomedical equipment has been introduced to facilitate reliable measurements when this type of technology is used [3]. Also, different advances have been made in the estimation of uncertainty in biomedical equipment by applying non-stochastic methods, such as the Guide to the Expression of Measurement Uncertainty (GUM) [4]. The GUM has been applied to the control of sphygmomanometers and their legal implications in the field of metrology [5].

    In Brazil, the National Institute of Metrology, Quality and Technology (INMETRO) has developed regulations for the mandatory verification of medical equipment used to weigh adults, pediatric scales, and sphygmomanometers [6]. In addition, the Institute of Metrology of Bosnia and Herzegovina (IMBIH) [7] highlighted the importance of applying metrology to clinical medicine, especially to standardize the norms for the inspection of medical devices. Other authors have evaluated and applied measurement uncertainty to determine emissions from fixed sources [8]. These studies are mentioned here because this paper presents a review of the methodologies most commonly used to estimate uncertainty (e.g., the non-stochastic methodology of the GUM) and their relationship with stochastic methodologies (e.g., the Monte Carlo method).

    Ophthalmologists use frontophotometers or lensometers to measure the dioptric power of ophthalmic lenses. These devices employ a composite lens system to converge a light beam and geometrical optics to indirectly measure the dioptric power of lenses [3],[5],[6]. Phoropters are used to measure the refractive error of the eye, and other studies have proposed a calibration method for this equipment [7]- [9].

    The existing literature in this field includes some reports of uncertainty in the measurement of the dioptric power of intraocular lenses [10]. However, regarding normal (or common) ophthalmic lenses, the reports of uncertainty have only been focused on the measurement of focal length and not dioptric power. There are many ways to measure the focal length, and all of them are well supported by physics, as well as good experimental results. Some methods reported in the literature include techniques based on the Talbot effect [10]-[18], while some others use moiré reflectometry [19]-[22], Fizeau interferometry [23], and techniques such as Lau interferometry [24], [25], digital Fourier transform [26], and Fresnel diffraction [27]. In nearly all these studies, the common denominator is the high cost of the required equipment, as well as the high instability of the experimental setups, and, hence, the need for a highly stable laboratory assembly. Thus, these techniques present experimental difficulties that are hard to overcome and whose solution, in many cases, involves additional economic costs.

    The specialized literature about the percentage of error or uncertainty obtained in different measurements of focal length includes the article by Nakano and Murata [14], who calculated an uncertainty in focal length in the order of 0.01 % using the Talbot effect. In turn, Glatt and Kafri [19]-[28] used moiré interferometry and obtained an uncertainty in the order of 2.5 %. Uncertainties of 0.8 % have also been obtained using Fresnel diffraction [29].

    In addition, another method used electromagnetic waves to measure the distance by applying Heisenberg’s uncertainty principle [30].

    However, the measurement of dioptric power is more accurate when there are fewer optical elements involved in the light’s trajectory, such as lenses, mirrors, or filters. Thus, the most appropriate way to measure dioptric power is by directly measuring the focal length of the lens. This is done by physically converging light rays from infinity in the case of positive lenses. In the case of negative lenses, an auxiliary positive lens is used to achieve the convergence of the rays and take an “indirect” measurement of the focal length by means of the experimental determination of the “back focal length.” This process can be performed using different techniques or physical approaches, such as Fresnel diffraction, which, unlike the method proposed in this study, considers an angle of incidence [31]. Thus, the method proposed here produces better measurable results.

    In this paper, we develop a method to experimentally determine dioptric power that is more precise than those found in the literature and does not depend on any additional factor other than the physical phenomenon itself. This study also presents an effective method to measure the uncertainty associated with the measurement of the dioptric power of ophthalmic lenses based on physical principles, clearly differentiating between positive and negative lenses. For this purpose, it was necessary to develop a mathematical model that accounts for the uncertainty associated with this measurement.

    2. EXPERIMENTAL SETUP AND METHOD

    There is no single method to measure focal length because, in the case of positive lenses, it can be measured directly by converging rays from infinity, while the focal length of negative lenses is measured indirectly, and an auxiliary lens is needed. In any case, the measurement of dioptric power is indirect in all methods because the value of the measurand is obtained by transforming, converting, or calculating other direct measurements. For all lenses, the dioptric power p is given by (1):

                                         (1)

    where f is the focal length of the lens, which can be positive or negative. In either case, the appropriate mathematical method and measurement must be implemented.

    2.1 Measurement method

    To determine the focal length of a positive lens, there must be a set of rays parallel to the optical axis coming from infinity. The physical property of a positive lens consists of making the rays coming from infinity converge to the focal point (Figure 1). A very precise measurement of the focal length of a lens can be obtained with an optical assembly that facilitates such an arrangement of rays. The precision of the measurement of the focal length will depend to a great extent on whether the incident rays are considered parallel to the optical axis, i.e., whether they can be considered to come from infinity (see Figure 1).

    Figure 1. Light rays coming from infinity onto a positive lens and converging to its focal point
    Source: created by the authors.

    The characteristics of a negative lens (also called diverging lens) are different from those of its positive counterpart because this type of lenses cause the rays coming from infinity to diverge. Therefore, their focal point cannot be directly observed, as illustrated in Figure 2. In this case, another auxiliary lens should be used to make the rays converge, thus establishing a precise optical criterion to measure the focal length of the diverging lens.

    Figure 2. Behavior of light rays coming from infinity onto a negative lens
    Source: created by the authors.

    According to geometrical optics, when two thin lenses one positive and one negative are aligned on the optical axis and separated by a d distance, as shown in Figure 3, there is not a single focal length, but rather two: a front focal length (FFL) and a back focal length (BFL), which will be abbreviated as B for the remainder of this paper.

    Figure 3. Optical system composed of a negative and an auxiliary positive lens to measure the focal length of the negative lens
    Source: Created by the authors.

    To measure B, we know from geometrical optics that it is given by (2):

                                       (2)

    Then, from (1), the focal length of the negative lens (f1) is solved in (3) in terms of other variables that can be measured:

                                 (3)

    2.2 Experimental setup

    To guarantee the accuracy of the calibration function, the measurements were made in an area with controlled relative humidity and temperature and positive pressure inside it. The environmental conditions used in this study to conduct the tests are: Temperature 22 °C ± 2 °C and relative humidity (RH) 40 % – 60 %.

    The relative humidity and temperature values are optimal for the measurement process because these conditions lie within the acceptable range to measure focal length.

    Conversely, if these conditions are not controlled, the measurements of the different focal lengths can be affected because thermal conditions have an impact on material expansion.

    To make the measurements, the experimental setup should include a system of thin lenses and a microscope objective on the same axis, as illustrated in Figure. 4. A microscope objective, which is located at the laser output, expands the spot of the beam by making it diverge. Then, the beam hits the L1 negative lens, which increases the divergence of the spot until it reaches a considerable diameter. Then, the L2 positive lens, which is located next, should be positioned with great precision because obtaining a widened spot of constant diameter depends on this step. Once the constant diameter spot has been obtained, it can be considered to have a set of rays coming from infinity. As a result, the focal length of any lens in the L3 position can be accurately measured directly using a previously calibrated tape measure. This setup and this spot are used to find the focal lengths of positive lenses directly.

    Figure 4. Experimental setup to determine a positive focal length.
    The microscope objective and lenses L1 and L2 are used to widen the beam
    Source: Created by the authors.

    To measure the focal length of negative lenses, an auxiliary positive lens is used to converge the rays. In this type of setups, the spot is widened between lenses L2 and L3, but this can be easily verified by placing a screen at least 5 m (infinity) away from L2 and checking if the diameter of the spot remains constant.

    Figure. 5 shows the device used for a negative lens such as L3. This setup requires adding an auxiliary positive lens (L4) to achieve the convergence of the rays at a B distance and then indirectly measuring its focal length (f1) using (3) since d, b, and the focal length of L4 are known. In this case, measuring the focal length requires the measurements of the distances of interest involved in (3), which are taken directly using a previously calibrated tape measure.

    Figure 5. Experimental setup to determine the B distance in order to obtain the focal length of a negative lens (L3).
    In this case, the L4 positive auxiliary lens is used to converge the beam since L3 is divergent
    Source: Created by the authors.

    2.3 Mathematical model for determining uncertainty and error in the calibration process of ophthalmic lenses

    As established in the measurement method, the assembly was prepared in a laboratory, which facilitated the application of the measurement schemes proposed in Figures 4 and 5 for positive and negative lenses, respectively. Figure 6 shows the assembly, which only requires a rail with lens mounts aligned with the table on which the focal length measurements are made using a calibrated tape measure.

    Figure 6. Laboratory setup. Left: laser and lenses L1 to L3. Right: mechanism of beam expansion in operation
    Source: Created by the authors.

    The mathematical expression used here to determine the uncertainty is based on stochastic models; therefore, it is based on the Guide for the Expression of Measurement Uncertainty (GUM) [4]. Note that the error of the measurement being tested must be calculated based on the following input variables given by (4):

                             (4)

    where f1 is the focal length of the lens; di, the deviation of the measurement at point j, which corresponds to each point taken from the focal length in relation to the reference value; CP, the correction due to the pattern, which corresponds to the tape measure; CA, the correction to compensate for the misalignment effect between the ruler and the table, according to the setup shown in Figure 6; CT, the correction due to differential thermal expansion; and CDE, the correction due to the scale of the tape measure. The deviation of the dj measurement is given by (5):

                                            (5)

    where xi is the measurement of each focal length of the lens x ̅, the mean of the measurement of the focal distances of the lens; and n, the number of data obtained. The CP correction is given by (6):

                                            (6)

    where ucoj denotes the product of the expanded uncertainty of the tape measure (Ueoj) multiplied by the 𝑘 coverage factor, which was previously calibrated at a laboratory at Instituto Tecnológico Metropolitano (ITM) in Medellín, Colombia. The CA correction is given by (7):

                                                    (7)

    where θ2max is the correction of the measured angle due to perpendicularity loss, and 𝑥𝑜𝑖 is the measurement of the focal distance of the lens. The CDE correction is given by (8):

                                                           (8)

    where e is the scale being used. In this study, it is millimeters (i.e., 0.001 meter).

    According to the law of propagation of uncertainty, the expression of the combined standard uncertainty, uc2 (P), , for positive lenses (assuming there is no correlation between the variables) is given by (9):

                  (9)

    Applying the expression above to the function of ophthalmic lens calibration and considering the sources of uncertainty established in (1) and (4), the output variable is given by (10):

                               (10)

    Applying (8) in (1) and (3), we obtain (11):

                                   (11)

    The estimate of the effective degrees of freedom of the standard uncertainty, Uc (P), associated with the output estimate is obtained using the Welch–Satterthwaite formula [4], given by (12):

                                      (12)

    The k coverage factor can be obtained using this equation, which is derived from a table of values and based on a Student’s t-distribution evaluated for a coverage probability of 95.45 %. A coverage factor of k = 2.0 is used in this study, which must be multiplied by the combined uncertainty, uc, to find the expanded uncertainty, Ue, using (10) and (11) with the result shown in (13).

                                             (13)

    For negative lenses, we obtain (14):

                     (14)

    where f2 is the focal distance of the negative lens, and the uncertainties are the same as those of the positive lens, with one difference, i.e., the term for the auxiliary positive lens, 𝐶𝐿𝑃, which is applied to the same procedures to determine its uncertainty.

    According to the law of propagation of uncertainty, the expression of the combined standard uncertainty, uc2 (P), of a negative lens (assuming there is no correlation between the variables) is the same as that presented in (3). Applying the previous step to the ophthalmic lens calibration function and taking into account the sources of uncertainty declared in models (1) and (14), the output variable is given by (15):

                                            (15)

    Applying (14) in (1) and (9), we obtain Expression (16):

                                      (16)

    The estimation of the effective degrees of freedom of the standard uncertainty, uc (p), associated with the output estimate is obtained using the previously mentioned Welch–Satterthwaite formula (12). Likewise, a k coverage factor of 2.0 is used here based on a student’s t-distribution for a coverage probability of 95.45 %. These values are taken into account to find the expanded uncertainty in a similar way to the previous case.

    3. EXPERIMENTAL RESULTS

    The results obtained with the proposed procedure show the estimation of uncertainty in the calibration of ophthalmic lenses in accordance with what is established in a non-stochastic methodology such as the Guide for the Expression of Measurement Uncertainty (GUM) [4 ], which has been widely used in different processes, e.g., chemical metrology [30], gamma ray spectrometry [32], and electromagnetic compatibility [33].

    The GUM establishes a general structure for evaluating and expressing uncertainty in measurements that can be applied to multiple measurement processes with different levels of accuracy and precision. Moreover, the principles in this guide are intended to be applicable to a wide range of measurements. The steps proposed in the GUM, which were widely used in this paper, are followed to identify and characterize the sources of uncertainty and estimate combined and expanded uncertainties.

    To implement the method proposed here and considering that the physical phenomena being intervened are diopter distances, the contributions in document DI-011 [34] by the Spanish Metrology Center were also taken as a reference regarding the concepts of dimensional metrology. Therefore, the measurements were taken using the method proposed in [34], and the results thus obtained are presented in Tables 1 and 2. The measurements were taken linearly for each diopter being verified, taking ten measurements from each one. These points were defined according to the experience of the laboratory, considering that the points included here cover a range of lenses normally used in the medical field. Table 1 shows that, in the entire set of experimental measurements, the deviations of each optical lens are very low, lower than one diopter. Since focal lengths are measured directly and dioptric power (DP) is calculated afterward, very precise and accurate reference measurements (RM) were found in these experimental data.

    Table 1. Focal lengths measured at the laboratory
    Source: Created by the authors.
    LensesRMm1m2m3m4m5m6m7m8m9m10Aver.DP
    l110.0000.1000.1010.0980.0990.1000.1000.1030.1030.1020.1030.1019.911
    l210.0000.1030.1030.1020.1000.1020.1020.1040.1060.1020.1010.1039.756
    l33.3300.3100.3050.3110.3090.3060.3020.3050.3040.3080.3120.3073.255
    l43.3300.3080.3050.3100.3080.3080.3040.3040.3010.3080.3120.3073.259
    l50.5002.1302.1402.1552.1602.1512.1552.1602.1502.1652.1702.1540.464
    l6-20.000-0.050-0.049-0.049-0.049-0.049-0.050-0.047-0.049-0.052-0.049-0.049-20.301
    l7-20.000-0.050-0.049-0.050-0.050-0.050-0.049-0.054-0.051-0.048-0.049-0.050-20.004
    l8-5.000-0.192-0.192-0.192-0.192-0.192-0.195-0.193-0.191-0.188-0.189-0.192-5.215
    l9-5.000-0.196-0.189-0.190-0.192-0.192-0.191-0.190-0.189-0.191-0.184-0.190-5.252
    l10-4.000-0.263-0.256-0.261-0.260-0.260-0.253-0.248-0.247-0.250-0.252-0.255-3.921
    lAUX10.0000.1000.1000.1010.1030.1030.1010.1000.1000.1030.1020.1019.872

    Table 2. Calibration corrections and associated uncertainties
    Source: Created by the authors.
    Lens(dj)ErrorCombined uncertaintyExpanded uncertainty
    l10.000600.089000.000620.00120
    l20.000550.244000.000570.00110
    l30.001100.075000.001120.00220
    l40.001090.071000.001110.00220
    l50.003920.036000.003950.00790
    l60.000440.301000.000460.00093
    l70.000530.004000.000550.00110
    l80.000690.215000.000720.00140
    l90.001060.252000.001080.00220
    l100.00197-0.079000.002000.00400
    lAUX0.000450.128000.000470.00094

    The uncertainty in the ophthalmic lens calibration function is estimated after obtaining the results in the measurement process, applying (9), (10), (11), (12), (13), (14), and (16) to the data in Table 1. Hence, an error and an expanded uncertainty, 𝑈𝑒, are obtained for each measurement point in Table 2.

    Table 1 lists the focal lens measurements with index 1, and subindices in numerical order. In addition, lAUX denotes the auxiliary lens used in the setup for measuring distance B. RM is the reference measurements; DP, dioptric power; and m, the measurement with an index that corresponds to the measurement order.

    Table 2 shows the data obtained after applying the mathematical model based on the law of propagation of uncertainty. This law is also used to express the calculation with a type A uncertainty (dj); a type B uncertainty provided by the reference (Cp; the uncertainty due to compensation for the sliding effect between the ruler and the table which is constant with a value of 0.00478 (Ca); the uncertainty due to differential thermal expansion correction and will always taken as zero (Ct); and the uncertainty provided by the scale division of the tape measure which is also taken as zero(Cde). Subsequently, the combined uncertainty is calculated using (16) resulting in a constant value of 0.00041. For the expanded uncertainty, the effective degrees of freedom should be estimated first by applying (12), using a k coverage factor of 2.0 based on a Student’s t-distribution for a coverage probability of 95.45 %. This coverage factor is multiplied by the combined uncertainty so that the expanded uncertainty can be found at each point, as shown in Table 2.

    According to the results obtained with the calibration function, a similar behavior can be expected in each repetition of the measurement executed for each nominally true value with the reference diopters, thus obtaining errors of around one diopter, whether for positive or negative lenses. The predominant uncertainty in most points is provided by the calibration of the standard being used (CP). Thus, the other sources that would contribute to the expanded uncertainty are the uncertainty provided by the repeatability (dj) and that provided by the scale division of the tape measure (Cde). This is because the sources due to compensation of the sliding effect between the ruler and the table (Ca), the source due to the differential thermal expansion correction (Ct), and the environmental conditions applied in the laboratory are considered null by the method.

    Finally, considering the recommendations provided by the Spanish Metrology Center [15], we obtained optimal results regarding deviations and a maximum expanded uncertainty of ± 0.0079 diopters in a 0.5-diopter lens.

    4. CONCLUSIONS

    This article presented a procedure for calibrating ophthalmic lenses. Additionally, it proposed a mathematical model for estimating uncertainty based on a non-stochastic methodology such as the Guide for the Expression of Measurement Uncertainty (GUM) [11] and the application of dimensional metrology concepts. This study also identified sources of uncertainty that had not been previously observed or analyzed in most research papers in the literature. The experimental results reported here show that, derived from the specifications of the equipment being tested, i.e., the uncertainty that most affects each of the measurements is that provided by the calibration of the standard being used. The proposed methodology is novel because it demonstrates the importance of calibrating ophthalmic lenses to obtain reliable measurements as well as estimating the uncertainty using a structure that combines a mathematical model with a non-stochastic method in accordance with the GUM.

    However, this method may present some limitations compared to existing ones (e.g., the interferometry method) because the measurements are made manually by the operator. In addition, not using a robotic or automated system can lead to human error, and, although this is considered in the estimation of uncertainty, it still is a limitation that would be easily overcome with financial resources. Nevertheless, this study is important because it investigates the reliability of measurements of biomedical equipment, specifically ophthalmic lenses. In addition, it provides relevant information for ophthalmic lens manufacturers because the maximum expanded uncertainty of the method proposed here was optimal: ± 0.0079 diopters in a 0-5 diopters lens. By contrast, other authors [1] have reported a maximum expanded uncertainty of ± 0.0086 D.

    The focal length measurement procedure used here offers two advantages: simplicity of the assembly and low cost. None of the papers reviewed in this study describes a method based on a cheaper assembly that also presents high stability and easy operation. The current disadvantages of this process are associated with its rudimentary and low-cost implementation that does not use any electronic sensors or measurement components. However, this can be overcome by obtaining financial resources to purchase more accurate measuring instruments.

    A future line of research is the application of this calibration function to other ophthalmic equipment based on physical principles to guarantee the validity of the results in order to obtain reliable measurements in eye diagnostics.

    5. ACKNOWLEDGMENTS

    This study was conducted in the framework of the project developed by Alejandro Salgar-Marin as a young researcher, whose general objective is “To define a calibration method for optical lenses applied to the metrological assurance of ophthalmic instruments.” This study was funded by the Instituto Tecnológico Metropolitano (ITM) in Medellín, Colombia, under the program called “Young researchers and innovators,” which is part of the institutional development plan: “ITM to another level.” The data in this paper were obtained at the physics laboratory and analyzed at the AMYSOD laboratory located in the Parque i research laboratory center at the ITM. Professor Javier Vargas thanks the Faculty of Exact and Applied Sciences, ITM, for their financial support.

    CONFLICTS OF INTEREST

    The authors declare that there is no conflict of interest.

    AUTHOR CONTRIBUTIONS

    Alejandro Salgar-Marin: Laboratory work, optical setup, manuscript drafting, research.

    Javier Alberto-Vargas: Laboratory work, optical setup, manuscript drafting and revision, research.

    Andrés Felipe Ramirez-Barrera: Laboratory work, manuscript drafting and revision, research, mathematical model development.

    6. REFERENCES