Optical Forces on an Oscillating Dipole Near VO2 Phase Transition

Abstract

We investigate optical forces on oscillating dipoles close to a phase change vanadium dioxide (VO$${}_2$$) film, which exhibits a metal-insulator transition around 340 K and low thermal hysteresis. This configuration emulates the interaction between an illuminated nanosphere and an interface and we employ a classical description to capture its important aspects. We consider both electric and magnetic dipoles for two different configurations, namely with the dipole moments parallel and perpendicular to the VO$${}_2$$ film. By using Bruggeman theory to describe the effective optical response of the material, we show that the thermal hysteresis present in the VO$${}_2$$ transition clearly shows up in the behavior of optical forces. In the near-field regime, the force on both dipoles can change from attractive to repulsive just by heating (or cooling) the film for a selected frequency range. We also verified that the optical forces are comparable to the Casimir-Polder force in a similar system, revealing the possibility of modulating or even changing the sign of the resultant force on an illuminated nano-object due to the presence of a thermochromic material. We hope that this work contributes to set the grounds for alternative approaches to control light-matter interactions using phase-change materials.

1. Introduction

Optical forces play a pivotal role in photonics with many applications. As interesting examples, we can mention radiation pressure forces [1,2], forces in optical tweezers [3,4,5,6], nanostructures [7,8], and waveguides [9,10], as well as interdisciplinary applications in biology [11] and atomic physics [1,5,12,13,14]. Hence, the possibility of harnessing light-matter interactions to tailor and control optical forces at the nanoscale is a sought-after goal in nanophotonics. At the nanoscale range, fluctuation-induced phenomena are non-negligible and Casimir forces between two bodies (or Casimir-Polder forces between a polarizable particle and a body) may be comparable to the classical electromagnetic forces [15,16,17]. In particular, the feasibility of switching on and off the repulsion between particles and surfaces in micro- and nano-mechanical devices can lead to new functionalities, such as the levitation of objects from surfaces to eliminate undesired adhesion and stiction of nanomechanical components [6,18].

Recent advances in plasmonics and metamaterials allow for the development of new material platforms to tune optical and dispersive forces at increasingly smaller scales. Remarkable examples are phase-change materials (PCM) [19], such as transition metal dichalcogenides [20], transparent conductive oxides [21], and liquid crystals [22,23], which have been integrated into metasurfaces and metadevices to allow for external control of their functionalities. Using this strategy, one can either progressively tune or abruptly switch the structural and/or optical properties of metadevices by externally varying an applied voltage [20], electric current [24], and incident light intensity [25]. Amid the PCMs for photonic applications, vanadium dioxide (VO${}_2$) may be singled out for its low thermal hysteresis, for exhibiting a metal-insulator transition (MIT) at low temperature (around 340 K) over a broad frequency range, and for its high refractive-index contrast in the visible range [26,27,28,29,30]. The modulation of plasmonic responses in VO${}_2$ has been reported [31,32,33,34]. In particular, hybrid structures of VO${}_2$ (in combination with metallic nanostrustures) may display localized surface plasmon resonances, a distinct characteristic for the field of active plasmonics [35]. Altogether, these features make VO${}_2$ a promising PCM as a tunable metamaterial [36], enabling the emergence of applications, such as smart windows [37,38]. Photonic applications of VO${}_2$-based structures, both in the infrared [39] and visible ranges [40], have been recently developed, including the temperature control of heat transfer [41] and quantum emission [42].

Due to its tunable optical properties, PCMs have been considered in investigations of the Casimir force in some situations [43,44,45]. In particular, Casimir-Polder forces between an atom and a VO${}_2$ film [46] have already been calculated. In the present work, we investigate optical forces on oscillating electric and magnetic dipoles close to a VO${}_2$ film. This system mimics the interaction of an illuminated nano-object with a neighbouring interface, as such an object would develop electric and/or magnetic dipole moments induced by the incoming radiation. In addition, this setup also relates to an excited two-level quantum emitter located near a surface (in the dipole approximation). Therefore, the classical description alone provides newsworthy results and may not only shed light upon the competition of light-induced and dispersive forces, but also predict important aspects of cavity effects on excited quantum states, a kind of approach that has been widely applied to study optical forces in such systems [47,48,49,50]. We consider two distinct configurations to compute optical forces: In one of them, the oscillating dipole (electric or magnetic) is parallel to the VO${}_2$ film, and, in the other one, it is perpendicular to the film. Remarkably, we show that it is possible to achieve a thermal control of the optical force on the dipole and change the attractive/repulsive character of the force in the near field on both kinds of dipoles just by varying the temperature. Hence, we conclude that thermal hysteresis clearly shows up in optical forces. Our calculations for optical forces are comparable to the Casimir-Polder force in a similar system, so that the resultant force on an illuminated nano-object can be regulated or even have its sign modified due to the presence of a thermochromic material. Our results expand the applicability of PCMs, VO${}_2$ in particular, to the external tuning of optical forces, and in general controlling light-matter interactions at the nanoscale.

This paper is organized as follows. In the next section, we present the methodology employed to model the VO${}_2$ medium and its metal-insulator phase transition and to calculate the classical expressions for the optical force near a planar surface for both configurations of dipoles. Section 3 comprises our main results, whereas Section 4 is dedicated to our final comments and main conclusions.

2. Methodology

In order to study optical forces, we shall consider a prescribed oscillating dipole near a VO${}_2$ film of thickness $d=200$ nm, supported by a saphire substrate (Al${}_2$O${}_3$). The oscillating dipole can be either electric or magnetic, and it is placed at a distance z from the surface, as shown in Figure 1.

In the following subsections, we introduce the theoretical approach employed to characterize the film and its metal-insulator phase transition. Next, we describe the main equations used to evaluate the optical forces on the dipole.

2.1. Modelling the VO${}_2$ Metal-Insulator Transition

In this work, we used the approach developed in Refs. [42,51] to model VO${}_2$ optical properties in terms of its dielectric function, applying the Bruggeman effective medium theory (BEMT) [52]. In short, VO${}_2$ presents an insulating behavior and monoclic structure [53] up to a critical temperature $T_{\mathrm{MIT}}\sim 340$ K, at which it undergoes a metal-insulator transition that consists of a structural phase transition to a rutile-type structure. In practice, however, the transition is smooth and at temperatures close to $T_{\mathrm{MIT}}$, the material presents a mixture of both metallic and insulating domains [54,55], which can be modeled as spheroidal metalic inclusions in a dielectric host medium. It is, therefore, possible to treat VO${}_2$ as a two-phase system with an effective permittivity depending on the permittivities of both metallic and insulating regions and their respective volume filling fractions f and $1-f$ ($0<f<1$) [56]. Here, the filling fractions are modeled as a logistic function of the temperature to emulate the histeretical behavior of the VO${}_2$ electric response as the temperature is ramped up or down [42].

