Field-Effect-Aware Fluctuation Tunnelling (FEAFIT)

• Field-Effect-Aware Fluctuation-Induced Tunnelling is a framework describing electron transport in disordered semiconductors and granular materials through quantum tunnelling across insulating gaps.
• It integrates field-effect modulation with stochastic thermal and electromagnetic fluctuations to extract material-specific transport properties directly from device measurements.
• The model explains the transition from low-temperature tunnelling plateaus to high-temperature Arrhenius conduction, aiding optimization in systems like IGZO and granular metals.

Field-Effect-Aware Fluctuation-Induced Tunnelling (FEAFIT) is a unified theoretical framework for understanding electronic transport in highly disordered semiconductors and granular materials. It rigorously accounts for tunnelling across nanoscale gaps, incorporating both stochastic thermal and electromagnetic fluctuations and field-induced modulation of the tunnelling barriers. FEAFIT enables extraction of material-specific transport properties directly from device measurements, providing quantitative links between experimental observations, first-principles calculations, and microscopic disorder. Its predictive capacity has been validated in amorphous oxide semiconductors such as IGZO and in granular metal systems subjected to high-frequency fields (Zhao et al., 26 Dec 2025, Hirlimann, 2013).

1. Physical Basis: Partially Coherent Domains and Tunnelling Junctions

FEAFIT is fundamentally built on three physical assumptions, with their realization most clearly elucidated in amorphous IGZO thin films (Zhao et al., 26 Dec 2025):

1. Partially coherent electronic states: Electrons retain spatial coherence only within discrete nanometer-scale domains or "puddles." Within a domain, electronic states are well described by plane waves; coherence is lost at domain boundaries, precluding globally delocalized transport.
2. Non-degenerate conduction: The Fermi level $E_F$ is positioned within a dense manifold of localized tail states below the threshold for fully delocalized states ($E_{dloc}$). Global free-electron conduction does not occur; instead, transport is dominated by transfer across energetic barriers.
3. Insulating gaps as tunnelling junctions: Electrical conduction proceeds by quantum tunnelling across narrow, energetically defined insulating gaps that separate coherent domains. These junctions have a width $d$ and barrier height $\Phi$, both modulated by external electric fields such as gate voltages in transistors.

This paradigm is supported by both Hall effect measurements, which reveal the coexistence of coherent and incoherent carriers, and first-principles simulations showing separated spatial pockets of conduction-band states whose separation matches the extracted tunnelling gap width (Zhao et al., 26 Dec 2025).

2. Tunnelling Action, Fluctuations, and Field Modulation

WKB Transmission and Fluctuation-Induced Enhancement

The tunnelling probability for an individual junction is governed by the WKB transmission action. At zero temperature, the action is given by

$$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$

where $m^*$ is the effective mass and $E$ the incident electron energy.

Thermal and electromagnetic field fluctuations transiently reduce the barrier height, thus boosting tunnelling rates. The junction can be mapped to a nanoscale parallel-plate capacitor (capacitance $C = \epsilon A/d$), so the variance of the voltage drop across the barrier, $\langle(\delta V)^2\rangle = k_B T / C$, characterizes the amplitude of thermal fluctuations. Applying Sheng’s theory, averaging over these fluctuations gives a smooth interpolation between elastic (low-$T$) tunnelling and high-$E_{dloc}$0 Arrhenius activation (Zhao et al., 26 Dec 2025).

Field-Effect Modulation

For field-effect devices, barrier parameters are modulated by the external field (e.g., gate voltage $E_{dloc}$1):

• Effective barrier height: $E_{dloc}$2
• Barrier width: $E_{dloc}$3

Here, $E_{dloc}$4 and $E_{dloc}$5 parameterize the sensitivity of the respective quantities to field, and $E_{dloc}$6 is the flat-band voltage.

3. Analytical Model: FEAFIT Conductivity and Resistance Expressions

The FEAFIT formalism provides closed-form analytical expressions for observables. For channel conductivity in a thin-film transistor:

$E_{dloc}$7

where:

• $E_{dloc}$8 is the domain-scale conductance prefactor,
• $E_{dloc}$9 is the gate-dependent effective activation energy,
• $d$0 encodes the characteristic energy scale set by barrier width,
• $d$1, $d$2 are crossover functions defined via a reduced fluctuation variable $d$3, with $d$4 and $d$5.

This form naturally recovers both the high-$d$6 Arrhenius-like behavior ($d$7 for $d$8) and the low-$d$9 tunnelling plateau ($\Phi$0 for $\Phi$1) without artificial regime splicing (Zhao et al., 26 Dec 2025).

For granular metallic contacts under electromagnetic stimulation, the resistance is modeled as

$\Phi$2

with $\Phi$3 the field amplitude, $\Phi$4 the frequency, $\Phi$5 the nominal barrier, and $\Phi$6 a frequency-domain cutoff controlled by the electron reflection time $\Phi$7 (Hirlimann, 2013).

4. Extraction and Correlation of FEAFIT Parameters

Comprehensive device data enable direct extraction of the key FEAFIT parameters:

• Barrier height and width ($\Phi$8, $\Phi$9): Obtained by global nonlinear least-squares fitting of conductivity data as a function of gate voltage and temperature. The dependence of $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$0 and $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$1 on $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$2 quantifies $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$3, $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$4, and their field couplings $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$5, $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$6.
• Coherence length $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$7: Inferred from the extracted effective area $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$8, using $$S_0 = 2\int_0^d \sqrt{2m^*[\,\Phi - E\,]} \, dx \approx 2d \sqrt{2m^*\Phi}/\hbar$$9.
• Energetic disorder: Tail-state density-of-states analyses yield disorder parameters (widths $m^*$0, $m^*$1), which are found to tightly correlate with the empirical barrier energies ($m^*$2).

These parameters track the material’s composition. For example, first-principles supercell calculations for IGZO with varying Ga content demonstrate that as Ga fraction increases, the mean separation between coherent domains ($m^*$34 to 8 nm) grows in parallel with fitted $m^*$4, and the energetic disorder extracted from density of states matches the FEAFIT-determined barriers (Zhao et al., 26 Dec 2025).

5. Universality and Applicability Across Materials Systems

The FEAFIT model generalizes beyond amorphous oxide semiconductors to describe granular metals and classical Branly-effect systems. For each metal–insulator–metal (MIM) junction, the application of field-dependent induced tunnelling and stochastic fluctuation terms directly recovers measured scaling laws:

• Low field/temperature: pure tunnelling regime with resistance $m^*$5.
• Intermediate field: field-induced enhancement, $m^*$6, matching observed quadratic-in-field resistance decreases up to the threshold for barrier suppression.
• High field/low frequency: saturation or irreversible drops due to contact welding (Hirlimann, 2013).

The framework thus covers continuous variation from pure quantum tunnelling, fluctuation-assisted activation, to strong field-driven nonlinearity. In field-effect transistors, it explains monotonic enhancement of conductivity with gate voltage, the crossover from Arrhenius to plateau forms, and composition-dependent suppression or enhancement of low-temperature conduction plateaus (Zhao et al., 26 Dec 2025).

6. Limitations, Boundary Conditions, and Materials Engineering Implications

FEAFIT, while predictive, is subject to several key boundary conditions:

• Assumes 1D square barriers; real systems present distributions in area and thickness, requiring statistical averaging or network representation for full accuracy (Hirlimann, 2013).
• Neglects heating and feedback below the threshold for contact welding; sustained field application may introduce additional nonlinearities.
• The induced-tunnelling enhancement is contingent on sub–$m^*$7 (fs-scale) coherent field spikes; long-pulse fields produce only WKB modulation, setting a frequency cutoff at $m^*$8.
• In low-temperature regimes, quantum (zero-$m^*$9) fluctuations rather than thermal noise set the variance, requiring appropriate modification of $E$0.
• The coherence domain framework is inapplicable if global (delocalized) band conduction is restored, i.e., $E$1 rises above $E$2.

From a materials engineering perspective, FEAFIT quantitatively demonstrates that raising carrier mobility and on-current in AOS or granular systems hinges on maximizing the size of coherent domains (e.g., by increasing In fraction or optimizing O content) and minimizing tail-state energetic disorder through compositional control (e.g., reducing Ga-induced fluctuation). All extracted parameters are directly interpretable in terms of microscopic disorder and coherence, providing an actionable roadmap for device optimization (Zhao et al., 26 Dec 2025).

7. Tables of Key FEAFIT Quantities and Relationships

The following table summarizes the key variable dependencies for FEAFIT transport in a-IGZO and granular systems:

Quantity	Physical Meaning	Functional Dependence
$E$3	Effective barrier height	$E$4
$E$5	Effective barrier width	$E$6
$E$7	Activation scale	$E$8
$E$9	Tunnelling scale	$C = \epsilon A/d$0
$C = \epsilon A/d$1	Crossover functions	$C = \epsilon A/d$2, $C = \epsilon A/d$3 with $C = \epsilon A/d$4
$C = \epsilon A/d$5	Contact resistance (Branly type)	See above; fields, frequency, and temperature

Empirical correspondence with first-principles calculations and observed device behavior underpins the universality and robustness of the FEAFIT approach (Zhao et al., 26 Dec 2025, Hirlimann, 2013).