1

Numerical Removal of Water-Vapor Effects from
THz-TDS Measurements

arXiv:0710.3621v1 [cs.CE] 19 Oct 2007

Withawat Withayachumnankul, Bernd M. Fischer, Samuel P. Mickan, Member, IEEE,
and Derek Abbott, Fellow, IEEE

Abstract— One source of disturbance in a pulsed T-ray signal
is attributed to ambient water vapor. Water molecules in the
gas phase selectively absorb T-rays at discrete frequencies corresponding to their molecular rotational transitions. This results
in prominent resonances spread over the T-ray spectrum, and in
the time domain the T-ray signal is observed as fluctuations after
the main pulse. These effects are generally undesired, since they
may mask critical spectroscopic data. So, ambient water vapor
is commonly removed from the T-ray path by using a closed
chamber during the measurement. Yet, in some applications
a closed chamber is not applicable. This situation, therefore,
motivates the need for another method to reduce these unwanted
artifacts. This paper presents a study on a computational means
to address the problem. Initially, a complex frequency response
of water vapor is modeled from a spectroscopic catalog. Using a
deconvolution technique, together with fine tuning of the strength
of each resonance, parts of the water-vapor response are removed
from a measured T-ray signal, with minimal signal distortion.
Index Terms— T-rays, Terahertz, THz-TDS, water vapor, rotational transitions, removal

I. I NTRODUCTION
Terahertz time-domain spectroscopy (THz-TDS) has received much attention from researchers due to its outstanding
performance [1], [2]. With the assistance of ultrafast femtosecond laser technology, a T-ray system generates and detects
a short pulse, and thus the corresponding frequency range
spanning from a few hundred gigahertz to a few terahertz or
more. This frequency range has a relatively high black-body
radiation background, prohibiting a high-SNR detection with
a conventional detector. However, this problem is suppressed
by adopting a coherent detection scheme, lifting the SNR to
as high as 60 dB [3].
In the T-ray frequency region, many polar gases of general
interest possess unique rotational transition energies, which
give rise to spectral resonances [4]. Because of these unique
fingerprints, THz-TDS proves useful for gas classification and
recognition. However, this property, on the other hand, affects
the THz-TDS capability in an open air setting, in which water
vapor is ubiquitous.
Manuscript received August 9, 2007; revised ..., 2007. This program was
supported by The Australian Research Council (ARC), the Sir Ross and Sir
Keith Smith Fund, the Defence Science and Technology Organisation (DSTO),
and NHEW P/L.
The authors are with the Centre for Biomedical Engineering (CBME) and
the School of Electrical & Electronic Engineering, The University of Adelaide,
SA 5005, Australia.
S. P. Mickan is also with Davies Collison Cave, 1 Nicholson Street,
Melbourne, VIC 3000, Australia.
Corresponding author: W. W. (withawat@eleceng.adelaide.edu.au).

Water vapor, the third most abundant gas in nature [5],
is known to have numerous rotational resonances in the Tray band [6]. THz-TDS of a sample, in open air, therefore
unavoidably results in a combination of the sample’s spectral
features and water vapor resonances in the frequency domain.
In the time domain, this results in field fluctuations after the
main T-ray pulse. Mostly, these effects are undesirable, since
they can obscure spectroscopic results of interest.
The water-vapor effects are removable during the measurement by purging an enclosed T-ray path with dry air
or such a non-polar gas as nitrogen, which does not have
transition energy levels in the T-ray regime [7]—alternatively
a vacuum is sometimes used. In some applications, it is not
always possible to enclose the entire T-ray beam path, for
example, in applications where stand-off detection is required
[8]. The effects can be partially canceled out by calculation,
if two measurements, one for the sample and the other for
the reference, are available. Again, not all situations can
provide both measurements. Performing narrowband T-ray
sensing in an atmospheric transmission window could avoid
the effects [9], at the expense of a full spectral fingerprint.
As a completely different approach, a numerical method to
alleviate the water vapor effects on the data would therefore
be beneficial.
Knowing the resonances’ characteristics in the T-ray frequency range allows a numerical estimation of the complex
response of water vapor. Theoretically, the estimated complex
response can be deconvolved directly from a received Tray signal. But pragmatically, fitting of a modeled complex
response to a measured response is complicated by many
factors, e.g. the limited dynamic range [10], limited frequency
resolution [11], and measurement uncertainties [12]. Besides,
the exact model requires precise knowledge of geometric and
atmospheric conditions during the measurement. Moreover, if
direct deconvolution is carried out, there is no measure to
validate the result.
Fine-tuning the strengths of complex resonances based on a
brute-force search is introduced in this work. Each resonance
is tuned in magnitude within a predefined range and then
deconvolved from a measured signal. A tuning criterion is
met when the ratio of the fluctuation energy to the main pulse
energy of an adjusted signal is minimized. Repeatedly tuning
the strength line-by-line ultimately results in the mitigation of
water-vapor-induced effects. Constrained by the condition of
minimum fluctuation energy, the algorithm always provides
the optimum results.
The proposed algorithm based on signal processing has

2

amplitude (a.u.)

nitrogen atmosphere
water-vapor atmosphere

0

5

10

103

15
time (ps)

20

In addition to the fluctuation and the resonance effects, the
analysis shows that the T-ray energy loss by water absorption
calculated between 0.0 and 4.0 THz is as high as 10% for the
spectrum in Fig. 1. The loss is likely due to the non-directional
energy re-emission of water molecules. In the time domain, the
ratio of the main pulse energy to the tail (fluctuation) energy
is calculated for both of the time-domain signals (for more
details on the calculation, see Section IV-C). The energy ratio
for the nitrogen measurement is 429.98, whereas the ratio for
the water vapor measurement is 18.82. This measure will be
used later on to construct a criterion for fluctuation removal.

25

nitrogen atmosphere
water-vapor atmosphere

2

10

101
noise floor

1
0

1

2

3
4
frequency (THz)

5

6

Fig. 1. Effects of water vapor on a T-ray pulse and its spectrum. (above) Tray signals recorded in nitrogen and water-vapor atmospheres. The sampling
interval is 0.067 ps, and the total time duration is 136 ps. (below) Their
corresponding spectra with the spectral resolution of 7.32 GHz. For this
measurement the noise floor marks the cutoff frequency at approximately
4.0 THz.

