US20100108888A1 - Reducing scattering related features in terhertz time domain spectroscopy bz averaging the impulse waveform over a plurality of sample points - Google Patents
Reducing scattering related features in terhertz time domain spectroscopy bz averaging the impulse waveform over a plurality of sample points Download PDFInfo
- Publication number
- US20100108888A1 US20100108888A1 US12/613,355 US61335509A US2010108888A1 US 20100108888 A1 US20100108888 A1 US 20100108888A1 US 61335509 A US61335509 A US 61335509A US 2010108888 A1 US2010108888 A1 US 2010108888A1
- Authority
- US
- United States
- Prior art keywords
- sample
- time domain
- radiation
- thz
- domain waveform
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 238000012935 Averaging Methods 0.000 title claims description 9
- 238000004611 spectroscopical analysis Methods 0.000 title description 7
- 230000005855 radiation Effects 0.000 claims abstract description 84
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000001678 irradiating effect Effects 0.000 claims abstract description 16
- 230000005670 electromagnetic radiation Effects 0.000 claims abstract description 11
- 230000003595 spectral effect Effects 0.000 claims description 24
- 230000000694 effects Effects 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 11
- 238000003384 imaging method Methods 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 5
- 239000000523 sample Substances 0.000 description 91
- 229920000573 polyethylene Polymers 0.000 description 25
- 239000004698 Polyethylene Substances 0.000 description 19
- 238000010521 absorption reaction Methods 0.000 description 18
- CZMRCDWAGMRECN-UGDNZRGBSA-N Sucrose Chemical compound O[C@H]1[C@H](O)[C@@H](CO)O[C@@]1(CO)O[C@@H]1[C@H](O)[C@@H](O)[C@H](O)[C@@H](CO)O1 CZMRCDWAGMRECN-UGDNZRGBSA-N 0.000 description 16
- 229930006000 Sucrose Natural products 0.000 description 16
- 230000005540 biological transmission Effects 0.000 description 16
- 239000005720 sucrose Substances 0.000 description 16
- GUBGYTABKSRVRQ-QKKXKWKRSA-N Lactose Natural products OC[C@H]1O[C@@H](O[C@H]2[C@H](O)[C@@H](O)C(O)O[C@@H]2CO)[C@H](O)[C@@H](O)[C@H]1O GUBGYTABKSRVRQ-QKKXKWKRSA-N 0.000 description 13
- 239000008101 lactose Substances 0.000 description 13
- 239000003570 air Substances 0.000 description 10
- 238000001514 detection method Methods 0.000 description 9
- 239000002245 particle Substances 0.000 description 9
- 239000002360 explosive Substances 0.000 description 8
- 239000000463 material Substances 0.000 description 8
- 230000005684 electric field Effects 0.000 description 6
- -1 masonry Substances 0.000 description 6
- GUBGYTABKSRVRQ-XLOQQCSPSA-N Alpha-Lactose Chemical compound O[C@@H]1[C@@H](O)[C@@H](O)[C@@H](CO)O[C@H]1O[C@@H]1[C@@H](CO)O[C@H](O)[C@H](O)[C@H]1O GUBGYTABKSRVRQ-XLOQQCSPSA-N 0.000 description 5
- 238000011835 investigation Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 239000000843 powder Substances 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 230000002123 temporal effect Effects 0.000 description 4
- 229910001218 Gallium arsenide Inorganic materials 0.000 description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical compound O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 239000013074 reference sample Substances 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 239000000123 paper Substances 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000001055 reflectance spectroscopy Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 229910052594 sapphire Inorganic materials 0.000 description 2
- 239000010980 sapphire Substances 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- XTWYTFMLZFPYCI-KQYNXXCUSA-N 5'-adenylphosphoric acid Chemical compound C1=NC=2C(N)=NC=NC=2N1[C@@H]1O[C@H](COP(O)(=O)OP(O)(O)=O)[C@@H](O)[C@H]1O XTWYTFMLZFPYCI-KQYNXXCUSA-N 0.000 description 1
- 229910005540 GaP Inorganic materials 0.000 description 1
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 239000007836 KH2PO4 Substances 0.000 description 1
- 229910003327 LiNbO3 Inorganic materials 0.000 description 1
- 229910012463 LiTaO3 Inorganic materials 0.000 description 1
- 229910007709 ZnTe Inorganic materials 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 229910000147 aluminium phosphate Inorganic materials 0.000 description 1
- LFVGISIMTYGQHF-UHFFFAOYSA-N ammonium dihydrogen phosphate Chemical compound [NH4+].OP(O)([O-])=O LFVGISIMTYGQHF-UHFFFAOYSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 229910002113 barium titanate Inorganic materials 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 229910052980 cadmium sulfide Inorganic materials 0.000 description 1
- UHYPYGJEEGLRJD-UHFFFAOYSA-N cadmium(2+);selenium(2-) Chemical compound [Se-2].[Cd+2] UHYPYGJEEGLRJD-UHFFFAOYSA-N 0.000 description 1
- 239000011111 cardboard Substances 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- CSJLBAMHHLJAAS-UHFFFAOYSA-N diethylaminosulfur trifluoride Substances CCN(CC)S(F)(F)F CSJLBAMHHLJAAS-UHFFFAOYSA-N 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000000116 mitigating effect Effects 0.000 description 1
- 229910000402 monopotassium phosphate Inorganic materials 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 239000008177 pharmaceutical agent Substances 0.000 description 1
- 239000005365 phosphate glass Substances 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- GNSKLFRGEWLPPA-UHFFFAOYSA-M potassium dihydrogen phosphate Chemical compound [K+].OP(O)([O-])=O GNSKLFRGEWLPPA-UHFFFAOYSA-M 0.000 description 1
- 229910052968 proustite Inorganic materials 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- SBIBMFFZSBJNJF-UHFFFAOYSA-N selenium;zinc Chemical compound [Se]=[Zn] SBIBMFFZSBJNJF-UHFFFAOYSA-N 0.000 description 1
- 239000000377 silicon dioxide Substances 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 229910052714 tellurium Inorganic materials 0.000 description 1
- 238000000411 transmission spectrum Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
- XLOMVQKBTHCTTD-UHFFFAOYSA-N zinc oxide Inorganic materials [Zn]=O XLOMVQKBTHCTTD-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
- G01N21/3586—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation by Terahertz time domain spectroscopy [THz-TDS]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
Definitions
- the present invention relates generally to the field of apparatus and methods for imaging and/or investigating samples in the far-infrared (far-IR)/Terahertz (THz) frequency range from 100 GHz to 100 THz.
- the radiation utilised is in the frequency range of 500 GHz to 100 THz and more preferably from 1 THz to 100 THz and most preferably from 700 GHz to 10 THz.
- US Patent Application Number 2001/0033636 discloses a method of detecting explosives in luggage which uses X-rays to determine an average atomic number of a material. The average atomic number is then compared with known average atomic numbers of explosive materials to determine if the material being examined is an explosive. This technique, however is unsuitable for routine-security screening of people as its radiation is ionising.
- THz radiation is suitable for screening people and their possessions as it is non-ionising and can pass through clothing, paper, cardboard, wood, masonry, plastic and ceramics. Therefore, it is safer than x-ray techniques.
- THz radiation could be used to obtain spectral information of explosive materials, as the organic molecular nature of explosives, and their crystalline structure are of the appropriate form. To date, however, this has not been effectively achieved.
- THz systems utilise the transmission of THz radiation through the sample under investigation to a detector.
- samples can be prepared to avoid attenuation effects.
- samples are often granulated which leads to scattering of the irradiating radiation.
- Terahertz pulse spectroscopy is a coherent detection system which means that it has a large dynamic range.
- the effects of scattering or attenuation means that a THz signal may be too weak to detect in certain samples.
- THz systems which utilise the reflection of THz radiation from a sample will not be limited by attenuation effects in the same way that transmission-based systems are.
- reflection-based systems experience different problems which in practice has meant that effective THz reflection spectroscopy is not possible.
- Thz radiation will impinge upon the surface of the object under investigation.
- the surface of the object represents an interface between a medium with a first refractive index (usually air) and the object itself (which has a second refractive index which is different to the first refractive index).
- the difference in the refractive indices means that a proportion of the radiation will be transmitted into the object (for an air/polythene interface this amounts to around 96% of the incident radiation) and a proportion will be reflected (for the air/PE example this therefore equates to around 4% of the incident radiation).
- the reflected beam will always have a weaker intensity than the transmitted beam. This means that the reflected beam is more susceptible to noise such as scattered radiation.
- the arrival times of radiation arriving at the detector from the surface of the object and scattering centres will be different.
- the time of arrival will depend on the depth of the scattering centre from the surface of the object.
- the time of arrival depends on the extra path length the scattered radiation has traveled relative to radiation that has directly passed through the object. This spread of arrival times causes artificial spectral features which can either mask actual features or create false features.
- One technique for mitigating the effects of scattering is to window the radiation detected at the detecting means. By gating the arrival time within which the system will detect photons any late arriving photons (from scattering events) will not be detected. Such a technique however reduces the efficiency of any spectroscopic analysis of the sample in question as it inevitably loses spectral information.
- the present invention provides a method of investigating an object, comprising the steps of
- the invention which is the subject of the present invention is suitable for use in both transmission based systems and also reflection based systems.
- the method comprises irradiating the object under investigation and detecting the radiation reflected/transmitted from a point on the object.
- the irradiation/detection steps are then repeated for a plurality of points such that an averaged time domain waveform can be calculated for the object.
- the method of the present invention acknowledges that signals due to scattering depend strongly on the geometry of the scattering centres within the object whereas directly reflected or transmitted signals do not. By taking repeated measurements the random signal elements introduced by scattering effects are removed and the true reflection/transmission signals, which have no phase dependence on the geometry of the object, remain.
- the intensity of the incident (I o ( ⁇ ) and the transmitted (I t ( ⁇ )) radiation is measured.
- the optical absorption coefficient, ⁇ ( ⁇ ) can then be calculated from the following equation:
- TPS terahertz pulsed spectroscopy
- the transmittance drops significantly from 100% to a few percent. Therefore the transmission measurement is very sensitive to absorption.
- the reflection loss at the air/sample interfaces which is mainly determined by the refractive index difference, only causes a baseline shift in the transmission spectrum.
- the reflectance can be calculated from the measured electric field of THz radiation as:
- steps are taken to remove features in the detected radiation that relate to the system components. Therefore, preferably a reference measurement is taken in the absence of a sample. The signal detected in the absence of a sample can then be subtracted from the signal detected in the presence of a sample.
- the reference measurement can be taken in the presence of a reference sample with well known THz properties.
- a reference sample could be air, water or polyethylene.
- an absorption profile of the object can be determined by transforming the averaged time domain waveform into the frequency domain by means of a Fourier transform.
- the composition of the sample can be derived from this absorption profile.
- any absorption features will produce maxima and minima in the absorption profile.
- the derivative of the absorption profile with respect to frequency can be determined.
- the method according to the first aspect of the present invention further comprises the step of progressively increasing the number of points that are irradiated until the scattering related effects are sufficiently reduced in order to determine the features of the object under investigation.
- the method according to the first aspect of the invention comprises the step of removing that portion of the averaged time domain waveform that corresponds to directly reflected/transmitted radiation and performing a Fourier transform on the remainder of the detected radiation in order to derive compositional information relating to the object.
- the “late time response” is analysed and there is therefore provided a method of investigating an object, comprising the steps of
- an apparatus for investigating an object comprising
- compositional information relating to a sample can be determined by investigating the sample according to the method of the first aspect of the present invention.
- a compositional image of the sample can then be derived by imaging the sample through a spectral filter that is arranged only to let radiation relating to one (or more) of the components of the sample through to the detector. Therefore, in a fourth aspect of the present invention there is provided of imaging a sample comprising the steps of
- An alternative method of reducing the effects of scatter is to use a diffuse irradiating beam of radiation. Accordingly, in a fifth aspect of the present invention there is provided a method of investigating an object, comprising the steps of:
- FIG. 1 illustrates a schematic of a pulsed terahertz reflection investigative technique utilised in a first embodiment of the present invention
- FIG. 2 illustrates a schematic of a pulsed terahertz transmission investigative technique utilised in a second embodiment of the present invention
- FIG. 3 shows a schematic representation of an air/polythene interface showing the various reflection, transmission and scattering components
- FIGS. 4 a and 4 b show time domain waveforms for samples composed of sucrose and PE powder.
- FIG. 5 shows a series of waveform traces measured from a mixture of fine lactose particles and large sucrose particles (500-800 ⁇ m size/diameter).
- FIGS. 6 a and 6 b show absorption plots at a single point ( FIG. 6 a ) and an averaged dataset ( FIG. 6 b ).
- FIGS. 7 a and 7 b show the first derivative of the absorption plots of FIGS. 6 a and 6 b.
- FIGS. 8 a and 8 b show refractive index plots of the sample imaged in FIGS. 6 a and 6 b.
- FIG. 9 shows various THz reflection spectra of the same sample for various sample areas.
- FIG. 10 shows a schematic of a layered PE/lactose/PE sample.
- FIG. 11 shows various temporal waveform traces obtained from the sample shown in FIG. 10
- FIG. 12 shows various spectral plots obtained from the sample shown in FIG. 10 .
- FIG. 13 shows a chemical mapping image of a sucrose/lactose/polyethylene sample
- FIG. 14 a shows a further chemical distribution image of the surface of a sample.
- FIG. 14 b shows temporal waveforms associated with the image shown in FIG. 14 a.
- FIG. 15 shows an imaging configuration for a THz system
- a terahertz pulsed investigating arrangement which comprises an ultra-short pulse laser 10 which may be, for example, Ti:sapphire, Yb:Er doped fibre, Cr:LiSAF, Yb:silica, Nd:YLF, Nd:Glass, Nd:YAG, Alexandrite Yb:Phosphate Glass QX, Yb:GdCOB, Yb:YAG, Yb:KG d(WO 4 ) or Yb:BOYS laser.
- This laser 10 emits pulses of radiation 11 , such as a collimated beam of pulses, each of which comprise a plurality of frequencies.
- the pulses generated by the laser preferably having a pulse duration of less than 200 fs.
- the beam of generated pulses is directed into beam splitter 12 .
- the beam splitter splits the beam into a pump beam 13 , which is used to irradiate the sample, and a probe beam 14 , which is used during detection.
- Scanning delay line 16 is a variable optical delay, which in its simplest form comprises two mirrors that serve to reflect the beam through a 180° angle. Using a computer as a controller, these mirrors can be quickly swept backwards and forwards in order to vary the path length of the probe beam 14 . In this way the scanning delay line 16 assists in matching the relative path lengths of the pump and probe beams.
- the probe beam is then focussed onto receiver 26 for combining with the Terahertz beam.
- the pump beam 13 is directed onto a source 17 .
- this source 17 preferably comprises a GaAs based photoconductive switch.
- GasAs based devices use the principle of photoconductive mixing to generate their THz output.
- the THz radiation 18 emitted by the emitter 17 is directed via a hyper-hemispherical lens (not shown) towards a first parabolic mirror 19 , which is preferably an off axis parabolic (OAP) mirror, as are all the parabolic mirrors referred to herein.
- the beam is then reflected off the first parabolic mirror 19 and onto second parabolic mirror 20 , which directs the radiation onto sample 21 .
- OAP off axis parabolic
- the sample 21 may be moved relative to the beam of radiation through the focal plane of the THz beam or the beam may be moved relative to the sample or both.
- the sample may be placed on a translation stage to appropriately move the sample. This translation stage could move the sample one dimensionally along one axis of movement or through two or three axes of movement.
- the THz radiation that is reflected from sample 21 is collected by third parabolic mirror 22 and onto a fourth parabolic mirror 24 via plain mirror 23 .
- the fourth parabolic mirror 24 directs the reflected radiation 25 onto a second hyper-hemispherical lens (not shown) and onto a detector 26 , such as an electro-optic detector or a photoconductive detector.
- Photoconductive detectors comprise a detection member which may be, for example, GaAs, InGaAs, Si on Sapphire etc.
- the detection member 26 is used to detect both the amplitude and phase of the radiation emitted from the sample 20 .
- the THz radiation 25 from the sample is incident on the back surface of the detection member 26 .
- the radiation is collected by a lens (not shown), which may be hemispherical or have another shape.
- the Terahertz radiation 25 incident on the detection member 26 induces a photocurrent through the region between electrodes 26 a and 26 b , located on the opposing side of the detection member 26 , which is being illuminated by the laser radiation.
- the radiation illuminating the region between electrodes 26 a and 26 b is preferably the probe beam 14 , which carries this information.
- the current that can then be detected by the electrodes is proportional to the strength of the THz field 25 .
- the electrodes 26 a , 26 b may be of a simple diode formation embedded in a transmission line. Alternatively, they may be triangular and arranged in the shape of a bow-tie to form a so-called bow-tie antenna. They may also be interdigitated electrodes at the centre of a bow-tie or spiral antenna.
- FIG. 2 illustrates an alternative pulsed arrangement, whereby the sample is investigated using transmitted radiation rather than reflected radiation and the detector is an EOS detector rather than a photoconductive detector.
- the detector is an EOS detector rather than a photoconductive detector.
- FIG. 2 it is apparent that a variable delay is introduced to the probe beam 14 and that the pump beam 13 is projected onto the sample in the same manner as was illustrated in FIG. 1 .
- the sample in FIG. 2 is now generally perpendicular to the incident pump beam, in order to maximise radiation transmission.
- Transmitted radiation 25 is combined with the probe beam 14 .
- One particularly popular way to do this is to use electro-optic sampling (EOS).
- EOS electro-optic sampling
- the transmitted THz beam 25 and the probe beam 14 co-linearly propagate through an EOS detector 28 .
- the transmitted radiation 25 passes through the detector 28 , which modulates the probe beam 14 .
- the EOS detector 28 may be comprised of any material which possesses good non-linear characteristics, such as GaAs or Si based semiconductors and NH 4 H 2 PO 4 , ADP, KH 2 PO 4 , KH 2 ASO 4 , Quartz, AlPO 4 , ZnO, CdS, GaP, BaTiO 3 , LiTaO 3 , LiNbO 3 , Te, Se, ZnTe, ZnSe, Ba 2 NaNb 5 O 15 , AgAsS 3 , proustite, CdSe, CdGeAs 2 , AgGaSe 2 , AgSbS 3 , ZnS, organic crystals such as DAST (4-N-methylstilbazolium).
- any material which possesses good non-linear characteristics such as GaAs or Si based semiconductors and NH 4 H 2 PO 4 , ADP, KH 2 PO 4 , KH 2 ASO 4 , Quartz, AlPO 4 , ZnO, CdS
- the modulated beam is then passed into quarter wave plate 30 . This serves to circularly polarise the emitted radiation.
- the circularly polarised light is then fed through a Wollaston prism 32 , which divides the polarization of the light onto two orthogonal components. These two orthogonal components are then directed onto balanced photodiode assembly 34 .
- the balanced photodiode assembly comprises two photo diodes to respectively detect each of the orthogonal components from the Wollaston prism 32 .
- the output of the photodiodes are then linked together such that the balanced photodiode assembly 34 only outputs an electrical signal if there is a difference between the readings of the two photodiodes. This output signal corresponds to the strength of the transmitted THz beam 25 .
- THz beam 25 serves to make the radiation exiting the detector 28 slightly elliptically polarised. This change in the polarization still remains after the radiation is passed through quarter waveplate 30 . Extracting the orthogonal components of this radiation using prism 32 causes a different signal to be measured at the two photodiodes, and hence balanced photodiode assembly 34 outputs a signal corresponding to the strength of the THz field.
- EOS detection enables the phase and amplitude of the transmitted radiation to be detected. It is to be appreciated that it will be apparent to those skilled in the art that this type of analysis could be performed for any type of detector.
- FIGS. 1 and 2 it is also possible to combine the arrangements of FIGS. 1 and 2 , whereby both reflected and transmitted radiation from the sample is measured.
- THz beam with another beam of radiation which has substantially the same wavelength or which differs in frequency by at most 10 GHz.
- Such combined radiation can be detected using a bolometer, Schottky diode etc.
- FIG. 3 shows a schematic representation of an air/polythene interface 40 showing the various reflection, transmission and scattering components. It can be seen that the transmitted signal 42 is much stronger than the reflected signal 44 —96% to 4%.
- Scattering centres 46 introduce additional scattered signals 48 .
- the reflected beam is much more susceptible to these scattered pulses as it has a weaker intensity. In both transmission and reflection however scattering events can give rise to false spectral features.
- FIG. 4 a shows four time domain waveforms for samples composed of 60 mg of sucrose and 260 mg of PE powder. Each sample was measured within a powder cell.
- Curve 1 shows the recorded waveform when shows a reference trace of an empty powder cell.
- Curve 2 shows a sucrose sample with particle sizes of 53-75 ⁇ m.
- Curve 3 shows a sucrose sample with particle sizes greater than 250 ⁇ m.
- Curve 4 is also for particles greater than 250 ⁇ m but at a different location in the powder cell.
- FIG. 4 b The presence of features due to scattering can clearly be seen in FIG. 4 b (arrows 50 , 52 and 54 ) and are more common with increasing size of the sucrose particles. These scattering features would give rise to false spectral features if the data of FIG. 4 b were analysed in the frequency domain (by Fourier transforming the data).
- FIG. 5 shows a series of waveform traces measured from a mixture of fine lactose particles and large sucrose particles (800 ⁇ m size/diameter).
- Curves 2 , 3 and 4 show readings taken from three different points on the sample. It can be seen that after the main initial pulse there is a series of ripples caused by the absorption/reflection characteristics of the sample.
- Three additional features are however visible in these curves. These additional features come from scattering events and will result in false spectral features following a Fourier transform of the data.
- Curve 1 shows an averaged waveform (including the three curves 2 , 3 and 4 ) taken over an area of 5 mm 2 . It can be seen that the scattering related features are now absent from the waveform.
- FIGS. 6 a and 6 b show a comparison in the frequency domain between a spectroscopic measurement at a single point ( FIG. 6 a ) and an averaged dataset ( FIG. 6 b ).
- FIG. 6 a a Fourier transform of a waveform trace (similar to those in Curves 2 - 4 of FIG. 5 ) has been plotted. This therefore corresponds to an absorption profile.
- FIG. 6 b the waveform data has been averaged first and the Fourier transform of the averaged data has been plotted. It can been seen that, in comparison to FIG. 6 b , the trace in FIG. 6 a comprises a number of false spectral features.
- FIG. 6 b shows a rising level of absorption with increasing wavenumber. This is due to absorption by the sample and/or the scattering of polyethylene (PE). The phase difference between reference and sample measurement will also cause a similar shift of the absorption spectrum.
- the so-called “First derivative” (d(absorption)/a(Frequency) has been plotted in FIGS. 7 a and 7 b . It should be noted that FIG. 7 a corresponds to FIG. 6 a and FIGS. 7 b to 6 b.
- FIGS. 7 a and 7 b It is clear from FIGS. 7 a and 7 b that the presence of spectral features in FIG. 7 a makes spectral analysis virtually impossible.
- FIGS. 8 a and 8 b show plots of refractive index versus wavenumber.
- the refractive index is calculated from time domain and phase information.
- FIG. 8 a corresponds to FIGS. 6 a and 7 a .
- FIG. 8 b corresponds to FIGS. 6 b and 7 b . It can be seen that compared to FIG. 8 b the presence of scattering events in FIG. 8 a introduces many errors into the derivation of the refractive index.
- FIG. 9 shows various THz reflection spectra of the same sample obtained by averaging over a sample area of different size.
- the sample comprised a mixture of sucrose and lactose particles.
- a sample that is investigated using reflection spectroscopy will return a temporal waveform that comprises an initial impulse feature corresponding to the reflection from the surface of the sample and subsequent waveform data.
- This type of analysis is often referred to as the “late time response” of a sample.
- FIG. 10 shows a layered sample 62 .
- a central layer 64 of a lactose and sucrose mixture is sandwiched between two layers of polythene 66 , 68 .
- Terahertz radiation 70 is incident on one of the polythene layers and the reflected radiation 72 detected.
- the sample comprises an air/PE interface 74 and two PE/LT interfaces 76 , 78 .
- FIG. 11 shows temporal waveform traces obtained from the sample shown in FIG. 10 .
- Curve 1 is a polythene reference trace obtained by focussing the incident radiation on a separate reference sample.
- Curve 2 shows the trace obtained by focussing incident radiation on the air/PE interface.
- Curves 3 and 4 show the trace obtained by focussing the incident radiation at the first PE/LT interface (Curve 4 is obtained at a single pixel whilst curve 3 is averaged over the whole area).
- the THz signal reflected from the first PE/LT interface is much weaker than that from air/PE, ‘ ⁇ 5 ’ means that the corresponding signal is multiplied by a factor of 5.
- FIG. 12 shows a comparison of the whole response from the sample of FIG. 10 with the late time response measured from the sample following a Fast Fourier Transforms. In each case the incident radiation was focussed on the LT/PE interface.
- Curves 1 and 2 relate to the whole response.
- Curve 1 represents the averaged response from the sample and curve 2 shows a typical response from one imaging point.
- Curves 3 and 4 show the late time response from the sample. It is clear that the spectral feature of lactose at 18 cm ⁇ 1 can now be clearly resolved.
- Curve 3 represents a combination of late time response and averaging.
- FIG. 13 shows a chemical mapping image of a sucrose/lactose/polyethylene sample using the spectral features of (a) sucrose and (b) lactose.
- FIG. 14 a shows a further chemical distribution image of the surface of a sample.
- the spectral features of sucrose, lactose and polyethylene have been used to identify the location of the components on the surface of the sample.
- Sucrose is denoted by 80
- lactose by 82
- polyethylene by 84
- the area labelled 86 comprises both lactose and sucrose.
- the corresponding THz reflection spectra are shown in FIG. 14 b.
- FIG. 15 An alternative method of reducing the effects of scatter is to use a diffuse irradiating beam of radiation. Such an arrangement is shown in FIG. 15 .
- the emitter optics 88 are arranged to focus the irradiating THz radiation 90 on a single point 92 on the sample 94 .
- the focus 95 could deliberately be moved such that a “diffuse” THz irradiating beam 96 probes the sample. This will have a similar effect to the “averaging” method of the first aspect of the invention and consequently scattering effects will tend to cancel each other out across the width of the beam.
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Toxicology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
A method of investigating an object, comprising the steps of: (a) irradiating the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz; (b) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform; (c) repeating steps (a) and (b) for a plurality of points on the object and (d) combining data from step (c) to produce a time domain waveform for the object which has been averaged over the plurality of points.
Description
- The present invention relates generally to the field of apparatus and methods for imaging and/or investigating samples in the far-infrared (far-IR)/Terahertz (THz) frequency range from 100 GHz to 100 THz. Preferably the radiation utilised is in the frequency range of 500 GHz to 100 THz and more preferably from 1 THz to 100 THz and most preferably from 700 GHz to 10 THz.
- It is well established that many chemicals and pharmaceutical agents have spectral signatures in the THz region. The reason for the spectral signals in this frequency range is considered to be associated with inter-molecular or intra-molecular crystalline vibrations or collective phonon oscillations. It is known that in THz spectroscopy the spectral information of transmitted and/or scattered radiation can be used to identify materials. Most materials interact with terahertz waves to some degree, and each material has its own frequency pattern, which can be considered as a kind of ‘fingerprint’.
- Of late security concerns have increased the need for a system that is able to identify hidden explosives, such as under a person's clothing, in a person's suitcase or in a postal package.
- US Patent Application Number 2001/0033636 discloses a method of detecting explosives in luggage which uses X-rays to determine an average atomic number of a material. The average atomic number is then compared with known average atomic numbers of explosive materials to determine if the material being examined is an explosive. This technique, however is unsuitable for routine-security screening of people as its radiation is ionising.
- In a paper entitled “Spectroscopy with electronic terahertz techniques” by D. W. van der Weide et al published in Proc Spie 3828 (99), electronic pulses in the microwave range, up to 450 GHz, are used to investigate explosives. These pulses are generated electronically using non-linear transmission lines, coupled with varactor diodes that are patterned on a semiconductor. A sine wave is applied to the transmission line and the result at the other end is a rapid (1-2 ps) voltage step, which is used to produce the microwave pulse from the semiconductor device. This paper demonstrates spectral features of explosives, but only in the sub-terahertz or microwave frequency range.
- THz radiation is suitable for screening people and their possessions as it is non-ionising and can pass through clothing, paper, cardboard, wood, masonry, plastic and ceramics. Therefore, it is safer than x-ray techniques.
- It has been speculated that THz radiation could be used to obtain spectral information of explosive materials, as the organic molecular nature of explosives, and their crystalline structure are of the appropriate form. To date, however, this has not been effectively achieved.
- Many practical embodiments of THz systems utilise the transmission of THz radiation through the sample under investigation to a detector. In laboratory conditions samples can be prepared to avoid attenuation effects. However, in practice, samples are often granulated which leads to scattering of the irradiating radiation. Terahertz pulse spectroscopy is a coherent detection system which means that it has a large dynamic range. However, the effects of scattering or attenuation means that a THz signal may be too weak to detect in certain samples.
- THz systems which utilise the reflection of THz radiation from a sample will not be limited by attenuation effects in the same way that transmission-based systems are. However, reflection-based systems experience different problems which in practice has meant that effective THz reflection spectroscopy is not possible.
- In any THz system emitted Thz radiation will impinge upon the surface of the object under investigation. The surface of the object represents an interface between a medium with a first refractive index (usually air) and the object itself (which has a second refractive index which is different to the first refractive index).
- The difference in the refractive indices means that a proportion of the radiation will be transmitted into the object (for an air/polythene interface this amounts to around 96% of the incident radiation) and a proportion will be reflected (for the air/PE example this therefore equates to around 4% of the incident radiation).
- It can therefore be seen that the reflected beam will always have a weaker intensity than the transmitted beam. This means that the reflected beam is more susceptible to noise such as scattered radiation.
- A further problem arises when the sample under investigation comprises scattering centres capable of scattering any radiation that penetrates the object. For both transmission based and reflection based systems the arrival times of radiation arriving at the detector from the surface of the object and scattering centres will be different. For reflection the time of arrival will depend on the depth of the scattering centre from the surface of the object. For transmission based systems the time of arrival depends on the extra path length the scattered radiation has traveled relative to radiation that has directly passed through the object. This spread of arrival times causes artificial spectral features which can either mask actual features or create false features.
- One technique for mitigating the effects of scattering is to window the radiation detected at the detecting means. By gating the arrival time within which the system will detect photons any late arriving photons (from scattering events) will not be detected. Such a technique however reduces the efficiency of any spectroscopic analysis of the sample in question as it inevitably loses spectral information.
- It is therefore an aim of the present invention to overcome or alleviate at least some of the problems associated with both transmission based and reflection based the prior art systems.
- According to first aspect the present invention provides a method of investigating an object, comprising the steps of
-
- (a) irradiating the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
- (b) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform;
- (c) repeating steps (a) and (b) for a plurality of points on the object and
- (d) combining data from step (c) to produce a time domain waveform for the object which has been averaged over the plurality of points.
- The invention which is the subject of the present invention is suitable for use in both transmission based systems and also reflection based systems. The method comprises irradiating the object under investigation and detecting the radiation reflected/transmitted from a point on the object. The irradiation/detection steps are then repeated for a plurality of points such that an averaged time domain waveform can be calculated for the object.
- The method of the present invention acknowledges that signals due to scattering depend strongly on the geometry of the scattering centres within the object whereas directly reflected or transmitted signals do not. By taking repeated measurements the random signal elements introduced by scattering effects are removed and the true reflection/transmission signals, which have no phase dependence on the geometry of the object, remain.
- In conventional optical spectroscopy methods, the intensity of the incident (Io(ν) and the transmitted (It(ν)) radiation is measured. The optical absorption coefficient, α(ν), can then be calculated from the following equation:
-
T int(ν)≡I t(ν)/I o(ν)=e −2α(ν)d - where ν is the frequency.
- In terahertz pulsed spectroscopy (TPS) measurement however, the electric field, instead of the intensity, of the (THz) radiation is measured. The electric field of the transmitted THz radiation can be described as:
-
T(ν)≡E t(ν)/E o(ν)=e −α(ν)d e jk(ν)d -
k(ν)=n(ν)2πν/c - The fact that the electric field is measured as opposed to the intensity of the radiation, as in other spectroscopic systems, allows the detected radiation to be “averaged” as claimed in this aspect of the invention. This is because the electric field, E, can be both positive and negative and so the scattering events will tend to cancel each other out over a large enough sampling of points whereas intensities will always add.
- It is clear that both absorption coefficient (α(ν)) and refractive index (n(ν)) can be calculated from a single measurement because the electric field measurements provide both amplitude and phase information.
- At absorption peaks, the transmittance drops significantly from 100% to a few percent. Therefore the transmission measurement is very sensitive to absorption. The reflection loss at the air/sample interfaces, which is mainly determined by the refractive index difference, only causes a baseline shift in the transmission spectrum.
- In reflection measurement, things are different. The reflectance can be calculated from the measured electric field of THz radiation as:
-
R(ν)≡E r(ν)/E o(ν)=(√{square root over (ε(ν))}−1)/(√{square root over (ε(ν))}+1) -
√{square root over (ε(ν))}=n(ν)+jα(ν)c/2πν≈n(ν) - Preferably steps are taken to remove features in the detected radiation that relate to the system components. Therefore, preferably a reference measurement is taken in the absence of a sample. The signal detected in the absence of a sample can then be subtracted from the signal detected in the presence of a sample.
- Alternatively the reference measurement can be taken in the presence of a reference sample with well known THz properties. Such a reference sample could be air, water or polyethylene.
- Preferably, an absorption profile of the object can be determined by transforming the averaged time domain waveform into the frequency domain by means of a Fourier transform. The composition of the sample can be derived from this absorption profile.
- Any absorption features will produce maxima and minima in the absorption profile. Conveniently, the derivative of the absorption profile with respect to frequency can be determined.
- It is clear that the greater the area over which the sample is averaged the better the removal of the scattering features will be. However, averaging over the maximum area possible may not be necessary in order to resolve the true absorption features of a sample. Therefore, preferably, the method according to the first aspect of the present invention further comprises the step of progressively increasing the number of points that are irradiated until the scattering related effects are sufficiently reduced in order to determine the features of the object under investigation.
- When an object is irradiated with an impulse of radiation then detected radiation will comprise an initial impulse feature and subsequently received signal. The signal that follows the initial impulse function is known as the “late time response”. It has been suggested that this “late time response” provides useful information about the sample and therefore preferably the method according to the first aspect of the invention comprises the step of removing that portion of the averaged time domain waveform that corresponds to directly reflected/transmitted radiation and performing a Fourier transform on the remainder of the detected radiation in order to derive compositional information relating to the object.
- In a second aspect of the present invention the “late time response” is analysed and there is therefore provided a method of investigating an object, comprising the steps of
-
- (a) irradiating the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
- (b) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform;
- (c) separating that portion of the detected radiation that corresponds to directly reflected or transmitted radiation and applying a Fourier transform to the remaining portion of the detected radiation.
- In a third aspect of the present invention there is provided an apparatus for investigating an object comprising
-
- (a) a source of electromagnetic radiation for irradiating a point on the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
- (b) a detector for detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform;
- (c) scanning means for sequentially irradiating a plurality of points on the object and
- (d) means for combining time domain waveforms from each of the plurality of points to produce a time domain waveform for the object which has been averaged over the plurality of points.
- Compositional information relating to a sample can be determined by investigating the sample according to the method of the first aspect of the present invention. A compositional image of the sample can then be derived by imaging the sample through a spectral filter that is arranged only to let radiation relating to one (or more) of the components of the sample through to the detector. Therefore, in a fourth aspect of the present invention there is provided of imaging a sample comprising the steps of
-
- a) investigating the sample according to a method according to the first aspect of the invention;
- b) deriving a spectral waveform from the time domain waveform obtained in step (a)
- c) identifying the components of the sample from the spectral waveform in step (b)
- d) re-irradiating the sample with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz
- e) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform, the detected radiation having first passed through at least one bandpass filter corresponding to at least one component of the sample
- f) generating an image from the radiation detected in step (e)
- An alternative method of reducing the effects of scatter is to use a diffuse irradiating beam of radiation. Accordingly, in a fifth aspect of the present invention there is provided a method of investigating an object, comprising the steps of:
-
- (a) irradiating the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
- (b) detecting radiation transmitted and/or reflected from a first point on or within the object to obtain a time domain waveform;
- wherein the irradiating pulse of radiation is focussed at a second point on or within the object.
- The present invention will now be described with reference to the accompanying drawings in which:
-
FIG. 1 illustrates a schematic of a pulsed terahertz reflection investigative technique utilised in a first embodiment of the present invention; -
FIG. 2 illustrates a schematic of a pulsed terahertz transmission investigative technique utilised in a second embodiment of the present invention; -
FIG. 3 shows a schematic representation of an air/polythene interface showing the various reflection, transmission and scattering components -
FIGS. 4 a and 4 b show time domain waveforms for samples composed of sucrose and PE powder. -
FIG. 5 shows a series of waveform traces measured from a mixture of fine lactose particles and large sucrose particles (500-800 μm size/diameter). -
FIGS. 6 a and 6 b show absorption plots at a single point (FIG. 6 a) and an averaged dataset (FIG. 6 b). -
FIGS. 7 a and 7 b show the first derivative of the absorption plots ofFIGS. 6 a and 6 b. -
FIGS. 8 a and 8 b show refractive index plots of the sample imaged inFIGS. 6 a and 6 b. -
FIG. 9 shows various THz reflection spectra of the same sample for various sample areas. -
FIG. 10 shows a schematic of a layered PE/lactose/PE sample. -
FIG. 11 shows various temporal waveform traces obtained from the sample shown inFIG. 10 -
FIG. 12 shows various spectral plots obtained from the sample shown inFIG. 10 . -
FIG. 13 shows a chemical mapping image of a sucrose/lactose/polyethylene sampleFIG. 14 a shows a further chemical distribution image of the surface of a sample. -
FIG. 14 b shows temporal waveforms associated with the image shown inFIG. 14 a. -
FIG. 15 shows an imaging configuration for a THz system - Referring to
FIG. 1 , a terahertz pulsed investigating arrangement is illustrated, which comprises anultra-short pulse laser 10 which may be, for example, Ti:sapphire, Yb:Er doped fibre, Cr:LiSAF, Yb:silica, Nd:YLF, Nd:Glass, Nd:YAG, Alexandrite Yb:Phosphate Glass QX, Yb:GdCOB, Yb:YAG, Yb:KG d(WO4) or Yb:BOYS laser. Thislaser 10 emits pulses ofradiation 11, such as a collimated beam of pulses, each of which comprise a plurality of frequencies. The pulses generated by the laser preferably having a pulse duration of less than 200 fs. - The beam of generated pulses is directed into
beam splitter 12. The beam splitter splits the beam into apump beam 13, which is used to irradiate the sample, and aprobe beam 14, which is used during detection. - The
probe beam 14 is directed, viaplain mirror 15, intoscanning delay line 16.Scanning delay line 16 is a variable optical delay, which in its simplest form comprises two mirrors that serve to reflect the beam through a 180° angle. Using a computer as a controller, these mirrors can be quickly swept backwards and forwards in order to vary the path length of theprobe beam 14. In this way thescanning delay line 16 assists in matching the relative path lengths of the pump and probe beams. The probe beam is then focussed ontoreceiver 26 for combining with the Terahertz beam. - The
pump beam 13 is directed onto asource 17. For pulsed approaches thissource 17 preferably comprises a GaAs based photoconductive switch. GasAs based devices use the principle of photoconductive mixing to generate their THz output. - The
THz radiation 18 emitted by theemitter 17 is directed via a hyper-hemispherical lens (not shown) towards a firstparabolic mirror 19, which is preferably an off axis parabolic (OAP) mirror, as are all the parabolic mirrors referred to herein. The beam is then reflected off the firstparabolic mirror 19 and onto secondparabolic mirror 20, which directs the radiation ontosample 21. - To analyse a particular sample in situ, the
sample 21 may be moved relative to the beam of radiation through the focal plane of the THz beam or the beam may be moved relative to the sample or both. As shown inFIG. 1 , the sample may be placed on a translation stage to appropriately move the sample. This translation stage could move the sample one dimensionally along one axis of movement or through two or three axes of movement. - The THz radiation that is reflected from
sample 21 is collected by thirdparabolic mirror 22 and onto a fourthparabolic mirror 24 viaplain mirror 23. The fourthparabolic mirror 24 directs the reflectedradiation 25 onto a second hyper-hemispherical lens (not shown) and onto adetector 26, such as an electro-optic detector or a photoconductive detector. - Photoconductive detectors comprise a detection member which may be, for example, GaAs, InGaAs, Si on Sapphire etc. The
detection member 26 is used to detect both the amplitude and phase of the radiation emitted from thesample 20. In these detectors, theTHz radiation 25 from the sample is incident on the back surface of thedetection member 26. The radiation is collected by a lens (not shown), which may be hemispherical or have another shape. TheTerahertz radiation 25 incident on thedetection member 26 induces a photocurrent through the region betweenelectrodes detection member 26, which is being illuminated by the laser radiation. As the detector needs to know information about the phase of the radiation emitted from thegenerator 17, the radiation illuminating the region betweenelectrodes probe beam 14, which carries this information. The current that can then be detected by the electrodes is proportional to the strength of theTHz field 25. - The
electrodes -
FIG. 2 illustrates an alternative pulsed arrangement, whereby the sample is investigated using transmitted radiation rather than reflected radiation and the detector is an EOS detector rather than a photoconductive detector. To avoid unnecessary repetition with respect toFIG. 1 , like reference numerals will be used to denote like features. - In
FIG. 2 , it is apparent that a variable delay is introduced to theprobe beam 14 and that thepump beam 13 is projected onto the sample in the same manner as was illustrated inFIG. 1 . The sample inFIG. 2 , however is now generally perpendicular to the incident pump beam, in order to maximise radiation transmission. - Transmitted
radiation 25 is combined with theprobe beam 14. One particularly popular way to do this is to use electro-optic sampling (EOS). In this technique, the transmittedTHz beam 25 and theprobe beam 14 co-linearly propagate through anEOS detector 28. The transmittedradiation 25 passes through thedetector 28, which modulates theprobe beam 14. - The
EOS detector 28 may be comprised of any material which possesses good non-linear characteristics, such as GaAs or Si based semiconductors and NH4H2PO4, ADP, KH2PO4, KH2ASO4, Quartz, AlPO4, ZnO, CdS, GaP, BaTiO3, LiTaO3, LiNbO3, Te, Se, ZnTe, ZnSe, Ba2NaNb5O15, AgAsS3, proustite, CdSe, CdGeAs2, AgGaSe2, AgSbS3, ZnS, organic crystals such as DAST (4-N-methylstilbazolium). - The modulated beam is then passed into
quarter wave plate 30. This serves to circularly polarise the emitted radiation. The circularly polarised light is then fed through aWollaston prism 32, which divides the polarization of the light onto two orthogonal components. These two orthogonal components are then directed ontobalanced photodiode assembly 34. The balanced photodiode assembly comprises two photo diodes to respectively detect each of the orthogonal components from theWollaston prism 32. The output of the photodiodes are then linked together such that thebalanced photodiode assembly 34 only outputs an electrical signal if there is a difference between the readings of the two photodiodes. This output signal corresponds to the strength of the transmittedTHz beam 25. - This is because, where there is no THz beam present, there is no difference between the two photodiode signals. However, where there is a
THz beam 25, theTHz beam 25 serves to make the radiation exiting thedetector 28 slightly elliptically polarised. This change in the polarization still remains after the radiation is passed throughquarter waveplate 30. Extracting the orthogonal components of thisradiation using prism 32 causes a different signal to be measured at the two photodiodes, and hencebalanced photodiode assembly 34 outputs a signal corresponding to the strength of the THz field. - Therefore EOS detection enables the phase and amplitude of the transmitted radiation to be detected. It is to be appreciated that it will be apparent to those skilled in the art that this type of analysis could be performed for any type of detector.
- It is also to be appreciated that it is also possible to combine the arrangements of
FIGS. 1 and 2 , whereby both reflected and transmitted radiation from the sample is measured. - Further, rather than combining the beam which has been reflected from or transmitted by the sample with the
probe beam 14, it is also possible to combine the THz beam with another beam of radiation which has substantially the same wavelength or which differs in frequency by at most 10 GHz. Such combined radiation can be detected using a bolometer, Schottky diode etc. -
FIG. 3 shows a schematic representation of an air/polythene interface 40 showing the various reflection, transmission and scattering components. It can be seen that the transmittedsignal 42 is much stronger than the reflectedsignal 44—96% to 4%. - Scattering centres 46 introduce additional scattered signals 48. The reflected beam is much more susceptible to these scattered pulses as it has a weaker intensity. In both transmission and reflection however scattering events can give rise to false spectral features.
-
FIG. 4 a shows four time domain waveforms for samples composed of 60 mg of sucrose and 260 mg of PE powder. Each sample was measured within a powder cell.Curve 1 shows the recorded waveform when shows a reference trace of an empty powder cell.Curve 2 shows a sucrose sample with particle sizes of 53-75 μm.Curve 3 shows a sucrose sample with particle sizes greater than 250 μm.Curve 4 is also for particles greater than 250 μm but at a different location in the powder cell. - In order to show the waveforms due to the sample only, the recorded waveforms from
FIG. 4 a were then de-convolved with the reference trace (curve 1). The de-convolved THz waveforms are shown inFIG. 4 b. It should be noted that in bothFIGS. 4 a and 4 b the curves have been offset for clarity. Also, the signal intensity ofcurve 1 has been reduced by a factor of 4. - The presence of features due to scattering can clearly be seen in
FIG. 4 b (arrows FIG. 4 b were analysed in the frequency domain (by Fourier transforming the data). -
FIG. 5 shows a series of waveform traces measured from a mixture of fine lactose particles and large sucrose particles (800 μm size/diameter).Curves arrow 56 incurve 2,arrow 58 incurve 3 andarrow 60 in curve 4) are however visible in these curves. These additional features come from scattering events and will result in false spectral features following a Fourier transform of the data. -
Curve 1 shows an averaged waveform (including the threecurves -
FIGS. 6 a and 6 b show a comparison in the frequency domain between a spectroscopic measurement at a single point (FIG. 6 a) and an averaged dataset (FIG. 6 b). InFIG. 6 a a Fourier transform of a waveform trace (similar to those in Curves 2-4 ofFIG. 5 ) has been plotted. This therefore corresponds to an absorption profile. - In
FIG. 6 b the waveform data has been averaged first and the Fourier transform of the averaged data has been plotted. It can been seen that, in comparison toFIG. 6 b, the trace inFIG. 6 a comprises a number of false spectral features. - The trace in
FIG. 6 b shows a rising level of absorption with increasing wavenumber. This is due to absorption by the sample and/or the scattering of polyethylene (PE). The phase difference between reference and sample measurement will also cause a similar shift of the absorption spectrum. In order to remove this aspect and aid interpretation of the spectral features the so-called “First derivative” (d(absorption)/a(Frequency) has been plotted inFIGS. 7 a and 7 b. It should be noted thatFIG. 7 a corresponds toFIG. 6 a andFIGS. 7 b to 6 b. - It is clear from
FIGS. 7 a and 7 b that the presence of spectral features inFIG. 7 a makes spectral analysis virtually impossible. -
FIGS. 8 a and 8 b show plots of refractive index versus wavenumber. The refractive index is calculated from time domain and phase information. Again,FIG. 8 a corresponds toFIGS. 6 a and 7 a.FIG. 8 b corresponds toFIGS. 6 b and 7 b. It can be seen that compared toFIG. 8 b the presence of scattering events inFIG. 8 a introduces many errors into the derivation of the refractive index. -
FIG. 9 shows various THz reflection spectra of the same sample obtained by averaging over a sample area of different size. The sample comprised a mixture of sucrose and lactose particles. - The number of points (pixels) used in averaging is shown next to each spectrum and the images have been off-set for clarity. It is clear that averaging over a larger sample area improves the results. However, in this example it can be seen that a sample area of 13×13 pixels is sufficient to resolve the absorption features of lactose and sucrose.
- A sample that is investigated using reflection spectroscopy will return a temporal waveform that comprises an initial impulse feature corresponding to the reflection from the surface of the sample and subsequent waveform data. By analysing the waveform trace following the initial impulse feature further information can be derived about a sample. This type of analysis is often referred to as the “late time response” of a sample.
-
FIG. 10 shows a layeredsample 62. Acentral layer 64 of a lactose and sucrose mixture is sandwiched between two layers ofpolythene Terahertz radiation 70 is incident on one of the polythene layers and the reflectedradiation 72 detected. The sample comprises an air/PE interface 74 and two PE/LT interfaces 76, 78. -
FIG. 11 shows temporal waveform traces obtained from the sample shown inFIG. 10 .Curve 1 is a polythene reference trace obtained by focussing the incident radiation on a separate reference sample.Curve 2 shows the trace obtained by focussing incident radiation on the air/PE interface.Curves Curve 4 is obtained at a single pixel whilstcurve 3 is averaged over the whole area). The THz signal reflected from the first PE/LT interface is much weaker than that from air/PE, ‘×5 ’ means that the corresponding signal is multiplied by a factor of 5. -
FIG. 12 shows a comparison of the whole response from the sample ofFIG. 10 with the late time response measured from the sample following a Fast Fourier Transforms. In each case the incident radiation was focussed on the LT/PE interface. -
Curves Curve 1 represents the averaged response from the sample andcurve 2 shows a typical response from one imaging point.Curves Curve 3 represents a combination of late time response and averaging. - By using a method relating to the present invention it is possible to identify the true spectral components of a sample. Once these features have been identified it is possible to image the sample at the wavelengths relating to the components of the sample. In this way a spectral image can be compiled of the sample.
FIG. 13 shows a chemical mapping image of a sucrose/lactose/polyethylene sample using the spectral features of (a) sucrose and (b) lactose. -
FIG. 14 a shows a further chemical distribution image of the surface of a sample. The spectral features of sucrose, lactose and polyethylene have been used to identify the location of the components on the surface of the sample. Sucrose is denoted by 80, lactose by 82, and polyethylene by 84. The area labelled 86 comprises both lactose and sucrose. The corresponding THz reflection spectra are shown inFIG. 14 b. - An alternative method of reducing the effects of scatter is to use a diffuse irradiating beam of radiation. Such an arrangement is shown in
FIG. 15 . In a typical system configuration theemitter optics 88 are arranged to focus the irradiatingTHz radiation 90 on asingle point 92 on thesample 94. However, as an alternative thefocus 95 could deliberately be moved such that a “diffuse”THz irradiating beam 96 probes the sample. This will have a similar effect to the “averaging” method of the first aspect of the invention and consequently scattering effects will tend to cancel each other out across the width of the beam.
Claims (10)
1. A method of investigating an object, comprising the steps of:
(a) irradiating the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
(b) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform;
(c) repeating steps (a) and (b) for a plurality of points on the object and
(d) sequentially combining data from step (c) to produce a time domain waveform for the object which has been averaged over the plurality of points, said sequential combining and averaging continuing until the effect of scattering in the produced time domain waveform has been reduced.
2. The method of claim 1 further comprising the steps of:
(b)(i) obtaining a reference measurement taken in the absence of an object;
(b)(ii) detecting radiation in the absence of an object; and
(b)(iii) subtracting the signal measured by the detector during step (b)(ii) from the signal measured by the detector during step (b).
3. The method of claim 1 wherein the Fourier transform of the averaged time domain waveform from step (d) is calculated.
4. The method of claim 3 further comprising obtaining a first derivative of the obtained spectra.
5. The method of claim 1 wherein the plurality of points in step (c) is iteratively increased until the waveform from step (d) is substantially unchanging.
6. (canceled)
7. (canceled)
8. An apparatus for investigating an object comprising:
a source of electromagnetic radiation for irradiating a point on the object with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
a detector for detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform;
scanning means for sequentially irradiating a plurality of points on the object and
means for sequentially combining time domain waveforms from each of the plurality of points to produce a time domain waveform for the object which has been averaged over the plurality of points, said sequential combining and averaging continuing until the effect of scattering in the produced time domain waveform has been reduced.
9. A method of imaging a sample comprising the steps of:
(a) investigating the sample according to a method as claimed in claim 1 ;
(b) deriving a spectral waveform from the time domain waveform obtained in step (a);
(c) identifying the components of the sample from the spectral waveform in step (b);
(d) re-irradiating the sample with an optically-generated pulse of electromagnetic radiation, said pulse having a plurality of frequencies in the range from 100 GHz to 100 THz;
(e) detecting radiation transmitted and/or reflected from the object to obtain a time domain waveform, the detected radiation having first passed through at least one bandpass filter corresponding to at least one component of the sample; and
(f) generating an image from the radiation detected in step (e).
10-12. (canceled)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/613,355 US20100108888A1 (en) | 2004-08-26 | 2009-11-05 | Reducing scattering related features in terhertz time domain spectroscopy bz averaging the impulse waveform over a plurality of sample points |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB0419098.9 | 2004-08-26 | ||
GB0419098A GB2417554B (en) | 2004-08-26 | 2004-08-26 | Spectroscopy apparatus and associated technique |
US66109007A | 2007-03-27 | 2007-03-27 | |
US12/613,355 US20100108888A1 (en) | 2004-08-26 | 2009-11-05 | Reducing scattering related features in terhertz time domain spectroscopy bz averaging the impulse waveform over a plurality of sample points |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US66109007A Continuation | 2004-08-26 | 2007-03-27 |
Publications (1)
Publication Number | Publication Date |
---|---|
US20100108888A1 true US20100108888A1 (en) | 2010-05-06 |
Family
ID=34307605
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/661,090 Active US7675036B2 (en) | 2004-08-26 | 2005-08-26 | Reducing scattering related features in terhertz time domain spectroscopy by averaging the impulse waveform over a plurality of sample points |
US12/613,355 Abandoned US20100108888A1 (en) | 2004-08-26 | 2009-11-05 | Reducing scattering related features in terhertz time domain spectroscopy bz averaging the impulse waveform over a plurality of sample points |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/661,090 Active US7675036B2 (en) | 2004-08-26 | 2005-08-26 | Reducing scattering related features in terhertz time domain spectroscopy by averaging the impulse waveform over a plurality of sample points |
Country Status (5)
Country | Link |
---|---|
US (2) | US7675036B2 (en) |
EP (2) | EP1792165A1 (en) |
JP (1) | JP2008510980A (en) |
GB (1) | GB2417554B (en) |
WO (1) | WO2006021799A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020104011A1 (en) * | 2018-11-19 | 2020-05-28 | Tera Group Ag | Method and apparatus for classifying a sample by means of terahertz spectroscopy |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9909986B2 (en) | 2003-10-15 | 2018-03-06 | Applied Research And Photonics, Inc. | Thickness determination and layer characterization using terahertz scanning reflectometry |
US8620132B2 (en) | 2006-09-27 | 2013-12-31 | Anis Rahman | Terahertz scanning reflectometer |
US7368292B2 (en) | 2006-03-10 | 2008-05-06 | University Of Florida Research Foundation, Inc. | Differential reflection spectroscopy system and method for detecting explosives and other target materials |
WO2007121598A1 (en) * | 2006-04-21 | 2007-11-01 | Eth Zurich | Broadband terahertz radiation generation and detection system and method |
DE102006042642B4 (en) * | 2006-09-12 | 2010-06-24 | Batop Gmbh | Terahertz time-domain spectrometer |
CN101210874B (en) * | 2006-12-31 | 2010-05-19 | 清华大学 | Method and apparatus for measuring terahertz time-domain spectroscopy |
GB2445957B (en) * | 2007-01-29 | 2009-08-19 | Tera View Ltd | A method and apparatus for analysing an LCD |
JP5147454B2 (en) * | 2007-03-07 | 2013-02-20 | キヤノン株式会社 | Image forming apparatus and image forming method |
US7767968B2 (en) * | 2007-09-18 | 2010-08-03 | Honeywell International Inc. | Correlated ghost imager |
EP2263101A2 (en) | 2008-03-18 | 2010-12-22 | Manchester Metropolitan University | Remote detection and measurement of objects |
JP2009300108A (en) | 2008-06-10 | 2009-12-24 | Sony Corp | Terahertz spectroscopic device |
JP2009300109A (en) * | 2008-06-10 | 2009-12-24 | Sony Corp | Terahertz wave measuring method and terahertz spectroscopic device |
US8113427B2 (en) * | 2008-12-18 | 2012-02-14 | Ncr Corporation | Methods and apparatus for automated product identification in point of sale applications |
GB0916300D0 (en) | 2009-09-17 | 2009-10-28 | Univ Manchester Metropolitan | Remote detection of bladed objects |
DE102010032382A1 (en) | 2010-07-27 | 2012-02-02 | Batop Gmbh | Fiber-coupled terahertz time domain spectrometer has pulse laser that uses transmitting antenna and receiving antenna, where both antennas are assigned with collimating terahertz optics in each case |
US20120273681A1 (en) * | 2011-04-26 | 2012-11-01 | Zomega Terahertz Corporation | Terahertz spectrometer |
US10269676B2 (en) * | 2012-10-04 | 2019-04-23 | Taiwan Semiconductor Manufacturing Company, Ltd. | Thermally enhanced package-on-package (PoP) |
US9513213B2 (en) * | 2013-08-08 | 2016-12-06 | Schlumberger Technology Corporation | System and method of determining rock properties using terahertz-band dielectric measurements |
KR101510228B1 (en) | 2014-08-12 | 2015-04-10 | 연세대학교 산학협력단 | Gas detection method and a gas sensor usinng terahertz, method for manufacturing the gas sensor and gas detection system |
JP6266719B1 (en) * | 2016-08-31 | 2018-01-24 | フェムトディプロイメンツ株式会社 | Terahertz time domain spectrometer |
CN107167771B (en) * | 2017-04-28 | 2018-10-26 | 深圳市无牙太赫兹科技有限公司 | A kind of the direct wave suppressing method and system of microwave imaging system |
DE102020113306A1 (en) | 2020-05-15 | 2021-11-18 | Helmut Fischer GmbH Institut für Elektronik und Messtechnik | Device for transmitting and / or receiving terahertz radiation and control device therefor |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4101773A (en) * | 1975-11-04 | 1978-07-18 | Emi Limited | Radiography |
US4517835A (en) * | 1983-09-19 | 1985-05-21 | Dresser Industries, Inc. | Method for determining formation dip |
US5473412A (en) * | 1993-03-08 | 1995-12-05 | Nikon Corporation | Energy amount controlling method |
US5710430A (en) * | 1995-02-15 | 1998-01-20 | Lucent Technologies Inc. | Method and apparatus for terahertz imaging |
US6601459B1 (en) * | 1999-10-29 | 2003-08-05 | Universitat Zurich | Method of volumetric blood flow measurement |
US6777684B1 (en) * | 1999-08-23 | 2004-08-17 | Rose Research L.L.C. | Systems and methods for millimeter and sub-millimeter wave imaging |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6078047A (en) * | 1997-03-14 | 2000-06-20 | Lucent Technologies Inc. | Method and apparatus for terahertz tomographic imaging |
DE19954663B4 (en) | 1999-11-13 | 2006-06-08 | Smiths Heimann Gmbh | Method and device for determining a material of a detected object |
WO2002017231A2 (en) * | 2000-08-23 | 2002-02-28 | Rose Research Llc | Systems and methods for millimeter and sub-millimeter wave imaging |
US6723991B1 (en) * | 2000-10-20 | 2004-04-20 | Imra America, Inc. | Single-shot differential spectroscopy and spectral-imaging at submillimeter wavelengths |
GB0122052D0 (en) * | 2001-09-12 | 2001-10-31 | Teraview Ltd | Imaging apparatus and method |
US7335883B2 (en) * | 2001-09-12 | 2008-02-26 | Teraview Limited | Imaging apparatus and method |
US20050087690A1 (en) * | 2001-12-28 | 2005-04-28 | Mamoru Usami | Spectral measurnig device |
US7368280B2 (en) * | 2002-05-23 | 2008-05-06 | Rensselaer Polytechnic Institute | Detection of biospecific interactions using amplified differential time domain spectroscopy signal |
GB2397207B (en) * | 2003-01-10 | 2005-04-13 | Teraview Ltd | Imaging techniques and associated apparatus |
GB2399626B (en) * | 2003-03-21 | 2006-04-05 | Teraview Ltd | Spectroscopy apparatus and associated technique |
GB2405200B (en) * | 2003-08-22 | 2005-09-07 | Teraview Ltd | A method and apparatus for investigating a sample |
CN101203742B (en) * | 2004-05-26 | 2011-10-19 | 派克米瑞斯有限责任公司 | Terahertz imaging in reflection and transmission mode for luggage and personnel inspection |
GB2416204B (en) * | 2004-07-16 | 2007-03-21 | Teraview Ltd | Apparatus and method for investigating a sample |
-
2004
- 2004-08-26 GB GB0419098A patent/GB2417554B/en not_active Expired - Fee Related
-
2005
- 2005-08-26 EP EP05775318A patent/EP1792165A1/en not_active Withdrawn
- 2005-08-26 US US11/661,090 patent/US7675036B2/en active Active
- 2005-08-26 WO PCT/GB2005/003333 patent/WO2006021799A1/en active Application Filing
- 2005-08-26 JP JP2007528982A patent/JP2008510980A/en active Pending
- 2005-08-26 EP EP08014345A patent/EP2023124A2/en not_active Withdrawn
-
2009
- 2009-11-05 US US12/613,355 patent/US20100108888A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4101773A (en) * | 1975-11-04 | 1978-07-18 | Emi Limited | Radiography |
US4517835A (en) * | 1983-09-19 | 1985-05-21 | Dresser Industries, Inc. | Method for determining formation dip |
US5473412A (en) * | 1993-03-08 | 1995-12-05 | Nikon Corporation | Energy amount controlling method |
US5710430A (en) * | 1995-02-15 | 1998-01-20 | Lucent Technologies Inc. | Method and apparatus for terahertz imaging |
US6777684B1 (en) * | 1999-08-23 | 2004-08-17 | Rose Research L.L.C. | Systems and methods for millimeter and sub-millimeter wave imaging |
US6601459B1 (en) * | 1999-10-29 | 2003-08-05 | Universitat Zurich | Method of volumetric blood flow measurement |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020104011A1 (en) * | 2018-11-19 | 2020-05-28 | Tera Group Ag | Method and apparatus for classifying a sample by means of terahertz spectroscopy |
Also Published As
Publication number | Publication date |
---|---|
WO2006021799A1 (en) | 2006-03-02 |
EP2023124A2 (en) | 2009-02-11 |
GB2417554B (en) | 2007-08-08 |
JP2008510980A (en) | 2008-04-10 |
US7675036B2 (en) | 2010-03-09 |
GB0419098D0 (en) | 2005-03-09 |
EP1792165A1 (en) | 2007-06-06 |
GB2417554A (en) | 2006-03-01 |
US20080006767A1 (en) | 2008-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7675036B2 (en) | Reducing scattering related features in terhertz time domain spectroscopy by averaging the impulse waveform over a plurality of sample points | |
US7728296B2 (en) | Spectroscopy apparatus and associated technique | |
US6865014B2 (en) | Apparatus and method for investigating a sample | |
JP5166024B2 (en) | Terahertz imaging in reflection and transmission modes for inspection of luggage and people | |
US20050082479A1 (en) | Imaging apparatus and method | |
US20090066948A1 (en) | Compact Terahertz Spectrometer Using Optical Beam Recycling and Heterodyne Detection | |
US6794670B1 (en) | Method and apparatus for spectrometric analysis of turbid, pharmaceutical samples | |
US9404854B2 (en) | Second and third order simultaneously non-linear optical processes and measurements for surface analysis | |
US20030155512A1 (en) | Apparatus and method for investigating a sample | |
Baker et al. | People screening using terahertz technology | |
GB2435509A (en) | Terahertz pulsed spectroscopy method | |
JP2003279412A (en) | Spectrometry device using white electromagnetic wave emitted from photonic crystal member generated by single pulse photoexcitation as light source | |
Lo et al. | Pulsed terahertz bi-directional reflection distribution function (BRDF) measurements of materials and obscurants | |
US9557267B2 (en) | Terahertz imaging via simultaneous surface and sub-surface evaluation via non-linear optical response | |
Jeunesse et al. | “On-the-Fly” Monitoring With a Single-Shot Terahertz Time-Domain Spectrometer | |
GB2350673A (en) | Three dimensional imaging using terahertz or Far IR radiation | |
Palka | Detection of covered materials in the TDS-THz setup | |
Kong et al. | Terahertz time-domain spectroscopy for explosive trace detection | |
Beigang | Industrial applications of terahertz technology: from layer thickness measurements to mail screening | |
Bernacki et al. | Differential spectroscopic imaging of particulate explosives residue | |
Kawase et al. | Terahertz rays to detect drugs of abuse | |
Kawase et al. | Terahertz spectral imaging for drug inspection | |
Zhou et al. | Impacts of terahertz scattering on the reflection spectrum for explosive detection | |
JP2010019647A (en) | Electromagnetic wave measuring instrument and method | |
Panzner et al. | THz-Time-Domain Spectroscopy-A new tool for the analysis of artwork |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |