<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1353658</article-id>
<article-id pub-id-type="doi">10.3389/fphy.2024.1353658</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Formation of field-induced breakdown precursors on metallic electrode surfaces</article-title>
<alt-title alt-title-type="left-running-head">Bagchi et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fphy.2024.1353658">10.3389/fphy.2024.1353658</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Bagchi</surname>
<given-names>Soumendu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<xref ref-type="fn" rid="fn1">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2601199/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Simakov</surname>
<given-names>Evgenya</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Perez</surname>
<given-names>Danny</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Theoretical Division</institution>, <institution>Los Alamos National Laboratory</institution>, <addr-line>Los Alamos</addr-line>, <addr-line>NM</addr-line>, <country>United States</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Accelerator Operations and Technology Division</institution>, <institution>Los Alamos National Laboratory</institution>, <addr-line>Los Alamos</addr-line>, <addr-line>NM</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1387884/overview">Andreas Kyritsakis</ext-link>, University of Tartu, Estonia</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1382147/overview">Kristinn Torfason</ext-link>, Reykjav&#xed;k University, Iceland</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2624739/overview">Chris Moore</ext-link>, Sandia National Laboratories (DOE), United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Soumendu Bagchi, <email>soumendubagchi@gmail.com</email>, <email>bagchis@ornl.gov</email>
</corresp>
<fn fn-type="present-address" id="fn1">
<label>
<sup>&#x2020;</sup>
</label>
<p>
<bold>Present address:</bold> Soumendu Bagchi, Center for Nanophase Materials Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, United States</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>11</day>
<month>03</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1353658</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>20</day>
<month>02</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Bagchi, Simakov and Perez.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Bagchi, Simakov and Perez</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Understanding the underlying factors responsible for higher-than-anticipated local field enhancements required to trigger vacuum breakdown on pristine metal surfaces is crucial for the development of devices capable of withstanding intense operational fields. In this study, we investigate the behavior of nominally flat copper electrode surfaces exposed to electric fields of hundreds of MV/m. Our novel approach considers curvature-driven diffusion processes to elucidate the formation of sharp breakdown precursors. To do so, we develop a mesoscale finite element model that accounts for driving forces arising from both electrostatic and surface-tension-induced contributions to the free energy. Our findings reveal a dual influence: surface tension tends to mitigate local curvature, while the electric field drives mass transport toward regions of high local field density. This phenomenon can trigger the growth of sharper protrusions, ultimately leading to a rapid enhancement of local fields and, consequently, to a runaway growth instability. We delineate supercritical and subcritical regimes across a range of initial surface roughness. Our numerical results are in qualitative agreement with experimentally reported data, indicating the potential practical relevance of field-driven diffusion in the formation of breakdown precursors.</p>
</abstract>
<kwd-group>
<kwd>breakdown</kwd>
<kwd>precursor formation</kwd>
<kwd>mesoscale surface diffusion</kwd>
<kwd>electrodiffusion</kwd>
<kwd>field induced behavior</kwd>
<kwd>field-enhancment mechanisms</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Low-Temperature Plasma Physics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Strong electric fields are often encountered on solid material surfaces of high voltage and field-emission devices used in a broad range of applications, including high gradient accelerators [<xref ref-type="bibr" rid="B1">1</xref>&#x2013;<xref ref-type="bibr" rid="B6">6</xref>]. One of the key challenges in maintaining and improving the functionality of such devices at high fields is related to the control of field-induced breakdown [<xref ref-type="bibr" rid="B7">7</xref>&#x2013;<xref ref-type="bibr" rid="B10">10</xref>] or vacuum arc generating events [<xref ref-type="bibr" rid="B11">11</xref>]. As breakdown is routinely observed at fields that are much lower than those required for strong emission at flat clean surface [<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B13">13</xref>], the formation of sharp surface features which can locally enhance the applied field by factors of 50&#x2013;100x is broadly assumed [<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B14">14</xref>]. Indeed, the development of breakdown from sharp emitter geometries is now well understood to follow from thermal runaway following from strong field emission currents [<xref ref-type="bibr" rid="B15">15</xref>]. The process of partial melting followed by plasma/vacuum arc formation have been described via a multiphysics and multiscale modeling approaches [<xref ref-type="bibr" rid="B16">16</xref>]. Overall, plasma-driven and ejection mechanisms [<xref ref-type="bibr" rid="B17">17</xref>, <xref ref-type="bibr" rid="B18">18</xref>] during localized high current-induced Joule heating explain formation of micro craters in cathodes during vacuum arcs, under nominal electric fields (e.g., in the range of 100&#xa0;s of MV/m). In these models, the presence of a sharp emitter tip has proved to be an essential trigger of breakdown [<xref ref-type="bibr" rid="B19">19</xref>, <xref ref-type="bibr" rid="B20">20</xref>].</p>
<p>Since the observed surface roughness of polished surfaces is generally not thought to be compatible with the field-enhancement factors required to trigger spontaneous breakdown, the formation of precursors through some surface deformation mechanisms have to be invoked. To date, unraveling the precise mechanisms governing the formation of such sharp field-enhancing precursors has remained elusive, these are often destroyed in the breakdown process Therefore, understanding the physical mechanisms through which the material surface couples and dynamically evolves under experimentally relevant field to form breakdown <italic>precursors</italic> is paramount to design ultra-high gradient systems including the next-generation of compact accelerators. Several hypotheses have been proposed to explain the formation of such precursors [<xref ref-type="bibr" rid="B15">15</xref>, <xref ref-type="bibr" rid="B21">21</xref>, <xref ref-type="bibr" rid="B22">22</xref>]. These hypotheses have been guided by the notion that real electrodes deviate significantly from ideal materials. Material imperfections, such as point defects, voids, dislocations, slip extrusions, preexisting surface features, are proposed to interact with electric fields, resulting in localized field-enhancement factors. For example, the role of dislocation mediated plasticity [<xref ref-type="bibr" rid="B21">21</xref>] has been explored to address the formation of possible stress-concentrators [<xref ref-type="bibr" rid="B23">23</xref>] through slips interacting with free metal surfaces. However, both fully-atomistic [<xref ref-type="bibr" rid="B24">24</xref>] and hybrid continuum-atomistic [<xref ref-type="bibr" rid="B21">21</xref>] approaches suggest that fields on the order tens of GV/m (i.e., 2 orders of magnitude larger than typical experimental breakdown fields) are required to trigger substantial defect activity, unless these processes are assisted by other driving forces (e.g., thermo-elastic stresses caused by Joule heating [<xref ref-type="bibr" rid="B25">25</xref>]). Effect of electric field gradients on diffusion rates of metal emitter surfaces has been atomistically probed [<xref ref-type="bibr" rid="B26">26</xref>]. However, the long-range, long-time consequences of such biased diffusion on breakdown precursor formation was not directly elucidated. As of now, these models remain difficult to directly validate experimentally, so the precise mechanisms leading to the formation of breakdown precursors at electric fields on the order of hundreds of MV/m remains incompletely understood.</p>
<p>In this paper, we introduce a simple surface evolution model driven by surface tension and applied electric fields. Our findings underscore the interplay between these forces: while surface tension tends to diminish local curvature, the electric field drives mass transport towards regions with heightened field strengths. This dynamic leads to the formation of sharper protrusions, further intensifying local fields and ultimately resulting in a tip-growth instability and potentially in breakdown. Our investigation reveals a crucial aspect: the critical fields for instability exhibit a rapid increase with the growing aspect ratio of the initial surface morphology. For instance, in the case of a typical copper surface featuring a micron-scale roughness [<xref ref-type="bibr" rid="B27">27</xref>], our numerical modeling suggests that diffusive surface evolution can indeed lead to the formation of breakdown precursors at experimentally-relevant fields without invoking the presence of very specific microstructural features. Our investigation also reveals a crucial aspect: the critical fields for instability exhibit a rapid decrease with the growing aspect ratio of the initial surface morphology, as well as a slow but consequential decrease with increasing perturbation wavelength.</p>
<p>The paper is organized as follows: first in <xref ref-type="sec" rid="s2">Section 2</xref> we describe the mesoscale finite element model of morphological evolution of surfaces under electric fields. We then analyze the subcritical and supercritical regimes and identify fields associated with unstable growth of surface perturbations into sharp protrusions in <xref ref-type="sec" rid="s3">Section 3</xref>. Finally, <xref ref-type="sec" rid="s4">Section 4</xref> puts these results in relation with previous experimental findings and discusses possible extensions of the model.</p>
</sec>
<sec id="s2">
<title>2 Mesoscale curvature-driven growth model</title>
<p>
<italic>Surface kinetics</italic>&#x2013;Our computational model follows conventional descent dynamics where the system evolves to minimize its free energy. Such models where first put forth by Herring [<xref ref-type="bibr" rid="B28">28</xref>, <xref ref-type="bibr" rid="B29">29</xref>] and Mullins [<xref ref-type="bibr" rid="B30">30</xref>], and have since proven successful at predicting several surface and interface evolution phenomena [<xref ref-type="bibr" rid="B31">31</xref>] e.g., grain boundary grooving [<xref ref-type="bibr" rid="B32">32</xref>], electromigration-driven void collapse [<xref ref-type="bibr" rid="B33">33</xref>] in electronic interconnects, diffusive surface crack nucleation [<xref ref-type="bibr" rid="B34">34</xref>] and instability in stressed solids [<xref ref-type="bibr" rid="B35">35</xref>] etc. In our setting, atomic migration is assumed to occur only at the free surface and to be governed by a driving force that is given by the surface-gradient of free energy density (<italic>per-unit volume</italic>, <italic>g</italic>) <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>g</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>g</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>Depending on its wavelength (<italic>&#x3bb;</italic>), a small-amplitude perturbation on the surface can either decay or grow based on competition between effective contribution to free-energy density from surface tension (<italic>W</italic>
<sub>
<italic>surf</italic>
</sub>&#x2013;which tends to stabilize flat surfaces) and any external driving mechanisms (<italic>W</italic>
<sub>
<italic>ext</italic>
</sub>&#x2014;potentially leading to growth of curvature). We further assume the flux of mass (<italic>j</italic>
<sub>
<italic>s</italic>
</sub>) to be proportional to the driving forces (<xref ref-type="fig" rid="F1">Figure 1</xref>), i.e.,<disp-formula id="e1">
<mml:math id="m2">
<mml:msub>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>M</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>g</mml:mi>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m3">
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>&#x3b4;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">atom</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>g</mml:mi>
</mml:math>
<label>(2)</label>
</disp-formula>
<disp-formula id="e3">
<mml:math id="m4">
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>&#x3b4;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">atom</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">surf</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ext</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>.</mml:mo>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Curvature-driven growth of electrode surfaces under fields: <bold>(A)</bold> A metal electrode (magenta) is exposed to an applied electric field imposed at some distance in the vacuum (blue). Surface tension and electrostatic energies drive surface diffusion fluxes which leads to the evolution of the surface. <bold>(B)</bold> Breakdown precursor formation at a nominally flat electro-polished copper surface perturbed at a wavelength (<italic>&#x3bb;</italic>) 4&#xa0;<italic>&#x3bc;</italic>m [<xref ref-type="bibr" rid="B27">27</xref>] with a low aspect ratio (<italic>h</italic>
<sub>0</sub>/<italic>&#x3bb;</italic> &#x223c; 0.1) under <italic>E</italic>
<sub>
<italic>applied</italic>
</sub> &#x3d; 775&#xa0;MV/m. <bold>(C&#x2013;E)</bold> illustrates the morphological evolution of a sinusoidal perturbation of aspect ratio 0.1 under various dynamical regimes. Under a subcritical field of 250&#xa0;MV/m, the perturbation decays back toward a flat surface <bold>(C)</bold> Under a critical field of 520&#xa0;MV/m, the surface perturbation remains stable. <bold>(D)</bold> Under a super-critical field (775&#xa0;MV/m), a sharp tip quickly grows <bold>(E)</bold>.</p>
</caption>
<graphic xlink:href="fphy-12-1353658-g001.tif"/>
</fig>
<p>Here, <italic>M</italic> is the mobility of surface atoms which is related to surface self-diffusivity (<inline-formula id="inf2">
<mml:math id="m5">
<mml:mi>D</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula> with <inline-formula id="inf3">
<mml:math id="m6">
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>3.61</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mspace width="0.22em"/>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mo>&#x30a;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mi>s</mml:mi>
</mml:math>
</inline-formula>) by <inline-formula id="inf4">
<mml:math id="m7">
<mml:mi>M</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mi>&#x3b4;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">atom</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>B</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:math>
</inline-formula>, where <italic>k</italic>
<sub>
<italic>B</italic>
</sub> is the Boltzmann constant, <italic>T</italic> the temperature, &#x3a9;<sub>
<italic>atom</italic>
</sub> the volume per atom, and <italic>&#x3b4;</italic> the thickness of the diffusive layer. As the current study is mainly focused on unravelling the early stages of precursor formation, the model is restricted to a constant temperature of 300&#xa0;K. The surface free-energy contribution <italic>W</italic>
<sub>
<italic>surf</italic>
</sub> is taken to be proportional to the surface curvature, i.e., <italic>W</italic>
<sub>
<italic>surf</italic>
</sub> &#x3d; <italic>&#x3b3;&#x3ba;</italic>, where <italic>&#x3ba;</italic> is the surface curvature<disp-formula id="e4">
<mml:math id="m8">
<mml:mi>&#x3ba;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:math>
<label>(4)</label>
</disp-formula>
<italic>&#x3b3;</italic> is the surface tension, and <italic>h</italic>(<italic>x</italic>) is the surface height at position <italic>x</italic>. In the following, <italic>W</italic>
<sub>
<italic>ext</italic>
</sub> accounts for the coupling of the system with an externally applied electric field (hence, <italic>W</italic>
<sub>
<italic>ext</italic>
</sub> &#x3d; <italic>W</italic>
<sub>
<italic>elec</italic>
</sub>).</p>
<p>Invoking conservation of mass, a continuity equation for the flux leads to the evolution of velocity (<italic>v</italic>
<sub>
<italic>n</italic>
</sub>) normal to the surface,<disp-formula id="e5">
<mml:math id="m9">
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mo>,</mml:mo>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>Summarizing from Eqs <xref ref-type="disp-formula" rid="e1">1</xref>&#x2013;<xref ref-type="disp-formula" rid="e5">5</xref>, the total velocity <inline-formula id="inf5">
<mml:math id="m10">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">total</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> can be decomposed into the surface-tension (<italic>v</italic>
<sub>
<italic>surf</italic>
</sub>) and electric-field (<italic>v</italic>
<sub>
<italic>elec</italic>
</sub>) contribution as<disp-formula id="equ1">
<mml:math id="m11">
<mml:mtable class="align-star" columnalign="left">
<mml:mtr>
<mml:mtd columnalign="right">
<mml:msubsup>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">total</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>M</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">surf</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">elec</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mtext>&#x2009;for&#x2009;homogeneous&#x2009;</mml:mtext>
<mml:mi>M</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mo>&#x3d;</mml:mo>
<mml:mi>M</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">surf</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">elec</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">surf</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">elec</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
</p>
<p>For the vertical (<italic>y</italic>) component of the total velocity,<disp-formula id="e6">
<mml:math id="m12">
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>.</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mo>,</mml:mo>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>Where <italic>v</italic>
<sub>
<italic>y</italic>
</sub> is the local vertical velocity of the surface (Eq. <xref ref-type="disp-formula" rid="e7">7</xref>) such that for a surface profile <italic>h</italic>(<italic>x</italic>),<disp-formula id="e7">
<mml:math id="m13">
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>h</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x2202;</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>
<italic>Electrostatics</italic>&#x2013;Consider a domain as shown in <xref ref-type="fig" rid="F1">Figure 1A</xref> which has a vacuum (&#x3a9;<sub>
<italic>v</italic>
</sub>) and an electrode material component (&#x3a9;<sub>
<italic>m</italic>
</sub>) which is assumed to be a perfect conductor. We obtain the electrostatic potential <italic>&#x3d5;</italic> and the electric field <italic>E</italic> &#x3d; &#x2212;&#x2207;<italic>&#x3d5;</italic>) through the solution of the Poisson equation:<disp-formula id="e8">
<mml:math id="m14">
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">&#x394;</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
<label>(8)</label>
</disp-formula>
<disp-formula id="e9">
<mml:math id="m15">
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:mi>&#x2202;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2261;</mml:mo>
<mml:mi>&#x2202;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">bot</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2261;</mml:mo>
<mml:mi>&#x2202;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2229;</mml:mo>
<mml:mi>&#x2202;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2261;</mml:mo>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:math>
<label>(9)</label>
</disp-formula>
<disp-formula id="e10">
<mml:math id="m16">
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x2207;</mml:mi>
<mml:mi>&#x3d5;</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">applied</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mspace width="1em"/>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mspace width="1em"/>
<mml:mi>&#x2202;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2261;</mml:mo>
<mml:mi>&#x2202;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">top</mml:mi>
</mml:mrow>
</mml:msup>
</mml:math>
<label>(10)</label>
</disp-formula>with <italic>E</italic>
<sub>
<italic>applied</italic>
</sub> (Eq. <xref ref-type="disp-formula" rid="e10">10</xref>) the macroscopically-applied field. The electrostatic contribution to the free energy density at the surface is then obtained as <inline-formula id="inf6">
<mml:math id="m17">
<mml:msub>
<mml:mrow>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">elec</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi>&#x3f5;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula>, where <italic>&#x3f5;</italic> is the vacuum permittivity, and <italic>E</italic>
<sub>
<italic>n</italic>
</sub> as shown in <xref ref-type="fig" rid="F1">Figure 1A</xref> is normal field acting on the surface &#x393;. Periodic boundary conditions are applied in the lateral (<italic>x</italic>) direction.<italic>Finite element formulation</italic>&#x2013;We solve the governing equations of surface kinetics (Eqs <xref ref-type="disp-formula" rid="e1">1</xref>, <xref ref-type="disp-formula" rid="e6">6</xref>) using a mixed finite element formulation in two dimension. For test functions <italic>u</italic> &#x2208; <italic>U</italic> and <italic>w</italic> &#x2208; <italic>W</italic> defined over the surface manifold (&#x393;), we obtain a variational weak form by integrating by parts:<disp-formula id="e11">
<mml:math id="m18">
<mml:msub>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>M</mml:mi>
<mml:mi>g</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>s</mml:mi>
<mml:mspace width="1em"/>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>&#x2207;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>w</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>x</mml:mi>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>The discretized solution space of Eq. <xref ref-type="disp-formula" rid="e11">11</xref> is a mixed space <inline-formula id="inf7">
<mml:math id="m19">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>U</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>. For both surface kinetics and electrostatics, we use Lagrangian elements of order 1. We currently model the surface as a 1D manifold (&#x393;<sup>
<italic>h</italic>
</sup>) mesh. The 2D vacuum domain <inline-formula id="inf8">
<mml:math id="m20">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> for electrostatics is discretized with triangular elements using a finer mesh closer to the vacuum/material interface &#x393;, matching the nodes of the surface mesh. We use an explicit forward Euler time-integration scheme (<italic>h</italic> (<italic>x</italic>, <italic>t</italic> &#x2b; &#x394;<italic>t</italic>) &#x3d; <italic>h</italic> (<italic>x</italic>, <italic>t</italic>) &#x2b; &#x394;<italic>tv</italic>
<sub>
<italic>y</italic>
</sub>). We have used mesh sizes ranging from <italic>&#x3bb;</italic>/100 (aspect ratio &#x3d; 0.1) to <italic>&#x3bb;</italic>/200 (aspect ratio &#x3d; 2) which leads to 100 (200) to 500 (1,000) mesh-points along x (y) direction (i.e., 20&#xa0;K to 500&#xa0;k mesh-points in the electrostatics domain). Mesh-sizes are experimented with and confirmed that the analytical (i.e., linear stability predicted (sub/super) critical field regimes matched with our numerical results for a lower aspect ratio periodic sinusoidal surface. Additionally, for the Poisson solver, the fact that local e-field enhancement at the tip is almost linearly increasing with aspect ratio (see <xref ref-type="sec" rid="s3">Section 3</xref>), a crude measure of local curvature, is assuring that results are not significantly affected by mesh size. The code adaptively reduces the timestep to maintain a predefine tolerance value for the maximum allowed nodal displacement. For most of the reported cases, it is conservatively capped at 1e-7 <italic>&#x3bb;</italic>, starting with an initial timestep (&#x394;<italic>t</italic>) of 1e-9 <italic>&#x3bb;</italic>
<sup>4</sup> in units of <italic>&#x3c4;</italic>.</p>
<p>In a staggered solution approach, the electrostatics (Eq. <xref ref-type="disp-formula" rid="e8">8</xref>) is first solved to evaluate <italic>W</italic>
<sub>
<italic>elec</italic>
</sub> on the coinciding nodes of the meshes of the domains &#x3a9;<sub>
<italic>v</italic>
</sub> and &#x393;, the surface profile <italic>h</italic> (<italic>x</italic>, <italic>t</italic>) is then evolved using a small timesteps &#x394;<italic>t</italic>. As the electrostatic solution changes slowly on the timescale &#x394;<italic>t</italic> required for numerical stability of the integration of the surface kinetic equations, the electrostatic solution is updated less frequently. A frequency ratio of 100&#x2013;200 was found to give good results.</p>
<p>In the resulting model, termed <italic>SurFE-XD Surface curvature-driven Finite Elements model for Diffusion under eXtreme conditions</italic>, the variational forms are specified in the FEniCS workflow [<xref ref-type="bibr" rid="B36">36</xref>] in the high-level Python-based Unified Form Language (UFL). These forms are then automatically compiled [<xref ref-type="bibr" rid="B37">37</xref>] and executed through high-performance computational kernels using the finite element library DOLFIN [<xref ref-type="bibr" rid="B38">38</xref>].</p>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<p>First, we demonstrate and verify our simulation results for a periodic sinusoidal profile (<inline-formula id="inf9">
<mml:math id="m21">
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> with a low aspect ratio (<italic>h</italic>
<sub>0</sub>/<italic>&#x3bb;</italic> &#x223c; 0.1). Based on experimental measurements of electropolished copper surfaces [<xref ref-type="bibr" rid="B27">27</xref>], we choose a wavelength of 4&#xa0;<italic>&#x3bc;</italic>m to present the results in terms of absolute physical units. We note that our continuum FEM framework only depends on relative geometry and is therefore scale-invariant, and the results can be rescaled for different wavelength, as we show in the following.</p>
<p>In <xref ref-type="fig" rid="F1">Figures 1A&#x2013;D</xref>, a critical field of 520&#xa0;MV/m is identified below which surface tension dominates, leading to the gradual decrease in amplitude of the initial surface profile. Above the critical field, the electric field drives mass transport towards the tip, leading to the emergence of a sharp surface feature that has the characteristics of an efficient breakdown precursor (<xref ref-type="fig" rid="F1">Figure 1E</xref>).</p>
<p>At low aspect ratio, one can compare the numerical results with growth/decay rates obtained from linear-stability analysis with an initial profile, say e.g., <inline-formula id="inf10">
<mml:math id="m22">
<mml:msub>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:math>
</inline-formula>. The time evolution of the surface amplitude is then given by <italic>A</italic>(<italic>t</italic>) &#x3d; <italic>A</italic>
<sub>0</sub>
<italic>e</italic>
<sup>
<italic>&#x3b1;t</italic>
</sup>, where the rate <italic>&#x3b1;</italic> determines growth or decay under a given applied E-field (<italic>E</italic>
<sub>
<italic>app</italic>
</sub>) [<xref ref-type="bibr" rid="B35">35</xref>]. One obtains <inline-formula id="inf11">
<mml:math id="m23">
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x221d;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3f5;</mml:mi>
<mml:msubsup>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">app</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> leading to decay (<italic>&#x3b1;</italic> &#x3c; 0) for <inline-formula id="inf12">
<mml:math id="m24">
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">app</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>&#x3b3;</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:math>
</inline-formula> and growth when <inline-formula id="inf13">
<mml:math id="m25">
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">app</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>&#x3b3;</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>&#x3f5;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:math>
</inline-formula>. Therefore, for a copper surface with <italic>&#x3bb;</italic> &#x3d; 4<italic>&#x3bc;</italic>m, <italic>&#x3b3;</italic> &#x3d; 1.5&#xa0;J/<italic>m</italic>
<sup>2</sup> [<xref ref-type="bibr" rid="B39">39</xref>] and <italic>&#x3f5;</italic> &#x3d; 8.85 &#xd7; 10<sup>&#x2212;12</sup>&#xa0;F/m for vacuum, the critical field is predicted to be around 520&#xa0;MV/m.</p>
<p>Our simulation results are in close agreement with the predictions of this linear stability analysis. One of the key limitations of this approximation is that is it only valid in the limit of small aspect ratios, i.e., <italic>A</italic>
<sub>0</sub>/<italic>&#x3bb;</italic> &#x226A; 1. Outside of this regime, explicit numerical solutions like that provided by (<italic>SurFE-XD</italic>) is essential. It is important to note that the critical field is inversely proportional to the square root of the wavelength <italic>&#x3bb;</italic>, and so could be higher or lower depending on the physical scale of the surface roughness. In particular, longer wavelengths lead to a lowering of the critical field. E.g., for copper, a wavelength of 100&#xa0;<italic>&#x3bc;</italic>m would correspond to a critical field to 100&#xa0;MV/m. The specific values cited should therefore be taken as indicative of the order of magnitude of the critical fields, but can be expected to vary with the specifics of each surface finish.</p>
<p>In order to understand the effect of strong electric fields on the growth dynamics of nominally flat (<italic>h</italic>
<sub>0</sub>/<italic>&#x3bb;</italic> &#x223c; 0.1) copper surfaces, we analyze their evolution under a near-critical (775&#xa0;MV/m) and stronger supercritcal (1.5&#xa0;GV/m) regimes. <xref ref-type="fig" rid="F2">Figure 2A</xref> shows that, at slightly super-critical fields, precursor formation occurs by long-range mass transport towards the hillock region. As the field are further increased into the super-critical regime, the growth becomes much more localized, leading to the formation of sharp emitter (<xref ref-type="fig" rid="F2">Figure 2B</xref>). The growth speed of the initial perturbation is also observed to be strongly dependent on the applied field. At 775&#xa0;MV/m, the timescale required for growth up to an amplitude of 0.5&#xa0;<italic>&#x3bc;</italic>m in the orders of <inline-formula id="inf14">
<mml:math id="m26">
<mml:mo>&#x223c;</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> in comparison with a timescale of about <inline-formula id="inf15">
<mml:math id="m27">
<mml:mo>&#x223c;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>&#x3c4;</mml:mi>
</mml:math>
</inline-formula> at 1.5&#xa0;GV/m, where <italic>&#x3c4;</italic> is the normalized time discussed in more detail in <xref ref-type="sec" rid="s4">Section 4</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Breakdown precursor formation under strong e-fields: Surface evolution at 775&#xa0;MV/m <bold>(A)</bold> and 1.5&#xa0;GV/m <bold>(B)</bold>. Field-enhancement at a particular time instance, e.g., <italic>t</italic> &#x3d; 10 &#xd7; 10<sup>&#x2212;4</sup>
<italic>&#x3c4;</italic> (cyan surface profiles in <bold>(A,B)</bold>) is compared through the corresponding spatial distribution of the local electric field magnitude under an applied field of 1&#xa0;GV/m in <bold>(C)</bold> (case a.) and <bold>(D)</bold> (case b).</p>
</caption>
<graphic xlink:href="fphy-12-1353658-g002.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F2">Figures 2C, D</xref> shows that the electric field localizes at the hillock regions while it is being partially shielded at the trough or valley regions. In the near critical regime (e.g., 775&#xa0;MV/m), the contribution of surface tension towards mass transport is of similar order with that of the electric field. Hence, while local enhancement at the hillock region drives the curvature to grow, active mass transport from the valley towards the hillock lead to long-range mass transport. On the other hand, at 1.5&#xa0;GV/m, the electrostatic driving force completely overwhelms surface tension in the hillock, leading to a localized and rapid growth of a very sharp tip. The rapid nonlinear growth towards instability is accelerated by local field enhancement in highly-curved regions.</p>
<p>To further probe the effect of local geometry on field induced breakdown precursor formation, we model the evolution of isolated Gaussian-shaped features (<inline-formula id="inf16">
<mml:math id="m28">
<mml:mi>h</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mspace width="1em"/>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>h</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula>, i.e., mean &#x3d; 0 and standard deviation &#x3d; <italic>&#x3bb;</italic>/2<italic>&#x3c0;</italic>) with aspect ratios (<italic>h</italic>
<sub>0</sub>/<italic>&#x3bb;</italic>) ranging from very low (0.1) to high (2.0). For a similar wavelength of 4&#xa0;<italic>&#x3bc;</italic>m, <xref ref-type="fig" rid="F3">Figure 3A</xref> depicts the sharp magnification of curvature at the tip of the Gaussian. This leads to a drastic change in the surface evolution velocity, leading to a surface evolution that is more localized, larger in magnitude, and faster, as the aspect ratio becomes larger.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Effect of the initial aspect ratio on the relative contributions of surface tension and electric fields. <bold>(A)</bold> Curvature (<italic>&#x3ba;</italic>(<italic>x</italic>)) of Gaussian surface profiles with aspect ratio 0.5 and 1.5. Middle panel <bold>(B,C)</bold> with aspect ratio 0.5 and lower panel <bold>(D,E)</bold> with 1.5 show spatial variation of individual and total velocity contributions from electric field and surface tension driven forces under fields 250&#xa0;MV/m <bold>(B,D)</bold> and 775&#xa0;MV/m at <italic>t</italic> &#x3d; 0 <bold>(C,E)</bold>.</p>
</caption>
<graphic xlink:href="fphy-12-1353658-g003.tif"/>
</fig>
<p>For a fixed width of 4<italic>&#x3bc;</italic>m, the individual contributions of surface tension (<italic>v</italic>
<sub>
<italic>&#x3b3;</italic>
</sub>) and electrostatics (<italic>v</italic>
<sub>
<italic>elec</italic>
</sub>) again act in opposite directions, as shown in <xref ref-type="fig" rid="F3">Figures 3B&#x2013;E</xref>. Under a field of 250&#xa0;MV/m, for a surface feature with a modest aspect-ratio (0.5), (<xref ref-type="fig" rid="F3">Figure 3B</xref>) surface tension slightly dominates over the electrostatic driving force at the tip location (<italic>x</italic> &#x3d; 0), leading to the decay of the initial protrusion. This trend is completely reversed at a field of 775&#xa0;MV/m for the same aspect ratio 0.5 (<xref ref-type="fig" rid="F3">Figure 3C</xref>); here the electric-field-induced velocity drives the growth of the tip resulting in the rapid growth of sharp feature that could act as a breakdown precursor.</p>
<p>At even higher aspect ratios (1.5, c.f. <xref ref-type="fig" rid="F3">Figure 3D</xref>) the local field enhancement at highly curved regions magnifies the electro-static driving forces for diffusion, inducing tip growth at fields that were previously sub-critical for smaller protrusions. E.g., at field of 250&#xa0;MV/m, the electrostatics contributions slightly exceed the restoring tendency of the surface tension, resulting in a net positive velocity at tip. This is in contrast with the much higher linear stability estimation of 520&#xa0;MV/m reported above. We note that two effects are at play: first, field enhancement is more efficient at isolated protrusions due to the lack of interactions between tips (however note that in our simulations, periodic boundary conditions remain in the lateral direction so tip-tip interactions are not totally absent, but they are smaller in magnitude); second, non-linear interactions between modes invalidate the predictions obtained in the linear regime.</p>
<p>The respective contributions of surface tension and electric fields is further analyzed in <xref ref-type="fig" rid="F4">Figure 4</xref>, which presents the ratio of the different velocity contributions (<inline-formula id="inf17">
<mml:math id="m29">
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">elec</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:math>
</inline-formula> at the tip of the Gaussian feature. Here, the sub-critical and super-critical regimes correspond to <italic>v</italic>
<sub>
<italic>elec</italic>
</sub>/<italic>v</italic>
<sub>
<italic>&#x3b3;</italic>
</sub> &#x3c; 1 and <italic>v</italic>
<sub>
<italic>elec</italic>
</sub>/<italic>v</italic>
<sub>
<italic>&#x3b3;</italic>
</sub> &#x3e; 1, respectively. The results show that the velocity ratio grows super-linearly as the aspect-ratio (<italic>h</italic>
<sub>0</sub>/<italic>&#x3bb;</italic>) and that the effect becomes even more prominent when the applied fields become stronger. At a low aspect-ratio of 0.1, we observe <italic>v</italic>
<sub>
<italic>elec</italic>
</sub>/<italic>v</italic>
<sub>
<italic>&#x3b3;</italic>
</sub> &#x223c; 1 at 546&#xa0;MV/m, consistent with linear stability predictions. However, with larger aspect-ratio (e.g., 1.5), the velocity ratio is around 5 at this same field. The velocity comparison due to fields of 250&#xa0;MV/m and 775&#xa0;MV/m in <xref ref-type="fig" rid="F3">Figure 3</xref> show almost an order of magnitude difference when aspect-ratio is even larger (1.5), demonstrating that the balance quickly shifts in favor of electro-static contributions as the aspect ratio increases, leading to tip-growth instabilities at progressively lower fields. As can be seen in <xref ref-type="fig" rid="F5">Figure 5B</xref>, the critical field for the tip to grow can reduce from 520&#xa0;MV/m to 250&#xa0;MV/m as the aspect-ratio increases from nominal to a higher value of 2.0. This is in line with the 4-fold increase in field-enhancement near tip (<xref ref-type="fig" rid="F5">Figure 5B</xref>).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Ratio of the vertical surface velocity contributions stemming from electrostatics and surface tension at the location of the tip as a function of aspect ratio for a Gaussian protrusion of 4&#xa0;<italic>&#x3bc;</italic>m in width.</p>
</caption>
<graphic xlink:href="fphy-12-1353658-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Time evolution of precursor growth for identifying critical electric fields: <bold>(A)</bold> For an aspect-ratio of 1.0, subcritical (260&#xa0;MV/m), near-critical (312&#xa0;MV/m) and supercritical (364&#xa0;MV/m) regimes of tip-displacement evolution. Inset Sharpening of tip (magenta fill) at supercrtical (364&#xa0;MV/m) compared to initial tip geometry in blue dashed line. <bold>(B)</bold> Critical e-fields (blue) decrease with increasing aspect ratio while field enhancement at the tip (green) linearly increases.</p>
</caption>
<graphic xlink:href="fphy-12-1353658-g005.tif"/>
</fig>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>A key factor that should be quantified in order to assess the efficiency of this mechanism is the timescale over which the growth of such surface features can be expected to occur. The time-unit in the simulation <italic>&#x3c4;</italic> is estimated as:<disp-formula id="e12">
<mml:math id="m30">
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3b3;</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>Considering Cu parameters (c.f. <xref ref-type="sec" rid="s2">Section 2</xref>) <inline-formula id="inf18">
<mml:math id="m31">
<mml:mi>D</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula> with <inline-formula id="inf19">
<mml:math id="m32">
<mml:msub>
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>3.61</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>5</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>1</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>12</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mspace width="0.22em"/>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mo>&#x30a;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>/</mml:mo>
<mml:mi>s</mml:mi>
</mml:math>
</inline-formula>; <inline-formula id="inf20">
<mml:math id="m33">
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>3.615</mml:mn>
<mml:mspace width="0.22em"/>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mo>&#x30a;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf21">
<mml:math id="m34">
<mml:mi mathvariant="normal">&#x3a9;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>11.81</mml:mn>
<mml:mspace width="0.22em"/>
<mml:msup>
<mml:mrow>
<mml:mover accent="true">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mo>&#x30a;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula> [<xref ref-type="bibr" rid="B40">40</xref>]), (Eq. <xref ref-type="disp-formula" rid="e12">12</xref>) can range from 475&#xa0;s to 13&#xa0;days. Such a wide range in physical timescales is due to the exponential dependence of surface adatom diffusivity on energy barrier (<italic>E</italic>
<sub>
<italic>b</italic>
</sub>) which is changes considerably with change in local environment; for instance, <italic>E</italic>
<sub>
<italic>b</italic>
</sub> could range from 0.1&#x2013;0.15&#xa0;eV for (111), 0.25&#x2013;0.30&#xa0;eV for (110), and 0.38&#x2013;0.69&#xa0;eV for (100) [<xref ref-type="bibr" rid="B41">41</xref>&#x2013;<xref ref-type="bibr" rid="B44">44</xref>]. It is difficult to lump complex atomistic details on a single effective diffusive characteristic without more sophisticated multiscale investigation which is beyond the scope of the current study. In comparison, the geometric factors <italic>&#x3b4;</italic> (a typical interatomic distance) and &#x3a9; (the volume per atom in the bulk) are easier to constrain, and only lead to linear variations in the estimated physical timescales. We also note that the characteristic timescale <italic>&#x3c4;</italic> is also exponentially dependent on temperature, which can dramatically affect the correspondence between reduced and physical time. Importantly however, this conversion is independent of the applied field and of the initial surface characteristics, so the relative timescales observed in different simulations can be compared to one another in both reduced and physical units. While quantitative predictions of time to breakdown remain difficult, the order of magnitude predicted by this model are compatible with experimental timescales.</p>
<p>Further, the simulations have shown that breakdown precursors can be expected to spontaneously form through surface transport at fields in low hundreds of MV/m, depending on the initial roughness of the surface. This mechanism can therefore potentially contribute to the spontaneous formation of breakdown-inducing features in application-relevant conditions. Indeed, the conditions investigated here are consistent with experiments carried on for breakdown with DC pulses with Copper electrodes [<xref ref-type="bibr" rid="B12">12</xref>]. These experiments have shown that i) the spatial location of the breakdown events is highly correlated, with a high fraction of craters overlapping, and ii) the breakdown events occurred preferentially close to the edge of the electrode. The first point is consistent with the observation that breakdown events leave craters and debris behind [<xref ref-type="bibr" rid="B12">12</xref>, <xref ref-type="bibr" rid="B45">45</xref>]. As shown above, these surface features can couple efficiently with the electric field, potentially leading to fast tip-growth instabilities and to further breakdown events occurring in close proximity and in rapid succession, consistent with experimental results [<xref ref-type="bibr" rid="B45">45</xref>, <xref ref-type="bibr" rid="B46">46</xref>]. Second, the edge of the cathode will be exposed to additional macroscopic field gradients which can also drive surface diffusion. This effect has been observed in nanoscale experiments on gold surfaces where electric field gradients at the edge of the electrode led to enhanced formation of ridge and tip-like features [<xref ref-type="bibr" rid="B47">47</xref>]. Both these observations are consistent with the hypothesis that breakdown was there mediated by surface diffusion driven by electric fields.</p>
<p>Finally, our simulations show that the mechanism investigated here can lead to very-fast runaway instabilities. Indeed, <xref ref-type="fig" rid="F5">Figure 5A</xref> shows that tip growth proceeds in two stages: an initial stage where the tip grows roughly linearly in time, and a second runaway stage where exponential growth is observed. This suggests that directly observing these sharp precursors in experiments can be expected to be very challenging, as the onset of runaway growth occurs at relatively modest aspect ratios and then unfolds very rapidly, leading to the destruction of the tip during the breakdown process. During short-time scanning probe investigations of surface roughness of electrodes in vacuum under moderate electric fields, sharp features leading to high field-enhancement factors were rarely observed [<xref ref-type="bibr" rid="B48">48</xref>]. This is consistent with a mechanism where sharp breakdown precursors form only shortly before breakdown occurs.</p>
<p>While field-driven surface diffusion appears to be a viable mechanism to explain the formation of breakdown precursors, other phenomena could contribute. For example, theoretical treatments typically assume clean surfaces free of contamination. Any contaminants could potentially affect surface evolution in complex fashion. Second, many high-field applications occur in high-frequency AC settings [<xref ref-type="bibr" rid="B2">2</xref>, <xref ref-type="bibr" rid="B13">13</xref>]. In this case, losses through Joule heating induce thermo-elastic stresses. These introduce an additional driving force for surface deformation [<xref ref-type="bibr" rid="B34">34</xref>] and can also lead to the formation of surface roughness through thermal fatigue [<xref ref-type="bibr" rid="B25">25</xref>]. Other mechanisms could therefore &#x201c;seed&#x201d; the tip-growth process by introducing surface features, which would then grow and sharpen through the diffusive mechanism discussed here. The contributions of thermo-elastic stresses on surface evolution are considered in a follow-on paper, where significant reductions in critical fields are observed for temperature rises on the order of tens of degrees [<xref ref-type="bibr" rid="B49">49</xref>]. We also ote that the present study mostly focuses with micron-scale roughnesses in metal electrodes as observed in electropolished surfaces. A different polishing approach (e.g., &#x201c;mirror-polish&#x201d;) might lead to submicron wavelengths. Furthermore, postmortem investigations of high-gradient breakdown surfaces in rf cavities revealed that pristine roughness conditions might be greatly modified by the presence of micron-scale debris or adsorbates, which could further affect the coupling with the electric field [<xref ref-type="bibr" rid="B50">50</xref>].</p>
<p>Finally, our study focuses only on early-stages of precursor formation, before runaway mechanisms that dissipate large amounts of energy are activated. Once they do so, several other mechanisms related to plasma formation e.g., ion-desorption [<xref ref-type="bibr" rid="B17">17</xref>], droplet ejection and emission-current induced heat transport [<xref ref-type="bibr" rid="B18">18</xref>] will become relevant in the intermediate stages of breakdown once field-induced surface instability has triggered localized emission [<xref ref-type="bibr" rid="B51">51</xref>, <xref ref-type="bibr" rid="B52">52</xref>]/melting [<xref ref-type="bibr" rid="B16">16</xref>]. We have also neglected any effects on mass loss due to vaporization, which are not expected to be significant in the early stages of precursor formation, but would become significant once the runaway process begins.</p>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>In this study, the evolution of metallic surfaces exposed to strong electric fields was modeled using a surface diffusion model driven by surface tension and electro-static driving forces. The numerical results are consistent with the prediction of linear stability analysis in the regime small aspect-ratio perturbations on the surface. Non-linearities in the growth process start to emerge with high aspect-ratio surface features as well as under stronger electric fields. Our results show that exposure to electric fields in the low hundreds of MV/m can trigger the spontaneous formation of sharp surface features that can act as precursors for breakdown. These fields are of the same order of magnitude as those experimentally observed to lead to breakdown, suggesting that field-driven transport is a plausible mechanism for precursor formation, acting either alone, or in concert with other surface deformation mechanisms.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary materials, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>SB: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing&#x2013;original draft, Writing&#x2013;review and editing. ES: Funding acquisition, Project administration, Resources, Supervision, Validation, Writing&#x2013;original draft, Writing&#x2013;review and editing. DP: Conceptualization, Investigation, Project administration, Resources, Supervision, Validation, Writing&#x2013;original draft, Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Laboratory Directed Research and Development program of Los Alamos National Laboratory under project number 20230011DR. Los Alamos National Laboratory is operated by Triad National Security, LLC, for the National Nuclear Security Administration of U.S. Department of Energy (Contract No. 89233218CNA000001).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sicking</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Str&#xf6;m</surname>
<given-names>R</given-names>
</name>
</person-group>. <article-title>From precision physics to the energy frontier with the Compact Linear Collider</article-title>. <source>Nat Phys</source> (<year>2020</year>) <volume>16</volume>:<fpage>386</fpage>&#x2013;<lpage>92</lpage>. <pub-id pub-id-type="doi">10.1038/s41567-020-0834-8</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Simakov</surname>
<given-names>EI</given-names>
</name>
<name>
<surname>Dolgashev</surname>
<given-names>VA</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>SG</given-names>
</name>
</person-group>. <article-title>Advances in high gradient normal conducting accelerator structures</article-title>. <source>Nucl Instr Methods Phys Res Section A: Acc Spectrometers, Detectors Associated Equipment</source> (<year>2018</year>) <volume>907</volume>:<fpage>221</fpage>&#x2013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.1016/j.nima.2018.02.085</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schneider</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Dolgashev</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Lewellen</surname>
<given-names>JW</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>SG</given-names>
</name>
<name>
<surname>Nanni</surname>
<given-names>EA</given-names>
</name>
<name>
<surname>Zuboraj</surname>
<given-names>M</given-names>
</name>
<etal/>
</person-group> <article-title>High gradient off-axis coupled C-band Cu and CuAg accelerating structures</article-title>. <source>Appl Phys Lett</source> (<year>2022</year>) <volume>121</volume>. <pub-id pub-id-type="doi">10.1063/5.0132706</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nanni</surname>
<given-names>EA</given-names>
</name>
<name>
<surname>Breidenbach</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Vernieri</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F</given-names>
</name>
<name>
<surname>White</surname>
<given-names>G</given-names>
</name>
<etal/>
</person-group> <article-title>Status and future plans for C<sup>3</sup> R&#x26;amp;D</article-title>. <source>J Instrumentation</source> (<year>2023</year>) <volume>18</volume>:<fpage>P09040</fpage>. <pub-id pub-id-type="doi">10.1088/1748-0221/18/09/p09040</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schulte</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Johnstone</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Boucher</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Esarey</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Geddes</surname>
<given-names>CG</given-names>
</name>
<name>
<surname>Kravchenko</surname>
<given-names>M</given-names>
</name>
<etal/>
</person-group> <article-title>Transformative technology for FLASH radiation therapy</article-title>. <source>Appl Sci</source> (<year>2023</year>) <volume>13</volume>:<fpage>5021</fpage>. <pub-id pub-id-type="doi">10.3390/app13085021</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kutsaev</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Novel technologies for compact electron linear accelerators (review)</article-title>. <source>Instr Exp Tech</source> (<year>2021</year>) <volume>64</volume>:<fpage>641</fpage>&#x2013;<lpage>56</lpage>. <pub-id pub-id-type="doi">10.1134/s0020441221050079</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dolgashev</surname>
<given-names>VA</given-names>
</name>
<name>
<surname>Higashi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Higo</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>SG</given-names>
</name>
</person-group>. <conf-name>Proc EPAC</conf-name>, <conf-loc>Genoa, Italy</conf-loc>, <conf-date>June, 2008</conf-date>, (<year>2008</year>) <fpage>742</fpage>. </citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dolgashev</surname>
<given-names>VA</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>S. G.</given-names>
</name>
<name>
<surname>Higashi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Tsukuba</surname>
<given-names>KEK</given-names>
</name>
</person-group>. <conf-name>Proc LINAC</conf-name> <conf-loc>Tsukuba, Japan</conf-loc>, <conf-date>September, 2010</conf-date>(<year>2010</year>) <volume>1043</volume>.</citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Dolgashev</surname>
<given-names>VA</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>SG</given-names>
</name>
<name>
<surname>Yeremian</surname>
<given-names>AD</given-names>
</name>
<name>
<surname>Higashi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Spataro</surname>
<given-names>B</given-names>
</name>
</person-group>. <source>Proc IPAC</source>, <conf-loc>Kyoto, Japan</conf-loc>, <conf-date>May, 2010</conf-date>, (<year>2010</year>) <fpage>3810</fpage>.</citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grudiev</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Calatroni</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Wuensch</surname>
<given-names>W</given-names>
</name>
</person-group>. <article-title>New local field quantity describing the high gradient limit of accelerating structures</article-title>. <source>Phys Rev Spec Topics-Accelerators Beams</source> (<year>2009</year>) <volume>12</volume>:<fpage>102001</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.12.102001</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hopkins</surname>
<given-names>MM</given-names>
</name>
<name>
<surname>Jindal</surname>
<given-names>AK</given-names>
</name>
<name>
<surname>Bussmann</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Ohta</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Berg</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Thomas</surname>
<given-names>C</given-names>
</name>
<etal/>
</person-group> <source>Employing nanoscale surface characterizations in a field emission model</source>. <publisher-loc>Albuquerque, NM, United States</publisher-loc>: <publisher-name>Sandia National Lab.SNL-NM</publisher-name> (<year>2020</year>).</citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Saressalo</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Profatilova</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Paszkiewicz</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Calatroni</surname>
<given-names>S</given-names>
</name>
<etal/>
</person-group> <article-title>Classification of vacuum arc breakdowns in a pulsed dc system</article-title>. <source>Phys Rev Acc Beams</source> (<year>2020</year>) <volume>23</volume>:<fpage>023101</fpage>. <pub-id pub-id-type="doi">10.1103/physrevaccelbeams.23.023101</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dal Forno</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Dolgashev</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Bowden</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Clarke</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Hogan</surname>
<given-names>M</given-names>
</name>
<name>
<surname>McCormick</surname>
<given-names>D</given-names>
</name>
<etal/>
</person-group> <article-title>Rf breakdown tests of mm-wave metallic accelerating structures</article-title>. <source>Phys Rev Acc Beams</source> (<year>2016</year>) <volume>19</volume>:<fpage>011301</fpage>. <pub-id pub-id-type="doi">10.1103/physrevaccelbeams.19.011301</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Descoeudres</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Ramsvik</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Calatroni</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Taborelli</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Wuensch</surname>
<given-names>W</given-names>
</name>
</person-group>. <article-title>Dc breakdown conditioning and breakdown rate of metals and metallic alloys under ultrahigh vacuum</article-title>. <source>Phys Rev Spec topics-accelerators beams</source> (<year>2009</year>) <volume>12</volume>:<fpage>032001</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.12.032001</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nordlund</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
</person-group>. <article-title>Defect model for the dependence of breakdown rate on external electric fields</article-title>. <source>Phys Rev Spec Topics-Accelerators Beams</source> (<year>2012</year>) <volume>15</volume>:<fpage>071002</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.15.071002</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Veske</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Eimre</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Zadin</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
</person-group>. <article-title>Thermal runaway of metal nano-tips during intense electron emission</article-title>. <source>J Phys D: Appl Phys</source> (<year>2018</year>) <volume>51</volume>:<fpage>225203</fpage>. <pub-id pub-id-type="doi">10.1088/1361-6463/aac03b</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Uimanov</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Shmelev</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Oks</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Yushkov</surname>
<given-names>GY</given-names>
</name>
<name>
<surname>Barengolts</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Cathode and plasma phenomena in vacuum-arc sources of hydrogen isotope ions: I. Desorption of hydrogen isotopes during the operation of vacuum arc cathode spots</article-title>. <source>Plasma Sourc Sci Tech</source> (<year>2020</year>) <volume>29</volume>:<fpage>015021</fpage>. <pub-id pub-id-type="doi">10.1088/1361-6595/ab62da</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaufmann</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Cunha</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Benilov</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Hartmann</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Wenzel</surname>
<given-names>N</given-names>
</name>
</person-group>. <article-title>Detailed numerical simulation of cathode spots in vacuum arcs: interplay of different mechanisms and ejection of droplets</article-title>. <source>J Appl Phys</source> (<year>2017</year>) <volume>122</volume>. <pub-id pub-id-type="doi">10.1063/1.4995368</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Eimre</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Parviainen</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Aabloo</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Zadin</surname>
<given-names>V</given-names>
</name>
</person-group>. <article-title>Application of the general thermal field model to simulate the behaviour of nanoscale Cu field emitters</article-title>. <source>J Appl Phys</source> (<year>2015</year>) <volume>118</volume>:<fpage>033303</fpage>. <pub-id pub-id-type="doi">10.1063/1.4926490</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jansson</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Baibuz</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Vigonski</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Zadin</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Parviainen</surname>
<given-names>S</given-names>
</name>
<etal/>
</person-group> <article-title>Growth mechanism for nanotips in high electric fields</article-title>. <source>Nanotechnology</source> (<year>2020</year>) <volume>31</volume>:<fpage>355301</fpage>. <pub-id pub-id-type="doi">10.1088/1361-6528/ab9327</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pohjonen</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Nordlund</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Kuronen</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Fitzgerald</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Dislocation nucleation from near surface void under static tensile stress in Cu</article-title>. <source>J Appl Phys</source> (<year>2011</year>) <volume>110</volume>. <pub-id pub-id-type="doi">10.1063/1.3606582</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Engelberg</surname>
<given-names>EZ</given-names>
</name>
<name>
<surname>Ashkenazy</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Assaf</surname>
<given-names>M</given-names>
</name>
</person-group>. <article-title>Stochastic model of breakdown nucleation under intense electric fields</article-title>. <source>Phys Rev Lett</source> (<year>2018</year>) <volume>120</volume>:<fpage>124801</fpage>. <pub-id pub-id-type="doi">10.1103/physrevlett.120.124801</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zadin</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Pohjonen</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Aabloo</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Nordlund</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
</person-group>. <article-title>Electrostatic-elastoplastic simulations of copper surface under high electric fields</article-title>. <source>Phys Rev Spec Topics-Accelerators Beams</source> (<year>2014</year>) <volume>17</volume>:<fpage>103501</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.17.103501</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bagchi</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Perez</surname>
<given-names>D</given-names>
</name>
</person-group>. <article-title>Atomistic modeling of the coupling between electric field and bulk plastic deformation in fcc metals</article-title>. <source>Phys Rev Acc Beams</source> (<year>2022</year>) <volume>25</volume>:<fpage>033101</fpage>. <pub-id pub-id-type="doi">10.1103/physrevaccelbeams.25.033101</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Laurent</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Tantawi</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Dolgashev</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Nantista</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Higashi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Aicheler</surname>
<given-names>M</given-names>
</name>
<etal/>
</person-group> <article-title>Experimental study of rf pulsed heating</article-title>. <source>Phys Rev Spec topics-accelerators beams</source> (<year>2011</year>) <volume>14</volume>:<fpage>041001</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.14.041001</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kimari</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Zadin</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
</person-group>. <article-title>Biased self-diffusion on Cu surface due to electric field gradients</article-title>. <source>J Phys D: Appl Phys</source> (<year>2022</year>) <volume>55</volume>:<fpage>465302</fpage>. <pub-id pub-id-type="doi">10.1088/1361-6463/ac91dd</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pendyala</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Bobji</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Madras</surname>
<given-names>G</given-names>
</name>
</person-group>. <article-title>Evolution of surface roughness during electropolishing</article-title>. <source>Tribology Lett</source> (<year>2014</year>) <volume>55</volume>:<fpage>93</fpage>&#x2013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.1007/s11249-014-0336-x</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Herring</surname>
<given-names>C</given-names>
</name>
</person-group>. <article-title>Effect of change of scale on sintering phenomena</article-title>. <source>J Appl Phys</source> (<year>1950</year>) <volume>21</volume>:<fpage>301</fpage>&#x2013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.1063/1.1699658</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Herring</surname>
<given-names>C</given-names>
</name>
</person-group>. <source>Fundamental contributions to the continuum theory of evolving phase interfaces in solids: a collection of reprints of 14 seminal papers</source>. <publisher-loc>Berlin, Germany</publisher-loc>: <publisher-name>Springer</publisher-name> (<year>1999</year>). p. <fpage>33</fpage>&#x2013;<lpage>69</lpage>.</citation>
</ref>
<ref id="B30">
<label>30.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mullins</surname>
<given-names>WW</given-names>
</name>
</person-group>. <article-title>Flattening of a nearly plane solid surface due to capillarity</article-title>. <source>J Appl Phys</source> (<year>1959</year>) <volume>30</volume>:<fpage>77</fpage>&#x2013;<lpage>83</lpage>. <pub-id pub-id-type="doi">10.1063/1.1734979</pub-id>
</citation>
</ref>
<ref id="B31">
<label>31.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Suo</surname>
<given-names>Z</given-names>
</name>
</person-group>. <article-title>Motions of Microscopic Surfaces in Materials</article-title>. <source>Adv Appl Mech</source> (<year>1997</year>) <volume>33</volume>:<fpage>193</fpage>.</citation>
</ref>
<ref id="B32">
<label>32.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Suo</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>W</given-names>
</name>
</person-group>. <article-title>A finite element method for simulating interface motion&#x2014;I. Migration of phase and grain boundaries</article-title>. <source>Acta materialia</source> (<year>1997</year>) <volume>45</volume>:<fpage>1907</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1016/s1359-6454(96)00323-0</pub-id>
</citation>
</ref>
<ref id="B33">
<label>33.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Bower</surname>
<given-names>AF</given-names>
</name>
<name>
<surname>Suo</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Shih</surname>
<given-names>C</given-names>
</name>
</person-group>. <article-title>A finite element analysis of the motion and evolution of voids due to strain and electromigration induced surface diffusion</article-title>. <source>J Mech Phys Sol</source> (<year>1997</year>) <volume>45</volume>:<fpage>1473</fpage>&#x2013;<lpage>93</lpage>. <pub-id pub-id-type="doi">10.1016/s0022-5096(97)00013-6</pub-id>
</citation>
</ref>
<ref id="B34">
<label>34.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>WH</given-names>
</name>
<name>
<surname>Srolovitz</surname>
<given-names>DJ</given-names>
</name>
</person-group>. <article-title>Surface morphology evolution in stressed solids: surface diffusion controlled crack initiation</article-title>. <source>J Mech Phys Sol</source> (<year>1994</year>) <volume>42</volume>:<fpage>1551</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1016/0022-5096(94)90087-6</pub-id>
</citation>
</ref>
<ref id="B35">
<label>35.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Srolovitz</surname>
<given-names>DJ</given-names>
</name>
</person-group>. <article-title>On the stability of surfaces of stressed solids</article-title>. <source>Acta metallurgica</source> (<year>1989</year>) <volume>37</volume>:<fpage>621</fpage>&#x2013;<lpage>5</lpage>. <pub-id pub-id-type="doi">10.1016/0001-6160(89)90246-0</pub-id>
</citation>
</ref>
<ref id="B36">
<label>36.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Langtangen</surname>
<given-names>HP</given-names>
</name>
<name>
<surname>Logg</surname>
<given-names>A</given-names>
</name>
</person-group>. <source>Solving PDEs in python: the FEniCS tutorial I</source>. <publisher-loc>Berlin, Germany</publisher-loc>: <publisher-name>Springer Nature</publisher-name> (<year>2017</year>).</citation>
</ref>
<ref id="B37">
<label>37.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kirby</surname>
<given-names>RC</given-names>
</name>
<name>
<surname>Logg</surname>
<given-names>A</given-names>
</name>
</person-group>. <article-title>A compiler for variational forms</article-title>. <source>ACM Trans Math Softw (Toms)</source> (<year>2006</year>) <volume>32</volume>:<fpage>417</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1145/1163641.1163644</pub-id>
</citation>
</ref>
<ref id="B38">
<label>38.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Logg</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Wells</surname>
<given-names>GN</given-names>
</name>
</person-group>. <article-title>DOLFIN: automated finite element computing</article-title>. <source>ACM Trans Math Softw (Toms)</source> (<year>2010</year>) <volume>37</volume>:<fpage>1</fpage>&#x2013;<lpage>28</lpage>. <pub-id pub-id-type="doi">10.1145/1731022.1731030</pub-id>
</citation>
</ref>
<ref id="B39">
<label>39.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kumykov</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Sergeev</surname>
<given-names>I</given-names>
</name>
<name>
<surname>Sozaev</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Gedgagova</surname>
<given-names>M</given-names>
</name>
</person-group>. <article-title>Surface tension of copper in solid phase</article-title>. <source>Bull Russ Acad Sci Phys</source> (<year>2017</year>) <volume>81</volume>:<fpage>357</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.3103/s1062873817030236</pub-id>
</citation>
</ref>
<ref id="B40">
<label>40.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Simon</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Drexler</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Reed</surname>
<given-names>RP</given-names>
</name>
</person-group>. <source>Properties of copper and copper alloys at cryogenic temperatures. Final report</source>. <publisher-loc>Boulder, CO, United States</publisher-loc>: <publisher-name>National Inst. of Standards and Technology MSEL</publisher-name> (<year>1992</year>).</citation>
</ref>
<ref id="B41">
<label>41.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Karimi</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Tomkowski</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Vidali</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Biham</surname>
<given-names>O</given-names>
</name>
</person-group>. <article-title>Diffusion of Cu on Cu surfaces</article-title>. <source>Phys Rev B</source> (<year>1995</year>) <volume>52</volume>:<fpage>5364</fpage>&#x2013;<lpage>74</lpage>. <pub-id pub-id-type="doi">10.1103/physrevb.52.5364</pub-id>
</citation>
</ref>
<ref id="B42">
<label>42.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schulze Icking-Konert</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Giesen</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Ibach</surname>
<given-names>H</given-names>
</name>
</person-group>. <article-title>Decay of Cu adatom islands on Cu(111)</article-title>. <source>Surf Sci</source> (<year>1998</year>) <volume>398</volume>:<fpage>37</fpage>&#x2013;<lpage>48</lpage>. <pub-id pub-id-type="doi">10.1016/s0039-6028(98)80009-5</pub-id>
</citation>
</ref>
<ref id="B43">
<label>43.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kallinteris</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Evangelakis</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Papanicolaou</surname>
<given-names>N</given-names>
</name>
</person-group>. <article-title>Molecular dynamics study of the vibrational and transport properties of copper adatoms on the (111) copper surface; comparison with the (001) face</article-title>. <source>Surf Sci</source> (<year>1996</year>) <volume>369</volume>:<fpage>185</fpage>&#x2013;<lpage>98</lpage>. <pub-id pub-id-type="doi">10.1016/s0039-6028(96)00920-x</pub-id>
</citation>
</ref>
<ref id="B44">
<label>44.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Evangelakis</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Papageorgiou</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Kallinteris</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Lekka</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Papanicolaou</surname>
<given-names>N</given-names>
</name>
</person-group>. <article-title>Self-diffusion processes of copper adatom on Cu(110) surface by molecular dynamics simulations</article-title>. <source>Vacuum</source> (<year>1998</year>) <volume>50</volume>:<fpage>165</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1016/s0042-207x(98)00039-6</pub-id>
</citation>
</ref>
<ref id="B45">
<label>45.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Timko</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Nordlund</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Costelle</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Matyash</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Schneider</surname>
<given-names>R</given-names>
</name>
<etal/>
</person-group> <article-title>Mechanism of surface modification in the plasma-surface interaction in electrical arcs</article-title>. <source>Phys Rev B</source> (<year>2010</year>) <volume>81</volume>:<fpage>184109</fpage>. <pub-id pub-id-type="doi">10.1103/physrevb.81.184109</pub-id>
</citation>
</ref>
<ref id="B46">
<label>46.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wuensch</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Degiovanni</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Calatroni</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Korsb&#xe4;ck</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Djurabekova</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Rajam&#xe4;ki</surname>
<given-names>R</given-names>
</name>
<etal/>
</person-group> <article-title>Statistics of vacuum breakdown in the high-gradient and low-rate regime</article-title>. <source>Phys Rev Accel Beams</source> (<year>2017</year>) <volume>20</volume>:<fpage>011007</fpage>. <pub-id pub-id-type="doi">10.1103/physrevaccelbeams.20.011007</pub-id>
</citation>
</ref>
<ref id="B47">
<label>47.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gill</surname>
<given-names>V</given-names>
</name>
<name>
<surname>Guduru</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Sheldon</surname>
<given-names>B</given-names>
</name>
</person-group>. <article-title>Electric field induced surface diffusion and micro/nano-scale island growth</article-title>. <source>Int J Sol Structures</source> (<year>2008</year>) <volume>45</volume>:<fpage>943</fpage>&#x2013;<lpage>58</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijsolstr.2007.09.010</pub-id>
</citation>
</ref>
<ref id="B48">
<label>48.</label>
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Hopkins</surname>
<given-names>MM</given-names>
</name>
<name>
<surname>Bussmann</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Smith</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Scrymgeour</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Brumbach</surname>
<given-names>MT</given-names>
</name>
<name>
<surname>Hjalmarson</surname>
<given-names>HP</given-names>
</name>
<etal/>
</person-group> <source>Vacuum field emission models measurements and simulations</source>. <publisher-loc>Albuquerque, NM, United States</publisher-loc>: <publisher-name>Sandia National Lab.SNL-NM</publisher-name> (<year>2017</year>).</citation>
</ref>
<ref id="B49">
<label>49.</label>
<citation citation-type="web">
<person-group person-group-type="author">
<name>
<surname>Shinohara</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Bagchi</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Simakov</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Baryshev</surname>
<given-names>SV</given-names>
</name>
<name>
<surname>Perez</surname>
<given-names>D</given-names>
</name>
</person-group>. <article-title>Thermal and electric field driven breakdown precursor formation on metal surfaces</article-title> (<year>2023</year>). <comment>
<ext-link ext-link-type="uri" xlink:href="https://arxiv.org/abs/2311.02519">https://arxiv.org/abs/2311.02519</ext-link>. (Accessed 29 February 2024)</comment>.</citation>
</ref>
<ref id="B50">
<label>50.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nusinovich</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Kashyn</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Antonsen</surname>
<given-names>T</given-names>
<suffix>Jr</suffix>
</name>
</person-group>. <article-title>Possible role of rf melted microparticles on the operation of high-gradient accelerating structures</article-title>. <source>Phys Rev Spec Topics-Accelerators Beams</source> (<year>2009</year>) <volume>12</volume>:<fpage>101001</fpage>. <pub-id pub-id-type="doi">10.1103/physrevstab.12.101001</pub-id>
</citation>
</ref>
<ref id="B51">
<label>51.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Xanthakis</surname>
<given-names>J</given-names>
</name>
</person-group>. <article-title>Extension of the general thermal field equation for nanosized emitters</article-title>. <source>J Appl Phys</source> (<year>2016</year>) <volume>119</volume>. <pub-id pub-id-type="doi">10.1063/1.4940721</pub-id>
</citation>
</ref>
<ref id="B52">
<label>52.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kyritsakis</surname>
<given-names>A</given-names>
</name>
</person-group>. <article-title>General form of the tunneling barrier for nanometrically sharp electron emitters</article-title>. <source>J Appl Phys</source> (<year>2023</year>) <volume>133</volume>. <pub-id pub-id-type="doi">10.1063/5.0144608</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>