Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Journal of Applied Physics

An exact tunneling model and its application to transmission and reflection delay times

Kevin L. Jensen, Jeanne M. Riga, Joel L. Lebowitz, Rebecca Seviour, Donald A. Shiffler

Published Sep 27, 2022• 7 citations
Publisher

Abstract

Electron emission into a nanogap is not instantaneous, which presents a difficulty in simulating ultra-fast behavior using particle models. A method of approximating the transmission and reflection delay (TARD) times of a wave packet interacting with barriers described by a delta function, a metal–insulator–metal (MIM, rectangular) barrier, and a Fowler Nordheim (FN, triangular) barrier is given and has application to simulation. It is based on the superposition of a finite number of exact basis states obtained from Schrödinger’s equation, analogous to how quantum carpets are simulated. As a result, it can exactly and uniquely follow exponentially small tunneling currents. A Bohm-like trajectory is obtained from the time evolution of the density: it shows delay in both the transmitted and reflected packets that can be simply evaluated. The relations to prior studies of the analytic [Formula: see text]-function barrier and the Wigner distribution function (WDF) methods are described. A comparison of the TARD times is contrasted to alternate times in the Büttiker–Landauer (BL) and McColl–Hartman (MH) times; the MH approach is further reformulated explicitly in terms of Gamow factors to consider how the McColl–Hartman effect is to be related, particularly in the case of the FN barrier of field emission.

DOI
10.1063/5.0096568

Related papers

Wigner wave packets: Transmission, reflection, and tunneling

Kevin L. Jensen, Joel L. Lebowitz, Jeanne M. Riga +2 more

A numerical solution to the time evolution equation of the Wigner distribution function (WDF) with an accuracy necessary to simulate the passage of a wave packet past a barrier is developed, where quantum effects require high accuracy and fine discretization. A wave packet incident on a barrier, a portion of which tunnels through, demonstrates behavior that can define various characteristic transmission and reflection delay (TARD) times useful in the simulation of electron emission. A model for the TARD times is proposed that relies only on the asymptotic maxima of the position $\ensuremath{\rho}(x,t)$ and wave number $\ensuremath{\rho}(k,t)$ densities given by the WDF and applied to a ballistic trajectory model for the (faster) transmitted and (slower) reflected parts. The dependence of the TARD times on barrier width, symmetry, and abruptness is analyzed. For symmetrical barriers with characteristics similar to field emission barrier heights and widths, TARD times are on the order of a fraction of a femtosecond. The TARD times for when tunneling predominates are contrasted to tunneling times in the literature. Use of the TARD times in simulations of field emission in nanogap devices or to model ultrashort pulses generated under rapidly changing conditions for electron sources are proposed.

Physical review. B./Physical review. B•Apr 26, 2021•15 citations

Reevaluating the Hartman effect for field emission

Kevin L. Jensen, Andrew Shabaev, Jeanne M. Riga +3 more

An electron wave packet tunneling through a barrier has a transmission (or ``group delay'') time ${\ensuremath{\tau}}_{g}$ that, for a rectangular barrier, is commonly held to become independent of the barrier width $L$ as the width increases (the McColl-Hartman effect). In the present study, it is shown that first, the McColl-Hartman effect for a rectangular barrier is dependent upon $L$ as the Gamow tunneling factor $\ensuremath{\theta}(k)$ vanishes, and ${\ensuremath{\tau}}_{g}$ is only independent of $L$ when $\ensuremath{\theta}(k)$ is large; and second, for a triangular barrier to model field emission, although ${\ensuremath{\tau}}_{g}$ can be large for small field, it vanishes when the energy matches the barrier height.

Physical Review A•Dec 9, 2021•9 citations

Dwell times for transmission and reflection

Nikolay Leonidovich Chuprikov

As was shown in quant-ph/0405028, the state of a tunneling particle can be uniquely presented as a coherent superposition of two states to describe alternative sub-processes, transmission and reflection. In this paper, on the basis of the stationary wave functions for these sub-processes, we give new definitions of the dwell times for transmission and reflection. In the case of rectangular potential barriers the dwell times are obtained explicitly. In contrast with the well-known Büttiker's dwell-time, our dwell time for transmission increases exponentially, for the under-barrier tunneling, with increasing the barrier's width. By our approach the well-known Hartman effect is rather an artifact resulted from an improper interpretation of the wave-packet tunneling and experimental data, but not a real physical effect accompanying the tunneling phenomenon.

