## A Data-Driven Model for the Field Emission from Broad-Area Electrodes

Moein Borghei, Robin Langtry a Avalanche Energy, 9100 E Marginal Way S, Tukwila, 98108, WA, USA

## Abstract

Electron emission from cathodes in high field gradients is a quantum tunneling effect. The 1928 FowlerNordheim field emission (FE) equation and the 1956 Murphy-Good FE equation have traditionally been key in describing cold field emissions, offering estimates for emitters for almost a century. Nevertheless, applying FE theory in practice is often constrained by the lack of data on the distribution and geometry of the emission sites. Predictions become more challenging with an uneven electric field distribution at the cathode surface. Consequently, FE formulations are frequently calibrated using current-voltage data after test, limiting their efficacy as true predictive models.

This study develops an alternative model for field emission using a data-driven predictive approach based on (1) vast experimental data, (2) electrostatic simulations of the cathode surface, and (3) detailed material and geometry properties, which together overcome these limitations. The objective of this work is to develop and harness this comprehensive dataset to train a machine learning model capable of providing precise predictions of the cathode current in order to further the understanding and application of field emission phenomena. More than 259 hours of experimental data have been processed to train and benchmark some of the well-known machine learning models. After two stages of optimization, a coefficient of determination &gt; 98% is achieved in the prediction total field emission current using ensemble models.

Keywords: field emission, Fowler-Nordheim theory, machine learning, quantum tunneling

## 1. Introduction

Conductors are often used to transport electrical charges. In most applications, these electrons are totally confined within the material by the presence of a surface potential barrier known as the material's work function. Until the early twentieth century, the emission of electrons from conductors was engineered through thermionic and photoelectric emission processes.

on the 1D Schr¨dinger equation, Fowler-Nordheim o derived an analytic formulation to approximate the tunneling rate of electrons through the potential barrier as a function of the local electric field and the work function of the material. Stern found and corrected a numerical error in the Fowler-Nordheim equation, and also introduced the concept of what is now called the Fowler-Nordheim plot [10].

Interest in field emission arose in the 1910s through J. E. Lillienfield's work developing portable X-ray machines [1], and later through W. D. Coolidge's observation of a new phenomenon of electron emission from the cold cathode due to the high electric field strength at the surface [2]. Following a series of experimental and theoretical works [3, 4, 5, 6, 7, 8], Fowler and Nordheim derived a closed-form equation for the local emission current density of an emission site in 1928 [9]. Based

Email address: mborghei@avalanchefusion.com (Moein Borghei)

Notable errors were found in the 1920s analytic solution of Fowler-Nordheim. These were corrected in 1956 by Murphy and Good [11] who proposed a corrected FE equation predicting considerably higher local current emission densities. Since then, the local electric field of the emitter is explained by field enhancement factors. Traditionally, sharp metallic whiskers that can be created artificially or developed through surface dislocation activity [12] or surface diffusion [13] under high field strength were thought to cause enhancement (see Figure 1). Over the years, efforts have been made to approximate the shape of field enhancing features [14, 15, 16, 17]. Despite these efforts, the observa-

Figure 1: Exemplary illustration of the field enhancement features.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-748e295a-3ea1-4874-adbf-1d36eeebc02c.jpg)

tion of field emission current on an extremely polished cathode surface indicates alternative sources of field emission, apart from whiskers. This has been confirmed in [18]. Additionally, an observed shift in electron emission spectra to energies below the Fermi level [19, 20] and the cathode electroluminescence [21], further point towards a nonmetallic origin.

Aside from the debates over the origin of the field emission current, nearly a century later after the discovery of the field emission, there is still no widely accepted model for the prediction of total field emission for practical applications involving broad-area electrodes, large gaps, and nonuniform electric field distribution where numerous emission sites make up the total field emission current, or dark current . In practice, it is nearly impossible to acquire all the geometric details of the microscopic emission sites on a broad-area electrode.

In this work, our aim is to train a data-driven model based on more than 259 hours of experimental data. We supplement the input dataset with electrostatic, geometric, and material information about the cathode. The model is agnostic to the individual field emission sites and rather aims to recognize their collective pattern on a broad-area electrode. To complement the cathode voltage and current measurements, we utilize the electric field distribution from an electrostatic simulation, which yields statistical features that explain the probability distribution of the electric field magnitude. Optical microscopy of a 1 µm × 1 µm section of a cathode data provides statistical features of the projection heights of the surface structure. Lastly, the work function of the material and the total cathode area conclude the input dataset.

Relaxing the modeling of broad-area emission from the traditional curve fitting to a more datainformed complex model allows one to capture a wider variation of data. Moreover, such a model can be constantly reinforced with more test results, provided that the aforementioned set of data is provided. Performance comparison of different supervised machine learning models on the dataset shows that we can achieve &gt; 98% accuracy in the prediction of the total field emission current with the gradient boosting ensemble model.

