Imaginary optical constants in near-infrared (NIR) spectral range for the separation and discrimination of adulterated diesel oil binary mixtures

Fuel adulteration is a major factor contributing to environmental pollution. The tendency to adulterate diesel oil with kerosene is becoming prevalent in many countries, which is particularly driven by the lower cost of kerosene than of diesel oil. The objective of this study is to develop a method of separating and discriminating between different compositions of kerosene in binary liquid fuel mixtures. In this study, we have utilized refractive indices of binary fuel mixtures measured using an Abbe refractometer, and the inversion of measured transmittance spectra in the visible and near-infrared (NIR) spectral range by a general method based on Kramers–Kronig dispersion analysis, where the singly subtractive Kramers–Kronig relations have been used to obtain the wavelength-dependent imaginary excess permittivity of a binary liquid mixture. The excess permittivity has considerable importance in the interpretation of different liquid molecule interactions in binary liquid mixtures. It is herein demonstrated that the problematic cases of relatively low adulterations of 5% and 10% can be detected and distinguished from 15% adulteration using imaginary optical properties in the NIR spectral range. We show that imaginary optical constants such as excess imaginary permittivity and the extinction coefficient are useful for not only screening for adulteration, but also discriminating and separating low from high kerosene compositions in mixtures. Moreover, these quantities can also be exploited for the development of practical sensors, especially for sensing diesel oil adulteration under field conditions.


Introduction
Liquid fuel adulteration is becoming a common problem in many countries. This unethical practice is carried out to increase profit. The underlying reason for this practice is the uneven taxation system for different liquid fuels, where kerosene is generally subsidized, while diesel and gasoline have the higher tax rates. This illegal practice leads to adverse effects such as damage to engines, loss of taxation for authorities, increased traffic accidents, and environmental pollution resulting from exhaust gas emissions [1,2].
The most common form of fuel adulteration in many parts of the world is to add kerosene to either diesel oils or gasoline. These fuels consist of hydrocarbons, which have spectral fingerprints in the near-infrared (NIR) spectral region. Nevertheless, it is challenging to screen diesel oils for adulteration with kerosene, especially when the content of the adulterants tends to be below 20%, using NIR spectra owing to the closely overlapping fingerprints of the two liquids. However, the development of optical measurement methods for adulterated fuels is feasible because of the different optical properties of these liquid fuels, which can usually be detected rapidly by different and relatively simple measurement devices. Furthermore, it is possible to perform these measurements in the field. Both spectral and non-spectral optical measurement methods have been proposed as methods for screening adulterated fuel [3][4][5][6].
For the rapid optical evaluation of fuel quality, an Abbe refractometer is a useful device for screening adulterated diesel oils via the refractive index [6]. However, there are cases, where the refractive index of diesel oil is very close to that of kerosene, and such a case is demonstrated in this study. Therefore, under such circumstances further data, in addition to the refractometer reading, is required to confirm sample adulteration. In our previous work, we have demonstrated the use of a prototype hand-held sensor for the detection of diesel oils adulterated with kerosene. The prototype is based on a modified commercial hand-held glossmeter. However, the measurement time with this type of sensor is still long compared with the almost instant reading with an Abbe refractometer or the relatively rapid NIR measurement with a (portable) spectrophotometer. Another type of hand-held sensor for screening adulterated diesel oils is based on an Abbe refractometer and a data analysis method using the real part of excess permittivity with a fixed probe wavelength [7]. Regarding field measurements, NIR spectrum measurement has also been tested for crude oil quality inspection under difficult field measurement conditions [8].
In our recent studies of adulterated diesel oils, we have learned that when using either an Abbe refractometer or a spectrophotometer, there are cases when one may need additional complementary methods to reliably confirm diesel oil adulteration. To circumvent this issue, we have proposed a combination of both refractive index and NIR transmittance measurements, along with exploitation of the so-called excess permittivity of a binary mixture to obtain a more comprehensive picture of adulterated diesel oils [9]. In other words, we have shown that the real part of the complex excess permittivity of adulterated diesel oil reveals fine structures of spectral fingerprints that can be used to screen adulterated diesel oil. The study reported in this paper is a continuation of [9], namely, we have investigated the imaginary part of the complex excess permittivity of authentic and adulterated diesel oils. It is shown that the combination of the real part of the refractive index, the imaginary part of the refractive index (obtained from the absorption coefficient of the fuel), and the imaginary part of the excess permittivity of adulterated diesel oils, provides a plausible means of deciding with confidence the relative amounts of kerosene adulteration in diesel oils.
To exploit the imaginary excess permittivity in adulterated diesel oil, we utilize the complex-valued excess relative permittivity, which is a useful measure for studying the difference between the measured and ideal permittivities of a liquid mixture. This relative excess permittivity is related to the complex refractive index of a binary liquid. Furthermore, the excess permittivity can be used as a measure of the interaction between molecules of different liquid components in binary liquid mixtures [10,11]. As examples of this type of application, it is possible to assess the validity of different mixing rules for the optical constants of binary liquids by studying perturbed hydrogen bonding in the IR region and the corresponding excess absorbance [12], and to study the excess absorbance of water-alcohol mixtures as well as ionic liquid mixtures in the terahertz (THz) gap [13,14]. The theoretical foundations of the complex excess permittivity which depends on the frequency of the incident electromagnetic field (the frequency-dependent complex excess permittivity is crucial in our case), have been derived relatively recently for binary and multiple liquid mixtures [15]. In this paper, we are concerned with the imaginary part of the wavelength-dependent excess permittivity. The complex permittivity of a binary liquid mixture depends on the displacement of electrons (light interaction) and molecules (IR interaction), and the molecular structure, size, and weight of the molecules. These factors are beyond the scope of these studies, because there are hundreds of different hydrocarbons in diesel oils.
Here, we report, to the best of our knowledge, the first ever study on liquid fuel adulteration based on imaginary excess permittivity focusing on diesel oil adulteration with different volume fractions of kerosene. The imaginary excess permittivity is calculated using NIR transmission spectra measured with a spectrophotometer, the refractive index measured with an Abbe refractometer, and using the singly subtractive Kramers-Kronig (SSKK) relation [16] to obtain an estimate of the wavelength-dependent refractive index, which is needed in the calculation of the imaginary excess permittivity. The basis of the SSKK relation is a rather fundamental property of physics and chemistry, namely, causality [16]. Therefore, SSKK analysis has general validity and can be used regardless of the spectral line shape model assumed for. The validity of the method is exemplified using diesel of two different grades adulterated with kerosene of the same grade.
To provide a more comprehensive picture of screening problematic adulterated diesel oils, here, we study the extinction coefficient and imaginary excess permittivity. The former is obtained directly by measurement, whereas the latter is calculated.