arXiv:quant-ph/0502073•Feb 14, 2005•2 citations

Tunneling of a Wave Packet

Thomas E. Hartman

Analytic expressions are given for the time spent by a particle tunneling through a potential barrier. The expressions are derived for an incident wave packet which is initially Gaussian, centered about a point an arbitrary distance away from a rectangular potential barrier and moving toward the barrier with constant average velocity. Upon collision with the barrier, the packet splits into a transmitted and a reflected packet. The resultant transmission time is positive, nonzero and in principle measurable. Although the transmission time becomes quite large as the incident kinetic energy becomes very small, in general, for nonzero incident momentum and finite potential barriers which are neither very thick nor very thin, the transmission times are less than the time that would be required for the incident particle to travel a distance equal to the barrier thickness. The transmission times for metal-insulator-metal thin film sandwiches, given approximately by δt3=ℏ/(Efφ)12, where Ef is the Fermi energy of the metal, and φ the vacuum work function, are of the order of 10−16 sec, compared to RC time constants of about 10−13 sec.

Journal of Applied Physics•Dec 1, 1962•821 citations

Wavelet solution of the Schrödinger equation for a triangular potential barrier and applied results on a MOM tunnel junction

Emilia-Simona Malureanu, Florea Craciunoiu

Quantum tunneling has many important applications in electronics. In thin layers structures, when the tunneling phenomenon occurs at the metal-insulator contact barrier, as for the MIM or MOS structures, we are dealing with the Fowler-Nordheim field emission. In these structures we have to calculate the currents densities that go through. The current density is calculated based on the transmission coefficient through the barrier, which is determined by solving the Schrödinger equation. The calculation methods currently used lead to values that get predictions of emission current density too low. The transmission coefficient is calculated by wavelet approximation and the results are validated experimentally on a lateral MOM tunnel junction.

•Oct 1, 2013•0 citations

A delta barrier in a well and the exact time evolution of its eigenstates

Kevin L. Jensen, Jeanne M. Riga, Andrew Shabaev +2 more

The analytic nature of the transmission coefficient for a δ-function barrier makes it a useful tool to examine a variety of technologically important applications, such as photoemission from semiconductors with an alkali coating, the examination of tunneling times for wave packets incident on a barrier, and for parameterizing tunneling through the narrow barrier of a normal-superconducting point contact. The analytic model of a δ-function barrier inside a confining well is extended to the finite height and width rectangular barrier (a delta-function sequence). Methods to exactly evaluate the eigenstates are given and their dependencies are examined. The time evolution of a superposition of the lowest eigenstates is considered for barriers having comparable Gamow tunneling factors so as to quantify the impact of barrier height and shape on time evolution in a simple and exact system and, therefore, serve as a proxy for tunneling time. Last, density profiles and associated quantum potentials are examined for coupled wells to show changes induced by weaker and wider barriers.