The structure of the remainder of this article is delineated as follows: In Section 2, we provide a brief review of the Murphy-Good theory of field emission with a focus on key aspects of the theory relevant to our work. Section 3 elucidates the framework of our approach to field emission prediction, which will be discussed in more detail in subsequent sections. The experimental setup is briefly discussed in Section 4. Section 5 details the preprocessing of data for application in machine learning models. The results of training and validation of various machine learning models are presented in Section 6. Finally, a discussion of the results and the conclusions are presented in Sections 7 and 8, respectively.

## 2. Overview of Murphy-Good Field Emission Theory

Murphy and Good [11] used quantum tunneling theory to calculate the probability that an electron escapes from a metal surface under a high electric field, which is then integrated into the energy states to obtain the local emission current density. The zero-temperature version of the Murphy-Good equation for steady electric fields is as follows:

$$J _ { L } = \frac { a \, E ^ { 2 } } { \phi \, t ^ { 2 } ( y ) } \ e ^ { - \frac { b \, \phi ^ { 3 / 2 } \ v ( y ) } { E } } \quad \quad ( 1 )$$

where a and b are constants, E is the intensity of the local electric field at the emitter, and ϕ is the work function of the material. Finally, t ( y ) and v y ( ) are slowly varying functions of ϕ and E . t ( y ) is close to unity, but v y ( ) takes values slightly less than unity.

The deviation of v y ( ) from unity is responsible for the large difference between the predictions of the Fowler-Nordheim and the Good-Murphy equation. Forbes and Deane provide a simple yet accurate approximation for v y ( ) in [22], and Miller reports the tabulated data in [23]. For additional review of recent contributions to the Nordheim parameter and more, see [24, 25].

The local electric field at the emission site, denoted E , frequently exceeds the expected magnitude derived from the estimation of the macroscopic field E 0 , by two to three orders of magnitude. Consequently, an effective enhancement factor, β eff = 1 I ∫ S β ( x ) J ( x ) dS , is typically used to approximate the microscopic electric field [26]. However, using E = β eff E 0 usually does not yield accurate results. A more practical approach is to use the dimensionless field enhancement factor at its apex γ a , which defines the local electrostatic field as E = γ E a 0 . -assuming a cylindrically symmetric protrusion [27].

To derive the total emission current from an emission site, the local emission current density can be associated with the total emission current via the relation I = AJ L , where A represents the notional emission area of the emitter which is often field dependent [22].

## 2.1. From Single Emitter to Multiple Emitters

In practical applications, the measured field emission current includes the contribution of not one but many emitters. Ideally, the field enhancement factor, work function, and emission area would be known for each emission site. However, it is merely impossible to do so except for the case of electrodes which have been specifically engineered for a particular enhancement factor.

Tomaschke and Alpert [28] examined the collective field emission of 100 emission sites, where each site is randomly assigned a field enhancement factor and an emission area. The total accumulation of the emitters' contributions to the total current is then given as:

$$I = \sum _ { i = 1 } ^ { N } A _ { i } \frac { a \, E _ { i } ^ { 2 } } { \phi _ { i } \ t ( y ) ^ { 2 } } \ e ^ { - \frac { b \, \phi _ { i } ^ { 3 / 2 } \ v ( y ) } { E _ { i } } } \quad ( 2 ) \quad \text{cr} \ \text{tic} \ \text{hc}$$

Their findings reveal that the total emission is closely aligned with a linear Fowler-Nordheim plot. Therefore, they can determine an effective field enhancement factor and a total field-emitting area. Several studies [29, 30] have sought to find an effective field enhancement factor for broad-area cathodes. The common observation of these studies is that as the applied electric field increases, β eff decreases. As the electric field increases, the lower β eff areas on the surface begin to contribute more to the local current density.

The inability of these models to adapt to more complex geometries and to provide a systematic approach to predict the field emission current from practical electrode geometries necessitates the development of a novel model based on experimental data, agnostic to the individual emitter characteristics while correctly capturing the collective behavior of broad-area electrodes in terms of field electron emission.

## 3. Predictive Modeling of Field Emission Current Using Data-Driven Approaches

Based on the theoretical overview of field emission, the key parameters of field emission are classified into these categories: (a) electrostatic parameters that describe the macroscopic field gradient at the cathode surface, (b) geometric parameters that influence the microscopic field gradient at the cathode surface and the emitting area, and (c) the chemical composition of the cathode surface. This study assumes that the vacuum level and the residual gas composition remain reasonably constant across different tests.

In this study, the data collected are processed and utilized to train a machine learning model to predict the field emission current of the cathode based on a given voltage. The objective is to transition from curve fitting the Fowler-Nordheim equation to a generalized, nonlinear model that more accurately captures the relationship between cathode voltage and current, and that better incorporates the described statistical features. This framework is depicted in Figure 2.

At the center of this process is machine learning (ML) model training, which integrates input from three primary groups: experimental measurements, electric field simulation, and optical microscopy. Additionally, factors such as work function and cathode area are considered, ensuring a holistic approach to training.

Starting with Experimental Measurements, the process begins by measuring the current-voltage ( I -V ) characteristics. This involves calculating the voltage drop across the series resistor and subsequently adjusting the measured voltage and current to obtain the actual cathode voltage, V cathode and current, I cathode . V cathode is used as one of the input parameters and I cathode is the target quantity for the approximation.

