Next Article in Journal
Three-Dimensional Imaging via Time-Correlated Single-Photon Counting
Next Article in Special Issue
Optical and Galvanomagnetic Properties of Bi1-xSbx Thin Films in the Terahertz Frequency Range
Previous Article in Journal
An Overview on the Latest Nature-Inspired and Metaheuristics-Based Image Registration Algorithms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

FEM Simulation of THz Detector Based on Sb and Bi88Sb12 Thermoelectric Thin Films

by
Anastasiia S. Tukmakova
1,*,
Alexei V. Asach
1,
Anna V. Novotelnova
1,
Ivan L. Tkhorzhevskiy
1,
Natallya S. Kablukova
2,3,
Petr S. Demchenko
4,
Anton D. Zaitsev
4 and
Mikhail K. Khodzitsky
4
1
Faculty of Cryogenic Engineering, ITMO University, 197101 Saint-Petersburg, Russia
2
International Scientific and Research Institute of Bioengineering, ITMO University, 197101 Saint-Petersburg, Russia
3
The Department of Physics, Saint Petersburg State University of Industrial Technologies and Design, 191186 Saint-Petersburg, Russia
4
Terahertz Biomedicine Laboratory, ITMO University, 197101 Saint-Petersburg, Russia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(6), 1929; https://doi.org/10.3390/app10061929
Submission received: 17 February 2020 / Revised: 6 March 2020 / Accepted: 6 March 2020 / Published: 11 March 2020
(This article belongs to the Special Issue Terahertz Sensing)

Abstract

:
A terahertz (THz) detector based on thermoelectric thin films was simulated using the finite elements method. The thermoelectric circuit consisted of S b and B i 88 S b 12 150-nm films on the mica substrate. S b , B i 88 S b 12 , and mica-substrate properties have been measured experimentally in the THz frequency range. The model of electromagnetic heating was used in order to estimate possible heating of S b - B i 88 S b 12 contact. THz radiation power varied from 1 μ W to 50 mW, and frequency varied in the range from 0.3 to 0.5 THz. The calculations showed a temperature difference of up to 1 K, voltage up to 0.1 mV, and responsivity of several mVW 1 . The results show that thin S b and B i S b thermoelectric films can be used for THz radiation detection at room temperatures.

Graphical Abstract

1. Introduction

