## Height distribution of equipotential lines in a region confined by a rough conducting boundary

## C. P. de Castro

Instituto de F´ ısica, Universidade Federal da Bahia, Campus Universit´rio da a Federa¸ c˜o, Rua Bar˜o de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil a a

E-mail:

caioporto@ufba.br

## T. A. de Assis

Instituto de F´ ısica, Universidade Federal da Bahia, Campus Universit´rio da a Federa¸ c˜o, Rua Bar˜o de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil a a

E-mail: thiagoaa@ufba.br

## C. M. C. de Castilho

Instituto de F´ ısica, Universidade Federal da Bahia, Campus Universit´rio da a Federa¸ c˜o, Rua Bar˜o de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil a a Instituto Nacional de Ciˆncia e Tecnologia em Energia e Ambiente, Universidade e Federal da Bahia Campus Universit´rio da Federa¸ c˜o, 40170-115, Salvador, BA, a a Brazil

E-mail:

caio@ufba.br

## R. F. S. Andrade

Instituto de F´ ısica, Universidade Federal da Bahia, Campus Universit´rio da a Federa¸ c˜o, Rua Bar˜o de Jeremoabo s/n, 40170-115, Salvador, BA, Brazil a a

E-mail:

randrade@ufba.br

## Abstract.

This work considers the behavior of the height distributions of the equipotential lines in a region confined by two interfaces: a cathode with an irregular interface and a distant flat anode. Both boundaries, which are maintained at distinct and constant potential values, are assumed to be conductors. The morphology of the cathode interface results from the deposit of 2 × 10 4 monolayers that are produced using a single competitive growth model based on the rules of the Restricted Solid on Solid and Ballistic Deposition models, both of which belong to the Kadar-Parisi-Zhang (KPZ) universality class. At each time step, these rules are selected with probability p and q = 1 -p . For several irregular profiles that depend on p , a family of equipotential lines is evaluated. The lines are characterized by the skewness and kurtosis of the height distribution. The results indicate that the skewness of the equipotential line increases when they approach the flat anode, and this increase has a non-trivial convergence to a delta distribution that characterizes the equipotential line in a uniform electric field. The morphology of the equipotential lines is discussed; the discussion emphasizes their features for different ranges of p that correspond to positive, null and negative values of the coefficient of the non-linear term in the KPZ equation.

Height distribution of equipotential lines in a region confined by a rough conducting boundary 2

PACS numbers: 68.55.-a, 68.35.Ct, 81.15.Aa , 05.40.-a

## 1. Introduction

Growth phenomena in non-equilibrium conditions [1, 2] are an important subject of condensed matter physics. Many properties of specific devices depend on rough surfaces that are produced under such non-equilibrium conditions. Technological applications such as conducting-field emitter devices that operate in Ultra High Vacuum (UHV) conditions, require the control of several properties such as the surface morphology. It is well known that large values of the current intensity on field emitter devices are connected to extremely small values of the effective emitting area [3]. One method to increase the effective area of emission is to use surfaces with an ending fractal boundary with self-affine scale invariance instead of considering structurally regular tips. In fact, the connection between the morphology of irregular surfaces including sharp conducting tips and the field emission properties has been a subject of intense research [4, 5]. In a recent work [6], Cabrera et al. measured experimentally the current-voltage (I-V) of a tunnel junction consisting of a sharp electron emitting metallic tip at a variable distance ('∆') from a planar collector and emitting electrons using electric field assisted emission. Their results showed scale invariance of the tunnel junction with respect to changes in ∆ (from few mm to several nm ), which means that the physical laws governing the flow of current are invariant with respect the changes of the length scale ∆. However, this behavior fails when ∆ was only a few nanometers as showed in reference [7]. The reported scale invariance for a single tip can be regarded as a preliminary study of the more general problem of field emission by an irregular conducting interface that will be considered in this work.

In the statistical mechanics framework, macroscopic laws that emerge in surface growth can be explained at different spatial and temporal scales considering a theoretical microscopic description that involves simple probabilistic laws. An interesting related point is the possibility of modeling these systems so that in the large-scale limit (i.e., where the scale invariance arises), they do not depend on the corresponding details but only on symmetries and the corresponding conservation laws.

Previous theoretical works related to the emitting properties of rough surfaces [8, 9, 10, 11, 12] have studied the behavior of the electric potential in a region bounded by a rough conducting profile/surface and a smooth conducting line/plane, which is held to a constant voltage bias. In particular, the references [9, 10] analyzed the roughness exponent of the equipotential lines (surfaces) in 2 (3)-dimensional systems. However, to our knowledge, no previous work discussed the properties of equipotential lines/surfaces under the perspective of its height distribution properties.

In this work, we investigate the height distribution properties of a family of equipotential lines in a region that is confined by a rough profile, which is regarded as an electron-emitting cathode, and a flat anode which is placed sufficiently away from it. The profile results from a numerical simulation of the deposition and growth surface models, whereas the cathode and the anode are subject to an electric potential difference. Our results are mainly based on the model that was introduced by Silva and Moreira [13, 14], which is compared to those obtained for the simpler Family model [15] in some limiting cases. We emphasize two consequences of our results for the following aspects: the height distribution presents a more detailed characterization of the equipotential lines than that provided by the roughness exponents and such investigation of equipotential lines offers a much more reliable method to characterize the surface. For particular applications, our work can contribute to elucidate the image