In parallel, Electric Field Simulation involves simulating the geometry explained in Section 4, fil-

Figure 2: Flowchart illustrating the feature extraction and machine learning model training process for field emission.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-9d62bbdf-2933-4ed3-8c47-a04022736a82.jpg)

tering fields above 5 MV/m to focus on significant emission areas, applying a Gaussian mixture model (GMM) for statistical analysis, and extracting relevant features that make up the input matrix.

includes empirical measurements, theoretical simulations, and microscopic surface characteristics, allowing accurate and reliable predictions of field emission behavior.

Optical microscopy is used to collect detailed surface information. This involves obtaining depth profiles, determining probability distributions, fitting skewed Gaussian models, and extracting roughness features. Although the emission area and field enhancement factor are not explicitly provided to the data-driven models, the statistical features derived from microscopy images supply the model with information about the density and profile of emitters. The machine learning models capture the variations in the emission area with the electric field through correlations in the provided data.

Lastly, two features are added to complete the input features vector. One is the cathode work function, which affects electron emission properties, and the second is the total cathode area which is used to extrapolate the emitter distribution from the localized microscopy images to the whole surface. Together, these diverse inputs ensure that the ML model is trained on a robust dataset, which

## 4. Experiment Setup

Figure 3 shows the experimental setup for field emission testing on different cathode materials. The electrode assembly involves a negatively charged cathode, a dielectric spacer (made from MACOR), and a grounded anode (fabricated from stainless steel). The assembly is housed within a vacuum chamber, where a roughing pump and a turbo pump maintain a pressure below 1 × 10 -8 Torr . The high voltage is transmitted from the power supply to the cathode using a cable along with a resistor limiting current 1 M Ω. The power supply voltage and current are recorded at a rate of one sample per second. The voltage is increased incrementally to follow a current conditioning method from 0 kV to 80 kV . During the experiment, measurements of the anode current, vacuum pressure,

Figure 3: Schematic representation of the apparatus for the measurement of field emission current.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-3c2779ec-3401-4810-ae24-b8313652925c.jpg)

and X-ray intensity are taken along with the voltage and current read at the power supply terminal. Observations of discharges are made visually through a chamber window using a high-speed camera equipped with wide-field optics to detect any discharge inside the vacuum chamber.

## 4.1. Sample Preparation

All specimens for testing are made from bulk materials. The selection of cathode materials included molybdenum (type 361) with a purity of 99.95%, titanium (grade 2), and copper (grade 101). The specimens were shaped using computer numerical control (CNC) machining. The desired surface roughness is achieved through manual sanding (2 µm or finer) or grit-blasting for increased surface roughness (4 -10 µm ). Figure 4 illustrates the characteristic copper and molybdenum specimen under 10x magnification. Following fabrication, the samples are subjected to a cleaning protocol to remove hydrocarbons and other residues that remain from the manufacturing process. Initially, each sample is washed with hot soapy water in a sink to remove bulk contaminants. Next, each sample is subjected to a 15-minute ultrasonic cleaning cycle with hot soapy water. After being rinsed, each sample is then transferred to another ultrasonic cleaner containing deionized water for a brief two-minute cleansing. Subsequently, the sample is immersed in an ultrasonic cleaner filled with isopropyl alcohol (IPA) for an additional 15 minutes. The final step includes drying the sample with a lint-free cloth in a class-10,000 cleanroom.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-2d067f80-243e-45be-90d9-1437339f9267.jpg)

(a) Polished molybdenum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-394af198-34e4-4160-8bce-ba325d384e8b.jpg)

(b) Grit blasted copper.

Figure 4: The surface morphology of two specimens magnified at 10x under the 3D surface profiler laser microscope.

## 5. Data Preprocessing

To prepare the raw experimental and simulation data for use in a data-driven model, the data must be cleaned and enriched with additional information. The flow chart in Figure 2 shows the steps discussed in the rest of the section.

## 5.1. Experimental Data

The first step is the removal of the data points where, due to electromagnetic interference and transient faults, the voltage or current has not been properly stored. In the next step, the tests are automatically separated by tracing the rise and fall of the cathode voltage in the dataset. Any test missing information on the cathode condition is deleted. Data points where the voltage setpoint and the voltage measurement deviate by greater than 10% are also removed.

The data collected from each test are acquired only once the cathode has been conditioned to a specific voltage level. The conditioning process effectively eliminates problematic emitters and adsorbates by gradually increasing the voltage until the current stabilizes. As a result, only stable field

emission current values are considered, facilitating the reproducibility of the results.

The experimental data used in this model are measurements of cathode voltage, V cathode , and cathode current, I cathode . However, the experimental setup measures the current of the power supply, I PS , and the voltage, V PS . Since there is a series resistor, R series , in line between the power supply and the vacuum chamber, the measured voltage ( V PS ) includes an additional voltage due to the drop across the resistor. Therefore, the cathode voltage is calculated as V cathode = V PS -R series × I PS .