Terahertz (THz) devices are of significant importance for such fields as disease diagnostics and medical imaging [1,2,3], security systems [4,5], and sensing and imaging [6]. Thermal THz sensors, such as bolometers, are used for THz detection [7]. The effect of electromagnetic heating due to THz radiation absorption can be used as the signal for thermoelectric (TE) detector. TE detectors work on the basis of the Seebeck effect, which generates electric voltage in TE circuits in case of temperature difference between two junctions (hot and cold junctions) of TE materials [8,9]. Basically, the temperature of cold junction is maintained at room temperature, while hot junction temperature is to be increased due to the electromagnetic heating. In literature, the effect combining optical and thermoelectric phenomena is called photothermoelectric, and it has been studied with the use of different advanced materials: graphene [10], single-layered M o S 2 [11], B i 2 S e 3 nanoribbons [12], and organic thermoelectrics [13]. A line of works about THz detectors working on the TE effect have been published. In [14], a responsivity of several hundred VW 1 was reported for a detector based on single-layer graphene on 300-nm S i O 2 substrate. A THz TE detector based on carbon nanotubes with responsivity of several VW 1 was reported in [15].
The increase of efficiency or responsivity of TE detector is possible due to the increase of temperature difference between junctions or due to the increase of TE figure of merit Z. Z is proportional to Seebeck coefficient squared S 2 and electrical conductivity σ ; Z is inversely proportional to thermal conductivity κ . In order to increase temperature difference, the cross-section of hot junction can be decreased. The hot junction can be thermally insulated, for example, in vacuum, and the device can be encapsulated in order to prevent convection.
Besides of temperature difference increase, the usage of low-dimensional TE materials interconnects with Z-value improvement. It was shown that 2D and 1D thermoelectrics could give a significant increase of TE figure of merit [16]. In particular, TE thin films are known for the possibility to produce a high Seebeck coefficient and power factor S 2 · σ . As all thermoelectric properties are the functions of temperature, not every material will be suitable for detecting purposes requiring near-room temperatures. TE films based on B i and T e alloys are good candidates for the applications at room temperature, their properties have been reported in [17]. Other candidates are TE films based on B i and S b solid solutions. In a B i ( 1 x ) S b ( x ) system with Sb concentration lower than 15%, the Seebeck coefficient value at room temperatures is around −100 μ V/K. This system was studied in the works of Volklein et al. in the 1990s [18] and Suslov et al. in 2017 [19]. Close values of Seebeck coefficient have been reported for the B i 0.91 S b 0.09 system in [20]. In an article published by Linseis et al. in 2018 [21], a Seebeck coefficient of B i 87 S b 13 was around 80–110 μ V·K 1 with comparatively low thermal conductivity (1.5–3.5 W·m 1 K 1 ) and high electrical conductivity (around 100–300 kS·m 1 ) at room temperatures. Sb thin film was used in our study as a possible p-type material, which can have a good agreement in electrical conductivity with Bi-Sb solid solutions. Seebeck coefficient of Sb is around 30 μ V/K [22]. Moreover, the presence of a narrow band gap makes it possible for Bi-Sb to absorb THz radiation. The band structure of Bi-Sb may be seen in [23], the results have been taken from the work published in 1970 [24]. Surface band structure of B i ( x 1 ) S b ( x ) was reported in [25]. Hence, from the point of view of the highest Seebeck coefficient at room temperature and band structure, Bi-Sb thin films potentially can be used as a sensitive element for THz detection.
The bismuth-antimony thermocouple has already been considered as a possible basis for far-infrared detection in 1982 and showed relatively low noise equivalent power. In [26], an antenna based on n-type B i 0.87 S b 0.13 and p-type Sb thin films on silicon substrate was proposed for 0.812-THz frequency detection. The simulations showed a possible temperature difference of up to 7.3 K in vacuum sealing under the input power of 7.5 μ W. The responsivity was 30 V / W in the air and 119 V / W in vacuum for a pixel with 8 thermocouples. In [27], an antenna made of 8 S b B i 87 S b 13 thermocouples was fabricated with 30- μ m length and 2- μ m width. It showed an output voltage of 0.1–0.7 μ V. The main advantage of this antenna was the possibility of room-temperature operation, low noise equivalent power (from 100 to 1000 pW·Hz 0.5 ), and a fast response of 22 μ s.
Finite elements simulation method showed itself as a reasonable and effective tool for THz antenna modelling, making it possible to combine thermal, electrical, and optical phenomena in one multiphysical model. For example, in [28], THz-resonant antenna was simulated to investigate the photocurrent value as the function on the antenna length. The simulation results were in good agreement with the experimental measurements. In [26], an electromagnetic heating model was used in order to evaluate temperature increase due to THz radiation absorption. A line of works about FEM simulation of thermoelectric detectors of infrared radiation, which is close to THz range, has been presented in the literature [29].
The aim of this work is to simulate a THz detector based on Sb and B i 88 S b 12 TE thin films in order to evaluate possible heating of hot junction due to the electromagnetic heating at 0.3, 0.4, and 0.5 THz. Such a detector can be potentially used as a pixel of THz imaging system.

2. Materials and Methods

2.1. Model Geometry and Mesh

3D model geometry is presented in Figure 1. Cross section in x-z plane is presented in Figure 2 and gives a detailed representation of antenna geometry. Position numbers (1–6) in both figures correspond to the same elements. Antenna had a square x-y cross-section 1 × 1 mm 2 . We used free tetrahedral user-controlled mesh, as the correlation between maximum mesh element size and wavelength must be controlled. Maximum element size was set equal to 40 μ m.

2.2. Mathematical Model

Interaction between detector and THz radiation was described by 2 interfaces: (1) electromagnetic heating and (2) heat transfer in solids. Radiation was excited from the port, passed through the vacuum, absorbed by TE films, and caused heating. The resulting temperature is the dependent variable to be found.

2.2.1. General Equations

Wave equation:
× μ r 1 ( × E ) k 0 2 ( ϵ r i σ / ( ω ϵ 0 ) ) E = 0 ,
where μ r is the relative permeability; E is the electric field; k 0 is the wave number in free space; ϵ r and ϵ 0 are the relative and free space permittivity, respectively; σ is the electrical conductivity; and ω is the angular frequency.
For the electric field displacement description, the dielectric loss model was used:
ϵ r = ϵ i ϵ ,
where ϵ and ϵ are real and imaginary parts of permittivity.
Heat transfer:
q = κ T = Q e ,
where q is the heat flux and Q e is the electromagnetic heat source, T is the temperature, κ —thermal conductivity.
Electromagnetic heat source or electromagnetic loss density includes resistive and magnetic losses:
Q e = Q r h + Q m l = 1 2 R e ( J · E * ) + 1 2 R e ( i ω B · H * ) ,
where J is the electric current density, E is the electric field strength, B is the magnetic induction, and H is the magnetic field strength.