distortions in probe microscopies such as Electrostatic Force Microscopy (EFM), where the knowledge of the electric field distribution along the tip probe is relevant to extract morphological information about the real irregular surface. This main conclusion from our work is particularly useful to study surfaces under the influence of non-linear mechanisms in the growth phase, where the surface width is not expressed as a power law with a single growth exponent for a large number of deposited monolayers.

This investigation is also motivated by the quoted experimental aspects related to the electronic properties of rough surfaces and the recent theoretical advances regarding the solution of the Kardar-Parisi-Zhang (KPZ) equation [16] in d = 1+1, where, in general, d = n + 1 indicates that a system composed by a n-dimensional substrate (here, 1-dimensional) and an extra dimension to where an important effect produced by the substrate is propagated.

The remainder of the paper is organized as follows. In Sec. 2, we describe the main features of the used models and emphasize the conditions under which they reproduce the typical features of the KPZ and Edwards-Wilkinson (EW) [17] growth dynamics. In Sec. 3, we describe the methodology to calculate and characterize of the equipotential lines. In Sec. 4, we discuss our results for the morphological properties of equipotential lines and analyze the corresponding behavior of the skewness and kurtosis with the electric potential. In Sec. 5, we summarize the results and present our conclusions.

## 2. Models

It is convenient to discuss the properties of the computational surface growth models that are used in this work with those of the KPZ equation in 1 + 1 dimensions

$$\partial _ { t } h ( x, t ) = \mu _ { 0 } + \nu \triangle h + \frac { \lambda } { 2 } ( \nabla h ) ^ { 2 } + \eta ( x, t ),$$

where h x, t ( ) represents the height of the profile in the y direction with respect to the line, h x, t ( = 0) = 0, that describes the flat profile at t = 0. Eq. 1 comprises several terms that consider the essential features of a stochastic surface growth. These features include irreversibility and locality, which enhance the surface roughness. In eq. 1 µ 0 corresponds to a constant driving force, ν is the diffusion coefficient of the deposited atoms on the surface, and η x, t ( ) is a white Gaussian noise that mimics the stochastic nature of the growth process. The non-linear term accounts for a considerable number of experimental results [18, 19] for the height distribution. Many of them cannot be described by a Gaussian function, which is the solution of the linear EW equation, which corresponds to setting λ = 0 in eq. 1.

The investigated model depends on a parameter p , which describes the probability that the deposited particle follows the Restricted Solid on Solid (RSOS) local rule, whereas q = 1 -p is the probability related to the Ballistic Deposition (BD) rule. Both RSOS and BD discrete models belong to the universality class that is defined by the KPZ equation but are characterized by negative and positive values of the parameter λ , respectively. Although the used model has not been proved to be equivalent to a KPZ equation, recent numerical investigations suggest that a change in the probability p changes the value of λ in a KPZ equation, which may be suitable to describe the competitive model [20, 21].

For almost all ranges of p , a growing profile first experiences a transient phase before attaining a stationary growth regime, where its width W scales with the

number of monolayers with the typical exponent β KPZ = 1 3 of the KPZ universality / class. For finite profiles, this regime is interrupted when W exponentially relaxes to a saturation regime and becomes constant. In this work, we analyze only rough interfaces in the growth regime, which corresponds to typical experimental conditions. In the first transient phase, depending on p , the growth exponent can reach values close to those typical of EW dynamics. Starting from p = 0, the length of the transient phase increases until a specific value p ∗ ≈ 0 83, where the EW typical value is always . valid. For p &gt; p ∗ , the model reaches KPZ scaling again after a decreasing transient time [22]. Because the condition p = p ∗ is very special, we find it wise to compare the properties of the profiles and equipotential lines with those obtained when the cathode is represented by Family-model produced profile [15], which shows the EW properties notably clearly.

The interesting behavior of the competitive model makes it suitable to investigate how the equipotential lines change when p is selected to reproduce either KPZ or EW dynamics, and for a fixed value of p , how the magnitude of the non-linear terms affects the EW-KPZ crossover as a function of time [22], which is particularly important for experimental interests. Indeed, the competition between physical mechanisms is present in many real processes in thin-film science. Then, the actual study clearly advances on the previously one where single models were considered, and the roughness exponent of equipotential lines were evaluated without a perspective for height distribution properties.

The distribution of height fluctuations of the KPZ equation, for a flat initial condition, can be expressed by Gaussian Orthogonal Ensemble (GOE) Tracy-Widom (TW) distributions [23]. Therefore our results can be better appreciated by comparing the skewness ( S ) and kurtosis ( K ) of the equipotential height distributions with the typical TW values for the corresponding quantities. Both S and K depend on the average distance from the equipotential lines to the cathode and anode. Such dependence has intrinsic features according to whether p &lt; p ∗ or p &gt; p ∗ . A third distinct behavior is observed at p = p ∗ , when the profile falls into an EW universality class.

We emphasize again that our investigation opens the possibility of analyzing the effects on the behavior of equipotential lines of distinct surfaces with growth exponents near those that feature an EW class. This possibility is important because of the emergence of local growth exponents in the first transient growth phase, which is near the phase that characterizes the EW dynamics when p ≈ p ∗ . As shown in Sec. 4, our results indicate that the height distributions of equipotential lines for such transient patterns are clearly different from those obtained for an actual EW profile , no matter whether they are produced using the competitive model at p = p ∗ or the Family model.

## 3. Methodology: Determination and Characterization of Equipotential Lines