merits in that (i) its generality allows application to T-ray
signals with reduced effect on samples’ temporal and spectral
features, and (ii) the exact measurement conditions, including
humidity, temperature, pressure, and propagation length, are
not crucial in order to remove the water-vapor response.
This paper is organized as follows: Section II elaborates
the effects of water vapor on T-ray signals and spectra. Section III, the model characteristics of each complex resonance is
determined and validated with a measurement. In Section IV,
the algorithm for removal of water-vapor effects is introduced.
This algorithm is verified by experiments in Section V, followed by discussions and a conclusion in Section VI.
II. E FFECTS OF WATER VAPOR ON PULSED T- RAY SIGNAL
Water vapor contains many modes and the details of its
interaction with an electromagnetic wave depends on the frequency. Water vapor exhibits pure rotational modes of energy
transitions spanning the millimeter wave to mid-infrared or
approximately 3 GHz to 19.5 THz [13], [14]. In the infrared
region, both pure rotational line transitions and rotationvibration bands are noticeable [15]. These transitions cause
absorption and re-emission of wave energy in a number of
narrow frequency bands, unique to water molecules.
Fig. 1 shows T-ray signals and spectra measured by THzTDS under nitrogen and water-vapor atmospheres at room
temperature and pressure. In the time domain, the signal
recorded in the presence of water vapor undergoes long
fluctuations after the main pulse. This is due to the energy
re-emission by rotational transitions of water molecules. In
the frequency domain, the water vapor causes several sharp
resonances at discrete frequencies, as a result of quantised
rotational transition energies. In terms of Fourier theory,
the pair of time-domain fluctuations and frequency-domain
resonances is based on the principle that a sharp feature in
one domain is related to a broad feature in the other domain.

III. M ODEL OF WATER VAPOR ABSORPTION AND
DISPERSION

It is required that the water-vapor response be modeled
appropriately before a further step in water-vapor removal is
taken. Modeling of absorption and dispersion for a rotational
resonance involves selecting a suitable lineshape and calculating the line position, line strength, and linewidth. Fortunately,
the line position and line strength are, only to a small extent,
dependent on atmospheric conditions, i.e. the line position is
constant at low pressures [16] and the peak absorption is not
a function of the pressure [16], [17]. In fact, the line position
and line strength have been well parameterized in many
publications. Thus these two parameters are readily available,
leaving only the linewidth and lineshape to be determined. In
the following subsections these parameters are discussed in
detail relative to water molecule.
A. Line strength and line position
Water vapor rotational resonances in the T-ray regime have
been extensively measured and analyzed via either conventional FTIR spectroscopy [18], [19], THz-TDS [6], [20], or
other techniques [21]–[23]—the extensive study is probably
due to the abundance of water vapor, which involves many
physical processes in nature. Subsequently, the computeraccessible spectroscopic parameters of water vapor, including
the line position, line strengths and linewidth, are available
from many research groups. Here in this work the catalog
published by the JPL group [24] is adopted. As this catalog has
been regularly updated since before 1985 [25], it is unlikely
that the high-intensity lines of water vapor are absent from the
list.
According to the JPL catalog, the line strength is reported in
terms of the integrated line intensity at 300 K, Ia (300 K), and
is scalable with the temperature, T . Despite the temperature
dependency, the line strength is independent of pressure [16].
Also included in the catalog is the line position. The shifting
of the position is an order of magnitude lower than the line
broadening [26], and thus the position is presumably pressureand temperature-dependent [16].
B. Linewidth
A rotational absorption resonance is possibly broadened
by a number of factors, for instance, self or foreign-gas
collisions, molecule-wall collisions, the Doppler effect, and

WITHAYACHUMNANKUL et al.: NUMERICAL REMOVAL OF WATER-VAPOR EFFECTS FROM THZ-TDS MEASUREMENTS

where ∆ω0 is the FWHM of the line at pressure p0 and temperature T0 , ∆ω is the FWHM at pressure p and temperature
T , and m is the temperature index, varying between 0.5 and
1.0. An index of m = 0.68 can be used in the calculation [30].
Determined from the experiment, the average FWHM,
∆ω0 /2π, of water vapor resonance at ambient temperature
and standard pressure is 6 GHz. This value is in agreement
with an estimate of 10 MHz per Torr, or 7.6 GHz at 1 atm [27],
[31] within the limit of measurement uncertainty.
C. Lineshape
The rotational absorption lineshape can be modeled either
by Lorentz [32], Vleck-Weisskopf [33], or Gross [13], [34]
profiles. The choice depends on the frequency range of interest
and the nature of line broadening. Notwithstanding, no significant difference is observable among these profiles at or nearby
a T-ray resonance, with the linewidth nearly three orders of
magnitude lower than the resonance frequency. In this paper,
a Lorentzian profile is adopted in modeling the water vapor
response in the T-ray frequency range.
A Lorentz absorption and dispersion profiles are given
by [13], respectively,


1
c∆ωa
, (2a)
κa (ω) =
2π
(ωa − ω)2 + ∆ωa2


c
ωa − ω
na (ω) − 1 =
,
(2b)
2π (ωa − ω)2 + ∆ωa2
where ωa denotes the ath transition frequency, and ∆ωa
denotes the HWHM of the profile at that frequency. The
absorption and dispersion profiles are linked together via
Kramers-Krönig relation.
D. Ensemble of rotational transition resonances
This subsection gives a complete model, in terms of optical
constants, as a result of several rotational transition resonances
in the frequency range of interest. The model is a function of
the lineshape, linewidth, line strength, and line position, which
are discussed in the previous subsections. Later on, the model
is compared to the measurement to verify the accuracy.
The model extinction coefficient and index of refraction are
summation of an offset and a set of Lorentzian profiles, or,
respectively,
X
κ(ω) = κ(0) +
ma κa (ω) ,
(3a)
n(ω) − 1

= n(0) +

a
X
a

ma [na (ω) − 1] ,

(3b)

0.3
absorption coefficient, α (cm-1)

natural lifetime broadening [16], [27]. Among these factors, at
a standard pressure the collisional broadening is predominant
in determining the spectral linewidth [17], [28].
More specifically, the collision broadened linewidth is dependent on pressure [27] and temperature [15], [20], [26]. This
dependency yields the Benedict and Kaplan relationship [13],
[29],
   m
T0
p
,
(1)
∆ω = ∆ω0
p0
T

0.2

3

measured absorption
modeled absorption
maximum measurable absorption

0.1
0
0
0.3

0.5

1

1.5

2

2.5

3
frequency (THz)

3.5

4

0.2
0.1
0
2

Fig. 2. Model and measurement of H2 O absorption. The lower graph is a
continuation of the upper graph, as indicated by the arrowheads. The measured
and modeled curves match well in the range of 0.0 to 1.6 THz, but match
poorly in the higher frequency range, where the absorption strengths become
higher than the dynamic range of the system. The maximum measurable
absorption coefficient of the system, represented by the dashed line, is
calculated using the method of Jepsen and Fischer [10].

