#### Abstract

A study of the electrostatic force between an Electrostatic Force Microscope tip and a dielectric thin film with finite conductivity is presented. By using the Thomas-Fermi approximation and the method of image charges, we calculate the electrostatic potential and force as a function of the thin film screening length, which is a magnitude related to the amount of free charge in the thin film and is defined as the maximum length that the electric field is able to penetrate in the sample. We show the microscope’s signal on dielectric films can change significantly in the presence of a finite conductivity even in the limit of large screening lengths. This is particularly relevant in determining the effective dielectric constant of thin films from Electrostatic Force Microscopy measurements. According to our model, for example, a small conductivity can induce an error of more than two orders of magnitude in the determination of the dielectric constant of a material. Finally, we suggest a method to discriminate between permittivity and conductivity effects by analyzing the dependence of the signal with the tip-sample distance.

#### 1. Introduction

In recent years, Electrostatic Force Microscopy (EFM) has become a useful tool in quantitative characterization of different properties at the nanoscale [1–6]. These technical improvements have been especially relevant in the characterization of thin films [7–10]. Most of the previous work has been focused on perfectly insulating dielectrics [11], but the possible effect of a small conductivity on the EFM signal has not been addressed before (i.e., previously, every layer included in the sample was only characterized by the thickness and dielectric constant). The conductivity, even with small values, is very important in the characterization of thin films in air ambient conditions where water and molecular adsorption and water condensation [12–14] can be detected as a finite surface conductivity due to the ions added on the otherwise insulating sample. Significant finite conductivity effects are also expected for thin films of conducting materials with a finite length since they cannot accumulate enough free charge in the required regions to compensate external electric fields [15]. The presence of defects on the surfaces is also a factor that implies the presence free charges in certain regions that could change the conductivity of the sample [16, 17].

For the discussion of the effects of the thin film conductivity, we are going to introduce the concept of thin film screening length , which is a magnitude directly related to the free charge in the film and can be defined as the maximum length that the electric field is able to penetrate in the sample. This concept is useful in this context because it is a length and it can be directly compared to the length of the thin film itself. For example, conducting thin films, with a screening length much smaller than the film thickness, are relatively easy to characterize since there is a perfect screening of the whole sample composed by the thin film and substrate. This system has been widely studied before [18, 19]. For this reason, in this paper we will focus on the effects of having thin films with small values of conductivity. In the present paper, we will include a theoretical model that includes the possibility of thin films with certain conductivity, which is a physical effect that was not taken into account in previous related articles [11]. This limit is much more difficult to analyse since, due to the relatively large values of the screening length, the electric field is able to penetrate the thin film and the substrate and relevant effects of the electric field distribution must be expected. With the electric field being the most important magnitude in Electrostatic Force Microscopy (EFM), we will analyse the effects of the thin film conductivity on the microscope signal. We will demonstrate that the presence of a small amount of free charge in thin films can be seen as a thin film with an extremely high dielectric constant. This is a very important result since we demonstrated previously that EFM has a very good resolution in those circumstances.

#### 2. Theoretical Model

In a typical EFM setup, we have a metallic tip connected to a battery that applies a constant electric potential . The tip is placed over a sample at a tip-sample distance . The tip usually takes the shape of a cone that is attached to the cantilever in its wider size. In that way, the sharp apex of the tip is placed in the closest size to the sample. It is well known that the macroscopic shape of the tip plays a relevant role in the EFM signal [20, 21], especially when the tip-sample distance is comparable to other magnitudes such as the tip radius. However, the tip cone shape does not affect the field distribution near the tip apex in the very proximity of the sample [22]. In this paper we will focus on tip-sample distances that are between and . Within this limit and in order to understand the main physics associated with sample’s finite conductivity, a simpler tip characterization is adequate. The tip in our simulations is going to be defined by the tip apex radius . The sample conductivity is included through a simple linearized Thomas-Fermi (Debby-Hückel) approximation. In this model, taking into account that , the electrostatic potential can be obtained by solving the following expression:where we are assuming that is uniform inside each region and is the screening length, which is related to the conductivity by , where is the diffusion constant of the material and is the thin film dielectric constant.

In the case of a thin film with finite nonzero conductivity, the sample is going to be composed of two layers, as shown in Figure 1: (1) a thin film with dielectric constant , thickness , and screening length and (2) a dielectric substrate with dielectric constant . When the sample does not have any conductivity (i.e., ) the electrostatic potential of a punctual charge is obtained by the image charge theory, where the image charges are usually obtained by solving the Laplace/Poisson equation in cylindrical coordinates .

**(a)**

**(b)**

**(c)**

**(d)**

From the geometry described in the previous paragraph, we only have inside the thin film. For any other region, and (1) becomes again the Poisson/Laplace equation. Following the scheme shown in Figure 1(d), for every region, the electrostatic potential , for a single point charge located in region at (, ), can be written aswhere , are the first-order Bessel functions, and and are coefficients obtained by applying the electrostatic boundary conditions ( and at ) to the sample interfaces. In the region where the punctual charge is placed (), we have . The other coefficients take the following form:where and are described in Figure 1(d) and can be defined aswhere and . Equation (2) is the general form of the electrostatic potential obtained in cylindrical coordinates when the Laplace/Poisson equation is being solved. This solution is the most adequate for geometries like the one in this paper (see Figure 1) since the cylindrical symmetry around the -axis allows us to easily obtained the specific solution of the electrostatic potential at any region by applying the electrostatic boundary conditions. Equations (3) and (4) give the exact values of the coefficients from regions 0 to 2, which are the only parameters that are unknown in (2). The same procedure can be done for any other geometry with the same axial symmetry. However, for geometries with more than 3 regions, the exact analytical values for the coefficients and become hard to obtain and manage. In these cases, numerical methods such as the transfer or scattering matrix methods, which define the relation between the coefficients of different interfaces by a numerical matrix, have been used successfully in the past [23]. It is worth noting that the forms of (3) and (4) for the particular case of have been previously described in [11]. The generalization of this procedure to the particular cases when there is certain conductivity in the thin film () is one of the main contributions of this paper. These results can be also useful to improve the quantitative description of the conductivity of samples such as graphene layers, where it has been demonstrated that their conductivity can be measured (or at least estimated) by EFM [15].