The task now is to calculate the VO${}_2$ effective dielectric constant ${\epsilon }_{{\mathrm{VO}}_2}(\lambda ,T)$ for temperatures in which the heteregeneous regime is present, i.e., close to $T_{\mathrm{MIT}}$. In the framework of BEMT, this may be be obtained from [42]

(1)$\begin{array}{cc}\hfill (1-& f)\left\{\frac{{\epsilon }_{\mathrm{d}}-{\epsilon }_{{\mathrm{VO}}_2}}{{\epsilon }_{{\mathrm{VO}}_2}+L({\epsilon }_{\mathrm{d}}-{\epsilon }_{{\mathrm{VO}}_2})}+\frac{4({\epsilon }_{\mathrm{d}}-{\epsilon }_{{\mathrm{VO}}_2})}{2{\epsilon }_{{\mathrm{VO}}_2}+(1-L)({\epsilon }_{\mathrm{d}}-{\epsilon }_{{\mathrm{VO}}_2})}\right\}\hfill \\ \hfill +& f\left\{\frac{{\epsilon }_{\mathrm{m}}-{\epsilon }_{{\mathrm{VO}}_2}}{{\epsilon }_{{\mathrm{VO}}_2}+L({\epsilon }_{\mathrm{m}}-{\epsilon }_{{\mathrm{VO}}_2})}+\frac{4({\epsilon }_{\mathrm{m}}-{\epsilon }_{{\mathrm{VO}}_2})}{2{\epsilon }_{{\mathrm{VO}}_2}+(1-L)({\epsilon }_{\mathrm{m}}-{\epsilon }_{{\mathrm{VO}}_2})}\right\}=0,\hfill \end{array}$

where L ($0\le L\le 1$) is the depolarization factor related to the shape of the metallic inclusions, and wavelength and temperature dependences have been omitted for simplicity. In addition, ${\epsilon }_{\mathrm{d}}$ and ${\epsilon }_{\mathrm{m}}$ denote the dielectric constants of VO${}_2$ in the purely insulating ($T\ll T_{\mathrm{MIT}}$, $f=0$) and in the purely metallic ($T\gg T_{\mathrm{MIT}}$, $f=1$) phases, respectively. These quantities are given by

(2)$\begin{array}{ccc}\hfill {\epsilon }_{\mathrm{m}}\left(\omega \right)& =& 1-\frac{{\omega }_{\mathrm{m}}^2}{{\omega }^2+i{\gamma }_{\mathrm{m}}\omega },\hfill \end{array}$

(3)$\begin{array}{ccc}\hfill {\epsilon }_{\mathrm{d}}\left(\omega \right)& =& 1+\frac{{\omega }_{\mathrm{d}}^2}{{\omega }_{\mathrm{R}}^2-{\omega }^2-i\omega {\gamma }_{\mathrm{d}}},\hfill \end{array}$

where ${\omega }_{\mathrm{R}}$ is the material resonance frequency, ${\omega }_{\mathrm{m}}\left({\omega }_{\mathrm{d}}\right)$ is the the plasma frequency (oscillating strength), and ${\gamma }_{\mathrm{m}}\left({\gamma }_{\mathrm{d}}\right)$ is the inverse of the relaxation time of the metallic phase (insulator phase). These parameters were obtained by numerically fitting the experimental data reported in Reference [57] in the range 1 μm $\le \lambda \le 10$ μm.

2.2. Optical Forces near a Planar Surface

In this subsection we shall briefly establish the expressions for the optical forces acting on a harmonically oscillating electric or magnetic dipole, when they are near a given planar surface whose reflection Fresnel coefficients for the transverse electric (TE) and transverse magnetic (TM) waves are known. Hence, the electromagnetic fields responsible for the optical force acting on the dipole will be the scattered part of the own dipole field. The optical force acting on a time-dependent electric dipole $\mathit{d}\left(t\right)$ located at a generic position ${\mathit{r}}_0$ is given by [58,59]

(4)${\mathit{F}}^e\left(t\right)=(\mathit{d}\left(t\right)·\nabla )\mathit{E}({\mathit{r}}_0,t)+\frac{1}{c}\dot{\mathit{d}}\left(t\right)\times \mathit{B}({\mathit{r}}_0,t).$

For the case at hand, the dipole and the electromagnetic fields have the same harmonic time-dependence, namely:

(5)$\mathit{d}\left(t\right)={\mathit{d}}_0{e}^{-i\omega t};\mathit{E}(\mathit{r},t)={\mathit{E}}_0\left(\mathit{r}\right)e^{-i\omega t};\mathrm{and}\mathit{B}(\mathit{r},t)={\mathit{B}}_0\left(\mathit{r}\right)e^{-i\omega t}.$

The time average of the force is given by

(6)$\begin{array}{ccc}\hfill ⟨{\mathit{F}}^e⟩& =& ⟨(\mathcal{R}e\mathit{d}·\nabla )\mathcal{R}e\mathit{E}⟩+\frac{1}{c}⟨\mathcal{R}e\dot{\mathit{d}}\times \mathcal{R}e\mathit{B}⟩\hfill \\ & =& \frac{1}{2}\mathcal{R}e\left\{({\mathit{d}}_0·\nabla ){\mathit{E}}_0^{*}-\frac{i\omega }{c}{\mathit{d}}_0\times {\mathit{B}}_0^{*}\right\}.\hfill \end{array}$

For the i-component of the force, we have

(7)$⟨F_i^e⟩=\frac{1}{2}\mathcal{R}e\left\{d_{0j}{\partial }_j{E}_{0i}^{*}-\frac{i\omega }{c}{ϵ}_{ijk}{d}_{0j}{B}_{0k}^{*}\right\}.$

Throughout this paper, we adopted the implict sum notation for repeted indices. Using Faraday’s law, ${\mathit{B}}_0=\frac{c}{i\omega }\nabla \times {\mathit{E}}_0$, so that $B_{ok}^{*}=\frac{ic}{\omega }{ϵ}_{lnk}{\partial }_l{E}_{0n}^{*}$, as well as the identity ${ϵ}_{ijk}{ϵ}_{lnk}={\delta }_{il}{\delta }_{jn}-{\delta }_{in}{\delta }_{jl}$, we finally obtain [59]

(8)$F_i^e\left({\mathit{r}}_0\right)=\frac{1}{2}\mathcal{R}e\left\{d_{0j}{\partial }_i{E}_{0j}^{*}\left({\mathit{r}}_0\right)\right\},$

