Next Article in Journal
A Novel Arc Fault Detector for Early Detection of Electrical Fires
Next Article in Special Issue
Low-Loss Coupling of Quantum Cascade Lasers into Hollow-Core Waveguides with Single-Mode Output in the 3.7–7.6 μm Spectral Range
Previous Article in Journal
Mechanical Coupling Error Suppression Technology for an Improved Decoupled Dual-Mass Micro-Gyroscope
Previous Article in Special Issue
Microfluidic Devices for Terahertz Spectroscopy of Live Cells Toward Lab-on-a-Chip Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Essential Limitations of the Standard THz TDS Method for Substance Detection and Identification and a Way of Overcoming Them

by
Vyacheslav A. Trofimov
* and
Svetlana A. Varentsova
Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, Leninskiye Gory, Moscow 119992, Russia
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(4), 502; https://doi.org/10.3390/s16040502
Submission received: 31 January 2016 / Revised: 26 March 2016 / Accepted: 29 March 2016 / Published: 8 April 2016
(This article belongs to the Special Issue Infrared and THz Sensing and Imaging)

Abstract

:
Low efficiency of the standard THz TDS method of the detection and identification of substances based on a comparison of the spectrum for the signal under investigation with a standard signal spectrum is demonstrated using the physical experiments conducted under real conditions with a thick paper bag as well as with Si-based semiconductors under laboratory conditions. In fact, standard THz spectroscopy leads to false detection of hazardous substances in neutral samples, which do not contain them. This disadvantage of the THz TDS method can be overcome by using time-dependent THz pulse spectrum analysis. For a quality assessment of the standard substance spectral features presence in the signal under analysis, one may use time-dependent integral correlation criteria.

Graphical Abstract

1. Introduction

As it is well-known, a detection and identification of explosives, drugs and other hazardous chemical and biological agents is an urgent security problem. For this problem solving, THz radiation has been widely used during the past twenty years [1,2,3,4,5,6,7,8,9,10,11,12,13,14]. Usually, the detection and identification of substances is based on a comparison of the absorption frequencies of the substance under consideration with the absorption frequencies of a standard substance from a database [5,6,7,8,9,10].
We will call this method the standard THz TDS method. It should be stressed that by standard THz TDS method we do not refer to the way for obtaining and registering the THz signal, but rather to the analysis of the spectral characteristics (Fourier spectra, absorbance or reflectance spectra) of the THz signal reflected from or transmitted through the object. Obviously, this method has well-known disadvantages. For example, many explosives have simulants—ordinary substances with similar absorption frequencies—which makes use of the TDS method not effective enough. Opaque packaging, inhomogeneity of the substance surface, high atmospheric humidity during the measurements also lead to low efficiency of this method [11,12,13,14]. We notice that in [10] the author reported the same disadvantages of THz spectroscopy.
In the present paper, we show the principal limitations of the standard THz TDS method for the detection and identification of substances under real conditions—at long distances of about 3.5 m under room conditions with a relative humidity of about 50%—which differ from the disadvantages mentioned above. For this purpose, we conducted a physical experiment with a neutral substance—a thick paper bag. Another example is the identification of Si-based semiconductors under laboratory conditions—at a short distance of about 30 cm and in dry air with a relative humidity less than 2%. We show the standard THz TDS method determines the spectral features of hazardous substances in the neutral substances in both physical experiment cases. This property is not inherent to any particular installation (for example, to that in the Lomonosov Moscow State University), because all our data were obtained from different laboratories. This fact proves the generality of the phenomenon and results in difficulty for construction a device based on the use of the standard THz TDS method, which can really detect dangerous substances in real conditions, as it will lead to a large number of false positives. A few words about why we chose to investigate semiconductors: currently, almost everyone has in their pocket or purse a cell phone or tablet PC that contains semiconductors. If false detection of explosives instead of cell phones occurs many times, the remote screening efficiency will tend to zero.
To overcome these negative features of THz TDS, we have developed a novel method for the detection and identification of substances by using the broadband THz signal with a few cycles. Its main feature consists in the spectral dynamics analysis of the substance response to THz radiation. We call it the SDA-method.
Earlier, the SDA-method was successfully applied in the transmission mode for the identification of various neutral and dangerous substances [15,16,17]. Later we showed the possibility of applying the SDA-method for this purpose in reflection mode [18]. In [19,20] the integral correlation and likeness criteria of probability assessment for the detection and identification of substance were proposed. In [21,22,23,24,25] the spectral properties of THz pulses measured at long distances of about 3.5 m were investigated by means of proposed modified integral criteria.
As our practice shows, to detect the presence of an explosive in a sample is a simpler task than proving its absence. Let us note that we have demonstrated the presence of explosives in the samples by means of the SDA-method and integral correlation criteria in a number of our previous papers, for example, in [15,20,24,25,26]. In particular, in [26] we showed the possibility to detect the explosive RDX in the plastic explosive PWM C4 with an inhomogeneous (rough and concave) surface in reflection mode using the integral correlation criteria (ICC).
By analogy with the “electronic nose”, the method discussed in [19,20,21,22,23,24,25,26] was called by us a “terahertz nose”. We recall that the “electronic nose” [27,28,29] is a multi-sensor system for the detection and analysis of multi-component gas mixtures. Identification in the modern sensor system occurs on the basis of chemical or physical sensor properties’ changes (e.g., a change in conductivity, the change of mass, fluorescence, etc.). Disadvantages of this technology are well-known and are associated with the identification of hidden objects, for example, if packed in polyethylene or paper, or hidden under clothes.
We emphasize that the SDA-method is based on processing of the reflected or transmitted THz signal. The method does not depend on the particular installation which was used for the THz signal measurement. In our opinion, it is a big advantage of our method and that is why we use in our article THz signals obtained in different laboratories. Let us stress once again, that the main aim of this paper is to demonstrate the principal limitations of the standard THz TDS method for the detection and identification of substances.
It should also be noted that besides the security screening applications, the demonstrated ultrafast THz spectroscopy approach also provides new opportunities to study technologically-related materials ranging from nanostructures to strongly correlated electron materials [30,31].
The measurements of the THz signal from thick paper bag at long distance were made in Lomonosov Moscow State University (Moscow, Russia). THz signals transmitted through RDX were registered in the Center for Terahertz Research, Rensselaer Polytechnic Institute (New York, NY, USA) and were provided to us by Xi-Cheng Zhang. The examples with the thick bag and semiconductors were previously partially discussed in [24,25].

2. Description of the Setup and Measurements

