Open access peer-reviewed chapter

Local Electric Fields in Dielectric and Semiconductors

By Dmitry E. Milovzorov

Submitted: June 6th 2017Reviewed: January 23rd 2018Published: May 23rd 2018

DOI: 10.5772/intechopen.74310

Downloaded: 476


Local electric fields are appeared in dielectric and semiconductors due to the destruction of symmetry, creating the vacancies, point defects and chemical impurities in material. By increasing in external electric field value there are numerous structural changes will be generated. Point defects in silicon films were characterized by using electron-paramagnetic resonance spectroscopy and laser picoseconds spectroscopy. Chemical bonding properties was investigated by means of Fourier-transformed infrared spectroscopy. The possible mechanism of phase destruction was proposed.


  • local field
  • point defect
  • dangling bonds
  • Raman spectroscopy
  • nanocrystals
  • second-harmonic generation
  • silicon films

1. Introduction

Local electric fields are appeared in dielectric and semiconductors due to the destruction of symmetry, creating the vacancies, point defects and chemical impurities in material. By increasing in external electric field value there are numerous structural changes will be generated. Some of them will produce such great local fields that will destroy all material or change its physical properties. The studying the nature of local electric fields will open new tendency in electronic device producing, from one side, and, help to change materials’ properties according to our needs, from another side.

Description of local electromagnetic fields is a continuously durable through the all history of physics and was began with publication of first articles written by Maxwell Garnett which were devoted to colors in metal glasses and metallic films [1], Lorentz [2], and later in works of Brugeman [3] was developed by Edmund Stoner from University of Leeds [4] and Osborn from Naval Research Laboratory [5]. For a complicated medium such as the binary system with components A and B the dielectric function can be estimated as following [5]:


In a case that one component is included in another dielectric component and polarized media with averaged value of polarizability of αw=εAwεBwεAw+2εBwthe Maxwell-Garnet formula is given by εw=εBw1+2CAεAwεBwεAw+2εBw1CAεAwεBwεAw+2εBw1.

Structural properties of material may be strictly different as for surface and thin films, as for nanostructures such as clusters or nanocrystals, as for bulk material. However, it is obviously that most amounts of media in the universe is nanostructured or even in nanocrystal phase. For example, interplanetary dust was observed charged coupled devices (CCD) detectors and infrared space telescope [6]. They observed a cometary coma of Hale-Bopp comet. The dust destruction in space (Cygnus Loop supernova) was observed by using a Spitzer Telescope tuned in infrared (IR) range from 22 to 36 μm [7] and shows us the chemical properties, such as dust chemical compositions. These observations show the great fraction of silicon in all space dust. Space dust destruction and ion formation was studied by Mann and Czechovsky [8], which results from model calculations in silicate grains, carbon and ice grains. Grain destruction in a supernova remnant shock wave was investigated by astronomers of Harvard University [9]. It was observed by Spitzer telescope in IR 24 μm range of wave length. The case of impact of nanoscopic dust grain with solar wind of spacecraft already was estimated by using the dimensionless parameter equals to ratio between Debye length and radius of dust cloud spherical shell with radius R [10]. It was shown that the dust particle with mass 10−20 kg produces by impact 107 charged particles.

From the other physical scientific trends we have an observations of local field by a nonlinear spectroscopic experiments with nonmaterial and nanocomposites, particularly, semiconductors. By using semiconductor materials have been made numerous types of devices, such as electronic devices and photon detectors, integrated circuits and thin film transistors, optoelectronic devices and others. Every time when the device is developed the problem of reproducibility of its work regimes and durability of their realization is appeared. The solution of this problem is very important for device manufacturing, and it depends on properties of used active semiconductor materials. The electrical properties often are not so transparent due to slightly nonlinear behavior of their characteristics. Figure 1 shows the current-voltage and resistance-voltage characteristics for two silicon films prepared by plasma-enhanced chemical vapor deposition technique with gas mixture of silane diluted by hydrogen and silicon tetra fluoride gas: amorphous and nanocrystallized [11]. It is seen, than their current-voltage characteristics are similar in this voltage range, but resistance-voltage characteristics are strictly different. Such difference can be explained by the disorder of amorphous phase and generating of numerous point defects by applying external electrical field. It is clear, that the voltage is varied in the range from the −10 to 10 V. Hydrogenated amorphous silicon was widely used in last decades in electronics. In recent years the nanocrystalline silicon are studying for many technological applications. The structural transformation from crystal to disordered materials, however, is investigated very poor, mainly resulted in Staebler-Wronski photo-stimulated effect. However, the electric field applied to the nanostructured silicon thin film gives the new possibility to change structural order. Such kind of structural transformation is caused because of there are numerous defects inside the silicon film.