where we brought back the dipole position ${\mathit{r}}_0$, and, for convenience of notation, we wrote the time average $⟨F_i^e⟩\left({\mathit{r}}_0\right)$ simply as $F_i^e\left({\mathit{r}}_0\right)$.

Analogously, it can be shown that the time average force acting on an oscillating magnetic dipole $\mathit{m}\left(t\right)=\mathcal{R}e\left({\mathit{m}}_0{e}^{-i\omega t}\right)$ is given by $F_i^m\left({\mathit{r}}_0\right)=\frac{1}{2}\mathcal{R}e\left\{{\mu }_0{m}_{0j}{\partial }_i{H}_j^{0*}\left({\mathit{r}}_0\right)\right\}$, so that, if the two dipole moments are present, the total force on the particle is

(9)$F_i^e\left({\mathit{r}}_0\right)+F_i^m\left({\mathit{r}}_0\right)=\frac{1}{2}\mathcal{R}e\left\{d_{0j}{\partial }_i{E}_{0j}^{*}\left({\mathit{r}}_0\right)+{\mu }_0{m}_{0j}{\partial }_i{H}_j^{0*}\left({\mathit{r}}_0\right)\right\}.$

Now, let us consider that the dipole is near a planar surface. In this case, in order to compute the optical force on the dipole, we need to take into account in the previous equations the electromagnetic field that is scattered by the neighboring surface and acts back in the dipole. As usual, these scattered electromagnetic fields can be calculated with the aid of the corresponding Green function $\mathbb{G}(\mathit{r},{\mathit{r}}^{\prime };\omega )$.

Let us first compute the electromagnetic force on an oscillating electric dipole at position ${\mathit{r}}_0$. Adopting $\mathit{m}=0$ in the previous equation the time average of the electric force reads

(10)${\mathit{F}}^e\left({\mathit{r}}_0\right)=\frac{1}{2}\mathcal{R}e\left\{d_j\nabla E_j^{*}\left({\mathit{r}}_0\right)\right\},$

where, for simplicity of notation, we are omitting the subscripts 0 in the electric dipole moment, as well as in the electric field. In order to compute the force on the electric dipole using the previous equation, we need the scattered electric field, which can be written in terms of the scattered Green function as ${\mathit{E}}^{\left(S\right)}\left(\mathit{r}\right)={\omega }^2{\mu }_0{\mathbb{G}}_e^{\left(S\right)}(\mathit{r},{\mathit{r}}_0;\omega )·\mathit{d}$. The scattered Green function for a planar geometry is well-known in the literature [47], so that last equation leads to

(11)$\begin{array}{cc}\hfill F_l^e& =\frac{1}{2}{\omega }^2{\mu }_0\mathcal{R}e{\left\{d_i^{*}{d}_j{\partial }_l{G}_{ij}^{\left(S\right)}(z,z_0;\omega )\right\}}_{z=z_0}\hfill \\ \hfill & =-\frac{1}{4}{\omega }^2{\mu }_0\mathcal{R}e\left\{d_i^{*}{d}_j{\int }_0^{\infty }d{k}_{\Vert }{e}^{2i{k}_{z0}{z}_0}{\int }_0^{2\pi }\frac{d\phi }{{\left(2\pi \right)}^2}\frac{k_{\Vert }{k}_l}{k_{z0}}{R}_{ij}\right\}\hfill \end{array}$

for the l-component of the force, with $k_{z0}=\sqrt{k_0^2-k_{\Vert }^2}$, $k_0=\omega /c$ and $\mathbb{R}$ given by

(12)$\mathbb{R}=\sum _{\mathrm{p},\mathrm{q}=\{\mathrm{TE},\mathrm{TM}\}}{r}^{\mathrm{p},\mathrm{q}}{\mathit{ϵ}}_{\mathrm{p}}^{+}\otimes {\mathit{ϵ}}_{\mathrm{q}}^{-},$

where $r^{\mathrm{p},\mathrm{q}}$ ($\mathrm{p},\mathrm{q}=$ TE, TM) are the usual Fresnel reflection coefficients for a p-polarized incident wave being scattered into a q-polarized reflected wave, and ${\mathit{ϵ}}_{\mathrm{TE}}^{\pm }$ and ${\mathit{ϵ}}_{\mathrm{TM}}^{\pm }$ denote the TE- and TM-polarization unitary vectors, respectively [47]. It can be shown that

(13)$\begin{array}{cc}\hfill F_{\mathrm{x}}^e\left(z\right)& =\frac{1}{8\pi {\epsilon }_0}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }\left[i{r}^{\mathrm{TM},\mathrm{TM}}\mathcal{I}m\left(d_{\mathrm{x}}^{*}{d}_{\mathrm{z}}\right)+\frac{k_0}{k_{z0}}{r}^{\mathrm{TE},\mathrm{TM}}\mathcal{R}e\left(d_{\mathrm{y}}^{*}{d}_{\mathrm{z}}\right)\right]k_{\Vert }^3{e}^{2i{k}_{z0}z},\hfill \end{array}$

(14)$\begin{array}{c}{F}_{\mathrm{z}}^e\left(z\right)=\begin{array}{cc}\hfill -\frac{1}{16\pi {\epsilon }_0}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }[& |d_{\Vert }{|}^2\left(k_0^2{r}^{\mathrm{TE},\mathrm{TE}}-k_{z0}^2{r}^{\mathrm{TM},\mathrm{TM}}\right)+2{\left|d_{\mathrm{z}}\right|}^2{k}_{\Vert }^2{r}^{\mathrm{TM},\mathrm{TM}}\hfill \\ \hfill & +4i\mathcal{I}m\left(d_{\mathrm{x}}^{*}{d}_{\mathrm{y}}\right)k_0{k}_{z0}{r}^{\mathrm{TE},\mathrm{TM}}]k_{\Vert }{e}^{2i{k}_{z0}z}.\hfill \end{array}\hfill \end{array}$

In the previous equations, we have already performed the angular integral in $d\varphi$ and $|d_{\Vert }{|}^2=|d_{\mathrm{x}}{|}^2+{\left|d_{\mathrm{y}}\right|}^2$; $F_{\mathrm{y}}^e$ is the analogous to $F_{\mathrm{x}}^e$ provided we replace the dipole moment components properly, meaning $d_{\mathrm{x}}\to d_{\mathrm{y}}$ and $d_{\mathrm{y}}\to -d_{\mathrm{x}}$.

For an isotropic material, like VO${}_2$, $r^{\mathrm{TE},\mathrm{TM}}=0$. In this work, we will only be concerned with cases in which the dipole components do not present a relative phase between them and can be considered real. Consequently, $\mathcal{I}m\left(d_{\mathrm{x}}^{*}{d}_{\mathrm{z}}\right)=\mathcal{I}m\left(d_{\mathrm{x}}^{*}{d}_{\mathrm{y}}\right)=\mathcal{I}m\left(d_{\mathrm{y}}^{*}{d}_{\mathrm{z}}\right)=0$ and the optical forces acting on the dipoles will have only a component perpendicular to the VO${}_2$ film.