Theory
The imaginary excess permittivity can be derived from the equation for the excess permittivity of binary mixtures. If we denote the complex relative permittivity of diesel oil by ε D and the corresponding permittivity of kerosene by ε K , we can express the complex excess relative permittivity (here after we use the expression for permittivity, meaning that it is the dimensionless relative permittivity) as [9] (1) Here, ε is the complex permittivity obtained from optical measurements, ε ideal is the complex ideal permittivity of the binary mixture without any interactions between unlike liquid molecules, f D is the volume fraction of diesel oil, and ε D /ε K is the permittivity ratio. The sign and magnitude of ε E are a measure of the strength of interactions between molecules in a binary mixture. There are hundreds of different hydrocarbons in diesel oil [17], and after adulteration by mixing with different percentages of kerosene, the dipole moments of the electric charges of the binary mixture will change, causing small but observable differences in the measured spectral features. This change also leads to differences in the real and imaginary parts of the excess permittivity, which is the key idea of this paper. The ideal permittivity in Eq. (1) corresponds to the upper Wiener bound [18] of an effective medium as follows: In general, the expression for the frequency-dependent complex permittivity of a homogeneous medium is given by where ω is the angular frequency of incident electromagnetic radiation and i is the imaginary unit. The complex refractive index is given by N(ω) = n(ω) + ik(ω), where n is the conventional refractive index and k is the extinction coefficient of the medium. The complex relative permittivity of an insulator such as a pure liquid or one of the binary mixtures in this study is given by ε = N 2 , which was used in the second equality in Eq. (3).
Using the first equality in Eq. (1) and second equality in Eq. (3), the imaginary excess permittivity is given by where λ is the wavelength of light, n(λ) and k(λ) are the wavelength-dependent refractive index and extinction coefficient of an adulterated sample, respectively. In our case, n ideal (λ) and k ideal (λ) are the wavelength-dependent refractive index and extinction coefficients of ideal binary mixtures, n D (λ) and k D (λ) are the wavelength-dependent refractive index and extinction coefficient of authentic diesel oil, and n K (λ) and k K (λ) are the wavelength-dependent refractive index and extinction coefficient of kerosene, respectively. The wavelength-dependent refractive indices and extinction coefficients of authentic diesel oil and kerosene, n D (λ), k D (λ), and n K (λ), k K (λ), as well as those of adulterated diesel oils n K (λ) and k K (λ), respectively, were obtained in several steps. First, the transmission spectra and refractive index of the sample were measured, then, Lambert's intensity law (3) (Eq. 5), the extinction coefficient (Eq. 6), and the SSKK relation (Eq. 7) were applied.
Lambert's intensity law is given as from which one obtains the wavelength-dependent absorption coefficient α of the medium for an a priori known sample thickness d. Hence, we can calculate k from the formula: In this work, we exploited the wavelength of Vis-NIR radiation rather than the angular frequency. Because we can measure the discrete refractive index for an anchor point with an Abbe refractometer, we further utilize the SSKK dispersion relation to obtain n using [16] Here, n(λ′) is the wavelength-dependent refractive index, n(λ″) is a priori known refractive index at an anchor point λ″, and P denotes the Cauchy principal value. To the best of our knowledge, this is the first time that Kramers-Kronig analysis has been applied to the wavelengthdependent imaginary excess permittivity of binary liquid mixtures. In this work, the fixed wavelength λ″ is the same as the wavelength of the Abbe refractometer, namely, 589 nm. In the practical calculation, the integral, especially in the case of SSKK analysis, can usually be truncated to correspond to the initial and final points of the measured transmittance. Here, we use a method to optimize the anchor point λ″ to minimize the possible data inversion error. This method is based on the properties of a Chebyshev polynomial shown by Palmer et al. [19]. Because of the fixed wavelength of the Abbe refractometer, we are not free to choose the anchor point. Nevertheless, we employ the optimal anchor point by fixing the wavelength of the final point (λ f ) of the transmittance spectrum and we calculate the initial wavelength (λ i ) that corresponds to the optimal anchor point. The wavelength λ i is then obtained from a zero of a Chebyshev polynomial of the first kind as follows [19]: Hence, the scanning of the transmittance spectrum was started from the calculated wavelength. The SSKK method has been used, for example, to analyze the properties of diesel oil droplets [20], but not to investigate the imaginary part of the excess permittivity of adulterated diesel oils. In the .
SSKK method, the location of the anchor point is not usually a critical factor for the successful inversion of extinction coefficient data.

