**Sung-Min Hong***

Gwangju Institute of Science and Technology 123 Cheomdan-gwagiro (Oryong-dong), Buk-gu, Gwangju 61005, Republic of KoreaNano-Micro Conference, 2017, 1, 010002

Published Online: 15 September 2017 (Article)

DOI:10.11605/cp.nmc2017.01002

*Corresponding author. E-mail: This email address is being protected from spambots. You need JavaScript enabled to view it.

**Abstract**

The transient simulation of semiconductor devices using a deterministic Boltzmann equations solver is presented. In order to avoid the numerical difficulties originated from the conventional H-transformation, the kinetic-energy-based scheme is adopted. Within the kinetic-energy-based scheme, the transport parameters become time-independent, therefore, implementing the transient simulation capability can be done easily. Preliminary results for a device are shown.

**Introduction**

The terahertz (THz) technology has gained much research interest worldwide. However, realization of a THz emitter with a small form factor is still a difficult task. The plasma instability in the two-dimensional (2D) electron gas has been proposed as a possible mechanism to enable an extremely fast semiconductor device operating in the THz frequency range [1].

Recently, we have reported the numerical simulation results of plasma oscillation in two-dimensional electron gas [2]. A periodic steady-state solver has been used to obtain the oscillating solution efficiently. The set of governing equations contains the 2D Poisson equation, the electron continuity equation, and the electron current density equation. The time derivative of the current density and the convective derivative are considered in the electron current density equation. It has been numerically verified that inclusion of those terms are mandatory to support sufficiently strong plasma instability in the 2D electron gas.However, it is noted that the set of governing equations adopted in [2] is based on the moments of the Boltzmann equation. It is not sufficiently accurate to describe the behaviour of plasma waves quantitatively [3]. A more advanced modelling approach is required. In order to solve this issue ultimately, the Boltzmann equation solver with the transient simulation capability is mandatory.

Although the stochastic Boltzmann equation solvers (the Monte Carlo simulators) have been used to simulate a device operating in the THz frequency range (for example, [4]), the stochastic nature of the Monte Carlo simulators introduces some numerical difficulties. Therefore, the deterministic Boltzmann equation solver is desirable. In this work, our motivation is to present some transient simulation results using a deterministic Boltzmann equation solver.

**Projected Boltzmann Equation**

The Boltzmann equation is a microscopic transport equation. It governs the spatiotemporal motion of the electron gas in a semiconductor device. Since it is defined on the phase space, evaluating its solution is a highly non-trivial task.

As discussed above, the stochastic Monte Carlo method can be applied to solve the Boltzmann equation efficiently. However, nowadays, the deterministic Boltzmann equation solvers [5] have gained popularity. In these solvers, the angle-dependent part of the electron distribution function is expanded with basis functions. Depending on the dimensionality of the momentum space, there can be a few variants of deterministic Boltzmann equation solvers.In this work, the three-dimensional momentum space is understood. Therefore, the spherical harmonic expansion [6] is the appropriate technique. A numerical stabilization scheme called the H-transformation [6] has played an important role in the development of this simulation technique. In the H-transformation, the total energy (which is a sum of the kinetic energy and the potential energy) instead of the kinetic energy is used as the energy variable

Unfortunately, when the transient simulation is concerned, the H-transformation introduces a serious difficulty [7]. The time-varying potential energy yields the time-varying transport parameters. Unavoidable interpolation procedure prevents the transient simulation capability eventually. Of course, the fundamental solution to the above problem is to invent yet another numerical stabilization scheme, which does not suffer from the time-varying transport parameters. Since it is currently a formidable task, in this work, we simply apply the kinetic-energy-based scheme [8] to get the preliminary solutions.

**Numerical Results **

A spherical harmonics expansion solver has been newly implemented into our in-house device simulation framework [2]. All parameters for Si are taken from [6]. The lowest expansion order is used. Various N^{+}NN^{+} devices are simulated. Fig. 1(a) shows the DC IV characteristics of those devices. The transient simulation results of the 1200-nm-long device are shown in Fig. 1(b). The DC bias voltage is 0.8 V. The frequency of the AC voltage is 1 THz. It is clearly demonstrated that our deterministic Boltzmann equation solver can perform the transient simulation at such a high frequency without any numerical problem.

**Conclusion**

In this work, the transient device simulation results using a deterministic Boltzmann equation solver based on the spherical harmonics expansion has been presented. In order to avoid the numerical difficulties originated from the H-transformation, the kinetic-energy-based scheme is adopted. The numerical results demonstrate that both of DC and AC analyses can be performed without any numerical problem. It is expected that this simulator can be applied to the THz devices based on the plasma instability.

**Acknowledgements**

This work was supported by Samsung Research Funding Center of Samsung Electronics under Project Number SRFC-IT1401-08.

**References**

[1] Taiichi Otsuji, Takayuki Watanabe, Stephane A. Boubanga Tombet, Akira Satou, Wojciech M. Knap, Vyacheslav V. Popov, Maxim Ryzhii, Victor Ryzhii. IEEE Transactions on Terahertz Science and Technology, 2013,3(1), 63-7. doi: 10.1109/TTHZ.2012.2235911

[2] Sung-Min Hong, Jae-Hyung Jang. IEEE Transactions on Electron Devices. 2015, 62(12), 4192-4198. doi: 10.1109/TED.2015.2489220

[3] Zeinab Kargar, Tobias Linn, Dino Ruić, Christoph Jungemann. IEEE Transactions on Electron Devices. 2016, 63(11),4402-4408. doi: 10.1109/TED.2016.2608422

[4] Javier Mateos, Toms Gonzalez. IEEE Transactions on Terahertz Science and Technology. 2012, 2(5), 562-569. doi: 10.1109/TTHZ.2012.2209970

[5] Sung-Min Hong, Anh-Tuan Pham, Christoph Jungemann. Deterministic Solvers for the Boltzmann Transport Equation, ISBN 978-3-7091-0777-5, Spring Verlag Wien/New York, 2011.

[6] A.Gnudi, D. Ventura, G. Baccarani, F. Odeh. Solid-State Electronics, 1993, 36, (4), 575-581. doi: 10.1016/0038-1101(93)90269-V

[7] K. Rupp, C. Jungemann, S.-M. Hong, M. Bina, T. Grasser, A. Jüngel. Journal of Computational Electronics. 2016, 15(3), 939-958. doi: 10.1007/s10825-016-0828-z

[8] C. Jungemann, A. T. Pham, B. Meinerzhagen, C. Ringhofer, M. Bollhöfer. Journal of Applied Physics, 2006, 100, 024502. doi: 10.1063/1.2212207

**Open Access**

This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

© The Author(s) 2017