Microconical Structure Formation and Field Emission From Atomically Heterogeneous Surfaces Created by Microwave Plasma–Based Low-Energy Ion Beams

We report the formation of self-organized microconical arrays on copper surface when exposed to high flux (5.4 × 1015 cm−2 s−1) of 2 keV argon ion beams at normal incidence. The created microconical arrays are explored for field emission properties. The surface morphologies are investigated by scanning electron microscopy and atomic force microscopy. The local work function variation is analyzed by Kelvin probe force microscopy, and the argon content in the irradiated layer is measured with X-ray Photoelectron Spectroscopy. The average aspect ratio (base width/height) of microstructures for individual irradiated samples is found to increase from 0.7 to 1.5 with a decrease in ion fluence. The ion concentration is highest (3.89 % ) for a fluence of 4.7 × 1018 cm−2, which asserts the formation of atomically heterogeneous surface due to subsurface ion implantation. An enhancement in the field emission properties of the argon ion–treated copper substrates at a fluence of 4.7 × 1018 cm−2 with a low turn-on voltage of 2.33 kV and with electron emission current 0.5 nA has been observed. From the Fowler–Nordheim equations, the field enhancement factor is calculated to be 5,561 for pristine copper, which gets enhanced by a factor of 2–8 times for irradiated substrates. A parametric model is considered, by taking into account the modified local work function caused due to structural undulations of the microstructures and presence of implanted argon ions, for explaining the experimental results on the field enhancement factor and emission current.


INTRODUCTION
Physical texturing on material surface by bombardment of energetic ion species is a contaminationfree and controlled technique of generating micro-or nanostructures of required geometry, composition, and concentration [1][2][3]. There are extensive experimental reports on growth of variety of self-organized patterns like fuzz, craters, dots, and ripples on carbon nanotubes, graphene sheets, and metallic compounds like ZnO, semiconductors like Si and Ge, and high-Z metals like W, Mo, and Ta, including Pd when irradiated with inert gas ions such as He and Ar [4][5][6][7]. The effective surface area gets enhanced due to the production of surface patterns, which in turn act as naturally selected micro-or nanofield electron emitters. Electron emitters with high current density and low turn-on field is desirable for their practical applications in cold electron sources like field emission cathodes (vacuum electronic devices), flat panel displays, gas sensors, magnetic storage devices, and microwave amplifiers including catalytic reactors [8][9][10][11][12]. The apex shape, size, distribution, and the local work function variation of emitters contribute to modify the field emission parameters. Hence, to improve field emission several methods like increasing carrier concentration by ion doping [13], modifying the structures by plasma treatment [14,15], or creating modified structures by ion bombardment [16] are employed. There are some reports that suggest deposition of metallic impurity atoms of low sputter yield during ion bombardment of metals, leading to the formation of disordered arrays of conical protrusions known as sputter cones [17][18][19]. On the other hand, in a wide range of devices involving metal surfaces exposed to high electric fields like in plasma facing components of fusion reactors and particle accelerators, field emission current leads to initiation of catastrophic vacuum breakdown (dark current) [20]. Hence, various fundamental studies are carried out to understand the intrinsic mechanism of electron emission under such a high applied field in metallic substrates [21,22]. Apart from topographical modifications, the surface also becomes heterogeneous in atomic length scale due to the introduction of irradiated ions in the interstitial layers (nm) [23]. The development of surface structures with compositional heterogeneities at the atomic level because of ion irradiation is an important field of research in material processing and thin film studies because of the controllability in the field emission properties that it can offer. The electrical, surface, and optical properties of such atomically heterogeneous metallic systems with patterned arrays of micro-or nanometric dimensions produced by irradiation of high fluence of energetic inert ions have been investigated earlier [23][24][25].
The present study focuses on the formation of naturally selected metallic field emission arrays. These microprotruding arrays are advantageous for rudimentary studies on various aspects of field emission since generation of electrons (leading to high emission current) is easier because of their low resistivity, high refractivity, and thermal effects [26,27]. Copper being an earth-abundant element with similar electronic configuration (σ Cu 5.96 × 10 7 Sm −1 at 20°C) of inert metals like gold and silver is considered for generating microstructures by low-energy inert gas ion beam. These modified copper substrates have variety of applications as sensing materials, in aeronautical industries, and as photocatalytic materials [28][29][30]. However, fundamental field emission experimental studies on such structurally modified plasma ion-irradiated copper substrate have remained unexplored.
In this article, the mechanism of formation of self-organized structures by low-energy ion irradiation and their structural contribution toward enhancement of emission current is explored. The conical structural growths on copper substrate in the absence of secondary metallic impurities are comparatively new and their electron emission response to strong electric fields has not been reported previously. This work presents feasibility of creation of microconical arrays at normal incidence. The experiment of irradiating samples for growth formation is carried over a wide range of fluences. Thereafter, from selectively chosen irradiated substrates, the field emission experiment is carried out individually on the samples.

