2 December 2013 Nondestructive evaluation of crystallized-particle size in lactose-powder by terahertz time-domain spectroscopy
Author Affiliations +
Optical Engineering, 53(3), 031203 (2013). doi:10.1117/1.OE.53.3.031203
Transmission-type terahertz time-domain spectroscopy is applied to evaluate crystallized lactose particle of size below 30 μm, which is far too small compared to the wavelength of incident terahertz (THz)-wave. The THz-absorption spectrum of lactose is successfully deconvoluted by Lorentzian to two spectra with peaks at 17.1  cm[sup]−1 (0.53 THz) and 45.6  cm 1 (1.37 THz) derived from α -lactose monohydrate, and a spectrum at 39.7  cm1 (1.19 THz) from anhydrous β-lactose after removal of the broad-band spectrum by polynomial cubic function. Lactose is mainly crystallized into α -lactose monohydrate from the supersaturated solution at room temperature with a small amount of anhydrous β-lactose below 4%. The absorption feature is dependent on the crystallized particle size and the integrated intensity ratio of the two absorptions due to α-lactose monohydrate is correlated in linear for the size.
Yamauchi, Hatakeyama, Imai, and Tonouchi: Nondestructive evaluation of crystallized-particle size in lactose-powder by terahertz time-domain spectroscopy



A terahertz (THz)-electromagnetic wave in the far infrared region is useful to characterize material properties based on the intra-vibration in relatively large organic molecules and/or the inter-molecular vibration incorporated hydrogen-bonding. The absorption feature of such organic material is significantly dependent on the molecular and crystal structures as demonstrated for L-phenylalanine comparing to L-tyrosine,1 three different retinal isomers with polyene chain,2 DL-alanine racemic compound comparing to D- and L-alanine,3 and so on. Therefore, the lines peculiar to material can be recognized as the spectral fingerprint. To date, a lot of organic materials have been characterized by the spectral fingerprint which is useful not only for the qualitative analysis but also for the quantitative evaluation. For example, we previously demonstrated that the anomer content of lactose in the milled powder can be precisely determined by the integrated intensity ratio of THz-absorptions due to the α-form and the β-form crystals.4 However, the crystallized particle size of organic materials has not been studied by THz-spectroscopy. In the case of saccharides, such as lactose used in food products, evaluation of the crystallized particle size is often required in sandiness which can be sensorically detected for sizes over 10 to 16 μm. In general, the size is evaluated by optical microscope or scanning electron microscope. Although such methods achieve a precise evaluation on the surface, the particle inside the products cannot be observed because the transmission length of the light or the electron is limited to regions within a submicrometer. In contrast, a THz-wave is notably absorbed in water but easily passed through most organic materials compared to UV–vis–infrared light or electron-beam. In concept, it is considered that the THz-imaging developed in some fields56.7 is useful to disclose the shape of particle inside materials. However, it is not easy to perform the imaging of a particle with a size below 20 μm, which is far too small compared to the wavelength of a THz-wave around 1 THz. On the other hand, if the THz-absorption is dependent on size, THz-spectroscopy is useful to estimate a particle with such a small size.

In this article, transmission-type terahertz time-domain spectroscopy (THz-TDS) is applied to evaluate a lactose powder, which consists of crystallized particles with a size below 30 μm, and precise analysis on the absorption feature is performed to estimate the particle size.



A transmission-type THz-TDS system was used to characterize lactose-powders. A detailed configuration of the system was shown elsewhere.4 In this system, a femto-second fiber laser (peak wavelength at 782 nm, half-width of 87 fs, repetition rate of 48 MHz) was used as a pump light and a probe light after being split in two ways by a half-mirror. The pump light chopped at 1 kHz was focused and irradiated on a THz emitter consisting of a dipole antenna with a 10 μm-gap space fabricated on low-temperature grown GaAs (LT-GaAs) layer and attached on a hemispherical Si-lens, on which the antenna was biased at 10 V to generate a transient current in the pico-second order. A THz-pulse generated on the emitter was radiated through the Si-lens, then focused normally incident to the sample using two off-axis parabolic metal mirrors. After passing through the sample, the THz-pulse was introduced to the detector with the same antenna configuration of the emitter by two off-axis parabolic mirrors and a hemispherical Si-lens. When the THz-pulse was introduced to the antenna in the detector, the probe laser light was simultaneously irradiated to detect the pulse by a sampling technique, where time-delay of the probe light was controlled by a retro-reflector and a micro-step stage controller with a step distance of 1 μm. The sampling data were recorded in a PC after signal-amplification, lock-in noise reduction, and A/D conversion. The recorded pulse data was processed by a discrete Fourier transform (DFT) after a Gaussian-window was superposed on the pulse data to remove aliasing. Temperature in the THz-TDS system was carefully controlled at 20°C.

