Over the past decades, advances in the generation of intense terahertz waves have enabled new researches on light–matter interaction. In particular, tabletop systems generating high-power terahertz pulses , , , , ,  provide broad opportunities to study ultrafast dynamics of matters in the nonlinear regime, through resonant and nonresonant coupling , , , , , , , , , , , . The low photon energy of terahertz waves allows resonant coupling with direct excitation of low-energy modes such as vibrations of crystal lattices, rotation of molecules, precessions of spins, and internal excitations of bound electron–hole pairs. Therefore, when the strong terahertz field is incident on the matter, it causes numerous low-energy excitations with large amplitudes beyond harmonic response. In addition to resonant manipulation, high electric field in a picosecond time scale can provide long-enough acceleration of free electrons and transient barrier distortion with strong voltage within a half cycle of driving field. This permits that the electron overcomes energy barriers of bound states. This light-field–driven tunneling can lead to interesting phenomena such as terahertz field emission , , , high harmonic generation , , , and nonlinear transmission , , ,  by making a dramatic change of material properties arising from the whole distortion of electronic potential without damaging the material due to the picosecond time scale of the terahertz pulses.
On the other hand, recent development of nanofabrication technology brings a large impact on various photonic researches including terahertz plasmonics , , , , , , , , , . Artificially designed metal structures in subwavelength scale are used to manipulate the response of the terahertz field, realizing a wide range of applications in sensors , , , , , , switches , , , , and filters , , , . Especially, a perfect conductor-like behavior of the noble metals in terahertz frequency range makes possible for strong light–matter interaction with nanoscale metal gaps or point geometries , , , , , . Electric field enhancement factor can reach a few thousands by coupling light with the metal structure, paving a way to study extreme nonlinear phenomena and providing ideal platform for the terahertz quantum plasmonics. Designing the metal nanostructures has therefore a significant impact on the terahertz responses in the quantum regime.
Terahertz quantum plasmonics is a research field that involves the studies of the quantum properties of light and its interaction with matter at nanoscale and angstrom scale in terahertz frequency regime. In this review, we focus on a nonresonant terahertz nonlinear phenomenon of electron tunneling driven by intense terahertz radiation that plasmonic metal nanostructures provide. We introduce a few basic building blocks of the metal nanostructures: a tip, coupled tips (dipole antenna), and metal gap for efficient generation of tunneling electrons driven by the terahertz waves. Meeting with the strong terahertz pulses, the metal nanostructures can enhance the electron tunneling, thereby leading to various extreme phenomena and useful applications.
2 Terahertz field enhancement in various metal structures
Metal structures such as nanotip, nanoantennas, and nanoslit and slot antennas have been of great research interest over the past two decades because of their potentials and functionalities for a broad range of applications such as chemistry, biology, and medical applications, as well as from basic research aspects , , , , , , , , , , , , , , . Strong field enhancement and confinement obtained at the subwavelength volume are the main resultant properties from interaction between electromagnetic waves and metal nanostructure, which are essential for nonlinear optics in nanoscale.
Figure 1 shows various metal nanostructures, nanotip (Figure 1A), nanoantenna (Figure 1B and C), nanoslit (Figure 1D), and nanoslot antenna (Figure 1E). Basically, the large field enhancement originates from the tight focusing of the terahertz electric field near the tip apex or at the nanogap induced by charge accumulation and capacitive coupling (Figure 1F) , . These two physical phenomena, charge accumulation and capacitive coupling, not necessarily of an orthogonal nature to each other, allow us to explain the efficiency of field enhancement in a broad spectral range from visible light to microwaves , , . Beyond the broadband spectral response of nanotip and nanoslit, dipole antennas and slot antennas are employed for the resonant excitation, furthering the capacities for manipulation and control of spatially and temporally localized electrons. In the next section, we will review recent progress in ultrafast control of electrons using terahertz plasmonics, which is one of the most studied terahertz nonlinear optics that requires intense electromagnetic field.
3 Terahertz field–induced quantum electron tunneling
Tunneling is one of the most fundamental quantum mechanical phenomena; wave function of matters can penetrate and eventually transmit through a classically forbidden potential barrier. Strong electrical bias on the metal surface applied by the terahertz pulse excitation modifies the barrier and can emit electrons from the metal. This tunneling current directly affects the experiments exploiting strong terahertz pulsed sources with metal nanostructures, and we therefore need to consider the tunneling phenomena quantitatively when analyzing the data.
When a terahertz electric field is incident on two metallic electrodes separated by a thin insulating film, electrons can flow through the potential barrier by tunneling effect (Figure 2). Tunneling current density J across a one-dimensional barrier by applying potential difference V at 0 K is given by , 
where m is the electron mass; e, the electron charge; h, the Plank constant; Ef, the Fermi level of the metal; E, the tunneling electron energy; D(E)=exp(−4πΔw(2m(η+ϕm−E))1/2/h), the tunneling probability factor calculated by WKB approximation; the mean value of the barrier height; ϕ(x)=ϕ−eVx/w, the barrier potential profile affected by the THz electric field; ϕi(x)=−1.15e2w2ln(2)/16πεwx(w−x), the approximated form of barrier modification by image force; w, the thickness of the insulating layer; ε, the dielectric constant of the insulating layer; and Δw=w2−w1 is the effective barrier width, where w1 and w2 are acquired from the zeros of ϕ(x)+ϕi(x)=0. Because of the exponential dependence for the applied electric field V/w, the tunneling current exhibits a highly nonlinear character. For example, even a 1-nm increase for a 1-nm gap with a constant voltage (or the field strength reduced by half while fixing the gap size) changes the tunneling current by several orders of magnitudes smaller.
The first term of (1) indicates tunneling current from electrode 1 (left) to electrode 2 (right) where the energy level of electrode 2 is lowered by amount of eV, and the second term indicates the current flowing in a reverse direction. When the applied field is strong enough to suppress the Fermi level of electrode 2 below the bottom of conduction band of electrode 1, i.e. eV>Ef, electrons cannot tunnel from electrode 2 to electrode 1 because of the absence of available states. In this case, the second term of (1) can be neglected and J can be written by
with a constant C, which is the Fowler–Nordheim tunneling formula for field emission from a single metallic electrode . In this case of field emission, the barrier potential denotes the work function of the metal.
3.1 Metal nanotip
3.1.1 Terahertz scanning tunneling microscope
Terahertz scanning tunneling microscope (STM) (THz-STM) is capable of imaging surfaces with atomic spatial resolution and the subpicosecond time resolution (<0.5 ps), much faster than conventional STM system , , , , , . The basic principle of THz-STM is quantum tunneling through the gap between the tip and sample resulting from a transient voltage induced by terahertz pulses. The total time-dependent voltage across the tunnel junction is the sum of the DC bias voltage and the terahertz pulse-induced transient voltage (left of Figure 3B). Because of the nonlinear characteristic of current–voltage relation in the tunnel junction, a time-dependent tunneling current with a rectified component is produced (right of Figure 3B). The terahertz-induced tunneling current pulses result in the time-averaged shift in STM current–voltage curves.
The spatial resolution limited by diffraction of terahertz waves (~1 mm) can be overcome by a scanning probe tip, which is an atomically sharp metallic tip with the resolution on the order of 10 nm. The free-space traveling terahertz pulses are focused on the STM tip acting as a broadband antenna ,  that couples the electric field of the terahertz pulses to the tip–sample junction via a Sommerfeld mode . The strongly enhanced terahertz electric field at the tip apex enables to produce the transient bias voltage with maintaining nanometer spatial resolution. A recent work reports imaging and tunneling spectroscopy of conductive surfaces with atomic resolution (~0.3 nm), which shows the image of ultrafast nonequilibrium tunneling dynamics for individual atoms on a silicon surface .
The antenna-like coupling of terahertz waves to the metal nanotip in STM system also results in polarization-dependent current–voltage characteristic. Figure 3C shows the STM current–voltage curves for a bare, flat, highly ordered pyrolytic graphite (HOPG) sample, revealing that the terahertz pulse generates an increased tunnel current only when the electric field of the terahertz pulse is polarized parallel to the metal tip , , .
The THz-STM system is used to study ultrafast carrier dynamics of single nanoparticles in subpicosecond time resolution by combining pump–probe technique. The terahertz-induced tunnel current with a different rectified component before and after optical excitation is measured by optical-pump THz-STM–probe system as seen in Figure 3D. Figure 3E is the change of tunnel current caused by photoexcitation in InAs nanodot as a function of pump–probe time delay. After photoexcitation, the electron trapping into InAs nanodots increases local bias voltage. resulting in the increase of rectified component in tunnel current. In the time-resolved THz-STM signal, the rise time of ~500 fs and the decay time of ~1 ps were measured with a nanoscale spatial resolution . Also, by introducing pump–probe technique where one terahertz pulse (pump) modifies transiently the rectified currents induced by the second one (probe), the ultrafast dynamics of an individual molecule such as a coherent vibration of the molecules was visualized in real space with an atomic scale .
In the THz-STM system, the single-cycle terahertz electric field transiently produces asymmetric bias voltage based on DC bias voltage, leading to the rectified component of tunnel current. In turn, the rectified tunnel current direction in the tip–sample junction is determined by the DC bias voltage because of the symmetric terahertz pulse shape in time domain. By using carrier–envelope phase (CEP)–controlled terahertz electric fields, the tunnel electrons in the tip–sample junction can be manipulated coherently. To experimentally control the CEP of single-cycle terahertz electric field, a pair of cylindrical lenses or spherical lenses for π/2 or π phase change, respectively (Figure 4A). The CEP-controlled terahertz electric field leads to the coherent distortion of potential barrier directly (Figure 4B). For the symmetric terahertz pulse, the direction of ultrafast electron tunneling is determined by the DC bias voltage, where the electrons tunnel from the metal tip (sample) to the sample (metal tip) under a positive (negative) DC bias voltage. In case of the asymmetric terahertz pulses, however, ultrafast tunnel current induced by strong terahertz field has unidirectional behavior while sweeping the DC bias voltage from negative to positive value. The direction of ultrafast tunneling current strongly depends on the CEP of terahertz pulses (Figure 4C, D).
This coherent manipulation of electron motion is enabled by terahertz-induced bias voltage, which is strong enough to distort barrier height in the tip–sample junction. The strong transient voltage is induced by the large field enhancement realized by the broadband antenna effect in the metal nanotip. Therefore, for the shaper metal tip accompanying stronger field enhancement and confinement, the larger tunnel current is measured in a given incident terahertz electric field (Figure 5A). It is worth to note that the rectified electrons show the same curves as a function of electric field scaled by enhancement factor, indicating that the tip geometry plays an important role to reach strong field regime at a given incident field (Figure 5B). Furthermore, the strong field enhancement by metal nanotip assists to reach a new regime, an extremely high-voltage regime where an additional space-charge potential plays an important role in the strong saturation of the tunneling electrons.
3.1.2 Terahertz field emission
Terahertz field emission is one of terahertz field–driven processes enabled by the high field enhancement in metal nanotip. Figure 6A shows the experimental setup of terahertz-induced field emission from a sharp tungsten tip in high-vacuum condition (10−8 mbar range), in which an incident single-cycle terahertz pulse is focused by an off-axis parabolic mirror and the emitted electrons are detected by a time-of-flight (TOF) electron spectrometer and a microchannel plate detector. The enhanced transient terahertz pulses created by antenna-like coupling of metal tip induce nonlinear electron emission from the nanotip, through ac-tunneling into vacuum and acceleration process by the locally enhanced terahertz near field. The emitted electron kinetic energy distributions are characterized by transient strong terahertz field. Figure 6C displays TOF spectra for two different terahertz pulse strengths. The spectral shape stems from nonlinear emission at maximum field strength and subsequent acceleration in the terahertz-induced momentary potential. Another recent work reported that large local field enhancement of terahertz pulses in metal nanotips is responsible for producing electrons with energies exceeding 5 keV, via terahertz field–driven electron emission process .
Tip-enhanced single-cycle terahertz pulse can be used for control of photoelectron dynamics in dual-frequency excitation of metal nanotip. For example, assisted by optical field enhancement, near-infrared pulses (800-nm wavelength) induce nonlinear photoelectron emission from the nanotip, which is gated and streaked by the enhanced terahertz electric field (Figure 6D, E). This near-field streaking technique can be applicable for nanoscopic version of terahertz streaking device .
3.2 Metal nanoantenna
3.2.1 Photon-assisted tunneling in single-molecule transistors
Electron transport through single molecules has been intensively studied because of its potential for device applications such as single-electron transistor. Beyond static properties of the single-electron transistor, its dynamical transport properties have been investigated by terahertz waves because typical energy scales in single-molecule transport lie mostly in the terahertz frequency regime such as tunneling times, vibron energies, orbital energy spacing, charging energies, and so on , . To couple the terahertz waves to a single molecule trapped between nanogap electrodes, however, it is inevitable to utilize the plasmonic effects of metal electrodes with a nanometer sized gap.
Bowtie-shaped antenna for the source and drain electrode can be employed for single-molecule transistors to achieve a good coupling efficiency between the terahertz waves and the tunneling electrons. By using the plasmonic effect, furthermore, the terahertz field in the gap is enhanced by a factor of ~105 in a nanometric volume, enabling strong terahertz field–induced nonlinear response such as photon-assisted tunneling and two-photon absorption , . As shown in Figure 7A–C, it was reported that the strong terahertz fields can create additional electron tunneling pathways in a single-molecule transistor in which a bowtie-shaped metal antenna plays a role as a metal electrode. Experimentally, the additional lines parallel to the ground-state lines in Coulomb stability diagrams were observed in the single-electron transport and Coulomb blockaded regions (Figure 7C). The observed energy separations ~10 meV (middle) corresponding to the photon energy of the incident terahertz waves indicate resonant terahertz-assisted tunneling effect (middle of Figure 7C), and the data for the higher terahertz irradiation show that two-photon absorption takes place (right of Figure 7C). Furthermore, a recent work demonstrated low-energy excitation in single fullerene molecules transistor by vibron-assisted tunneling promoted by the terahertz-induced oscillation of the molecule (Figure 7D) . The vibrational modes of molecule (~2 and ~4 meV)  create new tunneling paths for electrons observed in terahertz-induced photocurrent measured by using the time-domain terahertz autocorrelation measurement combined with bowtie antenna-based single-molecule transistor geometry (Figure 7E, F).
3.2.2 Terahertz field emission
As shown in Figure 8A, the locally enhanced terahertz electric field in the gap of a dipole antenna contributes to electric field emission into vacuum. The resonant excitation of the dipole antenna by an intense terahertz pulse of ~100 kV/cm enables the field strength in the dipole gap to reach up to ~10 MV/cm, which requires a field enhancement factor of ~100 (Figure 8B, C). The dipole antennas for the terahertz field–driven field emission experiment are fabricated on a thin SiNx film, and then the SiNx substrate in the vicinity of the gap is selectively etched away (Figure 8D, F). This allows the motion of free electrons emitted across the gap. Figure 8D–G show scanning electron microscopy (SEM) images of dipole antenna with a capacitive gap of around 2 μm before and after terahertz exposure for 1 h. In the vicinity of the gap, the gold has moved inward with reducing the gap to around 1 μm. The gold deformation occurring near the gap of dipole antennas results from electromigration due to high-energy field–emitted electrons. The mass transfer of gold produces current channel across the gap; therefore, the spectral change of dipole antennas by strong terahertz field irradiation can be modeled by a decrease in resistivity between the antennas (Figure 8H).
3.3 Metal nanoslit and slot antenna
3.3.1 Terahertz nonlinear transmission response
The vertically aligned metallic nanogaps can initiate the terahertz nonlinear response arising from the quantum tunneling effect. Noble metals such as gold or copper are good conductors in terahertz frequencies; therefore, they can induce extremely large field enhancement enough to reach the quantum regime with the few nanometer-sized gap structures . Figure 9A shows the experimental scheme for the terahertz quantum plasmonics using intense terahertz pulse ,  and metallic nanogap structures . Here, nanogaps are fabricated with the gap sizes of 5 and 1.5 nm by Al2O3 layers sandwiched by gold film, and copper single-layer graphene–copper for extreme terahertz field enhancement . Figure 9B presents transmitted amplitudes through a 1.5-nm gap resonant ring antenna for different incident terahertz field strengths, which are normalized by the maximum field amplitudes passing through a bare quartz substrate. As the incident field approaches 150 kV/cm, the transmitted field is considerably reduced by the terahertz field–induced tunneling across the nanogaps. Figure 9C shows the temporally integrated transmittance (T ≡ ∫Esample2(t)dt/∫Ereference2(t)dt, where Esample is field amplitude through the nanogap sample and Ereference is the field amplitude though the bare substrate), which is normalized by the transmittance T0 of the minimum incident field strength. The transmittance value eventually reduces by 50% as the field inside the gap Egap increases to 5 V/nm. It is worth to note that the similar field strength of V/nm in static field regime can damage the gap .
By pushing the limit of the gap size reducing down to the single atomic scale, we can pursue to a dramatic terahertz nonlinear response from the nanogap structures. The bottom image of Figure 9A shows the extreme, a metallic nanogap separated by a single-layer graphene, which is made by recently developed lithography method . These gaps with a single layer of carbon atoms in the middle have been approximated by a single 3-angstrom-wide gap with a dielectric constant of three , , across which strong quantum tunneling may occur , , . Figure 9D shows measured time traces of terahertz field amplitudes transmitted through the sample for various incident terahertz pulse strengths, normalized by the maximum field amplitudes transmitted through the quartz substrate only. By increasing the incident terahertz electric field up to E0,max ~192 kV/cm, the normalized transmitted amplitude in the time domain dramatically decreases, down to 0.16% at the maximum. This gives rise to a 97% reduction in normalized transmitted power. Figure 9E shows the normalized transmittance T and field enhancement factor (inset) at the maximum in time domain, as a function of the strength of the incident terahertz electric field.
When the applied electric field is sufficiently strong to flow tunneling current across the gap, the effective refractive index of the gap medium becomes larger, which in turn limits the field enhancement , . This self-limited transmission by quantum tunneling can be observed without damaging the gaps by using a single-cycle terahertz pulse (~1 ps) , . On the other hand, it is not easy to observe the nonlinear transmission by using optical pulses . This trend is summarized in Figure 10. It shows transmitted peak power vs. incident power for both terahertz and near-infrared frequencies. As we can see, the terahertz nonlinearities are obvious for both structures. As the incident power increases, the transmitted power deviates further away from the linear regime (dashed line) due to the quantum tunneling induced by the terahertz field. On the contrary, nonlinearity is hardly seen in the near-infrared regime for both nanogap samples. Before reaching the nonlinear transmission regime, the sample is damaged in the near-infrared. This behavior can be explained by estimating tunneling conductivity and displacement conductivity in both frequency ranges. The results indicate that the tunneling conductivity σt is the same for terahertz and near-infrared frequencies for the same field strength, while the frequency-dependent displacement conductivity (σd=−iωϵ0ϵgap) takes relatively large portion compared to the tunneling conductivity in the near-infrared . As a result, tunneling-induced transmission modulation is negligible in the near-infrared.
3.3.2 Terahertz tunneling transport
The nanogap structure in configuration of a metallic slit, composed of optically thick metal films extended far for many wavelengths in lateral direction and optically transparent insulating layer sandwiched between the metals, can be modeled by an equivalent circuit with appropriate parameters ,  as shown in Figure 11A. Here, the fields of the incident and transmitted light coincide with the applied current and voltage that are studied in impedance measurements. When an electromagnetic wave has its magnetic field aligned parallel to the slit (TM polarization), a current I0(t) is induced in the metal, which flows toward the gap. For normally incident plane waves, the current is expressed in terms of an effective surface current density K0=n×(2Hinc), where n is the unit vector normal to the metal surface and Hinc is the magnetic field vector of the incident wave. Therefore, the current I0 can be expressed with the incident field Hinc by the relation I0=|K0|×l with slit length scale l, which corresponds to the beam spot size illuminating the sample. The voltage across the gap V(t) simply equals the width of the gap times the (enhanced) electric field strength inside the gap. As derived from the Kirchhoff integral formalism , , , , the near-field enhancement is directly proportional to the transmitted field determined in the far field. In short, evaluating the incident and transmitted fields allows for a direct measure of the applied current and induced voltage at the gap, respectively. For a narrow slit, whose width is much narrower than the metal thickness, a capacitor-like behavior is expected , . Then, the voltage across gap V(t) is connected to the applied current I0(t) by the relation V(t)=(∫ I0(t) dt)/C, where C is the capacitance of the nanogap structure. In particular, if we consider a harmonic time dependence, given by e−iωt, the equation is rewritten as V=I0/(−iωC) in the frequency domain. This is illustrated in Figure 11B by the phase lag of π/2 between the induced voltage and the applied current. The great advantage of the circuit model analysis is on the ability to provide quantitative agreement with the experimental results without any fitting parameters but by considering only the amplitude profile of the incident pulse and the dimensions of the capacitor (such as the gap size, metal thickness, and refractive index of the insulating layer), which are known prior to the measurements. For example, Figure 11C shows the terahertz transmission through a slit that has gap widths of 1.5 and 5 nm and compared with the circuit modeling results.
As we have seen, the nonlinear effect can alter the terahertz transmission by the tunneling electrons across the gap, and the above simple capacitor picture for the slit structure may not be applied to the nonlinear response . It could be necessary to modify the optical properties of the insulating medium by including additional parameters such as resistance R. When intense terahertz pulse (peak value of E0=30 kV/cm) is incident on a nanogap whose gap dimension is only a few angstroms, a maximum of 2 V is induced across the atomic gap, as indicated in Figure 12A. Determined from the structural dimensions of the gap and an effective permittivity εgap of 3 , , , a capacitance of C=2.7 pF alone fails to correctly describe the features of the transmitted terahertz wave in this case. Electron tunneling through the gap prevents charge accumulation at the metal edges and causes a reduced voltage across the gap, and the shape of the transmitted pulse closely resembles the profile of the single-cycle incident pulse. Thus, a finite tunneling resistance Rt connected in parallel to the innate slit capacitance should be added to the total impedance, which relates the transient voltage to the current applied by the incident radiation . If a constant tunneling resistance is assumed over the studied voltage range, a value of 0.07 Ω offers a good fit to the experimental results. On the other hand, the circuit model is able to provide a differential equation that links the measured applied current I0(t) and induced voltage V(t), and the tunneling resistance can be deduced without assuming that the resistance value is fixed; the tunneling current I0(t)−C dV/dt is obtained from the source current I0(t) and measured voltage across the gap V(t). The extracted tunneling current shown in Figure 12B is proportional to the voltage across the gap V(t), suggesting an ohmic response of the tunneling electrons through the gap. Figure 12C is a current–voltage characteristic obtained from the mean value of the data points within every 0.1-V interval of the voltage across the gap; error bars indicate the standard deviation. The plot agrees with the slope of 0.07 Ω that was evaluated from the fit in Figure 12A.
Because the tunneling resistance and the capacitance C of the slit form a parallel circuit, the resistance being measured should be comparable to the impedance Z=1/(−iωC) or smaller to be detectable. Therefore, the factor 1/|Z| specifies the resolution of the circuit analysis when measuring the tunneling current. Smaller tunneling currents are measurable if the slit capacitance is decreased; presumably, the most efficient way to achieve this in the slit structure is to shorten the slit length l so that the slit reduces to a close-ended slot. Additionally, transmission of electromagnetic waves becomes challenging, as extremely narrow gaps are required to observe quantum effects. Resonance structures, such as slot antennas, allow an increase in transmittance for certain bandwidths, which can also be improved for broadband purposes using multiresonance structures , and exhibit a better signal-to-noise ratio than slits that suffer from a capacitor-like broad 1/f-type spectral response. Figure 12D and E show the time trace of the voltage across the gap and its Fourier analysis through ring-shaped slot antenna whose resonance is placed near 0.3 THz for an incident pulse strength of E0=20 V/cm. For this case, the circuit analysis is also possible to describe the terahertz response by an RLC (a resistor (R), an inductor (L), and a capacitor (C)) circuit , as depicted in the inset of Figure 12D. The parameters R, L, and C are determined from the complex impedance for a parallel RLC circuit to acquire the best fit with the transmitted pulse V(t) from the measured incident pulse I0(t). Moreover, a phase parameter θ is introduced to compensate for an additional phase shift in the transmitted wave due to coupling effects between nearby antennas, which is predicted from calculations based on modal expansion . Accordingly, an equivalent waveform is reproduced in Figure 12D by using the equation together with the values of R=5.83 Ω, L=1.72×10−12 H, C=1.69×10−13 F, and θ=20.6°. These four parameters are essential to fit the amplitude of the resonance peak, the resonant frequency, the line width, and the phase in the frequency domain. The same procedure is carried out for an incident pulse of E0=30 kV/cm. For such moderate intensities, where quantum effects are negligible, the same values of the four parameters also thoroughly duplicate the experimental data, as shown in Figure 12E.
3.3.3 Terahertz tunneling rectification
Figure 13A shows the strategy for the terahertz tunneling rectification by using lateral, ring-shaped tunneling barriers, encasing a metallic island that is surrounded by a metallic plane . By adjoining the one-dimensional tunneling junctions in two dimensions with triangle shape , the tunneling current flowing across the quantum barriers is determined not only by the phase of the external electromagnetic pulse , , , but also by the geometry of the barrier (i.e. lateral symmetry of the ring). Using the proposed method, it is possible to control the ultrafast tunneling process additionally by modifying the lateral shape of the ring and by simply changing the polarization of incoming pulses. The time-integrated total tunneling currents through the barriers are measured directly by attaching electrical probes on the sample surface as shown in Figure 13D. A THz polarizer is placed to control the direction of the surface current K. Figure 13E shows measured tunneling rectification currents across triangular (side length of 70 μm and perimeter of 210 μm shown in Figure 13B, C) and square (side length of 25 μm and perimeter of 100 μm) nanogap rings (gap size of 2 nm) as a function of the terahertz polarizer angle, i.e. changing incident terahertz polarization. The results show strikingly different behaviors as mentioned above. The asymmetric potential distribution along the equilateral triangle results in a net tunneling current through the contour, while the potential distribution at any point of a square is mostly counterbalanced by its corresponding point across the center, independent of the polarization of the incident pulse. Therefore, the triangle shape barrier can be used to accumulate more tunneling charges inside the metallic island under a single terahertz pulse excitation than the rectangle one, which enhances directly the experimental signal. Figure 13F displays polar plots of the tunneling currents for the triangular and square geometries. As expected, the total current vs. polarization angle shows the threefold rotational symmetry of an equilateral triangle.
Sending a femtosecond optical pulse with the terahertz pulse can provide another control parameter for the ultrafast modulation of the tunneling current. As shown in Figure 13G, the electric fields of the femtosecond optical pulse can be focused to a specific position in the ring and summed with picosecond terahertz electric fields. The optical pulse rapidly distorts the local potential barrier under the quasi-constant terahertz field, generating an additional, local tunneling current. The optically driven current is sensitively affected by the terahertz field strength at the specific position on the barrier and by the time delay between the terahertz and optical pulses, thus providing a way to visualize the spatiotemporal dynamics of the terahertz potential across the barriers. An interesting aspect of the optically modulated quantum barriers is the position-dependent ultrafast optical gating, as illustrated in Figure 13H.
Using the ring-shaped barriers and the time-resolved detection method of terahertz tunneling current by the optical pulse, terahertz full-wave rectification can be achieved, as shown in Figure 13I. Here, the spot size of the optical beam is expanded to cover the entire ring, which makes it possible to probe the total terahertz current response through the ring. Additionally, a DC bias is applied to generate a unidirectional (i.e. into or out of the loop) terahertz current across the ring barrier independent of the polarity of the incident terahertz pulse. As a result, the terahertz wave generates a fully rectified terahertz tunneling current depending on the DC bias strength, and this transient is visualized by the optical pulse (Figure 13J).
4 Discussions and outlook
In this review, we surveyed recent works concerning terahertz control of electrons using various metallic nanostructures. The intense terahertz pulses can manipulate tunneling electrons by ultrafast modification of the potential barriers in metal nano- and subnano-junctions. Here, three different metal nanostructures utilized for terahertz-induced quantum tunneling phenomena have been discussed: metal nanotip, nanoantenna, nanoslit, and slot antenna. The first one, terahertz-coupled metal nanotip, provides atomic scale imaging and picosecond time resolution in combination with STM system and can generate localized high-energy electron bursts. The second structure, dipole antenna system with various shapes, enables terahertz-assisted single-electron tunneling based on the resonant vibrational excitation of molecule, which is trapped inside the tips. In the case of metal nanoslit and slot antennas, terahertz-induced voltage and tunneling current can be extracted experimentally with a subpicosecond time resolution, leading to quantitative analysis of tunneling currents. Further developments of terahertz quantum plasmonics in terms of tunneling electron manipulation will require improvement of spatiotemporal resolutions and advances in the generation and detection of terahertz pulses, as well as easy-to-implement subatomic scale distance controls. Beyond terahertz field–driven electron tunneling, the coherent tunneling currents can be controlled by a few femtosecond phase–locked driving pulses using a single cycle of mid and near-infrared light.
Hirori H, Doi A, Blanchard F, Tanaka K. Single-cycle terahertz pulses with amplitudes exceeding 1 MV/cm generated by optical rectification in LiNbO3. Appl Phys Lett 2011;98:091106. CrossrefGoogle Scholar
Reimann K, Smith RP, Weiner AM, Elsaesser T, Woerner M. Direct field-resolved detection of terahertz transients with amplitudes of megavolts per centimeter. Opt Lett 2003;28:471–3. CrossrefPubMedGoogle Scholar
Hebling J, Yeh K-L, Hoffmann MC, Bartal B, Nelson KA. Generation of high-power terahertz pulses by tilted-pulse-front excitation and their application possibilities. J Opt Soc Am B 2008;25:B6–19. CrossrefGoogle Scholar
Sell A, Leitenstorfer A, Huber R. Phase-locked generation and field-resolved detection of widely tunable terahertz pulses with amplitudes exceeding 100 MV/cm. Opt Lett 2008;33:2767–9. PubMedCrossrefGoogle Scholar
Luo CW, Reimann K, Woerner M, Elsaesser T, Hey R, Ploog KH. Phase-resolved nonlinear response of a two-dimensional electron gas under femtosecond intersubband excitation. Phys Rev Lett 2004;92:047402. CrossrefPubMedGoogle Scholar
Yamaguchi K, Nakajima M, Suemoto T. Coherent control of spin precession motion with impulsive magnetic fields of half-cycle terahertz radiation. Phys Rev Lett 2010;105:237201. PubMedCrossrefGoogle Scholar
Katayama I, Aoki H, Takeda J, et al. Ferroelectric soft mode in a SrTiO3 thin film impulsively driven to the anharmonic regime using intense picosecond terahertz pulses. Phys Rev Lett 2012;108:097401. CrossrefGoogle Scholar
Yoshioka K, Katayama I, Minami Y, et al. Real-space coherent manipulation of electrons in a single tunnel junction by single-cycle terahertz electric fields. Nat Photon 2016;10:762–5. CrossrefGoogle Scholar
Tarekegne AT, Iwaszczuk K, Zalkovskij M, Strikwerda AC, Jepsen PU. Impact ionization in high resistivity silicon induced by an intense terahertz field enhanced by an antenna array. New J Phys 2015;17:043002. CrossrefGoogle Scholar
Strikwerda AC, Zalkovskij M, Iwaszczuk K, Lorenzen DL, Jepsen PU. Permanently reconfigured metamaterials due to terahertz induced mass transfer of gold. Opt Express 2015;23:11586–99. PubMedCrossrefGoogle Scholar
Keiser GR, Seren HR, Strikwerda AC, Zhang X, Averitt RD. Structural control of metamaterial oscillator strength and electric field enhancement at terahertz frequencies. Appl Phys Lett 2014;105:081112. CrossrefGoogle Scholar
Lange C, Maag T, Hohenleutner M, et al. Extremely nonperturbative nonlinearities in GaAs driven by atomically strong terahertz fields in gold metamaterials. Phys Rev Lett 2014;113:227401. PubMedCrossrefGoogle Scholar
Koo S, Kumar MS, Shin J, Kim D, Park N. Extraordinary magnetic field enhancement with metallic nanowire: role of surface impedance in Babinet’s principle for sub-skin-depth regime. Phys Rev Lett 2009;103:263901. CrossrefPubMedGoogle Scholar
Schertz F, Schmelzeisen M, Kreiter M, Elmers HJ, Schönhense G. Field emission of electrons generated by the near field of strongly coupled plasmons. Phys Rev Lett 2012;108:237602. CrossrefPubMedGoogle Scholar
Toma A, Tuccio S, Prato M, et al. Squeezing terahertz light into nanovolumes: nanoantenna enhanced terahertz spectroscopy (NETS) of semiconductor quantum dots. Nano Lett 2015;15:386–91. PubMedCrossrefGoogle Scholar
Shibata K, Umeno A, Cha KM, Hirakawa K. Photon-assisted tunneling through self-assembled InAs quantum dots in the terahertz frequency range. Phys Rev Lett 2012;109:077401. PubMedCrossrefGoogle Scholar
Yoshioka K, Katayama I, Arashida Y, et al. Tailoring single-cycle near field in a tunnel junction with carrier–envelope phase-controlled terahertz electric fields. Nano Lett 2018;18:5198–204. CrossrefGoogle Scholar
Du S, Yoshida K, Zhang Y, Hamada I, Hirakawa K. Terahertz dynamics of electron–vibron coupling in single molecules with tunable electrostatic potential. Nat Photon 2018;12:608–12. CrossrefGoogle Scholar
Baek IH, Kang BJ, Uk Jeong Y, Rotermund F. Diffraction-limited high-power single-cycle terahertz pulse generation in prism-cut LiNbO3 for precise terahertz applications. J Opt Soc Korea 2014;18:60–4. CrossrefGoogle Scholar
Marinica DC, Kazansky AK, Nordlander P, Aizpurua J, Borisov AG. Quantum plasmonics: nonlinear effects in the field enhancement of a plasmonic nanoparticle dimer. Nano Lett 2012;12:1333–9. CrossrefPubMedGoogle Scholar
Huber R, Tauser F, Brodschelm A, Bichler M, Abstreiter G, Leitenstorfer A. How many-particle interactions develop after ultrafast excitation of an electron-hole plasma. Nature 2001;414:286–9. CrossrefPubMedGoogle Scholar
Kyoung JS, Seo MA, Park HR, Ahn KJ, Kim DS. Far field detection of terahertz near field enhancement of sub-wavelength slits using Kirchhoff integral formalism. Opt Commun 2010;283:4907–10. CrossrefGoogle Scholar
About the article
Published Online: 2020-01-28
Published in Print: 2020-02-25
Funding Source: National Research Foundation of Korea
Award identifier / Grant number: NRF-2019R1G1A1007912
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIP: NRF-2015R1A3A2031768, Funder Id: http://dx.doi.org/10.13039/501100001321, NRF-2019R1G1A1007912) (MOE: BK21 Plus Program21A20131111123, Funder Id: http://dx.doi.org/10.13039/501100001321) and the 2019 Research Fund of UNIST.
Citation Information: Nanophotonics, Volume 9, Issue 2, Pages 435–451, ISSN (Online) 2192-8614, DOI: https://doi.org/10.1515/nanoph-2019-0436.
© 2020 Taehee Kang, Young-Mi Bahk and Dai-Sik Kim, published by De Gruyter, Berlin/Boston. This work is licensed under the Creative Commons Attribution 4.0 Public License. BY 4.0