In the physical experiment with a thick paper bag, we exploit an EKSPLA commercially available THz spectrometer developed by Teravil Inc. (Vilnius, Lithuania). It is made up of a femtosecond fiber laser, which generates laser pulses with an average power 1 W at 1030 nm centre wavelength, with 75 MHz repetition rate, and its pulse duration is equal to 80 fs. Low-temperature grown GaBiAs is used as photoconductor. The spectral range of the spectrometer is 0.1–5.0 THz, signal-to-noise is better than 103:1 (at 2 THz), 105:1 (at 1 THz) and 106:1 (at 0.4 THz). We use a parabolic mirror for the THz beam focusing on the object under investigation. Because the femtosecond fiber laser has low average power and the laser beam splits many times, we use additional flat mirror behind the object. Therefore, our setup operates in reflection-transmission mode simultaneously for each measurement.
The distance between the parabolic mirror and the object was about to 3.5 m. The principal scheme of the system is shown in Figure 1 for the transmission and reflection mode. A photo of the experimental setup for measurements at long distance can be found in [26].
A THz signal transmitted through the explosive RDX in transmission mode was measured using a TPS 3000 unit (Teraview, Cambridge, UK) in the standard configuration [26]. The main parameters of the system are: spectral range 0.06–3.6 THz, signal-to-noise better than 4000:1, dynamic range higher than 3 OD in the range 2 to 100 cm−1, spectral resolution 0.06 THz and rapid scan mode with 30 scans/s. Measurements were made at the distance of not more than 30 cm between a sample and the mirrors, as described in detail in [9]; the temperature in both cases was 22 °C and the only difference is that in our case the RDX signal was measured in air with a relative humidity of about 50%.

3. Spectral Analysis

3.1. Spectral Properties of the Thick Bag Sample

In Figure 2a, we present the sample investigated under real conditions at a long distance of about 3.5 m; it is a thick paper bag with a thickness of 5–7 mm. The measurement was carried out at a temperature of 18 °C and relative humidity of about 50%. We shall call this THz signal the Thick Bag signal.
Figure 2b shows the Thick Bag signal registered at the time interval t = [0, 110] ps. The pronounced sub-pulses, which are typical for reflected THz signals and caused by multiple reflections from the inner boundaries of the sample, are absent in (b). At the same time, high noise is observed in the signal. It should be stressed that the thick paper sample consists essentially of cellulose. In Figure 2c the chemical formula of cellulose is presented [32]. Note that reference [33] reports the absorption frequencies ν = 2.15 THz and ν = 3.03 THz for pure microcrystalline cellulose samples measured under laboratory conditions. In Figure 2d we show the formula of the explosive RDX [34], as we compare below the spectral properties of the thick paper bag and RDX.
The commonly used TDS deals with the spectrum of the main pulse of a THs signal transmitted or reflected from the object. As one can see from Figure 2b, the main pulse is observed in the time interval t = [0, 25] ps duration t = 25 ps. In order to get its spectrum and to avoid the influence of noise we cut the main pulse out from the long THz signal measured in this time interval.
Figure 3 shows the Fourier spectrum of the main pulse for the Thick Bag signal in the frequency ranges ν = [0, 1.5] THz (a) and ν = [1.5, 3.2] THz (b) calculated in the time interval t = [0, 25] ps with frequency resolution Δν = 0.04 THz. In Figure 3c,d the sample absorbance A is calculated in the same frequency ranges. Here:
A = log 10 ( | P ( ν ) | / | P R E F ( ν ) | )
where | P ( ν ) | , | P R E F ( ν ) | are the absolute spectral amplitude values of the measured and reference signals, respectively [35]. The Fourier spectrum of the whole signal can be found in [24]. In [24] we showed that the minima corresponding to frequencies ν = 0.56, 0.76 THz in Figure 3a may be caused by water vapor contained in the air and the superposition of reflected and transmitted signals. The absence of absorbance maxima at these frequencies in Figure 3c confirms our conclusion.
It is obvious that hazardous substances are absent in the thick paper bag sample. However, in the Fourier spectrum and absorbance of the signal measured under real conditions, it is possible to observe the spectral features of many dangerous substances. According to [6,7,9], the characteristic absorption frequencies of the explosives RDX, HMX, and PETN are: ν = 0.82, 1.05, 1.36, 1.54, 1.95, 2.19, 3.0 THz for RDX; ν = 1.78, 2.51, 2.82 THz for HMX; ν = 2.0, 2.16, 2.84 THz for PETN. In Figure 3, one can see minima in Figure 3a,b and maxima in Figure 3c,d, corresponding to the absorption frequencies of the explosives RDX, HMX, PETN and the illicit drug MDMA [24].
Some of these frequencies ν = 2.16, 2.88 THz are close to those of cellulose ν = 2.15 THz and ν = 3.03 THz [33], the others may be influenced by agents used during the manufacture and preserved in the samples.
This example shows that the main problem of identifying substances with the standard THz TDS method is not that the method determines the frequencies of these substance simulants, but in that the method is not able to show the absence of hazardous substances in the sample.

3.2. Spectral Properties of Si-Based Semiconductors

