Open Access
16 May 2017 Sensing of Streptococcus mutans by microscopic imaging ellipsometry
Mai Ibrahim Khaleel, Yu-Da Chen, Ching-Hang P. Chien, Yia-Chung Chang
Author Affiliations +
Abstract
Microscopic imaging ellipsometry is an optical technique that uses an objective and sensing procedure to measure the ellipsometric parameters Ψ and Δ in the form of microscopic maps. This technique is well known for being noninvasive and label-free. Therefore, it can be used to detect and characterize biological species without any impact. Microscopic imaging ellipsometry was used to measure the optical response of dried Streptococcus mutans cells on a glass substrate. The ellipsometric Ψ and Δ maps were obtained with the Optrel Multiskop system for specular reflection in the visible range (λ=450 to 750 nm). The Ψ and Δ images at 500, 600, and 700 nm were analyzed using three different theoretical models with single-bounce, two-bounce, and multibounce light paths to obtain the optical constants and height distribution. The obtained images of the optical constants show different aspects when comparing the single-bounce analysis with the two-bounce or multibounce analysis in detecting S. mutans samples. Furthermore, the height distributions estimated by two-bounce and multibounce analyses of S. mutans samples were in agreement with the thickness values measured by AFM, which implies that the two-bounce and multibounce analyses can provide information complementary to that obtained by a single-bounce light path.

1.

Introduction

Spectroscopic ellipsometry (SE) is a popular metrology technique used to determine thin film thickness and its optical constants. The technique allows us to reconstruct a nanoscale resolution image of the material surface by fitting the measured ellipsometric parameters psi (Ψ) and delta (Δ) with a suitable theoretical simulation.1 Here Ψ is a measure of the absolute value of the ratio of the complex Fresnel reflection coefficients (or reflectivities), rp and rs, for p- and s-polarized lights, and Δ is the phase difference between rp and rs. Combing SE with microscopy leads to a powerful noncontact imaging technique known as microscopic imaging ellipsometry (MIE).2 MIE enables a large area imaging with a spatial resolution of 1  μm and subnanometer-scale z-resolution.2,3 In addition to the subnanoscale thickness resolution, it does not require a vacuum ambient nor the labeling of molecules, as compared to other detection techniques that require labeling—usually fluorescent or enzymatic such as fluorescence imaging.46 The labeling itself in such techniques can be highly heterogeneous, leading to nonquantitative results.6 MIE enables faster real-time measurements (in seconds) in comparison with other slow scanning techniques such as AFM that requires a few minutes to scan one area. This is especially useful for biological processes that occur within seconds rather than minutes.7,8 These unique features show the versatility of this imaging technique which had increased its usage in sensing and analyzing biological systems.

One of the pioneering applications of imaging ellipsometry was to detect antigen–antibody interactions using a single wavelength light source and fitting the measured data before and after interaction happens.9 Other applications of imaging ellipsometry were reported, including studies of biological reaction in cancer cells,10 detection of insulin–antibody binding on solid substrate,11 detection of latent fingermarks,12 as well as mapping of multilamellar films, hydrated membranes, and fluid domains.13

Streptococcus mutans is an aerobic, gram-positive bacterium with coccus shape found in the human oral cavity.14,15 It has been reported as the main etiological agent of dental caries and normal inhabitant of dental plaque.16 S. mutans is known for its ability to synthesize extracellular polymeric substance (EPS), which works as a matrix holding microbial cells together to form a biofilm that attaches to teeth surface. This matrix serves as the colony of bacteria and contributes to structure maintenance and antibiotic resistance of the biofilm.17,18 In the metabolism process, S. mutans produces aggressive substances such as acids, which cause a low pH environment in the mouth and lead to corrosion of enamel and the formation of cavities.18,19

In recent years, several research groups have developed methods for detecting and studying S. mutans in either laboratory settings or field trials. To our knowledge, there have been no attempts to detect S. mutans cells in the visible region using MIE. In this work, we performed MIE measurements of S. mutans cells grown on a glass substrate and dried in air. We assumed a structure consisting of air/bacterial cells/glass substrate to study the dielectric properties and height distribution. The measured Ψ and Δ maps were fitted using three different optical models to detect bacterial cells and the surrounding dried culturing medium. In the first model, we considered only a single-bounce light path with the reflection from the interface between air and bacterial cells. In the second model, we considered the coherent superposition of the reflection from the interface between glass substrate and sample and the first reflection and called this the “two-bounce” model. In the third model, we included all reflections in the structure and called it the “multibounce” model. AFM Veeco Innova20 tapping mode was used in this study to map the surface topography of S. mutans cells and culturing medium. The observed height distributions of cells and culturing medium on glass substrate by AFM were compared to the deduced values of height distribution by two-bounce and multibounce analyses. The glass substrate was chosen to grow S. mutans for two main reasons. The first one is that a series of studies have been done by using a glass substrate to study and reveal many details about S. mutans adherence and other oral bacteria.21 The second reason is that glass substrate is easy to calibrate, which makes the analyses simpler.

2.

Materials and Methods

2.1.

Cell Cultures

An S. mutans strain (ATCC 25175) was purchased from Creative Life Sciences (CMP), Taipei, Taiwan.22 S. mutans was incubated in brain-heart infusion broth (BHI) purchased from CMP, Taipei, Taiwan,22 shaken at 37°C overnight. Another BHI solution was prepared to grow bacterial cells on a glass substrate. The solution consists of 7.4-g BHI, 20-g sucrose, and 200-ml Millipore water.23 Sucrose was added in order to enhance the rate of cell growth on a glass substrate.

2.2.

Cell Sample Preparation