Commercially available pure lactose α-D(+)-lactose monohydrate [Lα·H2O:Oβ-galactopyranosyl-(14)- α-D-glucopyranose monohydrate (C12H22O11·H2O)] powder (99% total lactose basis, 0.05% glucose Sigma-Aldrich, St. Louis, Missouri) was used to prepare a lactose solution and for seeding. The anomer (Lβ: anhydrous β-D(+)-lactose) content in the Lα·H2O powder was about 4% in the commercial specifications. Lactose particles were extracted from lactose supersaturated solution (1 g-lactose powder/3 cc-water) and the size was controlled by the extract period, where ultra-pure water with the resistivity above 18.2cm was used as a solvent. The crystallization was enhanced by seeding 10 mg-Lα·H2O powder added into the solution at 20°C with stirring (300 rpm), and then the extracted lactose was dried at 60°C for three days in an incubator after removal of the residual solution. The extracted lactose powders consist of the crystallized particles, which are then filled and lightly compressed in a metal aperture of 0.7 mm-thickness with a hole of 6 mm-diameter. A THz-wave was normally incident to the lactose-sample placed with the aperture in the THz-TDS system. The weight of the lactose-powder used for the THz-TDS was measured by an electronic weight-scale.


Result and Discussions


THz-Absorption Feature of Lactose Powder

Figure 1 shows the THz-spectra in the THz-TDS system without a lactose-sample, and the inset shows the THz-pulse and Gaussian-window superimposed on the pulse before DFT. Periodic fluctuation was observed in the THz-spectra (plain line) but successfully removed by using the Gaussian-window (bold-line). If the THz-pulse was sufficiently intense, the fluctuation due to aliasing was negligible, but the pulse-edge in relatively weak signal passed through samples should be taken into account for the precise analysis.

Fig. 1

Fourier transformed THz-spectra from the THz-pulse without Gaussian window (plain-line) and superposed Gaussian window (bold-line). The inset shows the THz-pulse and Gaussian window.


A typical prelimary absorption spectrum of the Lα·H2O powder consists of several narrow-bands and a broad-band as shown in Fig. 2, where the experimentally obtained data were shown by filled-circles with a fit (dash-line), by the deconvoluted spectra (solid-lines), and broad-band (dot-line). The spectral deconvolution was successfully achieved by Lorentzian as previously shown for the lowest lying spectrum8 after removal of a broad-band by a polynomial cubic function. The narrow-bands with peaks at 17.1cm1 (0.53 THz) and 45.6cm1 (1.37 THz) can be recognized as the fingerprint of Lα·H2O,9,10 which were denoted as α1 and α2, respectively. The absorptions derived from Lα·H2O are recognized to be associated to lactose molecular-rotations in the crystal as indicated by first principle calculations.10 A weak absorption with a peak at 39.7cm1 (1.19 THz) denoted as β1 can be suggested to be derived from Lβ because the absorption was dominant in Lβ powder.4

Fig. 2

Absorption spectrum of preliminary used α-lactose powder, where experimentally obtained data, broad band background, deconvoluted spectra and the fit are shown by filled-circle, bold dash-line, solid-line, and plain dash-line, respectively. Absorptions due to α-lactose monohydrate and anhydrous β-lactose are denoted as “α1,” “α2,” and “β1,” respectively.



Crystallization Behavior of Lactose in Supersaturated Solution

