Towards 4D modelisation of thermal-field emission from semiconductors

arXiv:2506.11927 · cond-mat.mtrl-sci, cond-mat.mes-hall · Submitted 2025-06-13 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: I'm Kai, and with me are Mira and Lev, guest researcher.

Mira: Today's paper: "Towards 4D modelisation of thermal-field emission from semiconductors".

Kai: The theoretical picture of thermal-field emission (TFE) from semiconductors has been limited to 1D and 2D models,

Mira: First, who's behind it and why it matters.

Title and authors: Kai: So, to summarize the core of "Towards 4D modelisation of thermal-field emission from semiconductors," the paper is about developing a comprehensive three-dimensional model for thermal-field emission that can take arbitrary geometries and doping levels into account <ref:2506.11927#pg0,Towards 4D modelisation of thermal-field emission from semiconductors>.

Mira: Essentially, it takes the previous limitations where TFE models were stuck in one or two dimensions and tries to bridge that gap by incorporating more complex physics.

Lev: What does this three dee modeling actually allow you to do that 1D models couldn't manage before <ref:2506.11927#pg0>? Does it solve some fundamental physical problem?

Kai: It allows them to solve for the interdependencies between the electric field, band structure, charge distribution, and temperature simultaneously through a set of coupled equations.

Mira: That self-consistent solution process is what’s important; it means the field affects the band bending, which changes how carriers are distributed based on doping, and that distribution then dictates where heat is generated.

Lev: So they are tackling the issue of how temperature influences the emission characteristics directly within their spatial model?

Kai: Precisely; they integrate thermal effects by calculating Nottingham heats for conduction and valence bands, Joule heating based on conductivity and field strength, and radiation loss at high temperatures.

Mira: That thermal integration is what elevates the work from a purely electronic structure study to a full thermo-electronic description of the process.

Lev: If we look at the final result, what's the most tangible physical insight they provide about TFE?

Kai: The main physical insight is that as the emitted current increases, the temperature rises because of internal heating and radiation loss, which then feeds back into modifying carrier concentrations via intrinsic density.

Mira: That feedback loop where current drives heat, which changes carriers, which changes emission characteristics is the key mechanism they are describing in this paper.

Lev: From a practical standpoint for error correction, this means we can't just assume a fixed emission rate; we have to account for thermal drift during operation.

Kai: So they successfully reproduce the non-linear I-V curves of semiconductors, showing the transition from quasilinear behavior in Region I to saturation in Region II.

Mira: That success is significant because it validates that their underlying physical assumptions about the material's response under bias are sound enough to model these characteristic features.

Lev: But we need to be careful; if this model doesn't handle transient events well, it won't be ready for real quantum computation environments.

Kai: They also found that saturation isn't just a simple drop in enhancement factor, but is linked to the Fermi level shifting, which provides a more detailed explanation.

Mira: That linkage between the macroscopic current behavior and microscopic electronic structure changes under high field is what gives this paper its depth.

Lev: We need to consider how complex that becomes when you try to map those three dee solutions onto a lattice structure for error correction simulations <ref:2506.11927#pg0>.

Kai: This paper provides a detailed look at how spatial charge density, field distribution, and thermal effects all work together within the semiconductor emitter.

Mira: It’s a lot of interconnected variables that need to be tracked together for any realistic simulation of device performance.

The paper's summary: Kai: Now we look at what improvements this paper suggests, which is essentially about moving from the current three dee model toward a full 4D framework <ref:2506.11927#pg0>.

Mira: They suggest extending this existing three-dimensional picture by incorporating temporal dynamics to include electron transport and heat dissipation in real time, which would take it into the fourth dimension.

Lev: Moving to 4D is ambitious; what kind of challenges do you foresee when you try to model time evolution alongside all that spatial complexity <ref:2506.11927#pg0>?

Kai: They state that their current three dee model serves as a solid platform for these future extensions, allowing for the incorporation of temporal dynamics <ref:2506.11927#pg0>.

Mira: The goal is to develop a system capable of simulating time-dependent phenomena, like how the temperature changes instantaneously during an emission pulse.

Lev: If you can do that, it would be incredibly useful for characterizing transient device behavior before you even deploy it on a quantum chip.

Kai: Future work will also focus on including other complex physics like photon enhanced field emission and quantum confinement effects as well.

Mira: Including those factors would really test the limits of the current model's ability to handle non-linear interactions, which is exactly what they need to do to fully validate it.

Lev: And addressing surface states explicitly in a quantitative way, rather than just qualitatively noting their presence, seems like a necessary step for any high-fidelity simulation.

Kai: They aim to quantify the relative contribution of different emission mechanisms—surface states, conduction band emission, and valence band emission—under varying field conditions.

Mira: Quantifying those contributions is crucial because it moves the analysis from just describing phenomena to predicting which mechanism dominates at any given operating point.

Lev: That level of quantitative detail would make it much more useful for designing devices where you need to understand precisely how surface effects influence the overall performance metrics.

Kai: Ultimately, this paper sets up a three dee picture of charge distribution, field, band structure, and temperature that is ready for those next steps into a full four-dimensional framework <ref:2506.11927#pg0>.

Mira: It’s a strong foundation because it provides the necessary spatial context to eventually tackle the time dynamics that are currently missing.

Lev: It gives us a clear roadmap on how to build up complexity systematically, which is helpful when you're trying to plan out long-term research projects.

Kai: So, essentially, they’re laying down a sophisticated three dee spatial and thermal picture before they start looking at the full time evolution of the system <ref:2506.11927#pg0>.

Mira: It’s a very logical progression to build up complexity this way, ensuring each layer is physically justified before adding the next dimension.

Lev: I just hope that when they get to those temporal components, you can actually manage the computational load without needing a supercomputer for every simulation run.

Kai: That’s the challenge ahead—making sure the future 4D model remains computationally tractable while incorporating all that physics <ref:2506.11927#pg0>.

The paper's improvements: Kai: So, to conclude this discussion on "Towards 4D modelisation of thermal-field emission from semiconductors," we see a paper that successfully develops a three-dimensional model capable of handling arbitrary geometries and doping levels for TFE <ref:2506.11927#pg0,Towards 4D modelisation of thermal-field emission from semiconductors>.

Mira: The key success here is the self-consistent solution process that links field determination to carrier concentration evaluation, which allows them to capture the essential non-linear behavior accurately.

Lev: It’s clear that this work provides a solid theoretical groundwork for understanding how temperature influences emission characteristics through internal heating and radiation loss in these devices.

Kai: The paper also highlights the finding that saturation is tied to the Fermi level shift, offering a deeper explanation than just an empirical observation.

Mira: Overall, it’s a very thorough description of the coupled physics involved in thermal-field emission within semiconductor emitters.

Lev: For us in error correction research, this model provides a tangible system to analyze how these complex factors interact before we try to apply them to actual hardware.

Kai: So, this paper is a valuable piece of work because it gives us a much better spatial understanding of the charge, field, band structure, and temperature dynamics.

Mira: We're excited about how this sets up the next phase for modeling these emitters with 4D temporal components <ref:2506.11927#pg0>.

Lev: I think this is where the real progress will happen in bridging theory and practical application for complex systems like this.

Kai: We’ve covered the development of the three dee model, its key findings on non-linear I-V curves, and how thermal effects are integrated into their analysis <ref:2506.11927#pg0>.

Mira: This paper really establishes a robust framework for modeling these emitters that is much more sophisticated than what was available previously.

Lev: We can now see exactly where the gaps are for future research, which helps us target our next research efforts effectively.

Kai: That's what we discussed regarding the three dee picture of charge distribution, field, band structure, and temperature dynamics in this paper <ref:2506.11927#pg0>.

Conclusion: Kai: So, we’ve walked through the development of this three dee model for thermal-field emission from semiconductors by looking at all those governing equations and how they self-consistently solve the field, band structure, charge distribution, and temperature together.

Mira: Exactly; it's a pretty solid framework because it forces you to account for those thermal effects like Joule heating and radiation loss directly into the electronic structure calculations.

Lev: From my side, I’m really interested in how this three dee spatial picture of charge and field would translate when we try to run this on real hardware, especially considering the need for time-dependent simulations.

Kai: Right, and that's where the future work is heading—they are looking to extend this into a full four-dimensional framework to include temporal dynamics for electron transport.

Mira: That extension is vital because it lets us see how temperature fluctuations impact the emission rate in real-time during device operation, which is something we can't do with a static three dee model.

Lev: I agree; if we can simulate transient events like that, it becomes much more useful for characterizing device behavior before deployment.

Kai: This paper, "Towards 4D modelisation of thermal-field emission from semiconductors," gives us a really detailed spatial picture of how all these variables interact in a semiconductor emitter <ref:2506.11927#pg0,Towards 4D modelisation of thermal-field emission from semiconductors>.

Mira: It’s a lot of interconnected physics they managed to tie together; the way they model saturation through the Fermi level shift is quite insightful.

Lev: I think that connection between macroscopic current behavior and microscopic electronic structure changes under high field is what really gives this paper its depth for error correction simulations.

