The transition from field emission to collisional space-charge limited current with nonzero initial velocity

Multiple electron emission mechanisms often contribute in electron devices, motivating theoretical studies characterizing the transitions between them. Previous studies unified thermionic and field emission, defined by the Richardson-Laue-Dushman (RLD) and Fowler–Nordheim (FN) equations, respectively, with the Child-Langmuir (CL) law for vacuum space-charge limited current (SCLC); another study unified FN and CL with the Mott-Gurney (MG) law for collisional SCLC. However, thermionic emission, which introduces a nonzero injection velocity, may also occur in gas, motivating this analysis to unify RLD, FN, CL, and MG. We exactly calculate the current density as a function of applied voltage over a range of injection velocity (i.e., temperature), mobility, and gap distance. This exact solution approaches RLD, FN, and generalized CL (GCL) and MG (GMG) for nonzero injection velocity under appropriate limits. For nonzero initial velocity, GMG approaches zero for sufficiently small applied voltage and mobility, making these gaps always space-charge limited by either GMG at low voltage or GCL at high voltage. The third-order nexus between FN, GMG, and GCL changes negligibly from the zero initial velocity calculation over ten orders of magnitude of applied voltage. These results provide a closed form solution for GMG and guidance on thermionic emission in a collisional gap.

www.nature.com/scientificreports/where E s is the electric field at the cathode, y = −1 √ 4QeE s where Q = e 2 /(16πǫ 0 ) , A FN = e 3 / 16π 2 ℏ� and B FN = 4 √ 2m� 3 /(3ℏe) are FN constants, is the electrode work function, and ℏ is the reduced Planck's constant.Table 1 summarizes the values of key physical parameters used in these equations and throughout this paper.
Recasting these equations in nondimensional variables yielded a universal (true for any diode geometry) set of equations and matching the resulting asymptotic solutions indicated the transition from a field emitting diode to a space-charge limited one 8 .Note that this matching invalidates the assumptions inherent in the fundamental equations since CL is commonly derived for E s = 0 , while FN requires nonzero E s = V /D .Thus, characterizing electron emission near this "nexus" necessitates solving the full solution that accounts for all relevant physics.Similar asymptotic approaches have been applied to characterize other emission and circuit phenomena.By introducing electron mobility µ into the electron force law to account for collisions, Benilov assessed the transi- tion from CL to SCLC with collisions 16 , or the Mott-Gurney (MG) law, given by 17 where µ is the electron mobility in the gas.Additional applications of this approach to assess field emission, SCLC, and resistive dissipation include the transitions between CL, FN, and Ohm's law 18 ; CL, FN, and MG 19 ; and CL, FN, MG, and Ohm's law 20 .This approach has also been applied to include the enhancement in quantum SCLC at small gaps with the transition to field emission driven microdischarges at large gaps 21 and the transitions between FN and SCLC in liquids and during the transition in phase from liquid to gas 22 .This approach of developing full analytic theories coupling relevant electron emission sources and assessing the transition between dominant mechanisms, referred to more broadly as "nexus theory" 4 , may be extended to other electron emission source mechanisms.For instance, Darr et al. 23 replaced FN with the general thermal-field (GTF) emission equation 2,[24][25][26] , which couples field and thermionic emission, to derive exact and asymptotic solutions linking FN, CL, and the Richardson-Laue-Dushman (RLD) solution for thermionic emission, given by 27 where T is the cathode temperature, k B is the Boltzmann constant, and A RLD = em/ 2π 2 ℏ 3 .Darr et al. extended previous studies examining the transition from RLD to CL using Miram curves [28][29][30][31][32][33] by incorporating FN and showing that an increasing contribution of field emission could soften the characteristic "knee" during the transition to CL. Lang et al. 34 replaced the GTF equation with the general thermo-field photoemission (GTFP) (2) equation 24,25 to include the transition to the Fowler-DuBridge (FD) equation for photoemission 2,35 .In addition to providing a detailed step-by-step process for developing nexus phase space plots showing the device conditions necessary for the transitions between FD, RLD, FN, MG, CL, and Ohm's law, Lang et al. derived exact solutions for the current density as a function of voltage and demonstrated the transitions to FD, RLD, FN, and CL under appropriate asymptotic conditions of temperature, gap distance, voltage, and laser frequency 34 .
Nexus theory may also be applied to assess variations in electron emission regime for rough electrodes.A recent study applied nexus theory to a fractional dimensional model to derive the conditions for transitions between FN, MG, and CL for different conditions of surface roughness 36 .Such a model may be extended to the condition here combining RLD with FN, MG, and CL by following the steps outlined in that study; however, it would also be necessary to extend the GTF representation of combined electron emission to fractional dimensions to derive the solution for the exact current density in the diode.
Although a given emitter may realistically exhibit multiple mechanisms due to field and temperature nonuniformities, the present study considers a one-dimensional (1D) system assuming all electron emission behavior along the cathode follows a common mechanism (or combination of mechanisms).Using 1D nexus theory simplifies the assessment of emitter behavior, elucidating inter-mechanism transitions based on operating system conditions.
While several studies have applied nexus theory to assess numerous transitions, none have explicitly considered the transition between RLD and MG.While this transition has been examined using nexus theory by matching the RLD and MG 23 , an exact solution fully coupling collisions and the GTF or GTFP source has not yet been conducted.While seemingly of academic interest, this assessment has practical implications in device physics.Although many vacuum electronics devices use cold cathodes (field emission with the electrons emitted with zero velocity) 35,37,38 , many use thermionic cathodes (electrons emitted with nonzero velocity) [38][39][40] .This motivated the earlier study linking RLD with FN and CL 23 .Less often considered theoretically, although motivating the previous study linking MG with FN and CL 19 , is that vacuum devices often do not operate in perfect vacuum.Pressures of 10 -5 Torr can dramatically degrade field emission 41 and the failure rate of emitters increases with increasing pressure 42 .Alternatively, at atmospheric pressure for microscale gaps, strong electric fields strip electrons off the cathode by field emission to reduce breakdown with decreasing gap size 5 .Reducing gap size to nanoscale causes space-charge to begin to limit the current prior to breakdown 7 .Since thermionic emission may also play a role in micro-or nano-scale breakdown 6 , this physics is also relevant for higher pressures; one example being high altitude discharge.For instance, when attempting to define "high pressure" and "low pressure" cutoffs for choosing the correct model to describe a device, it can be helpful to know the "midpoint pressure" as defined by nexus theory involving MG.This could be especially useful in thermionic diodes, to define a threshold for temperature-induced outgassing that may change device operation [43][44][45][46][47] .Hanquist et al. addressed high heating rates on hypersonic vehicles and developed a model combining SCLC with thermionically emitted electrons at the surface to include transpiration cooling 48 .
The other challenge is that when examining the nexuses of SCLC with RLD, one must also consider the implications of nonzero injection velocity on SCLC.The most widely-cited derivation of 1D, planar SCLC with nonzero electron injection velocity was derived by Jaffé [49][50][51] as where v 0 is electron injection velocity.Liu and Dougal derived an alternate solution, referred to as the bifurcation current density J B , given by 52 In (8), the electrons have zero velocity at the virtual cathode.In (7), the electric potential is shallower since the electrons maintain a nonzero velocity at the virtual cathode.From the asymptotic analysis of the transition from RLD to CL, Darr et al. derived a generalized CL law (GCL) J GCL as 23 which behaves identically to (7) in the limit of mv 2 0 (2eV ) −1 ≪ 1 .However, no readily available solution for an analogous general MG (GMG) solution for nonzero injection velocity exists to serve as an appropriate nexus equation for a thermionic cathode in a non-vacuum environment.
Thus, this paper will develop a first-principles based linkage of FN, RLD, GMG, and GCL to assess this behavior considering a 1D system.We first derive the theory starting from the electron force law coupled with the GTF and obtain the relevant asymptotic solution for GMG.We next apply the resulting asymptotic theories and exact solutions to assess nexuses under various conditions before making concluding remarks.

