Nanoantenna Enhanced Terahertz Interaction of Biomolecules
Nanoantenna Enhanced Terahertz Interaction of Biomolecules
Adak, Subham;Tripathi, Laxmi Narayan
1 1 , Subham Adak and Laxmi Narayan Tripathi Department of Physics, Birla Institute of Technology, Mesra, Ranchi-835215, Jharkhand, India. (Dated: March 11, 2019) Terahertz time-domain spectroscopy (THz-TDS) is a non-invasive, non-contact and label-free tech- nique for biological and chemical sensing as THz-spectra is less energetic and lies in the characteristic vibration frequency regime of proteins and DNA molecules. However, THz-TDS is less sensitive for detection of micro-organisms of size equal to or less than l/100 (where, l is wavelength of incident THz wave) and, molecules in extremely low concentrated solutions (like, a few femtomolar). Af- ter successful high-throughput fabrication of nanostructures, nanoantennas and metamaterials were found to be indispensable in enhancing the sensitivity of conventional THz-TDS. These nanostruc- tures lead to strong THz ﬁeld enhancement which when in resonance with absorption spectrum of ab- sorptive molecules, causing signiﬁcant changes in the magnitude of the transmission spectrum, there- fore, enhancing the sensitivity and allowing detection of molecules and biomaterials in extremely low concentrated solutions. Hereby, we review the recent developments in ultra-sensitive and se- lective nanogap biosensors. We have also provided an in-depth review of various high-throughput nanofabrication techniques. We also discussed the physics behind the ﬁeld enhancements in sub-skin depth as well as sub-nanometer sized nanogaps. We introduce ﬁnite-difference time-domain (FDTD) and molecular dynamics (MD) simulations tools to study THz biomolecular interactions. Finally, we provide a comprehensive account of nanoantenna enhanced sensing of viruses (like, H1N1) and biomolecules such as artiﬁcial sweeteners which are addictive and carcinogenic. CONTENTS I. INTRODUCTION 1,2 I. Introduction 1 Optical antennas, devices converting free optical ra- diations into localized ones, and vice versa are impor- II. High-Throughput Wafer-scale Fabrication tant components in modern nanophotonic devices. In Techniques 2 the recent decade, metallic nanoantenna has attracted much consideration with the signiﬁcant advancement III. Terahertz-Field Enhancement using of nanofabrication technologies, which provide the fea- Nanoantennas 4 sibility of fabrication at nanometre resolution, enabling A. THz Field Enhancement in Metal Nanogaps fascinating nano-optoelectronic devices. Since metals at Sub Skin-Depth Regime 5 exhibit ﬁnite conductivities and plasmonic effect at op- tical frequencies, therefore the optical ﬁeld response is B. THz Field Enhancement in Sub-Nanometer and Angstrom Metallic gaps: Quantum dictated by the resonance oscillation of electron clouds (surface plasmons, or SP) in the metal, coupling with Regime 7 the incident photon and sustaining surface plasmon po- C. The Finite-Difference Time-Domain (FDTD) 3,4 Algorithm : A Computational laritons (SSPs). Henceforth, at optical frequencies a Electromagnetism Approach 8 nanoantenna, depending upon the material properties, resonates for considerably shorter effective wavelength, 5,6 IV. THz-TDS and Sensing of Molecules and rather than half-wavelength of free-space. Biomolecules 9 Electromagnetic interaction of nanoantennas has lead A. Computational THz-TDS using Molecular to huge local ﬁeld enhancement and conﬁnement Dynamics (MD) Simulation 10 in nano-size volumes, enabling broadband applica- 7–15 B. THz-Nanoantenna and Metamaterials tion in nano-optics. Terahertz (THz) wave fun- Sensing of Micro-organisms, Organic- and nelling through nanogaps/nanoantennas has resulted Biomolecules 13 16–40 huge ﬁeld enhancement. Figure 1 summarizes the various investigations performed to study the THz V. Conclusion and Future Outlook 16 ﬁeld enhancement with different gap sizes. Such investigations have lead to the development of the VI. Acknowledgement 17 terahertz technology, which are now exploited in 41–46 47–49 50–52 electronics, photonics, medical sciences, VII. Conﬂict of Interest 18 53–55 56,57 military, security and conservation of cultural 58,59 heritages. These studies also lead to the improvisa- VIII. Keywords 18 60–76 tion of terahertz devices. 77–82 References 18 Terahertz time-domain spectroscopy (THz-TDS) arXiv:1903.03415v1 [physics.optics] 8 Mar 2019 2 > 100 nm 100 – 10 nm 10 – 1 nm < 1 nm w=20 nm w=1.5 nm E=1200 E=2000   w=1 𝜇 m E=300  w=2 nm w=70 nm E=1200 E=1000   w=10 𝜇 m E=30 w=9 nm  E=180 w=100 nm  E=1000  Decreasing Antenna Gap Width (w) Figure 1. THz-ﬁeld enhancement proﬁle with decreasing gap width. is a non-invasive and one of the most widely used 80,83–87 technique in biosensing because the THz-spectra are less energetic (about a few milli-electronvolts) lies in the characteristic vibration frequency regime of pro- teins and DNA molecules. Terahertz nanoantennas provide us label-free and ultra-sensitive sensing of 88,89 37,81,90–94 molecules and biomolecules as the nanoan- tennas focus the THz waves into an extremely con- ﬁned volume where we can place our analytes at very low concentrations. Signiﬁcant changes in the mag- nitude of the transmission spectrum were observed in the absorptive molecules when the resonance mode of 36,83 the waveguide matches its absorption spectrum, en- hancing the sensitivity of THz-resonators based sens- ing devices, therefore making it feasible for detection of biomolecules at extremely low concentration. Figure 2. Photolithography : A high resolution and high In this review work, we describe the high-throughput throughput technique for nanogap structures. Reprinted with fabrication techniques of THz nanoantenna via. pho- permission from Ref. . © 2018, De Gruyter. tolithography, atomic layer lithography, pattern and peel method, self-assembly lithography etc. Then we discuss the THz ﬁeld enhancement using metal nanogap have taken real examples from various literature that antenna where we discuss the ﬁeld enhancement be- how THz nanoantenna and metamaterials can be used yond the skin depth of the metal up to the quantum to enhance the terahertz biomolecular sensing. regime. Then we brieﬂy discuss ﬁnite-difference time- domain (FDTD) methods for calculation of ﬁeld en- hancement across the nanogap. At the last, we de- II. HIGH-THROUGHPUT WAFER-SCALE scribe the THz-TDS for sensing the biomolecules due FABRICATION TECHNIQUES to the ﬁeld enhancement across the nanogap. We also showed how a molecular dynamics (MD) simulation can be a versatile tool to estimate the terahertz absorp- A nanogap must have a reliable THz-ﬁeld interaction tion and vibrational density of states (VDOS). Then we and plasmonic characteristics (like, ﬁeld enhancement THz-Field Enhancement (E) Quantum Regime 3 Pattern and ALD Plug in second layer A B E G Adhesive C Apply adhesive D Figure 3. Atomic layer lithography. (A) Thin Al O layer coated on a patterened substrate using ALD. (B) Metal evaporation 2 3 plugging second metal layer. (C,D) Removing excess metal by using an adhesive. Scanning electron micrograph (SEM) of (E) 5 nm gap in Ag-ﬁlm (F) 10 nm gap between Au and Ag, and (G) 9.9 A Au - Al O -Au vertical nanogap. Reprinted with permission 2 3 95 21 from Ref. . © 2013, Nature Communications. (H) Adhesive tape planarization of nanogap. Adapted with permission from Ref. . © 2014, Springer Nature. (I) Planarization of half Si-wafer by using an adhesive. Reprinted with permission from Ref. . © 2013, Nature Communications. and non-linearity). But fabricating high throughput Polymers added to PR decides the solubility under UV wafer-scale nanogaps is one of the most challenging exposure. A PR which consist of an insoluble polymer part of the job. To achieve sufﬁciently intense far-ﬁeld which on UV exposure, changes into a soluble polymer, signals, nanogap arrays should be fabricated over is called a positive PR and a PR which consist of a sol- large sample area. Several methods like: fs-laser uble polymer which on UV exposure, polymerizes and 17,96–98 17,40,99–104 becomes insoluble, is called a negative PR. Therefore by beam machining, focussed ion beam, 17,105–108 using the desired PR, wafer-scale metal structures are photolithography, electron-beam (e-beam) 17,109–115 17,116 fabricated after performing a metal deposition process lithography, nanosphere lithography, on patterned resists and lift-off of the unwanted metal 117 17,118–122 sidewall lithography, nanoimprint, ﬁlm. pattern and peel method and atomic layer 17,41,95,124–126 The lithography method discussed earlier are lim- lithography have been developed for ited to fabrication of nanometer-gaps up to a few high resolution patterning over large area of sample. nanometers. To study the nano- and sub-nanoscale op- 17,23 tics, we need better lithography techniques to fabri- Photolithography (Fig. 2) is one of the most high- cate sub-10 nm wide gap. Many researchers have re- throughput methods for fabricating nanogap structures. ported the signiﬁcance of deposition of an atomically In photolithography, a desired polymer pattern is pat- thin insulating layer or sacriﬁcial layer in fabricating terned using ultraviolet (UV) light and a photomask. 127–132 small gaps between two metal structures. Chen This polymer pattering is formed by the physical change of a photo-sensitive material, known as photoresist (PR). et al. developed a pioneering method for fabricat- 4 Cr/Au Film Cr/Al Pattern Ion Milling on a clean and dried sapphire substrate. With the help of a thermal evaporator (evaporation rate = 1 A/sec), a chrome layer of thickness 10 nm was deposited over the PR pattern. Following the process, a silver ﬁlm of thickness 200 nm was evaporated over the chrome layer. After the evaporation of ﬁrst silver layer, the whole sam- ple was then dipped in acetone sonicated for 2-minutes at 150 W and 36 kHz, lifting off the PR pattern. After sonication, the sample was washed with isopropyl alco- Al/Cr Wet Etching Cr/Au Deposition Al O Coating 2 3 hol and dried using N gas. The dried sample was then dipped into a toluene solution of OT (1-8 octanedithiol) Figure 4. High-throughput fabrication technique of functionalized QDs, resulting in the formation of a self- nanometer-scale gaps introduced by Jeyong et al. Adapted assembled monolayer of QDs. The resulting sample with permission from Ref. . © 2014, Springer Nature. was again dried using N gas and a second silver layer of same/higher thickness was deposited over the self- assembled monolayer. Then by using a scotch tape , ing vertical nanogaps known as atomic layer lithogra- the second layer of silver was taped-off, leaving a verti- phy (Fig. 3 (A - D)). Figures 3 (E - G) are the scan- cal metal nanogap ﬁlled with a monolayer of QDs. ning electron micrographs (SEM) of 5 nm gap in Ag ﬁlm, 10 nm gap between Au and Ag, and 9.9 A Au - Al O - Au vertical nanogap, respectively, which are 2 3 fabricated using atomic layer lithography technique. III. TERAHERTZ-FIELD ENHANCEMENT USING This method also includes a new planarization scheme NANOANTENNAS which eliminates the background light transmission, enabling background-free transmission measurements. This lithography technique is the combination of atomic In THz frequency, a metal ﬁlm is not a per- layer deposition (ALD) and ‘plug-and-peel’ adhesive fect conductor as well as inﬁnitely ﬂat. Equiv- tape planarization and is a three-step process. The ﬁrst alence of THz-electric ﬁeld enhancement in metal 18,28,34,35,134–136 step includes the pre-patterning of metal using con- nanogap and magnetic ﬁeld enhance- 137–140 ventional lithography techniques, the second step in- ment around metallic nanoantennas can be 141–145 cludes the sidewall formation of materials like Al O or 2 3 shown by using Babinet's principle. This prin- SiO by ALD, and the third step includes the ‘plug-and- ciple is still used implicitly in different areas of peel’ taping process for planarization of metal plugs. 22,27,36,38,42,93,103,107,146–168 THz research. Sommerfeld's Therefore, atomic layer lithography is a high-resolution 169–171 half-plane problem is also used to study the patterning method with large scale uniformity. light-aperture interaction theoretically. Bethe and 173,174 Jeong et al. improvised the atomic layer lithogra- Bouwkamp reconsidered the Sommerfeld's half- phy technique which is demonstrated in Fig. 4. In this plane problem for small apertures in inﬁnitely thin lithography technique, conventional photolithography perfect conducting plates. Further studies were per- and liftoff process were used to pattern 30 nm Cr/150 formed on the light interaction with system of periodic nm Al layer on the 3 nm Cr/100 nm gold-sapphire sub- 175 apertures. Later, extensive studies were performed strate. The Cr-Al double layer acts as a sacriﬁcial layer on enhanced light transmission through the apertures which will be etched out later, removing excess metal 19,31,32,102,176–180 in various spectral regime. Theoretical and making the structure planar. As compared to Al, 135,181–188 approaches like: coupling-mode theory, mi- Au is less-resistant to ion beam therefore, excluding the 28,29,189 190,191 croscopic models, transfer-matrix method Au ﬁlm beneath Al-Cr layer, the exposed Au layer is and capacitor model, were combined with numerical removed using ion milling process. Uniform ALD of 76,192,193 methods to solve the enhanced light transmis- Al O is performed over the whole structure with a 2 3 sion and huge ﬁeld enhancement around metallic aper- thickness of nanometer accuracy. After that, the sec- tures. ond layer of Au with adhesive Cr is deposited, ﬁlling the trenches. Finally, Al/Cr wet etching is performed When an electromagnetic wave is incident normally which removes the overhanging Au and Al layers, ex- on a perfectly conducting metallic plane, an induced posing the dielectric gaps. This lithography technique current is developed on the surface, which reﬂects light can be used to fabricate large-scale nanogaps with ultra- back, with no charge accumulation anywhere. When high aspect ratio. this plane is divided into two perfectly conducting Sommerfeld half planes, macroscopic accumulation of Tripathi et al. developed self-assembled lithogra- charges take place at the edges with a length scale of one phy process by which they fabricated quantum dots wavelength, such that the surface charge density (s) de- (QDs) nanogap metamaterials (Fig. 5). In this lithogra- phy process, photolithography was used to pattern PR pends as a time function for small values of x l, 5 Sapphire Substrate Photolithography A B st Chrome Layer Deposition Silver 1 Layer Deposition C D Photoresist Chrome Silver Lift-off of Photoresist Self-assembly of Quantum Dots E F nd nd G Silver 2 Layer Deposition H Adhesive Taping of Silver 2 Layer CdSe Quantum Dots Mean size including capping molecule is 8 nm Figure 5. Self-assembly lithography technique. (A) Clean and dried sapphire substrate. (B) PR patterned using photolithography. (C) 10 m chrome layer deposition over PR pattern. (D) Evaporation of silver ﬁlm of thickness 200 nm over the chrome layer. (E) Sonication and lift-off of PR after ﬁrst silver layer deposition. (F) Self-assembly of QD monolayer over the ﬁrst layer of metal. (G) Second silver layer deposition after self-assembly of QD monolayer. (H) QDs ﬁlled nanogap metamaterial after taping-off second silver layer. Adapted with permission from Ref. . © 2015, Optical Society of America. given by and also in the quantum regime. We also discuss the classical limit of ﬁeld enhancement of single nanoslit e l 0 before entering quantum regime, and about ﬁnite- i(wt+p/4) s(x, t) = p E e (1) difference time-domain algorithm, which is a computa- 2p tional approach developed for the quantitative estima- Here, e , E , w and x describes the permittivity of vac- 0 0 tion of far- and near-ﬁeld enhancement measurements. uum, the incident electric ﬁeld, the angular frequency and distance from the edge, respectively. The charge singularity, at x = 0, for this half plane is very feeble and disappears with integration. When the two metallic half A. THz Field Enhancement in Metal Nanogaps at Sub planes are brought close together, an electrostatic force Skin-Depth Regime will be experienced by the charges due to their counter members across the gap, pulling the charges towards In this part of the section, we discuss the skin- the edge, and developing a strong electric ﬁeld across depth physics and the results reported in various the gap. As the gap shrinks, more charges gets accu- 16,34 literatures . Consider two perfectly conducting mulated at the edges as the light-induced currents be- 16,134 comes more and more stronger. This increases the sur- metallic planes of thickness h kept at a distance w face charge density at the edges as shown in Fig.6 B - C. from each other, forming a metal-air-metal nanogap. An This system can be portrayed as a line-capacitor driven electromagnetic wave is incident normally on the metal by the light-induced alternating current as shown in nanogap whose electric ﬁeld has a magnitude of E . Fig. 6 (A). Then the ultimate ﬁeld enhancement in a high-aspect ra- Numerous studies have reported the monotonic en- tio (w/h 1) is given by hancement of the electric ﬁeld with decreasing gap 34,93,95,145 size. Further studies have also reported the E l E = = (2) enhancement decrease in ﬁeld enhancement when the gap size de- E ph 24,194–199 creases and enters a quantum regime. Increase in ﬁeld enhancement beyond sub-skin depth is also where, E denotes the electric ﬁeld at the gap and l de- 16,23,34 reported. In this section of the review, we discuss notes the wavelength of the incident electromagnetic the THz ﬁeld enhancement for gaps beyond skin-depth wave in vacuum. If a dielectric material of permittivity 6 THz (E ) H = 2H . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . J ≈ 2H /h . . . . . . . . H = tH ; |t| << 1 B + + - C + - + - + - 𝛿 h + - J ≈ 2H /h + - + + - - + - 𝘀 H = 𝘀 E = 𝘀 E E = E 𝑚 𝑚 0 m t 0 𝜋ℎ + - + - + - + - Figure 7. Field enhancement in metal ﬁlm of thickness below + - skin-depth. (A) THz ﬁeld incident normally on a metal ﬁlm of thickness less than THz skin-depth. (B) Magnetic ﬁeld of light Figure 6. Capacitor-model representation of a nanogap. (A) (represented by blue encircled dots) and current density (rep- Nanogap gets charged just like a line-capacitor due to light- resented by red lines) near the thin metal ﬁlm. At incident sur- induced alternating current J. Hence, the electric ﬁeld en- face (up), the magnetic ﬁeld is twice that of the incident ﬁeld hances which is shown by gradual colour contour. Reprinted and a transmission side (down), the magnetic ﬁeld is much with permission from Ref. . © 2009, Springer Nature. smaller as compared to the incident side. (C) Displacement Schematics of charge accumulation near metal edges and cor- current and E around metal gap. Adapted with permission responding current density induced inside metal ﬁlm of two from Ref. . © 2018, De Gruyter. gap-sizes (B) w h( l) and (C) w h( l). Adapted with permission from Ref. . © 2017, American Physical Soci- ety. coefﬁcient (r) for such thin ﬁlms are given by E 2e c h h 0 0 0 t = = 1; r 1 1 (4) e is ﬁlled inside the gap, Eq. 2 modiﬁes to E s h h h 0 m E l e where, E denotes the transmitted electric ﬁeld. E = = (3) enhancement E ph e Considering a transverse magnetic polarized light (E , H ) is incident normally on a thin metal ﬁlm. Due x y where, e denotes the permittivity of vacuum. From the to reﬂection from the incident surface, the magnetic ﬁeld above equations, it can be said that the ﬁeld enhance- of the light near the incident surface becomes twice that ment (E ) is independent of metal characteris- enhancement of the incident magnetic ﬁeld. But the magnetic ﬁeld on tics. the transmission side is much smaller than the incident Considering an electromagnetic wave incident nor- ﬁeld (Fig. 7 (B)). Assuming that a constant electric ﬁeld- mally on a thin metal ﬁlm of thickness h which is smaller q /current density is developed inside the thin ﬁlm, us- than skin-depth d = , but larger the characteris- ing Ampere's law, the expression of the current density m s w 0 m 2H 2e c 0 inside the thin ﬁlm is J = (neglecting vacuum dis- tic thickness h = ; (h h d), shown in Fig. 7 h 0 0 placement current). Here, J denotes the current den- (A). Here, s denotes the conductivity of the metal, m sity developed inside the thin metal ﬁlm and H denotes denotes the magnetic permeability in vacuum and w de- the incident magnetic ﬁeld. The tangential component notes the angular frequency of electromagnetic wave. of electric ﬁelds at the air-metal interface are continuous The thickness at which the absorption loss by metal is and is given by 50%, is known as the characteristic thickness (h ). For a 7 1 1 metal with conductivity 10 W m , the characteristic J 2H 2e c 0 0 E = E = = = E (5) thickness is 0.53 nm. At 1 THz, the skin-depth is about t m 0 s s h s h m m m 100 nm or more for good metals. Therefore, the thick- ness range for transitional metal ﬁlm is 5-100 nm. The therefore, reproducing Eq. 4. At transmitting side, elec- electric ﬁeld amplitude transmission (t) and reﬂection tric ﬁeld just inside the metal surface is denoted as E . m 7 Figure 9. Electron tunnelling through thick and thin barriers. Due to the high gap electric ﬁeld E, separation of Fermi level h in two metals by potential energy V = eEw takes place. f represents the rectangular barrier height and the dotted line describes the barrier including image forces. Reprinted with permission from Ref. . © 2015, American Physical Society. Through this analogy, it can be said that with an in- crease in the conductivity of the metal, weaker electric ﬁelds are induced inside the metal. This is compen- sated by the metal's high dielectric constant when the Figure 8. THz ﬁeld enhancement in metal nanogaps at sub- displacement boundary condition is considered. There- skin depth regime. (A) SEM of a 70 nm wide gold nanogap fore, the ﬁeld enhancement is independent of the con- fabricated using focused ion beam. (B) Horizontal electric ﬁeld ductivity of the metal, which is mathematically proven enhancement in a 70 nm gap, obtained from FDTD simulation. above. Fig. 8 (A) is the SEM of a gold nanogap sample Reprinted with permission from Ref. . © 2009, Springer Na- ture. (C) Variation of electric ﬁeld enhancement with respect of gap width 70 nm, fabricated by Seo et al. using a fo- to gap thickness, indicating the 1/h dependence of the ﬁeld cused ion beam technique. Fig. 8 (B) shows the x-z plane enhancement. Adapted with permission from Ref. . © 2017, horizontal electric ﬁeld enhancement in 70 nm gap, ana- American Physical Society. lyzed from FDTD simulation at 0.1 THz. In the ﬁgure, it can also be seen that the electric ﬁeld does not penetrate into the metal, rather it is seen to be completely concen- Across the air nanogap, the normal component of the trated inside the nanogap, even if the gap size is much displacement current will be same (Fig. 7 (C)), therefore smaller than the skin depth (250 nm). Fig. 8 (C) shows electric ﬁeld inside the air nanogap will be the varying electric ﬁeld enhancement of 5 nm gap with respect to the gap thickness (h) at a frequency of 0.3 THz. e s 2e c l m m E = E E = E (6) m 0 0 The red dots signify the experimental data and the red e e w s h ph 0 0 m solid line shows the calculated data from the modal ex- where, e denotes the terahertz metal dielectric constant pansion. the ﬁgure also signiﬁes the 1/h dependence of given by the ﬁeld enhancement. s s m m e = e + i i (7) m ¥ w w w 1 i B. THz Field Enhancement in Sub-Nanometer and Angstrom Metallic gaps: Quantum Regime where, e denotes the dielectric constant of metal at high-frequency and g denotes the Drude damping con- Recently, Bahk et al. discussed about the ﬁeld en- stant. hancement in thin (w h < d) and thick (w < d < h) 8 current of THz-frequency. FDTD analysis of the elec- tric ﬁeld distribution around 1.5 nm wide gap in 150 nm thick gold ﬁlm at a frequency of 1.5 THz, as shown in Fig. 10 (A). As expected the horizontal electric ﬁeld is strongly concentrated and is enhanced by a factor of 2000 inside the gap. Fig. 10 (B) shows the electric ﬁeld enhancement in different gap size for various frequen- cies. The theoretical calculation of the modal expansion for the perfect electric conductor (PEC) model is shown by solid lines. C. The Finite-Difference Time-Domain (FDTD) Algorithm : A Computational Electromagnetism Approach The ﬁnite difference time domain (FDTD) method is used to solve Maxwell's equations in the time domain and generally used for computational electrodynamics. ~ ~ r D = r (8) ~ ~ r B = 0 (9) Figure 10. Classical limit of THz ﬁeld enhancement before en- tering quantum regime inside gold nanogap. (A) FDTD sim- ulated horizontal electric ﬁeld around gold nanogap of width ¶B ~ ~ r E = (10) 1.5 nm. (B) Variation of ﬁeld enhancement in different gap ¶t widths for various THz frequencies. Solid lines represents the theoretically calculated values from modal expansion formal- ism. Adapted with permission from Ref. . © 2017, American ¶D ~ ~ ~ Physical Society. r H = J + (11) ¶t are the four Maxwell's equations where, the electric narrow gaps, and in wide (d, h < w) gaps. In wide gaps, ﬁeld, magnetic ﬁeld, electric displacement ﬁeld, mag- the charges are mostly spread over the surface outside netic ﬂux density, free electric charge density, and free the gap rather than at metal edges (Fig. 6 (B)), which re- current density are denoted by E, H, D, B, r and J, 212,213 sults in a decrease of ﬁeld enhancement. In narrow gaps, respectively. These equations are solved numer- charges mostly accumulate at the metal edges of the gap ically on a discrete grid in both space and time, and (Fig. 6 (C)), resulting in stronger ﬁeld enhancement. As derivatives are handled with ﬁnite differences. No ap- the width of the gap (w) decreases further entering the proximations or assumptions are made about the sys- sub-nanometer and Angstrom regime, the charge dis- tem, making this method highly versatile and accurate. tribution becomes insensitive to gap size and electron It is a fully vectorial simulation method as it solves for tunnels through the potential barrier of the nanogap, all electric and magnetic ﬁeld vector components. Be- 197,199–210 showing quantum effects. When THz electro- ing a time-domain method, FDTD can be used to calcu- magnetic waves are incident on the nanogap, a transient late broadband results from a single simulation. FDTD voltage is developed in the dielectric gap which bends method is typically used when the feature size is of the the conduction band of the dielectric toward the Fermi order of the wavelength. It is general, versatile, accu- energy of metals (Fig. 9). This increases the chances for rate, broadband and fasts which makes it the most reli- electron tunneling through the potential barrier, causing able method in computational electrodynamics. non-linear transmissions. FDTD algorithm solves Maxwell's curl equations in Bahk et al. also investigated the classical limit of ﬁeld non-magnetic materials. enhancement before the gap-width entered the quan- ¶D tum regime. They reported that the ﬁeld enhancement = r H (12) exhibited saturation behavior rather than a monotonous ¶t increasing nature. This showed that the nanogap acts 34 ~ ~ ~ where, D = eE. Since the region is vacuum, J = 0. like a charged capacitor which is shown in the Fig. 6 (A), whose total induced charge is inversely propor- ~ ~ tional to the frequency of a light-induced alternating D(w) = e e (w)E(w) (13) 0 r 9 A B Figure 11. FDTD analysis of a gold nanogap and nanorod. (A) Horizontal electric ﬁeld around 100 nm gold nanogap. (B) Simulated magnetic ﬁeld around a gold nanorod. Yellow arrows represent the current ﬂow. Reprinted with permission from Ref. . © 2009, American Physical Society. ¶H 1 and for TE case, Maxwell's equations reduces to = r E (14) ¶t m ¶E ¶H 1 ¶E z y x = (22) with H, E, and D describing the magnetic ﬁeld, electric ¶t m ¶x ¶y ﬁeld and displacement ﬁeld, respectively, while e (w) is a complex relative dielectric constant given by D (w) = e e (w)E (w) (23) x r x e (w) = n (15) D (w) = e e (w)E (w) (24) y 0 r y , where n denotes refractive index of the material. In 3-D, the six electromagnetic components of Maxwell's equations are - E , E , E and H , H , H . Assuming that x y z x y z ¶D ¶H x z = (25) in z-dimension, the structure is inﬁnite and the ﬁelds are ¶t ¶y independent of z, such that e (w, x, y, z) = e (w, x, y) , and (16) r r ¶D y ¶H = (26) ¶t ¶x ~ ~ ¶E ¶H These equations are solved on a discrete, spatial and = = 0 (17) ¶z ¶z temporal grid. Each ﬁeld component is solved at a slightly different location within the grid cell (or, Yee then, two independent sets of Maxwell's equations will cell). be created. Each set will contain three vector quantities With this method we can simulate the electric and which can only be solved in x-y plane. The equation magnetic ﬁelds around any nanostructures using com- containing the components - E , E , H , is known as x y z mercially available softwares like, FDTD Solutions transverse electric (TE) equation, and the equation con- (Lumerical Inc., Canada), COMSOL, etc. Fig. 11 shows taining the components - H , H , E , is known as trans- x y z the FDTD analysis of gold nanogap (100 nm thick and verse magnetic (TM) equation. Therefore, in TM case, wide gap) and nanorod, respectively. Maxwell's equations reduces to ¶H ¶D ¶H z y x = (18) IV. THZ-TDS AND SENSING OF MOLECULES AND ¶t ¶x ¶y BIOMOLECULES 214–216 D (w) = e e (w)E (w) (19) z r z 0 In the late 1980s, Grischkowsky et al. intro- duced the terahertz time-domain spectroscopy tech- nique (Fig. 13). It is a spectroscopic technique in which ¶H 1 ¶E x z = (20) short-THz pulses are used to probe the properties of ¶t m ¶y matter. The system generally consisted of the femtosec- ond laser (fs-laser), which operates at a repetition rate ¶H y 1 ¶E z of 100 MHz, producing a 100 fs laser pulse train. Us- = (21) ing a beam splitter, the fs-pulse train is split into two ¶t m ¶x 0 10 A B C Figure 12. Spectral analysis for identiﬁcation of biomaterials. (A) Distinct spectral signatures can be used to discriminate a bacterial cell, spore and intracellular metabolite. (B) Bacterial species can be differentiated from their differences in THz opti- cal constants. (C) Resonant frequency shift analysis using THz nanoantenna biosensors can be used to study the presence of biomaterials. Adapted with permission from Ref. . © 2016, Springer Nature. 92–94 and metamaterials have proved to be effective de- vices for detection of micro-organisms, organic- and biomolecules. These devices show resonances with strong ﬁeld enhancement across the gap in the THz fre- quency range, which is highly sensitive to changes in dielectric constant of the gap region. Therefore, these devices are used as biosensors for ultra-sensitive and label-free detection of micro-organisms, organic- and biomolecules. In this section, we discuss the two widely used label-free THz-TDS methods : (i) THz-TDS using molecular dynamic (MD) simulations, and (ii) THz sens- ing by nanoantennas and metamaterials. Figure 13. Conventional THz time-domain spectroscopy setup developed by Grischkowski et al. Adapted with permission from Ref. . © 2016, Springer Nature. A. Computational THz-TDS using Molecular Dynamics (MD) Simulation beams: a pump beam and a probe beam. The pump beam is made to be incident on a THz emitter to emit MD simulation is a computational method used to THz pulses, which is collimated to the sample using a study the dynamics of a system of molecules in the con- pair of parabolic mirrors. Concurrently, the probe beam densed phase. To accurately study the dynamics of a is used in a time-gated manner for the detection of THz complex molecule, one needs to employ the quantum electric ﬁeld which contains time-domain information of model to study the wave function of each sub-atomic phase and amplitude. Coherently, the transmitted THz particle. But MD employs classical Newtonian mechan- electric ﬁeld is measured as a function of time to obtain ics to study the dynamics of the system, which makes it a time-domain signal which is converted to noise-free less accurate. Speciﬁcally, it integrates Newton's equa- frequency-domain signal using Fourier transformation tion of motion in discrete time-steps. Due to high THz function. absorption by water, THz-TDS becomes challenging for THz-TDS is a non-contact, non-invasive and label- water-based systems. Therefore, MD simulation is one free detection technique and is extensively applied in of the preferred methods to study the water dynamics biomedical imaging. However, THz-TDS has proved of a hydrated molecular system. In MD, water dynamics to be less sensitive for detection of micro-organisms has a great signiﬁcance in the structural arrangements of (like yeast, molds, and bacteria) of size equal to or molecules, which occurs because the molecule alters the less than l/100, due to their transparency under THz dynamics of the surrounding water molecules, adopting 88,89,217 frequencies. Recently, plasmonic nanoantennas a quasi-coherent character, caused by reorganized, loose 11 Figure 14. Flowchart of a general MD simulation process. Figure 15. THz time-domain analysis obtained using MD hydrogen-bond network. simulations. (A) Velocity autocorrelation fuction (VACF) and A signiﬁcant number of MD-based computational (B) vibrational density of states of oxygen (ox) and hydrogen analysis have been performed to interpret the THz (h) in ﬁrst solvation layer (3 A) of lysosome and bulk water. absorption spectra of protein solutions. Many stud- Reprinted with permission from Ref. . © 2013, American Chemical Society. (C) Absorption spectrum of water obtained ies showed the anomalous dynamics of water in the experimentally (black line) and from MD simulation (red line). hydration shell due to the fact that surrounding wa- Adapted with permission from Ref. . © 2014, American ter molecules were heterogeneously perturbed by the 220–224 Chemical Society. solute. Bandyopadhyay et al. have performed a signiﬁcant amount of MD simulations to study the ef- 223,225–232 fects of hydration of protein molecules. Var- ter dynamics, which were limited to the ﬁrst hydration ious methods have been used to study the thickness layer. Some similar results were obtained in which of villin headpiece sub-domain HP-36 which includes it was demonstrated a clear distinct power spectrum MSD, velocity autocorrelation function analysis, H- (VDOS) for water molecules bonded to different planes bond time correlation function analysis and reorienta- of antifreeze protein by hydrogen in the spectral do- tional correlation function analysis of water molecules. According to the simulation results, the different heli- main of 1-4 THz. While studying the vibrational spec- cal segments had a heterogeneous inﬂuence on the wa- trum of water in the villin headpiece subdomain hy- 12 dration layer, OOO bending mode showed a blue- shift in the ﬁrst hydration layer, which is noticed for water molecules bonded to protein by hydrogen. It also showed the possibility of the structural ﬂexibility of protein. Later an atomistic MD simulation of hen egg-white lysosome solvated in explicit water molecules at room temperature was performed, which reported that a few large-amplitude bistable motions exhibited by two coils controlled the overall ﬂexibility of protein molecules. A series of MD simulations were performed to compare the experimental data of far-infrared spec- troscopy used to study the dynamics of three aqueous peptides with varied helicity. Using the ﬁrst princi- ples, extensive studies on water, resolved in time and space, reported that the group motion of H-bonded molecules in the second solvation shell signiﬁcantly contributed to the absorption at about 2.4 THz, also showing the presence of third-shell effects. Heyden et al. in their future publication showed the capabil- ity of nonpolarizable water models in reproducing low- frequency, inter-molecular vibration of water since elec- tronic polarization is dominated by the static molecular dipoles. The depth of hydration shell in the hydrated lysosome, BPTI, TRP-cage, and TRP-tail was estimated using MD simulations. To perform a basic MD simulation of a molecule, ﬁrst and foremost a force ﬁeld must be applied to char- acterize the molecular interaction of the molecules. In Figure 16. Application of THz nanoantenna in developing a the next step, the molecule is positioned in a unit cell of highly-sensitive THz-TDS technique. desired shape and size followed by the addition of wa- ter molecules, known as solvation step. Then ions are added to neutralize the net charge of the system. After ulation steps are summed up in a ﬂowchart (Fig. 14). the addition of ions, the system is then relaxed through a process known as energy minimization. This process Fig. 15 (A) shows the MD simulation generated VACF also ensures that the system is free from steric clashes curves for oxygen and hydrogen atoms in bulk water or inappropriate geometry. For further simulation, the and in the water within 3 A of lysosome molecule (ﬁrst system needs to be brought to the desired temperature. hydration layer). The ﬁgure depicts two results: i) The After attaining the desired temperature, an adequate hydrogen dynamics get uncorrelated faster (0.85 ps ear- amount of pressure must be applied until it reaches a lier) than oxygen dynamics, and ii) The dynamics of proper density. This whole process is known as equili- the oxygen atoms are much more restricted in the ﬁrst bration and is conducted in two phases. The ﬁrst phase hydration layer than in bulk water (as compared from is conducted to attain the desired temperature under their extremums) due to caging-effect. The correspond- NVT (constant Number of particles, Volume, and Tem- ing VACF data was further Fourier transformed to get a perature) or isothermal-isochoric ensemble. The second set of VDOS data as shown in Fig. 15 (B). In the ﬁgure, phase is conducted to stabilize both pressure and den- at 1.1 THz a well-deﬁned peak is shown by the oxygen sity of the system under NPT (constant Number of par- atoms of bulk water signifying the bending motion of ticles, Pressure, and Temperature) or isothermal-isobaric O – O – O atoms (known as triplets of H-bonded oxy- ensemble. After attaining the desired temperature and gen atoms). But in the case of lysosome-bounded water pressure, the system is ﬁnally ready for the MD simu- molecules, this peak is blue-shifted by 0.4 THz signify- lation, known as MD production run. After performing ing that the H-bonds between O atoms gets stronger in the ﬁnal MD-run, the system is ready to be analyzed. the hydration layer. This peak has a lower amplitude From MD simulation, the time-dependent properties are compared to that of the oxygen atoms of bulk water, as obtained via. correlation function, like velocity autocor- the presence of lysosome does not have much inﬂuence relation function (VACF). The time-domain properties on the vibrational mode of water molecules. The VDOS are then converted to frequency-domain properties us- of hydrogen in both situations is very ﬂat as they are not ing Fourier transformation function, for example, the vi- much inﬂuenced by the protein molecule. The absorp- brational density of states (VDOS). The whole MD sim- tion spectrum can be solved by applying Fourier trans- 13 Figure 18. Effects of various gap-widths on THz-sensing of Figure 17. Interpretation of nanoantenna enhanced THz inter- RDX molecules. (A) Schematic of the cross section of a THz actions of biomaterials using FDTD method. (A) Schematic of single-slot nanoantenna containing RDX molecules. (B) Ab- the cross section view of a metallic nanogap cladded with a sorption spectrum of RDX molecules (1 mg) onto the bare absorptive homogeneous dielectric ﬁlm with complex refrac- quartz substrate. (C) Absorption coefﬁcients peaks and cross tive index. (B) Normalized transmittance simulated by vary- section of RDX molecules in various nanoslot gap-widths ing the imaginary part of the complex refractive index of the w. The black dashed line signiﬁes the 1/w dependence. homogeneous dielectric ﬁlm. (C) Transmittance peak values Reprinted with permission from Ref. . © 2013, American and resonance frequencies are plotted for various complex re- Chemical Society. fractive indices. Adapted with permission from Ref. . © 2018, De Gruyter. Figure 15 (C) is the absorption spectra of bulk water ob- tained experimentally (black line) and by MD simula- form on total dipole moment autocorrelation function of tion (red line), describing the inter-molecular hydro- the system computed from MD simulation, given by gen bond collective modes. 1 2pw iwt a(w) = e M(0)M(t) dt (27) 4pe 3Vk Tcn(w) 0 B B. THz-Nanoantenna and Metamaterials Sensing of Micro-organisms, Organic- and Biomolecules here, V denotes the volume of the system, k denotes the Boltzmann's constant, T denotes the absolute tem- perature of the system, c denotes the speed of light and THz nano-structures can be used for ultra-sensitive, n(w) denotes the frequency dependent refractive index. label-free and non-invasive sensing of molecules and 14 B C D E F Figure 19. Effects of different substrates on yeast detection. (A) Schematic of THz sensing of yeast cells using a THz slot antenna. (B) SEM image of yeast cells deposited on 2 mm wide and 100 mm long slot antenna. (C) THz transmission amplitudes for THz slot antenna on Si substrate showing a red-shift in the resonant frequency of the slot antenna due to the presence of yeast cells. THz transmission amplitudes for THz slot antenna on (D) Si- and (E) quartz-substrates at different yeast cell concentration (N ). av (F) Shifts in resonance frequency of THz slot antenna on Si- and quartz-substrates with respect to N . Reprinted with permission av from Ref. . © 2014, Springer Nature. biomaterials. As discussed earlier in this section, the refractive index is given by strong THz ﬁeld enhancement along with resonances in (f f )c the gap region of these nano-structures are sensitive to r s n(w) = 1 + (29) the changes in the dielectric constant of the insulating 2pd f gap material. If the gap is ﬁlled with desired bioma- terial sample, due to enhanced THz-interaction, the di- and applying Beer-Lambert law, absorption coefﬁcient electric (i.e., the biomaterial sample) will show an en- is given by hanced THz-absorption therefore, the THz optical prop- erties can be studied from the THz dielectric response. 2 2 A (w) 4p f k a(w) = ln(T) = ln = Thereby, making the necessary change in the THz-TDS d d A (w) c setup (Fig. 16), enables enhanced sensing of molecules (30) and micro-organisms which were not sensed by con- where, f and f are the phases of sample and reference r s ventional THz-TDS setup. Therefore, from the Fourier transmitted signals, respectively, T denotes the transmit- transformed spectra (frequency-domain amplitude and tance and f is the frequency of the input THz-signal. phase) obtained from the sample and reference trans- Through FDTD simulation, an intense study of this mitted signals, the complex optical properties ae com- biosensing process can be done. The biomaterial sam- puted from the following relation: ple ﬁlm is assumed to be a homogeneous dielectric ﬁlm of a thickness proportional to the molecular concentra- tion, is sandwiched between two gold ﬁlms (Fig. 17 (A)). da(w) i n(w)d 2 c A (w) = A (w) e e (28) s r The dielectric ﬁlm has a complex refractive index An + i Bk, where A and B are constants whose values range from 1.0 - 3.0 over a frequency band. The absorptive where, A and A are the amplitudes of sample and ref- dielectric ﬁlm was inserted in the gold nanogap by ap- s r erence transmitted signals, respectively. The real parts plying an auxiliary differential equation. To study the of the absorption coefﬁcient and the index of refraction electric ﬁeld interaction in the absorptive dielectric ﬁlm, are denoted by a(w) and n(w), respectively, and thick- a non-uniform mesh was applied over the whole ﬁlm ness of the dielectric is denoted by d. The real part of with the smallest step-size of 10 nm. To the nanoan- 15 A B C D E F Figure 20. Highly-sensitive and selective nanoantenna. (A) Schematic of a THz nanoantenna for detection of sugar molecules. (B) Fructose antenna with stains of 250 mg/dL fructose (top-right) and D-glucose (bottom-left). (C) THz transmittance through stained fructose antenna illustrating the selective nature of nanoantenna. (D) THz transmittance peaks of fructose and D-glucose at different concentrations using fructose antenna. (E) THz transmittances of added sugars and artiﬁcial sweeteners in various branded beverages and dietary sodas using a fructose antenna. (F) THz absorption coefﬁcients of low-concentrated artiﬁcial sweeteners: Acesulfame K and Aspartame. Reprinted with permission from Ref. . © 2015, Springer Nature. tenna, THz electromagnetic waves are made to be in- sensing, a single-slot THz nanoantennas of length 90 cident normally. As the THz waves pass through the mm and slot-widths of 50 nm (50 nm thick), 100, 200, dielectric ﬁlm, the transmitted THz-electric ﬁeld decays 500, 1000 and 5000 nm (100 nm thick) was fabricated on exponentially and is given by a thick quartz substrate and ﬁlled with 1 mg/ml RDX (1,3,5-trinitroperhydro-1,3,5-triazine) molecules. Ini- T (w) kk(w)h tially, the THz-absorption (a) of 1 mg RDX molecules = Ce (31) T (w) r placed over a bare quartz substrate was studied which is shown in Fig. 18 (B). Then the absorption of RDX 2 2 here, T = (E (w)) and T (w) = (E ) are the trans- s s r r molecules inside the various nano-slots was studied and mittances through dielectric and air nanogaps, respec- the data obtained was plotted against the slot-widths as tively. In the above relation, C is known as the trans- shown in the Fig. 18 (C). From Fig. 18 (B - C), it is evi- 2p f mittance ratio at the air-dielectric interface, k = is dent that nanoantenna enhances the THz-absorption by known as the incidence momentum and h is the thick- 3 a large factor of 10 . ness of the dielectric. Fig. 17 (B - C) shows the simu- lated transmittance for different complex refractive in- The effects of different substrates on biomaterial dices. It is evident from the ﬁgure, that the absorption detection were studied using THz-antenna (length = is dependent on the imaginary part of the complex re- 100 mm) with 10 10 slot-array each with a slot-width fractive index (k), and the resonance frequency is solely of 2 mm and periodicity of 200 mm (Fig. 19 (A - B)), was dependent on the real part of the complex refractive in- lithographed using electron-beam on two different sub- dex (n). However, the change in resonance frequency is strates, Si (undoped) and quartz, to sense yeast sam- quite small but stronger absorption can lead to an ap- ples. Using FFT (Fast Fourier Transform) algorithm, preciable change. Therefore, the variation in the trans- the normalized THz-transmission amplitudes for THz mission spectra will provide evident proofs for identiﬁ- slot-antenna, with (an average of 50 yeast molecules, cation of molecules and species contained in the bioma- or N = 50) and without yeast sample was analyzed av terial sample. (Fig. 19 (C)). From the ﬁgure, an evident 9 GHz red- To study the effects of various gap-widths on THz- shift was observed in the resonance frequency. To com- 16 A B C D E Figure 21. Detection of virus samples using THz nanoantenna sensing chip. (A) Schematic of a THz nano-slot antenna array sens- ing chip used to detect virus samples in liquid state. (B) Transmittance spectra measured using a multi-resonance nanoantenna, with and without H9N2. (C) Difference in transmitted intensity (DT) and resonance frequency shifts from each fundamental reso- nance peak (Df) for H9N2 virus sample. (D) THz transmission spectra through single-resonance nanoantenna, with and without virus samples (H1N1, H5N2 and H9N2). (E) Difference in transmitted intensity (DT) and resonance frequency shifts from the fundamental resonance peak (Df) of single-resonance nanoantenna for different virus samples. Reprinted with permission from Ref. . © 2017, Springer Nature. pare the effects of different substrate, the normalized ages of various brands in the THz range (0.5-2.5 THz) THz-transmission amplitudes were measured for Si- are shown in Fig. 20 (E). The low transmittances seen in and quartz-substrate THz antennas as shown in Fig. 19 the ﬁgure are due to the two artiﬁcially added sweet- (D - E). Comparing the measured amplitudes, quartz eners whose THz-absorptions is shown in Fig. 20 (F). showed stronger red-shifts compared to Si-substrate Hence, nanoantennas can provide highly sensitive de- (Fig. 19 (F)), showing 1/e dependence of the resonant tection even at low molecular concentration. e f f Recently, Lee et al. studied the THz-transmittances frequency-shift (Df/f ), where e = n is known as e f f e f f of various virus samples (H1N1, H5N2, AND H9N2) us- effective dielectric constant. Hence it can be concluded ing two THz nanoantennas: i) multi-resonance nanoan- that by using low dielectric constant substrate, the sen- tenna (with resonance frequencies of 0.63, 0.93 and 1.31 sitivity of a THz nanoantenna can be enhanced. THz), and ii) single-resonance nanoantenna (with reso- 88 nance frequency of 1.4 THz), which are shown in Fig. 21 THz nanoantennas are also selective in nature. This (B - E). Similar work has been reported by Park et al. in selective detection was shown using a nanoantenna of which they studied the transmission amplitude in both length 35 mm having a resonance frequency of 1.7 THz. presence and absence of virus samples (PRD1 and MS2 This antenna was specially designed to differentiate be- viruses) using a THz nanogap metamaterial (Fig. 22). tween fructose and D-glucose, and is also known as fructose antenna (Fig. 20 (B)). These molecules were dis- criminated from their measured transmittances using V. CONCLUSION AND FUTURE OUTLOOK the fructose antenna as shown in Fig. 20 (C - D). Using the same nanoantenna, the sugars and low-concentrated artiﬁcial sweeteners (acesulfame K and aspartame) con- THz waves are less penetrating and non-ionising elec- tained in beverages and dietary sodas, respectively, of tromagnetic waves. These waves are also known as popular brands were detected. The detection of these ’sub-millimeter waves’ and hence, is widely used in the artiﬁcial sweeteners is important as these are recently ﬁelds of astronomy and spectroscopy. Since the descrip- found to be addictive and toxic in nature. The different tion of electromagnetic interactions with metal is elu- transmittances of sugars contained in sweetened bever- cidated by using Maxwell's equations, researchers are 17 A C B D Figure 22. Detection of virus samples using THz nanogap metamaterial. (A) Schematic of THz nanogap metamaterial used to sense virus samples. (B) SEM image virus samples deposited on a THz nanogap metamaterial with 200 nm gap-width (w). THz transmission spectra of THz metamaterial (w = 2mm) containing (C) PRD1 and (D) MS2 virus samples. Noticeable red-shifts in resonance frequency are seen describing the sensitive detection of viruses using THz nanogap metamaterial. Reprinted with permission from Ref. . © 2017, Optical Society of America. very much fascinated to study the electric ﬁeld local- sitive and selective. Due to its enhanced detection in low ization and other plasmonic effects in various metal- molecular concentrations, in future nanoantennas can lic nanostructures. Numerous investigations reported be used to monitor blood sugar level. As an emerging that the electromagnetic interactions with metal ﬁlm in- future technology, nanoantennas can be used for early 243,244 crease the mobility of charges in the ﬁlm, tending it to detection of cancerous tumours which will lead to move towards the metal edges, leading to the enhance- the development of effective cancer treatment. Few re- ment of the electric ﬁeld inside different nanogap struc- cent works have shown that biosensing based on surface tures. To study the light-matter interactions, researchers plasmon resonance (SPR) has proved to be effective in 245–248 are fascinated towards THz waves as they can squeeze label-free tumor detection and can be used to de- through sub-nanometer metallic gaps (THz nanoanten- tect a single molecule of an early stage tumor. Currently, nas) and show non-linear optical responses. These research and development are in progress in making a nanoantennas are further integrated with different ma- nanoantenna based SPR biosensor which will be able to terials to study its novel plasmonic properties. Recently, detect cancer cells at an early stage. successful fabrication of graphene-integrated plasmonic 70,167,239–242 system has boosted further investigations on their plasmonic properties because graphene is atomi- VI. ACKNOWLEDGEMENT cally thin and is electrically tunable. As discussed earlier in this review, THz nanoanten- The authors thanks Birla Institute of Technology, nas have played a crucial role as biosensors for ultra- Mesra, Ranchi for providing research facilities and sensitive detection of various biomolecules and bioma- MHRD, government of India for support through terials, which is our prime focus in this review. We have TEQIP - III. The authors also thanks Pawan Kumar come across various literatures in which we have seen Dubey, Akriti Raj, Sameer Kumar Tiwari, Kamana the development of nanoantennas, becoming more sen- Mishra, Priyanshi Srivastava, Dhruv Sood and Devo- 18 tosh Ganguly for useful discussions and proof reading VIII. KEYWORDS of the manuscript. VII. CONFLICT OF INTEREST The authors declares no conﬂict of interest. firstname.lastname@example.org revb.95.075424. Y.-M. Bahk, D.-S. Kim, and H.-R. Park, “Large-area metal P. Bharadwaj, B. Deutsch, and L. Novotny, “Optical anten- gaps and their optical applications,” Advanced Optical nas,” Advances in Optics and Photonics 1, 438 (2009). Materials 7, 1800426 (2018). P. Muhlschlegel, “Resonant optical antennas,” Science 308, N. Bonod, E. Popov, D. Gerar ´ d, J. Wenger, and 1607 (2005). H. Rigneault, “Field enhancement in a circular aperture J. Alda, J. M. Rico-Garc´ ıa, J. M. Lopez-Alonso, ´ and G. Bore- surrounded by a single channel groove,” Optics Express man, “Optical antennas for nano-photonic applications,” 16, 2276 (2008). Nanotechnology 16, S230 (2005). 4 H. Cao and A. Nahata, “Resonantly enhanced transmis- B. M. Ross and L. P. Lee, “Comparison of near- and far- sion of terahertz radiation through a periodic array of sub- ﬁeld measures for plasmon resonance of metallic nanopar- wavelength apertures,” Optics Express 12, 1004 (2004). ticles,” Optics Letters 34, 896 (2009). 5 H.-T. Chen, W. J. Padilla, J. M. O. Zide, A. C. Gossard, A. J. A. Alu ` and N. Engheta, “Input impedance, nanocir- Taylor, and R. D. Averitt, “Active terahertz metamaterial cuit loading, and radiation tuning of optical nanoanten- devices,” Nature 444, 597 (2006). nas,” Physical Review Letters 101 (2008), 10.1103/phys- X. Chen, H.-R. Park, N. C. Lindquist, J. Shaver, M. Pelton, revlett.101.043901. 6 and S.-H. Oh, “Squeezing millimeter waves through a sin- L. Novotny, “Effective wavelength scaling for optical an- gle, nanometer-wide, centimeter-long slit,” Scientiﬁc Re- tennas,” Physical Review Letters 98 (2007), 10.1103/phys- ports 4 (2014), 10.1038/srep06722. revlett.98.266802. J. Jeong, J. Rhie, W. Jeon, C. S. Hwang, and D.-S. Kim, K. Kneipp, Y. Wang, H. Kneipp, L. T. Perelman, I. Itzkan, “High-throughput fabrication of inﬁnitely long 10 nm slit R. R. Dasari, and M. S. Feld, “Single molecule detection arrays for terahertz applications,” Journal of Infrared, Mil- using surface-enhanced raman scattering (SERS),” Physi- limeter, and Terahertz Waves 36, 262 (2014). cal Review Letters 78, 1667 (1997). J.-H. Kang, D.-S. Kim, and M. Seo, “Terahertz wave in- S. Nie, “Probing single molecules and single nanoparticles teraction with metallic nanostructures,” Nanophotonics 7, by surface-enhanced raman scattering,” Science 275, 1102 763 (2018). (1997). J.-Y. Kim, B. J. Kang, J. Park, Y.-M. Bahk, W. T. Kim, H. Xu, E. J. Bjerneld, M. Kll, and L. Brjesson, “Spec- J. Rhie, H. Jeon, F. Rotermund, and D.-S. Kim, “Tera- troscopy of single hemoglobin molecules by surface en- hertz quantum plasmonics of nanoslot antennas in nonlin- hanced raman scattering,” Physical Review Letters 83, ear regime,” Nano Letters 15, 6683 (2015). 4357 (1999). N. Kim, S. In, D. Lee, J. Rhie, J. Jeong, D.-S. Kim, and M. L. Juan, M. Righini, and R. Quidant, “Plasmon nano- N. Park, “Colossal terahertz ﬁeld enhancement using split- optical tweezers,” Nature Photonics 5, 349 (2011). ring resonators with a sub-10 nm gap,” ACS Photonics 5, O. M. Marago, ` P. H. Jones, P. G. Gucciardi, G. Volpe, 278 (2017). and A. C. Ferrari, “Optical trapping and manipulation of J. W. Lee, M. A. Seo, D. H. Kang, K. S. Khim, S. C. Jeoung, nanostructures,” Nature Nanotechnology 8, 807 (2013). and D. S. Kim, “Terahertz electromagnetic wave transmis- D. R. Ward, F. Hser, F. Pauly, J. C. Cuevas, and D. Natel- sion through random arrays of single rectangular holes son, “Optical rectiﬁcation and ﬁeld enhancement in a plas- and slits in thin metallic sheets,” Physical Review Letters monic nanogap,” Nature Nanotechnology 5, 732 (2010). 99 (2007), 10.1103/physrevlett.99.137401. M. Kauranen and A. V. Zayats, “Nonlinear plasmonics,” H. O. Moser, B. D. F. Casse, O. Wilhelmi, and B. T. Saw, Nature Photonics 6, 737 (2012). 14 “Terahertz response of a microfabricated rod–split-ring- Y.-M. Bahk, G. Ramakrishnan, J. Choi, H. Song, G. Choi, resonator electromagnetic metamaterial,” Physical Review Y. H. Kim, K. J. Ahn, D.-S. Kim, and P. C. M. Planken, Letters 94 (2005), 10.1103/physrevlett.94.063901. “Plasmon enhanced terahertz emission from single layer A. Novitsky, A. M. Ivinskaya, M. Zalkovskij, R. Malure- graphene,” ACS Nano 8, 9089 (2014). 15 anu, P. U. Jepsen, and A. V. Lavrinenko, “Non-resonant J. B. Lassiter, X. Chen, X. Liu, C. Cirac` ı, T. B. Hoang, terahertz ﬁeld enhancement in periodically arranged S. Larouche, S.-H. Oh, M. H. Mikkelsen, and D. R. nanoslits,” Journal of Applied Physics 112, 074318 (2012). Smith, “Third-harmonic generation enhancement by ﬁlm- A. Novitsky, M. Zalkovskij, R. Malureanu, and A. Lavri- coupled plasmonic stripe resonators,” ACS Photonics 1, nenko, “Microscopic model of the THz ﬁeld enhancement 1212 (2014). 16 in a metal nanoslit,” Optics Communications 284, 5495 Y.-M. Bahk, S. Han, J. Rhie, J. Park, H. Jeon, N. Park, and (2011). D.-S. Kim, “Ultimate terahertz ﬁeld enhancement of sin- gle nanoslits,” Physical Review B 95 (2017), 10.1103/phys- 19 H. R. Park, Y. M. Park, H. S. Kim, J. S. Kyoung, M. A. Seo, emission by nano-electrodes in a photoconductive pho- D. J. Park, Y. H. Ahn, K. J. Ahn, and D. S. Kim, “Ter- tomixer,” Nature Photonics 6, 121 (2012). ahertz nanoresonators: Giant ﬁeld enhancement and ul- L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, trabroadband performance,” Applied Physics Letters 96, H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, and F. Wang, 121106 (2010). “Graphene plasmonics for tunable terahertz metamateri- D. Qu, D. Grischkowsky, and W. Zhang, “Terahertz trans- als,” Nature Nanotechnology 6, 630 (2011). mission properties of thin, subwavelength metallic hole ar- M. Shur, “Plasma wave terahertz electronics,” Electronics rays,” Optics Letters 29, 896 (2004). Letters 46, S18 (2010). 32 46 J. G. Rivas, C. Schotsch, P. H. Bolivar, and H. Kurz, “En- P. Upadhyaya, S. Pramanik, and S. Bandyopadhyay, “Op- hanced transmission of THz radiation through subwave- tical transitions in a quantum wire with spin-orbit inter- length holes,” Physical Review B 68 (2003), 10.1103/phys- action and its applications in terahertz electronics: Be- revb.68.201306. yond zeroth-order theory,” Physical Review B 77 (2008), M. A. Seo, A. J. L. Adam, J. H. Kang, J. W. Lee, K. J. Ahn, 10.1103/physrevb.77.155439. Q. H. Park, P. C. M. Planken, and D. S. Kim, “Near ﬁeld M. Bahramipanah, M. S. Abrishamian, S. A. Mirtaheri, imaging of terahertz focusing onto rectangular apertures,” and J.-M. Liu, “Ultracompact plasmonic loop–stub notch Optics Express 16, 20484 (2008). ﬁlter and sensor,” Sensors and Actuators B: Chemical 194, M. A. Seo, H. R. Park, S. M. Koo, D. J. Park, J. H. Kang, 311 (2014). O. K. Suwal, S. S. Choi, P. C. M. Planken, G. S. Park, N. K. D. M. Mittleman, “Frontiers in terahertz sources and plas- Park, Q. H. Park, and D. S. Kim, “Terahertz ﬁeld enhance- monics,” Nature Photonics 7, 666 (2013). ment by a metallic nano slit operating beyond the skin- C. R. Williams, S. R. Andrews, S. A. Maier, A. I. Fernandez- ´ depth limit,” Nature Photonics 3, 152 (2009). Dom´ ınguez, L. Mart´ ın-Moreno, and F. J. Garc´ ıa-Vidal, “Highly conﬁned guiding of terahertz surface plasmon po- M. Shalaby, H. Merbold, M. Peccianti, L. Razzari, laritons on structured metal surfaces,” Nature Photonics 2, G. Sharma, T. Ozaki, R. Morandotti, T. Feurer, A. Weber, L. Heyderman, B. Patterson, and H. Sigg, “Concurrent 175 (2008). A. J. Fitzgerald, E. Berry, N. N. Zinovev, G. C. Walker, M. A. ﬁeld enhancement and high transmission of THz radia- Smith, and J. M. Chamberlain, “An introduction to medi- tion in nanoslit arrays,” Applied Physics Letters 99, 041110 cal imaging with coherent terahertz frequency radiation,” (2011). Physics in Medicine and Biology 47, R67 (2002). A. Toma, S. Tuccio, M. Prato, F. D. Donato, A. Perucchi, P. D. Pietro, S. Marras, C. Liberale, R. P. Zaccaria, F. D. An- E. P. J. Parrott, Y. Sun, and E. Pickwell-MacPherson, “Ter- gelis, L. Manna, S. Lupi, E. D. Fabrizio, and L. Razzari, ahertz spectroscopy: Its future role in medical diagnoses,” “Squeezing terahertz light into nanovolumes: Nanoan- Journal of Molecular Structure 1006, 66 (2011). tenna enhanced terahertz spectroscopy (NETS) of semi- J. W. Handley, A. J. Fitzgerald, E. Berry, and R. D. Boyle, conductor quantum dots,” Nano Letters 15, 386 (2014). “Wavelet compression in medical terahertz pulsed imag- ing,” Physics in Medicine and Biology 47, 3885 (2002). L.-N. Tripathi, Y.-M. Bahk, G. Choi, S. Han, N. Park, and D.-S. Kim, “Terahertz transmission through rings of quan- W. Tu, S. Zhong, Y. Shen, and A. Incecik, “Nondestruc- tum dots-nanogap,” Applied Physics Express 9, 032001 tive testing of marine protective coatings using terahertz (2016). waves with stationary wavelet transform,” Ocean Engi- C. A. Werley, K. Fan, A. C. Strikwerda, S. M. Teo, X. Zhang, neering 111, 582 (2016). R. D. Averitt, and K. A. Nelson, “Time-resolved imaging J. Li and Y. Pi, “Target detection for terahertz radar net- of near-ﬁelds in THz antennas and direct quantitative mea- works based on micro-doppler signatures,” International surement of ﬁeld enhancements,” Optics Express 20, 8551 Journal of Sensor Networks 17, 115 (2015). (2012). E. N. Semashkin and T. V. Artyushkina, “Operating range W. Park, Y. Lee, T. Kang, J. Jeong, and D.-S. Kim, and all-weather capability of terahertz (01 THz) and gi- “Terahertz-driven polymerization of resists in nanoan- gahertz (3–333 GHz) radars on horizontal and oblique tennas,” Scientiﬁc Reports 8 (2018), 10.1038/s41598-018- tracks,” Journal of Optical Technology 82, 430 (2015). 26214-w. D. Zimdars, J. S. White, G. Stuk, A. Chernovsky, G. Fichter, H.-R. Park, Y.-M. Bahk, K. J. Ahn, Q.-H. Park, D.-S. Kim, and S. Williamson, “Security and non destructive evalua- L. Mart´ ın-Moreno, F. J. Garc´ ıa-Vidal, and J. Bravo-Abad, tion application of high speed time domain terahertz imag- “Controlling terahertz radiation with nanoscale metal bar- ing,” in 2006 Conference on Lasers and Electro-Optics and riers embedded in nano slot antennas,” ACS Nano 5, 8340 2006 Quantum Electronics and Laser Science Conference (IEEE, (2011). 2006). N. C. Lindquist, P. Nagpal, K. M. McPeak, D. J. Norris, 57 S.-P. Han, N. Kim, W.-H. Lee, E. S. Lee, H. Ko, I.-M. Lee, and S.-H. Oh, “Engineering metallic nanostructures for K. Moon, D. H. Lee, and K. H. Park, “Real-time imaging of plasmonics and nanophotonics,” Reports on Progress in moving living objects using a compact terahertz scanner,” Physics 75, 036501 (2012). Applied Physics Express 9, 022501 (2016). M. C. Hoffmann, B. S. Monozon, D. Livshits, E. U. Rafailov, 58 K. Fukunaga and I. Hosako, “Innovative non-invasive and D. Turchinovich, “Terahertz electro-absorption effect analysis techniques for cultural heritage using terahertz enabling femtosecond all-optical switching in semicon- technology,” Comptes Rendus Physique 11, 519 (2010). ductor quantum dots,” Applied Physics Letters 97, 231108 59 J.-M. Manceau, A. Nevin, C. Fotakis, and S. Tzortzakis, (2010). “Terahertz time domain spectroscopy for the analysis of H. Tanoto, J. H. Teng, Q. Y. Wu, M. Sun, Z. N. Chen, S. A. cultural heritage related materials,” Applied Physics B 90, Maier, B. Wang, C. C. Chum, G. Y. Si, A. J. Danner, and 365 (2008). S. J. Chua, “Greatly enhanced continuous-wave terahertz 20 60 76 J. W. Lee, M. A. Seo, D. J. Park, S. C. Jeoung, Q. H. Park, K. Xu, “Integrated silicon directly modulated light source C. Lienau, and D. S. Kim, “Terahertz transparency at using p-well in standard CMOS technology,” IEEE Sensors fabry-perot resonances of periodic slit arrays in a metal Journal 16, 6184 (2016). plate: experiment and theory,” Optics Express 14, 12637 Y.-S. Lee, Principles of Terahertz Science and Technology (2006). (SPRINGER NATURE, 2008). 61 78 I. H. Libon, S. Baumgrtner, M. Hempel, N. E. Hecker, J. B. Baxter and G. W. Guglietta, “Terahertz spectroscopy,” J. Feldmann, M. Koch, and P. Dawson, “An optically con- Analytical Chemistry 83, 4342 (2011). trollable terahertz ﬁlter,” Applied Physics Letters 76, 2821 J. X. Xi-Cheng Zhang, Introduction to THz Wave Photonics (2000). (Springer-Verlag GmbH, 2009). R. Mendis, A. Nag, F. Chen, and D. M. Mittleman, “A S. Dexheimer, ed., Terahertz Spectroscopy (CRC Press, 2007). tunable universal terahertz ﬁlter using artiﬁcial dielectrics 81 M. Nagel, M. Frst, and H. Kurz, “THz biosensing devices: based on parallel-plate waveguides,” Applied Physics Let- fundamentals and technology,” Journal of Physics: Con- ters 97, 131106 (2010). densed Matter 18, S601 (2006). M. A. Kaliteevski, S. Brand, J. Garvie-Cook, R. A. Abram, 82 S. Kar, V. L. Nguyen, D. R. Mohapatra, Y. H. Lee, and and J. M. Chamberlain, “Terahertz ﬁlter based on refractive A. K. Sood, “Ultrafast spectral photoresponse of bilayer properties of metallic photonic crystal,” Optics Express 16, graphene: Optical pump–terahertz probe spectroscopy,” 7330 (2008). ACS Nano 12, 1785 (2018). F. Lan, Z. Yang, L. Qi, X. Gao, and Z. Shi, “Terahertz 83 B. M. Fischer, M. Walther, and P. U. Jepsen, “Far-infrared dual-resonance bandpass ﬁlter using bilayer reformative vibrational modes of DNA components studied by tera- complementary metamaterial structures,” Optics Letters hertz time-domain spectroscopy,” Physics in Medicine and 39, 1709 (2014). Biology 47, 3807 (2002). B. Heshmat, H. Pahlevaninezhad, P. Yuanjie, M. Masnadi- A. Markelz, S. Whitmire, J. Hillebrecht, and R. Birge, “THz Shirazi, R. B. Lewis, T. Tiedje, R. Gordon, and T. E. Dar- time domain spectroscopy of biomolecular conformational cie, “Nanoplasmonic terahertz photoconductive switch on modes,” Physics in Medicine and Biology 47, 3797 (2002). GaAs,” Nano Letters 12, 6255 (2012). A. Arora, T. Q. Luong, M. Krger, Y. J. Kim, C.-H. Nam, M. Rahm, J.-S. Li, and W. J. Padilla, “THz wave modula- A. Manz, and M. Havenith, “Terahertz-time domain spec- tors: A brief review on different modulation techniques,” troscopy for the detection of PCR ampliﬁed DNA in aque- Journal of Infrared, Millimeter, and Terahertz Waves 34, 1 ous solution,” The Analyst 137, 575 (2012). (2012). X. Yang, D. Wei, S. Yan, Y. Liu, S. Yu, M. Zhang, Z. Yang, S. B. Choi, J. S. Kyoung, H. S. Kim, H. R. Park, D. J. X. Zhu, Q. Huang, H.-L. Cui, and W. Fu, “Rapid and Park, B.-J. Kim, Y. H. Ahn, F. Rotermund, H.-T. Kim, K. J. label-free detection and assessment of bacteria by tera- Ahn, and D. S. Kim, “Nanopattern enabled terahertz all- hertz time-domain spectroscopy,” Journal of Biophotonics optical switching on vanadium dioxide thin ﬁlm,” Applied 9, 1050 (2016). Physics Letters 98, 071105 (2011). X. Yang, K. Yang, Y. Luo, and W. Fu, “Terahertz spec- M. Unlu and M. Jarrahi, “Miniature multi-contact MEMS troscopy for bacterial detection: opportunities and chal- switch for broadband terahertz modulation,” Optics Ex- lenges,” Applied Microbiology and Biotechnology 100, press 22, 32245 (2014). 5289 (2016). J. Li, J. He, and Z. Hong, “Terahertz wave switch based on D.-K. Lee, J.-H. Kang, J.-S. Lee, H.-S. Kim, C. Kim, J. H. silicon photonic crystals,” Applied Optics 46, 5034 (2007). Kim, T. Lee, J.-H. Son, Q.-H. Park, and M. Seo, “Highly W. Gao, J. Shu, K. Reichel, D. V. Nickel, X. He, G. Shi, sensitive and selective sugar detection by terahertz nano- R. Vajtai, P. M. Ajayan, J. Kono, D. M. Mittleman, and antennas,” Scientiﬁc Reports 5 (2015), 10.1038/srep15459. Q. Xu, “High-contrast terahertz wave modulation by H.-R. Park, K. J. Ahn, S. Han, Y.-M. Bahk, N. Park, and gated graphene enhanced by extraordinary transmission D.-S. Kim, “Colossal absorption of molecules inside single through ring apertures,” Nano Letters 14, 1242 (2014). terahertz nanoantennas,” Nano Letters 13, 1782 (2013). G. Liu, M. He, Z. Tian, J. Li, and J. Liu, “Terahertz sur- 90 O. Limaj, D. Etezadi, N. J. Wittenberg, D. Rodrigo, D. Yoo, face plasmon sensor for distinguishing gasolines,” Ap- S.-H. Oh, and H. Altug, “Infrared plasmonic biosensor for plied Optics 52, 5695 (2013). real-time and label-free monitoring of lipid membranes,” V. Astley, K. Reichel, R. Mendis, and D. M. Mittle- Nano Letters 16, 1502 (2016). man, “Terahertz microﬂuidic sensing using a parallel-plate D. Rodrigo, O. Limaj, D. Janner, D. Etezadi, F. J. G. waveguide sensor,” Journal of Visualized Experiments de Abajo, V. Pruneri, and H. Altug, “Mid-infrared plas- (2012), 10.3791/4304. monic biosensing with graphene,” Science 349, 165 (2015). F. Alves, D. Grbovic, B. Kearney, and G. Karunasiri, “Mi- 92 D.-K. Lee, J.-H. Kang, J. Kwon, J.-S. Lee, S. Lee, D. H. Woo, croelectromechanical systems bimaterial terahertz sensor J. H. Kim, C.-S. Song, Q.-H. Park, and M. Seo, “Nano with integrated metamaterial absorber,” Optics Letters 37, metamaterials for ultrasensitive terahertz biosensing,” Sci- 1886 (2012). entiﬁc Reports 7 (2017), 10.1038/s41598-017-08508-7. A. Hassani and M. Skorobogatiy, “Surface plasmon S. J. Park, J. T. Hong, S. J. Choi, H. S. Kim, W. K. Park, S. T. resonance-like integrated sensor at terahertz frequencies Han, J. Y. Park, S. Lee, D. S. Kim, and Y. H. Ahn, “De- for gaseous analytes,” Optics Express 16, 20206 (2008). tection of microorganisms using terahertz metamaterials,” F. Miyamaru, S. Hayashi, C. Otani, K. Kawase, Y. Ogawa, Scientiﬁc Reports 4 (2014), 10.1038/srep04988. H. Yoshida, and E. Kato, “Terahertz surface-wave reso- S. J. Park, S. H. Cha, G. A. Shin, and Y. H. Ahn, “Sensing nant sensor with a metal hole array,” Optics Letters 31, viruses using terahertz nano-gap metamaterials,” Biomed- 1118 (2006). ical Optics Express 8, 3551 (2017). 21 X. Chen, H.-R. Park, M. Pelton, X. Piao, N. C. Lindquist, measurements,” Applied Physics Letters 92, 113102 (2008). H. Im, Y. J. Kim, J. S. Ahn, K. J. Ahn, N. Park, D.-S. C. Vieu, F. Carcenac, A. Pepin, ´ Y. Chen, M. Mejias, Kim, and S.-H. Oh, “Atomic layer lithography of wafer- A. Lebib, L. Manin-Ferlazzo, L. Couraud, and H. Launois, scale nanogap arrays for extreme conﬁnement of elec- “Electron beam lithography: resolution limits and applica- tromagnetic waves,” Nature Communications 4 (2013), tions,” Applied Surface Science 164, 111 (2000). 10.1038/ncomms3361. A. Broers, A. Hoole, and J. Ryan, “Electron beam lithog- J. Lee, M. Seo, D. Park, D. Kim, S. Jeoung, C. Lienau, Q.-H. raphy—resolution limits,” Microelectronic Engineering 32, Park, and P. Planken, “Shape resonance omni-directional 131 (1996). terahertz ﬁlters with near-unity transmittance,” Optics Ex- Y. Chen, “Nanofabrication by electron beam lithography press 14, 1253 (2006). and its applications: A review,” Microelectronic Engineer- J. W. Lee, M. A. Seo, D. S. Kim, S. C. Jeoung, C. Lienau, J. H. ing 135, 57 (2015). Kang, and Q.-H. Park, “Fabry–perot effects in THz time- A. Tseng, K. Chen, C. Chen, and K. Ma, “Electron domain spectroscopy of plasmonic band-gap structures,” beam lithography in nanoscale fabrication: recent de- Applied Physics Letters 88, 071114 (2006). velopment,” IEEE Transactions on Electronics Packaging D. Ward, E. Statz, and K. Nelson, “Fabrication of polari- Manufacturing 26, 141 (2003). tonic structures in LiNbO3 and LiTaO3 using femtosecond J. C. Hulteen and R. P. V. Duyne, “Nanosphere lithography: laser machining,” Applied Physics A 86, 49 (2006). A materials general fabrication process for periodic parti- C. Enkrich, F. Per ´ ez-Willard, D. Gerthsen, J. F. Zhou, cle array surfaces,” Journal of Vacuum Science & Technol- T. Koschny, C. M. Soukoulis, M. Wegener, and S. Linden, ogy A: Vacuum, Surfaces, and Films 13, 1553 (1995). “Focused-ion-beam nanofabrication of near-infrared mag- P. Vettiger, “Nanometer sidewall lithography by resist sily- netic metamaterials,” Advanced Materials 17, 2547 (2005). lation,” Journal of Vacuum Science & Technology B: Micro- N. Ocelic and R. Hillenbrand, “Subwavelength-scale tai- electronics and Nanometer Structures 7, 1756 (1989). loring of surface phonon polaritons by focused ion-beam S. Y. Chou, P. R. Krauss, and P. J. Renstrom, “Imprint implantation,” Nature Materials 3, 606 (2004). lithography with 25-nanometer resolution,” Science 272, H. F. Ghaemi, T. Thio, D. E. Grupp, T. W. Ebbesen, and 85 (1996). H. J. Lezec, “Surface plasmons enhance optical transmis- S. Y. Chou, “Nanoimprint lithography,” Journal of Vacuum sion through subwavelength holes,” Physical Review B 58, Science & Technology B: Microelectronics and Nanometer 6779 (1998). Structures 14, 4129 (1996). T. W. Ebbesen, H. J. Lezec, H. F. Ghaemi, T. Thio, and M. D. Austin, H. Ge, W. Wu, M. Li, Z. Yu, D. Wasserman, P. A. Wolff, “Extraordinary optical transmission through S. A. Lyon, and S. Y. Chou, “Fabrication of 5nm linewidth sub-wavelength hole arrays,” Nature 391, 667 (1998). and 14nm pitch features by nanoimprint lithography,” Ap- H.-R. Park, Y.-M. Bahk, J. H. Choe, S. Han, S. S. Choi, K. J. plied Physics Letters 84, 5299 (2004). Ahn, N. Park, Q.-H. Park, and D.-S. Kim, “Terahertz pinch L. Guo, “Nanoimprint lithography: Methods and material harmonics enabled by single nano rods,” Optics Express requirements,” Advanced Materials 19, 495 (2007). 19, 24775 (2011). L. Guo, P. R. Krauss, and S. Y. Chou, “Nanoscale silicon D. C. Bell, M. C. Lemme, L. A. Stern, J. R. Williams, ﬁeld effect transistors fabricated using imprint lithogra- and C. M. Marcus, “Precision cutting and patterning of phy,” Applied Physics Letters 71, 1881 (1997). graphene with helium ions,” Nanotechnology 20, 455301 S. C. Ambhire, S. Palkhivala, A. Agrawal, A. Gupta, (2009). G. Rana, R. Mehta, D. Ghindani, A. Bhattacharya, V. G. H.-T. Chen, J. F. O'Hara, A. K. Azad, A. J. Taylor, R. D. Achanta, and S. S. Prabhu, ““pattern and peel” method Averitt, D. B. Shrekenhamer, and W. J. Padilla, “Experi- for fabricating mechanically tunable terahertz metasurface mental demonstration of frequency-agile terahertz meta- on an elastomeric substrate,” Optical Materials Express 8, materials,” Nature Photonics 2, 295 (2008). 3382 (2018). R. Singh, E. Plum, C. Menzel, C. Rockstuhl, A. K. Azad, C. Feuillet-Palma, Y. Todorov, A. Vasanelli, and C. Sirtori, R. A. Cheville, F. Lederer, W. Zhang, and N. I. Zheludev, “Strong near ﬁeld enhancement in THz nano-antenna ar- “Terahertz metamaterial with asymmetric transmission,” rays,” Scientiﬁc Reports 3 (2013), 10.1038/srep01361. Physical Review B 80 (2009), 10.1103/physrevb.80.153104. O. K. Suwal, J. Rhie, N. Kim, and D.-S. Kim, “Nonresonant M. Choi, S. H. Lee, Y. Kim, S. B. Kang, J. Shin, M. H. Kwak, 104 terahertz ﬁeld enhancement with 5-nm slits,” Scientiﬁc K.-Y. Kang, Y.-H. Lee, N. Park, and B. Min, “A terahertz Reports 7 (2017), 10.1038/srep45638. metamaterial with unnaturally high refractive index,” Na- W. Park, J. Rhie, N. Y. Kim, S. Hong, and D.-S. Kim, “Sub- ture 470, 369 (2011). 10 nm feature chromium photomasks for contact lithogra- C. G. Willson, R. R. Dammel, and A. Reiser, “Photoresist phy patterning of square metal ring arrays,” Scientiﬁc Re- materials: a historical perspective,” in Metrology, Inspec- ports 6 (2016), 10.1038/srep23823. tion, and Process Control for Microlithography XI (SPIE, 1997). H. Im, K. C. Bantz, N. C. Lindquist, C. L. Haynes, and S.- C. Lin, C. Chen, G. J. Schneider, P. Yao, S. Shi, A. Sharkawy, H. Oh, “Vertically oriented sub-10-nm plasmonic nanogap and D. W. Prather, “Wavelength scale terahertz two- arrays,” Nano Letters 10, 2231 (2010). dimensional photonic crystal waveguides,” Optics Ex- H. T. Miyazaki and Y. Kurokawa, “Squeezing visible light press 12, 5723 (2004). waves into a 3-nm-thick and 55-nm-long plasmon cav- H.-R. Park, X. Chen, N.-C. Nguyen, J. Peraire, and S.-H. ity,” Physical Review Letters 96 (2006), 10.1103/phys- Oh, “Nanogap-enhanced terahertz sensing of 1 nm thick revlett.96.097401. (l/106) dielectric ﬁlms,” ACS Photonics 2, 417 (2015). W. Zhu, M. G. Banaee, D. Wang, Y. Chu, and K. B. Crozier, A. Fursina, S. Lee, R. G. S. Soﬁn, I. V. Shvets, and D. Natel- “Lithographically fabricated optical antennas with gaps son, “Nanogaps with very large aspect ratios for electrical 22 well below 10 nm,” Small 7, 1761 (2011). hertz dipole nanoantenna arrays: Resonance characteris- W. Fan, S. Zhang, K. J. Malloy, and S. R. J. Brueck, tics,” Plasmonics 8, 133 (2012). “Enhanced mid-infrared transmission through nanoscale Y.-G. Jeong, M. J. Paul, S.-H. Kim, K.-J. Yee, D.-S. Kim, and metallic coaxial-aperture arrays,” Optics Express 13, 4406 Y.-S. Lee, “Large enhancement of nonlinear terahertz ab- (2005). sorption in intrinsic GaAs by plasmonic nano antennas,” D. Ji, A. Cheney, N. Zhang, H. Song, J. Gao, X. Zeng, Applied Physics Letters 103, 171109 (2013). H. Hu, S. Jiang, Z. Yu, and Q. Gan, “Efﬁcient mid- J. Shu, C. Qiu, V. Astley, D. Nickel, D. M. Mittleman, and infrared light conﬁnement within sub-5-nm gaps for ex- Q. Xu, “High-contrast terahertz modulator based on ex- treme ﬁeld enhancement,” Advanced Optical Materials 5, traordinary transmission through a ring aperture,” Optics 1700223 (2017). Express 19, 26666 (2011). T. Ding, D. Sigle, L. Zhang, J. Mertens, B. de Nijs, and H. Merbold, A. Bitzer, and T. Feurer, “Second harmonic J. Baumberg, “Controllable tuning plasmonic coupling generation based on strong ﬁeld enhancement in nanos- with nanoscale oxidation,” ACS Nano 9, 6110 (2015). tructured THz materials,” Optics Express 19, 7262 (2011). L. N. Tripathi, T. Kang, Y.-M. Bahk, S. Han, G. Choi, J. Rhie, 150 L. Dai and C. Jiang, “Anomalous near-perfect extraordi- J. Jeong, and D.-S. Kim, “Quantum dots-nanogap metama- nary optical absorption on subwavelength thin metal ﬁlm terials fabrication by self-assembly lithography and photo- grating.,” Optics Express 17, 20502 (2009). luminescence studies,” Optics Express 23, 14937 (2015). B. gt, R. Vogelgesang, W. Sigle, N. Talebi, C. T. Koch, and J. H. Kang, J.-H. Choe, D. S. Kim, and Q.-H. Park, “Sub- P. A. van Aken, “Hybridized metal slit eigenmodes as an il- strate effect on aperture resonances in a thin metal ﬁlm,” lustration of babinet’s principle,” ACS Nano 5, 6701 (2011). Optics Express 17, 15652 (2009). J. Yang, Q. Cao, and C. Zhou, “Theory for terahertz plas- J.-H. Choe, J.-H. Kang, D.-S. Kim, and Q.-H. Park, “Slot mons of metallic nanowires with sub-skin-depth diame- antenna as a bound charge oscillator,” Optics Express 20, ters,” Optics Express 18, 18550 (2010). 6521 (2012). D. Hu, X. Wang, S. Feng, J. Ye, W. Sun, Q. Kan, P. J. Klar, ´ ´ F. J. Garcıa-Vidal, L. Martın-Moreno, E. Moreno, L. K. S. and Y. Zhang, “Ultrathin terahertz planar elements,” Ad- Kumar, and R. Gordon, “Transmission of light through a vanced Optical Materials 1, 186 (2013). single rectangular hole in a real metal,” Physical Review B K. Iwaszczuk, A. Andryieuski, A. Lavrinenko, X.-C. 74 (2006), 10.1103/physrevb.74.153411. Zhang, and P. U. Jepsen, “Terahertz ﬁeld enhancement to S.-G. Park, Y. Choi, Y.-J. Oh, and K.-H. Jeong, “Terahertz the MV/cm regime in a tapered parallel plate waveguide,” photoconductive antenna with metal nanoislands,” Optics Optics Express 20, 8344 (2012). Express 20, 25530 (2012). D. S. Bulgarevich, M. Watanabe, and M. Shiwa, “Single Y. Yang, R. Singh, and W. Zhang, “Anomalous terahertz sub-wavelength aperture with greatly enhanced transmis- transmission in bow-tie plasmonic antenna apertures,” sion,” New Journal of Physics 14, 053001 (2012). Optics Letters 36, 2901 (2011). M. N. Gadalla, M. Abdel-Rahman, and A. Shamim, “De- I. S. Maksymov, A. E. Miroshnichenko, and Y. S. Kivshar, sign, optimization and fabrication of a 28.3 THz nano- “Actively tunable bistable optical yagi-uda nanoantenna,” rectenna for infrared detection and rectiﬁcation,” Scientiﬁc Optics Express 20, 8929 (2012). Reports 4 (2014), 10.1038/srep04270. T. H. Taminiau, F. D. Stefani, and N. F. van Hulst, “En- F. Fan, S.-T. Xu, X.-H. Wang, and S.-J. Chang, “Terahertz hanced directional excitation and emission of single emit- polarization converter and one-way transmission based on ters by a nano-optical yagi-uda antenna,” Optics Express double-layer magneto-plasmonics of magnetized InSb,” 16, 10858 (2008). Optics Express 24, 26431 (2016). G. Pelosi and S. Selleri, “Babinet's principle in electromag- T. Low and P. Avouris, “Graphene plasmonics for terahertz netics: Why does a slot radiate like a dipole? [historical to mid-infrared applications,” ACS Nano 8, 1086 (2014). corner],” IEEE Antennas and Propagation Magazine 59, F. Fan, W.-H. Gu, S. Chen, X.-H. Wang, and S.-J. Chang, 144 (2017). 142 “State conversion based on terahertz plasmonics with H. Booker, “Slot aerials and their relation to complemen- vanadium dioxide coating controlled by optical pump- tary wire aerials (babinet's principle),” Journal of the Insti- ing,” Optics Letters 38, 1582 (2013). tution of Electrical Engineers - Part IIIA: Radiolocation 93, A. H. Panaretos and D. H. Werner, “Spoof plasmon radi- 620 (1946). ation using sinusoidally modulated corrugated reactance Y. Mushiake, “A report on japanese development of an- surfaces,” Optics Express 24, 2443 (2016). tennas: from the yagi-uda antenna to self-complementary N. Yu, Q. J. Wang, M. A. Kats, J. A. Fan, S. P. Khanna, L. Li, antennas,” IEEE Antennas and Propagation Magazine 46, A. G. Davies, E. H. Linﬁeld, and F. Capasso, “Designer 47 (2004). spoof surface plasmon structures collimate terahertz laser M. Born and E. Wolf, Principles of Optics: Electromagnetic beams,” Nature Materials 9, 730 (2010). Theory of Propagation Interference and Diffraction of Light A. Ishikawa, S. Zhang, D. A. Genov, G. Bartal, and (Pergamon Pr, 1981). X. Zhang, “Deep subwavelength terahertz waveguides us- S. Koo, M. S. Kumar, J. Shin, D. Kim, and N. Park, ing gap magnetic plasmon,” Physical Review Letters 102 “Extraordinary magnetic ﬁeld enhancement with metallic (2009), 10.1103/physrevlett.102.043904. nanowire: Role of surface impedance in babinet’s principle R.-L. Chern, “Magnetic and surface plasmon resonances for sub-skin-depth regime,” Physical Review Letters 103 for periodic lattices of plasmonic split-ring resonators,” (2009), 10.1103/physrevlett.103.263901. Physical Review B 78 (2008), 10.1103/physrevb.78.085116. L. Razzari, A. Toma, M. Clerici, M. Shalaby, G. Das, C. Lib- D. Martin-Cano, O. Quevedo-Teruel, E. Moreno, L. Martin- erale, M. Chirumamilla, R. P. Zaccaria, F. D. Angelis, Moreno, and F. J. Garcia-Vidal, “Waveguided spoof sur- M. Peccianti, R. Morandotti, and E. D. Fabrizio, “Tera- 23 face plasmons with deep-subwavelength lateral conﬁne- R. Gordon and A. G. Brolo, “Increased cut-off wavelength ment,” Optics Letters 36, 4635 (2011). for a subwavelength hole in a real metal,” Optics Express 13, 1933 (2005). S. H. Lee, J. Choi, H.-D. Kim, H. Choi, and B. Min, “Ultra- J. Bravo-Abad, A. I. Fernandez-Dom ´ ´ ınguez, F. J. Garc´ ıa- fast refractive index control of a terahertz graphene meta- Vidal, and L. Mart´ ın-Moreno, “Theory of extraordinary material,” Scientiﬁc Reports 3 (2013), 10.1038/srep02135. transmission of light through quasiperiodic arrays of sub- D. Wang, Y. Gu, Y. Gong, C.-W. Qiu, and M. Hong, “An ul- wavelength holes,” Physical Review Letters 99 (2007), trathin terahertz quarter-wave plate using planar babinet- 10.1103/physrevlett.99.203905. inverted metasurface,” Optics Express 23, 11114 (2015). 167 Y. Takakura, “Optical resonance in a narrow slit in a thick S.-F. Shi, B. Zeng, H.-L. Han, X. Hong, H.-Z. Tsai, H. S. metallic screen,” Physical Review Letters 86, 5601 (2001). Jung, A. Zettl, M. F. Crommie, and F. Wang, “Op- V. Delgado and R. Marques, “Surface impedance model timizing broadband terahertz modulation with hybrid for extraordinary transmission in 1d metallic and dielec- graphene/metasurface structures,” Nano Letters 15, 372 tric screens,” Optics Express 19, 25290 (2011). (2014). 168 186 K. Takano, K. Shibuya, K. Akiyama, T. Nagashima, V. Galindo and C. Wu, “Numerical solutions for an inﬁnite F. Miyamaru, and M. Hangyo, “A metal-to-insulator tran- phased array of rectangular waveguides with thick walls,” sition in cut-wire-grid metamaterials in the terahertz re- IEEE Transactions on Antennas and Propagation 14, 149 gion,” Journal of Applied Physics 107, 024907 (2010). (1966). P. Sheng, R. S. Stepleman, and P. N. Sanda, “Exact eigen- A. Sommerfeld, “Mathematical theory of diffraction,” in functions for square-wave gratings: Application to diffrac- Mathematical Theory of Diffraction (Birkhuser Boston, 2004) tion and surface-plasmon calculations,” Physical Review B pp. 9–68. 26, 2907 (1982). C. J. R. Sheppard, J. Lin, and S. S. Kou, F. J. Garcia-Vidal, L. Martin-Moreno, T. W. Ebbesen, and “Rayleigh–sommerfeld diffraction formula in k space,” L. Kuipers, “Light passing through subwavelength aper- Journal of the Optical Society of America A 30, 1180 (2013). 171 tures,” Reviews of Modern Physics 82, 729 (2010). A. S. Marathay and J. F. McCalmont, “On the usual ap- H. Liu and P. Lalanne, “Microscopic theory of the extraor- proximation used in the rayleigh–sommerfeld diffraction dinary optical transmission,” Nature 452, 728 (2008). theory,” Journal of the Optical Society of America A 21, X.-Y. He, “Numerical analysis of the propagation proper- 510 (2004). ties of subwavelength semiconductor slit in the terahertz H. A. Bethe, “Theory of diffraction by small holes,” Physi- region,” Optics Express 17, 15359 (2009). cal Review 66, 163 (1944). 173 191 C. Bouwkamp and H. Casimir, “On multipole expansions P. Bell, J. Pendry, L. Moreno, and A. Ward, “A program for in the theory of electromagnetic radiation,” Physica 20, 539 calculating photonic band structures and transmission co- (1954). efﬁcients of complex structures,” Computer Physics Com- C. J. Bouwkamp, “Diffraction theory,” Reports on Progress munications 85, 306 (1995). in Physics 17, 35 (1954). L. Li, “New formulation of the fourier modal method for crossed surface-relief gratings,” Journal of the Optical So- R. Ulrich, “Far-infrared properties of metallic mesh and its ciety of America A 14, 2758 (1997). complementary structure,” Infrared Physics 7, 37 (1967). L. Salomon, F. Grillot, A. V. Zayats, and F. de Fornel, T. J. Kim, T. Thio, T. W. Ebbesen, D. E. Grupp, and H. J. “Near-ﬁeld distribution of optical transmission of periodic Lezec, “Control of optical transmission through metals subwavelength holes in a metal ﬁlm,” Physical Review perforated with subwavelength hole arrays,” Optics Let- Letters 86, 1110 (2001). ters 24, 256 (1999). J. Zuloaga, E. Prodan, and P. Nordlander, “Quantum D. S. Kim, S. C. Hohng, V. Malyarchuk, Y. C. Yoon, Y. H. description of the plasmon resonances of a nanoparticle Ahn, K. J. Yee, J. W. Park, J. Kim, Q. H. Park, and C. Lienau, dimer,” Nano Letters 9, 887 (2009). “Microscopic origin of surface-plasmon radiation in plas- R. Esteban, A. G. Borisov, P. Nordlander, and J. Aizpu- monic band-gap nanostructures,” Physical Review Letters rua, “Bridging quantum and classical plasmonics with 91 (2003), 10.1103/physrevlett.91.143901. a quantum-corrected model,” Nature Communications 3 B. K. Minhas, W. Fan, K. Agi, S. R. J. Brueck, and K. J. (2012), 10.1038/ncomms1806. Malloy, “Metallic inductive and capacitive grids: theory D. Marinica, A. Kazansky, P. Nordlander, J. Aizpurua, and and experiment,” Journal of the Optical Society of America A. G. Borisov, “Quantum plasmonics: Nonlinear effects in A 19, 1352 (2002). 179 the ﬁeld enhancement of a plasmonic nanoparticle dimer,” A. Naweed, F. Baumann, J. William A. Bailey, A. S. Nano Letters 12, 1333 (2012). Karakashian, and W. D. Goodhue, “Evidence for radiative C. Ciraci, R. T. Hill, J. J. Mock, Y. Urzhumov, damping in surface-plasmon-mediated light transmission A. I. Fernandez-Dominguez, S. A. Maier, J. B. Pendry, through perforated conducting ﬁlms,” Journal of the Opti- A. Chilkoti, and D. R. Smith, “Probing the ultimate lim- cal Society of America B 20, 2534 (2003). 180 its of plasmonic enhancement,” Science 337, 1072 (2012). D. Qu and D. Grischkowsky, “Observation of a new type M. S. Tame, K. R. McEnery, S ¸ . K. zdemir, J. Lee, S. A. Maier, of THz resonance of surface plasmons propagating on and M. S. Kim, “Quantum plasmonics,” Nature Physics 9, metal-ﬁlm hole arrays,” Physical Review Letters 93 (2004), 329 (2013). 10.1103/physrevlett.93.196804. 181 S. F. Tan, L. Wu, J. K. W. Yang, P. Bai, M. Bosman, and J. H. Kang, D. S. Kim, and Q.-H. Park, “Local capacitor C. A. Nijhuis, “Quantum plasmon resonances controlled model for plasmonic electric ﬁeld enhancement,” Physical by molecular tunnel junctions,” Science 343, 1496 (2014). Review Letters 102 (2009), 10.1103/physrevlett.102.093906. K. J. Savage, M. M. Hawkeye, R. Esteban, A. G. Borisov, J. Aizpurua, and J. J. Baumberg, “Revealing the quantum 24 regime in tunnelling plasmonics,” Nature 491, 574 (2012). tennas,” Optics Express 22, 30467 (2014). 201 218 J. A. Scholl, A. Garc´ ıa-Etxarri, A. L. Koh, and J. A. Dionne, O. Sushko, R. Dubrovka, and R. S. Donnan, “Terahertz “Observation of quantum tunneling between two plas- spectral domain computational analysis of hydration shell monic nanoparticles,” Nano Letters 13, 564 (2013). of proteins with increasingly complex tertiary structure,” The Journal of Physical Chemistry B 117, 16486 (2013). W. Zhu and K. B. Crozier, “Quantum mechanical limit to plasmonic enhancement as observed by surface-enhanced V. C. Nibali and M. Havenith, “New insights into the raman scattering,” Nature Communications 5 (2014), role of water in biological function: Studying solvated 10.1038/ncomms6228. biomolecules using terahertz absorption spectroscopy in W. Zhu, R. Esteban, A. G. Borisov, J. J. Baumberg, P. Nord- conjunction with molecular dynamics simulations,” Jour- lander, H. J. Lezec, J. Aizpurua, and K. B. Crozier, nal of the American Chemical Society 136, 12800 (2014). “Quantum mechanical effects in plasmonic structures with C. Rocchi, A. R. Bizzarri, and S. Cannistraro, “Water dy- subnanometre gaps,” Nature Communications 7 (2016), namical anomalies evidenced by molecular-dynamics sim- 10.1038/ncomms11495. ulations at the solvent-protein interface,” Physical Review B. de Nijs, F. Benz, S. J. Barrow, D. O. Sigle, R. Chikkaraddy, E 57, 3315 (1998). A. Palma, C. Carnegie, M. Kamp, R. Sundararaman, M. Marchi, F. Sterpone, and M. Ceccarelli, “Water ro- P. Narang, O. A. Scherman, and J. J. Baumberg, “Plas- tational relaxation and diffusion in hydrated lysozyme,” monic tunnel junctions for single-molecule redox chem- Journal of the American Chemical Society 124, 6787 (2002). istry,” Nature Communications 8 (2017), 10.1038/s41467- N. Sengupta, S. Jaud, and D. J. Tobias, “Hydration dynam- 017-00819-7. ics in a partially denatured ensemble of the globular pro- P. Cristofolini, G. Christmann, S. I. Tsintzos, G. Deligeor- tein human a-lactalbumin investigated with molecular dy- gis, G. Konstantinidis, Z. Hatzopoulos, P. G. Savvidis, and namics simulations,” Biophysical Journal 95, 5257 (2008). J. J. Baumberg, “Coupling quantum tunneling with cavity S. K. Sinha, S. Chakraborty, and S. Bandyopadhyay, photons,” Science 336, 704 (2012). “Thickness of the hydration layer of a protein from molec- J. D. Mar, X. L. Xu, J. J. Baumberg, A. C. Irvine, C. Stan- ular dynamics simulation,” The Journal of Physical Chem- ley, and D. A. Williams, “Voltage-controlled electron tun- istry B 112, 8203 (2008). neling from a single self-assembled quantum dot embed- Y. Xu, R. Gnanasekaran, and D. M. Leitner, “Analysis of ded in a two-dimensional-electron-gas-based photovoltaic water and hydrogen bond dynamics at the surface of an cell,” Journal of Applied Physics 110, 053110 (2011). antifreeze protein,” Journal of Atomic, Molecular, and Op- N. Cazier, M. Buret, A. V. Uskov, L. Markey, J. Arocas, tical Physics 2012, 1 (2012). G. C. D. Francs, and A. Bouhelier, “Electrical excitation of S. Chakraborty, S. K. Sinha, and S. Bandyopadhyay, “Low- waveguided surface plasmons by a light-emitting tunnel- frequency vibrational spectrum of water in the hydra- ing optical gap antenna,” Optics Express 24, 3873 (2016). tion layer of a protein: a molecular dynamics simulation A. V. Uskov, J. B. Khurgin, I. E. Protsenko, I. V. Smetanin, study,” The Journal of Physical Chemistry B 111, 13626 and A. Bouhelier, “Excitation of plasmonic nanoanten- (2007). nas by nonresonant and resonant electron tunnelling,” S. K. Sinha and S. Bandyopadhyay, “Differential ﬂexibility Nanoscale 8, 14573 (2016). of the secondary structures of lysozyme and the structure A. Stolz, J. Berthelot, M.-M. Mennemanteuil, G. C. des and ordering of surrounding water molecules,” The Jour- Francs, L. Markey, V. Meunier, and A. Bouhelier, “Non- nal of Chemical Physics 134, 115101 (2011). linear photon-assisted tunneling transport in optical gap 227 S. Bandyopadhyay, S. Chakraborty, and B. Bagchi, “Explo- antennas,” Nano Letters 14, 2330 (2014). ration of the secondary structure speciﬁc differential sol- V. Kravtsov, S. Berweger, J. M. Atkin, and M. B. Raschke, vation dynamics between the native and molten globule “Control of plasmon emission and dynamics at the transi- states of the protein HP-36,” The Journal of Physical Chem- tion from classical to quantum coupling,” Nano Letters 14, istry B 110, 20629 (2006). 5270 (2014). S. Bandyopadhyay, S. Chakraborty, S. Balasubramanian, S. C. H. Allen Taﬂove, Computational Electrodynamics: The S. Pal, and B. Bagchi, “Atomistic simulation study of Finite-Difference Time-Domain Method (ARTECH HOUSE the coupled motion of amino acid residues and water INC, 2005). molecules around protein HP-36: Fluctuations at and J. D. Jackson, Classical electrodynamics, Third Edition (John around the active sites,” The Journal of Physical Chemistry Wiley & Sons, 1998). B 108, 12608 (2004). D. J. Grifﬁths, Introduction to electrodynamics; 4th ed. (Pear- 229 S. Pal and S. Bandyopadhyay, “Effects of protein confor- son, Boston, MA, 2013) re-published by Cambridge Uni- mational ﬂexibilities and electrostatic interactions on the versity Press in 2017. low-frequency vibrational spectrum of hydration water,” C. Fattinger and D. Grischkowsky, “Point source terahertz The Journal of Physical Chemistry B 117, 5848 (2013). optics,” Applied Physics Letters 53, 1480 (1988). 230 S. Pal and S. Bandyopadhyay, “Importance of protein con- C. Fattinger and D. Grischkowsky, “Terahertz beams,” Ap- formational motions and electrostatic anchoring sites on plied Physics Letters 54, 490 (1989). the dynamics and hydrogen bond properties of hydration D. Grischkowsky, S. Keiding, M. van Exter, and C. Fat- water,” Langmuir 29, 1162 (2013). tinger, “Far-infrared time-domain spectroscopy with tera- S. K. Sinha and S. Bandyopadhyay, “Polar solvation dy- hertz beams of dielectrics and semiconductors,” Journal of namics of lysozyme from molecular dynamics studies,” the Optical Society of America B 7, 2006 (1990). The Journal of Chemical Physics 136, 185102 (2012). S. J. Park, B. H. Son, S. J. Choi, H. S. Kim, and Y. H. S. K. Sinha and S. Bandyopadhyay, “Local heterogeneous Ahn, “Sensitive detection of yeast using terahertz slot an- dynamics of water around lysozyme: a computer simula- 25 tion study,” Phys. Chem. Chem. Phys. 14, 899 (2012). H. liang Zhao, G. jun Ren, F. Liu, H. peng Xin, Y. bo Bai, and J. quan Yao, “Tunable terahertz source via liquid crys- T. Ding, T. Huber, A. P. Middelberg, and R. J. Fal- tal grating coated with electron beam excited graphene: coner, “Characterization of low-frequency modes in aque- A theoretical analysis,” Optics Communications 390, 137 ous peptides using far-infrared spectroscopy and molecu- (2017). lar dynamics simulation,” The Journal of Physical Chem- X. Gu, I.-T. Lin, and J.-M. Liu, “Extremely conﬁned tera- istry A 115, 11559 (2011). hertz surface plasmon-polaritons in graphene-metal struc- M. Heyden and M. Havenith, “Combining THz spec- tures,” Applied Physics Letters 103, 071103 (2013). troscopy and MD simulations to study protein-hydration H. Zhao, Q. Guo, F. Xia, and H. Wang, “Two-dimensional coupling,” Methods 52, 74 (2010). materials for nanophotonics application,” Nanophotonics M. Heyden, J. Sun, H. Forbert, G. Mathias, M. Havenith, 4 (2015), 10.1515/nanoph-2014-0022. and D. Marx, “Understanding the origins of dipolar cou- J. Conde, C. Bao, D. Cui, P. V. Baptista, and F. Tian, “Anti- plings and correlated motion in the vibrational spectrum body–drug gold nanoantennas with raman spectroscopic of water,” The Journal of Physical Chemistry Letters 3, ﬁngerprints for in vivo tumour theranostics,” Journal of 2135 (2012). Controlled Release 183, 87 (2014). M. P. Allen et al., “Introduction to molecular dynamics sim- X. Huang and M. A. El-Sayed, “Gold nanoparticles: Op- ulation,” Computational soft matter: from synthetic poly- tical properties and implementations in cancer diagnosis mers to proteins 23, 1 (2004). and photothermal therapy,” Journal of Advanced Research J. E. Bertie and Z. Lan, “Infrared intensities of liquids XX: 1, 13 (2010). The intensity of the OH stretching band of liquid water re- A. Abbas, M. J. Linman, and Q. Cheng, “New trends in visited, and the best current values of the optical constants instrumental design for surface plasmon resonance-based of h2o(l) at 25°c between 15,000 and 1 cm-1,” Applied Spec- biosensors,” Biosensors and Bioelectronics 26, 1815 (2011). troscopy 50, 1047 (1996). 238 Y. Yanase, T. Hiragun, K. Ishii, T. Kawaguchi, T. Yanase, M. Heyden, J. Sun, S. Funkner, G. Mathias, H. Forbert, M. Kawai, K. Sakamoto, and M. Hide, “Surface plasmon M. Havenith, and D. Marx, “Dissecting the THz spectrum resonance for cell-based clinical diagnosis,” Sensors 14, of liquid water from ﬁrst principles via correlations in time 4948 (2014). and space,” Proceedings of the National Academy of Sci- E. Helmerhorst, D. J. Chandler, M. Nussio, and C. D. ences 107, 12068 (2010). 239 Mamotte, “Real-time and label-free bio-sensing of molec- J. Jiang, M. Xiao, S. Wang, K. Liu, X. Wang, and T. Liu, ular interactions by surface plasmon resonance: A labo- “Polarized low-coherence interferometer based on a ma- ratory medicine perspective,” The Clinical biochemist. Re- trix CCD and birefringence crystal with a two-dimensional views 33, 161 (2012). angle,” Optics Express 25, 15977 (2017). S. L. Dodson, C. Cao, H. Zaribafzadeh, S. Li, and Q. Xiong, “Engineering plasmonic nanorod arrays for colon cancer marker detection,” Biosensors and Bioelectronics 63, 472 (2015).
arXiv (Cornell University)