Substrate and Fano Resonance Effects on the Reversal of Optical Binding Force between Plasmonic Cube Dimers.

The behavior of Fano resonance and the reversal of near field optical binding force of dimers over different substrates have not been studied so far. Notably, for particle clustering and aggregation, controlling the near filed binding force can be a key factor. In this work, we observe that if the closely located plasmonic cube homodimers over glass or high permittivity dielectric substrate are illuminated with plane wave, no reversal of lateral optical binding force occurs. But if we apply the same set-up over a plasmonic substrate, stable Fano resonance occurs along with the reversal of near field lateral binding force. It is observed that during such Fano resonance, stronger coupling occurs between the dimers and plasmonic substrate along with the strong enhancement of the substrate current. Such binding force reversals of plasmonic cube dimers have been explained based on the observed unusual behavior of optical Lorentz force during the induced stronger Fano resonance and the dipole-dipole resonance. Although previously reported reversals of near field optical binding forces were highly sensitive to particle size/shape (i.e. for heterodimers) and inter-particle distance, our configuration provides much relaxation of those parameters and hence could be verified experimentally with simpler experimental set-ups.

Introduction

Since its discovery, Fano resonance has been a typical feature of interacting quantum systems. However, recently the Fano resonance has been found in plasmonic nanoparticles, photonic crystals, and electromagnetic metamaterials. Fano resonance and the Fano line shape[1] cannot be described by the conventional and well-known Lorentz formula[2]. In plasmonic nanostructures, Fano resonance can happen due to the resonant destructive interference between a super-radiant (bright) mode and subradiant (dark) modes[3]. The promising applications of plasmonic hybridization and Fano resonances[4, 5] have been investigated in improved sensitivity of the resonance[6], bio sensing[7], surface-enhanced Raman scattering[8], photonic propagation and wave guiding[9], plasmon-induced transparency[10] and many others[11]. In contrast, much less attention is dedicated on near field optical force due to Fano resonance; especially for plasmonic dimers[4, 5] as discussed only in[12, 13]. Specially, plasmonic particle over substrate show remarkable properties: Fano resonances[14], broadband tunability in plasmonic resonance[15], modification of energies of the plasmon modes[16], changing the localized density of states[17], radiative enhancement[18], frequency shift of an electric dipole resonance[19] and so on. But the behavior of near field optical force for plasmonic dimers over different types of substrate and the effect of substrate on the reversal of near filed optical binding force have not been studied in literature. In fact, only few works have been reported on the reversal of near field optical binding force due to Fano resonance, i.e. for nanobar structures[12], for disk along with a ring structure[13] and those are highly dependent on inter-particle distance and particle size.

Notably, for particle clustering, particle aggregation and for improved sensors, controlling the near filed binding force is considered as a key factor. According to ref. 20: “The inter particle separations are usually comparable to the wavelength of the illuminating laser beam and, therefore, mainly long-range interactions mediated by the far-field scattered field are considered while the near-field coupling is usually omitted in the related studies.” But with the recent technology, the inter particle distance between the nano particles can be precisely controlled[21]. More importantly, the size and shape of the gap between the nanoparticle and film can also be controlled to sub nanometer precision bottom-up fabrication approaches[21]. The film-coupled nanoparticle geometry has recently been applied to enhancing optical fields, accessing the quantum regime of plasmonics[22], and the design of surfaces with controlled reflectance[23]. Still no investigation has been carried out regarding the behavior of both Fano resonance and the reversal of near field optical binding force with respect to the inter particle distance of two dimers over different substrates.

In general, reversal of near field lateral optical binding force for the same polarization of light is quite uncommon with the homodimers placed without substrates[24-27]. Even if the homodimers are placed over substrates, reversal of near field lateral optical binding force has not been observed in refs 28–30 for spherical shaped and rod shaped plasmonic homodimers. On the other hand, it is well known that Fano resonance is ultra-sensitive and induced reversal of optical binding force dies out very quickly if the inter particle gap of the heterodimers[12] increase even a little bit. Moreover, this Fano resonance induced lateral binding force reversal is highly dependent on particle size[12], which has not been achieved by homodimers so far. In this work, it is observed that if two plasmonic cubic homodimers (also applicable for heterodimers) are placed closely without any substrate or above glass or high refractive index substrate, no reversal of optical binding force occurs. But if we apply the same set-up over a plasmonic substrate, stable and stronger Fano resonance occurs along with the reversal of lateral binding force. Importantly, the inter-particle distance and the sizes of the cube homo-dimers are much less sensitive than the previously reported case in ref. 12.

In this article, we have investigated a fully alternative way in comparison with refs 12, 28–30 to achieve strong and stable reversal of optical lateral binding force based on cubic homodimer and substrate system, where stable and much stronger Fano resonance can be achieved due to the substrate coupling with larger area of the lower side of cube dimers. Though nanoparticles of sphere are easy to synthesize in experiment, the substrate-induced hybridization of the plasmon modes can be much larger[14] for a planar metallic nanoparticle (i.e. a cube) deposited on a plasmonic substrate[21, 23] than for a spherical particle[16]. This happens because the plasmon-induced steady state currents and surface charges will be located closer to the strong surface currents and plasmonic screening charges induced on the surface of the substrate respectively over a large area.

We have demonstrated that from the interplay between localized surface plasmon and propagating surface plasmon polariton along with the strong coupling between the two particles with the plasmonic substrate, the reversal of the optical binding force occurs during the strong Fano resonance. According to several full wave numerical results, our configuration provides much relaxation on some critical parameters (i.e. inter particle distance and particle size as stated previously) and hence can probably be verified experimentally with much simpler experimental set-ups.

Results and Discussions

We specify that throughout this paper we refer to ‘exterior’ or ‘outside’ forces as those evaluated outside the volume of the macroscopic particles, while ‘interior’ or ‘inside’ refer to those quantities inside this object volume. In order to consider the realistic effects, we have done all the numerical calculations using full wave simulations[31] in three dimensional (3D) structures.

The proposed simple set-up is illustrated in Fig. 1(a). The Silver cubes are placed near to each other above a silver substrate (whole set-up is embedded in water). Spacer (refractive index 1.4) of height 5 nm is always considered above all the substrates in this article. The real and imaginary part of the permittivity of silver is taken from the standard Palik data[31, 32]. Inter particle distance is d. The source is a simple x-polarized plane wave propagating towards −z direction and E 0 has been set to 1 V/m (so, the incident intensity is 0.001176 W/m2) in the simulation model. Our set-up is a symmetry broken system which later plays a vital role for the force reversal. Details on the calculation process of optical force (based on both Minkowski stress tensor and modified Lorentz force) has been discussed in the ‘METHODS’ section (at the end of this article). The lateral binding force is defined as, . Positive value of means attractive force and negative value means repulsive force. Here subscript (x), (1) and (2) represent: +x direction, left cube and right cube respectively.

Figure 1

(a) Two silver cubes (L = 120 nm) are placed over silver substrate and their inter-particle gap, d, is 100 nm (spacer height 5 nm from the bare substrate and the cubes are placed 5 nm away from the spacer). (b) Comparison of the extinction coefficients of two silver cubes [set-up of Fig. 1(a)] and a single silver cube placed over silver substrate. (c) The binding force of set-up (a). E 0 has been set to 1 V/m (so, the incident intensity is 0.001176 W/m2) in the simulation model. (d–f) Charge distributions at wavelengths: 966 (DD), 754 (Fano) and 627 (QQ) nm respectively. D and Q represent dipole and quadrupole respectively. (g–i) Electric field distribution for those same wavelengths respectively.

Plasmonic cubes over plasmonic substrate: Resonance modes and the reversal of binding force

In Fig. 1(a) we have considered two silver cubes of length 120 nm each (100 nm apart). These cubes are placed in water above a silver substrate. A spacer of 5 nm height is considered between the cubes and the substrate and the cubes are placed 5 nm above the spacer. This set-up is shined by a plane wave described in previous section. The extinction spectra in Fig. 1(b) reveals that Fano dip takes place at around 754 nm and the strength of extinction coefficient increases two times than a single cube placed over silver substrate. The quadropolar resonance (at wavelength 627 nm), dipolar resonance (at 966 nm) and Fano dip (at around 754 nm) do not experience shifting due to the addition of the identical cube.

It is demonstrated in Fig. 1(c) that reversal of lateral optical binding force takes place near the Fano dip at 754 nm wavelength. It should also be noted that if we decrease the inter particle distance of the cubes, usually the magnitude of the lateral binding force () increases. However, this does not mean that the magnitude of the repulsive force (reversed binding force) increases. For example- in Supplementary S1 we have shown the nature of binding force of two Ag cubes of 120 nm for the inter-particle gap of 10 nm, 25 nm, 50 nm, 175 nm and 200 nm. Notably, for extremely close inter-particle gap (i.e. 10 nm distance), the reversal of optical binding force vanishes but the magnitude of reversed binding force is quite strong for 25 nm inter-particle gap (and also for 50 nm gap) as shown in Supplementary S1 and for 100 nm gap in main article. Interestingly, with the increase of the inter-particle gap, the binding force magnitude decreases slowly [cf. Supplementary S1 for 120 nm cubes with inter-particle distances of 175 and 200 nm and in Supplementary S2 for 180 nm cubes with inter-particle distance of 150 nm] but still this overall set-up provides a very relaxed mechanism to verify the reversal of optical binding force in comparison with the previous case reported in ref. 12.

The distinct charge distributions of the plasmonic cubes have been shown in Fig. 1(d)–(f) for different resonance modes [i.e. dipole-dipole (DD), Fano and quadrupole-quadrupole (QQ) resonance respectively]. In addition, in Fig. 1(g)–(i) the electric field magnitudes of different modes [i.e. DD, Fano and QQ resonance respectively] have been demonstrated. The field distributions (field patterns) look similar (but certainly not identical) to ref. 14 where the substrate is considered as high refractive index dielectric.

At the DD resonance mode (i.e. 966 nm wavelength), field lines form closed loops around the separation between the particle and the metallic film[33]. Maximum value of the repulsive binding force occurs at this wavelength as shown in Fig. 1(c). This resonance mode corresponds to the localized surface plasmon (LSP) mode associated to the silver particle. According to another point of view[21], this resonance mode can be explained as a magnetic cavity mode that roughly satisfies the half wavelength criteria. In particular, the metal (silver)–insulator (spacer)–metal (silver) region supports a transmission line type of mode modified by the plasmonic dispersion of the metal. A large effective impedance mismatch occurs at the edges of the cube and most of the energy is reflected back under the cube [cf. the electric field distribution in Fig. 1(g)].

On the other hand, a resonance peak of extinction spectra occurs at around 627 nm as shown in Fig. 1(c), which corresponds to QQ resonance. Local maxima of the attractive binding force occur during this resonance. It corresponds to the surface plasmon polariton (SPP) mode propagating on the silver/air interface of the metallic film [cf. ref. 33], which is excited near this wavelength according to ref. 33. The corresponding enhancement of the particle excitation field increases the intensity of field scattered at infinity [cf. the electric field distribution in Fig. 1(i)]. From another point of view[21], this higher order mode is an example of waveguide modes and occurs at shorter wavelengths such that the cavity is longer than a half wavelength of the mode.

Plasmonic cubes over plasmonic substrate: Further discussion on Fano resonance mode

In Supplementary S2 it is shown that the reversal of binding force still occurs at the Fano dip position (around 1000 nm) even if we change the length of the cube from 120 to 180 nm. Such reversal of binding force is quite generic for other sized plasmonic cubes over plasmonic substrates. The only countable issue to achieve such force reversal is to achieve Fano dip at first with the single cube. However, the force reversal wavelength and the Fano resonance both red shift with the increase of the cube dimension.

In our optical set-up, the reason of stronger Fano resonance for the homodimers can be explained based on the destructive interference of stronger quadropolar and stronger dipolar mode of the whole system: plasmonic homodimers in presence of the plasmonic substrate. Though Fano resonance can be induced based on a single plasmonic cube and the plasmonic substrate as previously discussed considering a single cube over high refractive index substrate in ref. 14, the presence of the second dimer (homodimer set-up) makes both the DD and QQ resonances much stronger than the single cube and substrate system as shown in Fig. 1(b). Eventually, this leads to the stable and stronger Fano resonance, which is very important for the magnitude of the lateral binding force (especially when its sign reverses). On the other hand, if a second object is placed, which has bigger dimension than the first one (heterodimer set-up); the dipolar and quadrupole resonance of the first cube and the second cube usually do not match at the same wavelengths. Eventually this may affect the strength of the Fano resonance and also the magnitude of the lateral binding force (especially when its sign reverses [cf. Supplementary S3]).

The strength of electric field magnitude increases a lot (in comparison with all other substrates) due to the existence of plasmonic substrate, which is a necessary condition for stronger coupling and resonance that leads to binding force reversal. However, the mixed nature of electric field at the Fano dip clearly supports the mixing of dipolar and quadrupolar resonance[14, 15].

But the question is whether Fano resonance is the only factor to achieve such force reversal. The answer is: enhanced Fano resonance induced by strong plasmonic coupling is the key factor of such force reversal but it is not the only reason. For example- in the next section we shall demonstrate that Fano resonance also occurs for Silicon substrate with 120 nm lengthen cube dimers but no reversal of binding force takes place for such case. In addition, more details on the physical mechanism of such optical force reversal will be discussed in the next sections.

Plasmonic cubes over different substrates: effect on binding force

So far we have discussed regarding single or double plasmonic cubes over plasmonic substrates. Now we shall compare the behavior of optical binding force considering two plasmonic cubes without substrate and over the glass or Silicon (high refractive index) or silver (plasmonic substrate) embedded in water. In Fig. 2(a)–(c), no reversal of optical binding force is observed in Fig. 2(c). This same conclusion is true for the 2nd case [cf. Fig. 2(d)–(f)] where no reversal of optical binding force is observed. Albeit, in Fig. 2(g)–(i), the substrate is Silicon (modelled as Palik data[31, 32]), Fano dip spectral region is observed in Fig. 2(h) due to Fano resonance. But no reversal of optical binding force is observed in Fig. 2(i). This is an important observation. This observation suggests that Fano resonance is not the only criteria to achieve the reversal of optical binding force [which occurs for silver slab in Fig. 2(l)] rather some other factors (will be discussed next) significantly influence the origin of the reversal of optical binding force. Now we shall discuss the behavior of binding force distribution for these three resonances: (i) DD (ii) QQ and (iii) Fano resonance, when the substrate is silver.

Figure 2

Two silver cubes (L = 120 nm) are placed over different substrates and their inter-particle gap, d, is 100 nm. The spacer height is always 5 nm. Cubes are placed 5 nm (h = 5 nm) away from the spacer. x-polarized plane wave is propagating towards –z-direction. (a–c) No substrate is placed; the extinction coefficient and binding force for that configuration. (d–f) Glass substrate (refractive index 1.5) is placed; the extinction coefficient and binding force for that configuration. (g–i) Silicon substrate is placed; the extinction coefficient and binding force for that configuration. (j–l) Silver substrate is placed; the extinction coefficient and binding force for that configuration.

For the DD resonance mode, the resonance process with the plasmonic substrate is quite different in comparison with the other substrates. During the DD resonance mode, the accumulated charges are strongly coupled with the plasmonic substrate and the total power is mostly concentrated between the cube and the plasmonic substrate. The charge distribution essentially results in a current loop between the nano cube and metal film inducing a strongly enhanced magnetic field. This resonance is then mainly of magnetic nature[21, 23, 33]. But the most interesting part is that: though the resonance is of magnetic nature, the difference of the scattering part[34] (which originates from magnetic induction:  ×  in; where is the polarization induced current[34]) or bulk part of the total Lorentz force [cf. Eq. (6) in the ‘METHODS’ section] vanishes during the peak of LSP or DD resonance mode [cf. Fig. 3(a),(c) for both 120 and 180 nm lengthened cubes]:

Figure 3

For cube homodimers over silver substrate, plot of and [cf. Eq. (1) and (2) in main text]. Spacer height is always 5 nm and height of the cubes from the spacer is fixed 5 nm from the spacer. For first column: Length of each cube is 120 nm and they are placed 100 nm apart. For second column: Length of each cube is 180 nm and they are placed 150 nm apart.

Here and j (=1 or 2) represents the left (denoted as ‘1’) or right cube (denoted as ‘2’) respectively in the dimer set. ε is permittivity and μ is permeability of the scatterer. On the other hand, is permittivity and μ is permeability of the background (always water in this main article). So, the total binding force at LSP or DD resonance peak is fully due to the surface charges and the difference of surface force [cf. Eq. (5) in the ‘METHODS’ section] can be expressed as [cf. Figure 3 (b),(d) for both 120 and 180 nm lengthened cubes]:Here

The unit vector is the outward pointing normal to the surface. It should be noted that .

In Fig. 4(i) it can be seen that the flipped [and also Del FSurf (x); cf. Fig. 3(b) and (d)] reaches its maximum negative value near the DD resonance for plasmonic silver substrate [also cf. Fig. 2(l)]. From the first row of Fig. 4, it can also be seen that such reversal of does not happen for other substrates. This can be explained based on the electric field coupling between the cubes and also between cubes with plasmonic substrate. Electric field enhancement becomes much stronger for the case of plasmonic substrate in comparison with other substrates especially at this particular resonance. The surface force in Eq. (5) in the ‘METHODS’ section depends on the electric fields of the cube boundary/interface, which becomes much stronger during this DD resonance.

Figure 4

For the cases of different substrates, Lorentz force components and the steady state current (J y) of two silver cubes along with the substrate. First row represents the Lorentz force components: surface force, bulk force [cf. Eqs (5) and (6) in ‘METHODS’ section] and total force [cf. Eqs (3) and (4) in ‘METHODS’ section] only on cube-1 (left cube, denoted as P1). ‘°’ represents the chosen wavelengths for which the steady state current distributions are plotted later. Second, third and fourth row represent the steady state current (J y) from front view [in xz plane; setting the window very close to the cube surfaces from front view] for three different wavelengths (marked as ‘°’) respectively: (a–d) For glass substrate where wavelengths are chosen: 754, 622 and 550 nm respectively. (e–h) For Silicon substrate where resonance wavelengths are: 816 (DD), 679 (Fano) and 578 nm (QQ) respectively. (i–l) For Silver substrate where resonance wavelengths are: 966 (DD), 754 (Fano) and 627 nm (QQ) respectively.

Though the maximum value of repulsive binding force for DD or LSP mode is just due to surface charges, the local maximum value of attractive binding force for QQ or SPP resonance at wavelength 627 nm [cf. Fig. 2(l)] is the combined effect of static (opposite surface charges as shown in Fig. 1(f)) and dynamic (enhanced/stronger propagation of SPP as shown in Supplementary S4) process. If we give a close look at Lorentz force in Fig. 3, it is clearly observable that at this QQ resonance: though the resultant binding force is due to the dominance of Del FSurf(x), Del FBulk( also exists. This fact can be verified from Fig. 4(i), where the dominance of on cube-1 is clearly observable. The opposite charges of the cubes create the attractive force between the cubes placed over any substrate such as plasmonic or glass or silicon substrate. At the same time, the propagating surface plasmon of the plasmonic substrate also create force on the cubes in the opposite direction of its propagation, which ultimately causes the local maximum value of attractive force at this QQ or SPP resonance [cf. Fig. 2(l) in main article and Fig. S4 in Supplementary S4]. On the other hand, the strength of SPP becomes extremely weak at DD or LSP resonance as shown in Supplementary S4.

Now we shall consider the behavior of optical force for Fano resonance when the reversal of lateral binding force takes place [cf. Fig. 2(l)]. If we compare the extinction efficiency of the four cases [cf. the second column of Fig. 2], the magnitude is much stronger for the plasmonic (silver) substrate case in comparison with all other cases. In addition, if we consider the strength of induced current at resonant frequencies (i.e. Fano resonance and previously discussed DD resonance mode), the surface current for the case of silver slab is quite stronger in comparison with all other cases: no substrate, glass substrate and silicon substrate case. Now if we give a closer look to the individual cube (i.e. cube-1): it is clearly observable that the bulk force part [cf. Eq. (6) in ‘METHODS’ section] turns into negative force only near the Fano resonance regime for the case of plasmonic substrate. However, this bulk force (which is mainly connected with the scattering force of plasmonic objects[34, 35]) is never dominant, when the substrate is glass or silicon. The total force of cube-1 for glass and silicon substrate case is always dominated by surface force. In contrast, the scenario is fully reversed for the case of silver substrate especially during the Fano resonance. In Fig. 4(i) it can be observed that the reversal of optical binding force occurs due to the dominance of the reversed bulk Lorentz force. Our ultimate conclusion is that: (i) the reversal starts to occur from strong multiple scattering and (ii) though it is commonly believed that internal wave-field does not contribute much for optical force on plasmonic objects (cf. the dominance of surface force in refs 35, 36), during the Fano resonance the scenario is quite different.

Effect of height, size and background material on binding force

In previous section we have discussed the effect of different substrates on lateral binding force. In this section we shall focus our discussion on lateral binding force due to the change of dimer heights from the plasmonic substrate.

In Fig. 5(a)–(f), it is shown that when the cubic silver dimers are far away (45 nm away from the spacer) from the plasmonic substrate, Fano resonance does not occur. For these dimers, placed 45 nm away from the spacer, the dominant resonance in Fig. 5(a) is the DD resonance [but for the isolated dimer case in Fig. 2(a) the dominant resonance is not DD mode[14]]. It should be noted that the binding force at DD resonance mode behaves quite differently for these two cases: (i) Far from the substrate- strong mutual attractive force occurs and (ii) very close [i.e. 5 nm away from spacer as shown in Fig. 2 (j)–(l)] to the substrate- strong mutual repulsive force occurs [the reversal mainly occurs before this dipole resonance at the strong Fano resonance, when stronger coupling with the plasmonic substrate starts to occur].

Figure 5

Two silver cubes (L = 120 nm) are placed at different heights from the silver substrates (three different columns represent the position of the cubes: 50, 20 and 15 nm away from the substrate) and their inter-particle gap, d, is fixed 100 nm. The spacer height is 5 nm. x-polarized plane wave is propagating towards –z-direction. For different heights from the substrate: First row: extinction co-efficient; second row: binding force; ‘°’ represents the chosen wavelengths for which the steady state current distributions are plotted later. Third row: Lorentz force components along with the total force only on cube-1 (left cube, denoted as P1); fourth, fifth and sixth row: steady state current (J y) from front view [in xz plane; setting the window very close to the cube surfaces from front view] for two silver cubes along with the substrate in three different wavelengths (marked as ‘°’). For first column [50 nm height, (d–f)] wavelengths are chosen as: 855, 553 and 503 nm respectively. For second column [20 nm height, (j–l)] resonance wavelengths are: 872 (DD), 631 (Fano) and 567 nm (QQ) respectively. For third column [15 nm height, (p–r)] resonance wavelengths are: 890 (DD), 674 (Fano) and 590 nm (QQ) respectively.

Now, in Fig. 5(g)–(l), when the cubic silver dimers are closely placed (15 nm away from the spacer) from the plasmonic (silver) substrate, Fano resonance does not die out. However, such Fano resonance is not strong enough and optical binding force does not reverse. This configuration can better be comparable with the case: plasmonic homodimers (cube) over the previously discussed Si substrate. The current, electric fields and magnetic fields of such set-ups are also similar (but not identical) to that previously discussed case of dimers over Si substrate. Now, in Fig. 5(m)–(r), when the cubic silver dimers are closely placed (10 nm away from the spacer) from the plasmonic substrate, strong Fano resonance just starts to occur along with strong coupling with the substrate. During this resonance, the bulk part of Lorentz force reverses and just starts to dominate the total force on each dimer [cf. Fig. 5(o)]. Ultimately these lead to the reversal of optical binding force as shown in Fig. 4(n). The mechanism of binding force for this case is very similar to the case [5 nm away from the substrate] discussed previously in details in previous section.

The waveguide description[37] of resonance mode[21, 23] also aids in understanding the observed dependence of the resonance shift on both the sizes of the nano cube [as can be seen from the comparison of Fig. 1(b),(c) and Supplementary S2 Fig. S2(b),(c)] and the gap between the silver cube and the silver substrate [as shown in Fig. 5]. (1) For larger cube, the length of the cavity increases and will therefore support a resonance at a longer wavelength. As a result, the resonance redshifts. (2) Similarly, when the gap between the substrate and the cube gets smaller, the effective refractive index of the cavity mode increases. This effectively lengthens the cavity[21] and results in a resonance condition at a longer wavelength. For case (1) and (2), when the second cube is placed closely to the first cube, the induced stable Fano resonance also red shifts following the rapid red shift of DD resonance and hence the binding force reversal wavelength also red shifts following the Fano dip in Fig. 5 and Fig. 1(b),(c); Supplementary S2 Fig. S2(b),(c).

Throughout the main article, the background material has been considered water. The simplest way to blue shift all the resonances along with the lower wavelength of binding force reversal is to decrease the background refractive index [cf. the results shown in Supplementary S1 and S5 considering air background instead of water]. Last but not least, with the increase of the inter-particle gap, the binding force magnitude decreases slowly [cf. Supplementary S1 (120 nm lengthened cubes) and Supplementary S2 (180 nm lengthened cubes)] but still this overall set-up provides a very relaxed mechanism to verify the reversal of optical binding force. Notably, the conclusions drawn in this work may remain valid for other plasmonic substrates (i.e. Au substrate, as shown in Supplementary S5), too.

Conclusions

In summary, we have investigated a very simple possible configuration to demonstrate the reversal of lateral optical binding force with plasmonic homodimers based on strong Fano resonance. Among all the substrates (i.e. glass, Si, Ag, Au etc.), the closely placed plasmonic particles should remain very close to the plasmonic substrate so that the bulk part of the total Lorentz force dominates the total lateral force of each cube during the substrate mediated Fano resonance. The surface current (along with the strong electric and magnetic coupling) increases significantly only for the close presence of the plasmonic substrate, which bears significant influence on the reversal of binding force of the dimers. Though we have shown only the homodimer cases, our proposed idea should also work for the cube heterodimers providing more flexibility on particle size. We believe that our proposals can be verified experimentally due to the simplicity of the proposed set-ups; and thus the attractive and repulsive forces between two plasmonic objects can be robustly adjusted based on the idea proposed in this work.

Methods

The ‘outside optical force’[38] is calculated by the integration of time averaged Minkowski[38, 39] stress tensor at r = a + employing the background fields of the scatterer of radius a:where ‘out’ represents the exterior total field of the scatterer; , , and are the electric field, displacement vector, magnetic field and induction vectors respectively, 〈 〉 represents the time average and is the unity tensor.

On the other hand, based on the Lorentz force, the total force (surface force and the bulk force[40-42]) can be written as:where is the surface force density (the force which is felt by the bound electric and magnetic surface charges of a scatterer), which is calculated just at the boundary of a scatterer[40-42]. ‘in’ represents the interior fields of the scatterer; ‘avg’ represents the average of the field. σ and σ are the bound electric and magnetic surface charge densities of the scatterer respectively. The unit vector is the outward pointing normal to the surface. is permittivity and is permeability of the background. Bulk is the bulk force density, which is calculated from the interior of the scatterer by employing the inside field[40-42]. ε s is permittivity and μ s is permeability of the scatterer. As per we know, the Lorentz force dynamics of plasmonic particles and specially dimers have not been discussed previously. It is notable that the ‘external dipolar force’[43-48] (which has also been described as Lorentz force in ref. 43) is quite different than the Lorentz force[40–42, 49, 50] defined in our Eqs (4–6).

Supplement article