The current investigation starts by constructing the rough profile h x ( ), which acts as the cathode and results from one of the quoted deposition models. Here, x and y are integer numbers that represent distance measured in terms of some basic unit distance u , the precise size of which is not important in a theory of the kind under discussion. The growth process along the y direction of the considered 1 + 1 system (one-dimensional substrate + height, as mentioned in Sec.1) starts with a flat substrate

h x, t ( = 0) = 0, where x is an integer number, and x ∈ [1 , L ]. The column height h x ( ) is the y coordinate of the topmost adatom at position x . We have used profiles with lateral size up to L = 10 6 and up to T = 2 × 10 4 deposited monolayers (ML). For such parameter values, finite-size effects no longer influence the interpretation of our results and can be neglected. During the growth process, the number of monolayers is also referred to as the integration time t ; in this way, each time unit corresponds to the deposition of L particles. For the convenience of using a unified notation for all profiles and equipotential lines, let us switch from h x ( ) to h x, t ( ) = h p,φ ( x, t ), where we emphasize the dependence of h x ( ) on t . The subscript p stands for the probability of choosing the RSOS rule in the competitive model. Finally, the subscript φ is used to identify the constant value of the equipotential line. In this work, φ = 0 and φ = A correspond to the cathode and anode surfaces, respectively.

The flat anode, which is defined by h p,A , is placed at a distance 〈 d 〉 away from the cathode. 〈 d 〉 is measured by the number of vertical spacings between the average

$$\text{and cause} \, \underbrace { y } \, \text{in some or very one or even} \, \text{with} \, \text{with} \, \text{if} \, \text{the profile} \, \overline { h } _ { p, 0 } = \frac { 1 } { L } \sum _ { x = 1 } ^ { L } \left [ h _ { p, 0 } ( x, T ) \right ] \text{ and the flat anode, i.e:} \\ \langle d \rangle = h _ { p, A } - \overline { h } _ { p, 0 }.$$

$$\langle d \rangle = h _ { p, A } - \overline { h } _ { p, 0 }.$$

As previously indicated, h p, 0 is evaluated at t = T , except if explicitly indicated.

To compare the results for different rough profiles, we fixed the value of 〈 d 〉 . As a consequence h p,A , which is measured with respect to the substrate where the film is grown, changes with p . In this way, to select a proper and fixed value of 〈 d 〉 , such that for any value of p the anode is sufficiently away from the cathode, we first determined the value of p for which the roughness W (defined in eq. 4) attains its maximum, which is p = 0 when W ≈ 37 3. . Next, we identified the vertical coordinate of the highest peak h p =0 max for this value of p and set h 0 ,A = h p =0 max +∆, which result, using eq. 2, in:

$$\langle d \rangle = h _ { m a x } ^ { p = 0 } + \Delta - \overline { h } _ { 0, 0 },$$

where ∆ corresponds to the distance defined on the second paragraph of Sec.1. At this point, let us comment on the reasons to choose the particular value of ∆ used in this work. For ∆ &gt;&gt; h p =0 max , the results are expected to be largely independent of the precise value of ∆, but the computational work is largely increased. We have observed that, for all the rough profiles, the equipotential lines become flat and parallel to the x axis very rapidly after the value of y surpasses the highest peak. Therefore, we choose the value ∆ = 10, which corresponds to 〈 d 〉 = 198. We use this average distance for all p , because it is already sufficient large to have very smooth equipotential lines (with roughness exponent, α , near to unity [10]) that are observed at the flat anode. If we consider the basic unit distance u = 5 nm, we obtain ∆ = 50 nm for p = 0, and 〈 d 〉 /similarequal 1 µ m, which correspond to the typical scale used in the non-contact mode in atomic force microscopy (like EFM) where the electrostatic force is probed. Then, the results of this work can motivate experimental realizations in this regime.

Observed at this angle, the geometrical set up depends on two different scales (expressed by L and H ) along the horizontal and vertical directions. However, the rough profile has its own length scale in the vertical direction, which is the roughness

$$W _ { p, 0 } ( L, t ) = \left [ L ^ { - 1 } \sum _ { x = 1 } ^ { L } \left ( h _ { p, 0 } ( x, t ) - \overline { h } _ { p, 0 } ( t ) \right ) ^ { 2 } \right ] ^ { 1 / 2 }$$

Height distribution of equipotential lines in a region confined by a rough conducting boundary 7

$$\text{where} \, \overline { h } _ { p, 0 } ( t ) = L ^ { - 1 } \sum _ { \substack { x = 1 \\ x = 1 } } ^ { L } [ h _ { p, 0 } ( x, t ) ]. \\ \text{When} \, \ t \ = \ T, \ \text{the profiles that} \\ \text{are in the assumntotic growth response}$$

When t = T , the profiles that are grown using the competitive model are in the asymptotic growth regime ( W ∼ t β KPZ ) for p ∈ [0 , 0 75) . ⋃ (0 9 . , 1]. This regime corresponds to typical experimental conditions, in contrast to the exponential relaxation that is observed in steady-state (long-time) properties. For p ∈ [0 75 . , p ∗ ) ⋃ ( p , ∗ 0 90], the asymptotic growth KPZ regime has not yet been reached, . and the system is in the EW-KPZ crossover [24, 25]. For p = p ∗ ≈ 0 83, the system . is in the critical point and the corresponding growth regime is such that W ∼ t β EW , where β EW = 1 4 is the EW growth exponent. /

The next step is to compute the equipotential lines, h p,φ ( x ). Hence, we numerically solve Laplace's equation

$$\nabla ^ { 2 } \phi = 0,$$

/negationslash for the electric potential φ in the region between the two conductors, which are held at constant potential values of φ C and φ A for the rough cathode and the flat anode that is parallel to the x axis, respectively. The equation is solved using Liebmann's method and the appropriate Dirichlet conditions. The geometrical characterization of the equipotential lines is insensitive to the choice of φ C and φ A = φ C = 0. However, to directly compare some experimental devices, we choose the typical values φ C = 0 and φ A = 10 . 3

Then, the domain is divided into a two-dimensional grid, and the potential is iteratively calculated at each grid point for fixed boundary conditions at the rough cathode, which is described by the function h p, 0 ( x ) and the anode at h p,A . Periodic boundary conditions are imposed on the lateral borders. To determine the coordinates of the equipotential lines, we define the following difference,

$$D _ { \phi } = \frac { \phi _ { x, y + 1 } - \phi _ { x, y } } { \Delta y } > 0,$$

where x, y are integer numbers, and ∆ y =1 /lessmuch L . For any equipotential line h p,φ ( x ) where φ x,y &lt; φ &lt; φ x,y +1 , the corresponding coordinates are computed by evaluating h p,φ ( x ) = h x ( ) + dh φ , with

$$d h _ { \phi } = D _ { \phi } ^ { - 1 } [ \phi - \phi _ { x, y } ].$$

To illustrate the geometry of our problem, Figure 1 shows the substrate (or cathode) and a family of equipotential lines that were calculated using Eqs. (5), (6) and (7) for p = 0 (a) and p = 1 (b).

Next, we characterize the height distributions of the corresponding equipotential lines. In particular, the skewness, S φ ( ) = L -1 L =1 h p,φ ( x ) -h p,φ 3 /W 3 and kurtosis

$$\begin{array} { c } K ( \phi ) = \left [ L ^ { - 1 } \sum _ { x = 1 } ^ { L } \left ( h _ { p, \phi } ( x ) - \overline { h } _ { p, \phi } \right ) ^ { 4 } / W ^ { 4 } \right ] - 3 \, \text{are used throughout this work. Here, } W \\ \text{is a shorthand notation for } W _ { p, 0 } ( L, t ), \, \text{which is defined in } Eq. \, ( 4 ). \end{array}$$

∑ ( x )

## 4. Results and Discussion

Before discussing the results for the Family and the competitive models, we emphasize that in this last case, our results are based on a single profile for each distinct value of p . This usage is justified by the following facts: for each profile, our procedure amounts

/s32

/s32

Figure 1. Geometry of the problem for a substrate that was formed using an irregular profile, which was defined by (a) p=0 and (b) p=1. The equipotential lines h p,φ ( x ) for φ = 50 100 300 500 and 900 are shown. , , , The top panel shows the equipotential line h p, 900 ( x ) in a magnified scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-9e96c63e-8e4e-4710-b8c6-0352a9cd20a5.jpg)

to considering an entire set of equipotential lines, which demands a considerably numerical effort; we have found that, in the interval p /lessorsimilar 0 75 and . p /greaterorsimilar 0 9, . at time T , the height distribution of the cathode and the resulting equipotential lines are almost insensitive to p , which can be observed by comparing the values of S and K . Therefore, our procedure is equivalent to replacing several realizations for one value of p with a set of independent realizations at slightly different values of this parameter. This method can also be justified when we move to the close neighborhood of p = p ∗ . Here, S and K vary smoothly even for the results that are produced using 1 profile realization.

The first illustrated aspect in Figure 2 refers to the relative positions of the profile with respect to the anode, which can be expressed by the ratio 〈 d 〉 W ( 〈 d 〉 = 198) as a function of p . In units of the roughness W , the curve in Figure 2 represents the height where the anode is maintained above the average height of the corresponding rough

/s32

/s32

Figure 2. Ratio 〈 d 〉 W as function of probability p . 〈 d 〉 is the average distance between the anode and the irregular cathode (see Eq.2) while W ≡ W p, 0 ( L, t ) is the roughness of the irregular cathode (see Eq.4). The results are presented considering 〈 d 〉 =198, and it is possible to observe that the global roughness decreases at the interval 0 &lt; p &lt; 0 83 and increases for . p &gt; 0 83. . The global roughness is minimum at the point p = 0 83( . ≈ p ∗ )

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-b2e90e8d-035b-43bf-b033-8595d49d7727.jpg)

.

profile. It is possible to observe that the smoother rough profile occurs at p = 0 83, . i.e., near the critical point. Moreover, it is expected that in the infinite-size limit, the peak becomes more pronounced.

Next, we characterize the height distribution of the equipotential lines based on the behavior of S and K as a function of φ . These quantities vary in a wide interval, so that it becomes convenient to draw some of the graphs in the logarithmic scale. Because both quantities may assume negative values, we shift S and K by a fixed value 0 5, which leads to the functions . S ∗ ( φ ) = S φ ( ) + 0 5 and . K ∗ ( φ ) = K φ ( ) + 0 5 . illustrated in Figure 3. To better interpret the results, we draw the horizontal lines S ∗ = 0 5 and . K ∗ = 0 5, which delimitate the regions of negative and positive values . of S and K . The overall behavior of S ∗ ( φ ) is characterized by its monotonic increase with respect to the value of φ , as displayed in panel (a). This result is consistent with the convergence of S to a δ function, which is expected for a flat equipotential line. Furthermore, for p /lessorsimilar 0 75, the results in Figure 3 reveal that . S φ ( /lessorsimilar 10) ∼ 0 29 . (a) and K φ ( /lessorsimilar 10) ∼ 0 16. . Therefore, in the growth regime where the fluctuations in the surface height distribution is described using the GOE-TW distribution, the equipotential lines close to the profile follow a similar height distribution. The values of S ∗ and K ∗ do not significantly change until φ ≈ 10. However, for φ /greaterorsimilar 10, the derivative dS/dφ significantly grows, which indicates a non-uniform convergence to the δ function. This aspect can be observed by the analysis of Figure 3 (c), which shows the scaled height distribution P { [ h p,φ ( x, T ) -h T ( )] /W } of the equipotential lines for φ = 1 5 50 and 500. , , In fact, our results show that the large deviations of h T ( ) ≡ h p,φ ( T ) for each equipotential line determine the behavior of S . For p = 0 83, . the behavior of S φ ( ) has similar but less pronounced features, which indicates a weaker dependence of dS/dφ at this value of λ .