Results
Derivation.We consider a 1D, planar diode containing neutral gas with electron mobility µ , the cathode at x = 0 held at electric potential φ = 0 , and the anode at x = D held at φ = V with respect to the cathode.We assume the electron is emitted from the cathode at x(0) = 0 with initial velocity v(0) = v 0 and initial accelera- tion a(0) = eE s /m .Combining Poisson's equation with continuity, given by J = en e v , yields where n e , v, and J are electron number density, electron velocity, and current density, respectively.Assuming that the current is emitted due to a combination of thermionic and field emission, we define J using the GTF relation, , and β F , s , and N(n, s) are functions of F and T (see Supplementary Information) 23 .The force on an electron is given by The first term on the right-hand-side of (11) represents the force on the electron due to the electric field, while the second represents a friction term introduced by collisions.To reduce parameters and facilitate analysis, we nondimensionalize ( 10) and ( 11) by defining where the bars represent dimensionless parameters, terms with subscript 0 are scaling terms, and the FN coefficients are given by A FN = e 3 /(16π 2 �) and B FN = 4 √ 2m� 3 /(3 e) .Substituting ( 12) into ( 10) and ( 11) yields and respectively.Equations ( 13) and ( 14) are universal since all material dependence has been removed through nondimensionalization.
Differentiating (14) with respect to x , considering v = dx/dt to change variables, and combining with (13) gives Solving (15) for velocity using the initial conditions defined previously yields where v 0 = √ (k B T)/m e 23 , which may be written non-dimensionally as v 0 = T/26.1145 .Integrating ( 16) gives electron position as At low mobility (high pressure), exp(−t/µ) ≈ 0 , simplifying (16) and (17) to v(t) ≈ µ E + Jt + v 0 and x(t) ≈ µ Et + Jt 2 /2 + v 0 t , respectively.The critical current density can then be obtained by considering the condition x(τ ) = D .This gives the transit time as τ = χ E/J , with χ simplifying to for µ ≪ 1.
Instead of defining V using an energy balance equation 8 , we integrate (14) with respect to x and change variables to t to obtain 19 (10) www.nature.com/scientificreports/where v(t) is defined by (16).Because µ ≪ 1 , the second term on the right-hand side of ( 19) dominates.Substi- tuting the simplified velocity function, v(t) ≈ µ E + Jt + v 0 , and τ into (19) gives Neglecting higher order terms of µ , considering V ≫ 1, and incorporating the definition of χ from ( 18) simplifies (20) to give the general MG (GMG) equation for nonzero initial velocity as For v 0 = 0 , or T = 0 , J GMG simplifies to (5).Assuming µ ≫ 1 (low pressure) eliminates the collisional terms and simplifies ( 16) and ( 17) to v(t) ≈ (Jt 2 )/2 + Et and x(t) ≈ (Jt 3 )/6 + Et 2 /2 , respectively.Taking those equations and considering V ≫ 1 recovers the general Child-Langmuir (GCL) function 23 , given by [cf.(9)] FN defines the limit in which field emission dominates, and is recovered for V ≪ 1 to give RLD predicts the limit when thermionic emission dominates, T ≫ E , and is given by Theoretical analysis and results.Figure 1 compares the full solution from (19) with the asymptotes defined in ( 21), ( 22), ( 23) and (24) for GMG, GCL, FN, and RLD, respectively, for D = 250 .For µ = 5 and T = 0.04 , Fig. 1a shows that electron emission transitions from RLD to FN to GMG to GCL with increasing V .Reducing T to 0.002 for µ = 5 eliminates the RLD regime in Fig. 1b.In general, for sufficiently low µ in Fig. 1a  and b, electron emission first transitions to GMG before reaching GCL. Figure 1c shows that increasing µ to 100 for T = 0.04 causes electron emission to transition from RLD to FN to GCL with increasing V , bypassing the GMG regime.In addition to demonstrating the direct transition from FN to GCL, Fig. 1d, which considers µ = 100 and T = 0.002 , shows that emission bypasses the RLD regime at low V , just as for lower µ in Fig. 1b.In all cases, regardless of µ , electron emission follows GCL at sufficiently high V since the electrons will eventually have enough energy to exhibit vacuum-like behavior in the presence of gas 16,19 .At low V , electron emission will be driven by RLD at sufficiently high T (i.e., v 0 ) or FN at sufficiently low T (i.e., v 0 ).Electron emission will only transition to GMG before GCL for a sufficiently low µ .The presence or absence of these various transitions may be assessed by using nexus phase space plots.
To characterize the transition between two electron emission mechanisms, we equate J from their asymptotic solutions.To simplify GMG and reduce numerical errors, we first derive asymptotic equations in the limits of high ( v 0 ≫ 1 but still nonrelativistic) and low ( v 0 ≪ 1 , but nonzero) velocity as and respectively.Setting J FN = J GMG,v 0 ≪1 yields www.nature.com/scientificreports/and J GMG,v 0 ≪1 = J RLD gives where T ∼ = 26.1145v 2 0 .Since ( 27), ( 28), ( 29), ( 30), ( 31), (32) arise from matching asymptotic solutions, they will not perfectly match the exact solution from (19).This is a standard characteristic of matched asymptotic analyses since the assumptions used to obtain each asymptote will inherently conflict with those to determine the other asymptote (or asymptotes for a higher order nexus) at their intersection 5 .However, these nexuses specify a regime where the dominant electron emission mechanisms would be sensitive to any D , V , v 0 , and µ.
Figure 2 illustrates the transitions between electron emission mechanisms by plotting V as a function of D , µ , and v 0 .The second-order nexus curves in Fig. 2 demonstrate the conditions where the dominant electron mecha- nisms transition.RLD consistently dominates when V is small, except for µ 0.3 for D = 250 and T = 0.02 , where MG dominates at low V because the gap is always space-charge limited and J MG ≈ 0 for sufficiently low V and v 0 = 0 .This suggests that combinations of µ and V below this threshold (or, alternatively, any pressure above this threshold) prohibit the emission of electrons into the gap due to strong collisionality and weak electric field (and concomitant force), respectively, resulting in a near-zero SCLC.Mathematically, for nonzero v 0 , sufficiently low µ and V increases the magnitude of D 4 v 4 0 µ −4 V −4 such that J GMG → 0 in (21), making it less than both J RLD and J FN so that the gap becomes space-charge limited at any V below this threshold µ .Increasing V eventually causes GCL to dominate at high V independent of µ (since J GMG increases with increasing V , making it nonzero even at low µ ), implying that the gap behaves like vacuum for any electron emitted at sufficiently high V .Figure 3 (32) At least one third-order nexus, corresponding to the intersection of three asymptotic solutions, appears in each panel of Fig. 2. Figures 2a and c show the third-order nexus where J GMG = J GCL = J FN , analogous to our prior study demonstrating the intersection of MG, CL, and FN for zero injection velocity 19 .The third-order nexus between FN, GMG, and GCL informs whether electron emission bypasses the GMG regime during the transition from FN to SCLC with decreasing D or increasing µ .Figure 2a shows that for D < D GMG,GCL,FN for a fixed v 0 and µ , electron emission transitions directly from FN to GCL without undergoing GMG, indicating that the gap behaves essentially as vacuum.Figure 2c shows that for µ > µ GMG,GCL,FN for a fixed v 0 and D , the electrons encounter sufficiently few collisions crossing the gap such that electron emission transitions directly from FN to GCL without undergoing GMG.
We can derive this condition analytically by rewriting the second-order nexus for J GCL = J FN as The third-order nexus J GMG = J FN = J GCL can be obtained bv setting J GMG = J FN and J GMG = J GCL , adding them to obtain 2J GMG = J FN + J GCL , and solving to yield Equations ( 33) and ( 34) describe the third-order nexus between FN, GMG, and GCL for some fixed D, V , or µ for a given v 0 .Comparing the solution from (33) to (34) to the solution for zero injection velocity 19 shows that incorporating v 0 induces a negligible change in the third-order nexus (< 0.1%) for 10 5 ≤ V ≤ 10 15 ; hence, we may use the zero injection nexus relationships, given by and ( 33) www.nature.com/scientificreports/for the given voltage range with minimal loss of accuracy, as done previously when incorporating thermionic emission 23 .Figure 3 from 19 shows the behavior of the third-order nexus under this condition.The impact of injection velocity on the third-order does not become significant until V < 10 3 , as shown in Fig. 4, which shows how V at the third-order nexus behaves as a function of D and µ for different v 0 , respectively.The divergence due to nonzero v 0 becomes particularly pronounced at low V for both low D and µ.Two third-order nexuses arise when considering V as a function of v 0 in Fig. 2b: one nexus between RLD, FN, and GMG and a second between RLD, GMG, and GCL.This indicates that for a given v 0 , the full solution transitions from RLD to FN to GCL, which is demonstrated in Fig. 1, or transitions from RLD to GMG to GCL, bypassing the FN regime.For the third-order nexus J GMG = J FN = J RLD , the condition can be derived analyti- cally by considering J FN = J RLD to obtain where W is the Lambert W-function, or the product log function.The third-order nexus can then be recovered by setting J GMG = J FN and J FN = J RLD , adding them to obtain 2J FN = J RLD + J GMG , and solving to yield Similarly, the condition for the third-order nexus J GCL = J GMG = J RLD can be derived analytically by first considering J GCL = J RLD to obtain where T ∝ v 2 0 .As above, considering J GMG = J RLD and J GCL = J GMG , adding them to obtain 2J GMG = J GCL + J RLD , and solving for µ gives Although not visible in Fig. 2, a fourth-order nexus may also occur between RLD, FN, GMG, and GCL for an appropriate combination of V , D, µ, and v 0 .Selecting any one of these parameters uniquely defines the other three to achieve the fourth-order nexus.To predict the fourth-order nexus, we first determine D with respect to v 0 and V .Since J RLD = J FN = J GMG = J GCL , we can solve for D independent of µ by considering J RLD = J GCL , yielding D by (39).The fourth-order nexus can then be recovered by setting J GMG = J GCL , J GMG = J FN and J GMG = J RLD , adding them to obtain 3J GMG = J GCL + J RLD + J FN , and solving for µ to obtain www.nature.com/scientificreports/Equations ( 39) and ( 41) describe the fourth-order nexus between RLD, FN, GMG, and GCL for a combination of D, V , µ , or v 0 , shown in Fig. 5.
Figure 6 compares the two space-charge-limited regimes, GCL and GMG, by plotting each with respect to the corresponding equation with zero injection velocity.From above, GMG differs from MG by a recurring term χ 2 = [Dv 0 / µV ] 2 , which represents the ratio of initial velocity v 0 to the nominal drift velocity v D = µV /D ; therefore, we plot J GMG /J MG as a function of χ 2 in Fig. 6a.This demonstrates that J GMG decreases and approaches zero with increasing χ 2 , which corresponds to v d ≪ v 0 .This means that insufficient electric field and/or low mobility ( µ → 0 , which corresponds to a strongly collisional gap) makes it difficult for the electrons to move (41)  in the medium.Since J GMG → 0 , as χ 2 → ∞ , the gap will always be space-charge limited since its strong col- lisionality or weak electric field prohibits electron emission until the applied voltage becomes sufficiently high, at which point the gap becomes space-charge limited by CL [cf.Fig. 2c].With decreasing χ 2 (i.e., v d ≪ v 0 ), J GMG → J MG .Figure 6b shows J GCL /J CL as a function of U = v 2 0 / 2V , which is the standard scaling for CL with nonzero v 0 for planar [49][50][51] and nonplanar 53,54 diodes.As U → 0 , J GCL → J CL ; as U → ∞ , J GCL continues to increase.At some point, relativistic effects become important 55 , which is beyond the scope of the current study.

