Colossal Absorption of Molecules Inside Single Terahertz

Feb 26, 2013 - Brown , E. R.; Bjarnason , J. E.; Fedor , A. M.; Korter , T. M. Appl. Phys. Lett. 2007, 90 (6) ..... Daniel M. Mittleman. Nature Photon...
1 downloads 0 Views 326KB Size
Letter pubs.acs.org/NanoLett

Colossal Absorption of Molecules Inside Single Terahertz Nanoantennas Hyeong-Ryeol Park,† Kwang Jun Ahn,‡ Sanghoon Han,§ Young-Mi Bahk,† Namkyoo Park,§ and Dai-Sik Kim*,† †

Center for Subwavelength Optics, Department of Physics and Astronomy, Seoul National University, Seoul 151-747, Korea Global Frontier Center for Multiscale Energy Systems, Seoul National University, Seoul 151-742, Korea § Photonic Systems Laboratory, School of EECS, Seoul National University, Seoul 151-744, Korea ‡

S Supporting Information *

ABSTRACT: Molecules have extremely small absorption cross sections in the terahertz range even under resonant conditions, which severely limit their detectability, often requiring tens of milligrams. We demonstrate that nanoantennas tailored for the terahertz range resolves the small molecular cross section problem. The extremely asymmetric electromagnetic environment inside the slot antenna, which finds the electric field being enhanced by thousand times with the magnetic field changed little, forces the molecular cross section to be enhanced by >103 accompanied by a colossal absorption coefficient of ∼170 000 cm−1. Tens of nanograms of small molecules such as 1,3,5trinitroperhydro-1,3,5-triazine (RDX) and lactose drop-cast over an area of 10 mm2, with only tens of femtograms of molecules inside the single nanoslot, can readily be detected. Our work enables terahertz sensing of chemical and biological molecules in ultrasmall quantities. KEYWORDS: Terahertz nano, terahertz molecule sensing, terahertz absorption, terahertz absorption coefficient, terahertz cross section, terahertz nanoantenna electric field E, and magnetic field H, the Fermi’s golden rule states that the energy loss per unit time is proportional to μ2E2, which depletes from the incident energy flux S⃗ = E⃗ × H⃗ * (= EH). Therefore, the quantum mechanical absorption cross section σ can be written as follows:36

T

he demand for new chemical and biological sensing has been dramatically increasing in optical,1−3 infrared,4,5 terahertz (THz),6−8 or microwave9 regimes. Especially, the emergence of diverse materials, including organic- and biomolecules that respond resonantly to the THz waves, raises the possibility of sensitive THz detection.10−12 For improving sensitivities, various THz molecule sensing techniques utilizing parallel plate waveguide,13,14 microresonator,15 and frequency filter,16−19 have been developed, mostly by increasing the effective interaction length.20 Recently, plasmonic nanoantennas,21−23 which support collective resonances with strong local field in the infrared frequencies,24 have been suggested as a tool for the molecular detection, enhancing absorption signals of molecules around the antennas by 4−5 orders of magnitude.25−27 However, these tools have intrinsic limitations to measure the absolute absorptivity and cross section of their respective molecules due to the greater quantitative uncertainties of the reacting molecules in and around the positive-mold antennas and the lack of simple measure for energy and field enhancements. Meanwhile negative slot antennas28 allow simple and precise analysis of the electromagnetic fields over the whole aperture, both in terms of energy and field enhancements.29−35 When a molecule with a dipole moment μ is subject to an electromagnetic wave with intensity (Poynting vector) S, © 2013 American Chemical Society

σ=

ℏω 2π 2 2 2π 2 E2 μ E ρ(ℏω0) 0 = μ ρ(ℏω0)ℏω0 ℏ S ℏ S

(1)

where ρ(ℏω0) is the density of states, ω0 and S are the resonant angular frequency and the Poynting vector of the incident light, respectively (See Supporting Information for further details). In free space, (E2/S) = (E/H) = Z0 (Z0 = 377 Ω), and therefore σ = (2π/ℏ)μ2ρ(ℏω0)ℏω0Z0. Inside a terahertz nanoslot antenna (Figure 1a),37−39 however, the case is strikingly different although eq 1 is still valid. Dictated by the boundary condition of the metal slot, the electromagnetic wave is strongly perturbed, and the resultant electric field is enhanced by orders of magnitudes, in contrast to the magnetic field which deviates little from the incident value, within a factor of two.40,41 This extremely asymmetric electromagnetic environment in turn forces the cross section of the molecules to Received: January 29, 2013 Revised: February 20, 2013 Published: February 26, 2013 1782

dx.doi.org/10.1021/nl400374z | Nano Lett. 2013, 13, 1782−1786

Nano Letters

Letter

Figure 1. (a) The giant absorption cross section of RDX (1,3,5-trinitroperhydro-1,3,5-triazine) molecules inside a THz nanoslot antenna is induced by an asymmetric electromagnetic environment with E/Z0H ratio of about 103. The p-polarized terahertz pulses with polarization perpendicular to the long axis of the rectangle are normally incident on the sample. (b) Absorption spectrum using 1 mg of RDX molecules onto the bare quartz substrate. (c) Transmission spectra for RDX amounts ranging from 2 μg to 200 μg onto the bare quartz. (d) Transmission spectra for a single THz nanoslot antenna, having dimensions of the length 90 μm, the width 50 nm, and the metal thickness 50 nm, without and with RDX ranging from 40 ng to 2 μg.

Figure 2. (a, Left) Schematic of a single THz nanoslot antenna with a length of l = 90 μm, a width of w, and a thickness of h. (Right) SEM images of some of our slot samples with w = 50, 500, 1000, and 5000 nm from left to right. (b) Normalized resonant transmissions measured at 0.87 THz for all the slots: w = 50 nm (h = 50 nm), 100, 200, 500, 1000, and 5000 nm (h = 100 nm), with the RDX amounts ranging from 0 to 2 μg. (c) Peak absorption coefficients and cross sections of RDX molecules inside slots, plotted against the width w (red squares). The fit (black dashed line) indicates a 1/w dependence. (d) The schematic diagram of the molecular cross section characteristics for three types of asymmetric electromagnetic environments of E/Z0H = 1, 103, and 105.

1783

dx.doi.org/10.1021/nl400374z | Nano Lett. 2013, 13, 1782−1786

Nano Letters

Letter

Figure 3. (a) SEM images of two nanoslot antennas with different lengths of 90 and 150 μm, having resonances of 0.9 THz and 0.5 THz, respectively. (b) Transmission spectra through the two nanoslot antennas without and with RDX amounts ranged from 0 to 2 μg. (c) Absorption spectrum using 1.5 mg lactose molecules onto the bare substrate. The peak absorption at 0.54 THz is well-matched to the fundamental resonance of the 150 μm length nanoslot antenna. (d) Transmission spectra for the 150 μm length nanoslot antenna with lactose amounts ranged from 0 to 1 μg.

enhancement factor follows ∑ = (E2/Z0S) ∝ (l/w), as expected from an analytical calculation based on modal expansion26,29 (for more details, see Supporting Information), we vary the nanoantenna width. Single THz slot antennas with fixed l = 90 μm and varying widths (w) in thin gold films with the thickness of h were fabricated using a focused ion beam machine on a 500 μm thick quartz substrate, as shown in Figure 2a (left), with their resonance peaks at around 0.87 THz. Figure 2a (right) presents SEM images of some of our slot samples with w = 50, 500, 1000, and 5000 nm from left to right. The RDX molecules were filled inside the nanoslot by simple drop casting over an area of ∼8 mm2. The solutions were diluted from the RDX standard solution (Accustandard, Inc.) containing 1 mg/mL of RDX in acetonitrile/methanol.13 It is emphasized that although the RDX solutions from 0.02 to 2 μL were drop-cast over the said area, actually participating are only a few tens to hundreds of femtograms inside the slot. Figure 2b shows normalized resonant transmissions measured at 0.87 THz for the slots: w = 50 nm (h = 50 nm), 100, 200, 500, 1000, and 5000 nm (h = 100 nm), with the increasing RDX amounts up to 2 μg. For most of the samples, the transmission decrease saturates at ∼0.1 μg of RDX, consistent with the RDX film thickness t ∼ h. From these data, using the Beer’s law, we obtain the absorption coefficient α thereby the absorption cross section σ through α = σN = −(1/t)ln(T/T0). The number density of dried molecules is estimated as N = 4.93 × 1021cm−3, and T and T0 are the transmissions with and without RDX molecules, respectively. Figure 2c plots absorption coefficient and cross sections as a function of the antenna width w. The experimentally obtained α and σ are inversely proportional to w as expected, and are a few thousand times larger than the bulk values shown in Figure 1b, also consistent with what we predict from the ∑ ∝(l/w) dependence using our theoretical calculation.