It is convenient to split the vertical force into two contributions, namely: $F_{\perp }^e$, which is proportional to the component of the electric dipole moment perpendicular to the surface $d_{\mathrm{z}}$, and $F_{\Vert }^e$, which is proportional to the component of the electric dipole moment parallel to the surface $d_{\Vert }$. These contributions can be written as

(15)$\begin{array}{cc}\hfill \frac{F_{\perp }^e}{{\mathsf{\Gamma }}_0^e\hslash k_0}\left(z\right)& =-\frac{3}{8k_0^4}\frac{|d_z{|}^2}{{\left|\mathit{d}\right|}^2}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }{r}^{\mathrm{TM},\mathrm{TM}}{k}_{\Vert }^3{e}^{2i{k}_{z0}z},\hfill \end{array}$

(16)$\begin{array}{cc}\hfill \frac{F_{\Vert }^e}{{\mathsf{\Gamma }}_0^e\hslash k_0}\left(z\right)& =-\frac{3}{16k_0^4}\frac{|d_{\Vert }{|}^2}{{\left|\mathit{d}\right|}^2}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }\left(-k_{z0}^2{r}^{\mathrm{TM},\mathrm{TM}}+k_0^2{r}^{\mathrm{TE},\mathrm{TE}}\right)k_{\Vert }{e}^{2i{k}_{z0}z}.\hfill \end{array}$

In order to deal only with dimensionless quantities, we normalized the force by the quantity ${\mathsf{\Gamma }}_0^e\hslash k_0$, where ${\mathsf{\Gamma }}_0^e={\left|{\mathit{d}}_0\right|}^2{k}_0^3/(3\pi {\epsilon }_0\hslash )$ is the spontaneous emission rate of a two-level system in the empty space with transition dipole moment equal to the dipole moment ${\mathit{d}}_0$ (the quantity ${\mathsf{\Gamma }}_0^e\hslash k_0$ can be interpreted as the recoil force on a quantum emitter whose transition frequency is ${\omega }_0$). Recall that, if we take the average in all possible orientations, for an isotropic emitter, $|d_z{|}^2/{\left|\mathit{d}\right|}^2=1/3$ and $|d_{\Vert }{|}^2/{\left|\mathit{d}\right|}^2=2/3$.

Let us now turn our attention to the calculation of the optical force on a magnetic oscillating dipole. The procedure to compute this force follows the same steps as those for the electric dipole case. The main difference is that we need now the magnetic Green function ${\mathbb{G}}_m(\mathit{r},{\mathit{r}}_0;\omega )$, instead of the electric one. With this is mind, the optical force on the magnetic dipole reads

(17)${\mathit{F}}^m\left({\mathit{r}}_0\right)=\frac{{\mu }_0}{2}\mathcal{R}e\left\{m_j\nabla H_j^{*}\left({\mathit{r}}_0\right)\right\}.$

As in Equation (10), we are omitting the subscripts 0 in the magnetic dipole moment and in the magnetic field. The scattered part of the magnetic field at a generic position $\mathit{r}$ created by a magnetic dipole which is located at position ${\mathit{r}}_0$ is given by ${\mathit{H}}^{\left(S\right)}\left(\mathit{r}\right)={\mathbb{G}}_m^{\left(S\right)}(\mathit{r},{\mathit{r}}_0;\omega )·\mathit{m}$. The magnetic Green function can be written as ${\mathbb{G}}_m^{\left(S\right)}(\mathit{r},{\mathit{r}}^{\prime };\omega )={\mu }_0^{-1}\overrightarrow{\nabla }\times \mathbb{G}(\mathit{r},{\mathit{r}}^{\prime };\omega )\times {\overleftarrow{\nabla }}^{\prime }$ [60]. For a planar geometry, we find the expressions for the optical force on a magnetic oscillating dipole in terms of the Fresnel reflection coefficients, namely:

(18)$\begin{array}{cc}\hfill F_{\mathrm{x}}^m\left(z\right)& =\frac{{\mu }_0}{8\pi }\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }\left[i{r}^{\mathrm{TE},\mathrm{TE}}\mathcal{I}m\left(m_{\mathrm{x}}^{*}{m}_{\mathrm{z}}\right)-\frac{k_0}{k_{z0}}{r}^{\mathrm{TE},\mathrm{TM}}\mathcal{R}e\left(m_{\mathrm{y}}^{*}{m}_{\mathrm{z}}\right)\right]k_{\Vert }^3{e}^{2i{k}_{z0}z},\hfill \end{array}$

(19)$\begin{array}{cc}\hfill F_{\mathrm{z}}^m\left(z\right)& =\begin{array}{cc}\hfill -\frac{{\mu }_0}{16\pi }\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }[& |m_{\Vert }{|}^2\left(k_0^2{r}^{\mathrm{TM},\mathrm{TM}}-k_{z0}^2{r}^{\mathrm{TE},\mathrm{TE}}\right)+2{\left|m_{\mathrm{z}}\right|}^2{k}_{\Vert }^2{r}^{\mathrm{TE},\mathrm{TE}}\hfill \\ \hfill & -4i\mathcal{I}m\left(m_{\mathrm{x}}^{*}{m}_{\mathrm{y}}\right)k_0{k}_{z0}{r}^{\mathrm{TE},\mathrm{TM}}]k_{\Vert }{e}^{2i{k}_{z0}z}.\hfill \end{array}\hfill \end{array}$

Comparing the above formulas with those written in Equation (14), we see that they can be mapped one into another when we make the replacements $d_i\to m_i$, $r^{\mathrm{TM},\mathrm{TM}}\to r^{\mathrm{TE},\mathrm{TE}}$, $r^{\mathrm{TE},\mathrm{TE}}\to r^{\mathrm{TM},\mathrm{TM}}$, and $r^{\mathrm{TE},\mathrm{TM}}\to -r^{\mathrm{TE},\mathrm{TM}}$.

Finally, as we have done for the optical force on the oscillating electric dipole, it is also convenient to split the vertical force into two contributions, to wit: $F_{\perp }^m$ (proportional to the component of the electric dipole moment perpendicular to the surface $m_{\mathrm{z}}$) and $F_{\Vert }^m$ (proportional to the component of the electric dipole moment parallel to the surface $m_{\Vert }$). For an isotropic medium, they are given by

