Entropy is a thermodynamic state function. For a reversible heat transfer, ; the Clausius inequality generalizes this to irreversible processes. The combined entropy of a system and its surroundings does not decrease. A subsystem can lose entropy by exporting heat or matter, so the second-law statement must specify the whole isolated combination. In a simple compressible system,
which identifies temperature, pressure and chemical potential as derivatives of the energy representation of the state function. Thermodynamic equilibrium maximizes entropy under fixed total energy, volume and particle number. A heat reservoir instead leads to minimizing the Helmholtz free energy at fixed , or the Gibbs free energy at fixed .
Statistical mechanics explains the state count. In a microcanonical ensemble, equally probable compatible microstates give Boltzmann's entropy formula, . With general discrete probabilities, the Gibbs entropy is
The microcanonical formula follows by setting every nonzero probability to . For a canonical ensemble, , so
This turns microscopic energy levels and their populations into a measurable thermodynamic quantity. Classical phase-space counting needs a consistent microscopic normalization and indistinguishability convention; quantum state counting supplies that normalization. Independent systems have multiplicative state counts and therefore additive entropies.
For an ideal gas freely expanding from volume to in an insulated enclosure, energy remains fixed while each particle gains twice the accessible position volume. The state count increases by , giving . The actual irreversible process transfers no heat. A reversible isothermal comparison path transfers heat and computes the same entropy change because the initial and final thermodynamic equilibrium states are the same.
Calorimetry supplies another connection to experiment: at fixed pressure and within one phase, , with latent-heat contributions at first-order transitions. The third law of thermodynamics fixes for an ideal perfect crystal with a unique ground state. Degeneracy or frozen disorder can produce residual entropy. Entropy thus links heat-capacity measurements, phase thermodynamic equilibrium and the number of accessible microscopic configurations.
For a bulk liquid and vapor at common fixed temperature and pressure, transfer molecules from vapor to liquid while holding their sum fixed. The Gibbs free energy change is . Thermodynamic equilibrium with both phases present requires stationarity under either sign of transfer, hence . This is the bulk coexistence condition; an interface's curvature adds a surface-energy cost.
In the capillarity approximation, take an incompressible drop with constant liquid number density , bulk chemical potential at the specified reservoir conditions, and size-independent surface tension . The liquid bulk is extensive, so its Gibbs thermodynamic free energy is . Creating interface area reversibly at fixed temperature adds . Equivalently, along growth at fixed bulk parameters the stated differential integrates as , with the no-drop reference set to zero. For a sphere,
This approximation uses bulk properties for the interior and neglects the molecular interface thickness, curvature dependence of surface tension and higher compressibility corrections. It does not apply quantitatively to drops only a few molecules across.
The vapor contains molecules. Treat its chemical potential as approximately unchanged by forming the drop, as for a large surrounding reservoir or negligible vapor depletion. Then
The radius is limited by the available molecules in a finite sample; the constant-reservoir expression is a local nucleation model, not a description of arbitrarily complete vapor depletion. The excess over the vapor reference has derivative
If , it increases for every , leaving only the stable boundary minimum . If , is a metastable local minimum and the additional stationary radius is
At that radius , so it is an unstable maximum. The classical nucleation theory barrier is .
Figure 1.
Surface cost and bulk gain produce a critical droplet in supersaturated vapor
. Original dimensionless capillarity curves. The supersaturated case has a barrier at the critical radius; the other curve rises monotonically. These are reservoir-model free energies, before finite-sample depletion changes the chemical potentials. Details.
A fluctuation below tends to shrink because reducing its surface cost lowers thermodynamic free energy. A supercritical drop tends to grow, because its bulk free-energy gain outweighs the added surface cost. Growth continues while material supply and supersaturation allow it, eventually changing the vapor state or leading to macroscopic liquid. Curved-interface pressure and chemical thermodynamic equilibrium can also be described through the Young–Laplace equation and the Kelvin equation.
Semiconductor doping introduces electrically active impurities. In gallium arsenide, silicon has one more valence electron than gallium but one fewer than arsenic. Silicon on a gallium site can therefore supply a shallow donor electron, while silicon on an arsenic site can accept an electron and leave a mobile hole. It is an amphoteric dopant, so its net effect depends on site occupation and growth conditions.[ref]
For isotropic parabolic bands, shift the energy origin in the given free-electron density of states to the appropriate band edge and use its effective mass:
and, for valence electron states,
The stated normalization includes the two spin states and assumes one effective valley. In a nondegenerate semiconductor, the conduction occupation is approximately and the probability of a valence hole is approximately . Substituting in the first integral, and the corresponding hole energy in the second, gives
where the effective density of states values are and . The factor comes from the supplied Gaussian-weighted integral. Multiplying cancels the shared thermodynamic equilibrium chemical potential and proves the semiconductor mass-action law:
For the given sample and . With complete donor ionization and occupied acceptors, charge neutrality gives . At the given masses and gap give
Solving the compensated semiconductor equations, or using , yields
The very small minority density is consistent with the wide gap and the strongly n-type balance. Also and , checking the nondegenerate approximation. The acceptor levels are occupied in this n-type regime; complete acceptor ionization in the compensation equation is appropriate to the model.
At low temperature dopant freeze-out reduces the conduction-electron population; a quantitative curve would require dopant binding energies and occupation factors. Once donors are ionized while intrinsic excitation is negligible, gives an extrinsic plateau. At higher temperature intrinsic pairs eventually dominate, so if the gap and masses are treated as constant. Real gap and mass changes modify that last expression.
In a conductor, mobile electrons redistribute in response to an electric potential, producing an induced charge density that opposes its source. Finite electronic compressibility means a finite density change requires a change of local kinetic Fermi energy. Screening is therefore not perfect at every microscopic distance. On scales comparable to atomic structure or the Fermi wavelength, a continuum local-density model itself fails.
Thomas-Fermi screening treats electrons as a locally uniform degenerate free electron gas, replaces the ionic background by a smooth neutralizing charge, and neglects nonlocal response, exchange and correlation corrections. The electrochemical chemical potential is uniform in equilibrium. Because an electron has charge , its electrostatic energy is , giving
At zero potential . Linearizing at fixed , with , gives
Insert this in Poisson's equation to obtain , where the Thomas-Fermi screening length is
Away from a point source, direct radial differentiation shows . The decaying solution is therefore . Integrating Poisson's equation over a small sphere containing charge fixes . For a sufficiently large source, the linearized formula is only a far-field approximation, not a justified description of its singular core.
For monovalent sodium, the electron number density is its atomic number density,
This gives and
Room temperature is far below the corresponding Fermi temperature, so degeneracy is well justified, and sodium's simple conduction band makes a free-electron estimate useful. However and the mean atomic spacing is . The potential varies on a scale shorter than and the atomic spacing, so the slowly varying local-density and continuum assumptions are not quantitatively controlled at the screening length itself. A nonlocal treatment includes Friedel oscillations, absent from the pure exponential.
Nearly unchanged interlayer force constants in two-layer and twenty-layer films are consistent with very local screening and bonding forces: reducing the distant bulk need not strongly change the force between adjacent layers. They also indicate that metallic electronic response can survive in the thin film. The observation does not by itself validate the numerical Thomas-Fermi approximation, establish an identical full phonon spectrum, or rule out surface-induced electronic changes; it directly constrains the measured local force constant.
A positive electric potential difference between gate and substrate attracts electrons towards the insulating interface and repels the mobile electron holes of the p-type semiconductor substrate. The exposed, negatively charged acceptor impurities form a depletion region; the thin mobile electron layer at its surface is the semiconductor inversion layer. The conduction band and valence band bend downwards towards the positively biased gate, because the electron potential energy is .
Figure 1.
Band bending and charge separation at the conventional MOS inversion threshold
. Details.
To specify the conventional MOS inversion threshold, let be the ionized acceptor density, the intrinsic carrier concentration, and
In the neutral bulk the intrinsic level lies above the semiconductor chemical potential. A positive surface potential lowers both band edges by , giving, in the nondegenerate approximation,
At , the surface is intrinsic: . The usual onset of strong inversion instead means
Thus the surface electron density has grown to the bulk majority-carrier density. The chemical potential need not actually enter the conduction band. The phrase “just forms” can also refer to the earlier weak-inversion crossover; the diagram uses the conventional strong-inversion threshold and labels it explicitly.
Within the depletion approximation the negative space charge is approximately over a width
where is the absolute semiconductor permittivity. Its outer edge is where the bands become flat and mobile electron holes recover their bulk density. The electrical insulator blocks steady carrier transport while sustaining an electric field; the metal gate supplies positive surface charge balancing depletion and inversion charge. Under gate bias, the gate and substrate electrochemical potentials are different, so a single equilibrium Fermi level should not be drawn through the whole biased stack.
A solar cell uses the photovoltaic effect to turn photon energy into electrical work. In a semiconductor, above-gap photon absorption generates an electron in the conduction band and an electron hole in the valence band. A p-n junction, together with selective electrical contacts, collects the two carrier species at opposite terminals. Carrier diffusion and the built-in electric field help separate them before electron-hole recombination. Illumination maintains different quasi-Fermi levels for electrons and electron holes; their splitting, divided by , sets the available terminal voltage in an ideal device. The equilibrium built-in potential alone cannot continuously deliver work.
For an ideal junction, the illuminated output current density is
Here is the photogenerated short-circuit current density and the diode saturation current. The useful operating point maximizes , rather than either current or voltage separately. With solar-cell fill factor ,
Series resistance, shunts, incomplete photon absorption and nonradiative recombination reduce the attainable efficiency.
Figure 1.
Ideal illuminated solar cell: current, delivered power and the maximum-power point
. Details.
The Shockley-Queisser limit removes these avoidable losses but retains the following physical constraints. There is one semiconductor junction with one band gap; each absorbed above-gap photon yields at most one collected electron-hole pair; carriers thermalize to the cell temperature; and the inverse optical process produces unavoidable radiative recombination. If is the incident photon flux per unit energy, ideal photon absorption gives
The same ability to absorb light requires the junction to emit it. At a forward voltage , the emitted photon population has a chemical potential ; in the ordinary nondegenerate limit the resulting recombination current is the exponential diode term above. Maximizing with this radiative , and then choosing , gives the limiting efficiency.
A small band gap collects many photons but gives a low voltage and large thermalization losses. A large band gap permits a larger voltage but transmits too much of the incident spectrum. Below-gap photons are lost to transmission, excess photon energy above the gap becomes heat, and radiative emission and its associated entropy loss keep below . The optimum single-junction conversion is roughly thirty percent under ordinary unconcentrated sunlight; the precise ceiling depends on spectrum, cell temperature and optical geometry. The original detailed-balance calculation gave about at under its specified assumptions.[ref] Multiple junctions and optical concentration change the assumptions and consequently change the bound.
Use two plane waves, and , as a basis near their Bragg degeneracy. The potential has Fourier coefficients ; its zero Fourier coefficient vanishes. All other plane waves are off resonance and are neglected in the leading nearly-free electron model. With , the Schrödinger equation reduces to the Hermitian matrix problem
Setting the characteristic polynomial equal to zero gives
Solving this quadratic equation gives
For the real cosine potential its coefficient is real; the absolute-value form also allows an arbitrary choice of origin and corresponding Fourier phase.
At the Brillouin-zone boundary, , the free energies are equal. The two standing waves have different probability densities relative to the lattice potential and are split by
Repeating this avoided crossing at the reciprocal lattice boundaries joins smooth branches into energy bands. In one dimension a nonzero Fourier coefficient opens a gap at the corresponding crossing. In higher dimensions a direct gap at one wavevector need not give a common forbidden energy interval at every wavevector: different bands may still overlap in energy. The present two-dimensional example remains metallic even though the crossings are locally gapped.
For the rectangular lattice, choose real-space primitive vectors and . The reciprocal lattice vectors are
Define the half-widths of the first Brillouin zone by and . Their values and the zone are
A primitive reciprocal cell can be drawn with a corner at the origin; its centred translate is the first Brillouin zone here. Its area is .
Figure 1.
Real-space primitive cell and reciprocal-space first Brillouin zone
. Details.
Including both spin states, counting the occupied free-electron states per unit area gives
Consequently
The occupied disk has precisely the area of the first Brillouin zone, but their shapes differ: . It protrudes through both vertical zone boundaries and does not reach the horizontal ones. Thus two electrons per cell do not imply a filled first band in the absence of a separating band gap.
Fold the protruding caps by . The original circle and the translated circle intersect at , . The potential mixes these states and replaces each crossing by an avoided crossing. The lower branch reconnects into two open sheets running across successive zones in the direction; the upper branch forms a closed lens-shaped pocket about a vertical zone boundary. In a reduced-zone sketch its two halves sit at opposite edges, but periodic identification makes them one closed sheet. This is rectangular-lattice Fermi-surface reconstruction.
Figure 2.
Free and reconstructed Fermi contours: two open sheets and one closed pocket per reciprocal cell
. Details.
The largest pocket in this weak-potential construction is the lens-shaped intersection of two disks enclosed by the intersecting free circles. Each has radius and their centres are separated by . If , two sectors minus their two triangles give the lens bound on a reconstructed Fermi pocket:
Opening a gap raises the upper branch at fixed chemical potential and shrinks its pocket. The reported quantum oscillation area satisfies
It is therefore consistent with this closed sheet. Onsager quantization gives . An orbit on this closed contour can return to itself under the Lorentz force and occupy quantized Landau levels. The open sheets do not enclose a finite orbit area for a perpendicular field in this two-dimensional model. The observation is consistent with the reconstructed pocket, without uniquely determining a completely arbitrary band structure.
For a quantitative rough coupling estimate, use the supplied fixed-chemical-potential approximation and retain the two waves coupled by . Write
Their upper eigenvalue is
Setting it equal to gives the pocket boundary
Its intersection with occurs at , where
The two-wave estimate of a Bragg-gap coupling from a pocket area therefore uses
At this reproduces the lens bound on a reconstructed Fermi pocket; the pocket disappears at . Numerical integration gives and . Matching the measured area gives , so the justified rough result is
This estimate neglects remote plane-wave admixtures and other Fourier coefficients. It also holds the chemical potential fixed as requested; an exact fixed-density calculation would determine it self-consistently after reconstruction. The numerical root is a way to implement the approximate model, not a measurement of to four significant figures.
Dopant atoms and carrier type. Substitutional semiconductor doping places a small fraction of foreign atoms on host lattice sites. In a tetrahedrally bonded group-IV semiconductor, a group-V substitution such as phosphorus in silicon has an extra valence electron. A shallow donor level just below the conduction edge can release that electron, leaving a positively charged donor. In an n-type semiconductor, mobile electrons are the majority carriers and the chemical potential shifts toward the conduction band.
A group-III substitution such as boron in silicon has one fewer valence electron. An acceptor level just above the valence edge captures an electron from the valence states, leaving a mobile electron hole and a negatively charged acceptor. In a p-type semiconductor, electron holes are the majority carriers and the chemical potential shifts toward the valence band. The material can remain electrically neutral in both cases: mobile carrier electric charge is compensated by fixed ionized impurity electric charge. Dopant freeze-out, compensation and degeneracy qualify the simplest one-carrier-per-dopant picture.
Making the junction. Choose p material at and n material at . Electrons initially diffuse from n to p, and electron holes from p to n, recombining near the interface. Negative ionized acceptors are exposed on the p side and positive ionized donors on the n side. They produce an electric field directed from n to p. The field opposes further majority-carrier diffusion until drift and diffusion balance and the equilibrium electrochemical potential is uniform.
The depletion regions are depleted of mobile majority carriers, not of the fixed dopant charges that create the space charge. At equilibrium the conduction and valence edges bend together, while the chemical potential is flat. The specified is a electric potential difference: the corresponding initial chemical-potential energy separation is . This convention is required by the quoted units and width formulas. Write for the absolute permittivity; the paper's is relative permittivity.
Figure 1.
Substitutional donor/acceptor levels, equilibrium band bending, fixed space charge and the junction electric field
. Details.
Poisson equation and widths. Let the depleted intervals be and . Fully ionized number densities correspond to electric charge densities , respectively. Thus
The factor converts a dopant concentration to a space-charge density. Requiring zero field in both neutral regions gives charge neutrality,
With , direct integration gives the triangular field
Continuity at the interface is the same neutrality condition. The field is negative throughout the depleted interval. Set . Its integral gives
Potential continuity therefore requires
Combining it with neutrality yields the abrupt-junction depletion width and capacitance results
The more lightly doped side depletes farther, because equal and opposite space charges must be exposed.
Applied bias. Define positive as forward bias, raising the p-side terminal potential relative to the n-side terminal. It lowers the remaining barrier from to . Reverse bias has and raises that barrier. The same Poisson integration applies, with this replacement:
Forward bias narrows the depleted region; reverse bias widens it. Under bias there need not be one equilibrium chemical potential across the device, so the equilibrium flat-level argument must not be reapplied unchanged. The barrier/width change is illustrated in the next figure.
Differential capacitance and electric charge sign. Put and . The signed p-side fixed electric charge per unit area is
Differentiating with respect to forward electric potential difference gives the positive capacitance per unit area
The n-side electric charge is , so instead. If denotes a positive space-charge magnitude, the positive capacitance is for this forward-voltage convention. This resolves the sign implicit in the paper's positive . The total device capacitance is , where is its cross-sectional area.
Finite device and numerical sketch. For the symmetric device, and each side has thickness . Before full depletion,
Using , the zero-bias width on each side is , less than . The total zero-bias capacitance is .
Full depletion occurs when , so
At larger reverse bias the depleted thickness cannot exceed . In the ideal fully depleted slab, incremental terminal electric charge changes the field uniformly through that fixed thickness, giving the full-depletion capacitance of a symmetric junction slab
Thus the idealized piecewise sketch over the specified range is
It is continuous at full depletion, rises through at zero bias, and has a formal inverse-square-root divergence as from below.
Figure 2.
Bias-dependent depletion geometry and the numerical capacitance curve, including full depletion
. Details.
There are two physical limits to this extension. Once all dopants are exposed, their fixed space charge no longer increases, but the terminal electrodes can still acquire incremental electric charge; differentiating only the fixed ionized electric charge would incorrectly predict zero capacitance. Near , the depletion approximation itself fails as mobile carriers and forward injection become important, so the formal infinite endpoint is not a prediction of an actual infinite measured capacitance. The fully depleted plateau assumes ideal contacts and no breakdown, injection storage or extra series capacitance.
Why degenerate perturbation theory is needed. If several unperturbed states have equal or nearly equal energy, a perturbation can mix them at leading order. In that subspace, diagonalize the perturbation matrix rather than using nondegenerate formulas with small energy denominators. The resulting eigenvectors identify the combinations selected by the coupling, and the eigenvalues give the leading splitting. Other remote states can then be included perturbatively.
For the linear combination of atomic orbitals description of molecular hydrogen, equivalent atomic states on the two nuclei supply the nearly degenerate starting functions. In an orthonormal two-state model, equal diagonal energy and real hopping give and symmetric/antisymmetric combinations. For actual overlapping atomic functions with overlap and Hamiltonian coupling , one instead solves a generalized eigenvalue problem:
The bonding state has enhanced density between the nuclei and the antibonding state a separating node. Two opposite-spin electrons occupying the bonding molecular orbital describe the simplest closed-shell hydrogen-molecule bond. Electron-electron interaction and nuclear repulsion are needed to calculate a realistic total binding energy; the two-state splitting alone is not the full molecular energy.
In the nearly-free electron model, the analogous degenerate states are plane waves and at a Bragg scattering. A Fourier coefficient of the periodic potential couples them. At exact degeneracy , diagonalizing their two-by-two matrix gives
The associated standing waves place electronic density differently relative to the ionic potential. These are the same gap-opening combinations illustrated in Q1(b). The labels of the two starting basis states survive as mixed orbital or plane-wave character, rather than independent crossing eigenstates.
Origin of broad and narrow bands. Spatially extended s/p-derived orbitals have substantial overlap and therefore large hopping amplitudes and broad bands. More localized d orbitals overlap less strongly and produce narrower bands. A multiorbital tight-binding calculation assigns the on-site orbital energies and hopping matrices, forms Bloch sums, and diagonalizes the resulting matrix at each wavevector. Different hopping ranges and orbital symmetry determine the bandwidths and allowed interband coupling. The two-band description retains one broad combination and one nearly flat d-derived combination; a real transition metal contains several d orbitals and more complicated dispersions.
Hybridized eigenvalues and eigenstates. Let , in an orthonormal basis. Hermitian operator fixes the reverse coupling to be . Projection onto the two components gives
Solving this quadratic yields the band hybridization dispersion
For nonzero , one eigenvector representation is , followed by normalization. Thus the wavefunctions exchange broad-band and d-band character through the avoided crossing; the levels are equal mixtures at the uncoupled crossing.
At , the parabola crosses the horizontal level at , with . For small nonzero , the curves avoid crossing with minimum direct separation . The lower branch follows the parabola below the crossing and flattens toward above it; the upper does the opposite. The lower branch satisfies everywhere and tends to from below as , while the upper branch always lies above . This is the perfectly flat-band continuum approximation, with zone-boundary intersections excluded as requested.
Figure 1.
Uncoupled and weakly hybridized broad/flat bands, with the shifted Fermi radius and small lower-band slope
. Details.
Fixed particle number and the larger Fermi sphere. Denote the uncoupled broad-band Fermi radius by to distinguish it from the new one . With two spin states, the broad-band number density at is . The stipulated equal population in the flat sector makes the total density twice that value. At zero temperature and any nonzero hybridization, all occupied states fit into the monotonic lower branch below , while the upper branch is empty. Therefore
This is the stated factor increase. It uses fixed total particle number, equal spin counting, an exactly flat unhybridized band and enough available states below any Brillouin-zone restriction. In generic pairs of dispersive bands both hybridized pockets can remain occupied; then Fermi-volume conservation in two hybridized electron bands fixes a sum of volumes rather than this one radius. The exactly flat band is degenerate at , so its initial occupation must be specified, as it is here.
Fermi-velocity suppression. Set . The exact derivative of the lower branch is
At , and . For , the square bracket divided by two is . Since group velocity is and the broad-band reference velocity is ,
so the flat-band hybridization Fermi-velocity suppression is
The reference velocity is that of the original broad-band Fermi sphere, not the broad-band slope at the larger ; omitting the radius change loses the factor .
The low-temperature qualification becomes more stringent as : the new chemical potential approaches with , and the relevant thermal energy must be small compared with this hybridization scale. At fixed nonzero temperature the singular flat-band redistribution is not automatically resolved for arbitrarily small coupling. A small but finite narrow-band dispersion can also change the occupation topology.
Density of states and magnetism. For an isotropic three-dimensional branch with spin degeneracy two,
A small Fermi velocity therefore produces a large density of states. The localized, partially filled d-derived bands in iron, cobalt and nickel place many states close to the Fermi level; broad s/p-derived conduction bands in copper, silver and gold coexist with filled d bands lying below it. This explains the substantially larger Fermi-level density of states of the former group within the band picture. An exact factor of ten requires the actual band parameters and is not fixed by the two-band toy model.
A high density of states lowers the kinetic-energy cost of spin polarization. A repulsive exchange interaction can then satisfy the Stoner criterion and stabilize ferromagnetism. With a single-spin density of states and an on-site interaction normalization, the criterion is ; with spin-summed , it is . The large d-state contribution favors this instability in Fe, Co and Ni. The noble metals do not satisfy the corresponding instability condition in the same ordinary state; this does not imply that they have no weak magnetic response whatsoever.
Nondegenerate carrier product. Assume three-dimensional parabolic energy bands near their extrema, spin degeneracy two, and no additional uncounted valley multiplicity. Near the conduction edge, a Boltzmann approximation to the Fermi-Dirac distribution gives
Angular integration and the Gaussian radial integral yield
For valence electron holes, the empty-state probability is . Measuring electron hole kinetic energy downwards from similarly gives
The factors two count spin once in each band. Multiplying cancels the chemical potential and proves the semiconductor mass-action law:
In an intrinsic neutral semiconductor, . In a doped equilibrium semiconductor the same product holds under the stated nondegenerate assumptions, even when . Extra valleys or several valence branches must be included in the effective density-of-states factors rather than counted inconsistently.
Equilibrium junction. Put the p region on the left and the n region on the right. Before contact their majority carriers are electron holes and electrons, respectively. Diffusion initially sends electrons from n to p and electron holes from p to n. This leaves negative ionized acceptors on the p side and positive ionized donors on the n side, creating a depletion region and an electric field directed from n to p. At equilibrium this field balances diffusion, so the net electric current vanishes.
For a homojunction the gap remains constant while the conduction band and valence band edges bend downwards from p to n by . The equilibrium chemical potential is flat. Since , the electrostatic potential rises from p to n. The carrier profiles satisfy
With fully ionized nondegenerate dopants, the neutral p side has , , while the neutral n side has , . Thus the electron concentration rises across the junction and the electron hole concentration falls, with both small relative to the far-side majority densities in much of the depleted region. The ideal depletion approximation neglects these mobile populations when calculating the fixed space electric charge, but the actual equilibrium profiles are smooth.
Figure 1.
Equilibrium band bending, a flat chemical potential and complementary electron/hole profiles across a homojunction
. Details.
The built-in potential follows from these Boltzmann factors:
There is no applied terminal bias; the internal electrostatic potential variation does not contradict a constant electrochemical potential.
Opposing carrier fluxes under reverse bias. Use for the magnitude of reverse bias and for reverse electric current, as in the PDF. The junction barrier increases by . The Shockley diode equation in that orientation can be written
The reverse-directed contribution is supplied by minority electrons on the p side and minority electron holes on the n side. They reach the depletion edge by diffusion and are swept across by the built-in field. The opposing contribution is associated with majority carriers thermally crossing the barrier in the forward direction; reverse bias suppresses it by the Boltzmann factor. At zero bias detailed balance makes the two equal, while at large enough reverse bias the opposing flux becomes negligible.
In the ideal neutral-region diffusion model,
Here are minority-carrier diffusion coefficients, their diffusion lengths, and is junction area. The two terms inside this positive prefactor are the electron and electron hole contributions; both flow in the same reverse-current direction. This separates the two carrier-species contributions from the two opposing fluxes in the Shockley diode equation. Calling the saturation term a drift current refers to its passage across the depleted region, and does not eliminate the minority-carrier diffusion needed to supply it.
Reading the measured curves. At , the ideal equation gives , initially proportional to reverse electric potential difference. At , it approaches a plateau . The thermal electric potential differences at the three temperatures are about –, consistent with a rise at small reverse biases and near saturation at larger ones. Higher temperature produces a much larger plateau because thermal electron-hole excitation raises exponentially.
The lower-temperature curve has a noticeable residual electric potential difference dependence rather than a perfectly flat plateau. Depletion-region generation, surface leakage or other departures from the ideal diffusion model can contribute. The sharp upturn at the largest reverse electric potential differences is reverse breakdown, outside the unextended Shockley diode equation; the plot alone does not identify a unique microscopic breakdown mechanism. The original horizontal axis duplicates the final decade label, so a precise breakdown electric potential difference cannot be read reliably. Neither issue prevents estimating the gap from the intermediate-bias saturation region.
Band-gap estimate. Use plateau electric currents at a common intermediate reverse bias, away from both the small-bias rise and the breakdown upturn. Approximate readings are , –, and . Kelvin temperatures are required in the activation law.
For fixed effective masses, . If the remaining transport prefactor is approximately constant, extracting a semiconductor band gap from diode saturation current therefore uses
Using the widest temperature separation and the rounded readings gives
Neglecting all algebraic prefactors and fitting instead gives about . These differ by roughly , a material correction on the graph's precision scale. Thus an honest graph-based conclusion is
The nonflat low-temperature plateau, rough electric current readings, unprovided temperature dependence of , and possible gap variation prevent a more precise extraction. If depletion-region generation dominated rather than the specified mechanism, its activation slope could involve a different energy scale; the displayed estimate explicitly uses the diffusion-dominated saturation model requested here.
Assume three-dimensional parabolic band edges, with density-of-states masses and one valley per edge, each with spin degeneracy two. Their effective densities of states are
For the semiconductor mass-action law, integrating the square-root density of states against the Boltzmann tail gives and . Multiplication cancels the chemical potential, giving
If there are additional valley degeneracies, multiply by their product; anisotropic bands use the appropriate density-of-states masses. In the later equal-mass model, . No numerical mass or valley multiplicity has been invented.
When , the conduction-edge occupation is exactly , so the nondegenerate approximation for electrons fails at every positive temperature. Define the Normalized complete Fermi-Dirac integral
For the carrier densities with chemical potential pinned to a band edge, with the two specified equal masses, put . Direct Fermi-Dirac integration then gives
Thus exactly within the parabolic model. At low temperature, is thermally activated; at , from below. Both vanish at zero temperature and increase at fixed chemical potential. The other states that pin can supply the compensating charge; one must not additionally impose .
For the Pauli reduction of the pinned-edge mass-action product, the exact product relative to the classical law is
It approaches at low temperature and in the formal high-temperature limit. The low-temperature activation factor is unchanged but its prefactor is reduced; the classical mass-action formula overestimates the product. The high-temperature limit here refers to the given unbounded parabolic bands, not a claim that real semiconductor band structures remain unchanged at arbitrarily high temperature.[ref]
Finally, for a fixed positive , the zero-temperature Fermi sphere has radius , so
The low-temperature expansion at fixed gives the fixed-chemical-potential Sommerfeld density increase:
The monotonic density at fixed chemical potential in a square-root density of states initially rises quadratically from . For , its leading behaviour crosses to . The fixed chemical potential is what allows the particle number to change.
Figure 1.
Pinning the chemical potential to or above a band edge changes carrier populations and the mass-action product
. Source. Details.
Specify the spin counting first. Let denote the unpolarized single-spin density of states per atom per energy, and put . At fixed total occupation , the two Hubbard model mean-field energies can be written
The common shift can be absorbed into the chemical potential. Lowering the up-spin band by adds carriers, while raising the down-spin band removes the same number. Thus
If is the atom number density, the magnetization is . The Stoner-enhanced magnetic susceptibility in SI volume units is therefore
Here is the noninteracting Pauli paramagnetism result. If the paper's counts both spins, and the same result is
A susceptibility per atom omits . Stating the spin and volume conventions avoids an otherwise spurious factor of two or a dimensional mismatch.
The Stoner criterion is , equivalently with spin-summed . To see the instability, transfer electrons per atom from one spin band to the other. The band-filling cost is , whereas the interaction energy changes by . Their sum becomes negative at the threshold. The unpolarized state then ceases to be a local minimum, and nonlinear band filling determines the spontaneous ferromagnetism. The negative continuation of the displayed linear susceptibility beyond the threshold is not a stable paramagnetic response.
Band edges and optical transitions. A direct band gap has the conduction minimum and valence maximum at the same crystal momentum. A near-edge photon then connects the two almost vertically on an energy-versus-wavevector plot, because its momentum is small on the Brillouin zone scale. An indirect band gap has the extrema at different momenta, so an additional momentum transfer, commonly supplied by phonon absorption or emission, is required for a near-edge radiative transition. This reduces its probability relative to an allowed direct transition.
Near the zone centre, an isotropic direct-gap model has and downward valence dispersions . The heavy-hole and light-hole bands are the upper branches with large and small positive hole masses: the heavy branch is flatter. In common cubic semiconductors these arise from the manifold, while spin-orbit coupling places the split-off valence band lower by . Its thermally generated hole population is suppressed by approximately relative to upper-valence holes, apart from density-of-states factors. Neglect is justified when ; it is not a universal exclusion of this band at arbitrary temperature.
Figure 1.
Original schematic direct and indirect band edges, including heavy, light and split-off hole branches
. Details.
For an indirect material such as Si, the conduction minimum need not be at . In applying the carrier-density formulae, must mean the actual minimum of the included conduction valley. The supplied scalar effective masses and densities of states constitute a simplified band-edge model; they do not specify Si's full valley anisotropy or multiplicity.
Carrier populations. Let be the Fermi-Dirac distribution. In the nondegenerate conduction tail, . Set and integrate the given density of states:
The integral is . Defining the effective density of states gives
The factor two is already included in the supplied density of states; it must not be added a second time.
A hole is an unoccupied valence state. For sufficiently below , . Substitution of gives in exactly the same way
Their product eliminates the chemical potential:
This is the semiconductor mass-action law at thermal equilibrium. In an intrinsic semiconductor, charge neutrality additionally gives
In a doped semiconductor the product still applies in this nondegenerate equilibrium model, but the two populations need not be equal.
Intrinsic chemical potential. Taking logarithms of yields the intrinsic-semiconductor chemical potential
For unequal band degeneracies, use the first form with the full . For the supplied mass ratio , the shift from midgap at is approximately , toward the valence edge.
Three temperature regimes. The original plot runs from high temperature on the left to low temperature on the right. At high temperature, thermal excitation across the gap dominates and rises rapidly with temperature. In the middle, donors and acceptors are effectively ionized, the material is extrinsic, and makes the plateau. At low temperature, dopant freeze-out binds increasing numbers of electrons to donor levels, reducing the mobile-electron population as the sample cools. Its activation scale is a donor binding energy, not the full band gap. Compensation and incomplete ionization affect that regime.
Gap estimate from the original figure. On the dashed intrinsic branch, . Hence band-gap extraction from intrinsic carrier activation uses the slope of versus . With and expressed in eV/K, two points give
The very compressed intrinsic branch is only approximately readable: near and near . These give , rather than a precision determination. Neglecting the thermal prefactor gives about from those same rough points. Thus a gap of order – is justified by this schematic plot; its resolution does not support distinguishing a few hundredths of an electronvolt. The declining freeze-out branch must not be used to infer the full gap.
Room-temperature densities. At , the supplied masses give
In the fully ionized extrinsic plateau, put . Charge neutrality and the mass-action law imply
Since , is insensitive to the imprecise gap. The minority density is much more sensitive:
Using gives and . Using instead gives and . These are explicitly conditional estimates within the stated mass model. The plot alone cannot fix a precise minority density because the uncertain gap occurs inside an exponential.
Noninteracting spectral line. Treat as a wavevector, so physical crystal momentum is . Write ; an energy zero at the chemical potential can be chosen throughout. The temporal factor of a noninteracting state is . With the Fourier transform convention ,
Thus its normalized frequency spectral function is
One should not square the delta distribution as if it were an ordinary finite-time transition amplitude. More precisely, the fermionic retarded Green function is . Its Fourier transform is , whose negative imaginary part divided by is the same delta peak.
Interactions and lifetime. Interactions change the dispersion through the real part of the self-energy and permit decay into other excitations through its imaginary part. In the simplest quasiparticle-pole model, exponential decay of the retarded amplitude gives
Performing the causal integral,
The corresponding Lorentzian quasiparticle spectral peak is
Its total frequency integral is one. The exponential amplitude lifetime is ; the probability lifetime is . Describing as a scattering rate therefore requires the lifetime convention to be specified. The energy half-width is .
A general interacting electron is more accurately described by near the pole. The coherent peak then integrates to the quasiparticle residue and an incoherent continuum carries the remaining weight. In an energy-unit self-energy convention,
Linearizing the real denominator gives and at the pole. The requested unit-weight Lorentzian is the normalized quasiparticle model, not the most general exact interacting spectrum.
Peak position and width. At fixed , put and . Differentiation gives
It is positive for and negative for , so the only maximum is at , with height . Setting to half this height gives . Therefore
In energy units the maximum is and the half-width is . Frequency dependence of the self-energy would require reconsidering both the shape and peak position.
Fermi gas and Fermi liquid. A noninteracting Fermi gas has a delta-function spectral line at its bare dispersion for every momentum. A normal three-dimensional Fermi liquid has a shifted coherent peak plus incoherent spectral weight; close to the Fermi surface its width scales as . At zero temperature the peak becomes arbitrarily sharp at the Fermi surface, while it broadens farther away. Spectral functions include addition and removal weight; occupied-only photoemission contains an extra Fermi factor.
Figure 1.
A free spectral line and broadened quasiparticle peaks, with a separate incoherent background
. Details.
Photoemission measurement. A photon of known energy ejects an electron. Measure its kinetic energy and emission direction. For binding energy below the Fermi level and the appropriate calibrated work function , conservation of energy gives
Conservation of momentum parallel to a periodic surface gives
where is a surface reciprocal-lattice vector. At suitable photon energies the photon momentum can be neglected; at higher energies it must be retained. The surface breaks perpendicular translation symmetry, so the perpendicular crystal momentum generally needs a final-state model. A two-dimensional band can be reconstructed directly from the parallel components.
Within the sudden approximation and electric-dipole approximation, photoemission intensity and spectral function are related by
Here is measured relative to the chemical potential. Instrumental resolution convolves the product and inelastic scattering adds background. Thus ARPES maps the occupied electron removal spectral function, weighted by a matrix element; it does not measure the full unoccupied spectrum or an unweighted spectral function automatically. Bright ridges give dispersions, and intrinsic widths give decay information after instrumental effects are removed.
Reading the supplied image. Its bright ridges form a graphene-like cone. The crossing near is the Dirac point, not in general the Fermi level. In the paper's plotted energy coordinates its energy is roughly –. A velocity follows from a branch slope independently of the energy zero:
For example, the lower branch moves by about between and . This gives and
Reading other points changes the estimate by tens of percent. The original exam image has unequal energy increments on nearly equally spaced labelled ticks, so it does not justify a precise fitted velocity. The value is an estimate of the cone slope extrapolated to the Fermi crossing, rather than a measured derivative exactly at a visibly identified Fermi cutoff.
To determine a Fermi energy, first calibrate the Fermi edge, then subtract the chosen reference energy. Neither a Fermi-level marker nor an explicit energy-zero convention is supplied in this crop. If its upper termination near is intended to represent the Fermi cutoff, the conditional graph answer is on that axis, or relative to the Dirac crossing. Identifying the crossing itself as the Fermi energy would incorrectly assume undoped graphene.
The research article cited by the exam identifies this as linearly polarized graphene ARPES and independently states for the substrate-doped sample. That resolves the relative energy when this external calibration is admitted.[ref] The crop alone does not uniquely fix an absolute Fermi energy; the source-supported relative value and the conditional axis reading should be distinguished.
At equilibrium, majority electrons undergo diffusion from n to p, while majority electron holes diffuse from p to n. Their departure exposes positive donors on the n side and negative acceptors on the p side. The resulting electric field points from n to p and opposes further diffusion. Equilibrium is reached when the chemical potential is constant across the junction. The built-in potential is this electrostatic adjustment, not a discontinuity of the material's band gap.
Since electronic band energies shift by , both the conduction band and valence band undergo band bending down from p to n by at zero applied bias. Their separation remains the same gap. In the nondegenerate thermal equilibrium limit the built-in voltage can also be expressed using neutral-edge densities as .
Take the n-contact potential and its electrochemical potential as reference. Applying to the p side changes the total depletion-layer drop to
Thus is reverse bias: it increases the barrier and widens the depletion region. The bands have horizontal tangents at the neutral edges in this approximation and a smooth cubic bend within the linearly graded layer.
A biased junction is not in global thermal equilibrium. Its contact electron electrochemical potentials obey
Under reverse bias the p-contact level is higher by . A single horizontal equilibrium chemical potential cannot describe this state; electron and hole quasi-Fermi levels describe the nonequilibrium populations. The sketch marks the contact levels. Detailed spatial quasi-Fermi profiles require transport and recombination information beyond the electrostatic doping profile.
Figure 1.
Band bending, electric field and capacitance of a linearly graded junction
. Original depletion-model sketches with the p region on the left. The top panels compare equilibrium with reverse bias; under bias only the contact electrochemical levels are specified. The bottom panels show the negative electric field and the normalized positive differential capacitance. Energy and length scales are illustrative. Details.
Let the common uncoupled Fermi wavevector be . Equal carrier densities imply equal Fermi surface volumes, and hence equal radii for the spherical spin-degenerate electron pockets. The conserved total density is
If hybridisation leaves only one occupied electron-like band, state counting gives
This is the familiar single-large-Fermi-surface answer. It requires the other hybridised band to be empty.
For finite heavy-band mass, that conclusion does not follow from an arbitrary nonzero coupling. At sufficiently small , the upper band still has an occupied pocket, so there are two Fermi momenta, not one. The exact fixed-density statement is Fermi-volume conservation in two hybridized electron bands:
For an explicit weak-coupling check, write , , with , and enforce the uncoupled equality at : , . Put and . The Fermi-surface equation is
At fixed density , giving
For arbitrary coupling, put . The roots of the same characteristic equation, before discarding unoccupied pockets, are
Choose the chemical potential so that . This specifies the Fermi momentum of each pocket for any supplied curvatures and coupling.
The transition to one occupied pocket occurs at . Set ; at the transition the two squared radii are and zero. Their sum and product give
For there are two pockets; for the only nonzero radius is . The outer radius grows continuously and the inner one shrinks. The nearly flat band limit can make depopulation occur at a small coupling, but the mass ratio or a single-surface assumption is needed to turn the question into a unique finite-coupling numerical change.
In the nondegenerate conduction band, the Fermi-Dirac distribution reduces to its Boltzmann distribution tail. Thus
Here the substitution leaves . For electron holes use the empty-state probability , not the conduction-band Boltzmann factor. This gives
The chemical potential cancels in the product. With ,
This semiconductor mass-action law holds in thermal equilibrium for nondegenerate carriers, including a doped semiconductor. Only in an intrinsic semiconductor must . The factors are the effective density of states; the factor in the final expression follows from the two spin factors already included in (a).
Use the reduced attractive reduced BCS pairing Hamiltonian with volume V, and attraction :
A specified pairing shell or an ultraviolet renormalization is part of the model. The trace is grand canonical. Insert fermionic coherent states into its time-sliced trace. The resulting Grassmann variables are antiperiodic, and with inverse-energy imaginary time the dimensionless action is
This boundary condition is essential for Fermi-Dirac statistics; replacing it by a bosonic periodic trace changes the answer. A complex Hubbard–Stratonovich transformation uses the local identity
with time-slice factors understood. The pairing field renders the fermion action quadratic. For a uniform static field the independent Nambu spinor block has
Its determinant is . The phase-independent determinant sign cancels under consistent reference normalization. Integrating the fermions gives the BCS pairing-field effective action, , with the normal-ordering/reference term retained. The fermionic pairing-block determinant normalization uses one independent pairing block for every k; a fully doubled spin-particle-hole trace instead needs one half. The subtracted grand potential of the reduced BCS model is particularly unambiguous:
Performing the Matsubara frequency sum yields
up to a common constant. At zero gap this reduces exactly to the two-spin free-gas result. The Matsubara derivation of the BCS saddle follows by differentiating with respect to and using gives
The second equation fixes the chemical potential at specified number. This stationary action is BCS mean-field theory, not an exact integration of pairing-field fluctuations. With constant single-spin density of states per volume and a symmetric energy cutoff , weak coupling gives and . The zero-temperature condensation grand-potential density is to leading weak-coupling order; using a spin-summed density would change these displayed factors. The BCS coherence factors are and , and the associated BCS ground state is the normalized product of empty and occupied pair states. [ref]
Allowing a slowly varying complex pairing field generates phase stiffness and, near the transition, the Ginzburg-Landau theory of superconductivity. A neutral paired system has a gapless phase sound mode; in a charged bulk superconductor electromagnetic coupling produces the Anderson-Higgs mechanism in a superconductor and the Meissner effect. These collective fluctuations are distinct from the gapped fermionic excitations. A microscopic dynamical calculation is needed to fix collective-mode dispersion and damping; the static saddle alone does not determine them. [ref]
For a normal-ordered bosonic Hamiltonian, insert coherent-state resolutions between short imaginary-time evolution factors. The overlaps generate and the Hamiltonian symbol is obtained by replacing creation and annihilation operators by and . Including the chemical potential gives, in the stated convention,
The trace imposes periodic imaginary-time boundary conditions,
so bosonic Matsubara frequencies are . The two complex coherent-state integration variables are treated independently before imposing the physical saddle contour. There is no extra zero-point ordering term: the supplied interaction is already normal ordered.
The spatial confinement is not specified beyond the volume. The uniform bulk results below use a homogeneous box with periodic spatial boundaries, or the interior of a large vessel where boundary effects can be neglected. A hard wall instead requires a vanishing field at that wall and an inhomogeneous saddle near it. Temporal periodicity follows from the thermal trace and does not by itself impose periodic spatial boundaries.
The two-mode approximation for a Bose gas retains the real localized orbitals:
It is valid when the lowest doublet is well separated from higher energy levels, and the temperature and interaction-induced changes of the orbitals are small on that separation. More explicitly, if is the energy gap to discarded orbitals, sufficient low-energy conditions are , where . A deep double-well potential suppresses interwell interaction overlaps. Large alone does not justify the truncation: increasing at fixed coupling can violate the interaction condition.
Let and use the even and odd eigenfunctions to evaluate the projected one-particle Hamiltonian operator:
Thus
The same quantum tunnelling matrix element is . Reflection symmetry makes . Neglecting mixed contact interaction integrals leaves
The last step is the normal-ordered boson pair count, following from the canonical commutation relation. Consequently the two-site Bose-Hubbard model has precisely the requested coefficients
and
Here is the coefficient obtained by rewriting the fixed-particle number Hamiltonian operator, rather than an independently prescribed thermodynamic chemical potential. Adding a reservoir term would replace it by .
The negligible-overlap instruction applies to the mixed contact interaction integrals while the small quantum tunnelling matrix element is retained. Exactly disjoint smooth orbitals for a local would also give ; setting every overlap-related matrix element identically to zero would discard the physical splitting one is trying to describe.