become larger by a huge factor ∑ = (E2/Z0S) = (E/Z0H), at least 3 orders of magnitude, resulting in the absorption coefficient α reaching the colossal value of >100 000 cm−1. To probe absorption cross sections of molecules in THz region, we perform THz time-domain spectroscopy40,42,43 for the frequencies ranging from 0.3 to 2.0 THz. In Figure 1b, it was found that the drop-cast 1 mg RDX (1,3,5-trinitroperhydro-1,3,5-triazine) molecules over an area of ∼80 mm2 on the bare quartz have the most prominent absorption peak near 0.83 THz, in agreement with the previous free space THz measurements.44,45 It should be noted, as shown in Figure 1c, that we really do need 1 mg of RDX to get meaningful decrease in transmission- there is little change between 2 μg to 200 μg on the bare quartz. For observing small quantities of the molecules as promised by the unique electromagnetic environments of the slot antenna, we first designed a single THz nanoantenna with a large length/width ratio of l/w = 90 μm/50 nm =1800, to have 3 orders of magnitudes of electric field enhancement.29 The resonance wavelength of the antenna satisfies the condition λres = 2neffl, where neff is the effective index of refraction of the antenna−substrate composite.29,31 The length is designed to match one of the RDX absorption bands, about 0.8 THz, attributed to molecular conformations or a weak hydrogen bond between two RDX molecules.46,47 The contribution of the direct transmission through the unpatterned gold film is subtracted when obtaining the transmission spectra through the single nanoantennas. Strikingly, when we put RDX molecules into the metallic substrate with tailored nanoslot in the middle, as shown in Figure 1a, we start to see significant decrease in the resonant transmission intensity already for the amount of 40 ng spread over a sample area of 8 mm2, amounting to 22 fg inside the slot cavity (Figure 1d). For the amount of 2 μg RDX molecules, the resonance is almost gone. To see if the cross section 1784

dx.doi.org/10.1021/nl400374z | Nano Lett. 2013, 13, 1782−1786

Nano Letters

Letter

The maximum α and σ are, respectively, 170 000 cm−1 and 3.6 × 10−17 cm2 for the 50 nm width slot, about 2800 times larger than the values on the bare substrate. Hereby, it is possible to estimate how much more energy a single molecule inside the slot absorbs compared with the bare substrate. Since the absorbed power per molecule is σS, we need to approximate the Poynting vector enhancement as well: from the electric field enhancement of 2500 experimentally measured through the Kirchhoff formalism39,40,48 (see Supporting Information for further detail), we estimate the Poynting vector enhancement of ∼2900.29,49 The molecular absorption enhancement is therefore 2800 × 2900 ≈ 8 × 106. This is also in good quantitative agreement with the absorption enhancements estimated by PEC calculations and FDTD simulations (see Supporting Information). As made evident by these findings, because electric field E is much more enhanced than magnetic field Z0H inside the slot, it should be possible to amplify the molecular cross section inside the slot antenna further by closing the width even more (Figure 2d, middle). In this situation, only the ratio of E to Z0H, not the absolute power of the excited light is important to modulate the cross section. Using an ultimate nanoslot with a width of below 1 nm, it would be possible to reach 105 of the E/Z0H ratio necessary to induce the enhanced molecular cross section to approach the actual molecular size of ∼10−14 cm2 = 1 nm2 (Figure 2d, right). We now explore the possibility of extending our method to encompass any molecules with resonances in terahertz. We manufacture two single THz nanoslot antennas with the same width of 120 nm, same thickness of 100 nm, and different lengths l = 90 and 150 μm, having the resonances at 0.9 THz and 0.5 THz, respectively (Figure 3a). In stark contrast to the 90 μm length slot antenna, spectra for the 150 μm length slot antenna change little, as shown in Figure 3b. We put lactose molecules, which have substantial cross section near 0.5 THz (Figure 3c),50 into the 150 μm length slot antenna. The change in transmission is as dramatic as the RDX case (Figure 3d), and we deduce an enhanced absorption coefficient of 164,000 cm−1 demonstrating the generality of our terahertz nanoslot technology. In summary, we have demonstrated the molecular cross section and absorption coefficient enhancement created by the asymmetric electromagnetic environment with E/Z0H ratio of over 103 inside nanoslots. Up to 2800 enhancement of absorption cross section of RDX molecules, leading to the giant molecular absorption enhancement of 8 × 106, enables detection of extremely small quantities near 40 ng (22 fg inside the slot). With orders of magnitudes improved sensitivities and demonstrated generality, our method could pave the way toward THz intermolecular fingerprinting of target molecules, as well as ultrasmall quantity THz spectroscopy of chemicals and biological molecules in their developmental stages.