In the next example, we show that the application of the standard THz TDS method for the identification of substances can lead to inaccurate results even in the case of measurements under laboratory conditions—over a short distance of 30–40 cm, and at room temperature and low humidity.
Below we investigated the spectral properties of n-Si, p-Si semiconductors and silicon wafer with resistivity 40 Om·cm and measured with 12 mm aperture, in comparison with spectral properties of dangerous substances. We will call these THz signals n-Si, p-Si and Si-40-12 for brevity. Measurements of n-Si, p-Si signals were performed at Capital Normal University (Beijing, China) and of Si-40-12 signal at South China Normal University (Guangzhou, China).
The measurements of n-Si, p-Si signals were performed using a Ti:sapphire regenerative amplifier delivering ultrashort optical pulses with a duration of 100 fs and a central wavelength of 800 nm at a repetition rate of 1 kHz. The output of the laser was divided by beam splitters into three portions: a terahertz generation beam, probe beam, and pump beam. The generation pulse was incident on a 2-mm-thick (110) ZnTe emitter to generate the terahertz pulse via optical rectification [36]. The terahertz wave was normally incident to the sample. The transmitted THz signal was detected by free-space electro-optic sampling in a 1-mm-thick <110> ZnTe crystal with the sampling pulse [37]. Then the signal was collected by a lock-in amplifier with its frequency locked to an optical chopper. The system was purged by nitrogen to prevent absorptions by atmospheric humidity. The experiments were performed at room temperature.
Figure 4 shows the Fourier spectra of the main pulses of the n-Si and p-Si signals in the frequency range ν = [0.0, 3.0] THz (a) and ν = [3.0, 4.0] THz (b), in Figure 4c—absorbance in the frequency range ν = [0.0, 3.0] THz. One can see minima in Figure 4a and maxima in Figure 4c corresponding to the absorption frequencies of the explosives RDX, HMX, PETN [9].
In Figure 4d the absorbance of pure monocrystalline silicon is shown, using data from a THz database [38] (Teraphotonics Laboratory, RIKEN, Sendai, Japan). The paper [39] reported that the frequency ν = 3.6 THz is the well-resolved absorption feature of high-resistivity, float-zone silicon. Note that in (d) one can see maxima at frequencies ν = 2.24, 2.65, 3.61 THz. The first one is the same as the extremes in (a,c), marked by “1”. The second maximum in (d) coincides with the maximum n-Si absorbance at ν = 2.65 THz (c). For the p-Si signal the corresponding maximum in (c) is shifted to ν = 2.7 THz. The third maximum in (d) at a frequency ν = 3.61 THz is close to the minimum of the n-Si signal at frequency ν = 3.65 THz in (b) and close to the pure monocrystalline silicon absorption frequency ν = 3.6 THz [39]. At the same time, for the p-Si signal the corresponding minimum in (b) is shifted to ν = 3.7 THz. That is, the semiconductor samples n-Si and p-Si demonstrate spectral features, which are close to those of both of pure silicon and explosives.
The Si-40-12 signal was measured also in transmission mode, at a short distance of about 30–40 cm from the receiver at room temperature in the open air with non-zero humidity using the standard THz TDS unit. The measurement procedure is similar to that described in [40].
In Figure 5 the Fourier spectrum (a,b) and absorbance (c,d) of the main pulse of the Si-40-12 signal is presented in the frequency ranges ν = [0, 1.5] THz (a,c); [1.5, 3.7] THz (b,d). The Fourier spectrum and absorbance of the main pulse were calculated in the time interval t = [4, 24] ps length T = 20 ps. Note that in (b) one can see minima at frequencies ν = 2.25, 2.7, 3.55 THz, which are close to the absorption frequencies of pure silicon ν = 2.24, 2.65, 3.61 THz [38], ν = 3.6 THz [39]. The corresponding maxima in (d) are ν = 2.2, 2.7, 3.55 THz. That is, the silicon-based semiconductor samples from different laboratories (n-Si, p-Si and Si-40-12) demonstrate spectral features which are close to those of pure silicon.
One can see in Figure 5a,b minima and in Figure 5c,d maxima corresponding to the absorption frequencies of the explosives RDX, HMX and PETN. Note that in Figure 5 in the spectrum and absorbance there are more extremes that in Figure 4. This may be caused by the influence of the environment, including water vapor in the air. As in the previous case with n-Si and p-Si semiconductors, the Si-40-12 sample also demonstrates spectral features of explosives.
Therefore, the standard THz TDS method has fundamental limitations which do not allow its use for the detection and identification of substances under real conditions. However, integral correlation criteria together with the SDA-method allow one to show the absence of explosives in the samples of a thick paper bag and Si-based semiconductors.

4. SDA-Method and Modified Integral Correlation Criteria

4.1. Advantages of the Spectral Dynamics Analysis Method (SDA-Method)

It is obvious that the standard THz TDS method does not take into account the changes in instantaneous spectral intensities; it gives information about the spectrum averaged over the time of pulse registration. At the same time, the response of the medium to the action of the THz pulse with a few cycles is essentially non-stationary. The analysis of the evolution of the spectral intensities in time (spectral dynamics) at the chosen frequency ν allows one to get much more information about the substance than the analysis of the spectrum alone.
For the transmitted THz signals obtained under laboratory conditions, it is often sufficient to use the TDS method for identifying substances. However, in actual use the TDS method is inefficient because the THz signals detected under real conditions often have noisy Fourier spectra, which may be distorted by water vapor, the influence of the packaging material, and so on. In this case, it becomes necessary to develop new criteria that are not based on a comparison of the absorption coefficients for evaluating the detection and identification of substances. These effective criteria can be integral correlation and likeness criteria proposed by us on the basis of the application of SDA-method [19,20].
In order to get information about the characteristic absorption frequencies and relaxation times, we analyze the dynamics of the spectral intensities of the signal under investigation. For this purpose, we move the spectral dynamics of standard THz signal at characteristic absorption frequency along the spectral dynamics of the investigated signal in the chosen time intervals. As a standard signal, we use a THz signal transmitted through the sample with the desired substance and measured in laboratory conditions, or a signal measured in ambient air in view of water absorption frequencies. Analyzing the integral correlation between these spectral dynamics, we can conclude about the presence or absence of characteristic features of the standard substance in the sample.
Note that the frequency resolution is obviously dependent on the length of time interval T for which the measurement are performed. As a rule, the duration T was approximately 100 ps. Consequently, the minimum frequency difference Δ ν , which can be resolved by computer processing, is 10 GHz.

4.2. SDA-Method

In Section 4.2 and Section 4.3 we recall the basic steps in constructing the spectral line dynamics in the SDA-method and in calculating the integral correlation criteria, see for example, [23,24,25,26].
The proposed integral criteria are based on the analysis of total correlation characteristics over relatively short time intervals and taking into account the spectral brightness of frequencies ν 1 and ν 2 during these time intervals. Here ν 1 is a chosen frequency of the signal under investigation and ν 2 is a known absorption frequency of the standard signal. The correlation characteristics in turn, are based on dynamics of spectral amplitudes of THz signal under investigation at the chosen frequencies.
Let E ( t ) be the transmitted or reflected THz signal measured in the time interval [tb, te]. Information about the time evolution of the full spectrum or part of it can be obtained by using a time window with duration (length) T that slides along the signal. At each step, the time window is shifted on the chosen time interval Δ , and then the Fourier transform is applied to the function E ( t ) in this window. To avoid any “spreading” of the spectrum, we multiply the signal E ( t ) by function g(t), which tends to zero very quickly at the ends of the window. In order to construct the dynamics of the spectral line (or evolution of the modulus of spectral amplitude) P ( ν , t j ) of the function E ( t ) at chosen frequency ν , the Fourier-Gabor transform is carried out in each time window of length T:
P ( ν , t j ) = 1 T t j t j + T E ( t ) g ( t ) e i 2 π ν ( t t j ) d t ,   g ( t ) = e ( t t c 0.5 T ) k
where t j is a time of window beginning, j is a serial number of window, T is its length, ν is a frequency. The units of t j , T, Δ and ν are ps and THz, respectively. Then we calculate the absolute value of spectral amplitude P ( ν , t j ) in each time interval as its value at the end of the window in order to align the beginning of the physical pulse and its representation in the SDA method:
| P ν ( t j ) | = | P ( ν , t j + T ) |
In this paper, we choose parameters of window length, its shift and power k in the following way: the window length T = 2.8 ps, the window shift Δ = 0.2 ps and k = 20.