For p &gt; 0 90, . S φ ( ) becomes negative for h p,φ ( x ) near the rough conducting profiles h p, 0 ( x ). In this region, the absolute value of S is also near that typical of a GOE TW

/s41

/s40

/s42

/s49/s48

/s49

/s83

/s48/s46/s49

/s48/s46/s49

/s49

/s49/s48

/s49/s48/s48

/s49/s48/s48/s48

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-6db06489-1eef-4d60-a9e2-8c18842af796.jpg)

/s32

/s40/s120/s44/s84/s41/s45/s104/s40/s84/s41/s93/s47/s87/s125

/s112/s44

/s80/s123/s91/s104

/s48/s46/s49

/s48/s46/s48/s49

/s49/s69/s45/s51

/s91/s104

/s112/s61/s48/s46/s50/s48

/s45/s52

/s45/s50

/s48

/s50

/s40/s99/s41

/s52

/s61/s49

/s61/s53/s48

/s61/s53

/s61/s53/s48/s48

/s54

/s112/s44

/s40/s120/s44/s84/s41/s45/s104/s40/s84/s41/s93/s47/s87

Figure 3. Shifted skewness S ∗ (a) and kurtosis K ∗ (b) for a family of equipotential lines h p,φ ( x ) as a function of φ for different values of p . Three distinct regimes can be identified: for λ &gt; 0 ⇔ p = 0 0 20 and 0 50, , . . λ ≈ 0 ⇔ p = 0 83, and . λ &lt; 0 ⇔ p = 0 90 0 95 and . , . p = 1 0. . The inset details the behavior in a region where the rough profile morphologies are blurred because of finite time effects. The horizontal dashed line defines S ∗ ( φ ) = 0 5 . ⇒ S φ ( ) = 0 and K ∗ ( φ ) = 0 5 . ⇒ K φ ( ) = 0. (c) Scaled height distributions for the equipotential lines considering φ = 1 5 50 and 500 for a rough conducting profile grown with , , p = 0 20. . The vertical dashed line defines [ h p,φ ( x, T ) -h T ( ) ] /W = 0.

distribution. The important aspect is that S φ ( ) crosses the null value at φ ≈ 15, which coincides with the point where K ∗ ( φ ) exhibits a minimum value (see Figure

/s32

/s48/s46/s49

/s32/s112/s61/s48/s46/s56/s51

/s49

/s49/s48

/s32

/s49/s48/s48

/s32

/s49/s48/s48/s48

/s49/s48

/s49

/s48/s46/s49

/s32

/s40/s97/s41

/s32

/s32

/s32

/s32

/s32

/s112/s61/s48

/s32

/s112/s61/s48/s46/s50/s48

/s32

/s112/s61/s48/s46/s53/s48

/s32

/s112/s61/s48/s46/s56/s51

/s32

/s112/s61/s48/s46/s57/s48

/s32

/s112/s61/s48/s46/s57/s53

/s32

/s112/s61/s49

/s32

/s32

/s32

/s32

Figure 4. Dependence of the local growth exponent β t ( )( ≡ d ln W/d ln t ) for p = 0 8 and . p = 0 825. . The dashed and solid lines indicate the exact EW and KPZ values.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-8df7946e-63ae-44af-b879-b9f81b8b13bb.jpg)

3 (b)). This result indicates that, for conducting surfaces where λ &lt; 0, there is an equipotential line near the surface that presents a symmetric asymptotic behavior at the tails of the height distributions. This result suggests that the height distribution of the equipotential lines changes in a different manner (for positive and negative fluctuations) with respect to that of the actual conducting profile.

In the insets of Figures 3 (a) and 3 (b), we focus on the behavior of S ∗ ( φ ) and K ∗ ( φ ) for p = 0 80, 0 81, 0 87 and 0 89, where for . . . . L and t = T , the system presents morphologic features in the crossover EW-KPZ. Interestingly, the results indicate that similar features do not appear for the values of S ∗ ( φ ) and K ∗ ( φ ). We notice no appreciable differences for t = T , in the values of S ∗ ( φ ) and K ∗ ( φ ) for p ∈ [0, 0.75] and p ∈ (0.75, 0.83) or p ∈ (0.83, 0.9) and p ∈ [0.9, 1]. Our results suggest that the potential variation does not feel the consequences of these finite time effects (KPZEW crossover), following the typical values for models that obey the KPZ statistics. The analysis of equipotential lines appears as a better estimator to uncover the actual geometric features of the growing system.