simulations. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MEST) (SRC, No: 2008-0062255) (GRL, No: K20815000003) (Others, No: 2010-0029648, 2011-0019170, 2011-0020209), and Hi Seoul Science/Humanities Fellowship from Seoul Scholarship Foundation. We thank J. Y. Rhie for technical assistance with sample fabrication.



REFERENCES

(1) Zhu, J. G.; Ozdemir, S. K.; Xiao, Y. F.; Li, L.; He, L. N.; Chen, D. R.; Yang, L. Nat. Photon. 2010, 4 (1), 46−49. (2) Carmeli, I.; Lieberman, I.; Kraversky, L.; Fan, Z. Y.; Govorov, A. O.; Markovich, G.; Richter, S. Nano Lett. 2010, 10 (6), 2069−2074. (3) Fang, Z. Y.; Liu, Z.; Wang, Y. M.; Ajayan, P. M.; Nordlander, P.; Halas, N. J. Nano Lett. 2012, 12 (7), 3808−3813. (4) Osawa, M. B. Chem. Soc. Jpn. 1997, 70 (12), 2861−2880. (5) Cubukcu, E.; Zhang, S.; Park, Y. S.; Bartal, G.; Zhang, X. Appl. Phys. Lett. 2009, 95 (4), 043113. (6) Fitch, M. J.; Leahy-Hoppa, M. R.; Ott, E. W.; Osiander, R. Chem. Phys. Lett. 2007, 443 (4−6), 284−288. (7) Huber, A. J.; Keilmann, F.; Wittborn, J.; Aizpurua, J.; Hillenbrand, R. Nano Lett. 2008, 8 (11), 3766−3770. (8) Lim, M.; Choi, S. J.; Lee, G. S.; Seol, M. L.; Do, Y.; Choi, Y. K.; Han, H. Appl. Phys. Lett. 2012, 100 (21), 211102. (9) Liu, T. M.; Chen, H. P.; Wang, L. T.; Wang, J. R.; Luo, T. N.; Chen, Y. J.; Liu, S. I.; Sun, C. K. Appl. Phys. Lett. 2009, 94 (4), 043902. (10) Fischer, B. M.; Hoffmann, M.; Helm, H.; Wilk, R.; Rutz, F.; Kleine-Ostmann, T.; Koch, M.; Jepsen, P. U. Opt. Express 2005, 13 (14), 5205−5215. (11) Shi, Y. L.; Wang, L. J. Phys. D: Appl. Phys. 2005, 38 (19), 3741− 3745. (12) Kawase, K.; Ogawa, Y.; Watanabe, Y.; Inoue, H. Opt. Express 2003, 11 (20), 2549−2554. (13) Melinger, J. S.; Laman, N.; Grischkowsky, D. Appl. Phys. Lett. 2008, 93 (1), 011102. (14) Lee, E. S.; Lee, S. G.; Kee, C. S.; Jeon, T. I. Opt. Express 2011, 19 (16), 14852−14859. (15) Nagel, M.; Bolivar, P. H.; Brucherseifer, M.; Kurz, H.; Bosserhoff, A.; Buttner, R. Appl. Opt. 2002, 41 (10), 2074−2078. (16) Tiang, C. K.; Cunningham, J.; Wood, C.; Hunter, I. C.; Davies, A. G. J. Appl. Phys. 2006, 100 (6), 066105. (17) Ren, L.; Pint, C. L.; Booshenri, L. G.; Rice, W. D.; Wang, X. F.; Hilton, D. J.; Takeya, K.; Kawayama, I.; Tonouchi, M.; Hauge, R. H.; Kono, J. Nano Lett. 2009, 9 (7), 2610−2613. (18) Kyoung, J.; Jang, E. Y.; Lima, M. D.; Park, H. R.; Robles, R. O.; Lepro, X.; Kim, Y. H.; Baughman, R. H.; Kim, D. S. Nano Lett. 2011, 11 (10), 4227−4231. (19) Ren, L.; Pint, C. L.; Arikawa, T.; Takeya, K.; Kawayama, I.; Tonouchi, M.; Hauge, R. H.; Kono, J. Nano Lett. 2012, 12 (2), 787− 790. (20) Choi, M.; Lee, S. H.; Kim, Y.; Kang, S. B.; Shin, J.; Kwak, M. H.; Kang, K. Y.; Lee, Y. H.; Park, N.; Min, B. Nature 2011, 470 (7334), 369−373. (21) Heshmat, B.; Pahlevaninezhad, H.; Pang, Y. J.; Masnadi-Shirazi, M.; Lewis, R. B.; Tiedje, T.; Gordon, R.; Darcie, T. E. Nano Lett. 2012, 12 (12), 6255−6259.

ASSOCIATED CONTENT

S Supporting Information *

Derivation of quantum mechanical absorption cross section. Terahertz time-domain measurement and estimation of field enhancements. Analytical calculation based on modal expansion with slot antenna array. Molecular absorption enhancements inside slot antenna array in a PEC film estimated by the analytical model. Molecular absorption enhancements inside slot antenna array in a real metal film estimated by the FDTD 1785

dx.doi.org/10.1021/nl400374z | Nano Lett. 2013, 13, 1782−1786

Nano Letters

Letter

(50) Brown, E. R.; Bjarnason, J. E.; Fedor, A. M.; Korter, T. M. Appl. Phys. Lett. 2007, 90 (6), 061908.

(22) Knight, M. W.; Liu, L. F.; Wang, Y. M.; Brown, L.; Mukherjee, S.; King, N. S.; Everitt, H. O.; Nordlander, P.; Halas, N. J. Nano Lett. 2012, 12 (11), 6000−6004. (23) Seok, T. J.; Jamshidi, A.; Kim, M.; Dhuey, S.; Lakhani, A.; Choo, H.; Schuck, P. J.; Cabrini, S.; Schwartzberg, A. M.; Bokor, J.; Yablonovitch, E.; Wu, M. C. Nano Lett. 2011, 11 (7), 2606−2610. (24) Schuller, J. A.; Barnard, E. S.; Cai, W. S.; Jun, Y. C.; White, J. S.; Brongersma, M. L. Nat. Mater. 2010, 9 (3), 193−204. (25) Adato, R.; Yanik, A. A.; Amsden, J. J.; Kaplan, D. L.; Omenetto, F. G.; Hong, M. K.; Erramilli, S.; Altug, H. Proc. Natl. Acad. Sci. U.S.A. 2009, 106 (46), 19227−19232. (26) Wu, C. H.; Khanikaev, A. B.; Adato, R.; Arju, N.; Yanik, A. A.; Altug, H.; Shvets, G. Nat. Mater. 2012, 11 (1), 69−75. (27) Le, F.; Brandl, D. W.; Urzhumov, Y. A.; Wang, H.; Kundu, J.; Halas, N. J.; Aizpurua, J.; Nordlander, P. ACS Nano 2008, 2 (4), 707− 718. (28) Blumlein, A. D. Improvements in or relating to high frequency electrical conductors or radiators. British Patent 515683, 1938. (29) Garcia-Vidal, F. J.; Moreno, E.; Porto, J. A.; Martin-Moreno, L. Phys. Rev. Lett. 2005, 95 (10), 103901. (30) Seo, M. A.; Adam, A. J.; Kang, J. H.; Lee, J. W.; Ahn, K. J.; Park, Q. H.; Planken, P. C.; Kim, D. S. Opt. Express 2008, 16 (25), 20484− 20489. (31) Park, H. R.; Koo, S. M.; Suwal, O. K.; Park, Y. M.; Kyoung, J. S.; Seo, M. A.; Choi, S. S.; Park, N. K.; Kim, D. S.; Ahn, K. J. Appl. Phys. Lett. 2010, 96 (21), 211109. (32) Chen, H. T.; Lu, H.; Azad, A. K.; Averitt, R. D.; Gossard, A. C.; Trugman, S. A.; O’Hara, J. F.; Taylor, A. J. Opt. Express 2008, 16 (11), 7641−7648. (33) Seo, M.; Kyoung, J.; Park, H.; Koo, S.; Kim, H. S.; Bernien, H.; Kim, B. J.; Choe, J. H.; Ahn, Y. H.; Kim, H. T.; Park, N.; Park, Q. H.; Ahn, K.; Kim, D. S. Nano Lett. 2010, 10 (6), 2064−2068. (34) Park, S. G.; Jin, K. H.; Vi, M.; Ye, J. C.; Ahn, J.; Jeong, K. H. ACS Nano 2012, 6 (3), 2026−2031. (35) Kang, J. H.; Kim, D. S.; Park, Q. H. Phys. Rev. Lett. 2009, 102 (9), 093906. (36) Svelto, O. Principles of lasers, 5th ed.; Springer: New York, 2010. (37) Park, H. R.; Park, Y. M.; Kim, H. S.; Kyoung, J. S.; Seo, M. A.; Park, D. J.; Ahn, Y. H.; Ahn, K. J.; Kim, D. S. Appl. Phys. Lett. 2010, 96 (12), 121106. (38) Bahk, Y. M.; Park, H. R.; Ahn, K. J.; Kim, H. S.; Ahn, Y. H.; Kim, D. S.; Bravo-Abad, J.; Martin-Moreno, L.; Garcia-Vidal, F. J. Phys. Rev. Lett. 2011, 106 (1), 013902. (39) Park, H. R.; Bahk, Y. M.; Ahn, K. J.; Park, Q. H.; Kim, D. S.; Martin-Moreno, L.; Garcia-Vidal, F. J.; Bravo-Abad, J. ACS Nano 2011, 5 (10), 8340−8345. (40) Seo, M. A.; Park, H. R.; Koo, S. M.; Park, D. J.; Kang, J. H.; Suwal, O. K.; Choi, S. S.; Planken, P. C. M.; Park, G. S.; Park, N. K.; Park, Q. H.; Kim, D. S. Nat. Photon. 2009, 3 (3), 152−156. (41) Garcia-Vidal, F. J.; Martin-Moreno, L.; Ebbesen, T. W.; Kuipers, L. Rev. Mod. Phys. 2010, 82 (1), 729−787. (42) Vanexter, M.; Fattinger, C.; Grischkowsky, D. Opt. Lett. 1989, 14 (20), 1128−1130. (43) Wu, Q.; Litz, M.; Zhang, X. C. Appl. Phys. Lett. 1996, 68 (21), 2924−2926. (44) Chen, J.; Chen, Y. Q.; Zhao, H. W.; Bastiaans, G. J.; Zhang, X. C. Opt. Express 2007, 15 (19), 12060−12067. (45) Yamamoto, K.; Yamaguchi, M.; Miyamaru, F.; Tani, M.; Hangyo, M.; Ikeda, T.; Matsushita, A.; Koide, K.; Tatsuno, M.; Minami, Y. J. Appl. Phys. B 2004, 43 (3B), L414−L417. (46) Huang, F.; Schulkin, B.; Altan, H.; Federici, J. F.; Gary, D.; Barat, R.; Zimdars, D.; Chen, M.; Tanner, D. B. Appl. Phys. Lett. 2004, 85 (23), 5535−5537. (47) Boyd, S. G.; Boyd, K. J. J. Chem. Phys. 2008, 129 (13), 1. (48) Kyoung, J. S.; Seo, M. A.; Park, H. R.; Ahn, K. J.; Kim, D. S. Opt. Commun. 2010, 283 (24), 4907−4910. (49) Kang, J. H.; Choe, J. H.; Kim, D. S.; Park, Q. H. Opt. Express 2009, 17 (18), 15652−15658. 1786

dx.doi.org/10.1021/nl400374z | Nano Lett. 2013, 13, 1782−1786