Materials and methods
In this work, we studied summer-grade (sample A) and winter-grade [sample B (− 20 °C)] diesel oils that originate from the same crude oil, to take into account different climatic conditions throughout the year in different countries. The first very important additive in diesel oils (both summergrade and winter-grade) is wax, which is normally added to increase the combustion speed leading to more power owing to its higher cetane value. However, during winter (low operating temperatures), wax tends to solidify and cause problems during ignition; therefore, anti-wax additives are added to avoid solidification and enable smooth engine operation. Thus, the main difference between the summer and winter diesel compositions is the inclusion of anti-wax additives, which are present in winter-grade diesel oils but not in summer-grade diesel oils. The value of − 20 °C is the lowest temperature at which a diesel engine is expected to start. The kerosene sample in this study is a low-odor commercial product (sample C, Alfa Aesar). It is important to emphasize the variation of the optical properties of diesel oil and kerosene from different oil fields, because these are complex liquids with hundreds of hydrocarbons [17]. Owing to this fact, we chose a well-known laboratory quality kerosene as a reference that, in principle, is available to scientists worldwide. Obviously, the refractive index of the diesel oil and kerosene depends on the origin of the crude oil from which the fuel samples were refined. The refractive index of pure and adulterated diesel oil was measured at room temperature of 22.7 ± 0.02 °C using an Abbe refractometer (Atago RX5000) having an operating wavelength of 589 nm and a relatively high accuracy of ± 0.00004 refractive index units. Differences in the measurement temperature will affect the refractive index reading owing to changes in the thermo-optic coefficient of the fuel sample. Therefore, it is important in the field measurements to pay attention to the temperature variations to avoid possible misinterpretation of results.
The NIR transmittance spectra of the samples were measured at a room temperature of 23 °C using spectrophotometer (PerkinElmer Lambda 9) and a 10-mm-thick quartz cuvette. The thick cuvette was used to test the applicability to various field conditions. Before recording the initial spectra of the fuel samples, a baseline correction measurement was performed by recording the spectra of an empty cuvette: this baseline was recorded in the memory of the software of the spectrophotometer and was automatically considered in subsequent measurements. Furthermore, multiplicative scatter correction (MSC) pretreatment was performed on the spectral data using PLS-Toolbox (Eigenvector Research, Inc USA) in the MATLAB software, to minimize or eliminate Pathlength-related errors, baseline shifts, and interference. To ensure that our spectral data were suitable for further processing, we compared the pretreatment in MATLAB with MSC pretreatment carried out using UnscramblerX (CAMO, Norway). The spectrophotometer measures transmittance in the UV-Vis-NIR range with an accuracy of 0.07% for an absorbance value equal to 1. However, the target spectral range for this work was the Vis-NIR range. Our choice of wavelengths of considerably below 2000 nm for the spectra was based on the availability of such portable spectrometers on the market suitable for the detection of hydrocarbons in diesel oils which operate below 2000 nm. We remark that, during these studies, we also measured the transmittance using 10-mm-thick and 5-mm-thick cuvettes filled with fuel samples and observed the transmittance ratio from the obtained transmittance spectra to avoid any discrepancy of data due to an empty cuvette. However, we did not find any significant discrepancy for the samples of this study. For our target field measurement conditions, the utilization of a single cuvette, especially one with a thickness of 10 mm, is more practical than that of two cuvettes, because the former cuvettes are cheaper and can be easily rinsed than 5 mm and 1 mm cuvettes. Nevertheless, because of possibility of contamination of the cuvette during the field measurements, specially designed single cuvette with two different optical path lengths that incorporates a beam splitter to guide the probe beam via the two paths, and thus provides the transmittance ratio would be useful. Such a procedure would help to eliminate noise due to possible contamination of the cuvette.