Glass slides (BRAND GMBH + CO KG, Wertheim, Germany)24 were cut into 1×1  cm2 pieces. They were cleaned through sonication in acetone for 30 min, then for 20 min in isopropanol solution at room temperature. After that glass pieces were washed and sonicated in Millipore water for 10 min and finally dried in an oven at 60°C. S. mutans was cultivated on glass substrates as described by Hu et al.23 The cleaned glass substrates were immersed into a Petri dish containing 10 ml of freshly prepared BHI medium with sucrose. Then 2.5 ml of incubated S. mutans cells in BHI medium overnight was added to the Petri dish contacting the immersed substrates, sealed by parafilm and incubated for 72 h at 37°C. The glass samples were then washed three times with Millipore water and placed in a new Petri dish to dehydrate naturally for 1 day, then measured by MIE and AFM. Bare glass substrates kept in culturing medium (BHI and sucrose) without S. mutans inoculated were used in reference experiments.

2.3.

Multiskop Imaging Ellipsometry

The ellipsometric parameters Ψ and Δ of S. mutans cells on glass substrate were measured using an Optrel Multiskop system,25 as shown in Fig. 1 which is modified to a rotating compensator ellipsometry. Converting the system from null measurements (which requires rotating both the polarizer and analyzer) to RCE allowed us to avoid the image-shift problem caused by the rotating analyzer (a Glan–Taylor prism26) and to a faster extraction of the ellipsometric image.1

Fig. 1

Schematic diagram for rotating compensator ellipsometry.

JBO_22_5_056005_f001.png

The setup of our MIE consists of a laser arm, sample stage, and detector arm. The light beam comes from a Fianium WhiteLase short wavelength supercontinuum laser that covers wavelengths from 400 to 1000 nm with power >300  mW,27 then passes through a polarizer fixed at 45 deg and a rotating compensator. To enhance the intensity of the incident light, a 10× lens was inserted before the light hits the spot of interest on the sample. The reflected or scattered light from the sample then goes through an 80× objective lens28 that makes our system capable of probing an area as small as 1  μm×1  μm. The light beam then passes through an analyzer mounted at 90 deg with respect to the polarizer, and finally is collected by a charge-coupled device (CCD) camera (Pixelfly qe),29 with a pixel size equal to 164 nm. Measurements of ellipsometric parameters in this study were done using specular mode. In this mode, both the laser arm and detector arm are rotated in opposite directions with the use of a goniometer. The data collected by the CCD are the reflected light beam from the sample which is analyzed to determine the ellipsometric parameters Ψ and Δ. In RCE, the reflected light intensity as a function of compensator angle ϕc is given as

Eq. (1)

I(ϕc)=A0+A2cos(2ϕc)+B2sin(2ϕc)+A4cos(4ϕc)+B4sin(4ϕc).

We take measurements of I(ϕc) for 36 values of ϕc starting from 5 deg and ending at 355 deg with 10 deg spacing. These data allow us to determine the Fourier coefficients in Eq. (1), which can be used to evaluate Ψ and Δ.

2.4.

Theory and Numerical Calculations

The phase difference between rp and rs (where rp and rs are reflectivity for p- and s-polarized light beams) and the ratio of their amplitudes can be quantified by the ellipsometric measurements as Δ and tanΨ, respectively. The analysis of RCE signals was done to determine Ψ and Δ in terms of Fourier coefficients of the output light intensity as a function of the compensator angle.10

The fundamental equation of ellipsometry is then written as

Eq. (2)

ρ=tanΨeiΔ=  rprs  ,
where rp and rs are complex reflectivities for p and s polarizations. For the sample under consideration, the reflectivities can be calculated from Fig. 2 using three different models.

Fig. 2

Schematic diagram for light scattering paths from the sample under consideration: (a) single-bounce detection, (b) two-bounce detection, and (c) multibounce detection.

JBO_22_5_056005_f002.png

The light detected by an image sensor is composed of light paths after a single bounce, two bounces, or multiple bounces.30 For the single-bounce case, the coherent superposition of directly reflected light after a single reflection event with the sample is detected, while in the two- and multibounce cases the light is allowed to hit the interface between the sample and the substrate before reaching the sensor.31

In traditional methods of shape and reflection estimations such as time of flight and photonic mixer devices for opaque objects, the single-bounce measurement is commonly used, while the multibounce effect is ignored since it is considered as a source of noise.31 Recent studies of shape and reflection estimations have shown that two-bounce and multibounce light paths can provide information not achievable by using only single-bounce analysis. The two-bounce light paths were used for recovery of shape and reflection estimations,31 while multibounce light paths have been used to characterize flat and transparent thin films using a spectroscopic rotating polarizer-analyzer ellipsometer.32 To fully characterize the sample under consideration in this paper, we incorporated analyses of single-bounce, two-bounce, and multibounce models for the data obtained using RCE. The derivation of formulas relevant to the three models to be used is given in Appendix A.

2.4.1.

Single-bounce light path

The single-bounce light path was used to detect bumpy areas, including S. mutans cells and the culturing medium on glass substrate. Since each grid (with 5 pixel average) of bumpy area has a curved surface and size of 820 nm which is comparable to the wavelength, it is assumed that light going through the two-bounce and multibounce paths can be neglected [see Fig. 2(a)]. Therefore, light scattered from bumpy areas may be better described by the coherent sum of single-bounce light rays. By calculating rp and rs, we can determine the refractive index n1 and the extinction coefficient k1 of bumpy areas according to the following equation:

Eq. (3)

ϵ1=(n1+ik1)2=(1ρ)2(1+ρ)2sinθ02tanθ02+sinθ02,
where ϵ1 is the complex refractive index of bumpy area.

2.4.2.

Two-bounce light path

The model with a two-bounce light path is used to detect S. mutans cells and the dried culturing medium that represents a mixture of rough and nearly flat areas on the glass substrate. For the two-bounce light path shown in Fig. 2(b), the reflectivity is given as

Eq. (4)

rj=r01j+(1  r01j2)r12j  e2iβ,
where j stands for s-polarization and p-polarization, and β represents the phase factor for light going through the film given as

Eq. (5)

β=2πhn1cosθ1λ,
where h and n1 are the thickness and refractive index of the sample, respectively. Assuming the sample does not absorb light in the visible region (i.e., k1=0), then we have

Eq. (6)

ln|γ|=ln|e2iβ|=ln|r01pρr01s  ρ(1r01s2  )r12s  (1r01p2)r12p|=0.
An in-house MATLAB® program33 was used to solve Eq. (6) and find the refractive index and thickness of the cells and medium.

2.4.3.

Multibounce light paths

The model with multibounce light paths is also used to analyze the data. Based on the results of Figs. 4(a) and 4(b), we expect this model will give similar results as a two-bounce model (for our case with AOI=52  deg and height=80 and 350 nm), since the higher order terms in the multibounce model appear to be small. The reflectivity of multibounce model is given as

Eq. (7)

rj=r01j+r12j  e2iβ1+r01jr12j  e2iβ,
where r01j and r12j for s-polarization and p-polarization have been derived in Ref. 32. Assuming that k1=0 in the visible range, the final solution is written as

Eq. (8)

ln|γ|=0,
where γ is given in Eq. (20) of Appendix A. Our in-house MATLAB® program was also used to solve Eq. (8) and find the refractive index n1 of the sample. Given n1, the thickness of the sample can be obtained from Eq. (5) in terms of β, n1, and cosθ1 as

Eq. (9)

h=βλ2πn1cosθ1.

3.

Results and Discussion

For MIE measurements, three samples of S. mutans were investigated, while only the results for two samples are reported here, since the third sample gives qualitatively the same results. A carbon tape was attached to the backside of the glass substrate of these samples in order to eliminate backside reflections of the glass substrate in our measurements. For AFM measurements, two samples were investigated using tapping mode to map the surface morphology of cells and culturing medium. All the samples investigated by AFM and MIE were prepared under the same conditions.

3.1.

AFM Results

AFM tapping mode was used to morphologically characterize S. mutans cells cultured on a glass substrate for 72 h. All images were obtained with a commercial AFM (Veeco Innova) of ambient conditions. The OMCLO-AC160TS, R3 tip was used for tapping mode images.34 Figure 3 shows AFM images of S. mutans cells and heights distribution on a glass substrate with spatial resolution of 78  nm/pixel.

Fig. 3

S. mutans on a glass substrate after 72-h culture: (a) 20  μm×20  μm AFM topography image, (b) height profile of culturing medium (red, labeled 1) and S. mutans cells (green, labeled 2), and (c) height distribution of culturing medium and S. mutans cells.

JBO_22_5_056005_f003.png

The image shown in Fig. 3(a) shows the topography of S. mutans cells with different heights, while Fig. 3(b) shows the height profiles of S. mutans cells in red and culturing medium in green. Figure 3(c) shows the height distribution of S. mutans and culturing medium on glass substrate. The figure revealed that S. mutans cells are of height from 250 to 500 nm with maximum count at 350 nm, and a few stacked-up cells have heights of 600 to 650 nm. These heights correspond to the bright spots in Fig. 3(a). The height distribution of dried culturing medium (BHI and sucrose) on glass substrates in Fig. 3(c) shows height ranging from 10 to 200 nm with maximum count at 80 nm.

3.2.

MIE-RCE Results

Before performing ellipsometry measurements of S. mutans cells on glass substrate, the Ψ and Δ spectra were measured for a bare glass substrate with a carbon tape attached to the backside in order to block the backside reflections of the glass substrate. Cauchy model35 was used to fit Ψ and Δ spectra in the visible range (450 to 750 nm) and obtain the refractive index of glass substrate (1.549, 1.538, 1.532 at 500, 600, and 700 nm, respectively). MIE measurements of S. mutans samples cultured on a glass substrate for 72 h were done using specular mode at an angle of incidence (AOI) of 52 deg in the visible range (450 to 750 nm) and the analyses were done for three specific wavelengths. In order to overcome the signal noise in ellipsometry measurements, the values of Ψ and Δ of 5×5  pixels are averaged to produce a mean value which is registered for the center pixel. We shift the center pixel by 2 pixels (in both x and y directions) and then take the average of the 5×5  pixels around it, and then move on. This leads to a smoothed image with a total area of 58  grids×58  grids, with each grid being 328  nm×328  nm in size (which contains signals average over the surrounding 820  nm×820  nm area). The averaging procedure was done by using an in-house MATLAB® program to produce smoothed images of the sample.

To show the difference among single-bounce, two-bounce, and multibounce analyses, tanΨ and cosΔ obtained for the three models are plotted as functions of n1 in the interval (1, 3) as shown in Fig. 4. Here, tanΨ=|rp/rs| and cosΔ=Re(rp/rs)/tanΨ. In Fig. 4, tanΨ and cosΔ are plotted for three different wavelengths 500, 600, and 700 nm, except for the single-bounce case in Fig. 4(a) since the result is independent of wavelength. The height of the sample considered in Figs. 4(a) and 4(b) was set to 80 and 350 nm, where 80 nm is the average height of dried culturing medium and 350 nm is the average height of S. mutans revealed by AFM measurements in Fig. 3(c).

Fig. 4

tanΨ and cosΔ as functions of n1 at AOI of 52 deg for (a) two-bounce model (symbols) and multibounce model (continuous curves) at height=80  nm, and (b) two-bounce model (symbols) and multibounce model (continuous curves) at height=350  nm. Single-bounce results are also included in (a) to compare with two- and multibounce models.

JBO_22_5_056005_f004.png

Figures 4(a) and 4(b) show that results obtained by two-bounce and multibounce analyses overlap well for two different heights. This means that results of samples analyzed by two-bounce and multibounce models will be very close. As a calibration, in Fig. 5 we show the results of Ψ and Δ images for culturing medium sample (without S. mutans) at a wavelength equal to 600 nm analyzed using single-bounce, two-bounce, and multibounce models.

Fig. 5

CCD, Ψ, Δ, n1, k1, and h images of BHI sample using single-bounce, two-bounce, and multibounce models at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f005.png

We noticed that the topology of the n1 image is very similar to the Ψ image, while the k1 image is very similar to the Δ image. Comparing n1 images with the corresponding Ψ images in Fig. 5, we conclude that the cyan color areas in n1 images with n1 values ranging from 1.25 (1.39) to 1.3 (1.41) in single-bounce (two-bounce and multibounce) analyses coincide with the areas in dark blue of Ψ images. The yellow areas of n1 images match the light blue areas in Ψ images, while the dark blue spots in n1 images match the red spots in Ψ images. In single-bounce analysis, the yellow-color areas in the n1 image correspond to the blue areas in the k1 image which has negative (unphysical) values. This could be caused by depolarization of scattered light or it is an indication that those areas are flat and not suitable for single-bounce analysis. Even with unphysical k1 values in some areas, the n1 values obtained by the single-bounce model are still fairly close to the results obtained by two-bounce and multibounce methods. The blue spots in the n1 image correspond to the yellow and red spots in the k1 image which has positive values, indicating possible absorption there.

The n1 images obtained from two-bounce and multibounce analyses are very close (with difference less than 0.02 in all pixels). In the multibounce analysis, there are two roots for n1 corresponding to the “+” and “−” signs in Eq. (20). It is interesting to note that the physical results for n1 of the culturing medium are found only by adopting the “−” sign in Eq. (20) for this sample.

Comparing the k1 image in single-bounce analysis with h images in two-bounce and multibounce analyses, we found that the blue stripes in the k1 image match the red stripes in h images with height between 200 and 280 nm. Other areas correspond to k1 near zero and height less than 40 nm, which are non-absobing based on single-bounce analysis. It is interesting to note that the height of the culturing medium revealed by the two-bounce and multibounce analyses is either below 40 nm (flat areas) or above 200 nm (bumpy areas) and the topography of the h image is very close to the k1 image and Δ image.

3.2.1.

Single-bounce analysis

The single-bounce model was first used to analyze S. mutans with culturing medium on glass substrate. The left-side images of Figs. 6 and 7 show the CCD images of two measured areas (area 1 and area 2) of interest for S. mutans samples. The measured Ψ and Δ images at wavelengths of 500, 600, and 700 nm for the two different areas are shown in the second and third columns of Figs. 6 and 7.

Fig. 6

CCD, Ψ, Δ, refractive index n1, and extinction coefficient k1 images of S. mutans sample for area 1 using single-bounce analysis at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f006.png

Fig. 7

CCD, Ψ, Δ, refractive index n1, and extinction coefficient k1 images of S. mutans sample for area 2 using single-bounce analysis at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f007.png

Ψ and Δ images at different wavelengths were used to calculate ρ in Eq. (2). The obtained value of ρ is substituted in Eq. (3), and our in-house MATLAB® program was used to solve this equation and find values of refractive index n1 and extinction coefficient k1 for the three wavelengths and show their images in the fourth and fifth columns of Figs. 6 and 7.

We noticed again that the topology of the n1 image is very similar to the Ψ image, while the k1 image is very similar to the Δ image. Comparing n1 images with the corresponding Ψ images in Figs. 6 and 7 for the three wavelengths, we conclude that areas in blue color of n1 images with n1 ranging from 1.2 to 1.4 with a few (dark blue) spots less than 1.2 coincide with the blue areas in Ψ images.

The cyan area and green stripes (green and yellow stripes) in Fig. 6 (Fig. 7) for k1 images have positive values between 0 and 0.2 which are physically reasonable, and they correspond to blue areas in Δ images with values between 180  deg and 160  deg. The corresponding n1 values for these areas are between 1.6 and 2 which correspond to bumpy areas in culturing medium. There are a few spots in k1 images with values larger than 0.2 which indicates strong light trapping of the bumpy structure. The areas in k1 images with negative values (blue to cyan) are attributed to strong depolarization effect and the corresponding values of n1 there are between 1.2 and 1.6.

It should be noted that the depolarization effect due to scattering from S. mutans cells and rough areas of culturing medium has been neglected in this model. Thus, the n1 and k1 values obtained can be rather crude. The k1 value obtained can be misleading, since it is very sensitive to the phase and the depolarization effect. On the other hand, this one-bounce analysis is quite effective in locating which area is bumpy and the distribution of S. mutans cells can be related to the n1 distribution obtained.

3.2.2.

Two-bounce analysis

For the two-bounce analysis, the calculated ρ from measured Ψ and Δ images were substituted in Eq. (6). This equation was solved for different wavelengths to find the refractive index. The third column of Figs. 8 and 9 show the n1 images for the two samples obtained by solving Eq. (6) for wavelengths of 500, 600, and 700 nm.

Fig. 8

CCD, Ψ, Δ, refractive index n1, and height distribution h images of S. mutans sample for area 1 using two-bounce analysis at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f008.png

Fig. 9

CCD, Ψ, Δ, refractive index n1, and height distribution h images of S. mutans sample for area 2 using two-bounce analysis at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f009.png

Areas in cyan for n1 images in Figs. 8 and 9 show n1 values ranging from 1.2 to 1.4, and areas with light green and some grids with yellow, orange, and red colors show values ranging from 1.4 to 2.6. Comparing results of n1 images obtained from two-bounce and single-bounce analyses, we found that areas in cyan in Figs. 8 and 9 correspond to the same areas in blue color in Figs. 6 and 7, and n1 values obtained by both models are similar, although the k1 values there are negative (unphysical). This indicates that bacteria cells surrounded by culturing medium represented by the blue areas in n1 images in Figs. 6 and 7 can still be detected by single-bounce analysis, and the unphysical values of k1 there are due to depolarization effect. This is confirmed by the two-bounce analysis.

The two-bounce model can be used to analyze and sense the S. mutans cells in areas where the height distribution of cells is uniform, similar to height profile at 350 nm as shown in AFM results (green line) of Fig. 3(b). In addition, the refractive index values obtained from the two-bounce analysis can be used to calculate the height distribution of S. mutans cells surrounded by culturing medium according to Eq. (9). Thus, the two-bounce analysis is more versatile for the current application. The fifth columns of Figs. 8 and 9 show the obtained height (h) distribution for area 1 and area 2. The blue and cyan areas with h ranging from 40 to 200 nm in Figs. 8 and 9 correspond to areas filled with culturing medium. Areas with h above 200 nm (green, yellow, orange, and red) correspond to areas filled with cells, in agreement with the height distribution of S. mutans cells and dried culturing medium obtained by AFM in Fig. 3(c).

3.2.3.

Multibounce analysis

The first two columns of Fig. 10 show images of n1 and h distributions in area 1 for λ=500  nm obtained by solving Eq. (8) with both “−” and “+” signs in Eq. (20), respectively. In both calculations, the dark blue area represents the unphysical zone (since n1 obtained there is less than 1). We found that the two solutions are mutually exclusive with the unphysical zone obtained with “−” sign being the physical zone obtained with “+” sign, and vice versa. Furthermore, the physical values of n1 range from 1.3 to 1.5 (red zone) in the first column and from 1.5 to 3 (mixed-color zone) in the second column. The corresponding heights in the first column ranges from 200 to 300 nm, while in the second column it ranges from 20 to 180 nm. So, the two solutions of the multibounce analysis can be used to distinguish regions of low index (high h value) and high index (low h value), which correspond to regions filled with high and low concentrations of S. mutans cells, respectively. We can combine the physical solutions for n1 obtained with both signs, and the outcome is presented in the third column. The result is very similar to that obtained by the two-bounce analysis as shown in the first row of Fig. 8.

Fig. 10

Multibounce analysis for refractive index n1 obtained by using “+” and “−” signs in Eq. (20) and height (h) distribution images of S. mutans sample for area 1 at λ=500  nm at AOI=52  deg. The frame size is 19  μm×19  μm. The units of Ψ and Δ are degrees.

JBO_22_5_056005_f010.png

In the above analysis, we have neglected the effects of depolarization and stray light, which can lead to errors in our estimate of n1 and h. Thus, the n1 and height distributions deduced can be different for different wavelengths as shown in Figs. 6 and 7 as well as in Figs. 8 and 9. The height distribution analysis for different wavelengths of S. mutans and dried culturing medium in our study should be wavelength independent and the variation of n1 in the range of wavelengths considered should also be negligible. Thus, this variation in results for different wavelengths gives an estimate of the error caused by depolarization effect. The error estimate of n1 values obtained by single bounce and two bounce (same as multibounce) analyses are reported in Table 1, which shows a standard deviation in the 1% range. For the h distribution, since h varies significantly from pixel to pixel in the samples considered and the image resolution depends on wavelength, it is not straight-forward to estimate the error in h distribution by comparing results for different wavelengths. However, we can see that in areas where the height distribution is uniform, the results obtained become rather insensitive to the wavelength, while in bumpy areas the variation can be significant. The overall ranges and the topologies of h in Figs. 8 and 9 for different wavelengths are also similar. However, changes of distributions for wavelength from 500 to 700 nm are clearly visible.

Table 1

Standard deviation and mean value of n1 obtained by single-bounce and two-bounce models for area 1 and area 2 at AOI=52  deg.

Single bounceTwo bounce
n1n1
σ¯n1¯σ‾n1¯
Area10.1441.5760.1691.557
Area20.131.590.1611.567

Finally, we give some estimate of the statistical errors in our measurements. We determine the average standard deviation (σ¯) by repeating the same measurements five times for MIE measurements of a bare glass substrate with carbon tape attached to the backside. The mean value of Ψ and Δ maps (Ψ¯,Δ¯) and average standard deviation (σ¯) of this study are shown in Table 2. These statistical values were calculated for single-pixel maps and maps obtained by averaging over five surrounding pixels, since results reported in this study were based on Ψ and Δ maps averaged over five surrounding pixels explained in Sec. 3.2. As shown in Table 2, the average standard deviation (σ¯) for the 5-pixel average map is about one half of the single-pixel map, thus the resulting MIE image will be less noisy.

Table 2

Standard deviation and mean value of Ψ and Δ measured by MIE for glass substrate, area 1, and area 2 at AOI=52  deg.

Glass substrateArea1Area2
ΨΔΨΔΨΔ
σ‾Ψ‾σ‾|Δ|¯σ‾Ψ‾σ‾|Δ|‾σ‾Ψ‾σ‾|Δ|‾
500 nmSingle pixel0.15514.620.653178.20.339.8252.314167.20.39.1692.07175.9
5-pixel average0.11114.710.524178.90.13110.320.993178.10.1138.1020.726173.1
600 nmSingle pixel0.25114.320.563178.80.34610.772.1168.90.3378.32.03168.1
5-pixel average0.1514.360.468179.10.10610.330.708174.50.1078.6590.646177
700 nmSingle pixel0.26914.770.64178.80.4769.9512.533177.60.4427.4812.377178.3
5-pixel average0.12714.450.5261790.15710.340.7651730.13410.420.78174