Conclusion
This paper extends nexus theory to assess the transitions between RLD, FN, GMG, and GCL.Nexus theory provides a way to assess the dominant mechanisms before carrying out complicated simulations or experiments.Constructing the nexus phase space plot using the analytic equations for these mechanisms (including the GMG derived in this manuscript) with the desired operating conditions will demonstrate whether the equation for a single mechanism can be used (i.e., well away from one of the nexus curves) or if a more complicated equation combining multiple equations must be used (i.e., close to a nexus curve for two or more mechanisms).If this simple analysis shows that only one of these equations is dominant, this simplifies the design and simulation of the specific device; however, if operating near a nexus, then one must combine the mechanisms and use a more complete theory to predict behavior.The difference between using one of the simple equations can be significant, as illustrated by the 153.4% difference between GMG and the exact solution and 96.0%difference between FN and the exact solution in plot a of Fig. 1 at V = 25.73687near the nexus curve between FN and GMG.
The exact solution that accounts for temperature and mobility approaches the accepted equations for these mechanisms in the appropriate limits.Furthermore, we derive an analytic equation for GMG that includes the injection velocity, which is relevant for thermionic emitters under non-vacuum conditions.We also observe that applying a sufficiently small (but nonzero) bias voltage with a sufficient mobility yields a near-zero J GMG such that the gap is always space-charge limited, ultimately transitioning to GCL with increasing bias voltage.
We have recently experimentally and theoretically assessed the transitions between FN, MG, CL, and field emission-driven breakdown for nanoscale gaps at atmospheric pressure 7 and vacuum 56 near the third-order nexus between FN, MG, and CL.Other theories have examined thermo-field emission driven breakdown for microscale gaps at microwave frequencies 57,58 , indicating that the electrons undergoing these mechanisms need not only originate from cold cathodes.Thus, the results reported here, particularly for the fourth-order nexus between RLD, FN, GMG, and GCL, elucidate the contributions of the emission mechanisms (RLD and FN) and limits (GMG and GCL) for thermionic emitters in nano-and microscale devices that may undergo gas breakdown.Such theories may ultimately be extended to include additional emission mechanisms, such as photoemission 34 , or nonplanar geometries, which have been an ongoing area of study for both zero 59 and nonzero 53,54 injection velocity.Meadors and Poirier studied how to use a laser to heat a cathode to induce thermionic emission in vacuum and at atmospheric pressure without electromagnetic interference 60 .Another study addressed cathode heating and subsequent thermionic emission that play a critical role in arc formation 61 .Go 62 pointed out that ionenhanced thermo-field emission enabled the study of how slow-moving ions influenced thermionic emission in cathodic arcs [63][64][65][66][67] and then developed a theory for thermo-field emission driven microscale breakdown, extending the typi cal theories that considered strictly field emission driven breakdown 5 .As gap distances become smaller, the operating conditions approach the regime where the individual RLD, FN, MG, and CL may not capture the physics and a combined theory as derived in this paper becomes necessary.
Mobility may also be incorporated into recent work that unified RLD, FN, and SCLC in a vacuum crossedfield diode, where a magnetic field is applied perpendicular to the applied electric field 68 .All these applications of nexus theory demonstrate the importance of appropriately characterizing the dominant mechanism(s) to determine the current density.While perhaps not as important when considering two conditions (e.g., FN and CL), this becomes increasingly important as more mechanisms are added and the phase plot of contributing variables (e.g., voltage, gap distance, pressure, and temperature) increases.This complicates the transitions between the mechanisms and necessitates more care for ensuring that any theories or simulations properly account for the dominant mechanism.
Finally, we point out that the present paper considers the emission mechanism(s) predicted by the theory as coming from the full device.In other words, we do not consider that nonuniformities in temperature or electric field may result in different areas of the emitter undergoing emission mechanisms.Such a multidimensional model may be interesting to better understand behavior in these different regimes in future studies, but the approach here provides value for experimentalists to rapidly characterize the overall behavior of the measured current-voltage plots of the overall device. https://doi.org/10.1038/s41598-023-41615-2

Figure 1 .
Figure 1.Dimensionless current density J as a function of dimensionless voltage V for the exact solution represented by (19) and the asymptotic solutions GMG, GCL, FN, and RLD, respectively, with dimensionless gap distance D = 250 for (a) dimensionless mobility µ = 5 and dimensionless temperature T = 0.04 ; (b) µ = 5 and T = 0.002 ; (c) µ = 100 and T = 0.04 ; and (d) µ = 100 and T = 0.002 .For a low mobility (or high pressure), the full solution will follow GMG for sufficiently high V .For a sufficiently high mobility, electron emission bypasses GMG and transitions directly to GCL [cf.(c) and (d)].Low temperature eliminates the contribution of thermionic emission, removing the RLD regime from the full solution [cf.(b) and (d)].

Figure 2 .
Figure 2. (a) Dimensionless breakdown voltage V as a function of dimensionless gap distance D for dimensionless mobility µ = 5 and dimensionless temperature T = 0.02 .(b) Dimensionless breakdown voltage V as a function of dimensionless injection velocity v 0 , for dimensionless mobility µ = 5 and dimensionless gap distance D = 250 .(c) Dimensionless breakdown voltage V as a function of dimensionless mobility µ for dimensionless temperature T = 0.02 and dimensionless gap distance D = 250.

Figure 3 .
Figure 3. Dimensionless current density J as a function of dimensionless voltage V for the exact solution represented by (19) and the asymptotic solutions GMG and GCL, respectively, for dimensionless mobility µ = 0.1 , dimensionless gap distance D = 250 , and dimensionless temperature T = 0.02 .Asymptotic solutions FN and RLD are excluded because they are outside the relevant range.

Figure 4 .
Figure 4. Dimensionless voltage V as a function of (a) dimensionless gap distance D and (b) dimensionless mobility µ describing the third-order nexus between GCL, GMG, and FN.

Figure 5 .
Figure 5.The fourth-order nexus corresponding to J RLD = J FN = J GMG = J GCL uniquely defined by any one of the dimensionless voltage V , dimensionless injection velocity squared v 2 0 , dimensionless gap distance D , or dimensionless mobility µ.

Figure 6 .
Figure 6.(a) Ratio of dimensionless GMG current density J GMG to MG current density J MG as a function of χ 2 = [Dv 0 / µV ] 2 .(b) Ratio of dimensionless GCL current density J GCL to CL current density J CL as a function of the ratio of kinetic energy (KE) to potential energy (PE), given by U.

T , Table 1 .
Typical values of physical parameters.