DISCUSSION ON EXISTING THEORIES ON ION-IRRADIATED STRUCTURE FORMATION AND FIELD EMISSION STUDIES
A variety of theoretical models describe the growth phenomena as competitive processes between orientation-dependent surface diffusion process and sputtering. The Ehrlich-Schwoebel model describes that surface roughening is initiated by temperature-sensitive surface diffusion of adatoms in convex curvatures [31]. The nanoscale bubbles below the material surface coalesce to form craters and nanofibrils near the surface, which to some extent are explained by the surface adatom diffusion model by Martynenko and Nagel [32]. The Bradley-Harper surface instability theory [33] (or some extensions of this theory [34]) explains ripple formation (seen prominently on ion irradiated Si substrate) using the concept of preferential sputtering of concave surfaces with deposition on convex surfaces, especially for irradiation at an angle. The conical protrusions with impurity metal atom resting at the apex of the cones are believed to be formed on the surface to shield the underlying surface from incoming flux [35,36]. Robinson and Rossnagel's thermal diffusion model considers the formation and redeposition of impurity clusters as a result of thermal surface diffusion leading to migration of atoms in a random walk pattern [37]. Davidovitch et al. showed that the mass distribution along with erosive surface effects changes the surface height to a considerable degree [38]. Bradley tried to improve upon these models by considering early time dynamics of the initially planar solid surface seeded with impurities during ionic bombardment [39]. Perkinson's crater function theory incorporates both erosion and mass redistribution effects, neglecting the effects of local curvature at high incidence angles [40].
In the present experiment, no external impurities with differential sputtering yield strength are added. To understand the conical protrusion formation, we consider that at the impact point there is a volumetric removal of material and subsequent deposition after rearrangement, leading to the agglomeration of substrate atoms acting as seeding agents to aid in the growth process. Moreover, the theory for sputter cone formation in the nucleation and growth regime is highly complex and nonlinear and requires profound experimental evidence in order to properly correlate various conditions of the growth pattern formation with varying ion beam parameters and is beyond the scope of this study.
The field emission process (or Fowler-Nordheim tunneling) is quantum tunneling of electrons through electron rich surface-vacuum interface in the presence of a high electric field ( > 1 GV/m). The Fowler-Nordheim (FN) equation assumes a smooth planar metal surface with a given work function, free-electron model (Sommerfeld), and a semi-classical approximation (WKB) for the tunneling probability. The current density J from the Fowler-Nordheim equation is given by the following equation [41,42]: where E is the applied electric field, ϕ is the bulk work function, A and B are constants, and A 1.54 × 10 −6 AV −2 eV, and B 6.83 × 10 7 V cm −1 eV −3/2 . The field enhancement factor (β) is a morphology and work function-dependent parameter, affecting the local field around micro (nano)-scale protrusions on the metal surface. As seen from the FN equation (Eq. 1), the field emission current depends on ϕ and β which in turn depends on the aspect ratio and added ion impurities which can modify the work function locally. The curvature correction on field emitter tips has been recently studied in order to incorporate radius of curvature as an important parameter on which the field enhancement depends upon, reported independently by Kyritsakis and Biswas [43,44]. Their theoretical results are found to be consistent for tips having radius of curvature within 5-10 nm. The Fowler-Nordheim formula has been modified by considering the local spherical approximation of emitter surface, considering the effect on the image potential due to interactive electrons at a certain distance from the emitter [45]. Density functional theory (DFT) calculations by Djurabekova, Ruzibaev et. al., on copper (100) surface have shown a significant effect on the work function due to microstructures [46]. An ab initio computational method has been developed to calculate the field-emitted current from defective metal surfaces by combining the density functional theory and electronic structure calculations with quantum transport methods on copper (111) surface by Toijala [47]. However, from this model, it is difficult to justify the effective high value of the field enhancement factor in the field emission characteristics of a metal surface, as observed in their as well as our present experiments.