(20)$\begin{array}{cc}\hfill \frac{F_{\perp }^m}{{\mathsf{\Gamma }}_0^m\hslash k_0}\left(z\right)& =-\frac{3}{8k_0^4}\frac{|m_z{|}^2}{{\left|\mathit{m}\right|}^2}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }{r}^{\mathrm{TE},\mathrm{TE}}{k}_{\Vert }^3{e}^{2i{k}_{z0}z},\hfill \end{array}$

(21)$\begin{array}{cc}\hfill \frac{F_{\Vert }^m}{{\mathsf{\Gamma }}_0^m\hslash k_0}\left(z\right)& =-\frac{3}{16k_0^4}\frac{|m_{\Vert }{|}^2}{{\left|\mathit{m}\right|}^2}\mathcal{R}e{\int }_0^{\infty }d{k}_{\Vert }\left(-k_{z0}^2{r}^{\mathrm{TE},\mathrm{TE}}+k_0^2{r}^{\mathrm{TM},\mathrm{TM}}\right)k_{\Vert }{e}^{2i{k}_{z0}z},\hfill \end{array}$

where we applied an analogous normalization to the one used in Equations (15) and (16), but now with ${\mathsf{\Gamma }}_0^m={\mu }_0{k}_0^3{\left|\mathit{m}\right|}^2/(3\pi \hslash )$ being the spontaneous emission rate of a magnetic emitter in empty space.

3. Results and Discussions

We now proceed to detailed discussions of our results regarding the optical forces acting on oscillating dipoles close to a VO${}_2$ film. As previously mentioned, we performed our analysis for the cases of electric and magnetic dipoles, and each of the following subsections accounts for one of them.

3.1. Electric Dipole

In Figure 2, we plot the force on the oscillating electric dipole in the perpendicular configuration, due to the presence of VO${}_2$ medium near the MIT, as a function of temperature. We chose $z=50$ nm, so that we are in the near-field regime, and the electric dipole is oscillating perpendicularly to the film. As expected, the material thermal hysteresis is directly reflected in these curves. The most compelling feature that can be noted is the fact that there are some values of wavelengths $\lambda$ for which the attractive/repulsive character of the force may be interchanged just by heating or cooling the VO${}_2$ film.

In order to furtherly exploit this property, Figure 3 displays the force dependence with the wavelength $\lambda$ of the radiation emitted by a dipole at the same distance $z=50$ nm. Note that, due to the thermal hysteresis, two distinct temperatures can correspond to the same curve in Figure 3 depending on whether the system is being heated or cooled. It can be spotted that the change in the force sign with temperature occurs only for the short range $1.186$ μm $<\lambda <1.228$ μm (where the blue and the purple lines cross $F_{\perp }^e=0$). For higher wavelengths, the force is always attractive, whereas, for smaller wavelengths, always repulsive. It may also be noted from Figure 3a that, as the wavelength increases, the force does not significantly change with $\lambda$ for a fixed temperature. This behavior can also be inferred from Figure 2c, which reveals that the relative change in $F_{\perp }^e$ at a given phase reduces with $\lambda$. This remark should be considered carefully. We normalized Equations (15) and (16) by ${\mathsf{\Gamma }}_0^e\hslash k$ to keep our results as general as possible, so that they do not explicitly depend on any particular value of the dipole moment $\left|\mathit{d}\right|$ (just on the ratio $d_z/\left|\mathit{d}\right|$). However, in doing so, ${\mathsf{\Gamma }}_0^e\hslash k\sim 1/{\lambda }^4$ and our normalized results artificially grow with ${\lambda }^4$. To avoid any misleading conclusions, in Figure 2c and Figure 3a, we plot the non-normalized $F_{\perp }^e$ considering large $\lambda$ intervals so that this artificial $\lambda$-dependence is absent. In order to do so, we assumed a dipole moment $d_0\sim 3\times {10}^{-27}\mathrm{C}·\mathrm{m}$. This corresponds to a transition wavelength of 1 ¼m and a static polarizability of $4\pi {ϵ}_0{r}^3$, with $r=10$ nm. Indeed, the Casimir-Polder force between a metallic sphere with a few nanometers radius and a metallic wall at $z=50$ nm can be estimated as $F_{\mathrm{CP}}\sim {10}^{-15}$ N [17], comparable with our results from Figure 3.

The profile at Figure 3 can be depicted in terms of the near-field regime of Equation (15). The quasi-static limit, obtained by taking $c\to \infty$, leads us to

(22)$\begin{array}{ccc}\hfill F_{\perp }^e\simeq -\frac{1}{16\pi {\epsilon }_0}\frac{3d_{z,\Vert }^2}{4z^4}{\eta }^e,& \hfill & \hfill {\eta }^e=\frac{|{\epsilon }_{{\mathrm{VO}}_2}{|}^2-1}{|{\epsilon }_{{\mathrm{VO}}_2}{+1|}^2}.\end{array}$

Figure 4 shows the dependence on the factor ${\eta }^e$ with $\lambda$ for different temperatures. By comparing it with Figure 3, one realizes that it exhibits exactly the same behavior, except for the negative sign. In particular, the force changes sign at wavelengths for which $|{\epsilon }_{{\mathrm{VO}}_2}|\sim 1$, i.e., $\lambda \sim 1.2$ μm. Around this value, there is a peak on the factor ${\eta }^e$, responsible for the non-monotonic behavior with $\lambda$ appearing in Figure 2a. Note that, for wavelengths $\lambda \gtrsim 2$ μm, this factor is basically independent of $\lambda$ and its value in the dielectric and metallic phases differs only by a fixed value. These features can be explained due to the behavior of Drude and Drude-Lorenz permittivities ${\epsilon }_{\mathrm{m}}$ and ${\epsilon }_{\mathrm{d}}$ (Equations (2) and (3)) and it is also present in Figure 3b. The permittivities approach a fixed value as $\omega \to 0$. Similar conclusions occur for the case in which the electric dipole oscillates parallel to the film. In fact, the near-field regime of Equation (16) furnishes $F_{\perp }^e\simeq 2F_{\Vert }^e$, so that $F_{\Vert }^e$ only differs from the perpendicular contribution by a factor of two. This factor can be physically interpreted in terms of the field created by a static electric dipole. At a given distance from the dipole, this field has its modulus along the dipole direction twice as large as its modulus at the same distance along a direction perpendicular to dipole.

Figure 5 displays the normalized force $F_{\perp }^e/{\mathsf{\Gamma }}_0^e\hslash k$ as a function of z at the insulator and metal phases for different wavelengths $\lambda$. Note that the force oscillates around zero for large distances and diverges with $\sim 1/z^4$ in the near-field regime. Curiously, for certain wavelengths $\lambda \sim 1.2$ μm, the force may diverge positively or negatively depending on the phase in which the VO${}_2$ film is found. For instance, in the insulator phase (Figure 5a), $F_{\perp }^e$ diverges positively for $\lambda =1.2$ μm, while, in the metallic phase (Figure 5b), it diverges negatively. Conversely, the force $F_{\perp }^e$ for the others $\lambda$ shown in Figure 5 does not change qualitatively with the temperature.