Journal of Applied Physics•May 2, 2023•3 citations
Electron emission into a nanogap is not instantaneous, which presents a difficulty in simulating ultra-fast behavior using particle models. A method of approximating the transmission and reflection delay (TARD) times of a wave packet interacting with barriers described by a delta function, a metal–insulator–metal (MIM, rectangular) barrier, and a Fowler Nordheim (FN, triangular) barrier is given and has application to simulation. It is based on the superposition of a finite number of exact basis states obtained from Schrödinger’s equation, analogous to how quantum carpets are simulated. As a result, it can exactly and uniquely follow exponentially small tunneling currents. A Bohm-like trajectory is obtained from the time evolution of the density: it shows delay in both the transmitted and reflected packets that can be simply evaluated. The relations to prior studies of the analytic [Formula: see text]-function barrier and the Wigner distribution function (WDF) methods are described. A comparison of the TARD times is contrasted to alternate times in the Büttiker–Landauer (BL) and McColl–Hartman (MH) times; the MH approach is further reformulated explicitly in terms of Gamow factors to consider how the McColl–Hartman effect is to be related, particularly in the case of the FN barrier of field emission.
A numerical solution to the time evolution equation of the Wigner distribution function (WDF) with an accuracy necessary to simulate the passage of a wave packet past a barrier is developed, where quantum effects require high accuracy and fine discretization. A wave packet incident on a barrier, a portion of which tunnels through, demonstrates behavior that can define various characteristic transmission and reflection delay (TARD) times useful in the simulation of electron emission. A model for the TARD times is proposed that relies only on the asymptotic maxima of the position $\ensuremath{\rho}(x,t)$ and wave number $\ensuremath{\rho}(k,t)$ densities given by the WDF and applied to a ballistic trajectory model for the (faster) transmitted and (slower) reflected parts. The dependence of the TARD times on barrier width, symmetry, and abruptness is analyzed. For symmetrical barriers with characteristics similar to field emission barrier heights and widths, TARD times are on the order of a fraction of a femtosecond. The TARD times for when tunneling predominates are contrasted to tunneling times in the literature. Use of the TARD times in simulations of field emission in nanogap devices or to model ultrashort pulses generated under rapidly changing conditions for electron sources are proposed.
An electron wave packet tunneling through a barrier has a transmission (or ``group delay'') time ${\ensuremath{\tau}}_{g}$ that, for a rectangular barrier, is commonly held to become independent of the barrier width LLL as the width increases (the McColl-Hartman effect). In the present study, it is shown that first, the McColl-Hartman effect for a rectangular barrier is dependent upon LLL as the Gamow tunneling factor $\ensuremath{\theta}(k)$ vanishes, and ${\ensuremath{\tau}}_{g}$ is only independent of LLL when $\ensuremath{\theta}(k)$ is large; and second, for a triangular barrier to model field emission, although ${\ensuremath{\tau}}_{g}$ can be large for small field, it vanishes when the energy matches the barrier height.
As was shown in quant-ph/0405028, the state of a tunneling particle can be uniquely presented as a coherent superposition of two states to describe alternative sub-processes, transmission and reflection. In this paper, on the basis of the stationary wave functions for these sub-processes, we give new definitions of the dwell times for transmission and reflection. In the case of rectangular potential barriers the dwell times are obtained explicitly. In contrast with the well-known Büttiker's dwell-time, our dwell time for transmission increases exponentially, for the under-barrier tunneling, with increasing the barrier's width. By our approach the well-known Hartman effect is rather an artifact resulted from an improper interpretation of the wave-packet tunneling and experimental data, but not a real physical effect accompanying the tunneling phenomenon.
Analytic expressions are given for the time spent by a particle tunneling through a potential barrier. The expressions are derived for an incident wave packet which is initially Gaussian, centered about a point an arbitrary distance away from a rectangular potential barrier and moving toward the barrier with constant average velocity. Upon collision with the barrier, the packet splits into a transmitted and a reflected packet. The resultant transmission time is positive, nonzero and in principle measurable. Although the transmission time becomes quite large as the incident kinetic energy becomes very small, in general, for nonzero incident momentum and finite potential barriers which are neither very thick nor very thin, the transmission times are less than the time that would be required for the incident particle to travel a distance equal to the barrier thickness. The transmission times for metal-insulator-metal thin film sandwiches, given approximately by δt3=ℏ/(Efφ)12, where Ef is the Fermi energy of the metal, and φ the vacuum work function, are of the order of 10−16 sec, compared to RC time constants of about 10−13 sec.
Quantum tunneling has many important applications in electronics. In thin layers structures, when the tunneling phenomenon occurs at the metal-insulator contact barrier, as for the MIM or MOS structures, we are dealing with the Fowler-Nordheim field emission. In these structures we have to calculate the currents densities that go through. The current density is calculated based on the transmission coefficient through the barrier, which is determined by solving the Schrödinger equation. The calculation methods currently used lead to values that get predictions of emission current density too low. The transmission coefficient is calculated by wavelet approximation and the results are validated experimentally on a lateral MOM tunnel junction.
The analytic nature of the transmission coefficient for a δ-function barrier makes it a useful tool to examine a variety of technologically important applications, such as photoemission from semiconductors with an alkali coating, the examination of tunneling times for wave packets incident on a barrier, and for parameterizing tunneling through the narrow barrier of a normal-superconducting point contact. The analytic model of a δ-function barrier inside a confining well is extended to the finite height and width rectangular barrier (a delta-function sequence). Methods to exactly evaluate the eigenstates are given and their dependencies are examined. The time evolution of a superposition of the lowest eigenstates is considered for barriers having comparable Gamow tunneling factors so as to quantify the impact of barrier height and shape on time evolution in a simple and exact system and, therefore, serve as a proxy for tunneling time. Last, density profiles and associated quantum potentials are examined for coupled wells to show changes induced by weaker and wider barriers.