# A Nanowire with charge neutrality in band gap

Tags:
1. Mar 24, 2017

### Verdict

I'm working on a finite element simulation of the electrostatic potential $V(r)$ in and around semiconductor nanowires, based on solving Poisson's equation. While the details of the problem will follow shortly, the crux of where I run into trouble is that for nanowires it is important to include the effects of surface states. This is usually done in the form of a constant surface state density, but I want to include that the surface represents a capacitor with a variable charge as a function of the Fermi level.

To do so, I will include the so called charge neutrality level (CNL) $\Phi_{NL}$, following the paper Resolving ambiguities in nanowire field-effect transistor characterization by Heedt et al. (which should be open access). According to this, one can take the surface states and the charge neutrality level into account by assuming a finite density of surface states $D_s$ and a surface charge density given by

\sigma_s(r) = D_s\left[\Delta\Phi_{NL} - eV(r) - \chi_{InAs}\right]

where $\Delta\Phi_{NL} = \Phi_{NL} - E_c$ is the location of the charge neutrality level w.r.t. the conduction band edge $E_c$ and $\chi_{SE}$ is the electron affinity of the semiconductor. Note that this expression is what sets the model apart from many others.

Now, the above gives us the surface density. What remains is to define the the space charge density. This can be done with the Kane model in the context of $k \cdot p$ perturbation theory, where one writes the electron concentration in terms of its density of states integrated with the fermi-dirac distribution

n(r) = \int_{E_c}^{\infty} \frac{(2m_e^*)^{3/2}}{2 \pi^2 \hbar^3}\sqrt{E-E_c} \frac{1}{\exp(\frac{E-V(r)}{k_B T})}dE

and similarly for the hole concentration

p(r) = \int_{-\infty}^{E_c-E_g} \frac{(2m_p^*)^{3/2}}{2 \pi^2 \hbar^3}\sqrt{E-(E_c-E_g)} \left(1-\frac{1}{\exp(\frac{E-V(r)}{k_B T})}\right)dE

leading to a space charge density of

\rho(r) = e\left[p(r) - n(r)\right]

In the above $m_{\{e,p\}}$ are the effective electron and hole masses, and $E_g$ is the band gap.

Okay, great! At this stage we actually have all that is needed to set up a finite element simulation by drawing the hexagonal wire, defining the space charge and the surface charge density, surrounding the wire with something like air and filling in appropriate parameters.

But here we are getting to my problem. What the above model is supposed to do (if I understand it correctly) is find $V(r)$ so that it can match $\rho(r)$ and $\sigma(r)$, essentially bending the bands. This works for $\chi_{NL} > E_c$, where the CNL lies in the conduction band and results in an electron accumulation layer, and also for $\chi_{NL} < E_c - E_g$, where it lies in the valence band. But how should the above model work for $E_c - E_g <\Delta\chi_{NL} < E_c$? In this region the space charge density expression is simply zero, as the integrals are only evaluated outside of the bandgap. So the model can't exactly modify $V(r)$ to bend the bands into the region.

If I understand it correctly the model thus needs a modification that allows the bands to go into the gap around the surface, and I am not sure how to do this. What I've come to realise is that what one essentially needs is that $E_c \rightarrow E_c(V(r))$; this would make the limits of the integrals move with the potential, as they should. But how does the conduction band edge respond to the potential? That seems rather tricky..

[1]: http://pubs.rsc.org/en/Content/ArticleLanding/2015/NR/C5NR03608A#!divAbstract

2. Mar 29, 2017

### PF_Help_Bot

Thanks for the thread! This is an automated courtesy bump. Sorry you aren't generating responses at the moment. Do you have any further information, come to any new conclusions or is it possible to reword the post? The more details the better.