Figure 3 shows a photograph of crystallized lactose-particles observed by a Nomarski differential interference microscope, where the lactose was crystallized for (a) 24, (b) 72, and (c) 120 h and then dispersed in ethanol by ultra-sonic vibration for 20 s after being dried in an incubator. It has been recognized that lactose in the supersaturated solution around room-temperature can be mainly crystallized into Lα·H2O by seeding of the Lα·H2O powder. The crystallization behavior is also dependent on this condition because the β-lactose present in the aqueous solution selectively blocks the growth of the (020) face, which results in growth along the b direction (not the b direction) and the Tomahawk shape as shown in Fig. 4.11,12 Further, it was demonstrated by a dimethyl sulfoxide-lactose system with controlled β-lactose content in the solution that the Tomahawk shape is extended along the b direction, reducing the growth along the c direction with increasing β-lactose content in the lactose solution.12 Since β-lactose is present above 60% in aqueous lactose solutions around room-temperature by the mutarotation in equilibrium,13 the crystal is enlarged to the b direction as shown in Fig. 4. In this article, the crystal size was determined along the b axis in the anisotropic shape, and the average crystal size was uniquely increased with the crystallized period as shown in Fig. 5.

Fig. 3

Nomarski photographs of extracted lactose from the supersaturated solution for (a) 24, (b) 72, and (c) 120 h.


Fig. 4

Miller indicates of Tomahawk shaped α-lactose monohydrate crystal.12


Fig. 5

Dependence of average crystallized particle size of α-lactose monohydrate on the crystallized period, where the crystal size was estimated along the b-axis.



Dependence of THz-Absorption Feature on Crystallized Lactose Particle Size

Figure 6 shows the absorption spectra (solid-line) of lactose consisting of crystallized lactose particles with an average size of 4.3, 14.7, and 27.6 μm, where the lactose crystals were extracted for 24, 72, and 120 h from the lactose solution, respectively. The inset shows deconvoluted spectra of the absorptions due to Lα·H2O crystal, in which the absorption of anhydrous Lα was removed. It was mentioned that the broad-band was numerically removed by a polynomial cubic function, and the absorption coefficient is normalized by the weight of the lactose sample because the density of lactose was different in aperture for each sample. As expected from the particle shape, lactose was mainly crystallized into Lα·H2O because the spectra were dominated by the absorptions due to Lα·H2O (α1 and α2). Anhydrous Lβ was also crystallized. It was probably formed on the Lα·H2O crystal in the drying process, but the content in the extracted lactose was within 4%. The absorption features of α1 and α2 were dependent on the particle size as shown in the inset of Fig. 6. Details on the integrated absorption intensity and the full-width at half-maximum (FWHM) for the average particle size are shown in Fig. 7, in which the results on the average particle size of 3 μm is for the preliminarily used Lα·H2O powder. The α1 line width (filled-square) was independent of the particle size. A significantly narrow absorption spectrum of α1 with a FWHM of about 1.3cm1 has been reported by using a photomixer and continuous-wave, however, the obtained THz-TDS was broader than that,8 which is because of the resolution of the spectroscopy methods. Since the FWHM of 3.7cm1 obtained in this work was smaller than that of 4.91cm1 as demonstrated by THz-TDS,8 but much broader compared to the native feature with an FWHM less than 1.3cm1, it is considered that the spectrum width was not apparently changed in the results under the resolution limit of the THz-TDS system. Previously, the significantly sharp feature of α1 was discussed by the solid-state density functional theory and assigned to externally hindered rotational modes in the crystal b-axis, not to the internal modes.14 In such modes, the absorption intensity can be expected to increase with the crystal size much less than the wavelength of the incident THz-electromagnetic wave as shown in Fig. 7 (closed circles). On the other hand, the absorption of α2 was broader than that of α1 and the line width increased with the particle size (open-squares). However, the integrated intensity of α2 (filled-squares) was scarcely increased with the particle size, but the dependence was very small compared to α1. The vibration mode of α2 was identified as molecular-rotation in a lactose active mode by the first principles calculation.10 The significant difference of α2 feature is possibly due to the absorption is due to the internal vibrations associated with intermolecular hydrogen-bonding. The significant difference of α1 and α2 intensities resulted in assignment of the particle size. Figure 8 shows the intensity ratio of α2/α1 for the average particle size. A ratio as high as 3.98 was linearly decreased with the particles’ size and to 3.45 for the size of 27.6 μm, where the correlative square-factor was 99.2% on the least squares method. As a result, if the particles are contained inside samples, the average particle size can be nondestructively evaluated by the ratio of α2/α1 obtained by THz-TDS.