We can also estimate the error introduced in our Ψ and Δ values determined by using a finite number of compensator angles (ϕc) in our measurements. If there is no statistical error or systematic error, one should be able to determine the five Fourier coefficients in Eq. (1) with just five values of ϕc. Given the measurements at 36 values of ϕc, we can select them all (with 10 deg spacing), or adopt only 18 values at 20 deg spacing, or 12 values at 30 deg spacing. We then have three ways to determine the Ψ and Δ values. The deviation of these determined values from the mean gives us an estimate of the measurement error. The estimated average standard deviation σ¯ for areas 1 and 2 of our sample with S. mutans cells are also listed in Table 2.

For single-pixel Ψ maps, error ranges from 1.8% for glass substrate to 4.7% for area 1 and 5.9% for area 2. In single-pixel Δ maps, error ranges from 0.3% for glass substrate to 1.4% for area 1 and 1.3% for area 2, indicating that 5-pixel averaged maps are more accurate and reproducible by MIE.

4.

Conclusions

In this paper, we performed MIE measurements in the visible range to detect dried S. mutans cells and surrounding culturing medium on glass substrate. We found that two-bounce and multibounce lead to very close results for refractive index values and height distributions. Incorporating single-bounce, two-bounce, and multibounce analyses can reveal more information about the biological sample under investigation. The two-bounce and multibounce models provide the height distribution, while the single-bounce model provides additional information on the extinction coefficient or the depolarization effect. Results obtained in this work show the ability of MIE to detect biological structures and extract their average refractive index and height distribution. One of the major motivations of using MIE in biosensing is to gain more knowledge about the dielectric properties of biological systems for better understanding and further development of biomaterials. Moreover, it constitutes important driving forces in developments directed toward sensors and diagnostic applications for use in clinic medicine.4 Results obtained in this work show that MIE is a valuable detecting technique, which can differentiate bacteria types by examining their morphologies (including shapes and heights) and dielectric properties. Since it is a label-free, noninvasive sensing technique, it does not require any sample preparation, and it can perform measurements fairly fast (within seconds). Because of the above features, MIE can be used to study the interaction mechanism between a substrate and adsorbed biological layer via real-time measurements and analysis. The amount of information obtained from the analysis of MIE measurements can be extended to include the effects due to surface roughness and molecular orientation.4 Finally, comparing this sensing tool to other label-free and noninvasive imaging techniques such as Raman spectroscopy and IR/thermal imaging, MIE provides complementary information of the samples compared to Raman spectroscopy36 and IR/thermal imaging.37 Since Raman spectroscopy measures the vibrational modes of the sample, it is more sensitive to the chemical nature of the sample, while MIE is capable of detecting the distribution of cells and measuring their thickness and dielectric properties. The IR/thermal imaging technique detects the thermal radiation from the sample, which cannot achieve the submicron spatial resolution considered here, and it does not have the capability of extracting detailed structure information as MIE does. Although we report MIE technique only in the visible range, the same experimental setup and analysis can be applicable for other wavelengths such as ultraviolet or infrared.

Appendices

Appendix A:

Derivation of Formulas for the Three Models Used

A.1.

Single-Bounce Light Path

The Fresnel reflection coefficients rp and rs in z-direction for single-bounce light path from Fig. 2(a) are

Eq. (10)

r01p=ErxEix=1ϵ0Kiz1ϵ1Ktz1ϵ0Kiz+1ϵ1Ktz=1ϵ0ϵ0sinθ021ϵ1ϵ1sinθ021ϵ0ϵ0sinθ02+1ϵ1ϵ1sinθ02,

Eq. (11)

r01s=EryEiy=K0z  K1zK0z  +K1z=ϵ0sinθ02ϵ1sinθ02ϵ0sinθ02+ϵ1sinθ02,
where Kiz (i=0,1) denotes the z-component of the wave-vector for light in media i. θ0, ϵ0, and ϵ1 are the AOI, the complex refractive index of air, and the complex refractive index of S. mutans cells, respectively. Substituting for rp and rs in Eq. (2) in Sec. 2.4, we obtain

Eq. (12)

ρ=tanΨeiΔ=(1ϵ0ϵ0sinθ021ϵ1ϵ1sinθ02)(ϵ0sinθ02+ϵ1sinθ02)(1ϵ0ϵ0sinθ02+1ϵ1ϵ1sinθ02)(ϵ0sinθ02ϵ1sinθ02).

Solving Eq. (12) with respect to ϵ1, we have

Eq. (13)

ϵ1=(n1+ik1)2=(1ρ)2(1+ρ)2sinθ02tanθ02+sinθ02.  

A.2.

Two-Bounce Light Path

For the two-bounce light path shown in Fig. 2(b), the Fresnel reflection coefficients is given by

Eq. (14)

rj=r01j+(1r01j2)r12je2iβ,
where j stands for s-polarization and p-polarization, and β represents the phase film factor given as

Eq. (15)

β=2πhn1cosθ1λ,
where h and n1 are the thickness and refractive index of the sample, respectively.

Substituting Eq. (14) for s and p polarizations in Eq. (2), Sec. 2.4, we obtain

Eq. (16)

ρ=tanΨeiΔ=r01p+(1r01p2)r12pγr01s+(1r01s2)r12sγ,
where r01p, r01s, r12p, and r12s have been derived in a previous work by Taya et al.32

Solving Eq. (16) with respect to e2iβ we obtain

Eq. (17)