Additionally, the current, I PS , captures not only the cathode current but also any leakage current, I leakage , such as that of the corona in the exposed electrodes. Therefore, the baseline leakage current as a function of the power supply voltage is deducted from the actual current measurement at each corresponding voltage level. By deducting the leakage current from the measured current, we acquire the cathode current I cathode = I PS -I leakage . Ultimately, we obtain a matrix containing numerous pairs of V cathode and I cathode , with each pair representing the total field emission current of a specific cathode under a particular voltage.

## 5.2. Electrostatic Simulation

The field emission of electrons is strongly related to the strength of the electric field on the cathode surface, as discussed in Section 2. For a single emitter, the field at that emission site is needed to calculate the emission current based on the Fowler-Nordheim expression. However, in the case of broad-area emission, the electron current stems from a multitude of emission sites. The apex (or maximum) value of the local electric field at each of these emission sites depends on (1) the macroscopic electric field and (2) the microscopic field enhancement factor. In this study, the objective is to incorporate these parameters statistically. As a first step, the macroscopic electric field at the surface of the cathode is simulated electrostatically.

Figure 5 shows the distribution of the electric field in the periphery of the cathode. To feed the machine learning model, this study uses statistical features that describe the inhomogeneity of the electric field distribution. For a virgin electrode, there is a threshold field in the range of 5 -15 MV/m associated with the initial switch-on of the field emission current [31]. Therefore, the values of the electric field below 5 MV/m are excluded from the probability distribution of the electric field.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-91c7ce33-df4d-45d3-ad95-7d01e88928cb.jpg)

(a) Electric potential distribution.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-93881f89-60c9-45c0-9890-270ba4922c1a.jpg)

(b) Electric field distribution.

Figure 5: The electrostatic simulation of the high voltage assembly with cathode at -80kV.

AGaussian mixture model (GMM) is used to represent an inhomogeneous electric field distribution at the surface of a cathode by assuming that the field distribution is composed of multiple Gaussian components, each representing a distinct region with its own mean and variance. This model allows for the accommodation of complex and multi-modal field distributions, which are common in practical scenarios due to irregularities in material properties and geometric configurations. Each component in the GMM can be thought of as a local 'peak' or 'valley' in the field intensity, and the overall model provides a probabilistic framework for estimating

Figure 6: Probability distribution of electric field at the surface of cathode.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-ba32e5e6-14de-4500-9acb-dd4468a83800.jpg)

the likelihood of field strengths at different points on the cathode surface. The sensitivity analysis of model prediction with respect to the number of components are discussed in section 7.

Figure 6 demonstrates the distribution of the electric field intensity using four Gaussian distributions. Therefore, eight statistical features are used to explain the contribution of the electric field to the field emission in the model. After benchmarking different numerical models. A sensitivity analysis of the number of components in the GMM approximation is performed to find the optimal mixture model.

## 5.3. Optical Microscopy

For broad-area electrodes, it is well-established that field enhancement can exhibit significant spatial variability. Consequently, an effective approximation of surface microprotrusions is essential to capture the extensive variance in the currentvoltage measurements. In this study, we employ laser microscopy to characterize the surface profile of the sample, thereby enabling a detailed examination of the microscopic features of the cathode surface. The chemical composition has previously been analyzed by the authors in [32] and is not repeated here for brevity.

Surface roughness is determined with a Keyence VKX3050 3D surface profiler. Each specimen is examined using the laser microscope and the height array for an area of 1 µm × 1 µm at the center of the cathode surface. Due to the uniform surface preparation process, the 1 µm 2 area can provide a valid representation of the total cathode area. Figure 7 depicts the probability distributions of the surfaces of two specimens shown in Figure 4. If one plots the histogram of the height distribution, a skewed Gaussian distribution can approximate the projection heights of the probability distribution. Therefore, three statistical features incorporate a microscopical fingerprint.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-81428316-966c-48ef-ba54-5132f2af9a13.jpg)

(a) Polished molybdenum.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-1af3ca94-9de7-4110-8918-10a357cb94dc.jpg)

(b) Grit blasted copper.

Figure 7: The probability distribution of projection heights for two specimens.

## 6. Numerical Results

After the data have been processed and the datasets have been created, a collection of physical and environmental parameters is used as input features, and the cathode current vector serves as the target vector. The layout of a single row of input/output features can be seen in Figure 8. Only stable field emission current data points postconditioning are utilized in the modeling. The total count of data points is 920, which will be divided

Figure 8: The structure of single datapoint in the input and output arrays.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-87d70c32-dcd5-49b0-825b-454773892fa0.jpg)

into 80% for training and 20% for testing. This section benchmarks several well-known models. Given the widespread use of these models in different areas of machine learning, descriptions of each algorithm are avoided (for comprehensive information on these models, refer to [33]).

The hyperparameters of each model can be adjusted to improve the accuracy of the prediction. We begin with simple models such as linear regression, and then move toward more advanced models. Performance metrics used include: (1) mean absolute error (MAE) that calculates the average absolute difference, treating all errors equally, (2) mean squared error (MSE), which measures the average squared difference between predicted and actual values, penalizing larger errors more severely, and (3) the coefficient of determination ( R 2 score) that quantifies the proportion of variance in the dependent variable explained by the model, with 1 indicating perfect prediction and 0 indicating no explanatory power. To find the optimal hyperparameters, we utilize the Bayesian optimization method, specifically the tree-structured Parzen estimator, with the aim of maximizing the R 2 score. The models are developed in Python using open-source libraries, including Pandas , NumPy , TensorFlow , Scikit-learn , Optuna , and XGBoost .

## 6.1. Linear Regression

Linear regression is the simplest regression model if a linear correlation is anticipated between the input and the target variable. In the first attempt, we have employed Ridge and Lasso in linear regression as regularization strategies to attenuate overfitting.

The focus of hyperparameter tuning is on the regularization parameter α , exploring a logarithmic scale from 10 -4 to 10 4 . Through 1000 trials, the optimal α values for Ridge and Lasso were found to be approximately 26.58 and 1.07, respectively. With these optimal parameters, the models were trained on the training dataset and then evaluated on the test dataset. The Ridge model achieves an MSE of 1414.1, an MAE of 26.9, and an R 2 score of 43.8%.

In contrast, the Lasso model achieves an MSE of 1529.8, an MAE of 27.4, and an R 2 score of 39.2%. The R 2 score indicates that the Ridge regression model explains approximately 43. 8% and 39. 2% of the variance in the target variable, respectively. The performance of these models indicates a necessity for more sophisticated modeling techniques to better capture the highly nonlinear relationship between field emission current and its input features, thereby enhancing predictive accuracy.

## 6.2. Support Vector Regression

Support Vector Regression (SVR) has shown to be effective in high-dimensional spaces and when a certain tolerance margin within the model is required. It often achieves good performance with non-linear relationships when the suitable kernel is employed. In our implementation of SVR, hyperparameter optimization is carried out focusing on the regularization parameter, α , which ranged from 10 -4 to 10 , 4 ϵ spanning from 10 -4 to 10 2 defining the tolerance margin where no penalties are applied for errors, kernel types including linear, polynomial, radial basis function (RBF), and sigmoid. From more than 1000 trials, the optimal parameters determined were α = 4967 8, . ϵ = 0 27, RBF kernel, . and its tuning parameter γ = 3 06, which describes . the extent of influence of a single training example.

Using these optimal parameters, the SVR model is trained on the full training dataset and evaluated on the test dataset. The SVR model achieves an MSE of 1051 67, an MAE of 18 67, and an R . . 2 score of 58.2%. The improvement in predictive accuracy, compared to the Ridge and Lasso regression models, highlights the higher efficacy of SVR in capturing nonlinear relationships within this set of data. Nonetheless, the results suggest potential for further enhancements in predictive performance.

## 6.3. Random Forest Regression

In the next model, we explore random forest regression (RFR), an ensemble learning method that builds multiple decision trees and combines them to

achieve more accurate and stable predictions. This approach effectively handles complex interactions and nonlinearities, providing better generalization and feature importance scores, though it is less interpretable and computationally intensive. The key optimization parameters included the number of estimators (from 100 to 10000), maximum tree depth (from 10 to 1000), minimum samples for splitting (from 2 to 20) and leaf nodes (from 1 to 20), maximum number of features for splitting (from 4 to 14) and bootstrap sampling. Through 1000 optimization trials, the optimal parameters identified were 7965 estimators, a maximum depth of 316, a minimum of 2 samples for node splitting, 2 samples per leaf node, the maximum number of features of 14 without bootstrapping.

The RFR model achieved an MSE of 82 6, . an MAE of 4 4, and an R . 2 score of 96 7%. . This high R score indicates that the random forest model ac2 counts for approximately 97% of the variance in the target variable, demonstrating its superior performance in terms of capturing complex relationships within the dataset compared to previous models.

## 6.4. Gradient Boosting Regression

Further extending our analysis, we applied gradient boost regression (GBR), another ensemble technique that builds models sequentially to focus on the errors of previous models, often enhancing accuracy. GBRs are effective in achieving high predictive accuracy, but can be prone to overfitting, sensitive to noisy data and outliers, and require meticulous parameter tuning. The optimization process involved parameters such as the number of estimators (100 to 10,000), learning rate (0.01 to 0.5), maximum depth (1 to 10), minimum samples to split a node (2 to 20), minimum samples in a leaf node (1 to 20), subsample ratio (0.25 to 1.0), and the maximum features to split. After 1000 trials, the optimal parameters identified were 7895 estimators, a learning rate of 0.43, a maximum depth of 2, 2 samples to split a node, 10 samples per leaf node, a subsample ratio of 0.985, and the maximum number of features of 13.

The GBR model achieved an MSE of 42.1, an MAE of 3.2, and an R 2 score of 98 2%. . This high R 2 score indicates that the gradient boost model accounts for more than 98% of the variance in the target variable. The significant reduction in error metrics and the impressive R 2 score underscore the robustness of gradient boosting to capture complex data patterns and relationships.

## 6.5. Extreme Gradient Boosting Regression

