A new constituent of electrostatic energy in semiconductors An attempt to reformulate electrostatic energy in matter

. The concept of electric energy is revisited in detail for semiconductors. We come to the conclusion that the main relationship used to calculate the energy related to the penetration of the electric ﬁeld in semiconductors is missing a fundamental term. For instance, spatial derivate of the electrostatic energy using the traditional formula fails at giving the correct electrostatic force between semiconductor based capacitor plates, and reveals unambiguously the existence of an extra contribution to the standard electrostatic free energy. The additional term is found to be related to the generation of space charge regions which are predicted when combining electrostatics with semiconductor physics laws, such as for accumulation and inversion layers. On the contrary, no such energy is needed when relying on electrostatics only, as for instance when adopting the so-called full depletion approximation. The same holds for neutral and charged insulators that are still consistent with the customary deﬁnition, but these two examples are in fact singular cases. In semiconductors for instance, this additional energy can largely exceed the energy gained by the dipoles, thus becoming the dominant term. This unexpected result clearly asks for a generalization of electrostatic energy in matter in order to reconcile basic concepts of electrostatic energy in the framework of classical physics.


Introduction
Interpretation of electric energy in conjunction with thermodynamics has been widely investigated, with a special interest for dielectric bodies and ideal conductors [1][2][3][4][5]. The electric energy stored inside of a body can be expressed whether in terms of charges and potentials restricted to the volume of the body [5], or in terms of fields including contributions beyond the physical boundary of the system, see relations (1) and (2) for linear polarisable systems [5]. In addition, careful considerations and exhaustive criticisms about the validity of electrostatic energy formulation in conductors and insulators has been addressed in references [1,2]. Beyond the thorough literature, instructive and complete analysis of thermodynamics of electric and magnetic fields has also been developed in references [1,2,4,6].
While the topic of electric energy in matter seems to be well established and widely accepted, we will bring some evidence that this paradigm must be revisited in semiconductors to say the least, resulting in a new contribution to the well-known Helmholtz free energy arising from the influence of the electric field.
Before discussing an instructive virtual experiment, we revisit some fundamental relationships dealing with eleca e-mail: jean-michel.sallese@epfl.ch trostatic energy from a general thermodynamic point of view.
The electrical work U E that must be spent to gather charges from infinity into a volume Ω is given by [3,5,6]: where ρ is the local charge density and φ is the electric potential (note that the potential φ must vanish at infinity [5]). The integral is limited to the volume of the body containing the charges, and in this sense, relation (1) represents the electric energy of the content of Ω (assuming the system linear). The internal energy of electrical nature is then expected to be implicitly contained in (1), which is indeed how Frankl [7] analyzed the free energy stored in the depletion region of a silicon layer. This formulation attributes energy to electric charges. Alternatively, adopting the electric field and displacement vectors concepts (E = −grad(φ), div(D) = ρ, bold letters hold for vectors), the electrostatic energy can also be expressed from the electric field generated by the charges enclosed in the volume Ω, provided that the integration is performed over the whole space Ω ∞ , including matter: Again, the electric field and displacement vectors in (2) do not account for all sources in the universe; they are only assigned to the charges located in Ω. It can be shown that equations (1) and (2) are equivalent and represent the electrostatic energy of those charges [5]. A generalisation of relation (2) to nonlinear polarisable materials [1][2][3][4][5] is also proposed: Relation (3) is more general as it gives the incremental work spent upon creating the electric field in matter and in free space, but without any assumption on the relationship between E and D, i.e. without assuming that the medium is linear. From a thermodynamic point of view, relations (2) and (3) represent the Helmholtz free electric energy of the system when assuming an isothermal process under constant deformation [1][2][3].
For an isotropic medium [5], the energy belongings to the dipoles [1,2] is also included in (3) through the polarization vector P that satisfies D = ε 0 · E + P . Finally, the electrostatic energy is made up of two contributions, i.e. U E = U f + U P , where U f is the electric field energy, valid in matter as in free space [2]: whereas U P is an energy related to polarization processes experienced by the body Ω only: In dielectrics, U P can be thought as a transformation of electric energy in some internal energy that belongs to the body. Obviously, this term cancels in ideal conductors since no electric field penetrates inside. However, concerning semiconductors, we can wonder if U P is still the only contribution to the Helmholtz free energy of electric nature.
In this work, we propose to analyze how the electric energy is transferred to semiconductors, and if this still follows the same law as for dielectrics. To the best of our knowledge, a detailed transfer of electric energy in semiconductors has never been examined so far.

Virtual experiment with semiconductor based capacitors
This section analyzes the work spent upon moving a semiconductor based capacitor plate with respect to a counter conductor plate, and compares it with the variation of the Helmholtz free energy as predicted from relations (1) to (3). To this purpose, we assume two semi-infinite capacitors, C0 and C1, with their plates connected such as in Figure 1. Except for C1 where one electrode is a semiconductor (p-type doped, without loss of generality), others electrodes are ideal conductors. This special arrangement makes the thermodynamic analysis self consistent as there is no need for introducing any external voltage source. The case of capacitors biased by means of an external voltage source has been inspected in many details by Bobbio [1,2].

Regular derivation of the total electrostatic energy
Applying the definition of the electric energy given by relation (3) and assuming that electric charges on counter electrodes must compensate each other (we consider semiinfinite plates), we obtain the electric energy for each capacitor system C0 and C1 (scalars consistent with the axis orientation are used instead of vectors): Here Ω 0 and Ω 1 represent the free space between the capacitor plates, g 0 and g 1 are the electrodes gaps for C0 and C1, E g0 and E g1 are the uniform (scalar) electric fields between the electrodes (according to Fig. 1 E g0 , E g1 0) and "t" is the semiconductor thickness (note that t is negative given the origin of the axis). In addition, since we suppose that the polarisation of the semiconductor is linear, i.e. D = ε SC E (ε SC is the semiconductor dielectric constant), relation (7) becomes (according to Fig. 2, E g1 is positive while ψ is negative): where ψ 0 and ψ t are the potentials evaluated at the boundaries of the semiconductor body such as shown in Figure 2, which implicitly assumes the Fermi potential as reference.
We will come to this point later. Note that the last equality in relation (8) holds if the electric field takes a unique value for a given potential. However, this could be unfulfilled when the potential exhibits extrema across the semiconductor layer. In that case, the integral should be split in different domains wherein a unique correspondence exists between E and ψ. Then, the total electrostatic energy U E of the capacitors system is readily obtained as a function of the electric fields and potentials in the semiconductor: Before proceeding further, a discussion to clarify the meaning of the potential ψ is needed. When expressing the electric field as the gradient of a potential, the origin for that potential is a priori not relevant. However, when translating the integral over space into an integral over the potential as in relation (8), the potential ψ should be defined in such a way that it remains consistent with the expression of the electric field. A typical potential profile is shown in Figure 2 where the origin for potentials is taken at the constant Fermi potential (here ψ 0 < 0 and ψ t > 0). But this choice could be different and only after the link between the electric field and the potential has been established, the origin for the potential peculiar to relations (8) and (9) can be set.
After presenting a quite general treatment, the special case of a uniformly doped semiconductor in which a neutral region is recovered when moving away from the free surface will be analyzed in details.