Let us proceed with a closer discussion of our results for p ∈ (0 75 . , 0 83) and . p ∈ (0 83 . , 0 9) at times where the EW features are believed to exist. . It is well known [24, 25] that in this region, the duration of the transient phase EW-KPZ increases quickly. More specifically, it is claimed that the time to reach the asymptotic limit (i.e. the KPZ class) scales with λ -4 [22]. In competitive models, λ grows at least linearly with the difference p -p ∗ [26], and the crossover time is expect to depend on p with ( p -p ∗ ) -4 . This result is corroborated by evaluating the local growth exponent β t ( ) ≡ d ln W/d ln t (or effective growth exponent) illustrated in Figure 4. Let us first recall that, for any growing profile, a single growth exponent β can only be defined if a strict linear dependence between log( W L,t ( )) and log( ) is observed. t Because this dependence is not verified for a large interval of values of t when 0 75 . /lessorsimilar p /lessorsimilar 0 9, we . analyze the behavior of β t ( ), which is computed by numerically evaluating the local slope of this curve. We consider two conditions for which linear features should prevail in the morphology of conducting rough profiles as p = 0 8 and . p = 0 825. . Figure 4

shows the time evolution of β t ( ), which indicates the presence of a minimum that hints the value of t where the profile should be closer to one grown using a pure EW dynamics (this minimum was estimated at the times t = 347 and t = 1652 for p = 0 8 . and p = 0 825, respectively). .

To show that the analysis of equipotential lines helps elucidating the fine differences between profiles, we also explicitly consider a rough profile that is grown using the well-know linear Family model. Then, we compare the resulting behavior of the corresponding equipotential lines with those for the competitive model for p = 0 83 . ≈ p ∗ and p ∼ 0 83 ( . p = 0 8 and . p = 0 825 at times . t = 347 and t = 1652, respectively). In Figure 5 (a) and (b), we observe a similar behavior for S ∗ and K ∗ , respectively, as a function of φ for the Family model and as p = 0 83 . ≈ p ∗ . This result is interesting because modeling the condition p = 0 83 provides vacancies in . the volume, but these features are not reflected in the behavior of the equipotential lines, as we can observe by directly comparing with the results obtained using the Family model. However, for p ∼ 0 83, we do not observe a similar behavior for the . equipotential lines, which now exhibit different morphologies from the previous ones. This result indicates that β t ( ) of the competitive model close to those that characterize the EW dynamic (minimum of β t ( )), is not a reliable measure to confirm the electrical features of the rough devices that are similar to that grown using the EW dynamics. Moreover, we observe that for p = 0 8, . S ∗ exhibits low values for far equipotential lines, compared with the other models ( S ∗ ∼ 1 for φ = 900). Finally we can observe that a notable convergence for the limit p → p ∗ and for values of t that correspond to the minimum of β t ( ), to the behavior of S ∗ and K ∗ of equipotential lines, which are associated to EW features. At these values of t , the misinterpretations of small differences in β t ( ) compared to β EW can significantly affect the design of electrical devices, where roughness is an important tool.

At this point, let us remark that the height distribution of the equipotential lines becomes a reliable criterion to analyze the morphological differences that are produced by distinct rough cathodes. Indeed, the previous results, which were only based on the local roughness exponent (see reference [9]), have not produced such clear-cut differences in the local roughness exponent of the equipotential lines that are produced by rough cathodes generated by a ballistic ( λ &gt; 0) and random deposition with surface relaxation (RDSR) ( λ = 0). In the current work, such differences become clear when we consider the height distributions of equipotential lines for all three regimes, which are characterized by λ ( &gt; 0 = 0 , , &lt; 0) and p ( &lt; p , ∗ = p , &gt; p ∗ ∗ ).

Further morphological aspects of the profile have been revealed by our results. Let us define the set of peaks P p,φ in an equipotential line. A point x ∈ P p,φ , where its associated vertical coordinate is larger than those of its first lateral neighbors, i.e., h p,φ ( x -1 , T ) &lt; h p,φ ( x, T ) &gt; h p,φ ( x +1 , T ). In a similar way, let us define F + p,φ , the set of peaks of an equipotential line that are also characterized by positive fluctuations. A point x ∈ F + p,φ when x ∈ P p,φ and h p,φ ( x, T ) &gt; h p,φ ( T ). Finally, let us define f + as the fraction of the number of points x ∈ P p,φ that also satisfy x ∈ F + p,φ . It is possible to show that a relation h 3 p,φ ( T ) ∼ ( h + p,φ ( T )) 3 f p should be valid, where h + p,φ ( T ) is the average value of h p,φ ( x, T ) that was restricted to points x ∈ F + p,φ .

This regime identifies the regions for equipotential lines where the small fluctuations quickly disappear and large fluctuations with lower wavelengths prevail when the electric potential grows. In fact, for equipotential lines that are far from the rough profile, we find that h + p,φ ( T ) &gt;&gt; h p,φ ( T ) (For example, the main peaks in

/s32

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-db6f0f11-87aa-41b1-9687-e082ae2ab65a.jpg)

/s32

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-f3cef363-f422-4559-97d5-8666df357a3f.jpg)

/s32

Figure 5. Shifted (a) skewness S ∗ and (b) kurtosis K ∗ of equipotential lines as functions of the electric potential for the rough profiles that were grown according to the Family model and the competitive model at p = 0 8, . p = 0 825 and . p = 0 83 . ≈ p ∗ . For p = 0 8 and . p = 0 825, . the rough cathode surfaces were considered at times when the β t ( ) is minimum (See Figure 4). This minimum was estimated when t = 347 and t = 1652 for p = 0 8 and . p = 0 825, respectively. . It is possible to observe the convergence to EW behavior in the p → p ∗ limit.