Fig. 6

Absorption spectra of extracted lactose from the supersaturated solution with the average size of (a) 4.3, (b) 14.7, and (c) 27.6 μm. The inset shows the deconvoluted α-lactose monohydrate spectra.


Fig. 7

Variation of relative integrated absorption intensity of α1 and α2 (filled circles and open circles) and full-width at half-maximum of α1 and α2 (filled squares and open squares) for average crystallized particle size of α-lactose monohydrate.


Fig. 8

Dependence of integrated absorption ratio of α2/α1 on average crystallized particle size of α-lactose monohydrate.




THz-TDS was applied to evaluate the average particle size of crystallized α-lactose monohydrate. The dependence of THz-absorption on the particle size below 30 μm was different for the two spectra, with peaks at 17.5cm1 (0.53 THz) and 45.6cm1 (1.37 THz). The integrated absorption intensity ratio, which is linearly dependent on the particle size, was able to assign the average particle size. The results in this work showed a THz-absorption feature in organic crystal that incorporates hydrogen-bonding was dependent on the crystal size which was far smaller compared to the wavelength of the incident THz-wave and can present a new concept focused on the absorption intensities to evaluate the crystal size of organic crystal.


This research was partially supported by JSPS KAKENHI Grant-in-Aid for Exploratory Research Number 24655165.



F. Miyamaruet al., “THz-time-domain spectroscopy of amino acids in solid phase,” in Abst. Conf. Lasers and Electro-Optics, OSA, CMG3 (2003).Google Scholar


M. Waltheret al., “Far-infrared vibrational spectra of all-trans, 9-cis and 13-cis retinal measured by THz time-domain spectroscopy,” Chem. Phys. Lett. 332(3–4), 389–395 (2000).CHPLBC0009-2614http://dx.doi.org/10.1016/S0009-2614(00)01271-9Google Scholar


M. Yamaguchiet al., “Terahertz absorption spectra of L-, D-, and DL-alanine and their application to determination of enantiometric composition,” Appl. Phys. Lett. 86(5), 053903 (2005).APPLAB0003-6951http://dx.doi.org/10.1063/1.1857080Google Scholar


S. Yamauchiet al., “Determination of stereoisomer in sugars by THz-TDS,” Proc. SPIE 8011, 80118L (2011).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.902044Google Scholar


B. B. HuM. C. Nuss, “Imaging with terahertz waves,” Opt. Lett. 20(16), 1716–1718 (1995).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.20.001716Google Scholar


K. Kawaseet al., “Non-destructive terahertz imaging of illicit drugs using spectral fingerprints,” Opt. Exp. 11(20), 2549–2554 (2003).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.11.002549Google Scholar


M. Yamashitaet al, “Backside observation of large-scale integrated circuits with multilayered interconnections using laser terahertz emission microscope,” Appl. Phys. Lett. 94(19), 191104 (2009).APPLAB0003-6951http://dx.doi.org/10.1063/1.3133346Google Scholar


E. R. Brownet al., “On the strong and narrow absorption signature in lactose at 0.53 THz,” Appl. Phys. Lett. 90(6), 061908 (2007).APPLAB0003-6951http://dx.doi.org/10.1063/1.2437107Google Scholar


A. Roggenbucket al., “Coherent broadband continuous-wave terahertz spectroscopy on solid-state samples,” New J. Phys. 12(4), 043017 (2010).NJOPFM1367-2630http://dx.doi.org/10.1088/1367-2630/12/4/043017Google Scholar


S. Saitoet al., “First principles calculation of terahertz vibrational modes of a disaccharide monohydrate crystal of lactose,” Jpn. J. Appl. Phys. 45(43), L1156–L1158 (2006).JJAPA50021-4922http://dx.doi.org/10.1143/JJAP.45.L1156Google Scholar