Compared to GBR, Extreme Gradient Boosting Regression (XGBR) uses more advanced regularization techniques, handles missing values more efficiently, provides parallel and distributed computing, and offers more options for customization. The parameters considered include the column subsampling ratio (varying from 0.1 to 1.0), which determines the proportion of features used for each tree, the learning rate (from 0.01 to 0.3), the maximum depth of the trees (from 3 to 50), the L 1 regularization term on weights, α (varying from 0.1 to 100), the L 2 regularization term on weights, λ (from 10 -8 to 1), and the number of boosting rounds (from 100 to 1000). Additional parameters include the minimum loss reduction required for the partition of leaf nodes, γ (from 10 -8 to 1), the minimum weight required to create a new tree node (from 1 to 10), and the controlling factor for the balance of positive and negative weights (between 0.1 and 10).

The best tuning results for the XGBR model indicate that the optimal parameters include a tree feature selection rate of approximately 36%, a learning rate of about 29.4%, and a maximum tree depth of 29. Regularization is achieved with an α value of approximately 0.15 and a λ value near 0.002, while the model is trained using 604 estimators. The optimal model also has a γ value of around 0.003, a sample rate of approximately 64%, and a minimum child weight of 1. Additionally, the scale for positive class weights is set at about 7.56. The best model achieves an R 2 score of 90.9%, an MSE of 229.0 and an MAE of 7.5. The results demonstrate a relatively robust performance of the optimal XGBR model, but inferior to GBR and RF, due to the ineffectiveness of the extra regularization parameters given the size of the dataset and the number of features.

## 6.6. Neural Network Regression

Neural networks are highly flexible and powerful, capable of learning very complex patterns and relationships in large and complex datasets. They offer significant flexibility, but require careful tuning, are prone to overfitting, and generally lack interpretability. More data and computational resources might be needed.

The neural network architecture is optimized by fine-tuning several crucial parameters. These include the number of hidden layers (varying from

Figure 9: Performance comparison among various machine learning models.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00301v2-figure-267c7862-5bdf-4b11-8909-59ebaf87f96c.jpg)

1 to 8) and the number of units per layer (ranging from 16 to 128), the activation function (chosen from linear, ramp, softmax, sigmoid, hyperbolic tangent, exponential and scaled exponential linear units.), the dropout rate (fluctuated between 0.0 and 0.5), the learning rate ( η ) from 1 e -5 to 1 e -1 , and the optimizer including various versions of adaptive moment estimation (Adam, Adamax, NAdam), root mean square propagation (RMSProp) and gradient-based optimization methods (SGD, AdaDelta, AdaGrad, FTRL). In addition, the number of epochs can be up to 200 with early stopping on patience with 5 epochs, and the batch size is adjusted between 32 and 64.

The optimal hyperparameters identified were as follows: a network architecture consisting of 2 layers (first with 32 neurons and 12% dropout, and second with 128 neurons and 31% dropout), ramp activation function, RMSProp optimizer, a 2.4% learning rate and a batch size of 32. The best result from the neural network model achieves MSE, MAE and R 2 values of 922.9, 20.3 and 63.4%, respectively. Although this performance surpasses that of linear models, it does not match the effectiveness of ensemble methods. The low number of layers in the optimal shows the need for a larger dataset and possibly a higher number of features to optimize their performance.

## 7. Discussion

The performance of various regression models on the given dataset, indicated by their R 2 scores, is summarized in Figure 9. The results reveal that traditional linear models such as Ridge and Lasso exhibit relatively low R 2 scores of 43.8% and 39.2%, respectively, suggesting a limited ability to capture complex relationships in the data. In contrast, Support Vector Regression achieves a slightly higher score of 58.2%, indicating its ability to handle nonlinear patterns marginally better than Ridge and Lasso. However, it is the ensemble methods that significantly outperform the others: random forest, gradient boosting, and extreme gradient boosting show R 2 scores of 96.7%, 98.3%, and 90.0%, respectively. These high scores highlight the strength of ensemble learning in reducing variance and bias through the aggregation of multiple models. It is of note that the performance of the neural network has been subpar compared to the ensemble models, with an R 2 score of 63.3%.

The superior performance of ensemble methods such as random forest and gradient boosting can be attributed to their inherent ability to combine the strengths of multiple decision trees, thereby enhancing the overall robustness and accuracy of the model. Random forest achieves high performance through bagging, which reduces overfitting by averaging multiple decision trees trained on different subsets of the data. Gradient boosting, on the other hand, builds trees sequentially, each correcting the errors of its predecessor, leading to a highly accurate model, as evidenced by its highest R 2 score of 98.3. Extreme gradient boosting, a more efficient and regularized version of gradient boosting, also performs exceptionally well with a score of 90.0, balancing computational efficiency and accuracy. Neural networks, with their deep learning capabilities, perform better than traditional linear models but fall short of ensemble methods, likely due to their need for extensive tuning and larger datasets to achieve peak performance. The variation in these scores underscores the importance of model selection based on the specific characteristics and complexity of the dataset at hand.

A sensitivity analysis on the effect of the GMM components (1 to 9) on the electric field approximation is also performed to further enhance the accuracy of the GBR prediction. However, the R 2 score remains between 97.0% and 98.6% without any trend, indicating minimal impact from varying electric field approximation complexity. This is somewhat expected because of the identical cathode geometry in all data points. The impact is expected to be more significant with the varied cathode geometries.