Figure 1 (a) for φ = 300 and 500). Therefore, the corresponding global roughness can be written as

$$W ^ { 3 } \sim f _ { + } ^ { - 3 / 2 } ( \overline { h } _ { p, \phi } ^ { + } ( T ) ) ^ { 3 }.$$

Additionally, S should depend on f + as:

$$S \sim f _ { + } ^ { - 1 / 2 },$$

so that S &gt;&gt; 1 for equipotential lines that are far from the cathode. This result is observed in Figure 3 (a), where S ∗ ∼ 10 as φ = 900 for p &lt; 0 83. . It is worth noting that such behavior is observed in the large roughness limit ( p &lt;&lt; 1) of the conducting profile. In Figure 6, we show that these theoretical predictions are confirmed by numerical calculations for the rough conducting profile that were generated for p = 0 1 . and 0 2. . It is clearly observed that in the regime where h + p,φ ( T ) &gt;&gt; h p,φ ( T ), the exponent of f P tends to -1 / 2. In this region, dS/dφ quickly grows for 300 &lt; φ &lt; 500 immediately before the inflexion point (see Figure 3 (a)).

/s32

/s32

/s32

Figure 6. Skewness as a function of the fraction of peaks for equipotential lines that are calculated in the region bounded by a rough profile and the flat line for p = 0 10 0 20 and 0 95. . , . . In the inset, we show a magnified region where the scale relation S ∼ f -1 / 2 P holds.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0421v1-figure-28dcf5dd-f7bd-447f-999f-c80d745998c4.jpg)

If we only consider the cases where the equipotential lines have S &gt; 0, the predictions by Eqs. (8) and (9) fail for p = 0 95, i.e., in the case of smooth surfaces, . where the values of the skewness of distant equipotential lines are S ∗ ∼ 4 as p → 1 (For example, see Figure 3 (a) for p &gt; 0 83). . In that case, dS/dφ grows slowly before the inflexion point mainly because of the variations of fluctuations, which are associated to longer wavelengths (See Figure 1 (b)). For p = 0 95, the region before . the inflection point corresponds to -3 0 . /lessorsimilar log 10 ( f P ) /lessorsimilar -2 5. . However, Figure 6 shows that, in this region, the absolute value of the exponent of f + is greater than 1 / 2. This result explains the rapid disappearance of roughness with lower wavelengths (which do not corresponds to pronounced peaks). Then, we encounter another regime with the absolute value of the exponent of f + less than 1 / 2, which explains the limit dS/dφ → 0 in Figure 3 (a) for p = 0 95. .

## 5. Summary and Conclusions

In this paper, we study the behavior of the height distributions of the equipotential lines and focus the behavior of their skewness and kurtosis. We used rough profiles that were grown following a competitive model with KPZ components, which show opposite signals of λ . Our results show a non-trivial behavior of S and K as a function of φ . We conducted our investigation for profiles with a lateral extension L = 10 6 grown until a time t = 2 × 10 , which considerably reduces (but not eliminate) the presence 4 of finite time effects for any p = 0.83.

/negationslash

For all KPZ profiles with λ &gt; 0, an almost invariant behavior of S and K on φ was found. The same behavior occurs to families for which λ &lt; 0, where S and K for the height distribution of the equipotential lines have an invariant behavior with respect to φ , although they are different from the previous families with λ &gt; 0. For families with λ ≈ 0, our results show a third distinct behavior for S and K . These results suggest that the electric potential variation can identify the signal of λ ,

which provides the signal of the skewness of the rough conducting profile. However, the height distribution statistics of equipotential lines far from the conducting profile were shown as being notably different.

The particular behavior of S and K as a function of φ when p ≈ p ∗ motivated the comparison of these parameters with the corresponding ones that were obtained for substrates that are characterized by identical EW values of growth exponents. This comparison was conducted considering the equipotential lines of rough profiles: (i) which are defined by p = 0 83, (ii) for . the case of a Family model, and (iii) at the time in the competitive model (for p =0.80 and p =0.825) where the effective growth exponent nearly 1 / 4. The results show similar behaviors for p = 0 83 and the Family . model but clearly different results when the effective growth exponent is nearly 1 / 4. At these values of t , where thin-film technology arises, the misinterpretations of small differences in β t ( ) (if experimentally reliable) compared with β EW can significantly affect the design of the electrical devices.

Finally, our results can provide new insights for cases where the rough morphologies are an important tool to construct devices where electrical properties play a fundamental role. For example in the design of devices for the flexible stacked resistive random access memory (RRAM) applications the rough surface leads to a local electric field enhancement which accelerates the oxygen ion back-drift by the external positive bias (reset process) [27]. Radio-frequency absorption and electric field enhancements due to surface roughness are also subjects of theoretical research [28]. However, there are situations where high roughness is an important tool, such as when the deposited surface is combined with metallic and nonmetallic materials. In that case, understanding the evolution of the surface morphology and the electric potential distribution along an irregular shape during electrodeposition phenomena is greatly significant [29].

Moreover, in our viewpoint, the current study has called the attention to interesting characteristics that help one more systematically interpret of images from probe microscopies, such as the Electrostatic Force Microscopy. In that case, the tip must be maintained far enough from the rough surface to avoid unwanted short-range interactions (which may result from van der Waals forces) and close enough to not loose important morphologic information about the real profile.

## Acknowledgements