where the line strength ma is proportional to the integrated
line intensity, Ia , given in the JPL catalog, divided by the
line position, ωa , or, ma ∝ Ia /ωa . The offset is partially a
result of electronic and molecular vibrational resonances at
high frequencies [13], and also a contribution from the tails
of other rotational resonances unaccounted for by the model.
Typically, at the low frequency range, from DC to a few
hundred gigahertz, THz-TDS cannot produce sufficient energy
to overcome the noise. As a result, the resolved optical
constants in this frequency range are unreliable. A phase
extrapolation technique is usually introduced to correct the
unwrapping process, forcing the phase to start at zero [35].
Thus, it is not necessary in the model to consider the index
offset, n(0), which is derived from the phase.
The absorption of water vapor, shown in Fig. 2, is determined from the measured data (Fig. 1) and the Lorentzian
model in the frequency range between 0.0 and 4.0 THz. It
can be clearly seen that below 1.6 THz the model closely
resembles the measurement. This match is possible since the
resonances in this low frequency range have their strengths
beneath the maximum absorption coefficient measurable by
the system. However, as the frequency goes beyond 1.6 THz,
the model cannot track the measurement due to the dynamic
range of the system. In this situation, the system can no longer
measure the absorption coefficient correctly, and the clipped
absorption peaks are obvious, in particular, in the range from
2.0 to 4.0 THz. Furthermore, the clipped lines become less
distinct, when the T-ray power reaches the noise floor beyond
3.0 THz.
E. Continuum absorption
The continuum absorption, unaccounted for so far, is defined
as the excess measured absorption unable to be quantified
by the resonance spectrum. The mechanism underpinning the
continuum has not been fully understood [36], and most of the
mathematical models are fitted to experimental observations.

4

It is known that the continuum absorption is pressure and
frequency dependent. A continuum for H2 O collisions, derived
from an experiment in the T-ray frequency range, between 0.4
and 1.83 THz, is estimated at 4.22×10−8 (dB/km)/(hPa GHz)2
[36]. A normal atmospheric measurement in a relatively short
path length is insignificantly influenced by the continuum
absorption [36], [37].
IV. R EMOVAL OF H2 O RESPONSE BY STRENGTH TUNING
This section proposes a numerical means for removal of
water-vapor effects, as an alternative to the conventional
‘sealed chamber’ method. Subsection IV-A establishes the Tray spectroscopic system, and points out the difficulties in using a direct deconvolution to remove the water vapor response.
The strength-tuning algorithm is introduced in Subsection IVB, followed by the fluctuation ratio as a tuning criteria in
Subsection IV-C. Discussion on the algorithm’s generality can
be found in Subsection IV-D.
A. Water vapor as a black box
During a THz-TDS measurement, an emitted T-ray signal
evolves based on several factors, e.g. the sample response,
the water vapor response, the instrument response, and the
noise. These can be modelled as a system, as illustrated in
Fig. 3. Given that x(t) denotes the input T-ray pulse, which
is immediately deployed from a transmitter, s(t) the impulse
response of a sample, w(t) the impulse response of water
vapor, u(t) the impulse response of the instrument, n(t) the
noise, and r(t) the time windowing, the received pulse, y(t),
can be described as
y(t)

=

[x(t) ∗ s(t) ∗ w(t) ∗ u(t)] · r(t)
+[n(t) ∗ u(t)] · r(t) ,

(4)

where ∗ is the convolution operator. The noise is simply
additive and can thus be treated as an extra term in Equation 4.
Via Fourier transform, the above equation in the frequency
domain is
Y (ω) = [X(ω) · S(ω) · W (ω) · U (ω)] ∗ R(ω) + Ń (ω) ,

(5)

where Y (ω), X(ω), S(ω), W (ω), U (ω), R(ω), and Ń (ω)
represent complex frequency responses, as the Fourier pairs
of y(t), x(t), s(t), w(t), u(t), r(t), and [n(t) ∗ u(t)] · r(t)
respectively.
N(ω)

X(ω)

S(ω)

W(ω)

sample
response

water-vapor
response

∑

U(ω)

*

additive
noise

instrument
response

window
convolution

Here, L denotes the T-ray propagation length in free space,
excluding the space occupied by sample(s). The water vapor’s
complex index of refraction, n̂vap = nvap − jκvap , contains
the refractive index nvap and the extinction coefficient κvap ,
which are modeled in Equation 3. The water-vapor response
is also determined by the humidity [38], where ρ(T ) is the
vapor density at temperature T , ρ0 (T ) is the saturation vapor
density at the same temperature.
In order to remove the effects of water vapor from a
measurement, the water-vapor response, W (ω), in Equation 5
must be replaced by the vacuum response, V (ω), with the
same propagation length,
V (ω)

Y(ω)

Fig. 3.
Model of a spectroscopic system. The system response is the
combination of the frequency responses of the sample, water vapor, and
instrument, denoted by S(ω), W (ω), and U (ω), respectively. The symbols,
X(ω), Y (ω), N (ω), and R(ω) represent input T-ray spectrum, measured
spectrum, noise, and window (or limited recording duration) responses,
respectively. The water-vapor response in the dotted box needs to be removed.

=

exp [−jωL/c] .

(7)

From Equation 5 the noisy response can be modelled in two
parts,

 [X(ω)·S(ω)·W (ω)·U (ω)] ∗ R(ω)
, if |X ·S ·W | > |N |
(8)
Y (ω) ≈

Ń (ω)
, otherwise.
Assuming the noise level is unity and the model water-vapor
response matches well with the measurement, performing the
direct deconvolution would yield

 [X(ω) · S(ω) · V (ω) · U (ω)] ∗ R(ω)
, if |X · S · W | > 1
(9)
Y (ω) ≈

V (ω)/W (ω)
, otherwise.
Therefore, removal of the water vapor response fails at
|X · S · W | ≤ |N | or D(ω) ≤ |S · W |−1 , where D(ω) =
|X/N | is the system’s dynamic range. This usually occurs
at the peaks of absorption resonances, as shown in Fig. 2.
Consider the water vapor response closely. It can be separated into several components, and the separation is formulated
as
W (ω)

=

V (ω)W0 (ω)W1 (ω) . . . Wa (ω) .

(10)

Each component corresponds to a complex resonance, and has
its frequency response derived from Equations 3 and 6, or
Wa (ω)

R(ω)

+
+

The water-vapor frequency response can be expressed by


ρ(T ) ωL
W (ω) = exp −j n̂vap
·
.
(6)
ρ0 (T ) c

=

exp [−j ma (na − jκa )ω/c] .

(11)

The line strength, ma , now incorporates the humidity ratio,
ρ(T )/ρ0 (T ), and the propagation length, L. It can be inferred
from Equations 10 and 11 that all components or resonances
can be removed from the measured spectrum separately and
independently. Hereby, the deconvolution at the noisy parts of
the spectrum can be avoided, and also it is possible to tune
ma of each resonance to the measurement without accurate
information of the atmospheric and geometric parameters.
In order to achieve this parameter relaxation, somehow, a
criterion indicating the fit of tuning is essential.

WITHAYACHUMNANKUL et al.: NUMERICAL REMOVAL OF WATER-VAPOR EFFECTS FROM THZ-TDS MEASUREMENTS

B. Fluctuation-removal algorithm

