<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JCC</journal-id><journal-title-group><journal-title>Journal of Computer and Communications</journal-title></journal-title-group><issn pub-type="epub">2327-5219</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jcc.2021.912010</article-id><article-id pub-id-type="publisher-id">JCC-114424</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Terahertz Bessel Beam Applied to Thickness Measurement of Ellipsometry Methods
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Siyu</surname><given-names>Tu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Kejia</surname><given-names>Wang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Jinsong</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Zhengang</surname><given-names>Yang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan, China</addr-line></aff><pub-date pub-type="epub"><day>21</day><month>12</month><year>2021</year></pub-date><volume>09</volume><issue>12</issue><fpage>125</fpage><lpage>132</lpage><history><date date-type="received"><day>28,</day>	<month>August</month>	<year>2021</year></date><date date-type="rev-recd"><day>28,</day>	<month>December</month>	<year>2021</year>	</date><date date-type="accepted"><day>31,</day>	<month>December</month>	<year>2021</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  
    This paper introduces a Terahertz (THz) ellipsometer thickness measurement method based on Bessel beams. The ellipsometry method is used to measure the thickness of film in the THz band. And the thickness of film could be measured in the range of hundreds of microns which the measurement of film thickness by ellipsometer is usually only a few hundred nm in visible light. In addition, the photon energy of the THz-wave is very low and has little damage to the film. The THz Bessel beam has the characteristics of self-healing and diffraction-free. It can make the film thickness measurement within the diffraction-free distance, which is conducive to the flexibility of the ellipsometer system. We use a multi-frequency method to measure film with two-dimensional finite different time domain (FDTD) to numerical simulation analysis of light intensity successfully. 
  
 
</p></abstract><kwd-group><kwd>Bessel-Beam</kwd><kwd> Ellipsometry</kwd><kwd> Diffraction-Free</kwd><kwd> THz</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Terahertz (THz) waves are electromagnetic waves with a frequency from 0.1 THz to 10 THz and a wavelength from 30 um to 3 mm, which are between microwaves and infrared waves. THz waves have the characteristics of both microwaves and infrared waves. Macroscopically, THz waves can penetrate many dielectric materials and non-polar liquids. Microscopically, THz waves contain rich sample spectral information due to their broad spectrum [<xref ref-type="bibr" rid="scirp.114424-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.114424-ref2">2</xref>]. The unique characteristics of THz waves have brought far-reaching influence on communication (broadband communication), radar, electronic countermeasures, electromagnetic weapons, astronomy, medical imaging, non-destructive testing, safety inspection [<xref ref-type="bibr" rid="scirp.114424-ref3">3</xref>]-[<xref ref-type="bibr" rid="scirp.114424-ref8">8</xref>].</p><p>With thin film technology in information storage, electronic components, aerospace technology and optical instruments, etc. The accurate measurement of optical properties and surface morphology of film has gradually become an important direction of film research. Moreover, the physical properties of the film are closely related to the optical properties and thickness of the film. Therefore, the precise measurement of the optical properties and thickness of the film has become the primary research issue. The current methods for measuring the refractive index and thickness of the film include ellipsometry, heterodyne interferometry, equal thickness interferometry, reflection spectroscopy, transmission spectroscopy and other measurement methods [<xref ref-type="bibr" rid="scirp.114424-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.114424-ref10">10</xref>]. Ellipsometry is a non-contact and non-destructive measurement, which can measure samples under vacuum, atmosphere, moisture, and harmful environments with low environmental requirements [<xref ref-type="bibr" rid="scirp.114424-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.114424-ref12">12</xref>]. In addition, the ellipsometry is simple to operate and has a short measurement period. It can measure multiple physical quantities (film thickness, film refractive and absorption coefficient) at the same time. Hence, ellipsometry is widely used to measure films.</p><p>The Bessel beam has diffraction-free characteristics, the intensity distribution remains unchanged during propagation and has a high degree of localization [<xref ref-type="bibr" rid="scirp.114424-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.114424-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.114424-ref15">15</xref>]. It is widely used in precision measurement, laser processing and other fields in the optical waveband. Th diffraction-free and long focal depth of Bessel beams could be further used in THz there-dimensional (3D) tomography systems, such as THz CT systems [<xref ref-type="bibr" rid="scirp.114424-ref16">16</xref>]. The THz imaging system based on Bessel beam has unique advantages in applications such as long-distance standing human security inspection and non-destructive inspection of thick sample.</p><p>The film thickness measured by ellipsometer in the optical band is mainly in the range of hundreds of nanometers. For samples with a thickness of millimeters and hundreds of microns, the visible light waveband cannot be measured. Therefore, we proposed the application of THz waves to accurately measure the thickness of samples in the order of millimeters and hundreds of microns. In addition, the photon energy of the THz wave is relatively low, which hardly cause damage to the sample. The THz Bessel beam has the characteristics of diffraction-free and self-healing, and can realize precise measurement of samples within its diffraction-free distance. Thus, an ellipsometer thickness measurement method based on THz Bessel beam is proposed in this paper. This paper is organized as follows. In Section 2, the theory of ellipsometry calculation of film thickness and the generation of THz Bessel beams are showed. In Section 3, numerical simulation calculation of the thickness of film measured by THz Bessel beam is carried out by two-dimensional finite difference time domain (FDTD). Finally, our conclusion is summarized in Section 4.</p></sec><sec id="s2"><title>2. Theory</title><sec id="s2_1"><title>2.1. The Generation of THz Bessel Beam</title><p>A Gaussian beam with a waist radius of ω 0 is normally incident on the front surface of the axicon [<xref ref-type="bibr" rid="scirp.114424-ref17">17</xref>]. By geometric optics analysis, the Gaussian beam is regarded as a straight line, in <xref ref-type="fig" rid="fig1">Figure 1</xref>. According to refraction law:</p><p>n 0 sin ( θ ) = n sin ( γ ) ,</p><p>α 0 + γ = θ , τ + 2 γ = π ,</p><p>the refractive index of atmosphere is n 0 = 1 :</p><p>θ = arcsin [ n n 0 cos ( τ / 2 ) ] = arcsin [ n cos ( τ / 2 ) ] ,</p><p>the angle between the Gaussian beam and the z-axis after refractive:</p><p>α 0 = arcsin [ n cos ( τ / 2 ) ] + ( τ − π ) / 2 ,</p><p>when γ is small:</p><p>θ = n γ ,</p><p>α 0 = ( n − 1 ) γ .</p><p>Therefore, the Bessel beam with a diffraction-free distance of z max = ω 0 / tan ( α 0 ) is generated.</p></sec><sec id="s2_2"><title>2.2. Ellipsometry Measurement Theory</title><p>The basic principle of ellipsometer measurement is to incident a beam of elliptic polarization light onto the sample to be tested at a certain angle of incidence. Since the sample has different reflection and transmission coefficients for the p-light and s-light of the incident light, the polarization state of the light reflected from the sample changes with respect to the incident light which is showed in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The variation of polarization state is closely related to the structure of sample. The specific structure of sample can be determined according to the variation of polarization state. The variation of the polarization state of the light before and after reflection is observed to determine the reflection coefficient ratio ρ, which the information related to the optical parameters of the sample could be obtained.</p><p>The electric vector of incident light can be decomposed into p-light and s-light. E i p and E i s represents the p-light and s-light of the incident light respectively. E r p and E r s represents the p-light and s-light of the incident light respectively.</p><p>r p = E r p / E i p</p><p>r s = E r s / E i s</p><p>The Fresnel formula for the reflection of wave on the interface:</p><p>r p = n 2 cos φ 1 − n 2 cos φ 2 n 2 cos φ 1 + n 1 cos φ 2</p><p>r s = n 1 cos φ 1 − n 2 cos φ 2 n 1 cos φ 1 + n 2 cos φ 2</p><p>The relationship between the complex amplitude of the total reflected wave and the incident wave is:</p><p>E r p = r 1 p + r 2 p exp ( − i 2 δ ) 1 + r 1 p r 2 p exp ( − i 2 δ ) E i p</p><p>E r s = r 1 s + r 2 s exp ( − i 2 δ ) 1 + r 1 s r 2 s exp ( − i 2 δ ) E i s</p><p>The total reflection coefficients of the film are:</p><p>R p = E r p / E i p</p><p>R p = E r p / E i p</p><p>The reflectance coefficient:</p><p>ρ = R p / R s = E r p / E i p E r s / E i s = | E r p / E i p E r s / E i s | exp { i [ ( β r p − β i p ) − ( β r s − β i s ) ] } = tan ψ exp ( i Δ )</p><p>According to the Fresnel formula and the law of refraction:</p><p>{ r 1 p = ( n 2 cos φ 1 − n 1 cos φ 2 ) / ( n 2 cos φ 1 + n 1 cos φ 2 ) r 2 p = ( n 3 cos φ 2 − n 2 cos φ 3 ) / ( n 3 cos φ 2 + n 2 cos φ 3 ) r 1 s = ( n 1 cos φ 1 − n 2 cos φ 2 ) / ( n 1 cos φ 1 + n 2 cos φ 2 ) r 2 s = ( n 2 cos φ 2 − n 3 cos φ 3 ) / ( n 2 cos φ 2 + n 3 cos φ 3 ) n 1 sin φ 1 = n 2 sin φ 2 = n 3 sin φ 3 2 δ = ( 4 π d n 2 cos φ 2 ) / λ</p><p>Therefore, the relevant parameters of the film can be obtained by the variation of the polarization state of the incident beam before and after reflection.</p></sec></sec><sec id="s3"><title>3. Simulation Results</title><p>The 2D FDTD (Finite Difference Time Domain) numerical stimulation method is utilized to conduct a simulation on the model. In <xref ref-type="fig" rid="fig3">Figure 3</xref>, a Gaussian beam with a waist radius of 20 mm is incident on an axicon with base angle of 15˚. The Bessel beam is generated by the axicon. The sample is within the diffraction-free distance of the Bessel beam. And the Bessel beam is incident on the film obliquely at an incident angle of 45˚. A line monitor is set in the reflection direction of the beam.</p><p>The basic setup of the model is as follows. Since the thickness of the measured film is related to the polarization state of the incident beam, the polarization state is specially set when setting Gaussian beam source. We choose to use circularly polarization beam as the incident light. In FDTD, a Gaussian beam source can only set a certain angle of linear polarization light. Therefore, we set up two Gaussian beam sources. The amplitudes of two sources are same. Their polarization angles are 0˚and 90˚respectively as the p-light and s-light. The phase difference between the two Gaussian beams is 90˚. In Solidworks software, an isosceles triangle with a base angle of 15˚ and a bottom margin of three inches is created. The isosceles triangle is rotated 360˚around its center line to form axicon. The axicon is exported from Solidworks in STL format and imported it to FDTD. A sample with film thickness of 1 mm, refractive index n = 1.6 , substate thickness of 20 mm, refractive index n = 3.4 was established in FDTD. The position of the Gaussian sources, axicon, sample, monitor are respectively designed in the appropriate position.</p><p>In <xref ref-type="fig" rid="fig4">Figure 4</xref>, it is showed that the Gaussian beam is incident on the axicon to produce a Bessel beam. And the Bessel beam is incident on the film, part of the Bessel beam is reflected, part of the beam pass through the substrate. <xref ref-type="fig" rid="fig5">Figure 5</xref> are the light intensity measured by the line monitor. It can be seen that the strong middle peak is the reflected beam of the main lobe of the Bessel beam. <xref ref-type="fig" rid="fig5">Figure 5</xref>(a) &amp; <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) is the p-light and s-light of the Bessel beam generated by the axicon respectively, which is also the incident light. <xref ref-type="fig" rid="fig5">Figure 5</xref>(c) &amp; <xref ref-type="fig" rid="fig5">Figure 5</xref>(d) is</p><p>the p-light and s-light reflected by film. The ellipsometric parameters can be calculated by the light intensity of the incident and reflected p-light and s-light of simulation. The p-light to s-light ratio of incident light is 0.9959, the p-light to s-light of reflected light is 0.6524. According to calculation, the r<sub>1p</sub> is 0.1155, the r<sub>2p</sub> is 0.3218, the r<sub>1s</sub> is −0.3399, the r<sub>2s</sub> is −0.3971, the ellipsometric parameters is tan ψ = 0.6550 62567862263 . By:</p><p>2 δ = ( 4 π d n 2 cos φ 2 ) / λ ,</p><p>the Periodic film thickness is obtained.</p><p>The Periodic film thickness value is as follows:</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Part of the periodic film thickness</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >d/mm</th><th align="center" valign="middle" >0.7558</th><th align="center" valign="middle" >0.8322</th><th align="center" valign="middle" >1.0279</th><th align="center" valign="middle" >1.2897</th><th align="center" valign="middle" >1.5852</th><th align="center" valign="middle" >1.8987</th><th align="center" valign="middle" >2.2226</th></tr></thead></tbody></table></table-wrap><p><xref ref-type="table" rid="table1">Table 1</xref> gives part of the Periodic film thickness. When setting up model in simulation, the film thickness is 1 mm. Therefore, the value of the calculated periodic film thickness closest to 1 mm is the actual value of the film thickness. The 1.0279 mm is closest to 1 mm, and there is only 0.02 mm error which proves that the Bessel beam model is correct.</p></sec><sec id="s4"><title>4. Conclusion</title><p>This paper proposed a method of measuring film thickness by THz Bessel beam. It is proved that precise thickness measurement of thin films of millimeters and hundreds of microns can be achieved in the THz band. And the THz Bessel beam is a diffraction-free beam, which can realize the measurement system with a long depth of field. We have employed the ellipsometric calculation method in FDTD to verify the feasibility of this theory. In the future, we will further implement ellipsometric thickness measurement in the THz experimental system.</p></sec><sec id="s5"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s6"><title>Cite this paper</title><p>Tu, S.Y., Wang, K.J., Liu, J.S. and Yang, Z.G. (2021) Terahertz Bessel Beam Applied to Thickness Measurement of Ellipsometry Methods. Journal of Computer and Communications, 9, 125-132. https://doi.org/10.4236/jcc.2021.912010</p></sec></body><back><ref-list><title>References</title><ref id="scirp.114424-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Redo-Sanchez, A.S., et al. (2016) Terahertz Time-Gated Spectral Imaging for Content Extraction through Layered Structures. Nat Commun, 7, 12665.  
https://doi.org/10.1038/ncomms12665</mixed-citation></ref><ref id="scirp.114424-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Sun, W., Wang, X. and Zhang, Y. (2012) A Method to Monitor the Oil Pollution in Water with Reflective Pulsed Terahertz Tomography. Optik, 123, 1980-1984.  
https://doi.org/10.1016/j.ijleo.2011.10.002</mixed-citation></ref><ref id="scirp.114424-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Hu, B.B. and Nuss, M.C. (1995) Imaging with Terahertz Waves. Optics Letters, 20, 1716. https://doi.org/10.1364/OL.20.001716</mixed-citation></ref><ref id="scirp.114424-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Wu, Q., Hewitt, T.D. and Zhang, X.C. (1996) Two-Dimensional Electro-Optic Imaging of THz Beams. Applied Physics Letters, 69, 1026-1028.  
https://doi.org/10.1063/1.116920</mixed-citation></ref><ref id="scirp.114424-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Lui, H.S., Taimre, T., Bertling, X., et al. (2014) Terahertz Inverse Synthetic Aperture Radar Imaging Using Self-Mixing Interferometry with a Quantum Cascade Laser. Optics Letters, 39, 2629-2632. https://doi.org/10.1364/OL.39.002629</mixed-citation></ref><ref id="scirp.114424-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Song, H.J. and Nagatsuma, T. (2011) Present and Future of Terahertz Communications. IEEE Transactions on Terahertz Science and Technology, 1, 256-263.  
https://doi.org/10.1109/TTHZ.2011.2159552</mixed-citation></ref><ref id="scirp.114424-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Seeds, A.J., Shams, H., Fice, M.J., et al. (2015) Terahertz Photonics for Wireless Communications. Journal of Lightwave Technology, 33, 579-587.  
https://doi.org/10.1109/JLT.2014.2355137</mixed-citation></ref><ref id="scirp.114424-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Shen, Y.C. (2007) Terahertz Pulsed Spectroscopy and Imaging for Pharmaceutical Applications: A Review. J Pharm Pharmacol, 59, 209.</mixed-citation></ref><ref id="scirp.114424-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Paulick, T.C. (1986) Inversion of Normal-Incidence (R,T) Measurements to Obtain n+ik for Thin Films. Applied Optics, 25, 562-564.  
https://doi.org/10.1364/AO.25.000562</mixed-citation></ref><ref id="scirp.114424-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Chiu, M., Lee, J. and Su, D. (1999) Complex Refractive-Index Measurement Based on Fresnel’s Equations and the Use of Heterodyne Interferometry. Applied Optics, 38, 4047-4052. https://doi.org/10.1364/AO.38.004047</mixed-citation></ref><ref id="scirp.114424-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Kravets, V.G., Grigorenko, A.N., et al. (2010) Spectroscopic Ellipsometry of Graphene and an Exciton-Shifted van Hove Peak in Absorption. Physical Review B, 81, 1-6. https://doi.org/10.1103/PhysRevB.81.155413</mixed-citation></ref><ref id="scirp.114424-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Li, J., Ramanujam, B. and Collins, R.W. (2011) Dual Rotating Compensator Ellipsometry: Theory and Simulations. Thin Solid Films, 519, 2725-2729.  
https://doi.org/10.1016/j.tsf.2010.11.075</mixed-citation></ref><ref id="scirp.114424-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Cerda, S.C., Nava, M.M. and Hickmann, J.M. (1998) Interference of Traveling Nondiffracting Beams. Optics Letters, 23, 1871-1873.  
https://doi.org/10.1364/OL.23.001871</mixed-citation></ref><ref id="scirp.114424-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Durnin, J., Miceli Jr., J. and Eberly, J. (1987) Diffraction-Free Beams. Physical Review Letters, 58, 1499. https://doi.org/10.1103/PhysRevLett.58.1499</mixed-citation></ref><ref id="scirp.114424-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Scott, G. and Mcardle, N.J.O.E. (1992) Efficient Generation of Nearly Diffraction-Free Beams Using an Axicon. 31, 2640-2643. https://doi.org/10.1117/12.60017</mixed-citation></ref><ref id="scirp.114424-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Bitman, A., Goldring, S., Moshe, I., et al. (2014) Computed Tomography Using Broadband Bessel THz Beams and Phase Contrast. Optics Letters, 39, 1925-1928.  
https://doi.org/10.1364/OL.39.001925</mixed-citation></ref><ref id="scirp.114424-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Monk, S., Arlt, J., Robertson, D.A., et al. (1999) Generation of Bessel Beams at Millimeter Waves Frequencies by Use of an Axicon. Optics Communications, 170, 213-215. https://doi.org/10.1016/S0030-4018(99)00463-0</mixed-citation></ref></ref-list></back></article>