This work has the financial support of CNPq and FAPESB (Brazilian agencies). C. P. de Castro acknowledges specific support from Funda¸ c˜o CAPES. a

## References

- [1] Barabasi A L, Stanley H E 1995 Fractal Concepts in Surface Growth (Cambridge - Cambridge U. Press).
- [2] Krug J, 1997 Adv. Phys. 46 , 139.
- [3] Forbes R G, 2012 Nanotech. 23 , 095706.
- [4] Du J, Zhang Y, Deng S, Xu N, Xiao Z, She J, Wu Z and Cheng H, 2013 Carbon , 61 , 507.
- [5] de Assis T A, Benito R M, Losada J C, Andrade R F S, Miranda J G V, de Souza N C, de Castilho C M C, Mota F de b and Borondo F, 2013 J. Phys:. Cond. Matt. , 25 , 285106.
- [6] Cabrera H, Zanin D A, De Pietro L G, Michaels Th, Thalmann P, Ramsperger U, Vindigni A, Pescia D, Kyritsakis A, Xanthakis J P, Fuxiang Li and Abanov Ar, 2013 Phys. Rev. B , 87 , 115436.
- [7] Kyritsakis A, Xanthakis J P, Pescia D, 2014 Proc. R. Soc. A , 470 , 20130795.

Height distribution of equipotential lines in a region confined by a rough conducting boundary 16

- [8] Cajueiro D O, Sampaio V A de A, de Castilho C M C and Andrade R F S, 1999 J. Phys.: Condens. Matter , 18 , 4985.
- [9] Filho H de O D, de Castilho C M C, Miranda J G V and Andrade R F S, 2004 Physica A , 342 , 388.
- [10] de Assis T A, Mota F de B, Miranda J G V, Andrade R F S, Filho H de O D and de Castilho C M C, 2006 J. Phys.: Condens. Matter , 18 , 3393.
- [11] de Assis T A, Mota F de B, Miranda J G V, Andrade R F S and de Castilho C M C, 2007 J. Phys.: Condens. Matter , 19 , 476215.
- [12] Budaev V P and Yakovlev M, 2008 Plasma and Fusion Ressearch (R) , 3 , 1.
- [13] da Silva T J and Moreira J G, 2001 Phys. Rev. E , 63 , 041601.
- [14] da Silva T J and Moreira J G, 2002 Phys. Rev. E , 66 , 061604.
- [15] Family F, 1986 J. Phys. A , 19 , l441.
- [16] Kardar M, Parisi G and Zhang Y C, 1986 Phys. Rev. Lett. , 56 , 889.
- [17] Edwards S F and Wilkinson D R, 1982 Proc. R. Soc. London , 381 , 17.
- [18] Takeuchi K A, Sano M, Sasamoto T and Spohn H, 2011 Sci. Rep. , 1 , 34.
- [19] Yunker P J, Lohr M A, Still T, Borodin A, Durian D J and Yodh A G, 2013 Phys. Rev. Lett. , 110 , 035501.
- [20] Muraca D, Braunstein L A and Buceta R C, 2004 Phys. Rev. E , 69 , 065103(R).
- [21] Silveira F A and Aar˜o Reis F D A, 2012 a Phys. Rev. E , 85 , 011601.
- [22] Oliveira T J, Dechoum K, Redniz J A and Aar˜o Reis F D A, 2006 a Phys. Rev. E , 74 , 011604.
- [23] Sasamoto T and Spohn H, 2010 Phys. Rev. Lett. , 104 , 230602.
- [24] de Assis T A, de Castro C P, de Brito Mota F, de Castilho C M C and Andrate R F S, 2012 Phys. Rev. E , 86 051607.
- [25] Oliveira T J , 2013 Phys. Rev. E 87 , 034401.
- [26] Silveira F A and Aar˜o Reis F D A, 2012 a Phys. Rev. E , 85 , 011601.
- [27] Hu Young Jeong, Yong In Kim, Lee J Y and Sung-Yool Choi, 2010 Nanotech. , 21 , 115203.
- [28] Zhang P, Lau Y Y and Gilgenbach R M, 2009 J. Appl. Phys. 105 , 114908.
- [29] Gamburg Y D and Zangari G 2011 Theory and Practice of Metal Electrodeposition (Springer, New York, USA).
C. P. de Castro T. A. de Assis C. M. C. de Castilho R. F. S. Andrade
Height distribution of equipotential lines in a region confined by a rough conducting boundary
This work considers the behavior of the height distributions of the
equipotential lines in a region confined by two interfaces: a cathode with an
irregular interface and a distant flat anode. Both boundaries, which are
maintained at distinct and constant potential values, are assumed to be
conductors. The morphology of the cathode interface results from the deposit of
$2 \times 10^{4}$ monolayers that are produced using a single competitive
growth model based on the rules of the Restricted Solid on Solid and Ballistic
Deposition models, both of which belong to the Kadar-Parisi-Zhang (KPZ)
universality class. At each time step, these rules are selected with
probability $p$ and $q = 1 - p$. For several irregular profiles that depend on
$p$, a family of equipotential lines is evaluated. The lines are characterized
by the skewness and kurtosis of the height distribution. The results indicate
that the skewness of the equipotential line increases when they approach the
flat anode, and this increase has a non-trivial convergence to a delta
distribution that characterizes the equipotential line in a uniform electric
field. The morphology of the equipotential lines is discussed; the discussion
emphasizes their features for different ranges of $p$ that correspond to
positive, null and negative values of the coefficient of the non-linear term in
the KPZ equation.