Figure 1.

The current-voltage and ohm-voltage characteristics for two silicon films: amorphous and nanocrystallized.

The anomalous characteristic of resistance-voltage can be explained by random distributed the point defects inside the amorphous film along with the hydrogen atoms, and existing the dipoles Si-O which turn to compensate the external electric field. But, for the nanocrystalline silicon film, the point defects are incorporated into silicon nanostructured net and cannot move freely, because there is a stabile electric characteristic for nanocrystalline silicon film, and anomalous for amorphous.

The other new area of scientific interests is crystal-amorphous phase transformations by applying electric fields and role of local fields in phase transition from order to disorder. Because, it is important to investigate the point defects which can be responsible for local electrical fields generation in polarized media, such as dielectric silicon oxide media or semiconductor thin film of silicon. The main role plays here the silicon-oxide bonding in side thin film of silicon. Si-O dipoles play a dramatic role in crystal phase destruction by applying electric fields. The induced dipole moment by applied electric field can be written in the following form [12]:

d=α0E0j=1Nδαβ3rijαrijβrij2drij3;where E0is the field at the i-th monomer, rijis a relative vector between i-th and j-th monomer molecules, α0 is the dipole polarizability. However, the Hamiltonian of semiconductor cluster can be surely expressed by using donor and acceptor states in bulk material [13]: H=H0+HD+HA, where H0 is the Hamiltonian pure semiconductor for electron matrix elements of transitions between own conductive and valence bands, but HD and HA are the exchange energies matrix elements due to the donor and acceptor states, and for them it is clear be presented the following expressions:

HDij=kDki2Eεkσkδij; HAij=lAli2Eεlσlδij; where Dkiand Aliare coupling between k-th donor state and l-th acceptor state, but the value σk,l=12θk,li4γk,l2θk,l2are the self-energy corrections which depend on locality of sites, and θk,l=Eεk,l, and γk,lare the ionization energy for electrons in sites, and overlap integral between atomic orbital for the donor in k-site and acceptor in l-site, respectively. Therefore, different explanations of polarization properties of semiconductor media can describe only partially in their own borders the electric induced local fields’ appearance and various theoretical predictions based on them are not satisfied. The macroscopic description and molecular nanoscopic model are poor for investigating the mechanism of crystal semiconductor structures destruction by applied electric fields because the their dielectric functions is not so transparent for these complicated media, for example SiOx, from one side, and their ratio between covalent and ion fractions of inter atomic bond are not so homogeneous, from the other side. Because, there is a necessity of detail investigation of nanoscopic nature of local fields appearance and mechanisms of crystal-amorphous phases’ transformations.

The present work is devoted to the nature of local field appearance in silicon nanoscopic material and role point defects in phase transformation of material.

2. Experimental researches of field-assisted destruction of silicon nanocrystals

Nonlinear polarization associated with the phonons can be written as P=1Vi=1N<μ>i=NVαQ0<Q>E. Using these equations it is possible to obtain the wave equation for field: ΔE+η2c22Et2=μ02Pt2, where c2=1μ0ε0,η=εε0.