A. Van Krevald, “Growth rates of lactose crystals in solutions of stable anhydrous α-lactose,” Neth. Milk Dairy J. 23(4), 258–275 (1969).Google Scholar


P. WalstraR. Jennes, Dairy Chemistry and Physics, p. 35, Wiley, New York (1984).Google Scholar


T. D. Dinceret al., “Crystallisation of a-lactose monohydrate from dimethyl sulfoxide (DMSO) solutions: infuence of β-lactose,” J. Cryst. Growth 205(3), 368–374 (1999).JCRGAE0022-0248http://dx.doi.org/10.1016/S0022-0248(99)00238-9Google Scholar


D. G. Alliset al., “Assignment of the lowest-lying THz absorption signatures in biotin and lactose monohydrate by solid-state density functional theory,” Chem. Phys. Lett. 440(4–6), 203–209 (2007).CHPLBC0009-2614http://dx.doi.org/10.1016/j.cplett.2007.04.032Google Scholar



Satoshi Yamauchi received the BE, ME, and Dr. Eng. degrees in 1986, 1988, and 1991, respectively, all in electronic engineering, Tohoku University, Japan. In 1991 he joined the ULSI R&D Center, OKI Electric Industry, Japan. He was promoted to lecturer in the Department of Systems Engineering, Ibaraki University, in 1997 and now is an associate professor, Department of Biomolecular Functional Engineering, Ibaraki University. His current interests are electronic materials and electronic devices. He is a member of the Institute of Electrical and Electronics Engineers (IEEE) and the Japan Society of Applied Physics.


Sakura Hatakeyama was a student in Department of Biomolecular Functional Engineering, Ibaraki University, Japan, from 2009 to 2013. Her major in research was evaluation of biomolecules by terahertz spectroscopy. She received the BE degree in 2013 from Ibaraki University, Japan.


Yoh Imai received the BS, MS, and Dr. Eng. degrees from Hokkaido University, Sapporo, Japan, in 1977, 1979, and 1982, respectively, all in applied physics. He was a research associate in the Department of Engineering Science, Hokkaido University from 1982 to 1989. From 1987 to 1989 he was a visiting scholar at the University of Toronto, Canada. He was an associate professor in the Department of Computer Science and Electronics, Kyushu Institute of Technology, Japan, from 1989 to 2002. In 2002, he was promoted to a professor in the Department of Electrical and Electronic Engineering, Ibaraki University, Japan. His current interests are the terahertz wave technology and its applications on terahertz sensing and the optical nonlinearity in optical fiber and its applications on optical sensing. He is a member of the Optical Society of America (OSA), the Optical Society of Japan, and the Japan Society of Applied Physics.


Masayoshi Tonouchi received the BS, MS, and Dr. Eng. degrees form Osaka University, Japan, in 1983, 1985, and 1988, respectively. From 1988 to 1989 he worked at the Faculty of Engineering Science of the same university. From 1989 to 1994 he joined the Faculty of Computer Science and System Engineering, Kyushu Institute of Technology, Japan. In 1994, he moved to Kansai Advanced Research Center, Communications Research Laboratory, Japan. From 1996 to 2000, he was an associate professor at the Research Center for Superconducting Materials and Electronics, Osaka University. Currently he is a professor in the Institute of Laser Engineering, Osaka University and a concurrent professor of Nanjing University. His current research interests include ultrafast optical and terahertz science of strongly correlated electron systems, optical interfaces for single-flux-quantum circuits, and development-and-applications of terahertz systems such as the laser terahertz emission microscope. He is a member of the Optical Society of America, the Japan Society of Applied Physics, the Physical Society of Japan, and the Institute of Electronics, Information and Communication Engineers.

Satoshi Yamauchi, Sakura Hatakeyama, Yoh Imai, Masayoshi Tonouchi, "Nondestructive evaluation of crystallized-particle size in lactose-powder by terahertz time-domain spectroscopy," Optical Engineering 53(3), 031203 (2 December 2013). http://dx.doi.org/10.1117/1.OE.53.3.031203




Terahertz spectroscopy

Terahertz radiation


Nondestructive evaluation

Back to Top