C. Fluctuation ratio
A criterion related to the quality of a T-ray signal is
necessary in selecting the optimal resonance line strengths.
The ‘quality’ here is defined as having a high pulse energy
with low fluctuations in the time domain, which implies the
absence of water-vapor resonances in the frequency domain.
In order to quantify the quality one must be able to evaluate
the total energies of the main pulse and of the fluctuations.
One potential way is to window the time-domain signal with
a Gaussian profile to obtain the amplitude at a desired portion.
A Gaussian window is selected, because a T-ray main pulse
is essentially derived from an optical pumping pulse, which
takes on a Gaussian pulse profile [39].
Fig. 5(a) shows a Gaussian window, gσ (t), overlapping the
T-ray signal, y(t). A Gaussian window with appropriate width
and position would eliminate the fluctuating tail of a signal, as
shown in Fig. 5(b). On the other hand, the complement of a
Gaussian window could also be used to remove the main pulse,
as in Fig. 5(c). With the assistance of a Gaussian window,

START

reset a, b
a for line position
b for line strength
model
resonance profile
at ωa

H2O
line parameters
from database

multiply
resonance
with ma,b

fine-tuning of resonance strength

The fluctuation-removal algorithm is shown in Fig. 4. The
counters a and b for the line position and line strength are
reset. The H2 O spectral line parameters, including intrinsic
line positions and strengths, in the frequency range of interest, are then fetched from an existing database. A complex
resonance profile, n(ω) − jκ(ω), of a rotational transition
at ωa=0 is modeled according to the Lorentzian profile in
Equation 2. The modeled profile multiplied by the initial
strength factor ma=0,b=0 is then temporarily deconvolved
from the measured complex response, Y (ω). The deconvolved
time-domain signal, y(t), is estimated for its fluctuation ratio
(see Subsection IV-C). The procedure repeats to find the
fluctuation ratio at other predefined line strengths ma=0,b ;
b ∈ {1, 2, 3, . . .}.
Once the tuning range, ma=0,b , is covered, the algorithm
picks up an optimal strength within ma=0,b , which gives the
minimum fluctuation ratio, if any, and permanently removes
that optimal complex resonance from the measured signal.
The tuning procedure starts again, but with the next transition
resonance, ωa=1 , and is repeated until all resonances are optimized and possibly removed. It should be noted that a strong
resonance, which suffers from noise i.e. |X · S · W | < |N |
at and around the peak of resonance, will be avoided by the
tuning process, because performing deconvolution of such an
ill-defined resonance will give rise to a high fluctuation ratio.
Since there might be dependencies among the resonances, with
regard to time-domain fluctuations, the whole process starts
again until the fluctuation ratio no longer decreases.
It is advisable not to perform zero padding prior to Fourier
transform in the removal process, as the interpolated spectrum
does not exactly reflect the reality—this is especially the case
for points at resonances that have high variability. Otherwise,
discrepancy between the model and the spectrum would cause
a large remnant in the spectrum after deconvolution. Also note
that to speed up the algorithm, the resonances that have a
strength lower than the amplitude uncertainty can be skipped.

5

deconvolve
resonance
from Y(ω)
calculate
fluctuation ratio
of y(t)

strength
tuning done
?
YES

YES

increase stength
b = b+1

NO

reset strength
b=0
increase frequency
a = a+1

optimum line
found ?

remove
optimum resonance
from Y(ω)

NO
all resonances
checked ?

NO

YES
YES

fluctuation ratio
decreases ?
NO
END

Fig. 4. Fluctuation-removal algorithm. Each resonance in the T-ray frequency
range is modeled and fine-tuned in its strength with the criterion of minimum
fluctuation, and then the optimal resonance is removed from the measured
complex response. The procedure is repeated until all resonances within the
frequency range of interest are investigated, and until the fluctuation ratio no
longer decreases.

the fluctuation energy normalised by the main pulse’s total
energy—coined the fluctuation ratio—can be formulated as
R
[y(t) · (1 − gσ (t))]2 dt
t R
F (y, g ; σ) =
.
(12)
[y(t) · gσ (t)]2 dt
t
The integration is carried out over the time duration of a
recorded T-ray signal. A Gaussian window is given by
gσ (t)

=

exp[−(t − t0 )2 /2σ 2 ] ,

(13)

where t0 is the √
peak position of the T-ray signal, y(t), and σ
multiplied by 2 2 ln 2 is the FWHM of a Gaussian window.
D. Generality of the algorithm
The algorithm is plausibly general—in the sense that it can
be applied to any T-ray signal with minimal disturbance of
desired signal features. This generality is due to the following
facts:
First of all, a T-ray spectrum is altered by the algorithm
only the frequencies at which water-vapor absorption lines

6

1

y(t)

0.5

2.0

gσ(t)

1.5

0

nitrogen atmosphere

t0

(a) fluctuating T-ray signal and a Gaussian window

1

y(t)·gσ(t)

0.5
0

amplitude (a.u.)

amplitude (a.u.)

-0.5

-0.5
(b) T-ray signal weighted by the Gaussian window
0.1

0

10

20

30
time (ps)

40

50

water-vapor atmosphere

0.0

60

(c) T-ray signal weighted by the Gaussian window complement
Fig. 5. Gaussian window applied to the T-ray signal. A T-ray signal weighted
by a Gaussian window, which has its peak position set exactly at the main
pulse peak, yields the main pulse with suppressed fluctuations. A T-ray pulse
weighted by the Gaussian window complement yields only fluctuations. The
FWHM of the Gaussian window used here is 8.33 ps.

are situated. The spectral resonances of other polar gases,
occupying narrow frequency bands and usually not overlapping with water resonances, are not disturbed. In addition, the
broad spectral features of a time-domain transient is unlikely
to be significantly affected by narrow water line removal. The
reflections, for example, which exhibit fringes over a broad
frequency range, are not affected by narrow band resonance
removal.
Second, unrelated fluctuations are ignored. Because the
algorithm senses the fluctuation change, caused by strength
tuning, rather than the fluctuation itself, any unrelated fluctuations cannot deceive the algorithm.
Last, the causality of a signal is preserved. The absorption
and dispersion profiles, used in the model of water vapor
response, exactly comply with the Kramers-Krönig relation.
Deconvolution of the measured spectrum by this causal response would yield the causal spectrum.
V. R ESULTS
A. Free-path measurement
The first T-ray signal, subject to the water-vapor removal
algorithm, is measured in free space without the presence of
any material. The signal has the temporal resolution of 66.7 fs
and the total duration of 136.63 ps, providing the spectrum
with the spectral resolution of 7.3 GHz.
The strength-tuning algorithm is carried out with a set
of complex resonances in the frequency range between 0
and 4 THz, where the T-ray magnitude is relatively high.
The sequence, in which the resonances are interrogated and
removed, follows the numerical order, i.e. from low to high
frequencies. Only the resonances that have the strength higher
than 0.01 of the maximum strength are inspected. The FWHM
of the Gaussian window is 3 ps, and the iteration of the
algorithm is set to five.

time (ps)

5

10

-0.5
water-vapor atmosphere with fluctuation removal

-1.5
0

10

20

30

40
50
time (ps)

60

70

80