4.3. Integral Correlation Criteria

In order to investigate the integral correlation between the dynamics of spectral lines for the reflected or transmitted signal S ( t ) and the standard transmitted signal s ( t ) at chosen frequencies, we introduce the following notation: we denote the discrete set of absolute values of spectral amplitudes of the standard transmitted signal s ( t ) at the chosen frequency ν (see Equations (2) and (3)) as p ν = { | p ν ( t m ) | } , m = 1 , ... , M 1 . The corresponding set of absolute values of spectral amplitudes of the long reflected (or transmitted) THz signal S ( t ) at frequency ν is denoted as P ν = { | P ν ( t m ) | } , m = 1 , ... , M 2 , and its part with M 1 components, which begins at point t n , as P ν ( n ) = { | P ν ( n ) ( t n + m ) | } . Here M 1 and M 2 are the numbers of points in the corresponding dynamics; they depend on the dynamics construction parameters—the window length T and its shift Δ, see Equations (2) and (3). Both sets p ν = { | p ν ( t m ) | } and P ν ( n ) = { | P ν ( n ) ( t n + m ) | } must be averaged at each step t n to avoid the influence of constant components of sets p ν and P ν ( n ) on correlation. Moving the set p ν 1 along the set P ν 2 , we get in each point t n the correlation coefficient:
c p , P ( t n ) = m = 0 M 1 1 ( | p ν 1 ( t m ) | p ν 1 ¯ ) ( | P ν 2 ( t m + n ) | P ν 2 ¯ ) / | | p ν 1 p ν 1 ¯ | | | | P ν 2 ( n ) P ν 2 ¯ | |
where p ν 1 ¯ = m = 0 M 1 1 | p ν 1 ( t m ) | / M 1 , P ν 2 ¯ = m = 0 M 1 1 | P ν 2 ( t m + n ) | / M 1 .
Then using the correlation coefficient c p , P ( t n ) for two spectral dynamics, we consider the following integral characteristic for the detection and identification problem [20]:
C p , P ( t n ) = m = 0 n | c p , P ( t m ) | ,   n = 0 , ... , M 2 M 1
In the present paper, we use the modified criterion on the base of Equation (5), which takes into account the spectral brightness of each frequency ν 1 and ν 2 during the interval of correlation:
C W p , P ( t n ) = m = 0 n | c p , P ( t m ) | w 1 w 2 ,   n = 0 , ... , M 2 M 1
where w 1 = w ( | P ( ν 1 ) | ) , w 2 = w ( | P ( ν 2 ) | ) are the weight coefficients during the interval of correlation.
Along with the criterion (5) we deal with another criterion, in which the sets p ν 1 2 = { | p ν 1 ( t m ) | 2 } ,   P ν 2 2 = { | P ν 2 ( t m ) | 2 } will be used instead of sets p ν 1 and P ν 2 :
C W p , P S Q ( t n ) = m = 0 n | c p 2 , P 2 ( t m ) | w 1 2 w 2 2 ,   n = 0 , ... , M 2 M 1
It should be stressed that if w 1 = 1 and w 2 = 1 , we get integral characteristic Equation (5).
In [20] a criterion was introduced that assesses the similarity (or likeness) of two spectral line dynamics:
L p , P ( t n ) = m = 0 n l p , P ( t m ) ,   n = 0 , ... , M 2 M 1
where:
l p , P ( t n ) = 1 | | ( p ν 1 p ν 1 ¯ ) N ( P ν 2 ( n ) P ν 2 ¯ ) N | | | | ( p ν 1 p ν 1 ¯ ) N | | + | | ( P ν 2 ( n ) P ν 2 ¯ ) N | | ,   n = 0 , ... , M 2 M 1
The subscript N indicates that corresponding sets in Equation (9) are to be normalized, for example, in L 2 . In this paper, we give an example of application of this criterion for identification.

5. Identification Based on Integral Correlation Criteria

5.1. Detecting the Absence of RDX in the Thick Bag Sample

In order to show the absence of explosive RDX in the Thick Bag sample and Si-based semiconductors, we will use the THz signal transmitted through the tablet containing 10% RDX and 90% PE in the ambient air as a standard one (RDX_Air signal). Figure 6 shows the Fourier spectrum of the RDX_Air signal Figure 6a and absorbance in the frequency range ν = [0.6, 3.2] THz. One can see in Figure 6a minima at frequencies ν = 1.15, 1.4, 1.68 THz. At the same time, the maxima of absorbance at these frequencies are absent in Figure 6b, i.e., these minima are false absorption frequencies in the Fourier spectrum of the RDX_Air signal. Their appearance is caused by the presence of strong absorption of signal energy by water vapor. However, the minima of the Fourier spectrum of the RDX_Air signal at frequencies ν = 0.82, 1.95, 2.2, 2.42, 3.0 THz (a) coincide with maxima of the absorbance (b), and they can be used to identify the explosive RDX. Note that these frequencies are in a good agreement with the RDX absorption frequencies given in [6,7,9].
In the spectrum of the Thick Bag signal (a) there is no minimum at frequency, which is equal or close to the characteristic minimum of the RDX_Air spectrum at ν = 2.42 THz. In addition the minimum at ν = 3.0 THz is also absent in this spectrum, so in order to show the absence of RDX in the Thick Bag sample with the help of modified criteria (6) and (7), we will use the dynamics of the spectral lines of RDX_Air signal at frequencies ν = 2.42, 3.0 THz as the standard ones. The corresponding spectral lines dynamics | P ν ( t ) | for the RDX_Air signal are depicted in Figure 6c,d, and for the Thick Bag signal—in Figure 6e,f. Note that each dynamics have their own individual shape.
Here and below the frequency ν for the signal under investigation is detected if the corresponding characteristic C W p , P ( t n ) ( C p , P ( t n ) , L p , P ( t n ) or C W p , P S Q ( t n ) ) calculated for the pair (ν, ν1) lies above all other characteristics in the frequency detection range (FDR). Here ν1 is the chosen absorption frequency of the standard signal. As a rule, the boundaries of the FDR are extremes of the spectrum closest to the frequency under investigation. Vice versa, the frequency ν is not detected if there is at least one other characteristic lying above it in this range of frequencies.
It should be noted that earlier we investigated the influence of the various time intervals chosen for the calculations of the integral correlation criteria (ICC), on the detection and identification of substances. As usual, the main pulse of a THz signal is used for this purpose in transmission mode. In [21] the ICC were applied for the detection and identification of amphetamine-type stimulants using the main pulse of a transmitted THz signal; the measurements were carried out under laboratory conditions—at the short distance of 30 cm from the receiver and low humidity (less than 2%). In [24] we applied ICC for the investigation of the spectral properties of chocolate under real conditions—at a distance of about 3.5 m and high relative humidity of about 50%. It was shown that detection and identification are possible not only on the time interval containing the main pulse but also on the long time interval which follows it. As for the reflected THz signals, in many cases, the main pulse of the reflected signal does not allow one to carry out the identification. In [26] we showed the possibility to detect plastic explosive PWM C4 with inhomogeneous surface using not only the first sub-pulse following the main pulse of the reflected THz signal but the temporal response of the medium at longer intervals, which do not contain both the main pulse and the first sub-pulse. Since in the present paper we discuss the shortcomings of the standard THz TDS method, here and in the following examples, we apply the ICC on the time intervals, which contain only the main pulses of the measured THz signals.
In Figure 7 the integral characteristics CWp,P(tn) are shown in the time interval, containing the main pulse of the Thick Bag THz signal, calculated for the frequencies ν = 2.42 THz (a), 3.0 THz (b) for the Thick Bag signal with RDX_Air as a standard one. In both cases these frequencies are not detected in the ranges ν = [2.36, 2.52] THz (a) and ν = [2.96, 3.04] THz (b). Therefore, frequencies ν = 2.42, 3.0 THz (a,b) are not absorption frequencies of RDX in the Thick Bag signal and explosive RDX is absent in the sample with thick paper bag. In the same way, it is possible to show the absence of explosives HMX, PETN and illicit drug MDMA in this sample.
Note that integral characteristics C W p , P S Q ( t n ) give the same result as CWp,P(tn) but increase the contrast of detection of the corresponding pairs. They can be used in the case of close or coinciding lines for different pairs of frequencies, when detection by characteristics CWp,P(tn) alone is difficult, see [21].

5.2. Detecting Paper in the Sample

Below we detect paper in the Thick Bag sample by means of integral correlation Criterias (6) and (7). For this purpose we will use the standard transmitted THz signal Paper_phase(+80.68), which was measured at short distance of about 30 cm from a receiver at room temperature at South China Normal University (Guangzhou, China). The measurement procedure is similar to that described in [40].
Figure 8 shows the corresponding Fourier spectrum of the main pulse of the transmitted Paper_phase(+80.68) THz signal in the frequency ranges ν = [0, 2.0] THz (a) and ν = [1.6, 3.0] THz (b). Comparing the Fourier spectrum of the standard signal (Figure 8) with the spectrum of the Thick Bag signal (Figure 3), one can see the common or close minima at frequencies ν = 0.56, 0.76, 2.2, 2.88 THz. To find paper in the Thick Bag sample, we will use the spectral lines dynamics of the standard transmitted Paper_phase(+80.68) signal at frequencies ν = 2.16, 2.88 THz which are close to the absorption frequencies of cellulose ν = 2.15, 3.03 THz [29]. The shift of absorption frequencies of the Paper_phase(+80.68) signal may be caused by influence of additional substances in the paper (pigments, fillers, etc.).
In Figure 9 the integral characteristics CWp,P(tn) are presented calculated for frequencies ν = 2.2 THz (a), 2.88 THz (b) for the Thick Bag signal with the Paper_Phase(+80.68) signal as a standard one. Both frequencies in (a), (b) are detected in the frequency ranges ν = [2.0, 2.24] THz (a), ν = [2.84, 3.08] THz (b). That is, spectral features of the standard Paper_Phase(+80.68) signal at ν = 2.16, 2.88 THz are detected at frequencies ν = 2.2, 2.88 THz in the Thick Bag signal. Therefore, paper is found in the Thick Bag sample.

5.3. Detecting RDX Absence in the Si-Based Semiconductors

The n-Si and p-Si spectra in Figure 4a do not contain minima at the frequencies which are equal or close to the main characteristic absorption frequency of the RDX_Air signal at ν = 0.82 THz. The minimum at ν = 1.95 THz is also absent in these spectra, so in order to confirm the absence of RDX in the samples with semiconductors n-Si, p-Si with the help of modified Criteria (6) and (7), we can use the spectral line dynamics of RDX_Air signal at frequencies ν = 0.82, 1.95 THz as standard spectral lines dynamics. In Figure 10 we show the dynamics of spectral lines of the standard RDX_Air signal at frequencies ν = 0.82 (a), 1.95 (b) THz, and in Figure 10c,d—the corresponding dynamics of the n-Si signal.
Figure 11 shows the integral characteristics CWp,P(tn) calculated for the frequencies ν = 0.82 THz (a), 1.95 THz (b) THz for the n-Si signal by using RDX_Air as a standard one. In both cases, these frequencies are not detected—in each frequency range there are characteristics which lie higher than the corresponding values for the standard frequency. We obtain the same result for the p-Si signal.
Integral Criterias (5) and (8) allow one to show the absence of integral correlation even in the case, when the spectrum of the THz signal under investigation contains the minimum at frequency close or equal to characteristic absorption frequency of the standard signal. Therefore, the frequency ν = 2.25 THz is a minimum of the n-Si signal spectrum (Figure 4a), close to absorption frequency ν = 2.2 THz of the RDX_Air (Figure 6).
Figure 12a shows the integral correlation characteristic CWp,P(tn) calculated for the frequency ν = 2.25 THz. It lies above other characteristics in the frequency range ν = [2.15, 2.35] THz. However, using the integral Criteria Cp,P(tn)(5) and Lp,P(tn) (8) one can greatly improve the detection result.
In Figure 12b,c the integral characteristic Equations (5) Cp,P(tn) and characteristic Lp,P(tn) (8) do not detect the frequency ν = 2.25 THz in the frequency range ν = [2.15, 2.35] THz. In (d) the characteristic for ν = 2.25 THz is shown in the reduced time interval t = [10, 20] ps, where one can see another characteristic lying above it. Therefore, integral Criteria (5) and (8) show that the frequency ν = 2.25 THz is not absorption frequency of RDX in the n-Si signal. The same result is obtained for the p-Si signal. In the same way, it is possible to show the absence of the spectral features of HMX and PETN in the n-Si, p-Si signals.
In Figure 13 the integral characteristics CWp,P(tn) are shown for the frequencies ν = 0.82 THz (a), 2.42 THz (b) THz for the Si-40-12 signal by using RDX_Air as a standard one. In both cases, these frequencies are not detected—in each frequency range there are characteristics which lie higher than the corresponding values for the standard frequency. That is, the explosive RDX is absent in the semiconductor Si-40-12. In the same way, it is possible to show the absence of the explosives HMX, PETN in this sample.
We demonstrated that integral Criteria (5) and (8) allow one to show the absence of RDX spectral features even in the case, when the spectrum of the THz signal under investigation contains a minimum at a frequency close or equal to the characteristic absorption frequency of the standard signal. In order to enhance the detection reliability it is necessary to use different types of integral criteria simultaneously.