## 8. Conclusion

In pursuit of a robust yet practical model to predict total field emission from broad-area electrodes, this study approached the problem with an emphasis on extensive experimental data. We paired the experimental data with the electrostatic field distribution, surface profile, and material information to build a matrix of input features for prediction. Upon benchmarking various machine learning models, we found that the ensemble models, especially the gradient boosting regression, perform remarkably well, achieving over 98% accuracy. Neural neural networks of various depths did not offer improved results.

The findings presented herein should not be considered definitive. As the number of features and/or the size of the dataset expands, more sophisticated models may become necessary. These advanced models will potentially encapsulate a broader range of variability and, it is anticipated, will eventually evolve into a sufficiently robust framework capable of accurately predicting the field emission current from an electrode. One should also note that more statistical features may hinder the practicality and cost of using such models.

Future studies will be directed towards utilizing more extensive datasets that cover a variety of cathode geometries and involve multiple work functions for each material corresponding to different crystal faces. Additionally, integrating additional physical phenomena into the simulation models can broaden the statistical distribution related to the field emission of unconditioned cathodes.

## Acknowledgment

We express our gratitude to Alexander Fiore, Madeline Vorenkamp, Craig Miller and Rod Latham for their insightful discussions and thorough review of our manuscript.

## Data Availability

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

## References

- [1] J. Lilienfeld, The autoelectronic discharge and its application to the construction of a new form of x-ray tube, American Journal of Roentgenology (1922) 172-179.
- [2] W. D. Coolidge, Use of very high-voltage in vacuum tubes, Journal of the American Institute of Electrical Engineers 47 (3) (1928) 212-213. doi:10.1109/JAIEE. 1928.6538023 .
- [3] W. Schottky, ¨ ber U kalte und warme Elektronenentladungen, Zeitschrift f¨r Physik 14 (1) (1923) 63-106. u doi:10.1007/BF01340034 .
- [4] R. A. Millikan, C. F. Eyring, Laws governing the pulling of electrons out of metals by intense electrical fields, Physical Review 27 (1926) 51-67. doi: 10.1103/PhysRev.27.51 .
- [5] R. A. Millikan, C. C. Lauritsen, Relations of fieldcurrents to thermionic-currents, Proceedings of the National Academy of Sciences 14 (1) (1928) 45-49. doi: 10.1073/pnas.14.1.45 .
- [6] O. W. Richardson, On the extraction of electrons from cold conductors in intense electric fields, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 117 (778) (1928) 719-730. doi:10.1098/rspa.1928. 0031 .
- [7] J. R. Oppenheimer, On the quantum theory of the autoelectric field currents, Proceedings of the National Academy of Sciences 14 (5) (1928) 363-365. doi: 10.1073/pnas.14.5.363 .
- [8] C. C. Lauritsen, Electron emission from metals in intense electron fields, Ph.d. dissertation, California Institute of Technology (1929). doi:10.7907/J2T2-QM49 .
- [9] R. H. Fowler, L. Nordheim, Electron emission in intense electric fields, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 119 (781) (1928) 173-181. doi:10.1098/rspa.1928.0091 .
- [10] T. E. Stern, B. S. Gossling, R. H. Fowler, Further studies in the emission of electrons from cold metals, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character 124 (795) (1929) 699-723. doi:10.1098/rspa.1929. 0147 .
- [11] E. L. Murphy, R. H. Good, Thermionic emission, field emission, and the transition region, Physical Review 102 (1956) 1464-1473. doi:10.1103/PhysRev.102. 1464 .
- [12] E. Z. Engelberg, Y. Ashkenazy, M. Assaf, Stochastic model of breakdown nucleation under intense electric fields, Physical review letters 120 (12) (2018) 124801.
- [13] V. Jansson, E. Baibuz, A. Kyritsakis, S. Vigonski, V. Zadin, S. Parviainen, A. Aabloo, F. Djurabekova, Growth mechanism for nanotips in high electric fields, Nanotechnology 31 (35) (2020) 355301. doi:10.1088/ 1361-6528/ab9327 .
- [14] J. D. Zuber, K. L. Jensen, T. E. Sullivan, An analytical solution for microtip field emission current and effective emission area, Journal of Applied Physics 91 (11) (2002) 9379-9384. doi:10.1063/1.1474596 .
- [15] C. J. Edgcombe, Development of fowler-nordheim theory for a spherical field emitter, Physical Review B 72 (2005) 045420. doi:10.1103/PhysRevB.72.045420 .
- [16] C. J. Edgcombe, N. de Jonge, Deduction of work function of carbon nanotube field emitter by use of curvedsurface theory, Journal of Physics D: Applied Physics 40 (14) (2007) 4123. doi:10.1088/0022-3727/40/14/ 005 .
- [17] A. Kyritsakis, J. P. Xanthakis, Derivation of a generalized fowler-nordheim equation for nanoscopic field-

emitters, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 471 (2174) (2015) 20140811. doi:10.1098/rspa.2014.0811 .