EXPERIMENTAL DETAILS
The experimental set-up, shown in Figure 1, comprises of two regions the plasma production chamber and the ion extraction chamber. Microwaves of 2.45 GHz are launched from a magnetron (Alter-TMA20) connected with its power supply and controller (Richardson PM 740). The waves are passed through a metallic rectangular waveguide (WR 340) and a quartz window into the plasma chamber, and the power is maintained at 180 W for continuous mode of operation. A water-cooled isolator (ISO) protects the magnetron from the reflected power. The directional coupler (DC) measures the forward and reflected power of the microwaves. The triple stub tuner (TST) is used for matching the impedance of the magnetron to the plasma, so that maximum power can be transferred to the plasma. The plasma production region consists of a vacuum chamber (VC) of length 50 cm and diameter 20 cm and a magnetic multicusp (octupole) which is placed inside the VC for confinement of the plasma [24]. The ports in the VC are used for gas inlet, probe diagnostics, pressure gages, and venting. The amount of Ar gas entering VC is controlled by a mass flow controller (MKS1179A). The base pressure of VC is kept around ∼ 1.33 × 10 −6 mbar using a turbomolecular pump (Varian 301 Navigator) and a rotary pump, while the working pressure is kept at ∼ 1.33 × 10 −4 mbar. In the ion extraction region, ion beams are extracted from the microwave-induced plasma by employing metallic electrodes. The plasma electrode (PLE), which is the electrode that faces the plasma, is held at the ground potential and has a small aperture diameter of 1 cm. The collector plate on which samples are mounted is floated to a negative high voltage (for the extraction of ions from the plasma) through a vacuum Frontiers in Physics | www.frontiersin.org electrical feed-through. At the other end, the collector is attached to an XY stage manipulator. The curved potential lines which bend through the PLE aperture and penetrate into the plasma sheath region act as a lens for the beams and help them focus onto the substrate held on the collector plate [48]. All ions arrive at the PLE with Bohm velocity and have energy (eV p ) corresponding to the difference between the plasma potential (few tens of volts) and the ground potential. From the PLE, they are extracted with the potential applied at the collector (V c ∼ kV). Therefore, all ions are expected to be more or less monoenergetic (eV p ≪ eV c ), as the extraction is made through vacuum and there are no ion-neutral collisions.
The samples are copper sheets (purity > 99.99 % ) of dimensions 1 × 1 × 0.44 cm 3 . In order to remove the hydroxide layer, the substrates are immersed in a prepared solution of 2MOL/l HNO 3 for 45 s and then ultrasonically cleaned by acetone and ethanol. The samples are inserted in an ion extraction chamber and are exposed to energetic argon ions (2 keV) with a flux of 5.431 × 10 15 cm −2 s −1 at normal incidence angle (as shown in the inset image of Figure 1). It has been found in earlier works that singly charged ions are dominant in the plasma [49]. Experiments on surface modification is carried out at fluences of 1.3 × 10 18 cm −2 -4.7 × 10 18 cm −2 . The percentage error on the determined fluence, especially when using plasma to generate the ion beam, is calculated to be ∼ 4.2%. The beam diameter is kept at 1 cm to maintain a uniform irradiation of the exposed part of the sample [23]. The samples after irradiation are kept in vacuum desiccators to prevent any environmental contamination.
The surface morphologies of pristine and irradiated samples are investigated by Zeiss scanning electron microscopy (SEM, Model-SUPRA40VP). Atomic force microscopy (AFM, Model-ASYELEC-01, Asylum Research) provides the evidence of change in surface roughness due to microstructure formation. The MFP: 3D classic atomic force microscopy (AFM) measurements are performed in air with n-type silicon cantilever coated with Ti-Ir. AFM and Kelvin probe force microscopy (KPFM) scans an area of 25 μm 2 at different locations on the samples to account for all types of irregularities. The KPFM provides the local work function of the surface scanned.
X-ray photoelectron spectroscopy (XPS, Model-5000 Versa Probe II spectrometer, FEI Inc.) technique is employed for elemental analysis. The X-ray beam diameter is 200 μm and the average area over which the sample is probed is 1 mm 2 . The achieved energy resolution is ΔE ≤ 0.5 eV for Al K α excitation (1486.6 eV) in ultrahigh vacuum (10 −9 mbar). XPS spectra are calibrated using the reference energy value of carbon core level E (C 1s ) 284.8 eV. The implanted argon concentration (%) up to second decimal place is calculated from the survey spectral XPS analysis of the peaks using the Casa XPS software.
The digital picture and schematic diagram of the field emission setup is shown in Figure 2. The experiments are performed on both pristine and irradiated samples with developed microcones, at room temperature. The experiments and measurements which are performed at room temperature implied that no additional heating of the substrate is done externally. The substrate is copper, which is a good conductor and the heat generated gets transported away by conduction through the metallic rod handle and radiation. Moreover, we have also performed the experiment by breaking down the time period of irradiation in slots of four for 60 s each and still have observed no significant change in pattern formation. So, we have taken into account the fluence variation as the key factor in generating dense structures. The field emission experimental setup is a parallel plate configuration with gap distance of 2 mm, kept in high vacuum chamber having a base pressure of 1.33 × 10 −6 mbar. The anode is mirror-polished copper and the cathode is the substrate. The digital picture is the end section of a high vacuum-focused ion beam system, which is employed for the field emission measurements [50]. Figure 3 shows evenly distributed microstructures as observed from SEM analysis. In Figure 3A, we can see pristine copper with minute bright spots. The contrasting bright spots are initial random roughness generated due to HNO 3 etching to remove the copper hydroxide layer which acts as a precursor of nucleation sites for stabilizing the growth process. Copper substrate exhibits high electrical conductivity, thermal conductivity, malleability, and ductility, which aids in easier structure formation for comparative analysis of field emission effects for a chosen range of fluences, when exposed to high electric field. In Figure 3B at a fluence of 1.3 × 10 18 cm −2 , the structures are mound-like with well-defined rim. During ion bombardment, the surface gets destabilized by implanted ionic impurities for certain critical impurity flux leading to a microscale mound formation. It is observed from Figure 3C that the mound boundary becomes well defined with the onset  Figure 3D to Figure 3F, we see well-defined micro hillocks formation with individual base boundaries and certain protrusion height. It is observed that, with the increase in fluence, the structures become uniformly distributed with a prominent conical shape formation. At the maximum fluence of 4.7 × 10 18 cm −2 , as seen in Figure 3F, we observe production of self-organized microcones which are densely distributed. The conical mounds slowly become prominent hillocks with increase in ion beam fluence. The protrusion height variation is between 2 and 4 μm having a relatively constant base width of ∼ 3 μm. For fluences higher than 5 × 10 18 cm −2 , the sharpness of these structures reduces with deformation of the peaks. At higher fluences, sputtering erodes away these structures depending upon the strength of the ion beam and tenure of exposure [51]. Thus, optimization of fluences is required as it increases the number of emitters per unit area due to the production of multiple nucleation sites. In Figure 4, the AFM and KPFM results are reported. The surface roughness variation of pristine and two of the exposed substrates at fluences of 3.2 × 10 18 cm −2 and 4.7 × 10 18 cm −2 , respectively, are shown. Their corresponding work function distribution with respect to the local surrounding is obtained from KPFM analysis. In Figure 4, as one moves from left to right, the variation in surface topography and work function of the substrates show well-defined structure formation with an increase in fluence, similar to the trend obtained previously in the SEM analysis in Figure 3. KPFM is an atomic force microscope-based technique aiding in determining the changes in the work function and the surface potential in a nondestructive manner. In KPFM, the distributed contact potential difference (V C ) between the scanning tip and the sample is calculated as follows:

EXPERIMENTAL RESULTS AND DISCUSSIONS
where ϕ tip ( ∼ 5.0 eV) and ϕ sample are the work functions of the tip and sample, respectively, while e is the electronic charge. The work function of the sample is calculated from the local potential change due to the presence of undulations on the surface. The local surface potential is averaged over a finite effective area of the surface involved in the interaction with the tip. In KPFM, the sample potential thus contributes to a bigger volume of the tip than its nanometric apex. Hence, an effective work function is calculated over the region where the stylus is scanning. Since the samples are irradiated uniformly, the overall work function has decreased. In the valley regions surrounding the microconical tips, the work function is observed to be 4.85 eV, while in the tip region, there is a significant decrease to 4.69 eV for fluence of 4.7 × 10 18 cm −2 . For the pristine sample, the variation in work function is over a short range from 4.97 to 5.09 eV, which is because the pristine substrate has no implanted argon ions and more or less has uniform height distribution ( < 100 nm) as seen in Figure 4A. Figure 5A shows high resolution spectra of argon with peaks at Ar 2p 1/2 ( ∼ 244.1 eV) and Ar 2p 3/2 ( ∼ 241.8 eV) determined using XPS. A Shirley background fit is done to fit the XPS measurements. The Savitzky-Golay smoothing is done for multipeak analysis which is a standard smoothing procedure done to fit the XPS data. The peaks are resolved at ΔE 2.3 eV, indicating inertness of Ar embedded on the metal surface. After argon irradiation, the intensities of the two peaks (Ar 2p 1/2 and Ar 2p 3/2 ) have increased, as seen from their areal coverage. For a flux of 4.7 × 10 18 cm −2 the argon concentration is the highest ∼ 3.89% and for 2.5 × 10 18 cm −2 the argon concentration is ∼ 2.72%. The implantation depth of ions and their concentration can be varied by adjusting the ion beam parameters such as energy of ions, fluence, gas pressure, and wave power, including projectile target geometry. The results confirm argon ion implantation on irradiated copper surface by introducing atomic heterogeneity on the exposed surface and contributing significantly to the variation in the local work function. The pristine sample shows zero argon ion concentration from XPS analysis, which is consistent with the fact that the argon ions are not adsorbed on the exposed surface from environmental contamination, when exposed to external conditions while loading and unloading the samples. The XPS survey spectra are shown in Figure 5B. The as-prepared samples are preserved in a clean vacuum desiccator to avoid any exposure in the ambient air. The compositional studies revealed that an appreciable amount of oxygen and carbon is present in the nearsurface atomic layers of the pristine and irradiated samples. Although the as-prepared samples are preserved in a clean vacuum desiccator to avoid any exposure in the ambient air, there could be some possibilities to absorb some atmospheric gases while transferring for measurements. The irradiated samples have a considerable amount of argon implanted whose percentage implantation is under our experimental control.
In Figure 6, the decrease in variation of work function distributed over the hillock (microconical tips) and valley (regions surrounding the microconical tips) regions, respectively, with the concentration of the implanted argon ions are plotted. The implanted argon concentration should increase with increase in time of exposure of the ion flux as the number of argon ions per unit area increases. The results obtained are consistent with this statement as they can be seen in Figures 5,  6. Argon implantation influence on the work function of metallic substrate is a comparatively new topic. A recent article aims at investigating the effect of particle irradiation on the work function in simple metals like Al and Zn using Ar + ions at room temperature with 30 keV. The article uses similar KPFM technique and delves deep into understanding the theoretical aspects of the work function variation due to argon implantation [52].
The field emission experiment is performed using the micro conical structured substrates. The Fowler-Nordheim (FN) relation can be written in terms of the experimentally measured emission current (I) and applied voltage (V) as follows [53]: where ϕ is the bulk work function (4.99 eV), d is the gap distance, S is the effective area, and A and B are constants as defined earlier in Eq. 1.
The characteristics of emission current versus applied voltage plot are shown in Figure 7A and the corresponding FN plot in Figure 7B. It is evident from the experimental data tabulated in Table 1 that there is a decrease in the work function with an FIGURE 5 | High-resolution X-ray photoelectron spectroscopy (XPS) spectra in (A) and X-ray photoelectron spectroscopy (XPS) survey spectra in (B) for two chosen argon-bombarded copper surface compared to pristine. increase in fluence. This is observed due to the formation of the atomically heterogeneous surface with increased production of uniformly aligned distinct micro conical tips as seen in Figure 3. Thus, the emission of electrons from these irradiated surfaces at a higher fluence is easier at a comparatively low applied voltage (turn-on voltage), as seen in Figure 7A. The turn-on voltage for substrate irradiated with a fluence of 4.7 × 10 18 cm −2 is 2.33 kV, is significantly lower than that of the pristine sample having a turnon voltage of 6.51 kV. The FN plot in Figure 7B is fitted with a straight line and the field enhancement factor (β) is extracted from the slope.
The field enhancement factor is a morphology-and work function-dependent parameter, affecting the local field around micro (nano)-scale protrusions on the metal surface. There is a significant rise in steepness of the linear fitting on experimentally observed curves (and consequently the enhancement factor), with an increase in fluence.
As previously stated, the Fowler-Nordheim equation assumes a smooth planar metal surface with a given work function. But the observed work function gets modified due to both protrusion formation and subsurface implantation of argon ions causing a considerable change in the field emission properties of any ion irradiated substrate. Thus, assuming the bulk work function of substrate in calculating the field emission parameters particularly for irradiated samples may produce erroneous results. To explain our experimental results on field emission, we consider a parametric model to describe the field emission properties for atomically heterogeneous surfaces where the work function is taken to vary across the area encompassing both hilly and valley regions, implanted impurities, and structural deformities. We consider the effective work function to be a linear combination of the work function observed over the undulated surfaces. The parameters γ and δ are the effective proportions in which hillocks and valley areas are distributed for a particular substrate observed using images obtained after analysis from KPFM and SEM results. Using ImageJ software, the hilly regions obtained in digital scan are demarcated and isolated. The areal percentage is calculated from the pixel density and then, incorporated in their corresponding color coded regions to obtain the percentage of work function and aspect ratio contributions. ImageJ software is well-known software used to isolate digital image of cells in soft matter to calculate their areal percentage using pixel density difference of the image taken [54]. The effective work function thus becomes ϕ eff , and is given as follows: which is employed in Eq. 3 in the place of ϕ and the β obtained is labeled as β * which is obtained directly from the slope as discussed earlier. ϕ 1 and ϕ 2 are area-averaged work functions of hilly and valley regions, respectively.
In Figure 8A, the work function variation of both hillock-like and valley regions with aspect ratio is shown. The aspect ratio decreases with the increase in fluence with the formation of prominent and well-defined hillocks and valleys. The valleys are basically upturned conics where there is significant contribution from the nonuniform side of the slopes. Thus, the work function decreases for both hilly and valley regions leading to the extraction of electrons with ease and increase in the enhancement factor. In Figure 8B, the enhancement factor is plotted along with the aspect ratio as a function of ion beam fluence. The electron emission from these structures with low aspect ratio is relatively easy because of their pointed geometry and small radius of curvature. The generation of naturally selected absolute uniform sets of sub-microconical metallic arrays having each tip of similar height and base diameter is difficult. Hence, we have calculated the average aspect ratio for each of the irradiated samples. The aspect ratio is defined as ratio of the base width of the conical structure and its corresponding height, averaged over a number of such microprotrusions evaluated from the SEM analysis. The aspect ratio is found to be in the range of 0.7-1.5, implying low aspect ratio for pointed microstructures. The electrons can easily tunnel through such configuration with additional aid from the local surface charge developed in the vicinity of such structures. The presence of implanted argon ions may also influence the surface charge production significantly by causing variation in the enhancement factor. The modified enhancement factor β * obtained by employing modified work function (ϕ eff ) is compared with enhancement factor β and evaluated considering bulk work function ϕ in Table 2. The overall high values of the enhancement factor are due to the increase of the field emission current due to atomic-level surface defects such as implantation of atoms and locally attached surface adatoms apart from the sharp surface deformities. It may be understood that the change occurring at the valley regions is just because of the subsurface argon atoms which perhaps modify the density of states and lead to a lowering of the work function. The change happening at the hilly regions (conical tips) is a combination of the effect of the implanted argon impurities and the field enhancement. These contributions are not accounted for in the original Fowler-Nordheim equations. Moreover, vacuum electronics are sensitive to micro-and nanometric deformities. Hence, more precise evaluation of the enhancement factor is beneficial over the enhancement factor which has already been evaluated purely considering the bulk work function. Thus, the parametric modeling focuses on the need for considering individual topographical contribution of ion-irradiated surfaces as these micrometric surface effects infact influence the effectivity of field emitters.

SUMMARY
In conclusion, the present study revisits sputter cones like structure with well-defined arrays of conical hillocks of micrometer dimension, challenging their requirements for seeded impurity atoms along with exploring their possibilities as promising field emitters. The growth mechanism of microconical arrays needs to be explored with a proper support from the substrate-oriented and modified existing theories in a channelized manner. The high enhancement factor of pristine is due to copper being a highly conductive metal, but the subsequent increase in the enhancement factor with field emission current for irradiated surfaces is because of the formation of atomically heterogeneous system (by low energy plasma ion irradiation) with multiple microprotrusions. The field emission results have shown almost eight times increase in the field enhancement factor for high fluence irradiation (4.7 × 10 18 cm −2 ).
The local variation of work function is mainly due to the evolution of the structural deformities and implantation of argon ions as observed from the employed experimental techniques and reported in this article. In our parametric model, we have tried to tune the effects of the work function of modified atomically heterogeneous surface on their field emission parameters. From the results presented so far, we could satisfy the objective of presenting these newly formed naturally selected structures and side by side show the need for modified Fowler-Nordheim theory as a slight change in the work function brought about by the deformed atomically heterogeneous system produces significant change in the experimentally obtained field emission parameters.
Hence, there is a requirement for a relook at the existing FN theory for incorporating the effects of both the work function variation due to atomically heterogeneous surface and the corresponding field enhancement at the microconical tips. As a future work, it would be helpful to distinguish between the contribution of surface structure formation and heterogeneity introduction separately on the increase in field emission parameters. Experimental techniques aimed at masking the geometric effects of micro protrusions on field enhancement factor so that the effect of atomically heterogeneous surface layer influencing the field emission parameters that is discernible may be devised. A theoretical model based on understanding of the individual contribution of the parameters obtained from the experimental results and analysis, where a statistical distribution of arrays of microconics may be considered on a pure copper substrate. These microconical arrays with randomly distributed implanted inert gas ions will then represent the irradiated substrate for a comparative study with the experimentally obtained results.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article, and further inquiries can be directed to the corresponding authors.

AUTHOR CONTRIBUTIONS
JM and SB: designed and carried out the experiments, analyzed the experimental results, and wrote the manuscript. All authors contributed to the article and approved the submitted version.