Fig. 6. T-ray signals measured in free path. The sampling interval is 66.7 fs
with the total duration of 136.63 ps (only the first 80 ps is shown here). The
inset shows a zoom-in of signals between 5 and 10 ps. The fluctuation ratios
of the N2 -atmosphere, H2 O-atmosphere, and H2 O-removal measurements
are 2.2, 9.2, and 3.3, respectively. The fluctuation after the main pulse is
remarkably reduced.

As shown in Fig. 6, the algorithm can significantly reduce
the time-domain fluctuations, which are previously located
immediately after the main pulse, i.e. after 10 ps. The main
pulses in the inset clearly illustrate the similarity between
the fluctuation-free signal and the processed signal, despite
a temporal shift which is not accounted for by the algorithm.
In Fig. 7, the spectrum result demonstrates success of the algorithm in removing, or reducing the strength of, the resonances
in the frequency range from 0 to 2.5 THz. However, some
resonances still persist, where they are spectrally congested or
severely disturbed by noise.
B. DL-Phenylalanine measurement
The second experiment is performed with a signal measured
from a DL-Phenylalanine sample [40] at 220 K. The cryostat containing the sample is evacuated, so no water vapor
104

magnitude (a.u.)

-0.1

0.5

-1.0

y(t)·(1-gσ(t))

0

1.0

nitrogen atmosphere
water-vapor atmosphere
with fluctuation removal
water-vapor atmosphere

103

102

101

1
0

0.5

1

1.5
2
2.5
frequency (THz)

3

3.5

4

Fig. 7.
T-ray spectra measured in free path. The spectral resolution is
7.3 GHz. Up to 2.5 THz most of the water resonances are removed, or at
least reduced in the strength, by the proposed algorithm. The magnitudes are
offset for clarity.

WITHAYACHUMNANKUL et al.: NUMERICAL REMOVAL OF WATER-VAPOR EFFECTS FROM THZ-TDS MEASUREMENTS

2.0

7

104

nitrogen atmosphere
water-vapor atmosphere
with fluctuation removal
water-vapor atmosphere

nitrogen atmosphere

1.0
0.5

water-vapor atmosphere
0

magnitude (a.u.)

amplitude (a.u.)

1.5
103

102

101

-0.5
water-vapor atmosphere with fluctuation removal
-1.0
A

-1.5
0

10

20

30

40

50

B
60

1
0

0.5

1

1.5

time (ps)
Fig. 8. T-ray signals measured with DL-Phenylalanine sample in place.
The sampling interval is 66.7 fs with the total duration of 68.38 ps . The
reflections at 32 and 63 ps due to the cryostat’s windows can be recovered
from fluctuation by the algorithm. The fluctuation ratios of the N2 -atmosphere,
H2 O-atmosphere, and H2 O-removal measurements are 2.6, 7.9, and 3.8,
respectively.

at cryogenic temperature involves in the measurement. The
interacting water vapor is from normal air, surrounding the
cryostat. The total scan is 68.3 ps at the step size of 66.7 fs.
The spectral resolution is 14.6 GHz.
The recorded signal particularly contains the reflections,
caused by the cryostat’s windows, made from a cyclo-olefin
copolymer (Topas). Use of such a signal demonstrates how
well the proposed algorithm can deal with water-vapor fluctuations without disturbing unrelated reflection features.
The parameters set for the algorithm are similar to the
previous test, except for the interrogated frequency range
which is in between 0 and 2 THz.
The processed signal is shown in Fig. 8, along with the
original. Again, it is obvious that the fluctuations located after
the main pulse, beyond 15 ps, are remarkably reduced, in spite
of an introduction of small oscillation at the beginning. The
reflection at 63 ps (point B) is not disturbed by the algorithm,
and, surprisingly, the reflection at 32 ps (point A), which was
initially buried in fluctuations, is now recovered. In Fig. 9, the
algorithm can remove most of the H2 O resonances. Still a few
resonances cannot be completely removed. These persistent
resonances are around 1.1 and 1.7 THz, exactly the same
positions as the unremoved resonances in the previous case
(see Fig. 7). The recovery of reflections suggests a promising
application of the algorithm in T-ray range finding.
C. Evaluation of the algorithm’s performance
Table I shows the fluctuation ratio and mean squared error
(MSE) of the signals and the total energy of the spectra, investigated in the previous subsections. Two additional results,
not shown in this paper, are tested with lactose’s reference and
sample signals. The fluctuation ratio helps measure a change
in the fluctuation energy and main pulse energy, before and
after the signals are processed. The MSE evaluates the fitness
of the improved signals to the N2 -atmosphere measurement.

2

2.5

3

3.5

4

frequency (THz)
Fig. 9. T-ray spectra measured with DL-Phenylalanine sample in place. The
spectral resolution is 14.6 GHz. From 0 to 2 THz most of the resonances are
removed, or at least reduced in the strength, by the proposed algorithm. The
magnitudes are offset for clarity.
TABLE I
Q UANTITATIVE EVALUATION OF THE ALGORITHM ’ S PERFORMANCE . T HE
NOTATION H2 O, AND H2 O STAND FOR MEASUREMENTS IN WATER VAPOR
AND WATER VAPOR WITH FLUCTUATION REMOVAL , RESPECTIVELY.

Signal
Free path
Phenylalanine
Lactose reference
Lactose sample

Fluctuation ratio

MSE (%)

Energy (a.u.)

H2 O

H2 O

H2 O

H2 O

H2 O

H2 O

4.2
3.0
16.1
4.8

1.5
1.5
2.0
1.7

0.076
0.13
0.31
0.31

0.047
0.085
0.19
0.21

0.90
0.90
0.90
0.93

0.96
0.95
0.98
0.99

The MSE is given by
MSE

=

1 X
(ým − ym )2 ,
m m

(14)

where y is a signal measured in the nitrogen atmosphere and
ý is a compared signal. The error is summed and averaged
over m temporal points. The total spectrum energy shows the
recovery of the energy from resonance absorptions. Note that
the fluctuation ratio and the spectral energy are normalised to
those of the signal measured in nitrogen atmosphere.
From the table the fluctuation ratio indicates that the fluctuation energy is greatly reduced in all cases, once the removal
algorithm is implemented with the signal measured in the H2 O
path. However, the MSE still reflects a considerable difference
between the processed signals and the signals measured in
a nitrogen-filled chamber, in spite of the great improvement
visualized in the previous subsections. This would probably
be due to the whole signal shift in the time domain, caused
by a constant refractive index difference unaccounted for by
the algorithm.
The total spectral energy in Table I reveals the algorithm is
successful in recovering a part of the signal’s energy, which
is absorbed by water vapor resonances during propagation.
However, not all the energy is recoverable. With regard to the
processed spectra shown in the previous subsections, the algorithm cannot completely remove the water vapor resonances
in the region where the SNR is low and where the spectral

8