Kai: That’s exactly what we need to see when we think about mapping these three dee solutions onto a lattice structure for those future quantum simulations.

Mira: So, while it’s not a complete 4D simulation yet, it provides the necessary spatial context to build that next layer of complexity on top of <ref:2506.11927#pg0>.

Lev: I’m just hopeful that when they get to those temporal components, they can manage the computational load without requiring massive supercomputers for every single run.

Kai: We'll see if they can do that, but for now, this paper is a huge step forward in establishing the physical foundation for these kinds of simulations.

Mira: It’s a great piece of work because it lays down a really robust platform before we start looking at the more complex physics like photon enhancement later on.

Lev: That's where I hope we see some interesting results, though I always keep an eye out for those surface state contributions to be fully quantified.

Kai: We’ll keep an eye out for those next steps in the research as they try to incorporate things like quantum confinement into this model.

School of Engineering, University of Edinburgh · Institute of Technology, University of Tartu

cond-mat.mtrl-sci, cond-mat.mes-hall

Submitted: 2025-06-13

Updated: 2025-06-13

DOI: 10.1063/5.0302109

License: http://creativecommons.org/licenses/by-nc-sa/4.0/

Importance score: 80/100

The gist: The theoretical picture of thermal-field emission (TFE) from semiconductors has been limited to 1D and 2D models, which this work addresses by developing a 3D model capable of incorporating arbitrary

Key concepts

Poisson’s Equation: This equation describes how the electric field bends within a semiconductor material. The degree of this bending is directly controlled by both the external applied electric field and how heavily the material is doped with impurities, which dictates where charge accumulates.
Nottingham Heats: These are specific heat calculations used to account for temperature effects in the model. They are derived from integrals involving current densities to determine how heat is generated within the conduction band (QNC) and valence band (QNV), which is essential for modeling thermal changes.
Field Enhancement Factor ($eta$): This factor quantifies how much stronger the electric field becomes at the very tip of a sharp emitter compared to the average field. It is calculated by dividing the maximum field ($E_x$) by a reference field ($E_{cro}$), helping define emission characteristics.
Saturation Mechanism: The model found that current saturation in semiconductors is not caused by the electric field enhancement decreasing. Instead, it occurs because the Fermi level drops, which aligns with previously understood physical mechanisms for current limiting in these devices.

Terminology

APEX — 2× in this paper · explained in 1 episode
Core — 1× in this paper · explained in 2 episodes
Degeneracy — 1× in this paper · explained in 3 episodes
Fidelity — 1× in this paper · explained in 2 episodes

Summary

The theoretical picture of thermal-field emission (TFE) from semiconductors has been limited to 1D and 2D models, which this work addresses by developing a 3D model capable of incorporating arbitrary geometries and doping levels to provide a more realistic description of thermal-field emission.

Model Development and Governing Equations

The model is built upon solving the interdependencies between the electric field, band structure, charge distribution, and temperature. The foundation involves Poisson’s Equation to describe band bending:

“The degree of band bending depends on the applied field and the doping level of the semiconductor, and can be described using the Poisson’s Equation [8]:”

The total space-charge density is given by:

ρ = q(ND+−NA− + n − p) (2)

The model requires self-consistent solution through iterative loops. The process involves:

Solving Eq. 1 (blue loop) for a given applied potential to determine the electric field and band structure.
Evaluating Eq. 2 using formulae for respective densities: electron concentration, hole concentration, and ionized donor/acceptor concentrations (Equations 3 through 6).
Applying continuity conditions to calculate the electron current density inside the emitter (Equations 14 through 20).

Thermal Effects Integration

A crucial aspect of the model is incorporating thermal effects to account for temperature dependence, which is vital for a complete description. The heat components calculated are:

“Now, the Nottingham heat can be calculated as [9]:”

The Nottingham heats are defined by integrals involving current densities:

The conduction band heat: QNC(E) = ∫ Jc(E) / E0 (E − ER) dE (21).
The valence band heat: QNV(E) = ∫ JV(E) / E0 (ER − E) dE (22).

The Joule heating is calculated as:

QJ = σF2 (24)

Radiation loss is included for high temperatures:

QR = −εσ(T04 − T4) (25)

Field Emission Characteristics Analysis

The model successfully reproduces the characteristic non-linear I-V curves of semiconductors, showing a transition from a quasilinear behavior in Region I to saturation in Region II. Key findings regarding emission characteristics include:

The emission area is defined using the concept of notional area: A r = I / jr (28).
The field enhancement factor β is calculated as: β = E x / E cro (29), where Emax is the field at the apex and Emacro results from applying a potential V between anode and cathode (Equation 30).
Saturation is not attributed to a drop in the enhancement factor, but rather to a change in the position of the Fermi level:

**The saturation regime coincides with a drop of the Fermi level, which is consistent with previously reported mechanisms for current saturation [6]. **

The temperature dependence is significant; as temperature rises, intrinsic carrier density increases exponentially, leading to an exponential increase in emitted current until saturation occurs.

Experimental Validation and Interpretation

The model was validated by simulating a germanium emitter based on experimental data from Shepherd and Peria (S&P) [13]. The comparison showed:

**“Figure 10 shows a comparison between the experimental and simulated I-F (current-field) and it can be seen that they are in good agreement.” **

However, discrepancies in the high current regime were attributed to surface states. The paper proposes an alternative interpretation of S&P's results:

The large peak in the electron energy spectrum is attributed to electrons coming from the conduction band.
The smaller peak is attributed to the surface states of clean germanium, which contribute more as the potential barrier thins at higher fields, while bulk emission saturates due to a Fermi level drop.

Future Directions

The current model provides a 3D picture of the distribution of charges, field, band structure, and temperature, allowing for future extensions into the full 4D framework by incorporating temporal dynamics for electron transport and heat dissipation. Future work aims to include the physics of photon enhanced field emission, quantum confinement, and surface states.

The gist: A 3D model is developed that successfully reproduces characteristic non-linear I-V curves and temperature dependence for thermal-field emission from semiconductors by self-consistently solving coupled equations for electric field, band structure, charge distribution, and temperature.

How it works

The model first solves Eq. 1 (blue loop) for a given applied potential to determine the electric field and band structure. The process then evaluates Eq.

Improvements for AI systems

Here are the specific improvements that can be made to AI systems based on this scientific paper:

Improved Simulation Fidelity for Semiconductor Device Characterization: The model provides a 3D, self-consistent framework incorporating Poisson's equation, band structure (including temperature dependence), carrier concentrations, and thermal effects (Joule and Nottingham heating).
Enhanced Predictive Capabilities for Field Emission Saturation: The system can accurately predict the characteristic non-linear I-V curves of semiconductors across arbitrary geometries and doping levels. Specifically, it can model the physical origin of current saturation by tracking the dynamic change in the Fermi level position relative to band edges under high fields, rather than relying solely on empirical fitting.
Integrated Thermal Management Analysis: The AI can calculate internal emitter temperature distributions by solving a coupled heat transfer equation (including radiation loss). This allows for predictive modeling of how temperature fluctuations impact emission characteristics and material properties (like thermal conductivity) in real-time during device operation, which is crucial for reliability assessment.
Advanced Parameter Extraction and Sensitivity Mapping: The model enables the extraction and dynamic study of key field emission parameters—emission area, field enhancement factor, and work function—as functions of applied potential, doping type, geometry (shape/cone angle), and temperature. This allows AI systems to map the sensitivity of device performance to various operating conditions.
Surface State Contribution Modeling: By incorporating the analysis derived from comparing model results with experimental electron energy distributions (specifically addressing surface states vs. bulk conduction band emission), the AI can move beyond simplified models by quantifying the relative contribution of different emission mechanisms (surface states, conduction band, valence band) to the total current under varying field conditions.
Foundation for Full 4D Temporal Modeling: The established 3D spatial model serves as a robust platform for future extensions into a full 4D framework. AI can be used to develop temporal components (using the provided structure) to simulate time-dependent phenomena like electron transport dynamics and transient heat dissipation, leading to truly predictive, time-resolved device simulations.

These improvements allow an AI system to transition from simple curve fitting or linear approximations to a high-fidelity, physics-informed simulation engine capable of designing and predicting the behavior of advanced semiconductor field emitters.

Abstract

The theoretical picture of thermal field-emission (TFE) from semiconductors has been limited to 1D and 2D models. This can be attributed to the complex and interdependent phenomena that is involved in TFE from semiconductors which makes the calculations cumbersome. Such limitations result in a partial understanding of the underlying physics of semiconducting surfaces under high electrical fields, which requires the addition of the temporal dimension (4D) to yield a realistic model. Here we develop a 3D model of TFE from semiconductors that can take arbitrary geometries and doping levels. Our model successfully reproduces the characteristic saturation plateau of some semiconductors, as well as its dependence in temperature. The model is found to be in good agreement with experimental data from ntype Germanium at a qualitative level. We propose this model as a platform for future extensions into the full 4D framework, incorporating temporal dynamics for a more complete and predictive description of thermal-field emission from semiconductors.