Raman effect is result from the interaction of an electromagnetic field and optical phonon mode. The vibration wave <Q>=Q2exp(iωvtkvz+c.c.The polarization at the Stokes frequency PS=ε0χRωvEi2ES. Raman susceptibility


The microcrystal wave function is a superposition of Eigen functions with kvectors centered at k0. We suppose that k0=0. If we are using the weighting function as Gaussian we will have


the normalized first-order correlation function g1τ=<EtEt+τ><EtEt>. First-order coherence g1τ=expkτΓ02τ,where ωkis the phonon dispersion curve, Γ0is the natural line width. The spectrum of natural broadened value is Fω1ωωk2+Γ022. The phonon probability ρ=Ψ0Ψ0=u02rd3kC0k2. The first-order Raman spectrum


For a microcrystalline and nanocrystalline silicon with sizes of crystals L if the weight function is Gaussian the first-order Raman spectrum is following [14]: Iω=01expk^2L2/4a2ωωk^2+Γ0/22d3k^, k=2πak^, ais lattice constant, k^is dimensionless value,Γ0˜3.6cm1line width of the Si LO phonons in c-Si. The dispersion of the LO phonon in c-Si ω2k^=A+Bcosπk^/2, where A = 1.714 × 105 cm−2 and B = 105 cm−2 [15]:


where local field factor can be written as


where q is a vector of inverted lattice, 2π/a, a is a lattice constant. If the value λ0 is the bond length Si‐Si in bulk silicon, λ1 is the weak bond resulted the silicon-oxygen surrounding. By the value of density of bonds is NSiO = 21021 cm3 and for crystalline silicon the density of silicon bonds NSi−Si = 51022 cm−3 the lattice constants are following: a2 = 0.98 a0 аnd a1 = 0.996 a0.

For the silicon nanocrystalline and microcrystalline films the phonons can be generated in crystals by laser field or annealing. The wave of deformations can be generated by picoseconds laser pulse [16]. The acousto-electric effect was observed in n-type germanium [17]. The electric field which was appeared by ultrasound waves can be estimated by using the formula F=eE=q2Sc2KTω2τ1+ωτ2,where E is an acousto-electric field, τ is a relaxation time τ1=43τэ+k2D, D is a diffusion coefficient, The deformation potential causes the appearance of effective acoustical charge (Figure 2).

Figure 2.

Raman spectra from nc-Si film with nanocrystal mean size 24 nm.

The nanocrystalline film was made by me using CVD method of silane diluted by hydrogen (gaze flow rates ratio is 1:10) at low temperature of substrate (80°C). The RF power was 20 W. Working pressure was 0.2 Torr. The crystalline volume fraction was 66%. The crystal orientations for nanocrystals were determined by means of X-ray diffraction technique (111) and their average size was 24 nm. The thickness of silicon film was more than 300 nm. Figure 5 shows the changes in Raman scattering spectral data by applying the external electric field. It is seen, that there is phase destruction by the relatively high voltage. It is assumed that the nanocrystals which have grain boundary with oxygen atoms incorporated into silicon were destroyed in their crystal structure by Si-O dipoles reorientations caused by applied field. The initial crystal orientation was (111). The incorporated oxygen atoms are adsorbed in determined places. Their position results the appearance of numerous dangling bonds which are multiplied by the electric field and create the deep cracks in crystals. The crystal order is damaged along the axis that is perpendicular to (111). According to the Raman data for SiO2 [18] the Raman spectrum of SiO2 has the variation modes D1 (at 490 cm−1) with defects and activation energies 0.14 eV and pure mode w1 of Si‐O‐Si bridge. The sum dipole moment consists of dipole moment that is created due to the ellipsoidal shape and because of surface charges are appeared by silicon net deformation due to the oxygen incorporation in silicon SiO or hydrogen termination of dangling bonds and creation the SiH bonds: PSiM=83πab2NSiMδSiMba, where NSiM, δSiM are the density of the bonds and partial charge for SiM bond (O or H), respectively (Table 1).

NDB,cm−3P ellipsoid, DebyeΔE = μEext, μeV
by the Eext = 106 V/m
ΔE = μEext, μeV
by the Eext = 107 V/m

Table 1.

The polarization of elliptic silicon grains evaluated by fixed NDB and stark energy shifts for levels in electronic structure by external electromagnetic field.

It is necessary to note that the fractal structure of several kinds of nanocrystals may cause the dramatically changes (four orders of magnitude) in intensity of Raman scattering due to existing of plasmon resonance into the gaps between the fractals [19]. The Raman intensity by these conditions can be expressed as


where α is a polarizability and σ is a local conductivity of a fractal structure.

In addition, the light irradiation of amorphous silicon film causes the point defects generating and, mainly for amorphous hydrogenated silicon films, causes the appearance new dihydride configurations: (H‐Si Si‐H)2(H‐Si Si‐H) and SiH2 [20]. The two atoms of hydrogen in the SiH2 unit show an average proton separation of 2.39 Å. Because, for the hydrogenated silicon nanocrystalline films under influence of applying the external electric field the hydrogen diffusion increases and polysilane chains are created, surely.

For poly-Si films with nanocrystals the values of densities of SiO and SiH bonds varies in wide range from 1019 to 1021 cm−3. The density of bonds were estimated for the poly-Si films prepared by using PECVD as following: for Si-Si bonding the density of bonds is equals to 51022 cm−3, but densities of SiO and SiH bonds are 1021 cm−3. In these films there is an oxygen contamination on the 2% level. The S/V ratio is 1.25%. I suppose that all the oxygen is concentrated around crystals in their grains boundaries. By these values of densities the dipole moments causes by surface charges can be estimated as PSiO = 0.04 D and PSiH = 0.015 D. By applying external voltage it is clear that all Si-O dipoles move to compensate external field and destroy the crystal structure. Because, by applied electric field I observed the crystal phase destruction and SiOx creation by using the Raman scattering data which correspond to the results reported in work [21].

3. Matrix Hamiltonian by small perturbation of Si‐Si‐Si bridge

For very small nanocrystals with sufficient ration S/V the mechanism of three elements simultaneous interaction is important for precise calculations. The energy shift due to the stress appearance for crystal orientation (111) is less than 0.14 eV for the vacancy-oxygen (VO) complex by a stress 0.3 GPa. The Hamiltonian of such system of n atoms as for example,‐Si‐Si‐Si‐ and, particularly the Hamiltonian of interaction between atoms with indexes k−1 and k + 1 can be explained in matrix form is given by using the operators of creation and elimination of boson particles, such as phonons: H'=αck1+ck+1+e.c.The Hamiltonian matrix of interaction between three atoms (with upper indexes k−1, k and k + 1) in chain by a small perturbation α due to the VO appearance (for bonded Si1 and Si3 atoms without oxygen) can be written as Ĥ=A11k1A12αA21A22kA23αA32A33k+1. This matrix can be transformed into next triangular form:


Here, the matrix elements Â22=A22A21A111A12;Â23=A23αA111A21;Â33'=A33αA111αÂ32Â221Â23;Â32=A32αA111A12;and according to the Vilandt-Hoffman theorem of matrix perturbation theory [22] the following inequalities can be written.

i=1nλAiλAi21/2A+EAEor λ̂332λ332+α2λ11. Because, for the creation operator ĉ2k+1tis true the following expression ĉ2k+1t=c2k+1texpi0tα2E1. Therefore, the dispersion curve is written as ωphonon = ω0 ± ξ, where ω0 dispersion curve without applied field (for energy 0.14 eV the estimated additional frequency is less than 3.4 × 1013 Hz, and additional frequency ξ= α2λ11by stress due to the appearance of VO defects. According to the proposed model that is suitable for description Raman scattering phenomena caused by nano-sized cavities it is assumed that the results Raman frequency of radiation after scattering can be written as ωS=μE;where μ=dQ is dipole moment of nanocrystal, E is field, d is the size of crystal, Q is charge (Figure 3).

Figure 3.

Scheme of Raman scattering of photons on phonons (a) and Raman spectra for silicon film without external field (b) and by applied external electric field with additional spectral peak at 540 cm−1 which was created due to the defects generating by electric field (c).

4. Model of polaron state in silicon nanocrystals

The probability of changing the polarization state from one to another can be described by using Golden rule of Fermi


where ωpolar=PEis a polarization energy that can be expressed in frequency units. Matrix element for such changing of polarization state from i to state j, according to theoretical work of V. Lakhno [23], can be written as


where V is a nanocrystal volume (Figure 4).

Figure 4.

Silicon nanocrystals with different sizes by applying external electric field have various electronic structure: with polaron state in all volume of nanocrystal, and partially polarized nanocrystal according to relation between the values of size and “skin”-layer or depth of field penetration.

The second-harmonic generation is forbidden for center symmetric crystal such as bulk silicon because the sum dipole moment is zero, but is possible due to the surface breaking symmetry and quadruple terms contributions. The opposite situation is for nanostructured oxidized silicon film, the surface area for a great amount of nanocrystals is significant, the breaking symmetry is permanent and lateral isotropic. The oxygen atoms with concentration up to the values of 1020–1021 cm−3 show the sharp increase in SHG by increase in polarization properties of material, that have its properties as silicon nanocrystals, as silicon oxide inclusions. Figures 5 and 6 illustrate the SHG spectra of radiation reflected from silicon films.

Figure 5.

SHG spectra for χ(2)xxx component of susceptibility.

Figure 6.

SHG spectra for Pin➔Pout laser scheme and χ(2)xxx component of susceptibility.

The reflected SHG response was measured by using the radiation of optical parametric oscillator/amplifier pumped by the third harmonic (355 nm) of a Q-switched Nd: YAG laser (Spectra-Physics, MOPO 730) at a 10 Hz repetition rate with spectral range between 440 and 1700 nm. The bandwidth of radiation is 0.3 cm−1. The SHG response was detected by a photomultiplier tube and gated electronics with an average of 100 pulses. The linear polarized radiation was focused on the surface of the sample at the angle 45° and detected SHG signal was observed at the angle 45°, too. Such optical scheme arrangement was useful for surface contributors’ detection from the silicon surface (111). The diameter of irradiated spot was 0.5 mm. The energy of the primary laser beam was 4 mJ. The second-harmonic intensity can be written as


where the L(ω, ρ) value is a local field factor of film with crystalline volume fraction equals to ρ = 70%.

The SHG intensity as a function of the average grain size in poly-Si films, with crystalline volume fraction 70%, is presented in Figures 5 and 6 where


is normalized SHG signal, L2ωρand Lωρare the factors of the local field, where εc and εa are dielectric functions of crystalline and amorphous silicon, respectively. For sphere depolarization factor Λ is equal to 1/3, β is the Lorentz constant (for homogeneous spherical surrounding β =1/3). For calculations the following estimated values were λw = 1064 nm, λ2w = 532 nm, εc’(w) = 13, εc”(w) = 0.03, εa’(w) = 12, εa”(w) = 0.23, εc’(2w) = 18, εc”(2w) = 0.5, εa’(2w) = 18, εa”(2w) = 7.5. Resonance spectra consist of two sharp peaks. The peak at 1.6 eV is caused by SHG response due to E0’= Γ21−Γ15 transition. We suppose that the second peak can be recognized as SHG response due to E1 = L2’−L1 transition in silicon nanocrystallites.


The contribution of point defects as deviations of local fields and external applied electric field for a phonon generation in silicon nanocrystalline can be described by using the perturbation theory. The model Hamiltonian’s matrix for two-level system including the point defect as small perturbation ε that causes the violence of phonon energies of system E1 (for unperturbed state) and E2 (perturbed state):


The changes in Eigen values from the E1 and E2 by the field E = 0 to the new values of energies are following


It is assumed, that the field value ε determined as linear combinations of external applied field and all local deviations due to film structural disorder. We assume that there is no strict disorder media, but some small disorder is determined. εE1E2; ΔE2E1E2=2Δ. Therefore, by a small perturbation of system εE1E2, the energetic gap between the two levels which are located closed each to other increases from the zero value to 2Δ. The value of perturbation of atomic orbital for dimer molecular like assembly Si‐Si in a point defect as (VO) can be evaluated by using the energy of their interaction of atoms with dipole Si‐O, that has its polarized charge 0.2 е: F=14πε00.2e2r2; where r is a distance from the dimer Si‐Si to dipole Si‐O equals to 1–2 Ǻ.

From the other side, for drift of particle, such as hydrogen atom, by driving forces in condensed matter can be expressed by using formula for force F=HMkTkT, where HM is an enthalpy or energy of transport by heat, Tis a gradient of temperature [24]. For elastic medium the force that treats the defect is given by using stress value σ and displacement field of defect in its surrounding medium as u: F=Σuσσuds. David Emin studied the deformable lattice in 1972 [25] by using short-range electron-lattice interaction in dielectric or semiconductor materials which have weakly coupled electron-lattice interactions or small polaron states, the distortion related energy of which expressed in terms of distortion parameter x, E=12Mω2x2, the energy reduction due to the linear electron-lattice interaction as E=Ax.

5. Modular group translation model for crystal phase destruction by applied electric field

I propose the modular group translation (MGT) model for crystal phase destruction by applied electric field for explanation the Raman data which are on Figure 7 and show the dramatic changes in silicon crystal phase related spectral component at 520 cm−1 due to the applying electric field.

Figure 7.

Scheme of dipole’s and electric field vectors location on 2D plate.

It is assumed that the electric fields of external field and local polarized field can be written as E1=E11E12; E2=E21E22; X1=X11X12; X2=X21X22; for applied and local electric fields, and for dipoles vectors on 2D plate. For free energy it is easy to write the expression by using binding energy and energies of dipoles in external field and local field: F=12NNiEjδQjEjxj, where Nj is a number of neighbor atoms, Ejand xjare the vectors of electric field and dipoles along the j direction. For j components of free energy along the axis 1 and 2 on 2D plate it is surely can be presented other form of such formula:


For analysis of deformation by applying the external electrical field it is clear to use ratio between free energy components for different bonding and directions:

F1F2=E1δQ1E11x11δQ2E21x21E2δQ1E12x12δQ2E22x22;and, according to the conservation law for a system shown on Figures 8, and 9 charge neutrality for all elements of system influenced by electric field equals to nonzero, it can be possible to write the following formula:


Figure 8.

Scheme of arrangement of dipole vectors along the axis.

Figure 9.

Scheme of nanocrystal’s destruction by external electric field.

Accordingly, the ratio between the polarization charges for two dipoles inside the electric field is given by

δQ1δQ2=E11x11+E12x12E22x22+E21x21;or δQ2=δQ1E22x22+E21x21E11x11+E12x12=Q1.

The first dipole is devoted to the description of electrical polarization properties of Si‐Si bond, but the second – Si‐O bond. By substituting the expression for Si‐O polarization charge in expression for free energy we can easily to obtain the following expression


The relation between the deformation values along the 1 and 2 axis can be used for analysis the translations consequences of modular group.

W1W2=E11x11AE21x21E12x12AE22x22;or W1W2=p+rq+s;where the first pair p,q describe the translation result applied field, but the second pair r, s causes the translation due to the local and external fields interaction with Si‐Si dipoles. Crystal phase destruction can be generated by using significant values of electric field combining with sufficient density of Si‐O dipoles inside the silicon film (Figure 10).

Figure 10.

Scheme of silicon nanocrystal film with orientation (111) and precipitate of crystal growth on silicon dioxide substrate with nuclear silicon-oxide dipole in the top of pyramid.

Transformations which create the low dimensional model structure are illustrated on Figures 11 and 12 and are explained as combination of modular group transformation of two-dimensional nanocrystal and Knop triangular transformations into Kantor dust set or fractal structure with low dimension. The surface and interface point defects and impurities cause the local electric fields which can generate by applying external field great values of field.

Figure 11.

Scheme of geometrical transformation the primary triangle with L and H parameters into triangle 1 with parameters L1 and H1.

Figure 12.

Scheme of charges density for anisotropic and non-homogeneous media.

6. Possible scenario of nanocrystal destruction: from bulk silicon to Kantor dust

Model of phase destruction by modular group substitution [26] which consists of arc series and Knop transformation of two-dimension area under arc through the triangular decomposition [27] to one dimensional structure. The down picture illustrates the creation of Kantor dust by dividing the triangular angle on the top and neglecting the area of triangle in the middle of primary triangle area. Such nonlinear triangular transformations can be caused by a point defects and impurities which were included in bulk silicon net of nanocrystal.


For such transformation the equations for triangular quantities XN = 2XN + 1 and for areas SN + 1 = 0.38*SN. The estimated value of Hausdorf dimension for such mathematical set is d=ln2ln2.63=0.72.

7. Classical and quantum mechanical models of charge and current densities

We have to propose the new model that is more suitable to explain electric properties of nanometrical scale media with strong anisotropic and non-homogeneous properties (see Figure 12). It will be necessary to describe the further possibility to design new nanoelectronic devices based on quantum conductivity properties and atomic scale sizes.


For current density of homogeneous media with charge density ρ in classical theory we usually use the formula


For ρ value of non-homogeneous anisotropic media it is possible to use the expression:

ρXXρXYρXZρYXρYYρYZρZXρZYρZZ. For current density the following form is given:

J=eSρXXρXYρXZρYXρYYρYZρZXρZYρZZvXvYvZ; where the absolute value of vector e=1.To determine the charge density and current density on nanoscopic scale the quantum mechanical approach is applicable by the way


It is clear, that such approach is approximate and can be applicable to study the electric properties of point defects.

Nonlinear polarization concludes as linear as nonlinear terms: P=αE, where α is a linear polarizability. Canonical equation for electro-magnetic fields can be written by using the Maxwell equations by a simplification of model and assuming that the first order of derivates are much more than others. The second supposal is that the solution can be explained as harmonic function. It is clear, that the equation for the field for second harmonic can be presented in the following form


It is supposed that currents which was created due to the electromagnetic field of second harmonic generation and induced in nanocrystals dominate in surface layers and grain boundaries. Because, such currents can be explained by the first term in equation and relate to absorption and emission of photons. Equation (17) can be written in suitable form:


Because, the surface current can play a significant role in nanostructured silicon film and the current density can be written as following J:


By a symmetrical form of wave functions Ψexpikxwtthe current value is zero. By a bonding of p orbital of silicon atoms for bonding and antibonding cases the energy gap between them is estimated as 10 μeV and their energetic location is closed to a bottom of conductivity band Ес-0.17 eV.


For two energetic levels (A and B) which are situated closed to each other the expression for the current density is following:


The values of currents for two energy states А and В are different due to the difference in their energies, and their occupations are also varied because they depends on Boltzmann distribution for unperturbed case, and for the laser excitation of carriers they distributed according to the Gauss distribution. Because, the current of charges depends strictly on energetic location of defects levels which are closed to the bottom of conductivity band of silicon. Therefore, the nonzero current is appeared because the field of SHG is applied in silicon nanocrystals. The current spectrum has a resonant energetic peak by the electron energy became equal to the energy of defect level:


The dipoles-field interaction causes the appearance of oscillations on frequencies Ω1=μ1Ehand Ω2=μ2Eh, which by reemission result in radiation with various harmonics such as 2ω±Ω1, 2ω±Ω2, 2ω±Ω1±Ω2. The estimates for the values а and b as a levels’ widths are following a=ΔEaΔEa+ΔEaga;b=ΔEbΔEa+ΔEagb, and by using the Erginsoy formula for semiconductors impurities it is possible to determine the sizes of local area for such kind of effect ~R3, where R is a size of area of surface SHG assisted currents are generated. Therefore, there is a surface current in nanocrystals which are generated by applying the SHG fields. Such currents are caused by repolarization of pairs of atoms of silicon by their distribution of the surfaces of all nanocrystals irradiated by laser light. By this irradiation the times of repolarization are much more than period of oscillation of electromagnetic field of laser radiation, because there is absorption of radiation.

The free energy for nanocrystal with volume V can be expressed as following: F=F0κVE22.By using determination of deformation of solid, it is possible to use the other definition for free energy by using the Lame coefficients η and υ [28]: F=F0+η2uil2+υuik2,By equating of two right parts from expression ((11)) and (12), it is easy to obtain the result: κVE22=η2uil2+υuik2.

The estimated value of free energy to destroy the silicon crystal phase is following F < 1.23 *10−4 J/cm2. There are various preliminary states between stabile Si‐Si bonds and broken bonds. The bond length distortion is varied from 0–15% of initial bond length. There are minima in potential energy for hydrogen which are varied from 1.3 eV to 2.3 eV. The Coulomb force driving the migration of charged impurity or defect can be written as F=ZeE;where the value of Z is an effective charge of defect, E is applied electric field. According to [28] the force due to the appeared interstitial atom is written as an Eshelby formula F=α43πr03Trσ,where Tr(σ) = σ11 + σ22+ σ33, σ is a stress, r0 is a radius of point defect. By this way it is clear that the electric field necessary to make stress can be evaluated as following: E=4απr033Zeσ11+σ22+σ33.

8. MGT model and mechanism of crystal phase destruction by local fields


Free energy is determined as energy Ejper one bond, Njis a quantity of neighbor atoms, Ejis electric field, and vector xjof dipole that is due to the inter atomic bond, such as Si‐O, for example. Inter atomic bond length (for Si‐O bond) can be calculated by using the following formula proposed by S. Hasegawa and co-workers [29] dSiO=dSiOnonpolarBSiδqSiBOδqO, where δqSiand δqOare partial charges for atoms, and BSi,andBSiare empirical factor that has its positive value. Therefore, for different directions is


For estimation the ratio between the energy of deformation and weak bond length which can be appeared by applying the external electric field we can use the Einstein relation for relation between the drift velocity and applying force vdrift=DkTF:μEd=kTDvN; where N is quantity of atoms that was locations were deformed, d is bond length of polarized silicon. The N value can be evaluated as following N=DkTvμEd. Here D is a diffusion coefficient, T is temperature, μis a dipole moment of bond. The model of kinetic behavior of densities of dangling bonds and weak bonds which are generated by local electric field can be described by using the following


where NH is the density of silicon-hydrogen bonds, NDBis the density of dangling bonds, and NWBis a density of weak bonds, W1is a rate of dangling bonds elimination and weak bonds creation due to the annealing of silicon, mainly, but the W2is a rate of generating of dangling bonds due to the hydrogen diffusion, W3is a rate of floating bonds generating due to the decreasing of density of weak bonds. It is assumed, that the values W2and W3are the same order of magnitude. For atomic diffusion coefficient the following formula is D=D0expEAkT; where EAis activation energy value and theoretical equation for D0can be written as follows D0=expΔSa2[30]; where n is a number of neighbor interstitial places of location, ν is frequency of vibrations in interstitial position, a is a lattice constant, αis a coefficient that value depends on the interstitial position, ΔSis an entropy ΔS=∂ΔFT, ΔFdifference in free energy that equals to the energy of activation. However, the total fraction of neutral interstitial is low compare with density of bulk atoms and according to W. Harrison studying [31], their number is approximately 1014 cm−3 by the melting temperature. Therefore, the diffusion as interstitial as dopant diffusion plays mainly if we observe only the hydrogen diffusion with energy of its activation 1.5 eV [32] and investigate the hydrogen moving through the sites with different potential caused by variation in electric charges.

Analytical solution of system of differential equations results in the following expression for density of dangling bonds

NDB=NDB0expA2texpA24Bt+expA24BtCB;where the values are given.

A=W2+2W1+W3;B=W2W1+W3W1W3;C=W1+W3W1NH+W1W3NFB. By the relation between values W2 and W3 as W2 = 0.001* W3, or W2 ➔0, the coefficients A, B, C have view.

AW2=0=2W1+W3;BW2=0=W1W3;CW2=0=W1+W3W1NH+W1W3NFBand the solution in this case can be written as.

NW2=0=NW2=00expW12+W324W112W3t1+exp4W12+W32t+W1+W3W3NH+NFB;where W1>W3/3should be realized for increasing the density of dangling bonds, and in this case the density of dangling bonds has a slow rising over the time (as it is shown in Figure 13). By using the Runge–Kutta method of 4-th order it is easy to calculate the evolution of density of dangling bonds by the local electric field in several regimes for NDBt, NWBtas normalized values of density of dangling and weak bonds.

Figure 13.

Evolution of density of dangling and weak bonds for different degrees of atomic sites order apart from the central of local electric field. It is assumed, that the values are W 2 = W 3 . It is seen, that there is a reconstruction of the order by initial increasing of weak bond density, and redistribution of particles’ local places according to the surrounding density values. It is assumed that this is a creation of order of amorphous phase. The third curve illustrates the evolution of density of perturbed bonds of next degree of order from the distance of a local electric field center. Oscillations of density of dangling bonds reflect the reconstruction of bonds and destruction by the annealing, hydrogen migration, according to the initial order of fixed atomic places of crystal structure and their new amorphous local placing. Over the period of time T the oscillations are suppressed by redestribution of initial impulse created by local field center and caused the increasing of weal bond density value. By this way the deformation penetrates in the surrounding of local area that will be larger than before.

It is known the model of defects generating by light irradiation in amorphous silicon which was proposed by the scientists of Ames Laboratory [33] which calculated the evolution of density of dangling bonds according to their proposed model. It is seen, that the evolution which was shown in Figure 14 has the same increasing tendency as evolution stimulated by light irradiation. By the values W2=0.001* W3the evolution changes its sharp increasing to the slow behavior of density of dangling bonds function as it is seen in Figure 14 .

Figure 14.

Evolution of density of dangling bonds for different degrees of atomic sites order apart from the central of local electric field, by the ratio between the values W 2 =0.001* W 3 . This scenario can describe the poor hydrogen contamination of silicon network.

By the ratio between rates W2and W3equals to 100 the calculating evolution of density of dangling bonds can be presented as that was shown in Figure 15.

Figure 15.

Evolution of density of dangling bonds for different degrees of atomic sites order apart from the central of local electric field, by the ratio between the values W 2 =0.01* W 3 . It is clear, that there is amorphization of the solid because the atomic sites order was damaged. It is assumed that such scenario can be explained the difference in quantity of dangling bonds near the local field center NDB Local and surrounding, NWB Surrounding, NWB Surrounding> > NDB Local. The conservation law can explain that the small number of particles which are placed near local field center cause the changes in impulses of numerous atoms in surrounding area.

By the ratio between rates W2and W3equals to 10 the calculating evolution of density of dangling bonds can be presented as that was shown in Figure 16.

Figure 16.

Evolution of density of dangling bonds for different degrees of atomic sites order apart from the central of local electric field, by the ratio between the values W 2 = 0.1* W 3 .

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Dmitry E. Milovzorov (May 23rd 2018). Local Electric Fields in Dielectric and Semiconductors, Electric Field, Mohsen Sheikholeslami Kandelousi, IntechOpen, DOI: 10.5772/intechopen.74310. Available from:

chapter statistics

476total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Local Electric Fields in Dielectric and Semiconductors: Part II

By Dmitry E. Milovzorov

Related Book

First chapter

Numerical Analysis of Three‐Dimensional MHD Nanofluid Flow over a Stretching Sheet with Convective Boundary Conditions through a Porous Medium

By Stanford Shateyi

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us