lines are crowded. The causes of these two cases are analyzed
separately as follows.
At a low SNR region, distortions in the absorption and
dispersion resonances by noise causes deviation of the measurement from the model. This situation is demonstrated in
Fig. 2. Applying the removal algorithm to a low SNR region
would be effective, provided that the distortion is not severe.
Otherwise, the large deviation inhibits the removal via a suboptimal value of the fluctuation ratio. A possible means to
alleviate the noise issue is to implement a signal enhancement
methodology such as wavelet denoising [41] prior to the
removal process; this might better recover water resonances
and subsequently improve the removal efficiency.
Crowded resonance features result in overlapping and merging between two or more resonances. If the spectral resolution
is too low, these blended resonances are rendered incorrectly,
resulting in distorted shape. The algorithm would see them
as a single stronger and wider resonance, and as a result, the
removal does not perform well. Improving the spectral resolution would allow better removal performance. For example, the
free space measurement in Fig. 7 shows a better result than the
DL-Phenylalanine measurement in Fig. 9 in the regions around
1.1 and 1.7 THz; the former and latter measurements have the
spectral resolutions of 7.3 GHz and 14.6 GHz, respectively. In
addition to the spectral resolution problem, blended resonances
cannot be represented exactly by a simple sum between two
Lorentzian lines [42].
VI. C ONCLUSION
A numerical algorithm for removal of water-vapor effects in
T-ray measurements is investigated in this paper. Given only a
sample signal with no reference, theoretically, the removal can
possibly be carried out by simple deconvolution of the model
complex water-vapor response from the signal. However, many
factors limit the usability of the simple deconvolution. An
exact model of water vapor resonance requires many physical data, including the pressure, temperature, humidity, and
propagation length. In addition, if the noise level is sufficiently
high such that some strong resonances are distorted, the whole
deconvolution process cannot be performed.
The proposed algorithm fine-tunes the strength of each
model spectral resonance. A criterion for strength tuning is met
when the fluctuation ratio of the processed signal reaches the
minimum. Once an optimal resonance is attained, it is removed
from the signal. The algorithm then proceeds to removal of the
next resonance. This tuning scheme relaxes the requirement for
precise information of the delicate physical conditions during
the measurement. Furthermore, the fluctuation ratio criterion
inhibits any fault deconvolution that might occur when the
noise disturbs the quality of a measured signal.
In the experiments, the algorithm produces promising results. It can reduce a significant amount of the signal fluctuations and spectral resonances with small disturbances to other
non-related features, such as reflections or sample-induced resonances. Moreover, some features, which are initially masked
by fluctuations, can be recovered by the algorithm. Optical
constants extracted from improved signals have an acceptable
quality.

The proposed algorithm is general in the sense that any
other polar gas response could be targeted, in principle, and
hence removed from a measured spectrum if desired. It is not
necessary that the parameters of a targeted gas be present in a
spectroscopic catalog, as long as its pure response in the frequency range of interest is fully estimable. This scheme might
have benefits in some particular situations, for instance, where
molecular contaminant(s) are unavoidable in measurement.
Finally, it should be noted that the algorithm is not meant to
serve as a full replacement to the ‘sealed chamber’ procedure
executed during the measurement stage. In fact, it assists the
measurement in cases where a sealed chamber is not feasible.
Still the algorithm cannot remove some traces of resonances
due to the ubiquitous noise and resonance crowding. Further
improvements of the algorithm remains an attractive challenge.
ACKNOWLEDGMENT
The authors gratefully acknowledge Morten Franz at the
Department of Molecular and Optical Physics, the University
of Freiburg, for his support in T-ray measurements, and Brian
W.-H. Ng at the School of Electrical & Electronic Engineering,
the University of Adelaide, for useful technical discussion.
R EFERENCES
[1] D. M. Mittleman, M. Gupta, R. Neelamani, R. G. Baraniuk, J. V. Rudd,
and M. Koch, “Recent advances in terahertz imaging,” Applied Physics
B: Lasers and Optics, vol. 68, no. 6, pp. 1085–1094, 1999.
[2] W. Withayachumnankul, G. Png, X. X. Yin, S. Atakaramians, I. Jones,
B. U. H. Lin, J. Balakrishnan, B. W.-H. Ng, B. Ferguson, S. P. Mickan,
B. Fischer, and D. Abbott, “Terahertz sensing and imaging,” Proceedings
of the IEEE, 2007, accepted for publication.
[3] F. Huang, B. Schulkin, H. Altan, J. F. Federici, D. Gary, R. Barat,
D. Zimdars, M. Chen, and D. B. Tanner, “Terahertz study of 1,3,5trinitro-s-triazine by time-domain and Fourier transform infrared spectroscopy,” Applied Physics Letters, vol. 85, no. 23, pp. 5535–5537, 2004.
[4] D. M. Mittleman, R. H. Jacobsen, R. Neelamani, R. G. Baraniuk, and
M. C. Nuss, “Gas sensing using terahertz time-domain spectroscopy,”
Applied Physics B: Lasers and Optics, vol. 67, no. 3, pp. 379–390, 1998.
[5] G. Tripoli, Handbook of Weather, Climate, and Water. John Wiley &
Sons, 2003, ch. 15 - Physical Atmospheric Science, pp. 177–206.
[6] M. van Exter, C. Fattinger, and D. Grischkowsky, “Terahertz timedomain spectroscopy of water vapor,” Optics Letters, vol. 14, no. 20,
pp. 1128–1130, 1989.
[7] D. Grischkowsky, S. Keiding, M. van Exter, and C. Fattinger, “Farinfrared time-domain spectroscopy with terahertz beams of dielectrics
and semiconductors,” Journal of the Optical Society of America B:
Optical Physics, vol. 7, no. 10, pp. 2006–2015, 1990.
[8] J. F. Federici, B. Schulkin, F. Huang, D. Gary, R. Barat, F. Oliveira,
and D. Zimdars, “THz imaging and sensing for security applications–
explosives, weapons and drugs,” Semiconductor Science and Technology,
vol. 20, pp. S266–S280, 2005.
[9] H.-B. Liu, H. Zhong, N. Karpowicz, Y. Chen, and X.-C. Zhang, “Terahertz spectroscopy and imaging for defense and security applications,”
Proceedings of the IEEE, 2007, (In Press).
[10] P. U. Jepsen and B. M. Fischer, “Dynamic range in terahertz timedomain transmission and reflection spectroscopy,” Optics Letters,
vol. 30, no. 1, pp. 29–31, 2005.
[11] J. Xu, T. Yuan, S. Mickan, and X. C. Zhang, “Limit of spectral resolution
in terahertz time-domain spectroscopy,” Chinese Physics Letters, vol. 20,
no. 8, pp. 1266–1268, 2003.
[12] W. Withayachumnankul, H. Lin, S. P. Mickan, B. M. Fischer, and
D. Abbott, “Analysis of measurement uncertainty in THz-TDS,” in Proc.
SPIE Photonic Materials, Devices, and Applications II, vol. 6593, 2007,
p. art. no. 659326.
[13] A. J. Kemp, J. R. Birch, and M. N. Afsar, “The refractive index of water
vapour: a comparison of measurement and theory,” Infrared Physics,
vol. 18, pp. 827–833, 1978.

WITHAYACHUMNANKUL et al.: NUMERICAL REMOVAL OF WATER-VAPOR EFFECTS FROM THZ-TDS MEASUREMENTS