In order to unveil such behavior, Figure 6 illustrates how the force $F_{\perp }^e$ varies as a function of the distance z for different temperatures at fixed wavelengths $\lambda \sim 1.2$ μm. For example, considering $z\lesssim 50$ nm in Figure 6b,c, $F_{\perp }^e$ can be either attractive or repulsive depending on the temperature, but the same behavior is not verified in the other panels. One may notice that, as $\lambda$ increases, the force at very small distances changes gradually from repulsive to attractive. Hence, the force may diverge negatively or positively as $z\to 0$ depending on $\lambda$ and also on the temperature. Actually, this type of behavior was already expected from Figure 2b, when we verified that at a fixed distance $z=50$ nm the force could change sign just through heating for wavelenghts around $\lambda \sim 1.2$ μm. Nevertheless, in studying this change of behavior as a function of z, one may find results that are robust with respect to the distance z. In particular, Figure 6 reveals a new property, to wit, the fact that the first equilibrium point of the system changes its stability character with the temperature.

To further investigate this feature, Figure 7 shows the normalized position $z_{\mathrm{eq}}/\lambda$ of the first equilibrium point as a function of $\lambda$ at fixed temperatures. Note that, for $1.17$ μm $<\lambda <1.21$ μm, $z_{\mathrm{eq}}/\lambda$ varies discontinuously with $\lambda$, and the position of the discontinuity depends on the temperature. In fact, this discontinuous behavior is a consequence of the fact that the equilibrium point changes from unstable to stable as $\lambda$ increases. Moreover, for any $\lambda$ in between the aforementioned range, it is always possible to place the dipole at a distance z at which the force will change sign just by heating (or cooling) the VO${}_2$ film. For instance, the force $F_{\perp }^e$ on the dipole is repulsive at a distance $z=50$ nm for $\lambda =1.18$ μm, regardless of the temperature (see Figure 2b and Figure 6b), but it changes sign for smaller z. In contrast, the force is always attractive for $\lambda =1.25$ μm, regardless of the distance z and the temperature of the film.

We also investigate the electric field distributions of the scattered field (Figure 8) at a given instant of time. Surely, the field distribution is qualitatively altered when we compare insulator and metallic phases, regardless of the $\lambda$ considered. Recall that, at $z=50$ nm, the force $F_{\perp }^e$ on the dipole changes from repulsive to attractive through heating for $\lambda =1.2$ μm, but it is always attractive for $\lambda =1.3$ μm (Figure 2a). In terms of the field distributions, this fact can be inferred from the relative field orientations at the dipole position. The orientation changes from insulator to the metallic phase for $\lambda =1.2$ μm, but it is essentially the same for $\lambda =1.3$ μm. We stress that it is not the absolute orientation of the field that matters for unveiling the character of the force—after all, the fields oscillate in time—but the fact that there is a relative orientation between the two phases that are flipped.

It is also interesting to explore the changes of the electric field intensities exhibited in Figure 8. Recalling that the electric dipole is pushed to regions where the field intensiy is greater, it can be realized that, in the situations depicted in Figure 8b–d, the dipole is attracted to the VO${}_2$ film in agreement with Figure 6. Although the dipole is repelled by the VO${}_2$ film in the situation described in Figure 8a, this is not seen so clearly in this panel.

In Figure 9, we calculated the scattered electric field distribution of an oscillating electric dipole perpendicular to the VO${}_2$ film in the insulator phase with $\lambda =1.2$ μm, but now, in the far-field regime. In this region, we checked that there is no qualitative changes between the distributions in metallic and insulator phases, as expected, since the more distant the dipole is from the VO${}_2$ film, the less important are the boundary conditions and, consequently, the less important is the difference between metals and insulators. Note the transverse character of the field and its oscillatory behavior with the distance, as expected. Observe that local minima are separated by a distance $\lambda /2$, and the field is attenuated for large distances.

3.2. Magnetic Dipole

Figure 10 represents the magnetic force on the magnetic dipole as a function of temperature for different values of $\lambda$. Besides the thermal hysteresis, note that for $\lambda \sim 1$ μm the force behavior is significantly different for perpendicular and parallel configurations (Figure 10a,c), but it is not for $\lambda \gtrsim 2.5$ μm (Figure 10b,d). Moreover, likewise the electric dipole case, there are wavelengths for which the force changes sign just by heating the VO${}_2$ film. However, the wavelength interval for which this may occur is much larger here.

Figure 11 shows the magnetic force as a function of $\lambda$ in both parallel and perpedicular configurations. Comparing Figure 11b,c, one may note that, in the interval 1 μm $<\lambda \lesssim 2.5$ μm, the forces for the perpendicular and parallel configurations of the dipole are qualitatively different. Particularly, the force may change from repulsive to attractive with heating in the interval $1.22$ μm $<\lambda <1.31$ μm for a magnetic dipole parallel to the surface (see Figure 11d), but it is always repulsive in this range for a perpendicular configuration. Nevertheless, heating changes the force from attractive to repulsive in a much broader range than for the electric dipole case, i.e., $\lambda \gtrsim 2.5$ μm, regardless of the dipole orientation (specifically, for a parallel configuration, this happens for $\lambda \gtrsim 2.3$ μm).

In order to further clarify these results, we also performed an extreme near-field ($c\to \infty$) approximation on Equations (20) and (21), obtaining

(23)$\begin{array}{cccc}\hfill F_{\perp }^m& \simeq -\frac{{\mu }_0}{16\pi }\frac{2m_z^2{k}^2}{4z^2}{\eta }_{\perp }^m,\hfill & \hfill & {\eta }_{\perp }^m=\frac{\mathcal{R}e{\epsilon }_{{\mathrm{VO}}_2}-1}{4},\hfill \end{array}$

(24)$\begin{array}{cccc}\hfill F_{\Vert }^m& \simeq -\frac{{\mu }_0}{16\pi }\frac{m_{\Vert }^2{k}^2}{4z^2}{\eta }_{\Vert }^m,\hfill & \hfill & {\eta }_{\Vert }^m=\frac{\mathcal{R}e{\epsilon }_{{\mathrm{VO}}_2}-1}{4}+\frac{|{\epsilon }_{{\mathrm{VO}}_2}{|}^2-1}{|{\epsilon }_{{\mathrm{VO}}_2}{+1|}^2}.\hfill \end{array}$