6. Conclusions

We showed the essential limitations of the standard THz TDS method for the detection and identification of substances, which are based on comparing their stationary spectra. As examples of neutral substances under analysis, we used a thick paper bag and the silicon-based semiconductors n-Si, p-Si and a silicon wafer. We demonstrated that the standard THz TDS method detects the spectral features of the explosives RDX, HMX, PETN and illicit drug MDMA in these samples. At the same time, based only on this method, it is impossible to show the actual absence of hazardous substances in the samples. This fact makes the THz TDS method insufficient for the reliable identification of substances not only under real conditions (at a long distance and high relative humidity), but also under laboratory conditions (at a short distance and low relative humidity). It should be stressed that this also makes it difficult to construct a device that can detect dangerous substances in the real world by means of the standard THz TDS method, as it will lead to a large number of false positives—for example, instead of a cell phone such a device will detect RDX in the pocket of a person.
At the same time, the proposed integral criteria and SDA-method allow us to detect the absence of dangerous substances in the samples under investigation and to find paper in the Thick Paper Bag sample. In order to enhance the detection reliability it is necessary to use the different types of integral criteria simultaneously.
We also showed that the Thick Paper Bag and Paper_phase(+86.68) samples demonstrate spectral features similar to those of pure microcrystalline cellulose, and the silicon-based semiconductor samples from different laboratories (n-Si, p-Si and Si-40-12) demonstrate similar spectral features of pure silicon.
Thus, the discussed method is a promising and competitive one for the effective detection and identification of various substances both under laboratory and real conditions, in comparison with the THz TDS method, based on comparison of the substances’ spectra. The method can be used with success for security problems solving, for non-destructive testing, as well as for quality control in the pharmaceutical and food industry.

Acknowledgments

The investigation was made using support of the Russian Science Foundation (Grant No. 14-21-00081). We are grateful to Xi-Cheng Zhang (Center for Terahertz Research, Rensselaer Polytechnic Institute, NY, USA), Sheng Lan (South China Normal University, Guangzhou, China), Qingli Zhou and Cunlin Zhang (Capital Normal University, Beijing, China), for providing us the data used in the present paper.

Author Contributions

Vyacheslav A. Trofimov provided the idea of the SDA-method and integral correlation criteria for the detection and identification of substances with the THz signal with a few cycles, performed the experiments at Moscow State University, provided data for the paper and wrote the paper. Svetlana A. Varentsova analyzed the data, made all figures, prepared the literature survey and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
THz TDSTerahertz time-domain spectroscopy
SDA-methodSpectral Dynamics Analysis method
FDRFrequency detection range