γ=e2iβ=r01pρr01sρ(1r01s2)r12s(1r01p2)r12p.

Assuming the sample does not absorb light in the visible region (i.e., k1=0), then

Eq. (18)

ln|γ|=ln|r01pρr01sρ(1r01s2)r12s(1r01p2)r12p|=0.

A.3.

Multibounce Light Path

The Fresnel reflection coefficients of multibounce light paths have been derived in Ref. 32 and given as

Eq. (19)

ρ=tanΨeiΔ=(r01p+r12pγ)(r01s+r12sγ)(1+r01sr12sγ)(1+r01pr12pγ).

Assuming that k1=0 in the visible range, then the final solution is written as

Eq. (20)

ln|γ|=ln|e2iβ|=ln|B±B24AC  2A|=0,
where

Eq. (21)

A=(ρr01pr01s)(r12pr12s),

Eq. (22)

B=r01pr01s(ρr12pr12s)+ρr12sr12p,

Eq. (23)

C=ρr01sr01p.

The thickness of the sample can be obtained from Eq. (15) in terms of β, n1, and cosθ1 as

Eq. (24)

h=βλ2πn1cosθ1.

Disclosures

No conflicts of interest, financial or otherwise, are declared by the authors.

Acknowledgments

This work was supported in part by the Ministry of Science and Technology of Taiwan under Contract No. MOST 104-2112-M-001-009-MY2.

References

1. 

Y. C. Chang et al., “Optical sensing of nanoscale objects via computation-aided microscopic imaging ellipsometry,” Sci. Lett. J., 5 229 (2016). Google Scholar

2. 

M. C. Howland et al., “Characterization of physical properties of supported phospholipid membranes using imaging ellipsometry at optical wavelengths,” Biophys. J., 92 (4), 1306 –1317 (2007). http://dx.doi.org/10.1529/biophysj.106.097071 BIOJAU 0006-3495 Google Scholar

3. 

Z. H. Wang and G. Jin, “A label-free multisensing immunosensor based on imaging ellipsometry,” Anal. Chem., 75 (22), 6119 –6123 (2003). http://dx.doi.org/10.1021/ac0347258 ANCHAM 0003-2700 Google Scholar

4. 

H. Arwin, “Spectroscopic ellipsometry and biology: recent developments and challenges,” Thin Solid Films, 313 764 –774 (1998). http://dx.doi.org/10.1016/S0040-6090(97)00993-0 THSFAP 0040-6090 Google Scholar

5. 

H. Arwin, “Ellipsometry on thin organic layers of biological interest: characterization and applications,” Thin Solid Films, 377 48 –56 (2000). http://dx.doi.org/10.1016/S0040-6090(00)01385-7 THSFAP 0040-6090 Google Scholar

6. 

A. J. Qavi et al., “Label-free technologies for quantitative multiparameter biological analysis,” Anal. Bioanal. Chem., 394 (1), 121 –135 (2009). http://dx.doi.org/10.1007/s00216-009-2637-8 ABCNBP 1618-2642 Google Scholar

7. 

M. Losurdo et al., “Spectroscopic ellipsometry and polarimetry for materials and systems analysis at the nanometer scale: state-of-the-art, potential, and perspectives,” J. Nanopart. Res., 11 (7), 1521 –1554 (2009). http://dx.doi.org/10.1007/s11051-009-9662-6 JNARFA 1388-0764 Google Scholar

8. 

P. P. Lehenkari et al., “Adapting atomic force microscopy for cell biology,” Ultramicroscopy, 82 (1–4), 289 –295 (2000). http://dx.doi.org/10.1016/S0304-3991(99)00138-2 ULTRD6 0304-3991 Google Scholar

9. 

A. Rothen, “The ellipsometer, an apparatus to measure thicknesses of thin surface films,” Rev. Sci. Instrum., 16 (2), 26 –30 (1945). http://dx.doi.org/10.1063/1.1770315 RSINAK 0034-6748 Google Scholar

10. 

Y.-D. Chen et al., “Study of biological reaction in cancer cell with spectroscopic imaging ellipsometry,” Proc. SPIE, 9925 992505 (2016). http://dx.doi.org/10.1117/12.2236288 Google Scholar

11. 

Y. M. Bae et al., “Detection of insulin-antibody binding on a solid surface using imaging ellipsometry,” Biosens. Bioelectron., 20 (4), 895 –902 (2004). http://dx.doi.org/10.1016/j.bios.2004.03.032 BBIOE4 0956-5663 Google Scholar

12. 

I. An, “Application of imaging ellipsometry to the detection of latent fingermarks,” Forensic Sci. Int., 253 28 –32 (2015). http://dx.doi.org/10.1016/j.forsciint.2015.05.009 FSINDR 0379-0738 Google Scholar

13. 

M. M. B. Nielsen and A. C. Simonsen, “Imaging ellipsometry of spin-coated membranes: mapping of multilamellar films, hydrated membranes, and fluid domains,” Langmuir, 29 (5), 1525 –1532 (2013). http://dx.doi.org/10.1021/la3046675 LANGD5 0743-7463 Google Scholar

14. 

K. H. Metwalli et al., “Streptococcus mutans, Candida albicans, and the human mouth: a sticky situation,” PLoS Pathog., 9 (10), e1003616 (2013). http://dx.doi.org/10.1371/journal.ppat.1003616 1553-7374 Google Scholar

15. 

D. Ajdic et al., “Genome sequence of Streptococcus mutans UA159, a cariogenic dental pathogen,” Proc. Natl. Acad. Sci. U. S. A., 99 (22), 14434 (2002). http://dx.doi.org/10.1073/pnas.172501299 Google Scholar

16. 