By comparison of these expressions with Equation (22), one can realize that, besides the factor 2 also present in the electric dipole case, the difference between parallel and perpendicular cases have now an additional term ${\eta }_{\perp }^m$. To clear up these differences, the factors ${\eta }_{\perp }^m$ and ${\eta }_{\Vert }^m$ are plotted in Figure 12. Note that, for $\lambda \gtrsim 2.5$ μm, the contribution ${\eta }_e=\left(\right|{\epsilon }_{{\mathrm{VO}}_2}{|}^2-1)/|{\epsilon }_{{\mathrm{VO}}_2}{+1|}^2$ only shifts up the force in the parallel configuration. This behavior is in agreement with the plot of ${\eta }^e$ in Figure 4b. In fact, a comparison between Figure 10b,d determines that, apart from a shift and a factor, the terms in parentheses in Equations (23) and (24) have qualitatively similar behaviors. Despite that, for 1 μm $<\lambda \lesssim 2.5$ μm, the effect of ${\eta }^e$ is non-negligible. This accounts for the differences in the predictions of parallel and perpendicular configurations in this interval (see Figure 10a,c). Specifically, $F_{\Vert }^m$ is non-monotonic with $\lambda$ in the metallic phase and, more importantly, it may change sign with heating in the interval $1.22$ μm $<\lambda <1.31$ μm (Figure 10c), which does not happen for $F_{\perp }^m$ (Figure 10a).

The dependence of the force on the magnetic dipole with the term involving ${\eta }_{\perp }^m$ explains why the effect of changing the force sign is more robust in frequency for a magnetic dipole than for an electric dipole. Strictly speaking, the lower the frequency, the better is the distinction between metal and dielectric [42]. The distinction between the force in the metallic and insulator phases is clearer for a magnetic dipole than for an electric dipole since in the former the dependence on $\mathcal{R}e{\epsilon }_{{\mathrm{VO}}_2}$ is stronger (see Equations (22)–(24)). The only exception is found for higher frequencies where the distinction between metal and dielectric is blurred. Even though, in this region, we were still able to find a short $\lambda$-interval in which a change in the sign of the force on the magnetic dipole occurs. However, the main reason for that is different, relying on the contribution ${\eta }^e$ with $\mathcal{R}e{\epsilon }_{{\mathrm{VO}}_2}\sim -1$ and small $\mathcal{I}m{\epsilon }_{{\mathrm{VO}}_2}$.

In addition, Equations (23) and (24) reveal a scaling law with the distance in the form $z^{-2}$. Comparison with $z^{-4}$ scaling law for the electric case, given by Equation (22), evidences a two powers difference between the electric and magnetic setups. It may be attributed to the fact that the electric field of an oscillating electric dipole $\mathit{d}={\mathit{d}}_0{e}^{-i\omega t}$ has three terms ($1/r$, $1/r^2$ and $1/r^3$), whereas the electric field of an oscillating magnetic dipole $\mathit{m}={\mathit{m}}_0{e}^{-i\omega t}$ has only two ($1/r$ and $1/r^2$) [61]. The reason for that is very simple: a static magnetic dipole does not create an electric field, so that there is no term proportional to $1/r^3$ in the expression of the electric field produced by an oscillating magnetic dipole. This fact also has interesting consequences in the dispersion force between an electrically polarizable atom and a magnetically polarizable one, namely: while, in the former case, the non-retarded force between them is proportional to $1/r^7$, in the latter, it is proportional to $1/r^5$ [62,63,64]. Likewise, the magnetic field of an oscillating magnetic dipole has three terms, while the magnetic field of an oscillating electric dipole has only two. In other words, the ultimate reason for such two powers difference is related to the fact that the magnetic field of an electric oscillating dipole does not contain the static-like term. Moreover, the scaling law $z^{-2}$ of the force on a magnetic dipole compared to $z^{-4}$ law of force on an electric dipole accounts for the main reason why the electric contribution is dominant in this near-field regime.

Figure 13 shows the force on a magnetic dipole as a function of z for two distinct wavelengths: $\lambda =2$ μm, for which the force is repulsive in the near-field; and $\lambda =4$ μm, for which the force changes sign in the near-field (recall Figure 10b). Besides the expected oscillatory behavior for large distances, note that, for $\lambda =2$ μm, the force $F_{\perp }^m$ diverges positively as $z\to 0$, regardless of the temperature, while, for $\lambda =4$ μm, there are temperatures for which the force diverges negatively as $z\to 0$. In fact, the latter behavior is present for $\lambda \gtrsim 2.5$ μm. In contrast with the electric case, there are no wavelengths for which the force on the magnetic dipole diverges only negatively as the distance z decreases for different temperatures.

These conclusions about whether the force diverges positively or negatively as $z\to 0$ can also be driven from Figure 14, that shows the position of the first equilibrium point as a function of the wavelength $\lambda$. Similarly to the electric case, the discontinuity on the plot informs about the stability of the equilibrium point. For $\lambda <2.3$ μm, the first equilibrium point is always unstable. As $\lambda$ increases, its stability depends on the temperature. Interestingly, for $T>342$ K, it is always unstable. This result is in agreement with the prediction that there is a $\lambda$ above which the force on the magnetic dipole in the perpendicular configuration always changes its sign with heating, i.e., $\lambda \sim 2.3$ μm. From the previous discussion, we can verify that, in order to have the possibility of changing the attractive/repulsive character of the force by varying the temperature, the position of the dipole and its oscillating frequency need to be properly chosen.

Lastly, we computed in Figure 15 the magnetic field distributions for a perpendicular magnetic dipole at a given instant of time. Our considerations are very similar to the ones made before about the electric field distributions for an electric dipole (Figure 8). When the force does not change its sign, the field at the dipole position also does not alter its orientation, as occurs, for example, for $\lambda =2$ μm. In contrast, when the force changes from attraction to repulsion with heating, it can be seen through the modification on the field orientation on the particle position (for $\lambda \gtrsim 2.5$ μm). In particular, for $\lambda =2$ μm, the magnetic field at the dipole position points downwards in both insulator and metallic phases (Figure 15a,b), while, for $\lambda =4$ μm, it points upwards in the insulator phase (Figure 15c) and downwards in the metallic phase (Figure 15d). Regarding the analysis of the field intensities, we should point out that, differently from the electric case, it is not possible to discern the sign of the force as the gradients are too smooth within the panel’s resolution.

In Figure 16, we calculated the magnetic field distribution in the far-field regime (compare with Figure 9). In this region, there is no qualitative changes between the behavior of the distributions in metallic and insulator phases. Note that the separation between two local minima is again given by $\lambda /2$, and the field is weakened and oscillates with the distance.

4. Final Remarks and Conclusions