- [18] B. M. Cox, The nature of field emission sites, Journal of Physics D: Applied Physics 8 (17) (1975) 2065. doi: 10.1088/0022-3727/8/17/008 .
- [19] C. Athwal, R. Latham, The effect of the applied field on the energy spectra of electrons field emitted from microscopic sites on broad-area copper electrodes, Physica B+C 104 (1) (1981) 189-195. doi:https://doi.org/ 10.1016/0378-4363(81)90050-4 .
- [20] C. S. Athwal, R. V. Latham, Switching and other nonlinear phenomena associated with prebreakdown electron emission currents, Journal of Physics D: Applied Physics 17 (5) (1984) 1029. doi:10.1088/0022-3727/ 17/5/017 .
- [21] R. E. Hurley, Electrical phenomena at the surface of electrically stressed metal cathodes. i. electroluminescence and breakdown phenomena with medium gap spacings (2-8 mm), Journal of Physics D: Applied Physics 12 (12) (1979) 2229. doi:10.1088/0022-3727/ 12/12/026 .
- [22] R. G. Forbes, J. H. Deane, Reformulation of the standard theory of fowler-nordheim tunnelling and cold field electron emission, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 463 (2087) (2007) 2907-2927. doi:10.1098/rspa. 2007.0030 .
- [23] H. C. Miller, Obtaining numerical values for the elliptic emission functions used in calculating electron emission from surfaces, Journal of Vacuum Science and Technology 17 (3) (1980) 691-695. doi:10.1116/1.570543 .
- [24] R. G. Forbes, Renewing the Mainstream Theory of Field and Thermal Electron Emission, Springer International Publishing, Cham, 2020, pp. 387-447. doi: 10.1007/978-3-030-47291-79 .
- [25] K. L. Jensen, Tutorial: The equations of electron emission and their evaluation, Journal of Applied Physics 135 (11) (2024) 111101. doi:10.1063/5.0196481 .
- [26] D. A. Kirkpatrick, A. Mankofsky, K. T. Tsang, Analysis of field emission from three-dimensional structures, Applied Physics Letters 60 (17) (1992) 2065-2067. doi: 10.1063/1.107090 .
- [27] T. A. de Assis, F. F. Dall'Agnol, Derivation of a currentvoltage-type plot beyond the fowler-nordheim one: The role of the voltage-dependency on the emission area, Journal of Applied Physics 126 (7) (2019).
- [28] H. Tomaschke, D. Alpert, Field Emission from a Multiplicity of Emitters on a Broad-Area Cathode, Journal of Applied Physics 38 (2) (1967) 881-883. doi: 10.1063/1.1709430 .
- [29] Y. Feng, J. P. Verboncoeur, A model for effective field enhancement for Fowler-Nordheim field emission, Physics of Plasmas 12 (10) (2005) 103301. doi:10. 1063/1.2103567 .
- [30] G. A. Farrall, Numerical Analysis of Field Emission and Thermally Enhanced Emission from Broad-Area Electrodes in Vacuum, Journal of Applied Physics 41 (2) (1970) 563-571. doi:10.1063/1.1658714 .
- [31] R. Latham, N. Xu, 'electron pin-holes': the limiting defect for insulating high voltages by vacuum, a basis for new cold cathode electron sources, Vacuum 42 (18) (1991) 1173-1181. doi:https://doi.org/10. 1016/0042-207X(91)90127-5 .
- [32] M. Borghei, D. Velazquez, R. McMullen, G. Latchford,
- B. Riordan, R. Langtry, Impact of direct current conditioning on cathode dark current in high vacuum, in: 2023 30th International Symposium on Discharges and Electrical Insulation in Vacuum (ISDEIV), 2023, pp. 74-77. doi:10.23919/ISDEIV55268.2023.10199791 .
- [33] C. M. Bishop, N. M. Nasrabadi, Pattern recognition and machine learning, Vol. 4, Springer, 2006.
Moein Borghei Robin Langtry
A Data-Driven Model for the Field Emission from Broad-Area Electrodes
Electron emission from cathodes in high field gradients is a quantum
tunneling effect. The 1928 Fowler-Nordheim field emission (FE) equation and the
1956 Murphy-Good FE equation have traditionally been key in describing cold
field emissions, offering estimates for emitters for almost a century.
Nevertheless, applying FE theory in practice is often constrained by the lack
of data on the distribution and geometry of the emission sites. Predictions
become more challenging with an uneven electric field distribution at the
cathode surface. Consequently, FE formulations are frequently calibrated using
current-voltage data after test, limiting their efficacy as true predictive
models.
  This study develops an alternative model for field emission using a
data-driven predictive approach based on (1) vast experimental data, (2)
electrostatic simulations of the cathode surface, and (3) detailed material and
geometry properties, which together overcome these limitations. The objective
of this work is to develop and harness this comprehensive dataset to train a
machine learning model capable of providing precise predictions of the cathode
current in order to further the understanding and application of field emission
phenomena. More than 259 hours of experimental data have been processed to
train and benchmark some of the well-known machine learning models. After two
stages of optimization, a coefficient of determination >98% is achieved in the
prediction total field emission current using ensemble models.