[14] P. F. Bernath, “The spectroscopy of water vapour: Experiment, theory
and applications,” Physical Chemistry Chemical Physics, vol. 4, pp.
1501–1509, 2002.
[15] S. D. Gasster, C. H. Townes, D. Goorvitch, and F. P. J. Valero, “Foreigngas collision broadening of the far-infrared spectrum of water vapor,”
Journal of the Optical Society of America B: Optical Physics, vol. 5,
pp. 593–601, 1988.
[16] C. H. Townes and A. L. Schawlow, Microwave Spectroscopy. McGrawHill, 1955.
[17] H. Harde, R. Cheville, and D. Grischkowsky, “Terahertz studies of
collision-broadened rotational lines,” Journal of Physical Chemistry A,
vol. 101, no. 20, pp. 3646–3660, 1997.
[18] J. Kauppinen, T. Kärkkäinen, and E. Kyrö, “High-resolution spectrum
of water vapor between 30 and 720 cm−1 ,” Journal of Molecular
Spectroscopy, vol. 71, no. 1-3, pp. 15–45, 1978.
[19] J. K. Messer, F. C. D. Lucia, and P. Helminger, “The pure rotational
spectrum of water vapor—a millimeter, submillimeter, and far infrared
analysis,” International Journal of Infrared and Millimeter Waves, vol. 4,
no. 4, pp. 505–539, 1983.
[20] R. A. Cheville and D. Grischkowsky, “Far-infrared foreign and selfbroadened rotational linewidths of high-temperature water vapor,” Journal of the Optical Society of America B: Optical Physics, vol. 16, no. 2,
pp. 317–322, 1999.
[21] P. Helminger, J. K. Messer, and F. C. D. Lucia, “Continuously
tunable coherent spectroscopy for the 0.1–1.0-THz region,” Applied
Physics Letters, vol. 42, no. 4, pp. 309–310, 1983. [Online]. Available:
http://link.aip.org/link/?APL/42/309/1
[22] F. Matsushima, H. Odashima, T. Iwasaki, S. Tsunekawa, and K. Takagi,
“Frequency measurement of pure rotational transitions of H2 O from 0.5
to 5 THz,” Journal of Molecular Structure, vol. 352-353, pp. 371–378,
1995.
[23] P. Chen, J. C. Pearson, H. M. Pickett, S. Matsuura, and G. A. Blake,
“Submillimeter-wave measurements and analysis of the ground and ν2
= 1 states of water,” The Astrophysical Journal Supplement Series, vol.
128, pp. 371–385, 2000.
[24] H. M. Pickett, R. L. Poynter, E. A. Cohen, M. L. Delitsky, J. C. Pearson,
and H. S. P. Müller, “Submillimeter, millimeter, and microwave spectral line catalog,” Journal of Quantitative Spectroscopy and Radiative
Transfer, vol. 60, no. 5, pp. 883–890, 1998.
[25] R. L. Poynter and H. M. Pickett, “Submillimeter, millimeter, and
microwave spectral line catalog,” Applied Optics, vol. 24, no. 14, pp.
2235–2240, 1985.
[26] V. B. Podobedov, D. F. Plusquellic, and G. T. Fraser, “THz laser study
of self-pressure and temperature broadening and shifts of water vapor
lines for pressures up to 1.4 kPa,” Journal of Quantitative Spectroscopy
and Radiative Transfer, vol. 87, pp. 377–385, Sept. 2004.
[27] P. F. Bernath, Spectra of Atoms and Molecules, 2nd ed.
Oxford
University Press, 2005.
[28] N. Gopalsami, S. Bakhtiari, A. C. Raptis, S. L. Dieckman, and F. C.
DeLucia, “Millimeter-wave measurements of molecular spectra with
application to environmental monitoring,” IEEE Transactions on Instrumentation and Measurement, vol. 45, no. 1, pp. 225–230, 1996.
[29] W. S. Benedict and L. D. Kaplan, “Calculation of line widths in H2 OH2 O and H2 O-O2 collisions,” Journal of Quantitative Spectroscopy and
Radiative Transfer, vol. 4, no. 3, pp. 453–469, 1964.
[30] L. S. Rothman, D. Jacquemart, A. Barbe, D. C. Benner, M. Birk,
L. R. Brown, M. R. Carleer, C. Chackerian, K. Chance, L. H. Coudert,
V. Dana, V. M. Devi, J. M. Flaud, R. R. Gamache, A. Goldman, J. M.
Hartmann, K. W. Jucks, A. G. Maki, J. Y. Mandin, S. T. Massie,
J. Orphal, A. Perrin, C. P. Rinsland, M. A. H. Smith, J. Tennyson,
R. N. Tolchenov, R. A. Toth, J. Vander Auwera, P. Varanasi, and
G. Wagner, “The HITRAN 2004 molecular spectroscopic database,”
Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 96,
pp. 139–204, 2005.
[31] F. C. De Lucia, Sensing with Terahertz Radiation, ser. Springer Series
in Optical Sciences. Springer, 2003, vol. 85, ch. Spectroscopy in the
terahertz spectral region, pp. 39–115.
[32] H. A. Lorentz, “The absorption and emission lines of gaseous bodies,”
Proceedings of Royal Academy of Amsterdam, vol. 8, pp. 591–611, 1906.
[33] J. H. Van Vleck and V. F. Weisskopf, “On the shape of collisionbroadened lines,” Reviews of Modern Physics, vol. 17, no. 2-3, pp. 227–
236, 1945.
[34] R. Emery, “Atmospheric absorption measurements in the region of 1
mm wavelength,” Infrared Physics, vol. 12, no. 1, pp. 65–79, 1972.
[35] W. Withayachumnankul, B. Ferguson, T. Rainsford, S. P. Mickan, and
D. Abbott, “Material parameter extraction for terahertz time-domain
spectroscopy using fixed-point iteration,” in Proc. SPIE Photonic

[36]

[37]
[38]

[39]

[40]

[41]
[42]

9