References

  1. Federici, J.F.; Schulkin, B.; Huang, F.; Gary, D.; Barat, R.; Oliveira, F.; Zimdars, D. THz imaging and sensing for security applications—Explosives, weapons and drugs. Semicond. Sci. Technol. 2005, 20, S266–S280. [Google Scholar] [CrossRef]
  2. Baker, C.; Lo, T.; Tribe, W.R.; Cole, B.E.; Hogbin, M.R.; Kemp, M.C. Detection of concealed explosives at a distance using terahertz technology. IEEE Proc. 2007, 95, 1559–1565. [Google Scholar] [CrossRef]
  3. Kemp, M.C. Screening mail for powders using terahertz technology. Proc. SPIE 2011, 8189. [Google Scholar] [CrossRef]
  4. Beigang, R.; Biedron, S.G.; Dyjak, S.; Ellrich, F.; Magnus, W.; Haakestad, M.W.; Hübsch, D.; Kartaloglu, T.; Ozbay, E.; Ospald, F.; et al. Comparison of terahertz technologies for detection and identification of explosives. Proc. SPIE 2014, 9102. [Google Scholar] [CrossRef]
  5. Katz, G.; Zybin, S.; Goddard, W.A.; Zeiri, Y.; Kosloff, R. Direct MD simulations of terahertz absorption and 2D spectroscopy applied to explosive crystals. J. Phys. Chem. Lett. 2014, 5, 772–776. [Google Scholar] [CrossRef] [PubMed]
  6. Leahy-Hoppa, M.R.; Fitch, M.J.; Zheng, X.; Hayden, L.M.; Osiander, R. Wideband terahertz spectroscopy of explosives. Chem. Phys. Lett. 2007, 434, 227–230. [Google Scholar] [CrossRef]
  7. Davies, A.G.; Burnett, A.D.; Fan, W.; Linfield, E.H.; Cunningham, J.E. Terahertz spectroscopy of explosives and drugs. Mater. Today 2008, 11, 18–26. [Google Scholar] [CrossRef]
  8. Choi, K.; Hong, T.; Sim, K.I.; Ha, T.; Park, B.C.; Chung, J.H.; Cho, S.G.; Kim, J.H. Reflection terahertz time-domain spectroscopy of RDX and HMX explosives. J. Appl. Phys. 2014, 115, 023105. [Google Scholar] [CrossRef]
  9. Chen, J.; Chen, Y.; Zhao, H.; Bastiaans, G.J.; Zhang, X.C. Absorption coefficients of selected explosives and related compounds in the range of 0.1–2.8 THz. Opt. Express 2007, 15, 12060–12067. [Google Scholar] [CrossRef] [PubMed]
  10. Kemp, M.C. Explosives detection by terahertz spectroscopy—A bridge too far? IEEE Trans. Terahertz Sci. Technol. 2011, 1, 282–292. [Google Scholar] [CrossRef]
  11. Xiong, W.; Shen, J. Fingerprint extraction from interference destruction terahertz spectrum. Opt. Express 2010, 18, 21798–21803. [Google Scholar] [CrossRef] [PubMed]
  12. Palka, N. Identification of concealed materials, including explosives, by terahertz reflection spectroscopy. Opt. Eng. 2014, 53. [Google Scholar] [CrossRef]
  13. Ortolani, M.; Lee, J.S.; Schade, U.; Hübers, H.W. Surface roughness effects on the terahertz reflectance of pure explosive materials. Appl. Phys. Lett. 2008, 93. [Google Scholar] [CrossRef]
  14. Withayachumnankul, W.; Fischer, B.M.; Abbott, D. Numerical removal of water vapor effects from terahertz time-domain spectroscopy measurements. Proc. R. Soc. A 2008, 464, 2435–2456. [Google Scholar] [CrossRef]
  15. Trofimov, V.A.; Varentsova, S.A.; Chen, J.; Zhang, X.-C. Identification of explosive media using their spectrum dynamics under the action of THz pulse. Proc. SPIE 2009, 7486. [Google Scholar] [CrossRef]
  16. Trofimov, V.A.; Varentsova, S.A. 2D THz signature for substance identification. Pros. SPIE 2010, 7687. [Google Scholar] [CrossRef]
  17. Trofimov, V.A.; Varentsova, S.A.; Shen, J.; Zhang, C.; Zhou, Q.; Shi, Y. 2D signature for identification of drugs. Proc. SPIE 2011, 8040. [Google Scholar] [CrossRef]
  18. Trofimov, V.A.; Varentsova, S.A.; Palka, N.; Szustakowski, M.; Trzcinski, T. The method of the spectral dynamics analysis of reflected signal for problem of identification of substance. Proc. SPIE 2011, 8023. [Google Scholar] [CrossRef]
  19. Trofimov, V.A.; Peskov, N.V.; Kirillov, D.A. Efficiency of using correlation function for estimation of probability of substance detection on the base of THz spectral dynamics. Proc. SPIE 2012, 8496. [Google Scholar] [CrossRef]
  20. Trofimov, V.A.; Varentsova, S.A.; Szustakowski, M.; Palka, N. Influence of surface of explosive on its detection and identification using the SDA method for analysis of the reflected THz signal. Proc. SPIE 2013, 8734. [Google Scholar] [CrossRef]
  21. Trofimov, V.A.; Varentsova, S.A.; Trofimov, V.V.; Tikhomirov, V.V. Peculiarities of the detection and identification of substance at long distance. Proc. SPIE 2014, 9102. [Google Scholar] [CrossRef]
  22. Trofimov, V.A.; Varentsova, S.A.; Trofimov, V.V. Transmission of THz pulse with a few circles through opaque samples placed at long distance (4–6 m). Proc. SPIE 2014, 9199. [Google Scholar] [CrossRef]
  23. Trofimov, V.A.; Varentsova, S.A. Real-time criteria based on spectral dynamics of medium response for the detection and identification of substance using THz signal. Proc. SPIE 2014, 9253. [Google Scholar] [CrossRef]
  24. Trofimov, V.A.; Varentsova, S.A.; Trofimov, V.V. Possibility of the detection and identification of substance at long distance at using broad THz pulse. Proc. SPIE 2014, 9253. [Google Scholar] [CrossRef]
  25. Trofimov, V.A.; Varentsova, S.A. Efficiency of using the spectral dynamics for pulsed THz spectroscopy of both explosives and other materials. Proc. SPIE 2015, 9454. [Google Scholar] [CrossRef]
  26. Trofimov, V.A.; Varentsova, S.A. An effective method for substance detection using the broad spectrum THz signal: A “terahertz nose”. Sensors 2015, 15, 12103–12132. [Google Scholar] [CrossRef] [PubMed]
  27. Singh, S. Sensors—An effective approach for the detection of explosives. J. Hazard. Mater. 2007, 144, 15–28. [Google Scholar] [CrossRef] [PubMed]
  28. Buryakov, I.A.; Buryakov, T.I.; Matsaev, V.T. Optical chemical sensors for the detection of explosives and associated substances. J. Anal. Chem. 2014, 69, 616–631. [Google Scholar] [CrossRef]
  29. Onodera, T.; Toko, K. Towards an electronic dog nose: Surface plasmon resonance immunosensor for security and safety. Sensors 2014, 14, 16586–16616. [Google Scholar] [CrossRef] [PubMed]
  30. Li, T.; Patz, A.; Mouchliadis, L.; Yan, J.; Lograsso, T.A.; Perakis, I.E.; Wang, J. Femtosecond switching of magnetism via strongly correlated spin-charge quantum excitations. Nature 2013, 496, 69–73. [Google Scholar] [CrossRef] [PubMed]
  31. Luo, L.; Chatzakis, I.; Patz, A.; Wang, J. Ultrafast terahertz probes of interacting dark excitons in chirality-specific semiconducting single-walled carbon nanotubes. Phys. Rev. Lett. 2015, 114. [Google Scholar] [CrossRef] [PubMed]
  32. A-level Applied Science/Colour Chemistry/Fibres/Cellulose. Available online: https://www.en.wikibooks.org/wiki/A-level_Applied_Science/Colour_Chemistry/Fibres/Cellulose) (accessed on 31 March 2016).
  33. Vieira, F.S.; Pasquini, C. Determination of cellulose crystallinity by terahertz-time domain spectroscopy. Anal. Chem. 2014, 86, 3780–3786. [Google Scholar] [CrossRef] [PubMed]
  34. RDX. Available online: https://www.en.wikipedia.org/wiki/RDX (accessed on 31 March 2016).
  35. Kemp, M.C.; Baker, C.; Gregory, I. NATO Security through Science Series; IOS Press Books Online: Brussels, Belgium, 2006. [Google Scholar]
  36. Rice, A.; Jin, Y.; Ma, X.F.; Zhang, X.-C.; Bliss, D.; Larkin, J.; Alexander, M. Terahertz optical rectification from <110> zinc-blende crystals. Appl. Phys. Lett. 1994, 64, 1324–1326. [Google Scholar] [CrossRef]
  37. Wu, Q.; Litz, M.; Zhang, X.-C. Broadband detection capability of ZnTe electro-optic field detectors. Appl. Phys. Lett. 1996, 68, 2924–2926. [Google Scholar] [CrossRef]
  38. THz Database 2.0. Available online: http://www.riken.jp/THzdatabase/ (accessed on 31 March 2016).
  39. Dai, J.; Zhang, J.; Zhang, W.; Grischkowsky, D. Terahertz time-domain spectroscopy characterization of the far-infrared absorption and index of refraction of high-resistivity, float-zone silicon. JOSA B 2004, 21, 1379–1386. [Google Scholar] [CrossRef]
  40. Liu, H.Y.; Dai, Q.F.; Wu, L.J.; Lan, S.; Trofimov, V.A.; Varentsova, S.A. Effects of p-type doping on the optical properties of InAs/GaAs quantum dots. Solid State Commun. 2012, 152, 435–439. [Google Scholar] [CrossRef]
