2. Two-Dimensional Positioning of a Nanoblock
The proposed 3D plasmonic ruler is composed of three silver nanoblocks. In order to readily illustrate the plasmonic interactions between the three blocks, we first investigated the optical properties of simpler structures that consist of two silver nanoblocks, called a 2D plasmonic ruler. After the determination of the lateral movement of two blocks, we extend the concept of the ruler to three silver nanoblocks.
Figure 1 shows a schematic diagram of a 2D plasmonic ruler that consists of a large silver nanoblock and a small nanoblock, where the blocks are spaced with an air gap of 10 nm. The upper small nanoblock and the lower large nanoblock have dimensions of 200 nm × 150 nm × 100 nm and 260 nm × 210 nm × 100 nm, respectively. In this structure, the lateral movement along the x–y plane of the upper block is estimated by observing the absorption spectra, and the lower block with a fixed position is declared as the reference block. Silver is an optimal material because of its lower absorption losses in noble metals for the spectral range of interest. The center block size was chosen to locate the spectral positions of plasmon resonant mode 1 (M1) and mode 2 (M2) near 1000 nm. The peak wavelengths are insensitive to the upper and lower block sizes. The given differences in
x and
y dimensions of the center block were determined for two peaks, M1 and M2, which were placed with enough spectral spacing to be larger than linewidths of the peaks.
In this paper, we calculate higher-order plasmonic resonant modes and their absorption cross-sections using a 3D finite-difference time-domain (FDTD) method. Silver was modeled by the Drude model, which is represented as follows:
The background dielectric constant (
), plasma frequency (
ωp), and collision frequency (
γ) were set to 3.14, 9.13, and 0.021 eV, respectively, by fitting experimentally-determined dielectric functions of silver [
14]. The spatial grid along the
x,
y and
z axes were 2.5, 2.5, and 1 nm, respectively.
When an incident light, linearly polarized along the
z-axis (
Ez), is injected into the proposed structure, a strongly confined surface plasmon cavity mode can be excited in the air gaps between nanoblocks. In this case, there are two orthogonal, higher-order plasmon resonant modes with one intensity node along the
x- or
y-axis. Moreover, compared with a fundamental gap cavity mode, the proposed cavity modes have one intensity node that suppresses radiation loss along the node direction [
15,
16]. Suppression of radiation loss provides a narrow linewidth, which can enhance the sensitivity of the gap size sensors [
15,
16]. Depending on the direction of the incident light, one of two modes can be selected by mode matching with an incident plane wave.
When the linearly polarized light is injected into the proposed structure along the
y-axis or
x-axis, a plasmon resonant mode 1 (M1) with an intensity node along the
x-axis or a mode 2 (M2) with a node along the
y-axis can be excited in the air gap between the two nanoblocks, as shown in
Figure 2a,b. To estimate the lateral movement (Δ
x and Δ
y) of the upper block, we need to exploit the optical properties of the two plasmonic modes, M1 and M2, simultaneously. In order to excite M1 and M2 at the same time, the incident angle of light was set to 40° from the
y-axis (see inset of
Figure 1). Our modes M1 and M2 have strong electric field in the air gap, not around sharp block edges. Additionally, as shown in the mode profiles of
Figure 2, significant enhancements at the edges are not observed. In fact, the resonant wavelengths of M1 and M2 are maintained, except for slight blue-shift due to the decrease of the effective cavity size when the corner of the blocks are rounded. In other words, the proposed plasmonic ruler is robust to the rounding of edges in the realistic experimental nanoparticles.
In spite of the similarity in the mode profiles of M1 and M2, the resonant wavelengths of the two modes are different because of the different
x- and
y-directional sizes of the rectangular cavity: 200 and 150 nm.
Figure 2a,b present a top-view of the dominant electric field (
Ez) mode profiles of M1 and M2 when the upper block is placed just above the center of the lower large block, where the profiles are obtained at the center of the air gap. The resonant wavelengths of M1 and M2 are 882.2 and 1081.3 nm, respectively. M2, with an intensity node along the
y-axis, primarily oscillates along the
x-axis (the longer axis) and has a longer wavelength than that of M1.
Next, we investigate changes in the properties of each mode for a lateral movement of the upper block when the air gap between the upper and lower blocks is fixed. As shown in
Figure 2c,d, when the upper block moves on the
x- and
y-axes by amounts of Δ
x = 30 nm and Δ
y = 30 nm from the original position, the position of the strongly confined field also changes with respect to the upper block position. Therefore, an effective cavity size that determines a resonant wavelength of a cavity mode can be changed because the electric field of the cavity mode partially experiences the placement of the air region out of the lower block. In this regard, M1 and M2 experience distinct resonant wavelength shifts depending on the direction and size of the lateral displacements, Δ
x or Δ
y, of the upper block. If the changes of the resonant wavelengths of M1 and M2 are measured, the direction and size of the lateral movement of the upper block are also obtained. As a result, the proposed structure can operate as a plasmonic ruler, which enables measurement of the horizontal relative position of the upper block to the lower block.
In order to quantify the change in the resonant wavelength of each mode for the lateral movement of the upper block, we calculated the resonant wavelengths of M1 and M2 as functions of the movement of the block along the
x- or
y-direction from 0 to 30 nm at 5 nm increments.
Figure 3a shows a schematic diagram of the upper block movement in the
x- and
y-directions. In this case, the air gap thickness between the upper and lower blocks was fixed to 10 nm. First, in the case of M1, as shown in
Figure 3b, distinctive resonant wavelength shifts are observed depending on the moving direction,
x- (black) or
y- (red). When the block moves along the
x-direction from 0 to 30 nm, the resonant wavelength of the M1 slightly increases from 882.2 to 885.1 nm. In contrast, when the block moves along the
y-direction from 0 to 30 nm, the resonant wavelength of M1 largely decreases from 882.2 to 865.1 nm. Hence, the wavelength of M1, with an intensity node along the
x-axis, has a larger wavelength shift for the movement of the upper block along the
y-axis.
On the other hand, the resonant wavelength of M2 shows opposite dependences on the movement of the block. As demonstrated in
Figure 3c, when the upper block moves along the
x-direction from 0 nm to 30 nm (black), the resonant wavelength largely blue-shifts from 1081.3 to 1064.7 nm; when the upper block moves along the
y-direction (red), the wavelength slightly red-shifts from 1081.3 to 1084.2 nm. The wavelength of M2 with a node along the
y-axis has a larger wavelength shift for the movement of the block along the
x-axis. That is to say, when the upper block is moved in the node direction, the resonant wavelength is slightly red-shifted; when it is moved perpendicular to the node direction, the resonant wavelength is largely blue-shifted. Based on these different dependences of the changes in wavelengths of the two modes, M1 and M2, we expect that the lateral position of the upper block can be estimated by observing the resonant wavelengths of M1 and M2 simultaneously.
The absorption cross-sections of the cavity structure are strongly enhanced at the resonant wavelength of the plasmonic cavity mode. In order to exploit the absorption cross-section spectra to obtain the resonant wavelength shifts of M1 and M2 for the lateral movement of the upper block, a monochromatic plane wave source with vertical polarization (
Ez) was injected into the proposed structure to excite the plasmonic modes M1 and M2 simultaneously. As shown in the inset of
Figure 1, the incidence angle of the plane wave was set to 40° and the absorption cross-section was calculated in the spectral range from 800 to 1250 nm, which includes the resonant wavelengths of M1 and M2.
Figure 4a shows the absorption cross-section spectrum for the original position of the upper block, indicated by (0-0) meaning (Δ
x = 0 nm and Δ
y = 0 nm). In this absorption cross-section spectrum, each peak position corresponds to the resonant wavelengths of M1 and M2 (882.2 and 1081.3 nm, respectively). Similarly,
Figure 4b–d show the absorption cross-section spectra for the lateral movement of the upper block: (30-0), (0-30), and (30-30), respectively. Black curves represent the absorption cross-section spectrum at the reference position, as shown in
Figure 4a, and red curves represent the absorption spectra with respect to the corresponding position of the upper block.
In the case of the movement of the block only along the
x-axis by 30 nm (denoted (30-0) in
Figure 4b), the corresponding absorption peaks in the red curve are 885.1 and 1064.7 nm, respectively, which correspond to M1 and M2. The shorter wavelength mode, M1, has a small red-shift of 2.9 nm and the longer wavelength mode, M2, has a large blue shift of 16.6 nm. In contrast, for the movement only along the
y-axis by 30 nm—(0-30) in
Figure 4c—M1 has a large blue-shift of 17.1 nm and M2 has a small red-shift of 2.9 nm. In the case of movement of Δ
x = 30 nm and Δ
y = 30 nm, (30-30), the two modes have large blue-shifts of 16.9 and 16.3 nm. Compared to the initial wavelengths, the shifted absorption resonant wavelengths have similar tendencies with the resonant wavelength shifts characteristic of M1 and M2, as depicted in
Figure 3b,c. A large wavelength shift of M1 (M2) along the
y- (
x-)axis exploits the resonant wavelength of M1 (M2) as a ruler for the
y- (
x-)directional position of the upper block.
The intensity of the mode which has a large wavelength shift for a given movement of the block changes considerably. For example, in
Figure 4b, the peak intensity of M2 which shows large wavelength shift decreases from 0.012 to 0.08 in contrast to the no significant change of M1 intensity because light coupling into each mode changes depending on the block positions. In order to calculate the absorption spectra, the absorbed energy in the block is calculated by subtracting the flux of the scattered light from the integrated flux of the incident light at each wavelength.
When the angle of incidence changes, the peak positions and linewidths of M1 and M2 are maintained, while the peak intensities change. For example, in the case of the incidence angle of 0°, along the y-axis, only M1 can be observed and M2 disappears. In contrast, for the incidence angle of 90°, along the x-axis, only M2 is observable. Therefore, in changing the incidence angle, one can excite M1 or M2 separately in the experimental measurement, which allows measurement of the spectral position of each peak more accurately by removing the overlap of peaks, resulting in higher sensitivity. In the experiment, a thin oxide layer could be observed on the surfaces of the silver metal blocks. The layer in the air gap induced a slight red-shift of the wavelength due to an increase of the effective optical path length in the gap; however, the operation mechanism of the 3D plasmonic ruler was maintained, except for the shift.
In order to determine the lateral position of the upper block more precisely via the two absorption peaks M1 and M2, we need to calculate the absorption cross-section peak wavelengths of M1 and M2 for 2D movement of the upper block using 2D color mapping of the wavelengths, as shown in
Figure 5a,b. In these 2D maps, the absorption peak wavelengths are displayed as different colors for Δ
x and Δ
y, ranging from 0 to 30 nm. In the wavelength map of M1 (
Figure 5a), the color of the corresponding wavelength is almost maintained along the
x-axis, Δ
x. However, the color changes dramatically along the
y-axis,
y-directional movement, Δ
y, which is demonstrated in
Figure 3b and
Figure 4. In the case of M2 in
Figure 5b, the dependencies are opposite, as shown in
Figure 3c and
Figure 4. Based on two-color mapping of M1 and M2 for the lateral movement of the upper block, it is simple to determine the position of the block if measuring the absorption cross-section spectrum of the incident light with an incidence angle of 40° and with vertical polarization. The method is as follows:
For example, in the case of Δ
x = 0 nm and Δ
y = 0 nm, the absorption peak wavelengths of M1 and M2 are 882.2 and 1081.3 nm, respectively. The wavelengths are depicted as cyan lines on
Figure 5a,b, M1 and M2 mapping. A point of intersection of the two lines is represented as a cyan circle that is almost at Δ
x = 0 nm and Δ
y = 0 nm. Similarly, in the case of Δ
x = 30 nm and Δ
y = 0 nm, the absorption cross-section peak wavelengths of M1 and M2—885.1 and 1064.7 nm—are drawn as black lines in
Figure 5a,b. The intersection, represented by a black circle, is Δ
x = 29.3 nm and Δ
y = 0.7 nm, which is the estimated lateral displacement of the upper block. Therefore, the lateral position of the block is precisely determined except for a small discrepancy that results from the finite grid size and the discrete scan step in 2D mapping. Consequently, from a single spectral measurement of two absorption peaks, one can determine the relative lateral displacement of the upper block relative to the reference lower block.
3. Three-Dimensional Positioning of a Nanoblock
Until now, we assumed the air gap thickness to be fixed at 10 nm. However, due to the strong plasmonic coupling in close proximity to metal nanoparticles, the resonant wavelength of the proposed structure has strong gap dependence [
1,
2,
3,
4,
5,
15,
16]. For the same structure as is in
Figure 3, we calculated the resonant wavelengths of M1 and M2 by varying the vertical position of the upper block along the
z-direction, as shown in
Figure 6a. The resonant wavelengths of M1 (black) and M2 (red) are plotted as functions of the gap thickness in
Figure 6b, with a thickness range from 2 to 20 nm. The wavelengths increase as the gap decreases, and the wavelengths dramatically increase at small gap size. For example, the resonance of M1 red-shifts from 882.2 to 949.2 nm with a change of 66.8 nm; for a decrease of the gap from 10 to 8 nm, there is a change of 2 nm. In contrast, for the same 2 nm change of the gap from 4 to 2 nm, a large red-shift of 359.8 nm is observed from 1164 to 1523.8 nm. This nonlinear dependence of the wavelength on the gap is widely observed in the gap plasmon resonances [
1,
2,
3,
4,
5]. The wavelength of the two modes, M1 and M2, have similar dependences on gap thickness.
As shown in
Figure 5, the lateral displacement of the upper block can be measured precisely using a two-nanoblock structure. However, because of the degrees of freedom in the 3D positioning of a nanoblock relative to the reference block, a combination of lateral and vertical movement cannot be estimated simultaneously in the two-block scheme. Thus, we introduce an extended nanoblock structure that can measure the position of a block three-dimensionally using the same mechanism: the measurement of the absorption peaks of orthogonal plasmonic modes in two silver nanoblocks.
Figure 7a shows the extended nanoblock structure, which consists of three silver nanoblocks. In this structure, we added a small silver nanoblock to the larger block of the two-block structure with the same air gap thickness of 10 nm. Top and bottom small blocks have a size of 200 nm × 150 nm × 100 nm, and the center large block has a size of 260 nm × 210 nm × 100 nm. Contrary to the two-block structure with one air gap, the extended structure consists of two air gaps: between the top small and the center large block, and between the center large and the bottom small blocks. Here, the top and bottom blocks are reference blocks and are assumed to have fixed positions. A center large block is moved three-dimensionally with values of Δ
x (nm), Δ
y (nm), and Δ
z (nm), which is measured by the absorption peaks due to the plasmonic modes in the two air gaps.
First, we calculate the absorption spectrum when the center block is placed at the origin with no movement, Δ
x = 0 nm, Δ
y = 0 nm, and Δ
z = 0 nm, denoted by (0-0-0), as shown in
Figure 7b. Because the thicknesses of the two air gaps are the same and both gaps are 10 nm, the absorption cross-section spectrum of the structure shows the same absorption peak wavelengths as those of the two-block structure in
Figure 4a: 882.2 and 1081.3 nm, the wavelengths of M1 and M2. However, when the center block moves along the
z-axis to 2 nm from the origin, the upper and lower air gaps become 8 and 12 nm, respectively. Due to the strong gap thickness dependence of the resonant wavelengths in the gap plasmonic cavity, as shown in
Figure 6b, there are large resonant wavelength differences between the plasmonic modes in the upper and lower air gaps.
Figure 7c shows the absorption spectrum of the center block with Δ
x = 0 nm, Δ
y = 0 nm, and Δ
z = 2 nm (the 0-0-2 position). In this figure, there are four absorption peaks at 849.0, 949.2, 1032.3, and 1163.2 nm, which are marked P1, P2, P3, and P4, respectively, from shortest to longest wavelength. P1 (849.0 nm) and P3 (1032.2 nm) indicate the absorption peaks of M1 and M2 excited in the air gap with a thickness of 12 nm, whilst P2 (949.2 nm) and P4 (1163.2 nm) indicate those in the gap with a thickness of 8 nm. The resonant wavelengths can also be confirmed by the curves of the gap dependences of M1 and M2 in
Figure 6b.
Next, we investigate the absorption spectrum when the lateral displacement of the center block has a vertical movement of 2 nm.
Figure 7d shows the absorption spectrum (red) corresponding to the position of the center block, Δ
x = 30 nm, Δ
y = 0 nm, and Δ
z = 2 nm, denoted by (30-0-2). In comparison with the absorption spectrum of Δ
x = 0 nm, Δ
y = 0 nm, and Δ
z = 2 nm, (0-0-2) (black), P1 and P2 corresponding to the M1 peaks of the air gap of 12 and 8 nm, respectively, are slightly red-shifted, and P3 and P4 corresponding to the M2 peaks of the air gap of 12 and 8 nm are largely blue-shifted. Wavelength shifts in M1 and M2 in the
x-directional movement of 30 nm are identical to the wavelength dependences of M1 and M2 in the two-block structure, as shown in
Figure 4b. Therefore, we can precisely estimate the position of the center block three-dimensionally using the single absorption spectrum of the three-block structure. On the other hand, four peaks—M1 and M2 modes for upper and lower air gaps—exist independently without any hybridization of each mode. Because M1 and M2 in the same air gap are orthogonal, the two M1 or M2 modes excited in different air gaps cannot couple each other due to large the block height of 100 nm.
The method of estimating the position of the center block from the absorption spectrum of the three-block structure is analogous to that of the two-block structure. From the absorption spectra (e.g., those in
Figure 7c,d), the wavelengths of the four peaks are determined as P1, P2, P3, and P4. Similar to the 2D color map of the wavelengths of M1 and M2 in
Figure 5, the absorption peak wavelengths of M1 and M2 in an air gap can be calculated for all of the 3D space of Δ
x and Δ
y and for Δ
z (the air gap thickness). For 3D mapping of the absorption peak wavelengths, the wavelength surfaces for P1 and P2 are plotted in the M1 mapping. In addition, the wavelength surfaces for P3 and P4 are plotted in the M2 mapping.
For example, in the absorption spectrum for the position of the center block, Δ
x = 30 nm, Δ
y = 0 nm, and Δ
z = 2 nm (30-0-2), P1 (852.1 nm), P2 (953.3 nm), P3 (1017.5 nm), and P4 (1148.3 nm) are observed. In
Figure 8a, the wavelength surfaces of P1 (red) and P2 (yellow) are plotted on the M1 mapping, and the surfaces of P3 (cyan) and P4 (blue) are plotted on the M2 mapping. Once the surfaces of P1 and P3 are obtained for the same air gap, a blue line—the line of the intersection of P1 and P3 surfaces—can be displayed, which is on the plane of z = 12 nm. Additionally, a red line—the intersection of the surfaces P2 and P4—can be drawn in the plane of z = 8 nm. Based on the planes of the red and blue lines, one can conclude that the vertical position of the center block is moved 2 nm from the original position.
Next, these lines are projected onto the x–y plane and the point of intersection of the two lines indicates the lateral position of the center block.
Figure 8 represents an orthographic projection of the blue and red lines onto the x–y plane. In this figure, the intersection of the two lines is indicated by a black circle and is at Δ
x = 28.7 nm and Δ
y = 2.2 nm. Finally, the position of the center block is estimated to be (Δ
x = 28.7 nm, Δ
y = 2.2 nm, Δ
z = 2 nm). As a result, the position of the center block can be estimated three-dimensionally using the spectral positions of the four peaks in the absorption spectrum and the calculated 3D mappings of the resonant wavelengths of M1 and M2 in the proposed structure. As mentioned previously, the small discrepancy results from the finite grid size and the discrete scan step in 3D mappings.
The method to determine 3D position of the block is as follows:
One prepares two 3D wavelength maps of M1 and M2 in 3D spaces of Δ
x, Δ
y, and Δ
z by calculating the wavelengths of M1 and M2 for all positions of the center blocks. The wavelengths of M1 and M2 are 882.2 and 1081.3 nm at origin, Δ
x = 0 nm, Δ
y = 0 nm, and
z = 10 nm. One wavelength of M1 or M2 is plotted as an equi-wavelength surface in 3D space, such as the equi-wavelength lines of
Figure 5. For a certain wavelength of a mode, only one equi-wavelength surface exists in a M1 or M2 3D map. At each
z position, the cross-sectional maps in the
x–
y plane correspond similarly to the 2D color mapping of the wavelengths shown in
Figure 5.
By measuring the absorption spectrum, the spectral positions of four peaks that are two short wavelength peaks are obtained: M1 modes of upper and lower gaps and two long wavelength peaks, and M2 modes of upper and lower gaps. Based on the four wavelengths and the two M1 and M2 wavelength maps in step 1, two M1 equi-wavelength surfaces are plotted on the M1 3D surfaces, and two M2 surfaces are plotted on the M2 3D surfaces similar to the two wavelength lines in the 2D color maps of
Figure 5.
One overlaps four equi-wavelength surfaces on a 3D space and obtains two lines from two intersections of M1 and M2 equi-wavelength surfaces. One set of M1 and M2 surfaces can make only a single intersection line, as shown in the red and blue lines in
Figure 8a.
Determine Δx and Δy of the center block by finding the intersection of two lines in 2D space of Δx and Δy. Based on the determined values of Δx and Δy and the four equi-wavelength surfaces in 3D spaces obtained in step 3, the vertical position, z, of the center block can be determined.