In this work, we investigated optical forces acting on oscillating electric and magnetic dipoles close to a VO${}_2$ phase-change film, in order to explore the effects of its well-known metal-insulator transition and thermal hysteresis at relatively low temperatures. Additionally, we considered two distinct configurations of the dipoles in our setup, to wit, perpendicular and parallel to the film. We verified the remarkable possibility of thermally controlling the optical force on both dipoles in the near-field regime, once the appropriate dipole frequency and the distance from the VO${}_2$ surface were chosen. More specifically, we show that the thermal hysteresis allows for a change in the attractive/repulsive character of the force just by heating (or cooling) the VO${}_2$ film. Interestingly, the force on electric and magnetic dipoles shows different power laws with the distance from the film. Furthermore, for magnetic dipoles, there are more wavelength intervals for which the force may change its sign depending on the temperature. Our results also show that the optical forces are comparable to the Casimir-Polder force on a nanosphere located at the same distance in a similar setup, so that the resultant force on an illuminated nano-object may be controlled or even reversed due to the presence of the VO${}_2$ film. Altogether, we hope that our results further expand the broad spectrum of applications of these materials, providing alternative ways to tune light-matter interactions using phase-change materials.

Author Contributions

Conceptualization, F.A.P. and W.J.M.K.-K.; methodology, D.S. and W.J.M.K.-K.; software, D.S. and W.J.M.K.-K.; validation, P.P.A.; formal analysis, D.S., F.S.S.R., C.F. and W.J.M.K.-K.; data curation, D.S.; writing—original draft preparation, D.S., F.S.S.R. and C.F.; writing—review and editing, P.P.A., F.A.P. and C.F.; visualization, D.S. and P.P.A.; supervision, C.F. and W.J.M.K.-K.; project administration, F.S.S.R. and C.F.; funding acquisition, C.F., F.S.S.R. and W.J.M.K.-K. All authors have read and agreed to the published version of the manuscript.

Funding

C.F. and F.S.S.R. acknowledge Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for financial support (grant numbers 310365/2018-0 9 and 309622/2018-2). F.S.S.R. (grant number E26/203.300/2017) and P.P.A. acknowledge Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ). D.S. and F.A.P. also acknowledge the funding agencies. W.K.-K. acknowledges the Laboratory Directed Research and Development program of Los Alamos National Laboratory for funding under Project No. 20210327ER.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:

MIT	metal–insulator transition
PCM	phase-change material
BEMT	Bruggeman effective medium theory

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. Schematic representation of the system. An oscillating dipole in vacuum is placed at a distance z of a VO2 film of thickness d suported by a Al2O3 substrate.

Figure 2. (a,b): Force on the electric dipole in the perpendicular configuration normalized by Γ0eℏk as a function of temperature T for different values of wavelengths λ. The dipole is assumed to be at a distance z=50 nm from the VO2 film. In (c), we considered a dipole moment d0∼3×10−27C·m.

Figure 3. (a) Force on the electric dipole (with dipole moment d0∼3×10−27C·m) in the perpendicular configuration as a function of wavelength λ for different temperatures T. The dipole is assumed to be at a distance z=50 nm from the VO2 film. (b) Details of the previous panel in the range 1 μm <λ<1.4 μm. Dashed lines indicate the λ-interval for which the force may change its sign depending on the temperature.

Figure 4. (a) Factor ηe as a function of wavelengths for different temperatures. (b) Details of the previous panel in the range 1 μm <λ<1.5 μm.

Figure 5. Force F⊥e in the perpendicular configuration normalized by Γ0eℏk as a function of the distance at the insulator (a) and metallic phases (b) for different wavelengths.

Figure 6. Force F⊥e in the perpendicular configuration as a function of the distance for different temperatures and fixed wavelengths λ=1.10 μm (a), λ=1.18 μm (b), λ=1.20 μm (c), and λ=1.25 μm (d). The insets in (a,d) show the intercept of these plots with x-axis.

Figure 7. The normalized position zeq/λ of the first equilibrium point as a function of λ for different temperatures. Dotted lines indicate discontinuities in zeq(λ) for each temperature.

Figure 8. Scattered electric field distribution in the near-field regime for an oscillating electric dipole at x=0 and z=50 nm (white point in the plots) perpendicular to the VO2 film in insulator (a,c) and metallic phases (b,d), considering the dipole wavelength λ=1.2 μm (a,b) and λ=1.3 μm (c,d). The VO2 film lies in the region −200 nm <z<0 nm and the sapphire substrate lies in the region z<−200 nm. The electric field is normalized by its maximum value.

Figure 9. Scattered electric field distribution in the far-field regime for an oscillating electric dipole at x=0 and z=50 nm perpendicular to the VO2 film in insulator phase for λ=1.2 μm. The electric field is normalized by its maximum value in this region.

Figure 10. Force on an oscillating magnetic dipole as a function of temperature for different wavelengths in both perpendicular (a,b) and parallel configurations (c,d) for distinct wavelengths. The dipole is assumed to be at z=50 nm from the VO2 film.

Figure 11. Force on an oscillating magnetic dipole as a function of λ in both perpendicular (a,b) and parallel configurations (c,d) for distinct temperatures. (b,d) show detailed intervals of (a,c), respectively. The dipole is assumed to be at a distance z=50 nm from the VO2 film. Dashed lines indicate λ-intervals for which the force may change its sign depending on the temperature.

Figure 12. (a) Factors η⊥m (dashed lines) and η‖m (solid lines) as functions of λ for different temperatures. (b) Details of the previous panel in the range 1 μm <λ<3 μm.

Figure 13. Force on the magnetic dipole F⊥m as a function of distance z, for different temperatures and for (a) λ=2 μm and (b) λ=4 μm.

Figure 14. The normalized position zeq/λ of the first equilibrium point of the magnetic dipole and VO2 film system as a function of λ for different temperatures. Dotted lines indicate discontinuities in zeq(λ) for each temperature. Note that, for T>342 K, there are no discontinuities in the range considered, and the first equilibrium point is always unstable.

Figure 15. Scattered magnetic field distribution in the near-field for an oscillating magnetic dipole at x=0 and z=50 nm (white point in the plots) perpendicular to the VO2 film in insulator (a,c) and metallic phases (b,d), considering the dipole wavelength λ=2 μm (a,b) and λ=4 μm (c,d). The VO2 film lies in the region −200 nm <z<0 nm, and the sapphire substrate lies in the region z<−200 nm. The magnetic field is normalized by its maximum value in this region.

Figure 16. Scattered magnetic field distribution in the far-field regime for an oscillating magnetic dipole at x=0 and z=50 nm perpendicular to the VO2 film in insulator phase for λ=3 μm. The magnetic field is normalized by its maximum value in this region.