An Accessible Integrated Nanoparticle in a Metallic Hole Structure for Efficient Plasmonic Applications

Addressing the severe deterioration of gap mode properties in spherical-shaped nanoparticles (NPs) becomes necessary due to their utilization in a wide range of multi-disciplinary applications. In this work, we report an integrated plasmonic nanostructure based on a spherical-shaped nanoparticle (NP) in a metallic hole as an alternative to a NP-only structure. With the help of three-dimensional (3D) electromagnetic simulations, we reveal that when a NP is positioned on the top of a metallic hole, it can exhibit superior gap-mode-based local-field intensity enhancement. The integrated nanostructure displayed a ~22-times increase in near-field enhancement characteristics, similar to cube- or disk-shaped nanostructure’s plasmonic properties. From an experimental perspective, the NP positioning on top of the metallic hole can be realized more easily, facilitating a simple fabrication meriting our design approach. In addition to the above advantages, a good geometrical tolerance (metallic hole-gap size error of ~20 nm) supported by gap mode characteristics enhances flexibility in fabrication. These combined advantages from an integrated plasmonic nanostructure can resolve spherical-shaped NP disadvantages as an individual nanostructure and enhance its utilization in multi-disciplinary applications.


Introduction
A methodology to create and design nanostructures or devices before fabrication will save time and cost in developing efficient applications in various fields [1][2][3][4][5]. One such field that merited from this modeling approach is optics, where attractive applications in plasmonics, photonics, non-classical light sources, quantum dots, semiconductors, and so on can be realized [4,[6][7][8][9][10]. Notably from those applications, plasmonic nanostructures were studied extensively, which reveals unique optical properties such as the negative index of refraction, complex sub-wavelength characteristics, extraordinary transmission, tunable resonances, photocatalytic activity, and color printing [10][11][12][13][14][15]. In particular, exploiting interactions between light and matter with the help of surface plasmon resonance (SPR) results in variety of properties and functions. This SPR can be categorized into two parts: surface plasmon polariton (SPP)-the propagation of electron oscillations along the planar interface; and localized SPR (LSPR)-the confinement of electron oscillations on a subwavelength structure. These two forms of SPR, once excited, can lead to remarkable enhancement of local field or hot-spots and allows manipulation of light below the diffraction limit. This remarkable property of SPR helped in yielding a diverse range of applications such as artificial magnetism, imaging, sensing, energy, optical switching, and photodetection [1,4,6,[16][17][18][19][20][21].
Numerical modeling approaches based on three-dimensional electromagnetic simulations created a positive impact in achieving better plasmonic devices by exploiting a variety of geometrical designs [1,4,6,10,[14][15][16]. Of many features, the following will be key goals of the modeling approach: to develop efficient plasmonic nanostructures; the availability of geometrical tolerance with minimal loss in plasmonic properties; and non-complex low-cost fabrication [1,8,9,15]. Metallic nanoparticles (NPs) are one such example where there is a possibility of realizing attractive plasmonic characteristics through a simplistic approach [4,10,[22][23][24][25][26]. In particular, spherically shaped NPs have been studied extensively (experiment and theory) owing to their simplicity in fabrication, and ease of use in versatile applications. Despite the advantages mentioned above, rapidly deteriorating plasmonic characteristics from spherical NPs were observed on the presence of minor nanoscale geometrical errors [27][28][29][30]. Due to these issues, alternative geometries with respect to sphere-like disk and cube were preferred due to their excellent optical nature. Even though a NP-on-mirror (NPOM) design can enhance spherical NPs' optical properties, gap plasmonic properties become absent even with gap sizes exceeding ≥6 nm [27]. It is necessary to address this issue, where a simple and effective nanostructure involving spherical NPs can be realized. With the aid of numerical modeling, it is possible to envision efficient spherical NP plasmonic nanostructures facilitating simple fabrication [10,15,[18][19][20].
In this work, we report an integrated nanostructure (NP in a metallic hole on a SiO 2 substrate) which can be realized practically. The primary purpose of this simulation study is to reveal its plasmonic characteristics and to discuss its optical properties and possible advantages by utilizing this integrated nanostructure. We believe our work can act as a design platform, exploiting meaningful insights and thereby coalescing this integrated structure for multi-disciplinary applications.

Materials and Methods
Optical modeling is carried out using the three-dimensional (3D) finite-difference timedomain (FDTD) method (Lumerical FDTD solutions, ANSYS Inc. Vancouver, Canada). The modeled structure consists of a metallic hole with a nanoparticle (NP) on a SiO 2 substrate (Figure 1a). To begin with, we simplified the model consisting of a single hole with a single NP. For both the hole and NP, we used gold as a material. The geometrical parameters are as follows: hole's width and thickness are given by "h" and "t," respectively; NP diameter is given by "d" (Figure 1b). For most of this study, we used the following conditions: d = 100 nm; h = 100 nm; and t = 100 nm. In the case of geometrical tolerance analysis, the condition is set to h = t. The 3D structure is surrounded by perfectly matched layer (PML) boundary conditions in XYZ directions. A broadband plane-wave source is used to excite the structure in a normal direction from the top (+Z) with an incident field of "E 0 " and wave vector "k". Meshing conditions: "5 nm mesh covering the total simulation area and a 0.3 nm mesh-override size surrounding nanostructure of our main concern". The following refractive indices are used: Johnson and Christy database for gold and Palik database for SiO 2 [31,32]. Air surrounds the nanostructure on top with a refractive index of n = 1. A box-shaped power monitor is placed close to the nanostructure to record the nearfield enhancement. Approximately, the fourth power of near field or |E/E0| 4 can be compared to the electromagnetic enhancement factor of SERS (for clarity in linear scale graphing, we used |E/E0|). By assuming the Raman probe molecules are distributed randomly and uniformly on the NP surface, an averaged electromagnetic enhancement factor can be calculated by averaging the volume integral of |E/E0| [1,4,6,8,10,27,29,30,[33][34][35][36][37][38][39]: As seen from Equation (1), E0 is the modulus of the incident electric field (|E0| = 1 V/m), E represents the generated local field (E = Ex, Ey, Ez), and V is the volume at a certain distance (here we considered 2 nm) above the metallic surface [1,29,40,41].
To identify and interpret complex plasmonic modes, we utilized three-dimensional surface charge density mappings (3DSCDM) using COMSOL Multiphysics software (Wave optics module). By considering skin effect, an integration of Gauss's law, outward normal vector, and local electric field, 3D surface charge density " " is calculated [1,29,40,41].
From Equation (2), radius of NP is given as "R", outward normal vector "n", skin depth "δ" and permittivity of vacuum "ε0". Figure 2 displays the fundamental issue concerning individual nanostructures either in the form of NP or a metallic hole on a glass or SiO2 substrate. We did not use a metallic mirror as a substrate, to clearly evaluate and distinguish the individual nanostructure's problems. Broadband near-field spectra reveals a EFmax of ~8 can be extracted from individual nanostructures (Figure 2a). Cross-sectional electric field profiles reveal no trace of gap mode characteristics from these individual nanostructures (Figure 2b,c) extracted at 550 nm and 535 nm wavelength positions for NP and hole structures, respectively. NP mode or hole modes are solely identified from these individual nanostructures on a SiO2 substrate (unless there is a possibility of an increased number of NPs or with extremely smaller metallic hole width "h" size). In case of NPs, even in the presence of a metallic mirror substrate (NPOM design), gap mode deteriorates rapidly and diminishes when A box-shaped power monitor is placed close to the nanostructure to record the nearfield enhancement. Approximately, the fourth power of near field or |E/E 0 | 4 can be compared to the electromagnetic enhancement factor of SERS (for clarity in linear scale graphing, we used |E/E 0 |). By assuming the Raman probe molecules are distributed randomly and uniformly on the NP surface, an averaged electromagnetic enhancement factor can be calculated by averaging the volume integral of |E/E 0 | [1,4,6,8,10,27,29,30,[33][34][35][36][37][38][39]:

Problems with Individual Nanostructure(s)
As seen from Equation (1), E 0 is the modulus of the incident electric field (|E 0 | = 1 V/m), E represents the generated local field (E = E x , E y , E z ), and V is the volume at a certain distance (here we considered 2 nm) above the metallic surface [1,29,40,41].
To identify and interpret complex plasmonic modes, we utilized three-dimensional surface charge density mappings (3DSCDM) using COMSOL Multiphysics software (Wave optics module). By considering skin effect, an integration of Gauss's law, outward normal vector, and local electric field, 3D surface charge density "ρ" is calculated [1,29,40,41].
From Equation (2), radius of NP is given as "R", outward normal vector "n", skin depth "δ" and permittivity of vacuum "ε 0 ". Figure 2 displays the fundamental issue concerning individual nanostructures either in the form of NP or a metallic hole on a glass or SiO 2 substrate. We did not use a metallic mirror as a substrate, to clearly evaluate and distinguish the individual nanostructure's problems. Broadband near-field spectra reveals a EF max of~8 can be extracted from individual nanostructures (Figure 2a). Cross-sectional electric field profiles reveal no trace of gap mode characteristics from these individual nanostructures (Figure 2b,c) extracted at 550 nm and 535 nm wavelength positions for NP and hole structures, respectively. NP mode or hole modes are solely identified from these individual nanostructures on a SiO 2 substrate (unless there is a possibility of an increased number of NPs or with extremely smaller metallic hole width "h" size). In case of NPs, even in the presence of a metallic mirror substrate (NPOM design), gap mode deteriorates rapidly and diminishes when spacer thickness becomes ≥6 nm [27]. Presenting a plasmonic nanostructure design with good geometrical tolerance facilitating flexible fabrication and enhanced gap mode properties will be necessary in developing efficient applications. spacer thickness becomes ≥6 nm [27]. Presenting a plasmonic nanostructure design with good geometrical tolerance facilitating flexible fabrication and enhanced gap mode properties will be necessary in developing efficient applications.

Integrated Nanostructure: Nanoparticle in a Metallic Hole
Here, we introduce an integrated plasmonic nanostructure (NP in a metallic hole) on a SiO2 substrate as an alternative geometrical design to enhance plasmonic properties. In the first part of the optical analysis, we evaluated how the position of the NP in the hole determines the plasmonic properties. Following geometrical conditions used in this section: h = t = d. We applied these geometrical conditions for the following requirements: position of NP on the top of the hole (suitable with a state of h = d); simplification of a model (h = t = d); and flexible understanding of plasmonic properties from our structure. The inset figure of 3a displays the three modeled NP positions: bottom, middle, and top. In practical realization, a NP positioned on the top will be easier to fabricate but, at the same time, there can be a minor probability of the NP's movement towards the bottom. By considering the possibilities of NP locations, we compared top, middle, and bottom conditions. Figure 3a displays the near field enhancement EFmax spectra from these nanostructures. We marked three dominant resonance peak positions from "i" to "iii" to interpret the origin of plasmonic modes. Superior EFmax properties are revealed when NP is positioned at the top of the metallic hole. The bottom position displays comparatively lower EFmax values from its resonance wavelength marked as "i". The near field EFmax properties deteriorated from 173, 125, and 87 when the NP position moved from top to the bottom. The wavelength position marked at "i" is chosen to distinguish primary plasmonic properties due to its dominant gap mode characteristics, which are revealed from the cross-sectional XZ electric field profiles (Figure 3b). Highly concentrated, enhanced near-field spots can be seen at metal-metal contact positions (with the better case being the NP positioned at top) in the case of "i," explaining its dominant gap mode origin when compared to "ii" and "iii". In comparison to "i", weakly enhanced near-field spots can be seen at metal-metal contact positions at "ii", thereby assigning it to "supporting gap mode". NP mode is observed at "iii" as there is no trace of enhanced local field at metalmetal contact [27,29].

Integrated Nanostructure: Nanoparticle in a Metallic Hole
Here, we introduce an integrated plasmonic nanostructure (NP in a metallic hole) on a SiO 2 substrate as an alternative geometrical design to enhance plasmonic properties. In the first part of the optical analysis, we evaluated how the position of the NP in the hole determines the plasmonic properties. Following geometrical conditions used in this section: h = t = d. We applied these geometrical conditions for the following requirements: position of NP on the top of the hole (suitable with a state of h = d); simplification of a model (h = t = d); and flexible understanding of plasmonic properties from our structure. The inset Figure 3a displays the three modeled NP positions: bottom, middle, and top. In practical realization, a NP positioned on the top will be easier to fabricate but, at the same time, there can be a minor probability of the NP's movement towards the bottom. By considering the possibilities of NP locations, we compared top, middle, and bottom conditions. Figure 3a displays the near field enhancement EF max spectra from these nanostructures. We marked three dominant resonance peak positions from "i" to "iii" to interpret the origin of plasmonic modes. Superior EF max properties are revealed when NP is positioned at the top of the metallic hole. The bottom position displays comparatively lower EF max values from its resonance wavelength marked as "i". The near field EF max properties deteriorated from 173, 125, and 87 when the NP position moved from top to the bottom. The wavelength position marked at "i" is chosen to distinguish primary plasmonic properties due to its dominant gap mode characteristics, which are revealed from the cross-sectional XZ electric field profiles (Figure 3b). Highly concentrated, enhanced near-field spots can be seen at metal-metal contact positions (with the better case being the NP positioned at top) in the case of "i," explaining its dominant gap mode origin when compared to "ii" and "iii". In comparison to "i", weakly enhanced near-field spots can be seen at metal-metal contact positions at "ii", thereby assigning it to "supporting gap mode". NP mode is observed at "iii" as there is no trace of enhanced local field at metal-metal contact [27,29]. Materials 2022, 15, x FOR PEER REVIEW 5 of 12 It will be difficult to judge plasmonic properties solely based on cross-sectional electric field profiles. Utilizing 3DSCDM calculations will be helpful in understanding the optical properties in detail. 3DSCDM profiles for integrated nanostructures taken at resonance wavelength positions marked from "i" to "iii" are displayed in Figure 4. To explain plasmonic modes in an easier way, we used the following symbols: a solid square symbol for presence of the dipole mode and a solid diamond symbol for the quadrupole mode [1]. Figure 4a-c shows 3DSCDM profiles for integrated nanostructures with the NP positioned from bottom, middle to top, respectively (green dotted inset figure displays surface charge mappings in NP for the respective integrated nanostructures). For "i" wavelength position, dipole mode is clearly visible for all the three integrated nanostructures. Differences in plasmonic modes appear for "ii" and "iii" positions. For both bottom and middle NP positions, the quadrupole mode is noted as "ii" and "iii". However, for the top NP position, it differed: the dipole mode appears at "ii" and quadrupole mode at "iii". In addition, for the top NP position, it is significant to note that whispering gallery modelike surface charge interactions (Figure 4c,d, observed at resonance wavelengths "i" and "ii"; absent for "iii") could be seen at the NP-hole interface, a probable reason for higher EFmax contribution as compared to middle and bottom NP positions. No such whispering gallery mode-like surface charge interactions could be seen with the middle and bottom NP positions ( Figure 4e); thus, we can conclude that when a strong coupling between a brighter (dipole) gap mode and whispering gallery mode-like surface charge interactions (at NP-hole interface) happens, it can generate a highly localized field from the integrated nanostructure. It will be difficult to judge plasmonic properties solely based on cross-sectional electric field profiles. Utilizing 3DSCDM calculations will be helpful in understanding the optical properties in detail. 3DSCDM profiles for integrated nanostructures taken at resonance wavelength positions marked from "i" to "iii" are displayed in Figure 4. To explain plasmonic modes in an easier way, we used the following symbols: a solid square symbol for presence of the dipole mode and a solid diamond symbol for the quadrupole mode [1]. Figure 4a-c shows 3DSCDM profiles for integrated nanostructures with the NP positioned from bottom, middle to top, respectively (green dotted inset figure displays surface charge mappings in NP for the respective integrated nanostructures). For "i" wavelength position, dipole mode is clearly visible for all the three integrated nanostructures. Differences in plasmonic modes appear for "ii" and "iii" positions. For both bottom and middle NP positions, the quadrupole mode is noted as "ii" and "iii". However, for the top NP position, it differed: the dipole mode appears at "ii" and quadrupole mode at "iii". In addition, for the top NP position, it is significant to note that whispering gallery mode-like surface charge interactions (Figure 4c,d, observed at resonance wavelengths "i" and "ii"; absent for "iii") could be seen at the NP-hole interface, a probable reason for higher EF max contribution as compared to middle and bottom NP positions. No such whispering gallery mode-like surface charge interactions could be seen with the middle and bottom NP positions ( Figure 4e); thus, we can conclude that when a strong coupling between a brighter (dipole) gap mode and whispering gallery mode-like surface charge interactions (at NP-hole interface) happens, it can generate a highly localized field from the integrated nanostructure.
Surface charge interactions at the center of NP (XY direction) in the integrated nanostructures happen because of NP-hole contact (orange dotted line from the inset Figure 4a  The assignment of gap modes for marked positions "i" and "ii" from integrated nanostructures in Figure 3a can be justified by comparing it to individual nanostructures ( Figure 2) with the following characteristics: highly enhanced EFmax values; brighter localfield spots originating at metal-metal contacts; and red-shifted wavelength positions. Critically, at the same time, an integrated nanostructure with the top NP position displayed better EFmax properties based upon two brighter gap modes (dipole mode) in combination with whispering gallery mode-like charge properties at the NP-hole interface. The assignment of gap modes for marked positions "i" and "ii" from integrated nanostructures in Figure 3a can be justified by comparing it to individual nanostructures ( Figure 2) with the following characteristics: highly enhanced EF max values; brighter local-field spots originating at metal-metal contacts; and red-shifted wavelength positions. Critically, at the same time, an integrated nanostructure with the top NP position displayed better EF max properties based upon two brighter gap modes (dipole mode) in combination with whispering gallery mode-like charge properties at the NP-hole interface. Significantly, integrated nanostructures exhibited a maximum of a~22-times increase in near-field enhancement EF max when compared with individual nanostructure(s), which is critically advantageous to the nano-structural design.
To further strengthen the importance of the integrated nanostructure, we calculated the Purcell factor in understanding the local intensity enhancement. We had extracted the Purcell factor from integrated nanostructures and compared it with individual NP structure ( Figure 5). For this purpose, we had placed the emitter close to the NP's EF max area as schematically shown in Figure 5a [42,43]. The Purcell enhancement is~10 times better with NP positioned at the top (integrated nanostructure) as compared to the individual nanostructure. Combination of higher Purcell factor and local field enhancement properties from the integrated nanostructure (with NP on the top) will benefit surface-enhanced applications, cavity nanostructures, sensing, and so on [44].
To further strengthen the importance of the integrated nanostructure, we calculated the Purcell factor in understanding the local intensity enhancement. We had extracted the Purcell factor from integrated nanostructures and compared it with individual NP structure ( Figure 5). For this purpose, we had placed the emitter close to the NP's EFmax area as schematically shown in Figure 5a [42,43]. The Purcell enhancement is ~10 times better with NP positioned at the top (integrated nanostructure) as compared to the individual nanostructure. Combination of higher Purcell factor and local field enhancement properties from the integrated nanostructure (with NP on the top) will benefit surface-enhanced applications, cavity nanostructures, sensing, and so on [44].

Geometrical Tolerance of Metallic Hole
We also considered a possible situation where hole width could be larger than NP diameter (h > d) based upon fabrication errors. To analyze the integrated nanostructure's geometrical tolerance, a geometrical parameter "hole gap" (Figure 6) was introduced and its influence towards gap-mode based near-field enhancement was evaluated. For this study, a NP positioned at the bottom will be a practical possibility. The hole gap is varied from 5 nm to 30 nm in 5 nm steps, and its influence over broadband near field spectra is shown in Figure 5a. Extracted EFmax values from its respective resonance wavelength positions as a function of "hole gap" sizes are plotted in Figure 6b. A decrease in near-field strength is seen when the "hole gap" size increases. The resonance wavelength position of gap mode blue-shifted from 862 nm to 691 nm (171 nm span) with a decrease in EFmax from 72 to 17 when the hole-gap size varied from 5 nm to 30 nm. The corresponding crosssectional XZ electric field profiles taken at their resonance wavelength positions were displayed in Figure 6c-h (5 nm to 30 nm in order). Gap mode properties are observed till the hole-gap size of ~20 nm. Once the plasmonic gap mode became negligibly small or severely deteriorated, EFmax resonance wavelength shift and intensity became ~ constant (hole-gap size > 20 nm). To distinguish these characteristics clearly, we extracted 3DSCDM profiles for hole-gap sizes of 5 nm (Figure 6i) and 25 nm (Figure 6i). Even though the surface charge interactions look identical on the NP surface in both cases, a difference is clearly seen on the substrate surface. A dipole mode-like interaction is observed when the

Geometrical Tolerance of Metallic Hole
We also considered a possible situation where hole width could be larger than NP diameter (h > d) based upon fabrication errors. To analyze the integrated nanostructure's geometrical tolerance, a geometrical parameter "hole gap" (Figure 6) was introduced and its influence towards gap-mode based near-field enhancement was evaluated. For this study, a NP positioned at the bottom will be a practical possibility. The hole gap is varied from 5 nm to 30 nm in 5 nm steps, and its influence over broadband near field spectra is shown in Figure 5a. Extracted EF max values from its respective resonance wavelength positions as a function of "hole gap" sizes are plotted in Figure 6b. A decrease in near-field strength is seen when the "hole gap" size increases. The resonance wavelength position of gap mode blue-shifted from 862 nm to 691 nm (171 nm span) with a decrease in EF max from 72 to 17 when the hole-gap size varied from 5 nm to 30 nm. The corresponding cross-sectional XZ electric field profiles taken at their resonance wavelength positions were displayed in Figure 6c-h (5 nm to 30 nm in order). Gap mode properties are observed till the hole-gap size of~20 nm. Once the plasmonic gap mode became negligibly small or severely deteriorated, EF max resonance wavelength shift and intensity became~constant (hole-gap size > 20 nm). To distinguish these characteristics clearly, we extracted 3DSCDM profiles for hole-gap sizes of 5 nm (Figure 6i) and 25 nm (Figure 6i). Even though the surface charge interactions look identical on the NP surface in both cases, a difference is clearly seen on the substrate surface. A dipole mode-like interaction is observed when the hole gap size is 5 nm (Figure 6i). For a hole-gap size of 25 nm, a quadrupole mode appeared (Figure 6j). This explains the why at a 5 nm hole-gap size we can extract better EF max , as compared to a 25 nm hole-gap size. A geometrical tolerance of~20 nm is significantly vital with recent advances in fabrication methods. hole gap size is 5 nm (Figure 6i). For a hole-gap size of 25 nm, a quadrupole mode appeared (Figure 6j). This explains the why at a 5 nm hole-gap size we can extract better EFmax, as compared to a 25 nm hole-gap size. A geometrical tolerance of ~20 nm is significantly vital with recent advances in fabrication methods.

Discussion
In general, NP positioned at the top (the best structure from this work) can be easily assembled upon the success of the hole nanostructure's identical width (w = d) fabrication.

Discussion
In general, NP positioned at the top (the best structure from this work) can be easily assembled upon the success of the hole nanostructure's identical width (w = d) fabrication. Lithography (nanoparticle lithography or electron-beam lithography) or self-assemblybased methods can facilitate the metallic hole fabrication [45][46][47][48]. After the successful formation of holes, plasmonic NPs can be spin coated. Depending upon hole and NP geometry, top or bottom NP positions can be realized. Enhanced Purcell factor, better EF max , and facile fabrication can be realized from the integrated nanostructure. Even though free movements of NP inside the hole can be useful in interesting sensing and trapping studies, NP positioned on the top of the hole will yield better near-field enhancement properties supported by brighter gap modes [49][50][51]. The benefits of our integrated nanostructure can be applied to various applications including catalytic devices. A strong coupling between the plasmon and optical modes will be attractive in facilitating catalytic reactions. As seen from our best-integrated nanostructure (top NP position), superior EF max are extracted as enhanced plasmonic (brighter) mode interactions between the NP and hole occured. Furthermore, enhancement in quantum efficiency will be possible on the basis of strong coupling of plasmonic mode(s) over that of decoupled conditions [52][53][54][55][56][57]. Precisely designed plasmonic structures are required to effectively promote catalytic reactions. The plasmonic structure concentrates the absorbed energy on far-field or near-field enhancement depending on the morphological characteristics of the structure. When we achieve dominant near-field enhancement properties, the energy absorbed by the plasmonic structure contributes to the catalytic reaction of its surroundings rather than emitting it to the outside. The root cause of this difference is based on the plasmonic mode. However, calculating and interpreting three-dimensional plasmonic modes is not yet widespread. Precise design examples (as reported in this work) help catalytic researchers to make the right choice(s) for using plasmonic structures. In addition to the above benefits, a good geometrical tolerance with a presence of gap-mode characteristics facilitates flexible fabrication. Considering these merits, it is possible to realize an efficient plasmonic nanostructure in an integrated form similar to that of a disk-or cube-based NPOM design [8,30]. Our optical design strategy based on an integrated nanostructure can open various attractive applications supported by a simple fabrication.

Conclusions
In summary, we had reported a three-dimensional finite-domain time-difference simulation study on a NP-hole on a SiO 2 substrate (integrated nanostructure) and discussed its plasmonic properties. NP positioned at the top of the hole showed better plasmonic characteristics with an enhanced gap mode capability (a maximum of a~22-times increase in near field strength) and Purcell factor resolving the individual nanostructure's disadvantage. Effective geometrical tolerance with the availability of gap mode properties till a hole-gap size of~20 nm can be realized, facilitating flexible fabrication. In addition to all the above advantages, the probability of positioning an NP on top of a hole can be more straightforward, resulting in a simple fabrication. An efficient plasmonic nanostructure or devices based on this integrated (NP-hole) nanostructure can foster a variety of applications in the field of photocatalysis, energy, sensors, surface-enhanced spectroscopy, and ultrafast optoelectronic applications.

Conflicts of Interest:
The authors declare no conflict of interest.