Results and discussion
We show the refractive index data for pure fuel samples A, B, and C, measured by the Abbe refractometer at 589 nm in Table 1. It is clear from Table 1 that the summer-grade diesel oil (sample A) has a higher refractive index than the winter-grade diesel oil (sample B) and the kerosene (sample C). Furthermore, the refractive indices of samples B and C differ only in the third decimal place. The refractive index of kerosene sample C is very close to that of winter-grade diesel oil sample B, which indicates the problem that one can expect only a minute difference in the refractive index, between the authentic and adulterated diesel oils. Table 1 also gives the refractive indices of the adulterated diesel oils. In all cases, the refractive index of the binary mixture decreases with increasing kerosene volume fraction. For the binary mixtures for sample A, the refractive index readings for 5%, 10%, and 15% kerosene become increasingly close to that of diesel oil sample A. These readings are always between those of sample A and kerosene, in line with the conventional binary mixture model. However, the mixing rules, such as those in [12], for an ideal binary liquid are not valid for these samples. Indeed, we have shown, for example, that the Lorentz-Lorenz (L-L) formula is invalid for the present samples. Thus, the permittivity ratio, especially with increasing amount of kerosene, does not linearly depend on the kerosene fill fraction. Hence, the quantification of the kerosene fill fraction using the L-L formula is problematic.
For sample B adulterated with 5% and 10% kerosene, the refractive indices are slightly above the reading of pure diesel oil, whereas that for sample B adulterated with 15% kerosene is slightly below it. In fact, this behavior for sample B is contrary to that usually expected for binary mixtures, namely, under ideal conditions, the refractive indices of mixtures are expected to lie between those of the two mixture constituents. Sample B adulterated with 15% kerosene could be erroneously concluded to be authentic diesel oil if and only if the refractive index reading is considered, namely, 1.44731 (15%) and 1.44775 (0%) are close to each other. When the amount of authentic diesel oil in the 15% adulterated sample is increased to reduce the adulteration to 10% and finally 5%, the abnormal behavior of the refractive index becomes evident, and hence, such a dilution process for the originally 15%-adulterated diesel oil increases the confidence in the judgement that the diesel has been adulterated. In our recent article [7], we demonstrated that sometimes when the refractive indices values for both authentic and adulterated samples are very similar, it might lead to authentic diesel erroneously being interprated as adulterated. Therefore, a more rigorous measure is required to ascertain the differences between samples. Moreover, in the same article, the insignificant effect of small temperature variations on the refractive index was explained. The article placed more emphasis on the role of chemical interactions in invalidating the L-L formula if one attempts to quantify the amount of kerosene. These interactions are the cause of nonlinear variations of the refractive index, especially for the present case of diesel oils adulterated with kerosene [7]. Figure 1 shows the transmittance curves of pure fuel sample A together with the corresponding adulterated samples in the spectral range of 431-1600 nm. The wavelength of 1600 nm was fixed, and that of 431 nm was calculated from Eq. (8). Diesel oil and kerosene also have a strong absorption band in the range of approximately 1700-2100 nm [21], but we have omitted this band as well as bands in the midinfrared spectral range. As shown in Fig. 1, the curves for authentic diesel oil (sample A) and kerosene (sample C) overlap. However, it is possible to clearly differentiate the authentic diesel oil from the adulterated samples, whereas it is difficult to differentiate the adulterated samples. For the sample with 5% adulteration, it is possible to differentiate the sample from other adulterated samples, especially at visible wavelengths. Moreover, the corresponding transmittance curves for 10% and 15% adulteration show even greater overlap in Fig. 1. For sample B, the screening of adulterated diesel oils is much more difficult regarding the interpretation of the spectral data, and its refractive index is close to that of sample C. This again shows the difficulty of screening adulterated diesel oils using their Vis-NIR spectra and refractive index. From the measured transmittance spectra of the pure and adulterated samples, we calculated the wavelength-dependent extinction coefficient of the samples using Eq. (6). Figure 2a, b shows magnifications of extinction coefficient curves for kerosene, and authentic and adulterated samples A and B in the vicinity of the strong absorption band ranging from 1190 to 1220 nm. It is evident from Fig. 2a, b that the shoulder of the fingerprint of the pure diesel oils at 1195 nm is higher than that at 1214 nm, whereas the opposite holds for the kerosene. Moreover, pure kerosene and diesel oil can already be identified from the relative heights of the two spectral peaks of these samples. However, when one mixes kerosene with diesel oils A and B, the presence of kerosene cannot be judged from the relative heights in this spectral region, because the left shoulder is always higher, thus emphasizing the presence of diesel oil. It is also clear from Fig. 2a, b that for 10%, and 15%-adulterated samples, the extinction coefficient curves overlap, and it is difficult to find differences because of the relative heights of the two nearby bands of diesel oil and kerosene. Hence, there is no clear rule that enables one to use the extinction coefficient to predict the adulteration percentage. Figure 3a, b shows magnifications of the extinction coefficient curves for kerosene, authentic sample A, and authentic sample B with their corresponding adulteration percentages in the spectral range 1390-1427 nm. In contrast to the observation made for the fingerprints between 1190 and 1220 nm, where the difference in the relative height exists between the nearby bands for diesel oil and kerosene, such difference in the relative height is not present for the two bands in the vicinity of the strong absorption ranging from 1395 to 1415 nm. In this spectral range, the maximum extinction coefficient curves at approximately 1395 nm and 1415 nm increase with increasing adulteration percentage. This is clear for both adulterated samples A and B in this spectral region. Furthermore, the curves of the extinction coefficient of authentic sample B and kerosene in Fig. 3b almost completely overlap, but the curves for the adulterated samples are much higher. For authentic sample A and kerosene, in Fig. 3a, both curves are well resolved, but after mixing kerosene and diesel oil, the extinction coefficient curves for the adulterated samples are higher. The maximum extinction coefficient increases with increasing adulteration percentage, which is clearest in the vicinity of 1415 nm. Hence, from the extinction coefficient curves in Fig. 3a, b, we can see an excess value of k, namely, the interaction of the molecules of the two fuels. From Figs. 2 and 3, it is evident that 1390-1420 nm is optimal for screening the adulterated fuels.
According to refractive index readings in Table 1, if we have a relatively highly adulterated sample 15% in this study, we can dilute it first to 10% and then to 5%, with authentic diesel to obtain more information about it. Note that the operating wavelength of the Abbe refractometer is 589 nm, whereas we monitor the extinction coefficient curves, such as those presented in Fig. 3a, b, at much higher wavelengths. Then, the extinction coefficient curves complement the prediction based on the readings of the Abbe refractometer, and we obtain the adulterated samples in the correct order of adulteration from the spectrophotometer readings at 1390-1420 nm. Although the Abbe refractometer reading for sample B with 15% adulteration, which is 1.44731, could lead to it being accepted as authentic diesel oil, because it is close to the refractive index of authentic diesel B, namely, 1.44775, the extinction coefficient data clearly show that it is an adulterated sample. We can thus use both Abbe refractometer data and extinction coefficient data (at 1400 nm) to screen adulterated diesel oils. To determine the excess extinction coefficient using the expression k(λ) − k ideal (λ), the extinction coefficient is applied. However, to clearly distinguish and discriminate between samples with different adulteration percentages, we employ a rigorous theory, given in [15], associated with excess permittivity. Therefore, in this work, in addition to the excess extinction coefficient, we introduce one more measure, which is based on the calculation of the imaginary excess permittivity of diesel oil and kerosene mixtures. As already concluded from the magnified extinction coefficient curves, the optimum spectral range for examining imaginary excess permittivity is in the vicinity of 1395-1415 nm. The imaginary excess permittivity was calculated using Eqs. (4), (6), and (7) within this spectral range. Figure 4a, b shows estimates for the refractive index of the adulterated fuel samples calculated using Eq. (7) and by finite-range integration. These data were used to calculate the imaginary excess permittivity. From Fig. 4a, b, it is clear that the greatest change in the refractive index is in the region of anomalous dispersion in the NIR absorption bands. The SSKK relation gives stronger convergence of the integral than the conventional Kramers-Kronig relation; hence, the accuracy of data inversion is high even for finite-range spectral data. Our aim is to develop an optical sensor that only uses measured finite spectral data in the analysis method to screen adulterated diesel oils. Successful data inversion using the SSKK relation in the case of a narrow spectral range and omission of the nearby spectral fingerprint was demonstrated in [22]. Recently, the validity of the method using the finite-range SSKK relation has been demonstrated even for saturated absorption bands, by simulations of absorption with a Lorentzian in [23].
Next, Fig. 5a shows the wavelength dependence of the imaginary excess permittivity of adulterated sample A in the spectral range of 1000-1600 nm, where the curves were obtained as the difference between the imaginary permittivity of the adulterated samples and that of the ideal binary mixtures. It is clear from Fig. 5a that 5% adulteration has the lowest imaginary excess permittivity, followed by 10% adulteration. Furthermore, the sample adulterated with 15% kerosene has a much higher imaginary excess permittivity than the samples with 5% and 10% adulteration. This demonstrates that using the imaginary excess permittivity, it is possible to separate and discriminate sample A adulterated with different volume fractions of kerosene. Furthermore, Fig. 5b shows the wavelength dependence of the imaginary excess permittivity of adulterated sample B in the spectral range of 1000-1600 nm. In this case, we observe a similar behavior for 5% and 10% adulteration to that in Fig. 5a: one can distinguish two nearby peaks of excess imaginary permittivity that reflect the two absorption bands. However, a single strong peak is observed for 15% adulteration in Fig. 5a, b, which implies that this is not merely a coincidence, because it is present for both diesel oils. The baseline of 10% in Fig. 5b is strikingly high compared with that in Fig. 5a. However, the 15% adulterated sample in Fig. 5a, b can be judged as an adulterated diesel oil owing to its higher peak. The imaginary excess permittivity can be used in decision making to screen diesel oil adulteration by kerosene.
As mentioned already, sample B has a refractive index very close to that of kerosene, and this further complicates the detection of adulteration using only Abbe refractometer readings, because the similar values increase the uncertainty. Hence, from the information of the refractive index measured by the Abbe refractometer, the extinction coefficient obtained from spectrophotometer measurements, as well as the imaginary excess permittivity obtained by a combination of these data through SSKK analysis, we can not only decide with certainty the presence of adulterated diesel oil, but also identify and discriminate low from high levels of adulteration.
In Fig. 5a, b, a strong peak is observed for 15% adulteration, whereas for 5% and 10% adulteration, a strong peak is not present. We suggest that this type of curve provides evidence for the abnormal dependence of the imaginary excess permittivity on the kerosene fill fraction and hence supports the observation that the L-L formula is invalid for adulterated diesel samples, especially when the amount of kerosene is relatively large [7]. The L-L formula for permittivity, indicating a linear dependence on the fill fraction of a foreign liquid, is valid in the absence of chemical activity of the liquids involved in the binary mixture. On the other hand, the interaction between different hydrocarbons, as well as other molecules that have been added to the diesel oil in the refinery, may cause various interactions with the molecules present. Therefore, the abnormal behaviors of the extinction coefficient and imaginary permittivity are due to the chemical activity of the mixed liquids. This, in turn, may explain the abnormal dependence of the imaginary excess permittivity on the kerosene fill fraction.
For practical field measurements, the most important task is to detect adulterated diesel oil. Hence, accurate knowledge of the kerosene fill fraction is a secondary issue. Since the origin of the crude oil, its refining to obtain diesel oil, and the additives used may result in chemical activity with kerosene, it is difficult to find a general rule for predicting the kerosene fill fraction, and the present study gives only a partial solution. Both experimental and theoretical investigations are needed to find a rule to predict the kerosene percentage from the imaginary excess permittivity.
In our recent article, we have shown that conventional kerosene and diesel oils from different origins show optical properties such as excess refractive index and excess permittivity [7]. However, excess permittivity, being a more fundamental quantity, is usually preferable. Furthermore, the wavelength dependence of imaginary excess permittivity can be used to estimate, for instance, the contact angle of binary mixtures, which is important in the development of a handheld sensor [6] for screening adulterated diesel oil.
In the case that markers are used, such as polar markers or impurities, for the identification of adulterated diesel oils, as suggested in [24], these markers contribute to the permittivity of diesel oil and kerosene and thus to the excess permittivity. Hence, the excess permittivity is a useful and relatively accessible measure for screening adulterated diesel oil containing fingerprints in the presence of polar markers.
The volume of liquid needed for the measurement of the liquid fuel is small: only a drop is required for measurement with an Abbe refractometer or the volume of a thin spectroscopic cuvette is required. In principle, it is possible to combine a refractometer and a spectrophotometer in a single measurement unit, as described in [25], and utilize the concept of the imaginary excess permittivity for the data. In such a case, only one liquid measurement compartment is required.

Conclusion
In this work, we have studied the wavelength dependence of the imaginary optical properties of diesel oils adulterated with different volume fractions of kerosene in the NIR spectral range using general dispersion analysis based on Kramers-Kronig analysis. We have demonstrated that using both the extinction coefficient and imaginary excess permittivity, it is possible to obtain extra information about diesel oils adulterated with different volume fractions of kerosene, particularly in situations, where the refractive indices of diesel oil and kerosene are very close for a fixed measurement wavelength of the probe light. In addition, when the volume of kerosene was 15% which is below the typical adulteration level of 20%, a strong peak of imaginary excess permittivity was observed, which was very different from the cases of 5% and 10% adulteration. Hence, the extinction coefficient and excess imaginary permittivity complement each other in the screening of diesel oil for adulteration. Moreover, the imaginary permittivity can also be exploited for the development of practical sensors for the screening of diesel oil adulteration.