2.2.2. Boundary Conditions

The excitation port 1 (excitation port on Figure 2) was applied to the upper surface of the simulated system (1 × 1 mm 2 ) and was described by the following equation:
S = δ Ω ( E E 1 ) · E 1 δ Ω E 1 · E 1 ,
where E 1 and E 2 are electric fields on port 1 and port 2, respectively.
Port 2 (“off-type” port on Figure 2) was applied to the whole lower surface of the detector (1 × 1 mm 2 ) corresponding to the mica-substrate downside. The “off-type” port was described by the following equation:
S = δ Ω E · E 2 δ Ω E 2 · E 2 .
The scattering boundary condition was applied to the lateral surfaces:
n × ( × ( E ) ) i k n × ( E × n ) = 0 ,
where n-normal vector.

2.2.3. Mathematical Description of Thermal and Electromagnetic Phenomena in Thin Films

Four thin-layer objects were in the model: TE thin films of S b and B i 88 S b 12 , and two copper films. The thickness of all films was 150 nm. In order to reduce the model size and computational time, analytical expressions were used to described thin objects.
In electromagnetic interfaces, a dielectric loss model is required. For this purpose, a transition boundary condition was used:
J s . up = Z s E t . u p Z t E t . d o w n Z s 2 Z t 2 ,
J s . down = Z s E t . d o w n Z t E t . u p Z s 2 Z t 2 ,
Z s = i ω μ k 1 t a n ( k d ) ,
Z t = i ω μ k 1 s i n ( k d ) ,
k = ω ϵ r + ( σ / ( i ω ) ) ) μ r ,
where indices “ u p ” and “ d o w n ” correspond to upsides and downsides of the layer, Z s —surface impedance, Z t —tangential impedance, J s —surface current density, and E t —tangential electric field.
In the heat transfer interface, a thin layer description was used in order to describe heat transfer in thin objects:
t q s = d s Q s + q 0 ,
q s = d s κ t T ,
where d s is film thickness, Q s is a density-distributed heat source, q 0 —heat flux received out-of-plane, q s —heat flux in the film, and t stands for tangential derivative.
Thermally thin approximation model was used for the layer description. This model is valid for films with high thermal conductivity, or conductivity which is higher than the surroundings.

2.3. Materials

The films were prepared by the thermal vacuum deposition with the subsequent annealing. This method was used for B i x 1 S b x -films preparation and was described in [30]. The properties of S b and B i 88 S b 12 films of 150-nm thickness have been chosen for simulation. Electrical conductivity and relative permittivity of mica substrate, S b , and B i 88 S b 12 thin films in terahertz frequency range were measured using a terahertz time-domain spectrometer (THz-TDS) in thin-film approximation [31] (when a wavelength of THz radiation is by orders of magnitude higher than the film thickness). The waveforms of THz single-pulse (temporal dependencies of THz signal electric field amplitude) were measured for an air (reference signal), mica substrate, and a film. The corresponding amplitude and phase spectra were obtained in the 0.3–0.8 THz frequency range through the Fourier transform method. The complex sheet conductivity of a film was calculated as
σ ^ ( f ) = [ ( n ^ s u b ( f ) + 1 ) E ^ s u b ( f ) / E ^ s a m ( f ) n ^ s u b ( f ) 1 ] / Z 0 ,
where f is the frequency of THz radiation, Z 0 is the free-space impedance, E ^ s u b and E ^ s a m are complex amplitudes of THz wave transmitted through the substrate and the sample, and n ^ s u b is the substrate complex refractive index—which is calculated as
n ^ s u b ( f ) = c [ ϕ s u b ( f ) ϕ a i r ( f ) ] / ( 2 π f d s u b ) + 1 i c l n ( | E ^ s u b 2 ( f ) | / | E ^ a i r 2 ( f ) | ) / ( 4 π f d s u b ) ,
where c is the speed of light, ϕ s u b and ϕ a i r are phases of substrate and air signals, d s u b is the substrate thickness, and E ^ a i r is the complex amplitude of the THz wave transmitted through the air. The complex film permittivity was extracted from the complex conductivity as
ε ^ ( f ) = 1 + i σ ^ ( f ) / ( 2 π f d f i l m ε 0 ) ,
where d f i l m is the film thickness and ε 0 is a permittivity of free space.
The results of measurements are shown in Figure 3, Figure 4 and Figure 5. It should be mentioned that they are in good enough agreement with the classical Drude model (for conductivity σ ^ and permittivity ε ^ ) due to the semimetallic nature of the films:
σ ^ ( ω ) = σ 0 1 i ω τ ,
ε ^ ( ω ) = ε ω p 2 ω 2 + i γ ω ,
where σ 0 is the DC conductivity, ω = 2 π f is the angular frequency of THz radiation, τ is the relaxation time of charge carriers, ε is high-frequency dielectric constant, ω p is the plasma frequency, and γ is the damping rate.
Material properties used in the heat transfer interface are listed in Table 1.
The real part of relative complex permittivity is strongly negative for both the films in the studied frequency range, in accordance with the Drude model, for metals/semimetals below the plasma frequency. The behaviour of electrical conductivity also follows this model; generally, it decreases with increasing frequency of THz radiation. The electrical conductivity is higher for pure stibium film, as well as absolute values of real and imaginary parts of its permittivity. The mica has dielectric properties, so it’s permittivity does not follow the Drude theory.
In order to evaluate the dispersion of copper-relative permittivity, the Drude–Lorentz model for metals was used [32]. Calculated real (Re) and imaginary (Im) parts of relative permittivity are shown in Table 2.

3. Results

3.1. Temperature Distribution Along the Antenna Length

The distribution of temperature difference between absolute temperature and ambient temperature T a m b = 293.15 K is presented in Figure 6. Hot junction corresponds to x-coordinate at 500 μ m. The B i 88 S b 12 layer is heated up to a higher temperature in comparison with the S b layer. Temperature in the B i 88 S b 12 layer is higher than in the hot junction. The impact of heating in the C u layer is almost negligible. The results are presented for different h p o r t values. The same picture of temperature distribution along the antenna length was true for all simulation cases (with different power input and frequencies).

3.2. Temperature Difference and Responsivity

One of the main goals of this paper is to evaluate the responsivity of the detector. For this purpose, one should know the amount of power scattered and transferred to the detector working surface. The amount of power absorbed will differ from the power on the excitation port. The difference between excitation power and power on the detector surface can be expressed in terms of S-parameter (scattering parameter). The part of S-parameter expressing the transmission is S 21 . The S 21 value differs within the change of h p o r t and has a dependence on frequency. Parameter S 21 was evaluated using the following expression:
S 21 = p o r t 2 ( E c · E 2 * ) dA 2 p o r t 2 ( E 2 · E 2 * ) d A 2 ,
where E c —the computed electric field which consists of the excitation and reflected field, A 2 is the port 2 surface area, and * symbol means that the calculation is to be made at a specific mode number.
The amount of power P transmitted to the material can be found in terms of the following expression:
P = P e x · S 21 2 ,
where S 21 —transmission parameter, P e x —emitted power.
In Figure 7, the square of S 21 absolute value is presented. The results for different P e x values were equal. In our case, distances between wave excitation port and antenna from several tens of microns to several millimetres are allowable. At the distance equal to 2 mm, the part of scattered radiation is from 0.81 to 0.97. The wave with the frequency of 0.3 THz scatters more intensively than with 0.4 and 0.5 THz. The lowest scattering rate can be seen at the frequency of 0.5 THz.
The oscillations in transmission parameter dependence on the distance between THz radiation source and detector are related to wave interference between source and detector. The troughs in the transmission occur at distances given by h p o r t = ( c / ( 2 f ) ) · N , where c is the speed of light, f is the frequency, and the integer N is the order of the trough.
The transmission parameter decreases with the increase of the distance between wave excitation port and antenna. Hence, the resulting hot junction temperature T h must be inversely proportional to the h p o r t . In Figure 8, Figure 9 and Figure 10, the dependence of temperature difference Δ T h between hot junction and ambient temperature on h p o r t is presented. Hot junction temperature was calculated as the average value within the hot contact length.
The responsivity of the antenna at different frequencies is presented in Figure 11. The dependence is nonlinear for all the frequencies and changes its value from 1 to 2.5 mVW 1 . The responsivity at 0.3 THz has the lowest values, no higher than 1.5 mVW 1 . The responsivity at 0.5 THz is the highest—from 1.25 to 2.5 mVW 1 . No dependence of responsivity on the excitation power was observed.

4. Discussion

Temperature difference between hot and cold junctions in TE circuit was investigated. Voltage in TE circuit is found as Δ V = ( | S n | + | S p | ) · Δ T . S b - B i 88 S b 12 circuit can produce ( | S n | + | S p | ) of around 100 μ V/K. For an excitation port power of 500 μ W, a possible response is expected to be around 1 μ V. At lower power (e.g., P e x =1 μ W) the expected response is lower. At P e x from 1 to 50 mW, the voltage is expected to be from 10 μ V to 0.1 mV, which is possible to detect.
The calculated responsivity is from 1 to 2.5 mVW 1 . The increase in frequency leads to the increase in responsivity. Higher values of TE Seebeck coefficients of thermocouple (up to 300 V K 1 ) can increase the responsivity by 2–3 times. The dependence of responsivity on the distance between excitation port and detector has smoother behaviour when the h p o r t -value is higher than 1 mm.
Thin layer of B i 88 S b 12 showed higher sensitivity to THz radiation with higher resulting temperatures. The temperature of hot junction is lower than in the B i 88 S b 12 layer. Hence, it is possible to do a thermoelectric circuit containing one TE material— B i 88 S b 12 and metal, for example N i . Due to the rather high temperature increase (1-2 K at single B i 88 S b 12 ) and the significant Seebeck coefficient (around 100 μ VK 1 ), a voltage of around 0.1 mV can be expected in such a circuit.
Moreover, other charge-density-wave compounds with high thermoelectric properties at room temperature may be potentially used for THz-detecting purposes. Thin films of well-known thermoelectrics based on Bi, Te, Sb, Sn, and their compounds are of great interest for such applications.

Author Contributions

Conceptualization, M.K.K., N.S.K., and A.S.T.; methodology, A.V.A., A.S.T., N.S.K., and P.S.D.; formal analysis, A.S.T., A.V.A., and A.V.N.; investigation, A.S.T., P.S.D.; data curation, A.D.Z. and A.V.N.; writing—original draft preparation, A.S.T. and A.V.A.; writing—review and editing, M.K.K.; visualization, I.L.T.; supervision, M.K.K.; funding acquisition, M.K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Russian Science Foundation, grant number 19-72-10141.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
TEThermoelectric
THzTerahertz
SbAntimony
BiBismuth
NiNickel

References

  1. Zaytsev, K.I.; Kudrin, K.G.; Koroleva, S.A.; Fokina, I.N.; Volodarskaya, S.I.; Novitskaya, E.V.; Yurchenko, S.O. Medical diagnostics using terahertz pulsed spectroscopy. Phys. Conf. Ser. 2014, 486, 012014. [Google Scholar] [CrossRef] [Green Version]
  2. Yang, X.; Zhao, X.; Yang, K.; Liu, Y.; Liu, Y.; Fu, W.; Luo, Y. Biomedical applications of terahertz spectroscopy and imaging. Trends Biotechnol. 2016, 34, 810–824. [Google Scholar] [CrossRef] [PubMed]
  3. Choi, H.; Son, J.-H. Terahertz Imaging and Tomography Techniques. In Terahertz Biomedical Science and Technology; Son, J.-H., Ed.; CRC Press: Boka Raton, FL, USA, 2014; p. 47. [Google Scholar]
  4. Hoshina, H.; Sasaki, Y.; Hayashi, A.; Otani, C.; Kawase, K. Noninvasive mail inspection system with terahertz radiation. Appl. Spectrosc. 2009, 63, 81–86. [Google Scholar] [CrossRef] [PubMed]
  5. Federici, J.F.; Schulkin, B.; Huang, F.; Gary, D.; Barat, R.; Oliveira, F.; Zimdars, D. THz imaging and sensing for security applications—explosives, weapons and drugs. Semicond. Sci. Tech. 2005, 20, 266–280. [Google Scholar] [CrossRef]
  6. Malhotra, I.; Jha, K.R.; Singh, G. Terahertz antenna technology for imaging applications: A technical review. Int. J. Microw. Wirel. T. 2018, 10, 271–290. [Google Scholar] [CrossRef]
  7. Sizov, F.; Rogalski, A. THz detectors. Prog. Quant. Electron. 2010, 34, 278–347. [Google Scholar] [CrossRef]
  8. Rowe, D.M. General Principles and Basic Considerations. In Thermoelectrics Handbook; Rowe, D.M., Ed.; CRC Press: Boka Raton, Fl, USA, 2006; pp. 1-1–1-2. [Google Scholar]
  9. Yan, J.; Liao, X.; Yan, D.; Chen, Y. Review of Micro Thermoelectric Generator. J. Microelectromechanical Syst. 2018, 27, 1–18. [Google Scholar] [CrossRef]
  10. Basko, D. A photothermoelectric effect in graphene. Science 2011, 334, 610–611. [Google Scholar] [CrossRef] [Green Version]
  11. Buscema, M.; Barkelid, M.; Zwiller, V.; van der Zant, H.S.; Steele, G.A.; Castellanos-Gomez, A. Large and tunable photothermoelectric effect in single-layer MoS2. Nano Lett. 2013, 13, 358–363. [Google Scholar] [CrossRef] [Green Version]
  12. Yan, Y.; Liao, Z.M.; Ke, X.; Van Tendeloo, G.; Wang, Q.; Sun, D.; Yu, D.P. Topological surface state enhanced photothermoelectric effect in Bi2Se3 nanoribbons. Nano Lett. 2014, 14, 4389–4394. [Google Scholar] [CrossRef]
  13. Huang, D.; Zou, Y.; Jiao, F.; Zhang, F.; Zang, Y.; Di, C.A.; Zhu, D. Interface-Located Photothermoelectric Effect of Organic Thermoelectric Materials in Enabling NIR Detection. ACS Appl. Mater. Inter. 2015, 7, 8968–8973. [Google Scholar] [CrossRef] [PubMed]
  14. Cai, X.; Cai, X.; Sushkov, A.B.; Suess, R.J.; Jadidi, M.M.; Jenkins, G.S.; Nyakiti, L.O.; Murphy, T.E. Sensitive room-temperature terahertz detection via the photothermoelectric effect in graphene. Nat. Nanotechnol. 2014, 9, 814–819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. He, X.; Fujimura, N.; Lloyd, J.M.; Erickson, K.J.; Talin, A.A.; Zhang, Q.; Léonard, F. Carbon nanotube terahertz detector. Nano Lett. 2014, 14, 3953–3958. [Google Scholar] [CrossRef] [PubMed]
  16. Mao, J.; Liu, Z.; Ren, Z. Size effect in thermoelectric materials. NPJ Quantum Mater. 2016, 1, 1–9. [Google Scholar] [CrossRef]
  17. Mele, P.; Narducci, D.; Ohta, M.; Biswas, K.; Morante, J.; Saini, S.; Endo, T. Masayuki Takashiri Thin Films of Bismuth-Telluride-Based Alloys. In Thermoelectric thin Films; Springer Nature Switzerland: Cham, Switzerland, 2019; pp. 1–29. [Google Scholar]
  18. Volklein, F.; Kessler, E. Thermoelectric properties of Bi1–xSbx films with 0< x < 0.3. Thin Solid Film. 1990, 155, 197–208. [Google Scholar]
  19. Suslov, M.V.; Komarov, V.A.; Suslov, A.V. Thermoelectric power of Bi92 Sb8 and Bi85 Sb15 thin films. Semiconductors 2017, 51, 862–865. [Google Scholar] [CrossRef]
  20. Cho, S.; DiVenere, A.; Wong, K.G.; Ketterson, B.J. Thermoelectric transport properties of n-doped and p-doped Bi0.91Sb0.09 alloy thin films. J. Appl. Phys. 1999, 85, 3655–3660. [Google Scholar] [CrossRef]
  21. Linseis, V.; Volklein, F.; Reith, H.; Huhne, R.; Schnatmann, L.; Nielsch, K.; Woias, P. Thickness and temperature dependent thermoelectric properties of Bi87Sb13 nanofilms measured with a novel measurement platform. Semicond. Sci. Technol. 2018, 33, 085014. [Google Scholar] [CrossRef]
  22. Völklein, F.; Kessler, E. Thermal Conductivity and Thermoelectric Figure of Merit of Thin Antimony Films. Phys. Status Solidi B 1990, 158, 521–529. [Google Scholar] [CrossRef]
  23. Zemskov, V.S.; Belaya, A.D.; Beluy, U.S.; Kozhemyakin, G.N. Growth and investigation of thermoelectric propertiesof Bi-Sb alloy single crystals. J. Cryst. Growth 2000, 212, 161–166. [Google Scholar] [CrossRef]
  24. Brandt, N.B.; Svistova, E.A. Uspekhi Phizich. Nauk 1970, 101, 249–272. (In Russian) [Google Scholar]
  25. Benia, H.M.; Straser, C.; Kern, K.; Ast, C.R. Surface band structure of Bi(1-x)Sbx (111). Phys. Rev. B 2015, 91, 161406. [Google Scholar] [CrossRef] [Green Version]
  26. Ihring, A.; Kessler, E.; Dillner, U.; Haenschke, F.; Schinkel, U.; Schubert, M.; Meyer, H.G. High performance uncooled THz sensing structures based on antenna-coupled air-bridges. Microelectron. Eng. 2012, 98, 512–515. [Google Scholar] [CrossRef]
  27. Huhn, A.K.; Spickermann, G.; Ihring, A.; Schinkel, U.; Meyer, H.G.; Haring Bolivar, P. Uncooled antenna-coupled terahertz detectors with 22 μs response time based on BiSb/Sb thermocouples. Appl. Phys. Lett. 2013, 102, 121102. [Google Scholar] [CrossRef] [Green Version]
  28. Mastel, S.; Lundeberg, B.M.; Alonso-Gonzalez, P.; Yuanda, G.; Watanabe, K.; Taniguchi, T.; Hone, J.; Koppens, H.L.F.; Nikitin, A. Terahertz Nanofocusing with Cantilevered Terahertz-Resonant Antenna Tips. Nano Lett. 2017, 17, 6526–6533. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Xu, D.; Wang, Y.; Xiong, B.; Li, T. MEMS-based thermoelectric infrared sensors: A review. Front. Mech. Eng. 2017, 12, 557–566. [Google Scholar] [CrossRef]
  30. Grabov, V.M.; Komarov, V.A.; Kablukova, N.S. Galvanomagnetic Properties of Thin Films of Bismuth and Bismuth–Antimony Alloys on Substrates with Different Thermal Expansions. Phys. Solid State 2016, 58, 605–611. [Google Scholar] [CrossRef]
  31. Tinkham, M. Energy gap interpretation of experiments on infrared transmission through superconducting films. Phys. Rev. 1956, 104, 245–246. [Google Scholar] [CrossRef]
  32. Stefan, A. Maier Plasmonics: Fundamential and Applications; Springer: Berlin, Germany, 2007; p. 223. [Google Scholar]
Figure 1. 3D view of simulated antenna. 1—mica substrate; 2,5—copper films, 3— B i 88 S b 12 film, 4— S b film, 6—vacuum chamber.
Figure 1. 3D view of simulated antenna. 1—mica substrate; 2,5—copper films, 3— B i 88 S b 12 film, 4— S b film, 6—vacuum chamber.
Applsci 10 01929 g001
Figure 2. Antenna cross-section in z-x plane; l p o r t = 1000 μ m—the length of the wave excitation port; l T E f i l m = 800 μ m–the length of thermoelectric couple containing two films, l a n t = l p o r t —the full antenna length; h f i l m = 150 nm—the thickness of the films; h m i c a = 40 μ m—the thickness of mica substrate; h p o r t —the distance between wave excitation and antenna, the value of h p o r t was changing within the simulations.
Figure 2. Antenna cross-section in z-x plane; l p o r t = 1000 μ m—the length of the wave excitation port; l T E f i l m = 800 μ m–the length of thermoelectric couple containing two films, l a n t = l p o r t —the full antenna length; h f i l m = 150 nm—the thickness of the films; h m i c a = 40 μ m—the thickness of mica substrate; h p o r t —the distance between wave excitation and antenna, the value of h p o r t was changing within the simulations.
Applsci 10 01929 g002
Figure 3. The real and imaginary parts of a relative permittivity of thin-film Sb and Bi 88 Sb 12 on mica substrate in the frequency range of 0.3–0.8 THz.
Figure 3. The real and imaginary parts of a relative permittivity of thin-film Sb and Bi 88 Sb 12 on mica substrate in the frequency range of 0.3–0.8 THz.
Applsci 10 01929 g003
Figure 4. Electrical conductivity of antimony and bismuth antimonide 150-nm thin films on mica substrate in the frequency range of 0.3–0.8 THz.
Figure 4. Electrical conductivity of antimony and bismuth antimonide 150-nm thin films on mica substrate in the frequency range of 0.3–0.8 THz.
Applsci 10 01929 g004
Figure 5. The properties of mica substrate used in the simulation in the frequency range of 0.3–0.8 THz.
Figure 5. The properties of mica substrate used in the simulation in the frequency range of 0.3–0.8 THz.
Applsci 10 01929 g005
Figure 6. Example of temperature distribution along the x-axis (antenna length). Copper films correspond to the x-coordinates: 0–100 μ m, B i 88 S b 12 to 100–500 μ m, S b to 500–800 μ m; hot junction corresponds to the x-coordinate at 500 μ m.
Figure 6. Example of temperature distribution along the x-axis (antenna length). Copper films correspond to the x-coordinates: 0–100 μ m, B i 88 S b 12 to 100–500 μ m, S b to 500–800 μ m; hot junction corresponds to the x-coordinate at 500 μ m.
Applsci 10 01929 g006
Figure 7. Transmission parameter dependence on the distance between THz radiation source and detector. The results are obtained for the power input of 50 mW.
Figure 7. Transmission parameter dependence on the distance between THz radiation source and detector. The results are obtained for the power input of 50 mW.
Applsci 10 01929 g007
Figure 8. Temperature difference on h p o r t at 0.3 THz; P e x values are listed in the legend.
Figure 8. Temperature difference on h p o r t at 0.3 THz; P e x values are listed in the legend.
Applsci 10 01929 g008
Figure 9. Temperature difference on h p o r t at 0.4 THz; P e x values are listed in the legend.
Figure 9. Temperature difference on h p o r t at 0.4 THz; P e x values are listed in the legend.
Applsci 10 01929 g009
Figure 10. Temperature difference on h p o r t at 0.5 THz; P e x values are listed in the legend.
Figure 10. Temperature difference on h p o r t at 0.5 THz; P e x values are listed in the legend.
Applsci 10 01929 g010
Figure 11. Responsivity of antenna on h p o r t in frequency range from 0.3–0.5 THz.
Figure 11. Responsivity of antenna on h p o r t in frequency range from 0.3–0.5 THz.
Applsci 10 01929 g011
Table 1. Thermophysical properties of the materials.
Table 1. Thermophysical properties of the materials.
MaterialThermal Conductivity, in-Plane,Thermal Conductivity, Cross-Plane,Heat Capacity,Density,
W·m 1 · K 1 W·m 1 · K 1 J·kg 1 · K 1 kg·m 3
Sb63.2424.432076691
B i 88 S b 12 13.215.11249790
mica5.10.518802900
Cu4014013848960
Table 2. Calculated relative permittivity of copper.
Table 2. Calculated relative permittivity of copper.
Frequency, THzReIm
0.1−6219.672.70 × 10 2
0.2−6219.421.35 × 10 2
0.3−6218.999.00 × 10
0.4−6218.396.75 × 10
0.5−6217.625.40 × 10
0.6−6216.684.50 × 10
0.7−6215.583.86 × 10
0.8−6214.33.37 × 10

Share and Cite

MDPI and ACS Style

Tukmakova, A.S.; Asach, A.V.; Novotelnova, A.V.; Tkhorzhevskiy, I.L.; Kablukova, N.S.; Demchenko, P.S.; Zaitsev, A.D.; Khodzitsky, M.K. FEM Simulation of THz Detector Based on Sb and Bi88Sb12 Thermoelectric Thin Films. Appl. Sci. 2020, 10, 1929. https://doi.org/10.3390/app10061929

AMA Style

Tukmakova AS, Asach AV, Novotelnova AV, Tkhorzhevskiy IL, Kablukova NS, Demchenko PS, Zaitsev AD, Khodzitsky MK. FEM Simulation of THz Detector Based on Sb and Bi88Sb12 Thermoelectric Thin Films. Applied Sciences. 2020; 10(6):1929. https://doi.org/10.3390/app10061929

Chicago/Turabian Style

Tukmakova, Anastasiia S., Alexei V. Asach, Anna V. Novotelnova, Ivan L. Tkhorzhevskiy, Natallya S. Kablukova, Petr S. Demchenko, Anton D. Zaitsev, and Mikhail K. Khodzitsky. 2020. "FEM Simulation of THz Detector Based on Sb and Bi88Sb12 Thermoelectric Thin Films" Applied Sciences 10, no. 6: 1929. https://doi.org/10.3390/app10061929

APA Style

Tukmakova, A. S., Asach, A. V., Novotelnova, A. V., Tkhorzhevskiy, I. L., Kablukova, N. S., Demchenko, P. S., Zaitsev, A. D., & Khodzitsky, M. K. (2020). FEM Simulation of THz Detector Based on Sb and Bi88Sb12 Thermoelectric Thin Films. Applied Sciences, 10(6), 1929. https://doi.org/10.3390/app10061929

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop