Casimir forces exerted by epsilon-near-zero hyperbolic materials.

The Casimir force exerted on a gold dipolar nanoparticle by a finite-thickness slab of the natural hyperbolic material namely, the ortorhombic crystalline modification of boron nitride, is investigated. The main contribution to the force originates from the TM-polarized waves, for frequencies at which the parallel and perpendicular components of the dielectric tensor reach minimal values. These frequencies differ from those corresponding to the Lorentzian resonances for the permittivity components. We show that when the slab is made of an isotropic epsilon-near-zero absorbing material the force on the nanoparticle is larger than that induced by a hyperbolic material, for similar values of the characteristic parameters. This fact makes these materials optimal in the use of Casimir's forces for nanotechnology applications.

Introduction

The presence of electromagnetic fluctuations is the origin of important phenomena such as thermal emission, radiative heat transfer, van der Waals interactions, Casimir effect and van der Waals (contact-free) friction between bodies[1] which play an important role in the behavior of matter at very short distances with important implications in nanoscience and nanotechnology. That is why the study of such fluctuations is currently the subject of numerous investigations.

One of the most intriguing effects confirming the foundations of the quantum field theory was predicted by H. Casimir in 1948[2] and was referred in[3] to as “the driving force from nothing”. In the same year, H. Casimir and D. Polder proposed a theory for dipole interactions taking into accound retardation effects[4]. Experimental confirmation of the existence of Casimir forces was carried out M.J. Sparnaay in 1958[5]. The theory of the Casimir forces for real materials at finite temperatures was proposed by E.M. Lifshitz[6] and the general theory of Van-der-Waals forces was developed by Dzyaloshinskii, Lifshith and Pitaevskii[7]. It was subsequently shown that the presence of a liquid in between the two interacting bodies may induce a change in the sign of the force that may become repulsive instead of attractive. The methods used to compute Casimir forces and the value of the force for different geommetries has been reviewed in[8-12].

Futher progress in the study of the Casimir effect and in the Van-der-Waals forces was due to the discovery of new materials such as metamaterials, and to investigations on how micro- and nanoparticles interact with electromagnetic fields. A theory of the Casimir forces under non-equilibrium conditions, for example when the objects are at different temperatures was also developed. These results have been reviewed in[12].

In the last decade, a number of papers on Casimir forces in hyperbolic metamaterials (HMMs) on the basis of arrays of metal nanowires were published. Metamaterials have shown to be very useful in near-field radiative transfer[13,14] and in other areas of fluctuation electrodynamics due to the fact that they exhibit stricking properties such as a high density of modes and an ability to support propagating waves with very large wave-vectors. This effect is due to an increase of the density of evanescent fields of plasmonic modes within the gap.

Casimir forces (attractive and repulsive) can act at large distances[15,16] due to the fact that HMMs may support propagating rather than evanescent waves at large values of the transverse component of the wave vector. The forces are not only perpendicular but also parallel to the surface of the material[17]. Conservative lateral forces induced by corrugations of the surface are local acting only within one of the corrugation periods and a consequence, they can only exert local lateral displacements. On the contrary, non-conservative lateral forces as the ones induced by fluctuating currents in an anisotropic HMM (boron nitride) can move a particle persistently along a flat homogeneous boundary[20]. The lateral movement of anisotropic particles along a surface was analyzed in[18]. The motion, however, ceases when the particle adopts an orientation for which its energy is a minimum. Lateral Casimir forces can also cause particle rotation[19].

The interest in metamaterials with near-zero (NZ) parameters (refractive index, permittivity and permeability) is due to the fact that structures made of these materials offer enormous possibilities for applications[21]. Near-zero-index photonics is currently an area of rapid growth[22]. Structures with near-zero parameters are: epsilon-near-zero (ENZ), , mu-near-zero (MNZ), , and epsilon-mu-near-zero (EMNZ)[22]. All of these cases exhibit a near-zero index of refraction: .

Ziolkovski proposed zero index of refraction metamaterials for the decoupling of spatial and temporal variations of the field[23]. This idea is based on the fact that the field in the region with near-zero parameters exhibits a static-like character even if it oscillates in time. Tunneling of electromagnetic waves through a narrow two-dimensional channel filled with an ENZ material was predicted in[24]. Near-zero-index media can enhance nonlinear processes[25,26], optical activity in one-dimensional epsilon-near-zero pseudochiral metamaterials[27], electric levitation[28] and other field-matter interaction processes.

In this article, we compute the Casimir force on a gold dipolar nanoparticle induced by a slab of an -near-zero hyperbolic and isotropic material and show that its main contribution comes from the frequency domain where .

The article is organized as follows. In Section II, we introduce the model and analyze the different contributions to the radiative force on a gold particle close to a slab of an absorptive anisotropic material. In Section III, we present our results of the Casimir force and in Section IV we summarize our main conclusions.

The model

We consider a slab made of an absorptive anisotropic material which is infinite in the x- and y-directions and has thickness h in the z-direction. The anisotropy axis is directed along the z-axis. We will compute the Casimir forces acting on a gold nanoparticle placed nearby the slab. The relative permittivity tensor of the material has the diagonal formwhere are the coordinate unit vectors.

This particular geometry allows us to analyze separately the propagation of TM and TE waves in the slab. Let us consider the TM modes and find fields excited by point-like fluctuating currents within the slab in the frequency domain. For the tangential field components , excited by the fluctuating currents located within the absorptive layer , the Maxwell equations reduce to the system of two ordinary differential equations:with the matrix elements of [A] given bywhere and Ohm are the wavenumber and wave impedance in vacuum, respectively. The components of the vector F(z) = areThe elementary bulk current source has the form: .

The solution of Eq. (2) for points in the interval is[34]:with the transfer matrix:where and are the transverse and normal components of the wave vector in the slab, respectively, andis the transverse wave impedance for the TM mode. The boundary conditions are: , where is the transverse wave impedance in vacuum. We can then express the tangential field components at the interface , created by a current located at in the formwhereThe average values of the fluctuating currents vanish, only their correlations contribute to the energy flux. These correlations are given through the fluctuation-dissipation theorem[35]. The ensemble-averaged correlator in the plane , for the mode, induced by fluctuating currents located within the slab, , readswhere the correlation is obtained through the fluctuation-dissipation theorem[35]with andthe Planck’s oscillator energy. In Eq. (11), , is the permittivity of vacuum, the reduced Planck constant, T the temperature, and the Boltzmann constant.

We will consider radiative forces on a small nanoparticle in the dipole approximation. The dipolar force acting on the particle can be written asThis formula was derived from the Maxwell stress tensor for the dipolar particle[36]. In it, is the polarizability of the particle given bywith r and its radius and permittivity, respectively, and .

The first term in (13), related to the gradient forces, causes attraction of the particle toward the slab interface due to the z-dependence of the fields through the factor , for . This term properly defines the Casimir or the van-der-Waals forces on a nanoparticle[1]. The explicit expression for iswhere , if , and , if [20]. Only evanescent waves contribute to this force.

In the second contribution, the x- and z-component of the Poynting vector corresponds to pulling forces along the corresponding directions. The x-component of the Poynting vector integrated over is zero in the case of a symmetric geometry and different from zero in the asymmetric case, as shown in[20]. At small |z|, the attractive gradient force is dominant whereas at larger |z| the dominant force is the repulsive force proportional to the z component of the Poynting vector. Very often, this contribution is considered as a part of the Casimir force. The third term in (13) does not contribute to radiative forces[20].

Using the fluctuation-dissipation theorem (11), and expression (15) and integrating over z, as done in[20] [see formulas (17)], we obtainwhere we have defined the coefficientswithThe non-conservative force, determined by the second term of (13), proportional to the Poynting vector, is negligible at small distances from the boundary of an absorptive medium.

The Casimir force exerted by the TM waves is obtained by integrating the force over and :

Results and discussion

Cubic, hexagonal, rhombohedral and orthorhombic cristalline forms of boron nitride exhibit hyperbolic dispersion in the infrared frequency range[29-32]. As a particular case, we will consider the orthorhombic form. The components of the permittivity tensor are given by the Lorentz model[32]:where and are, respectively, the transverse phonon frequency and the oscillator strength of the lattice vibration for the parallel and perpendicular polarizations, and is the damping constant. The constants are the components of the permittivity tensor at frequencies that greatly exceed the phonon resonance frequency . The values of the parameters of (20) used are: , ,  rad/s,  rad/s, , ,  rad/s,  rad/s. For these values of the parameters, the Lorentzian resonances of and take place at frequencies  THz and  THz, respectively. In the vicinity of these resonances, the real parts of and change their signs and the imaginary parts are very large. One can expect an increase of the Casimir forces near the resonance due to the singularity of if . Similarly, an increase of the Casimir force per unit of frequency is expected near the resonance when increases due to the increase of , if [see Eq. (9)].

Figure 1 illustrates the frequency dependence of Re() and Im() on the frequency range where the corresponding permittivity component experiences the Lorentzian resonance. Because of the losses, near the frequency 24.9 THz. Figure 2 shows similar dependencies for the perpendicular component of the permittivity. Here, we see that near 48.5 THz.

Figure 1

Re() (black solid line), Im() (red solid line), and (blue dashed line).

Figure 2

Re() (black solid line), Im() (red solid line), and (blue dashed line).

As an example of particle, we consider a spherical gold nanoparticle whose complex permittivity in the infrared range is, according to the Drude model, given bywhere  rad/s and  rad/s are the plasma frequency and the damping frequency, respectively and [33]. The radius of the particle is 10 nm.

As expected, the main contribution to the Casimir forces comes from the regions where and . Figs. 3 and 4 show the Casimir force per frequency unit in the vicinity of these frequencies. Oscillations in Fig. 4 are caused by the excitation of plasmon-polaritons in the vicinity of -near-zero frequencies[38,39]. For hyperbolic materials becomes very large at frequencies , where . At these frequencies in a finite-thickness slab exists a dense (countable in the lossless limit) spectrum of modes (see[40], Fig. 14) which manifest itself as the ‘fringes’ in Fig. 4.

Figure 3

The Casimir force per frequency unit versus frequency in the vicinity of the -near-zero region.  nm,  K,  nm.

Figure 4

Casimir force per frequency unit versus frequency in the vicinity of the -near-zero region calculated for the thickness  nm (red),  nm (black), and  nm (blue). Dashed blue line shows the Casimir force per frequency unit calculated for an isotropic material with permittivity , at  nm.

Figure 5 shows the overall Casimir force. The lower integration limit over frequency is 10 THz since at smaller frequencies contributions to the Casimir force are very small. The upper limit corresponds to a frequency at the abscisa axis. The result of the integration increases dramatically in the frequency domains where and are minimal. The repulsive force due to the second term in (13) (proportional to the Poynting vector) is of order  N. In the figure , we compare the Casimir forces exerted by the slab of hyperbolic material (boron nitride) to the ones obtained from a hypothetical isotropic material with permittivities and . The greatest Casimir force is induced by the isotropic material with the permittivity undergoing the Lorentzian resonance for in boron nitride.

Figure 5

The Casimir force [N] versus frequency calculated for boron nitride (black), for a material with permittivity (red), and for a material with (blue).

Conclusions

In summary, we have shown that the main contributions to the Casimir forces on a dipolar particle results from the TM-polarized waves and takes place at regions where and are minimal which differ from the regions at which the Lorentzian resonances for the corresponding permittivity components take place. The leading contribution comes from the -near-zero region. Hyperbolicity itself (i.e. different signs of the parallel and the perpendicular components of the permittivity) does not guarantee a high force value, compared to the one obtained for an -near-zero isotropic absorbing material corresponding to the TE-waves and excluding the term with .

Our result of the Casimir force differs from that obtained in the case that both objects are made of the same HMM[37] in which the force is much greater than that obtained between dielectric materials. This fact indicates that the nature of the materials and the frequency dependence of the permittivity could play a role in the value of the force. The effect of TE-waves in the Casimir force was also analyzed arriving at the conclusion that the force is three orders of magnitude smaller than the one resulting from the TM-waves.

The enhancement of the Casimir force found when ENZ hyperbolic materials are used shows that these materials could be advantageous for the use of Casimir’s forces in nanotechnology.