To solve electrostatic problems with this geometry, an algorithm called the Generalized Image Charge Method (GICM) has been developed [20]. The GICM replaces the surface charge density by a set of charges inside the metallic tip and the sample by a set of image charges. The charge values are obtained by a standard least-squares minimization (LSM). Once the charge distribution is obtained by the tip, the electrostatic force can be directly obtained by the interaction between the tip charges and their images.

#### 3. Effects of Small Thin Film Conductivities

In Figure 1 we show three different equipotential distributions that only differ in the value of the thin film. The EFM setup has been simulated by nm, nm, nm, and V. For simplicity, we have fixed the dielectric constant of both film and substrate to (i.e., the only difference between film and substrate in Figure 1 is the screening length ). When , the thin film-substrate system becomes a semi-infinite dielectric sample with , which corresponds to the equipotential distribution shown in Figure 1(b). As it can be seen, smaller values imply a smaller electric field inside the sample, as it should be. This effect is also found when the thin film dielectric constant increases, as shown in Figure 2(a) where we plot the electrostatic potential drop along the -axis (from the tip apex to the substrate below the film) for the equipotential distributions shown in Figure 1. To compare the effect of having a nonzero conductivity (i.e., finite ) and having a higher value in the thin film, we also show the electrostatic potential drop for (, ) and (, ). These results are particularly relevant since they show that different combinations of dielectric constant and screening length can give rise to the same electrostatic force; for example, , nm) = , ) at nm and , nm) = , ) at nm, as we show in Figure 2(b). Comparing the potential drop of these three samples, we see that finite and high values reduce the electric field inside the sample. However, they have very different effects on the electrostatic potential and on the electric field distribution. Focusing on these three curves, we see that when nm, the electric field (potential gradient) is lower in the region near to the tip (air or vacuum) and higher inside the sample, compared to the cases with high . It is worth noting that as we can see in Figure 2(b), the effects of conductivity can be also distinguishable for values higher than 1 nm. This fact makes the method useful even in air condition, where the minimum working tip-sample distance is limited by the jump into contact distance. In this paper we have assumed nm in most of the cases for clarity in the figures, but higher values could be also used.

**(a)**

**(b)**

In Figure 3 we show as a function of (Figure 3(a)) and (Figure 3(b)) for two different tip-sample distances ( nm and nm). Focusing on nm, we see that does not converge until or . This implies that when working at very small tip-sample distances and small film thicknesses, the EFM is able to distinguish huge values, as it has been reported before [24]. Due to the small thin film thickness, the electric field is able to go through the film and penetrate in the substrate. For higher or values, the electric field vanishes inside the film and no difference can be effectively found between a thin film and a semi-infinite metallic sample. The values of the screening length that can be distinguished also imply that small amounts of conductivity are able to change the electrostatic signal significantly. This fact makes extremely important having under control any possible presence of conductivity in the sample, especially when quantitative magnitudes are being measured. Another parameter that must be taken into account is the tip-sample distance. As we can see in Figure 3, nm) converges much faster than nm). This fact makes the problem even harder since the effect of having a finite screening length is different for different tip-sample distances. Moreover, as we can see in Figure 2(b), nm) = , ) when nm. However, both curves diverge when becomes higher. On the upper limit under study ( nm), the value that corresponds to nm) is . In other words, the presence of conductivity can be confused with insulating thin films characterized by dielectric constants that differ around an order of magnitude, depending on the working tip-sample distance.

**(a)**

**(b)**

In Figure 4 we show the combinations of and that gives the same EFM signals for and 10 nm. As we can see, there is a big window of values that, depending on the working tip-sample distance, are equivalent to a certain value. The window of possible values becomes wider when becomes smaller. For example, the values go from to when nm and from to when nm. Although this fact can be assumed as inconvenient at the first sight, it can be also used to distinguish the physical phenomena involved in the EFM signal by measuring the electrostatic force at two different tip-sample distances. If the apparent dielectric constant of the material differs significantly, we may think that the sample is having a conductive behavior.

It is worth noting that when working at very small tip-sample distances, van der Waals and water meniscus forces may have an influence in the final results and must be taken into account or avoided. In the case of van der Waals force, increasing the tip-sample distance is a good method to avoid its effect since it has a stronger dependence with the distance than the electrostatic forces. In the case of water meniscus, it is better to reduce the tip voltage since their formation is directly related to the electrostatic energy (proportional to the tip voltage) [25].

#### 4. Conclusions

In conclusion, we have developed a method to simulate and study the effects of the conductivity in thin film samples. We have demonstrated that very small values of conductivity can change the electrostatic force in EFM to values similar to those obtained from thin films with extremely high dielectric constants. For example, a thin film with nm and can be compared to thin films without any conductivity and between 110 and 500. We have also demonstrated that the equivalent value can change more than one order of magnitude depending on the working tip-sample distance. The strong dependence with the tip-sample distance can be used to distinguish the origin of the electrostatic signal since, only in the case of having conductivity, different values would be measured at different distances.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors acknowledge J. J. Sáenz, C. Gómez-Navarro, and J. Gómez-Herrero for insightful discussions. This work has been partially funded by the Banco Santander-UAM Research Program. G. M. Sacha acknowledges support from the Spanish Ramón y Cajal Program.