Electric work upon electrode displacement: need for a new energy term
In this virtual experiment, before being isolated, tied electrodes are pre-charged with a total average charge density (per unit surface) Q T = Q 0 +Q 1 (see Fig. 1). Any displacement dx of the semiconductor counter electrode (others are supposed fixed) will induce a variation of the electrostatic energy for C0 and C1; as well as a mechanical work δW F arising from the attractive electrostatic force F that exists between the plates.
Invoking the fundamental law of thermodynamics, when the displacement is performed at constant total charge, i.e. no connection to any voltage source, we can write: As in reference [2], we introduce the Helmholtz free energy for each of the bodies. Relation (10) becomes: where dU f is the variation of the electric field energy in the whole space (including bodies), while dA C and dA SC are the variations of the Helmholtz free energy of the conductor and semiconductor plates induced by the displacement dx while maintaining constant the total charge Q T . Finally, S C and S SC represent the entropy of each subsystem. Note that for dielectrics, the Helmholtz free energy reverts to U P [2]. Still as in reference [2], we assume that the temperature of the system is maintained fixed (T C and T SC are equal and constant). Assuming that there is no deformation, and since the temperature is constant, no variation of the Helmholtz free energy can arise in conducting bodies [2], i.e. dA C = 0. Then, relation (11) simplifies into: Bobbio [2] demonstrated that the term in the right hand side of (12) is the electric energy as defined from relation (3): The link between the mechanical work and energies of electrical nature is now in order. The way Q T will redistribute among the connected electrodes will depend on the total charge density itself and on the electrodes gaps, i.e. g 0 and g 1 ; as well as on the physical nature of the plates, i.e. conductors or semiconductors.
As discussed before, we start considering a very general situation where the semiconductor plate is not necessarily neutral at x = t (in this case, there could also be a charge sheet layer at the contact-semiconductor interface at x = t since the electric field must vanish for x < t). To be consistent with further developments, we will call ψ 0 the surface potential and use the notation ψ s instead.
Given that connected plates share the same potential, we have (E g0 , E g1 0 and ψ S − ψ t 0 (see Fig. 2): Since the sum of the charge densities on each pair of tied electrodes is fixed, we call it Q T , Gauss theorem imposes that the sum of the electric fields E g1 and E g0 is invariant: Next, from the continuity of the displacement vector, without presuming for any fixed charge sheet layer on the semiconductor surface (x = 0), the electric fields across the semiconductor/free space interface satisfy: where E S is the surface electric field evaluated inside the semiconductor at x = 0. Differentiating (9) gives the variation of the electric energy for the capacitors system: Noting that the electric field takes a unique value for a given potential, the last term in (17) simplifies: where E S and E t are the electric fields at x = 0 and x = t, still inside the semiconductor. Next, merging relations (14), (15) and (16) with (17), the change in the electric energy can be expressed in terms of C1 capacitor quantities only: Finally, the force per unit area acting on the semiconductor plate as derived from the electric energy is: This is what we call the 'field energy force' since it is derived from the electric field and dipole energies that are already included in relation (3).
On the other hand, the electric charge density on C1 creates an attractive force between the plates through the Coulomb force F C . This force per unit area is the product of the total charge density in the semiconductor times the electric field generated by the counter electrode, which is half the electric field in the gap since we must exclude the contribution to the field induced by the charge itself (F C is positive according to Fig. 1): Relation (21) defines the "genuine" electrostatic force [5]; as such it should be regarded as the actual force. We expect that the force derived from the free energy (relation (20)) and the Coulomb force should be strictly equal. However, relations (20) and (21) are actually not equivalent.
The unexpected inconsistency between these two formula is a major result per se and suggests that the definition of the electric energy as given by relation (3) does not represent the total electrostatic energy gained by the system.
Without loss of generality, a new contribution is introduced in the Helmholtz free energy. We define this extra quantity by A Extra . Then, the new electric energy U new E for the system of coupled capacitors writes: and relation (13) becomes: Now, if we impose F C and F U to be equal, the extra energy must satisfy the following differential equation: Energy representation of the conductor-vacuumsemiconductor capacitor structure when a neutral region is recovered away from the free surface. The surface potential is the potential drop inside the semiconductor (EF and EFi are respectively the Fermi and intrinsic Fermi energies).
Integrating over g 1 from −∞ to some finite value, and noting that at infinity the electric field E S must vanish since relation (14) cannot diverge, the extra energy is given by: where ψ ∞ S and ψ ∞ t are the potentials evaluated at x = 0 and x = t when the gap separation g 1 goes to infinity. These can be calculated once the physical parameters of the semiconductor system are known. In addition, we implicitly assume that A Extra = 0 when g 1 → ∞. This extra energy involves only physical quantities belonging to the semiconductor (electric fields, potentials and dielectric constant). Then, likewise for the energy of the dipoles U p , A Extra is also part of the semiconductor free energy.
Relation (25) expresses the additional electric energy in a semiconductor layer. This result was not anticipated if we concede that relation (3) was hold to account for the energy of electric nature in matter and in free space. In the next section, a special case where a neutral region is recovered inside the semiconductor layer is illustrated in details.

The case of a partially depleted semiconductor
Without loss of generality, we assume a p-type doped semiconductor layer with a doping density N A . An energy representation is shown in Figure 3 (the semiconductor is set to a lower potential with respect to the counter electrode). In addition, we consider that the semiconductor is neutral at x = t, and therefore E t = 0 (note that for doping densities greater than 10 15 cm −3 , neutrality is recovered after some micrometers only). In this case, it is quite common to define the intrinsic Fermi potential in the neutral body as the origin for the potentials. Therefore we have ψ t = 0 and the extra energy simplifies: Introducing the susceptibility χ = ε SC /ε 0 − 1 and assuming the medium linear, the polarization vector can be written P = ε 0 · χ · E. From relation (5) the dipoles free energy is: Combining (26) with (27) gives: It is instructive to define the energy contribution restricted to the volume of the semiconductor, U SC (note that this is not the internal energy, see [1] for clarification). According to relation (2), we have: In addition to the field and dipoles energies, the new energy U new SC still restricted to the semiconductor plate must now incorporate the extra energy A Extra : Noting that the link between E and ψ is unique and independent of the coordinate, using E s or E in relation (30) makes no difference for the integrals and relation (30) simplifies: As expected, U SC and U new SC are not equivalent. Adding to U new SC the contribution of the remaining electric energy stored in the free space surrounding the semiconductor (Ω −Ω SC ), the total electrostatic energy for the system of capacitors U new E is readily obtained: This electric energy is different from what we used initially, see relation (9). We can gain more insight about the meaning of U SC and U new SC when relying on the standard E-Ψ plot used to model the Metal-Insulator-Semiconductor capacitors system [8] that links the electric field E to the potential ψ (assuming non degenerate semiconductors): where n i is the intrinsic carrier density and U T is the thermodynamic potential ( This representation is very adequate to explain the distinction between the 'typical' and the 'new' formulations of electric energy restricted to the semiconductor body, U SC and U new SC . Basically, the new expression stated by relation (31) represents the area (blue) above the curve in the E-Ψ plot shown in Figure 4a, whereas the contribution given by the standard relation (29) is representing half of the lower area (green). When the surface potential is below the onset of strong inversion as shown in Figure 4a, the upper area is almost half of the bottom one, and it is likely that both definitions converge to the same value. However, when the surface potential is beyond the strong inversion limit as in Figure 4b, the upper area dominates, indicating that relations (29) and (31) are not equivalent.

The case of neutral and charged insulators
There are situations where the electrostatic force derived from the electrostatic energy given from relations (1) to (3) works correctly. We will provide some evidence that this does happen for insulators, and more generally for charged insulators, which are de facto traditional systems coming along with electric energy considerations.
Substituting the semiconductor for an insulator in the system of coupled capacitors sketched in Figure 1, we calculate the Coulomb F C and energy based F U forces and compare them. Regarding the energy stored in the capacitor C1, relation (7) can be rewritten in terms of the electric field in the insulator E i (ε i is the insulator dielectric constant): Introducing the local charge density in the insulator ρ(x) and assuming that ρ(x) does not depend on the local potential, integrating Poisson equation in the insulator gives: Since the charge density depends on the coordinate only, we can express the electric field in the form: where f (x) is a function that is independent of the electric potential.
Similarly, integration of (36) gives the potential distribution in the insulator: After rearrangement, the total electric energy U E = U C1 + U C0 can be written as: Noting that integrals in the bracket of (38) do not depend on the local potential, and thus on the value of the gap g 1 , the force arising from the derivative of the electric energy is: Additionally, connected plates must share the same potential, which from relation (37) gives: Introducing the identity (40) in relation (39), we obtain: Again, charge conservation writes: Similarly, the continuity of the displacement vector at the dielectric-gap boundary gives: Finally, using relations (42) and (43) in (41), we find that the electrostatic force based on the electric energy variation upon the electrode displacement reverts to the Coulomb force between the charged capacitor plates: We can then conclude that as far as insulators are concerned, there is no need to introduce any new free energy of electric nature and the usual definition is leading to the correct result. This apparently 'trivial' finding for insulators could explain why the inconsistency pointed out in semiconductors has been concealed.

A sufficient condition for the existence of the additional free enerxgy
In the quest for a more general criterium, we can analyze when the extra energy given by (26) vanishes, i.e. when the electric energy of the capacitor system can still be obtained from relation (3). Following former analysis, this condition is verified as soon as F U = F C in relations (20) and (21). Imposing this identity links surface electric fields to potentials at front and back interfaces.
Additionally, when the semiconductor layer is neutral at x = t, this condition simplifies further: Since there is a one-to-one correspondence between the local electric field and the local potential, see relation (33) for instance, this condition can also be extended to any coordinate inside the body, leading to the differential equation: Solutions of (47) are readily obtained: where C is a positive valued integration constant. Without loss of generality, we assume ψ S 0 as in Figure 3 so that E(x) = C −ψ(x).
Making use of the Poisson equation, we find that the charge density must be constant (and negative in our case) in the semiconductor where the electric field is not null: (when ψ S 0 the same conclusion applies with a positive charge). Basically, for semiconductors this condition reverts to the so-called full depletion approximation in a uniformly doped material [8]. It imposes that the body is fully depleted down to a given coordinate until it changes for neutrality in a step-like transition. In fact, this analysis makes use of the Poisson equation only, discarding the Fermi-Dirac or Boltzmann statistics, which is a quite common approximation for semiconductors operating in depletion mode.
We can illustrate this result by considering that the semiconductor layer (doped N A ) in Figure 1 gets fully depleted from the surface down to a coordinate x 0 , beyond which it recovers neutrality, meaning that ψ(x 0 ) = 0 and E(x 0 ) = 0. Solving the Poisson equation while assuming a charge density −qN A (−q is the electron charge) between [x 0 ,0] and 0 elsewhere, the potential and electric field at the surface are readily obtained: Therefore the surface electric field and the surface potential are linked as follows: Evaluation of integrals in relation (26) gives: Apart for the sign, these quantities are equal, implying that the extra energy does vanishes in this case. However, as it will be illustrated in the next section, adopting a more consistent approach involving Fermi-Dirac or Boltzmann statistics, which rules out this quite crude full depletion approximation, will highlight a new energy term. The additional energy gained by the semiconductor is then likely to have its origin in statistical physics.

Some tangible quantities: energies and forces in an ideal semiconductor based capacitor
In order to estimate the magnitude and impact of the free energy given by relation (26), we rely on the analytical expression that links the surface potential to the surface electric field in a semiconductor plate. Considering a nondegenerate p-type doped silicon, the surface potential and the surface electric field satisfy the following well-known relationship, valid for depletion, inversion and accumulation [8] (signs are consistent with Fig. 3): where U T is the thermal voltage, n i is the intrinsic carrier density and N A is the doping concentration, other symbols having their usual meaning.

Computation of the extra energy in non-degenerate silicon
As already discussed, adopting the conventional representation E S (ψ S ) and considering only depletion-inversion modes (positive values of the surface potential in this case), relation (26) translates into a simple picture: A Extra is merely half of the lower area (green) minus the upper one (blue) (see Figs. 4a and 4b). This graphical interpretation of A Extra is useful to appraise the meaning of the 'extra energy'. For instance, according to relations (52) and (53), the area in the upper side of the curve (blue) is about half the area below (green) as long as the surface potential Ψ S remains below 0.7 volt (see Fig. 4a), which is the threshold of strong inversion in our case. This is where the potential and the electric field satisfy the full depletion approximation to a good extent. Indeed, while Ψ S remains below 0.7 volt, holes can be neglected while keeping only the linear term in ψ s under the square root in relation (54), leading to: which is formally the same as relation (48b). However, increasing the surface potential beyond that limit (0.7 volt) increases the weight of the blue area with respect to the green one (see Fig. 4b), and predicts that the two definitions of electric energy given by relations (29) and (31) will start to deviate from one another. In fact the amount of the additional free energy A Extra will increase exponentially with the surface potential.
Similarly, for negative values of Ψ S , i.e. accumulation mode, we can show that the additional free energy A Extra is never negligible whatever the value of the surface potential (the surface 'above' the curve is always greater than half of the 'bottom' counterpart). Therefore, we anticipate that as soon as the surface potential matches with inversion or accumulation, A Extra will never be negligible. All the more, it will be the dominant contribution to the total electrostatic free energy.
To perceive how the Helmholtz free electric energy is shared between the two 'components' A Extra and U p , we evaluate the ratio A Extra /U P versus the surface potential for various doping densities (see Fig. 5a). The additional electric energy upon a change in the surface potential dominates when the semiconductor is set in accumulation, i.e. ψ S < 0 V, or above the threshold of inversion (i.e. ψ th S > 0.7 V for N A = 10 16 cm −3 ). In inversion mode, this 'threshold' ψ th S is shifted towards higher values when the doping density is increased, whereas in accumulation the doping density has almost no effect when illustrated with respect to the surface potential. Note that before reaching the strong inversion limit (which depends on the doping density), the ratio A Extra /U P remains lower than unity, meaning that the extra energy is minimized, which is expected since this is when the full depletion approximation is satisfied to a good extent.
A somehow different representation of A Extra /U P in terms of the charge density stored in the semiconductor is shown in Figure 5b. Except in a limited region where  5. Ratio between the extra free energy AExtra and the dipoles energy UP versus the surface potential (a) and versus the charge density (b) in a silicon layer, for different doping concentrations. In (b), we note that in accumulation, AExtra is always dominant, in contrast to depletion-inversion mode where it can have a minor contribution. At flat band, the ratio takes a well defined value independent of the doping density.
-Q SC is lower than about 10 −3 C/cm 2 (depletion), for a given space charge density the energy ratio increases with the doping level. In addition, A Extra exceeds always U P in accumulation, which is not the case in depletioninversion mode where it can be relatively small, but still never cancels. Interestingly, when the surface potential is close to the flat band (ψ S ≈ 0 V), the extra energy tends to an asymptotic value independent of the doping density, A Extra /U P ≈ 3ε SC /ε 0 · χ (about 1.09 for silicon). Then, the additional electrostatic energy stored in a quasi neutral silicon layer still represents about half of the electric field and dipoles energies restricted to that silicon body.

Appraisal of Coulomb and energy based electrostatic forces
The aim of this section is to illustrate the extra energy contribution through the tangible quantity of electrostatic force evaluated in a nanometer scale capacitor systems.
To that purpose, two values for the electrode gap of C1 are chosen, namely 100 and 10 nm, as well as different doping densities. Without loss of generality, for each case we assume identical electrodes separations for g 1 and g 0 , respectively for C1 and C0 capacitors. Regarding the total charge density per unit surface Q T , it has been assigned values consistent with what is found in MOS and MEMS devices under normal operation (note that Q T is a dummy variable used to generate charge densities, related potentials and electric fields in C1 and C0). Solving the set of equations involving electrostatics and semiconductors physics gives the charge densities on each capacitor plate, as well as the surface potential in the semiconductor. Next, the ratio of the Coulomb to the energy-based forces F C /F U (see relations (20) and (21)) are displayed in Figures 6a and 6b for the two values of the capacitors gap.
Considering the case where the electrodes gap is 100 nm, Figure 6a confirms that F U does not equal F C , which was indeed the argument revealing some inconsistency in the "regular" expression of electrostatic energy. This mismatch gets smaller for substrate doping densities greater than 10 16 cm −3 . However, it is worth noticing that the link with the extra energy is not straightforward when comparing Figures 5a and 5b with Figures 6a and 6b since a large amount of extra energy does not necessarily mean a large discrepancy in the Coulomb and energy derived forces. In the capacitor system, the dependence of the force with the doping and charge density is less evident as it depends on how this energy is affected by the electrodes separation, i.e. the derivative of the energy with the gap g 1 . In addition, part of the total electric energy is also stored between the plates.
Decreasing the gap down to 10 nm is even more significant in terms of forces deviation. Figure 6b reveals that F C and F U can differ by nearly one order of magnitude for a doping density of 10 15 cm −3 , and by a factor close to 4 for highly doped substrates. The reason why F C /F U is enhanced in smaller gap capacitors can be understood as follows: as far as the energy stored in free space between the plates exceeds the energy stored inside the semiconductor, A Extra will have a limited influence and F U will not depart so much from F C . This happens for relatively 'large' systems with a gap laying in the micrometer range. But in nanometer scale capacitors, the electric energy stored inside the semiconductor becomes comparable with the electric field energy stored between the plates (in free space). In this case, energy-based and Coulomb forces may deviate significantly from one another. As rule of thumb, we can say that as far as the gap between the plates and the extension of non-neutral regions in the semiconductor are in the same range, A Extra might not be negligible. This will happen for deep sub-micrometer semiconductor based capacitor systems.  . 6. Ratio between the electric force calculated using the Coulomb force (FC ) over the force obtained from the standard definition of electric energy (FU ), as a function of the charge density, and for different doping concentrations. Two values of electrode separation are addressed, namely 100 nm (a) and 10 nm (b). These calculations reveal that the usual definition of electrostatic energy is not giving the correct value. These discrepancies are even more obvious when decreasing the gap from 100 to 10 nm since the impact of the energy stored in the semiconductor layer becomes enhanced in regard to the energy stored in the free space between electrodes.

Generalisation to arbitrary geometries
The simple picture of an ideal semi-infinite capacitor system revealed unambiguously the need for an additional electric energy in matter, which has been straightforwardly illustrated for ideal semiconductors layers. Here, we propose to generalize this analysis by considering a three dimensional semiconductor body which is separated from an ideal conductor 'C' as shown in Figure 7. Likewise in the former analysis, we connect this system to an ideal capacitor C0, but here the idea is to transfer some electric energy from C0 to the Electrode-Semiconductor system Fig. 7. Representation of the system used to derive the electric energy in a semiconductor body of arbitrary shape. The virtual capacitor is used to generate a continuously varying voltage by means of the work done by an external force on one electrode. The energy transfer between the ideal capacitor and the semiconductor-outer electrode system still reveals that an energy term is missing in the standard formulation of the electrostatic energy.
by increasing gradually the gap between the pre-charged capacitor plates of C0 until reaching a preset value. Basically, such a process will raise continuously the potential V P at node P in the semiconductor while engendering a transfer of charges between the two systems.
As in the first experiment, this method will involve a mechanical work generated upon moving the plates against the Coulomb force that should still satisfy the principle of energy conservation for the whole system. The advantage for using this setup to generate a potential is that it relies on similar thermodynamics arguments that have been introduced in Section 2 given that both experiments share equal concepts.

Analysis of the energy transfer process
The electric energy stored in the capacitor C0 can be calculated from relation (1) since we assume ideal conductors for the variable capacitor plates (no penetration of the electric field in conductors): where Q is the total charge (assumed positive without loss of generality) on the fixed electrode of C0, which should not be mistaken with the charge density per unit area, and V P and V C are the inner and outer electrode potentials (see Fig. 7). Since the three dimensional semiconductor has a limited volume, we must also assign a finite value to C0 otherwise any finite transfer of charges from C0 to the semiconductor will affect neither the charge density, nor the potential drop across C0. The electric field between the plates and the potential difference are given by (the sign of E is consistent with the orientation of the x axis): where S is the area of the capacitor plates which are separated by x. Using (56), the electric field energy stored between the capacitor electrodes is: As in Section 2, we assume rigid bodies maintained at a fixed temperature and we invoke the same set of hypothesis. The mechanical work spent by applying a force F on one electrode to increase the gap between the plates by dx is related to the change in the electric energy of the whole system dU E (note that the work increases the energy of the system when increasing the gap between the electrodes since F · dx is then positive): Here U Ω holds for the electric energy stored in the semiconductor body and in free space, excluding the capacitor C0 which is already accounted in U C0 . From relation (59), we obtain: On the other hand, an attractive Coulomb force exists between the plates. Its module is given by: Under quasi-static conditions, the applied force F and the Coulomb force F C should compensate each other. Therefore the work spent upon moving the plate by dx is: Then, using (61) with (63), we get: where Q SC represents the charge in the semiconductor body and, since the total charge on connected electrodes is invariant (i.e. δ(Q SC + Q) = 0), dQ SC = −dQ.
The way dQ and dx are interrelated depends on the interaction between the semiconductor body and the conductor C. Here, what matters is that V P −V C can take any value (provided Q is not null) by varying the distance between the electrodes. In particular, when electrodes are "in contact" (but still "isolated" electrically), the potential drops to zero. Therefore, using the same kind of topology as in Section 2, we create a virtual variable voltage source while satisfying the principle of energy conservation.
Next, relation (64) can be formally written as a function of the charge transferred to the semiconductor body: where V P C = V P − V C depends implicitly on Q SC .

Derivation of the total electrostatic energy in the semiconductor body
According to Gauss theorem, a variation in the semiconductor charge reverts to a variation in the flux of the surface displacement vector dD through the surface S SC , where S SC consists of the combination of the outer surface and inner semiconductor-electrode interface (see Fig. 7): Here dS is the unitary surface vector oriented outward for the outer surface, and inward for the inner semiconductorcontact interface (note that a virtual "wire" should be created inside the semiconductor body to access the inner contact. However, this can be made infinitesimally narrow so that its shape will not contribute to the surface and volume integrals that come into play with the divergence theorem). Since the potential V P C is constant, it can be moved inside the integral: (67) Figure 7 shows a line L joining the point P in the semiconductor to the external conductor C. This line crosses the surface S of the semiconductor at N . Since the potential drop between P and C does not depend on the path (the voltage drop is zero in a closed loop), we can write: In relation (68), the quantity (V N − V P ) represents the potential drop in the semiconductor between the inner contact and the point N at the surface (which reverts to ψ S − ψ t in the case of the semiconductor layer discussed in Sect. 2): The last integral on the right hand side of (69) can be transformed into a field volume integral by using the divergence theorem. To do so, we consider the volume in the whole space (without C0), excluding the semiconductor, i.e. Ω − Ω SC . This consists in the volume of the semiconductor Ω SC , the volume of the electrode Ω C and the volume of the connecting "tube" between them Ω Connect introduced to ensure the continuity for the inner surface of Ω − Ω SC (see Fig. 8).
In order to simplify relation (69), we start using the divergence theorem: leading to: Since there is no charge in Ω − Ω SC , the second integral in (71) vanishes, leading to: Noting that Ω − Ω SC on one side, and {Ω SC , Ω C } on the other side have opposite surface orientations, using the divergence theorem again and noting that at infinity the displacement vector must vanish, we have: Next, the connecting element is made infinitesimally small, so that the last integral over S Connect in (73) can be omitted: Noting that V C is constant on the conductor surface and that the integral of dD over the conductor surface S C is simply −dQ SC , as imposed by the charge neutrality principle, we can write: (75) Using the above identity in relation (69), the incremental energy transfer from the capacitor to the rest of the system, including free space (except C0) is: (76) Next, the total work exchanged is obtained by integrating (76) over the displacement vector while noting that there is no charge in Ω − Ω SC : Permuting the integration sequence, we obtain: The electric energy supplied by the virtual voltage source is divided in some energy stored in free space and in the semiconductor body. On the other hand, the electrostatic energy in the free space Ω (excluding C0, but including the semiconductor ) is still given by relation (3).
Likewise for the analysis of the capacitor plate, the difference in these quantities represents the additional free energy for the body: Therefore, the new definition of the electric energy of a system including an arbitrary semiconductor body is proposed: It is straightforward to see that relation (80) is the generalization of relation (32) that was developed for a planar geometry. Therefore, even though we followed a quite different approach, we conclude that an additional free energy created by an electric field is predicted in semiconductors, and very likely in matter in general, and only in some special cases this energy is not engendered.
We can proceed further these developments and propose an expression in terms of the local charge density. Using the divergence theorem for the surface integral of (80), we obtain: Introducing this identity in relation (79), we find that the energy contribution coming from the semiconductor body takes a very simple form: At least in semiconductors, A Extra is always positive, i.e. the internal energy is increased upon generation of nonneutral regions. In addition, using relation (1) (valid only if the system is linear [5]) for the standard expression of electric energy, we can write the new energy as a function of the potential and charge density only: Two clarifications need to be done. Firstly, the volume Ω in (80) is used to identify the system of interest (semiconductor body with outer electrode) excluding the 'capacitor voltage source'. But this still represents the whole space termed as Ω ∞ in relation (1), reason why it is used in place of Ω in relation (83). Next, the meaning of the potentials φ and ψ should be well understood. Whereas φ(x) is the potential created at x by all the charges in the sense given by relation (1) and should therefore vanish at infinity, ψ(x) is the difference in the potentials between the inner electrode P and a coordinate x inside the semiconductor body (Fig. 7). Finally, when considering isotropic media and introducing the polarization vector as in reference [2], the different contributions to the electrostatic energy are highlighted: The first term represents the total energy stored in the electric field in the whole space, including the semiconductor body. The second term represents the energy supplied to polarize matter [1,2]. Lastly, the third term represents the new energy A Extra . Different expressions for the electrostatic energy with the new definition are then obtained. Note that whereas relations (80) and (84) generalize relation (3), relation (83) is somewhat more restrictive as it is makes use of relation (1) and hence should be used with the same restriction, i.e. for linear polarizable systems (in the sense discussed in Ref. [5]).

The full depletion approximation in a semiconductor sphere as a case of study
We illustrate the use of generalized expressions of the extra energy in the case of a p-type doped semiconductor sphere of radius R where the full depletion approximation is invoked. We consider a kind of semiconductorconductor capacitor system as shown in Figure 9. Likewise for the planar geometry discussed formerly, applying a lower potential to the silicon with respect to the surrounding conductive sphere will deplete the surface of the Fig. 9. Cross section of the spherical semiconductor body surrounded by an ideal counter electrode that creates a depletion in the inner sphere (electrical potentials are applied at the center and on the outer conductive shell). The radius of the semiconductor sphere is R and the radius where the transition between full depletion-neutral semiconductor takes place is labeled r0.
semiconductor sphere from mobile carriers. In addition, as the semiconductor is uniformly doped, there is no angular dependence of any physical quantity, implying that the Poisson equation in spherical coordinate satisfies: The full depletion approximation supposes that the core of the semiconductor sphere remains neutral up to a radius r 0 , then is fully depleted from r 0 to R(ψ(r 0 ) = 0, E(r 0 ) = 0). Imposing these conditions to relation (85), the potential and the radial electric field depend on the radius r and on the parameter r 0 such as: We are now able to evaluate the extra energy defined from relation (79) for instance. Concerning the surface integral, since we are considering a sphere with no angular dependence, the integral over the surface for a given surface displacement vector D S simplifies as: Concerning the surface potential, this is straightforwardly obtained from relation (86): Next, the incremental surface displacement vector is recovered from relation (87): (note that there is no dependence on the radius since D S is defined at the surface only). After calculations, the surface integral becomes: Next, the second contribution is evaluated. Still assuming that there is no angular dependence and that the medium is linear, we have (dΩ = 4 × πr 2 dr): It is worth noticing that in (91), integration is performed over the parameter r 0 , which is the radius where depletion starts, whereas in (92) the integration is done over the variable r from r 0 to R. Comparing both terms reveals that these integrals are equal and will cancel each others in relation (79). Therefore, the regular definition of the electrostatic energy remains valid when the full depletion approximation is used even in case of the semiconductor spheres. So far we have carried out the analysis making of use of fields, which was instructive to understand how these relationships can be used in practical cases. In the next section, we will adopt the dual charge-potential approach.

The full depletion approximation in arbitrary geometries
As stated in relation (82), an alternative expression of A Extra in terms of charges and potentials is possible. Here, we can gain more insight when integrating by parts the