Materials, Devices, and Applications, vol. 5840, no. 1, 2005, pp.
221–231. [Online]. Available: http://link.aip.org/link/?PSI/5840/221/1
V. B. Podobedov, D. F. Plusquellic, and G. T. Fraser, “Investigation of
the water-vapor continuum in the THz region using a multipass cell,”
Journal of Quantitative Spectroscopy and Radiative Transfer, vol. 91,
pp. 287–295, Mar. 2005.
S. A. Clough, F. X. Kneizys, and R. W. Davies, “Line shape and the
water vapor continuum,” Atmospheric Research, vol. 23, no. 3-4, pp.
229–241, 1989.
T. Yuan, H. Liu, J. Xu, F. Al-Douseri, Y. Hu, and X.-C. Zhang, “Terahertz time-domain spectroscopy of atmosphere with different humidity,”
in Proc. SPIE Terahertz for Military and Security Applications, R. J.
Hwu and D. L. Woolard, Eds., vol. 5070, 2003, pp. 28–37.
L. Duvillaret, F. Garet, J.-F. Roux, and J.-L. Coutaz, “Analytical modeling and optimization of terahertz time-domain spectroscopy experiments
using photoswitches as antennas,” IEEE Journal of Selected Topics in
Quantum Electronics, vol. 7, no. 4, pp. 615–623, 2001.
M. Yamaguchi, F. Miyamaru, K. Yamamoto, M. Tani, and M. Hangyo,
“Terahertz absorption spectra of L-, D-, and DL-alanine and their
application to determination of enantiometric composition,” Applied
Physics Letters, vol. 86, no. 5, p. art. no. 053903, 2005.
B. Ferguson and D. Abbott, “De-noising techniques for terahertz responses of biological samples,” Microelectronics Journal, vol. 32, no. 12,
pp. 943–953, 2001.
H. M. Pickett, “Determination of collisional linewidths and shifts by a
convolution method,” Applied Optics, vol. 19, no. 16, pp. 2745–2749,
1980.

Withawat Withayachumnankul was born in
Bangkok, Thailand on December 20, 1980. From
1996-2000, he was educated at King Mongkut’s
Institute of Technology at Ladkrabang (KMITL),
in his hometown, where he obtained a Bachelor’s
Degree in Electronics Engineering with honours.
Working at the Biomedical Signal & Image Processing Laboratory at KMITL in 2000, as a postgraduate student, he performed multi-disciplinary
research involving digital signal processing, medical
image reconstruction, and computer graphics. Here,
he played a major role in developing core programming codes for medical
software currently deployed in many leading Thai institutions and hospitals.
After receiving a Master’s degree in Electronics Engineering in 2002, he
served at his alma mater as a lecturer, a position that he currently retains.
In this role, he conducted the computer programming and computer graphics
courses in the Department of Information Engineering.
In 2005, he was a visiting scholar at the University of Adelaide, under Derek
Abbott. In 2006, Withawat Withayachumnankul was granted an Australian
Endeavour International Postgraduate Research Scholarship (EIPRS) to study
toward his PhD under Derek Abbott, Bernd Fischer, and Samuel Mickan,
within the Adelaide T-ray group, the School of Electrical & Electronic
Engineering, The University of Adelaide. Mr Withayachumnankul’s current
research involves T-ray signal processing, T-ray components, and superluminal
(group delay) electromagnetic propagation.

Bernd M. Fischer was born in Waldkirch, Germany.
In 1997/98 he won an ERASMUS bursary and
studied for one year in the Université Paris-sud 11,
France, where he completed a Stage de Maitrise in
biophysics and he received his Diplom Degree (with
distinction) in 2001 from the University of Freiburg,
Germany. He received his PhD degree (summa cum
laude) in physics from the University of Freiburg,
Germany, in 2006, under Peter Uhd Jepsen and
Hanspeter Helm in the area of T-ray spectroscopy
of biomolecules.
Dr Fischer joined the Adelaide T-ray Group, Australia, in 2006, and in
2007 he was awarded the prestigious Australian Research Council (ARC)
Postdoctoral (APD) Fellowship.

10

Samuel P. Mickan (S’99-M’03) was born on March
26, 1976 in Adelaide, Australia. He graduated with
a BEng in Electrical and Electronic Engineering,
(1997), and a with BA in German & History (2002),
both at the University of Adelaide. He received his
PhD (with Special Commendation) in 2004 under
Derek Abbott and Jesper Munch at the University
of Adelaide, in the area of T-ray biosensing. He
has won a number of awards and honors including, the Schneider Award for Excellence (1997),
AFUW Barbara Crase Bursary (2001), George Murray Scholarship (2002), the Australian Academy of Sciences (AAS) Young
Researcher’s Award (2004), the University Medal from the University of
Adelaide (2004), and the Tall Poppy Award for science (2005). Under a
Fulbright Scholarship, he collaborated with X.-C. Zhang’s group, Center for
Terahertz Technology, Rensselaer Polytechnic Institute in New York, USA,
from 2001 to 2002. In 2003, he began lecturing in the School of Electrical
and Electronic Engineering at the University of Adelaide. He is currently both
an adjunct senior lecturer at the University of Adelaide and a patent attorney
at Davies Collison Cave, Melbourne, Australia. He was Treasurer/Secretary
of Microwave Theory & Techniques Chapter of LEOS and was founder of the
IEEE student branch at the University of Adelaide. Dr Mickan is a member
of the IEEE, IEEE LEOS society, OSA, APS, and SPIE.

Derek Abbott (M’85-SM’99-F’05) was born on
May 3, 1960, in South Kensington, London, UK,
and he received a BSc (Hons) in physics (1982)
from Loughborough University of Technology, UK.
He completed his PhD (with commendation), in
electrical & electronic engineering (1995) from The
University of Adelaide, Australia, under Kamran
Eshraghian and Bruce R. Davis. He has led a number
of research programs in the imaging arena, ranging
from the optical to infrared to millimeter wave to
T-ray (terahertz gap) regimes. From 1978 to 1986,
he worked at the GEC Hirst Research Centre, London, UK, in the area of
visible and infrared image sensors. His expertise also spans VLSI design,
optoelectronics, device physics, and noise; where he has worked with nMOS,
CMOS, SOS, CCD, GaAs, and vacuum microelectronic technologies. On
migration to Australia, he worked for Austek Microsystems, Technology Park,
South Australia, in 1986. Since 1987, he has been with The University of
Adelaide, where he is presently a full Professor in the School of Electrical
& Electronic Engineering and the Director of the Centre for Biomedical
Engineering (CBME). He has appeared on national and international television
and radio and has also received scientific reportage in New Scientist, The
Sciences, Scientific American, Nature, The New York Times, and Sciences
et Avenir. He holds over 300 publications/patents and has been an invited
speaker at over 80 institutions, including Princeton, NJ; MIT, MA; Santa Fe
Institute, NM; Los Alamos National Laboratories, NM; Cambridge, UK; and
EPFL, Lausanne, Switzerland. He won the GEC Bursary (1977), the Stephen
Cole the Elder Prize (1998), the E.R.H. Tiekink Memorial Award (2002),
SPIE Scholarship Award for Optical Engineering and Science (2003), the
South Australian Tall Poppy Award for Science (2004) and the Premier’s
SA Great Award in Science and Technology for outstanding contributions to
South Australia (2004). He has served as an editor and/or guest editor for a
number of journals including IEEE J OURNAL OF S OLID -S TATE C IRCUITS,
Chaos (AIP), Smart Structures and Materials (IOP), Journal of Optics B
(IOP), Microelectronics Journal (Elsevier), Fluctuation Noise Letters (World
Scientific), and is currently on the Editorial Board of P ROCEEDINGS OF THE
IEEE. He has served on a number of IEEE technical program committees,
including the IEEE APCCS and the IEEE GaAs IC Symposium. Prof Abbott
is a Fellow of the Institute of Physics (IOP), with honorary life membership,
and is a Fellow of the IEEE.