Figure 1. The principal scheme of the system for the transmission (a) and reflection (b) mode.
Figure 1. The principal scheme of the system for the transmission (a) and reflection (b) mode.
Sensors 16 00502 g001
Figure 2. Thick paper bag (a) investigated in real conditions at long distance of about 3.5 m; the Thick Bag signal (b) at the time interval t = [0, 110] ps; chemical formulae of cellulose (c) and RDX (d).
Figure 2. Thick paper bag (a) investigated in real conditions at long distance of about 3.5 m; the Thick Bag signal (b) at the time interval t = [0, 110] ps; chemical formulae of cellulose (c) and RDX (d).
Sensors 16 00502 g002
Figure 3. Fourier spectrum of the main pulse of the Thick Bag signal (a,b) and the absorbance of the sample (c,d) in the frequency ranges ν = [0, 1.5] THz (a,c); ν = [1.5, 3.2] THz (b,d) calculated with frequency resolution Δν = 0.04 THz.
Figure 3. Fourier spectrum of the main pulse of the Thick Bag signal (a,b) and the absorbance of the sample (c,d) in the frequency ranges ν = [0, 1.5] THz (a,c); ν = [1.5, 3.2] THz (b,d) calculated with frequency resolution Δν = 0.04 THz.
Sensors 16 00502 g003
Figure 4. Fourier spectrum (a,b) and absorbance (c) of the n-Si and p-Si THz signals in the frequency ranges ν = [0.0, 3.0] (a,с); [3.0, 4.0] THz (b); absorbance of monocrystalline silicon in the frequency range ν = [0.0, 3.0] THz (d).
Figure 4. Fourier spectrum (a,b) and absorbance (c) of the n-Si and p-Si THz signals in the frequency ranges ν = [0.0, 3.0] (a,с); [3.0, 4.0] THz (b); absorbance of monocrystalline silicon in the frequency range ν = [0.0, 3.0] THz (d).
Sensors 16 00502 g004
Figure 5. Fourier spectrum (a,b) and absorbance (c,d) of the main pulse of the Si-40-12 signal in the frequency ranges ν = [0, 1.6] THz (a,c); [1.5, 3.2] THz (b,d).
Figure 5. Fourier spectrum (a,b) and absorbance (c,d) of the main pulse of the Si-40-12 signal in the frequency ranges ν = [0, 1.6] THz (a,c); [1.5, 3.2] THz (b,d).
Sensors 16 00502 g005
Figure 6. Fourier spectrum of the RDX_Air THz signal (a) and absorbance (b) in the frequency range ν = [0.6, 3.2] THz; dynamics of spectral lines of the RDX_Air (c,d) and Thick Bag (e,f) signals at frequencies ν = 2.42 THz (c,e); 3.0 THz (d,f).
Figure 6. Fourier spectrum of the RDX_Air THz signal (a) and absorbance (b) in the frequency range ν = [0.6, 3.2] THz; dynamics of spectral lines of the RDX_Air (c,d) and Thick Bag (e,f) signals at frequencies ν = 2.42 THz (c,e); 3.0 THz (d,f).
Sensors 16 00502 g006
Figure 7. Time-dependent integral characteristic CWp,P(tn) calculated for the frequencies ν = 2.42 THz (a); 3.0 THz (b) if the modified integral criterion is applied for the Thick Bag signal and RDX_Air as a standard signal.
Figure 7. Time-dependent integral characteristic CWp,P(tn) calculated for the frequencies ν = 2.42 THz (a); 3.0 THz (b) if the modified integral criterion is applied for the Thick Bag signal and RDX_Air as a standard signal.
Sensors 16 00502 g007
Figure 8. Fourier spectrum of the main pulse of the transmitted Paper_phase(+80.68) THz signal in the frequency ranges ν = [0, 2.0] THz (a) and ν = [1.6, 3.0] THz (b).
Figure 8. Fourier spectrum of the main pulse of the transmitted Paper_phase(+80.68) THz signal in the frequency ranges ν = [0, 2.0] THz (a) and ν = [1.6, 3.0] THz (b).
Sensors 16 00502 g008
Figure 9. Time-dependent integral characteristics CWp,P(tn) detection frequencies ν = 2.2 THz (a); 2.88 THz (b) for the Thick Bag signal and Paper_Phase(+80.68) as a standard signal.
Figure 9. Time-dependent integral characteristics CWp,P(tn) detection frequencies ν = 2.2 THz (a); 2.88 THz (b) for the Thick Bag signal and Paper_Phase(+80.68) as a standard signal.
Sensors 16 00502 g009
Figure 10. Dynamics of spectral lines of the RDX_Air (a,b) and n-Si (c,d) signals at frequencies ν = 0.82 THz (a,c), 1.95 THz (b,c).
Figure 10. Dynamics of spectral lines of the RDX_Air (a,b) and n-Si (c,d) signals at frequencies ν = 0.82 THz (a,c), 1.95 THz (b,c).
Sensors 16 00502 g010
Figure 11. Integral characteristics CWp,P(tn) calculated for the frequencies ν = 0.82 THz (a); 1.95 THz (b) THz for n-Si signal with RDX_Air as a standard one.
Figure 11. Integral characteristics CWp,P(tn) calculated for the frequencies ν = 0.82 THz (a); 1.95 THz (b) THz for n-Si signal with RDX_Air as a standard one.
Sensors 16 00502 g011
Figure 12. Integral characteristics CWp,P(tn) (a); Cp,P(tn) (b); Lp,P(tn) (c,d) calculated for the frequency ν = 2.25 THz for n-Si signal with RDX_Air as a standard one.
Figure 12. Integral characteristics CWp,P(tn) (a); Cp,P(tn) (b); Lp,P(tn) (c,d) calculated for the frequency ν = 2.25 THz for n-Si signal with RDX_Air as a standard one.
Sensors 16 00502 g012
Figure 13. Integral characteristics CWp,P(tn) calculated for the frequencies ν = 0.82 THz (a); 2.42 THz (b) THz for Si-40-12 signal with RDX_Air as a standard one.
Figure 13. Integral characteristics CWp,P(tn) calculated for the frequencies ν = 0.82 THz (a); 2.42 THz (b) THz for Si-40-12 signal with RDX_Air as a standard one.
Sensors 16 00502 g013

Share and Cite

MDPI and ACS Style

Trofimov, V.A.; Varentsova, S.A. Essential Limitations of the Standard THz TDS Method for Substance Detection and Identification and a Way of Overcoming Them. Sensors 2016, 16, 502. https://doi.org/10.3390/s16040502

AMA Style

Trofimov VA, Varentsova SA. Essential Limitations of the Standard THz TDS Method for Substance Detection and Identification and a Way of Overcoming Them. Sensors. 2016; 16(4):502. https://doi.org/10.3390/s16040502

Chicago/Turabian Style

Trofimov, Vyacheslav A., and Svetlana A. Varentsova. 2016. "Essential Limitations of the Standard THz TDS Method for Substance Detection and Identification and a Way of Overcoming Them" Sensors 16, no. 4: 502. https://doi.org/10.3390/s16040502

APA Style

Trofimov, V. A., & Varentsova, S. A. (2016). Essential Limitations of the Standard THz TDS Method for Substance Detection and Identification and a Way of Overcoming Them. Sensors, 16(4), 502. https://doi.org/10.3390/s16040502

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop