太赫兹大孔径折射透镜的焦移特性研究

李惟帆; 祁峰; 汪业龙; 刘鹏翔; 刘朝阳; LI Wei-Fan; QI Feng; WANG Ye-Long; LIU Peng-Xiang; LIU Zhao-Yang

网刊加载中。。。

使用Chrome浏览器效果最佳，继续浏览，你可能不会看到最佳的展示效果，

确定继续浏览么?

复制成功，请在其他浏览器进行阅读

Focal shift effect of terahertz wide-aperture refractive lens PDF

- ORCID：
LI Wei-Fan ^1,2,3
✉
- ORCID：
QI Feng ^1,2,3
✉
- ORCID：
WANG Ye-Long ^1,2,3
- ORCID：
LIU Peng-Xiang ^1,2,3
- ORCID：
LIU Zhao-Yang ^1,2,3

1. Key Laboratory of Opto-Electronic Information Processing， Chinese Academy of Sciences，Shenyang 110169， China； 2. Key Laboratory in Terahertz Imaging and Sensing， Liaoning Province， Shenyang 110169， China； 3. Shenyang Institute of Automation， Chinese Academy of Sciences， Shenyang 110169， China

CLC： O43

Updated：2023-06-18

DOI：10.11972/j.issn.1001-9014.2023.03.008

Abstract

The focal shift effect of terahertz （THz） beam focusing when using a wide-aperture refractive lens has been investigated. The deviation of focus position caused by the focal shift effect can adversely affect the imaging or measurement quality of a THz system. In this study， reference values of relative focal shifts and combinations of different lens apertures， focal lengths， and working frequencies were analyzed and discussed through theoretical calculation and finite element analysis simulation. When using the commercial lens， the actual focus was determined based on the focal shift effect to ensure the working efficiency of a terahertz system. Concerning the customized lens design， the focal shift distance was compensated in the focal length optimization according to the working frequency. These two approaches guaranteed the good performance of a THz system.

Keywords

terahertz lens; focal shift; finite element analysis; aspherics

Introduction

The performance of a terahertz （THz） imaging or spectroscopy system is closely related to the quality of the THz beam^［

1-2］. Due to the diffractive divergence caused by the sub-millimeter scale wavelength， THz wave usually needs reshaping whether it is generated from electrical equipment or an optical system^{［Reference 3-4}3-4］. THz wide-aperture refractive lens （WARL） is a commonly used THz beam-shaping device with the advantages of broadband applicability， simple structure， high stability， and low cost. Sub-wavelength scale focusing spot and high spatial resolution THz image can be obtained using THz WARL with a high diameter to focal length ratio^{［Reference 5

Baidu Scholar}5］. Most THz systems require an accurate focus position rather than a small spot because the focus is where test samples are usually placed and determines where to image or measure. One of the factors influencing the accuracy of the focus position is spherical aberration. In our previous work， the necessity of eliminating spherical aberration in the THz WARL design was discussed^{［Reference 6

Baidu Scholar}6］. Another crucial and frequently overlooked factor is the focal shift effect， a phenomenon first known in Gauss optics. When a Gaussian beam is focused by a lens， with the beam waist in the plane of the lens and much smaller than the lens aperture， the maximum intensity is located at the waist of the focused beam that is closer to the lens rather than the geometric focus. The focal shift effect can be generalized to uniform pupil illumination in aberration-free systems as the incident beam waist becomes much larger than the lens diameter^{［Reference 7-8}7-8］. The focal shift effect can be explained by the diffraction theory and is related to the Fresnel number of a THz beam shaping device. For the same THz WARL， the lower the working frequency， the more severe the adverse effect of focal shift on a THz system^{［Reference 9

Baidu Scholar}9］. Therefore， the focal shift effect should be fully considered and compensated with the design and use of THz WARLs. In this paper， the focal shift effect of the THz beam focused by THz WARLs was theoretically studied. Finally， the design and use of THz WARLs under different research or application requirements were proposed after the analysis of different apertures， focal lengths， and working wavelengths.

1 Theoretical calculation of focal shift effect

The focal shift theory was first derived from the Huygens-Fresnel principle by Li and Wolf^［

8］. Fig. 1 illustrates the schematic of THz beam focusing and focal shift effect. The field at any point P along the optical axis that is far away from the diffraction aperture^{［Reference 10

Baidu Scholar}10］ is expressed as：

U (P) = - \frac{i k}{2 π} \frac{A e x p (- i k f)}{f} \iint_{S}^{} \frac{e x p (i k s)}{s} d S

（1）

Fig. 1 The schematic of THz beam focusing and focal shift effect

图1 太赫兹波束聚焦及焦移效应示意图

Eq. （1） can be simplified as：

U [P (z)] = A \frac{e x p (i k z)}{z} [e x p (- \frac{1}{2} i k \frac{z}{f} \frac{a^{2}}{f + z}) - 1]

（2）

where a and f denote the semi-aperture and focal length of THz WARL， respectively. The intensity of point P is：

I (P) = U (P) U^{*} (P) = I_{0} {(1 - \frac{u_{N}}{2 π N})}^{2} {[\frac{s i n (u_{N} / 4)}{u_{N} / 4}]}^{2}

（3a）

u_{N} = \frac{k a^{2} z}{z (f + z)} = 2 π N \frac{z}{f + z}

，

（3b）

where $I_{0} = {(π a^{2} |A| / λ f^{2})}^{2}$ indicates the intensity of geometrical focus， and $N = a^{2} / λ f$ represents the Fresnel number.

Fig. 2（a） exhibits the intensity distributions along the optical axis given by Eq. （3）. The $z = 0$ on the horizontal axis is the position of geometrical focus， and the intensity is set to unity. As N decreases， the peak intensity increases， and the focus position moves towards the negative direction of the optical axis where the lens is located. Fig. 2（b） demonstrates the relative focal shift （ $Δ f / f$ ） and relative peak intensity （ ${I / I}_{0}$ ） with N. As N increases， $Δ f / f$ increases to zero， and ${I / I}_{0}$ approaches unity. $Δ f / f$ is -1%， -10%， and -40% when $N = 11$ ， $N = 3$ ， and $N = 1$ ， respectively. The values of the focal shift are all negative， reflecting that the focus always moves toward the lens.

Fig. 2 Theoretical calculations of focal shift effect （a） intensity distributions along the optical axis at different Fresnel numbers，（b） relative focal shift and peak intensity with Fresnel number

图2 焦移效应的理论计算（a）不同菲涅尔数下沿光轴方向的强度分布曲线，（b）相对焦移和峰值强度与菲涅尔数的关系曲线

In the simplification of Eq. （2）， the conditions of $a ≫ λ$ and ${(a / f)}^{2} ≪ 1$ are used to neglect the high-order terms of the expansion. Regarding a THz WARL， the above two inequations cannot be satisfied in many cases. As the focal length decreases， the center thickness of the lens increases， and the thin lens approximation no longer holds. All of these factors may induce deviations between practical and theoretical focal shift values.

2 Simulation of focal shift effect of THz WARL

Electromagnetic simulation based on finite element analysis （FEA） was performed to analyze the intensity distribution of THz beam focusing， so as to explore the focal shift effect of THz WARLs. Considering broadband THz research and applications， the focusing characteristics of the same THz WARL at different working frequencies were first studied. The incident wave was chosen to be a Gaussian plane wave. The boundary was set to be an absorbing condition to simulate the actual infinite space. The specific methods for THz WARL design and optimization were analyzed in our previous work^［

6］. The refractive index of THz WARL material was set to 1.41， and the absorption factor was set to 0 since the weak absorption characteristic of polymer material was below 1 THz. The aperture and focal length of the selected THz WARL were both 20 mm， and the working frequencies were 0.1-0.6 THz with an interval of 0.1 THz. The collimated THz beam passed through the plano surface and then was refracted and focused by the aspherical convex surface， as illustrated in Fig. 1. Theoretically， the plano-convex surfaces can eliminate spherical aberration. Figures 3（a） to 3（c） present the amplitude distributions obtained by FEA simulations at working frequencies of 0.1， 0.3， and 0.5 THz， respectively. In Fig. 3（a）， the white solid line and dashed line indicate the positions of the physical focal plane and geometrical focal plane， respectively. Fig. 3（d） depicts the intensity distributions along the optical axis. The peak of each curve represents the focus and is marked by a solid circle. Fig. 3（e） exhibits

Δ f / f

and

{I / I}_{0}

with the working frequency and N.

Δ f / f

is -11.1% at 0.2 THz and -1.4% at 0.5 THz. These figures verify that the focal shift effect is weaker at higher working frequencies. The differences between the simulation and calculation results in Fig. 3（e） are reasonable since the fixed lens is thick and does not meet

{(a / f)}^{2} ≪ 1

and partly meets

a ≫ λ

Fig. 3 Focal shift analysis of the THz WARL with a 20 mm aperture and focal length （a） to （c） FEA simulations at 0.1， 0.3， and 0.5 THz， respectively，（d） intensity distributions along the optical axis，（e） relative focal shift and peak intensity with working frequency and Fresnel number

图3 孔径和焦距均为20 mm的太赫兹WARL的焦移特性分析（a）至（c）0.1、0.3和0.5 THz下的FEA仿真结果，（d）沿光轴方向的强度分布曲线，（e）相对焦移和峰值强度与工作频率和菲涅尔数的关系

Similarly， FEA simulations were conducted using different apertures （2a）， focal lengths （f）， and working frequencies （ν）， with N values ranging from 2 to 32. The results are rendered in Fig. 4， which are consistent with the theoretical calculation. The error was derived from two aspects： the neglect of high-order terms in the derivation of the theoretical equation； the invalidation of thin lens approximation. Nevertheless， the theoretical curve can still provide a good reference for the design and use of a THz WARL.

Fig. 4 FEA simulation results obtained using different lens apertures， focal lengths， and working frequencies

图4 不同透镜孔径、焦距和工作频率下的FEA仿真结果

3 Minimize focal shift effect in design and use

According to the FEA simulations， the focus position of a THz WARL is closer to the theoretical position when N is large. With respect to a customized THz WARL， the aperture and focal length can be freely designed following the working frequency and desired N value. Although N is relatively small， the focal shift can be compensated by optimizing the surface parameters. For example， the THz WARL with a 20 mm aperture and the focal length was discussed in section 3. When the working frequency is 0.1 THz， the THz WARL has an N of 1.67 and $Δ f / f$ of -23.4%， and the physical focal length is calculated to be 15.32 mm. If the physical focal length is set to 20 mm， the geometrical focal length is calculated to be 28.75 mm. When the working frequency is 0.1 THz， the physical focal length is 14.7 mm， and the focal shift value is 5.3 mm （ $Δ f / f$ = -26.5%）， as demonstrated in Fig. 5（a）. Fig. 5（b） renders the simulation with focal shift compensation. The central thickness of the THz WARL remained the same， and the curvature of the convex surface was modified based on the theoretical calculation in Fig. 2（b）. The optimized physical focal length is 19.0 mm， and the focal shift value is 1.0 mm （ $Δ f / f$ = -5.0%）， which is an improvement of more than 20%.

Fig. 5 FEA simulations at 0.1 THz of the THz WARL with a 20 mm aperture and focal length，（a） without and （b） with focal shift compensation

图5 孔径和焦距均为20 mm的太赫兹WARL在0.1 THz频率下的FEA仿真，（a）有焦移补偿，（b）无焦移补偿

However， the choice of aperture and focal length is relatively limited for commercial THz WARLs. The commercial THz WARLs with plano-convex surfaces of three companies were investigated. The parameters and analyses of some commercial lenses are listed in Table 1. The focal length is greater than the aperture in most cases. Particularly， the real lens aperture should be smaller than the lens diameter since the THz WARL needs to be clamped when in use. Additionally， 90% of the diameter is taken as the aperture （2a） in the last column of Table 1.

Table 1 Commercial THz WARLs with plano-convex surfaces

表1 具有平凸表面的商用太赫兹WARL

Brand	Surface type	Diameter （mm）	Focal length （mm）	$a^{2} / f$ （mm）
Batop	Aspherical	25.4	10-67	1.9-13.0
	Aspherical	50.8	35-100	5.2-14.9
	Spherical	25.4	100-150	0.9-1.3
	Spherical	50.8	200-250	2.1-2.6
Thorlabs	Spherical	50.8	75-500	1.0-7.0
		76.2	115-150	7.8-10.2
		101.6	151.5-200	10.5-13.8
Tydex	Spherical	25.4	50-250	0.5-2.6
		50.8	75-2000	0.3-7.0
		76.2	100-450	2.6-11.8
		101.6	200-5 000	0.4-10.5

Fig. 6 presents the focal shift （ $Δ f / f$ ） in relation to lens parameters （ $a^{2} / f$ ） and working frequency. The area corresponding to the situation where commercial THz WARLs could be used is marked with a light-yellow background. With the curve of $Δ f / f = 1 %$ as an example， any combination of the $a^{2} / f$ and working frequency above this curve satisfies $Δ f / f < 1 %$ . Considering that the convex surface of most lenses in Table 1 is spherical， the actual shift of focus position is larger due to spherical aberration. The spherical aberration causes the focal plane to form a diffuse spot， which is also approximated by a shift in the focus position. Thus， the actual focus position should be determined according to the focal shift distance and spherical aberration when using commercial THz WARLs to weaken the adverse effect of focal shift， such as reducing the accuracy of imaging or measurement. Regarding customized THz WARLs， the focal length should be increased following the working frequency in the design to compensate for the focal shift effect， if a specific aperture and focus position are required.

Fig. 6 Focal shift （ $Δ f / f$ ） in relation to lens parameters （ $a^{2} / f$ ） and working frequency， the area corresponding to the situation where commercial THz WARLs could be used is marked with a light-yellow background

图6 焦移（Δf/f）与透镜参数（a2/f）和工作频率的关系，浅黄色背景标记了使用商用太赫兹WARL时涉及的参数所对应的区域

4 Conclusions

In conclusion， the focal shift effect of THz beam focusing when using THz WARLs was studied. As revealed by the calculation of the theoretical expression， the actual focus position always shifted towards the lens direction. The relative focal shift is inversely proportional to the Fresnel number of a THz WARL. FEA simulations of different lens apertures， focal lengths， and working frequencies were conducted. The focal shift values were consistent with the theoretical calculation. The inaccurate focus position will impact the working efficiency of the components of a THz system， leading to curtailed imaging or measurement quality. Hence， the focal shift effect should be minimized in the practical use of THz WARLs. Concerning commercial THz WARLs， especially the configurations of focal length greater than the aperture， the actual focus position should be determined under the focal shift effect. For customized THz WARLs， if there are specific requirements for the lens size and working frequency， the focal shift distance can be offset by increasing the theoretical focal length in the design. These two approaches can guarantee the good performance of a THz system， such as the high spatial resolution of a THz image and large dynamic range of the THz spectrum.

References

Mittleman D M. Twenty years of terahertz imaging ［Invited］［J］. Optics Express， 2018， 26（8）：9417–9431. 10.1364/oe.26.009417 [Baidu Scholar]

Neu J， Schmuttenmaer C A. Tutorial： An introduction to terahertz time domain spectroscopy （THz-TDS）［J］. Journal of Applied Physics， 2018， 124（23）：231101. 10.1063/1.5047659 [Baidu Scholar]

Ferguson B， Zhang X C. Materials for terahertz science and technology［J］. Nature Materials， 2002， 1（1）：26–33. 10.1038/nmat708 [Baidu Scholar]

Liu Z， Qi F， Wang Y， et al. A 220- to 299-GHz CMOS terahertz detector［J］. Journal of Infrared Millimeter and Terahertz Waves， 2019， 40（6）：606–619. 10.1007/s10762-019-00592-2 [Baidu Scholar]

Lo Y H， Leonhardt R. Aspheric lenses for terahertz imaging［J］. Optics Express， 2008， 16（20）：15991–15998. 10.1364/oe.16.015991 [Baidu Scholar]

Li W， Qi F， Wang Y， et al. Refractive aspherical lens for terahertz imaging［J］. Optics Communications， 2018， 433：14-17. 10.1016/j.optcom.2018.09.071 [Baidu Scholar]

Li Y. Focal shift formula for focused， apertured Gaussian beams［J］. Journal of Modern Optics， 1992， 39（8）： 1761–1764. 10.1080/713823585 [Baidu Scholar]

Li Y， Wolf E. Focal shifts in diffracted converging spherical waves［J］. Optics Communications， 1981， 39（4）：211-215. 10.1016/0030-4018(81)90108-5 [Baidu Scholar]

Ferlet M. Imaging at low Fresnel number： some challenges and applications［J］. Proceedings of SPIE - The International Society for Optical Engineering， 2006， 6289. 10.1117/12.680509 [Baidu Scholar]

Born M， Wolf E. Principles of optics： electromagnetic theory of propagation， interference and diffraction of light， seventh edition［M］. Cambridge： Cambridge University Press， 1999.杨葭荪译，光学原理：光的传播、干涉和衍射的电磁理论，第7版，北京，电子工业出版社， 2016. [Baidu Scholar]

You are the first Visitors

Governing Body：中国科学院

Organizer：中国科学院上海技术物理研究所，中国光学学会

Address：上海市玉田路500号 Phone：021-25051553

51La

Home

Introduction

Editorial Board

Authors' Guidelines

Open Access

Publishing Ethics

Downloads

Copyright

Contact Us

中文版

Focal shift effect of terahertz wide-aperture refractive lens PDF

Abstract

Keywords

Introduction

1 Theoretical calculation of focal shift effect

2 Simulation of focal shift effect of THz WARL

3 Minimize focal shift effect in design and use

4 Conclusions

References