N. Hakimiha et al., “The susceptibility of Streptococcus mutans to antibacterial photodynamic therapy: a comparison of two different photosensitizers and light sources,” J. Appl. Oral Sci., 22 (2), 80 –84 (2014). http://dx.doi.org/10.1590/1678-775720130038 Google Scholar

17. 

Q. Wei and L. Y. Z. Ma, “Biofilm matrix and its regulation in Pseudomonas aeruginosa,” Int. J. Mol. Sci., 14 (10), 20983 (2013). http://dx.doi.org/10.3390/ijms141020983 1422-0067 Google Scholar

18. 

B. H. Liu et al., “In situ biosensing of the nanomechanical property and electrochemical spectroscopy of Streptococcus mutans-containing biofilms,” J. Phys. D Appl. Phys., 46 (27), 275401 (2013). http://dx.doi.org/10.1088/0022-3727/46/27/275401 Google Scholar

19. 

S. E. Cross et al., “Atomic force microscopy study of the structure-function relationships of the biofilm-forming bacterium Streptococcus mutans,” Nanotechnology, 17 (4), S1 –S7 (2006). http://dx.doi.org/10.1088/0957-4484/17/4/001 NNOTER 0957-4484 Google Scholar

21. 

E. H. Beachey, Bacterial Adherence, 6 109 Chapman and Hall Ltd., London (1980). Google Scholar

22. 

“Creative life science,” http://www.cmp-micro.com Google Scholar

23. 

Y. F. Hu, J. D. Zhang and J. Ulstrup, “Interfacial electrochemical electron transfer processes in bacterial biofilm environments on Au(111),” Langmuir, 26 (11), 9094 –9103 (2010). http://dx.doi.org/10.1021/la9047853 LANGD5 0743-7463 Google Scholar

26. 

J.-Y. Fan et al., “A study on transmitted intensity of disturbance for air-spaced Glan-type polarizing prisms,” Opt. Commun., 223 11 –16 (2003). http://dx.doi.org/10.1016/S0030-4018(03)01618-3 OPCOB8 0030-4018 Google Scholar

27. 

Fianium, “WhiteLase SC series: blue & UV enhanced supercontinuum fiber lasers,” (2017) http://www.fianium.com/pdf/WhiteLase_SC400_UV_v1.pdf January ). 2017). Google Scholar

28. 

Edmund Optics Inc, “Mitutoyo infinity corrected long working distance objectives,” 2017, (2017) https://www.edmundoptics.com.tw/microscopy/infinity-corrected-objectives/Mitutoyo-Infinity-Corrected-Long-Working-Distance-Objectives/ March ). 2017). Google Scholar

29. 

PCO, “Pixelfly qe,” (2017) https://www.pco.de/sensitive-cameras/pixelfly-qe/ March 2017). Google Scholar

30. 

A. Velten et al., “Recovering three-dimensional shape around a corner using ultrafast time-of-flight imaging,” Nature Commun., 3 745 (2012). http://dx.doi.org/10.1038/ncomms1747 NCAOBW 2041-1723 Google Scholar

31. 

C.-Y. Tsai, A. Veeraraghavan and A. C. Sankaranarayanan, “Shape and reflectance from two-bounce light transients—IEEE Xplore document,” (2016) http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=7492882 Google Scholar

32. 

S. A. Taya, T. M. El-Agez and A. A. Alkanoo, “Thin film characterization using rotating polarizer analyzer ellipsometer with a speed ratio 1:3,” J. Electromagn. Anal. Appl., 3 351 –358 (2011). http://dx.doi.org/10.4236/jemaa.2011.39056 Google Scholar

33. 

MathWorks, “Matlab,” (2017) https://www.mathworks.com/products/matlab.html March ). 2017). Google Scholar

34. 

Olympos, “Micro cantilevers,” OMCL-AC160TS, (2017) http://probe.olympus-global.com/en/product/omcl_ac160ts_r3/ March ). 2017). Google Scholar

35. 

R. Synowicki, “Optical characterization by spectroscopic ellipsometry,” 96 (2010) http://grover.mirc.gatech.edu/data/gatech_ellipsometry_seminar-dec_2010.pdf Google Scholar

36. 

R. Lewandowska et al., “What can Raman spectroscopy and spectroscopic ellipsometry bring for the characterisation of thin films and materials surface?,” Surf. Interface Anal., 40 (3–4), 588 –592 (2008). http://dx.doi.org/10.1002/(ISSN)1096-9918 Google Scholar

37. 

A. C. Jones and M. B. Raschke, “Thermal infrared near-field spectroscopy,” Nano Lett., 12 (3), 1475 –1481 (2012). http://dx.doi.org/10.1021/nl204201g NALEFD 1530-6984 Google Scholar

Biography

Mai Ibrahim Khaleel is a PhD student in the nanoscience and technology program—Taiwan International Graduate Program—in collaboration with the Department of Engineering and System Science at National Tsing Hua University. She received her BS and MS degrees in physics from Al-Quds University in 2007 and 2009, respectively. Her current research interests include bio-optics, imaging ellipsometry, and photonics.

Biographies for the other authors are not available.

© 2017 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2017/$25.00 © 2017 SPIE
Mai Ibrahim Khaleel, Yu-Da Chen, Ching-Hang P. Chien, and Yia-Chung Chang "Sensing of Streptococcus mutans by microscopic imaging ellipsometry," Journal of Biomedical Optics 22(5), 056005 (16 May 2017). https://doi.org/10.1117/1.JBO.22.5.056005
Received: 29 January 2017; Accepted: 21 April 2017; Published: 16 May 2017
Lens.org Logo
CITATIONS
Cited by 1 scholarly publication.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Ellipsometry

Statistical analysis

Atomic force microscopy

Glasses

Optical testing

Specular reflections

Back to Top