Local Electric Fields in Dielectric and Semiconductors Local Electric Fields in Dielectric and Semiconductors

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.


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 ð Þ h i¼ ε A w ð ÞÀεB w ð Þ ε A w ð Þþ2εB w ð Þ the Maxwell-Garnet formula is 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 10 7 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.
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]: ð Þ is the field at the i-th monomer, r ij is 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]: where H 0 is the Hamiltonian pure semiconductor for electron matrix elements of transitions between own conductive and valence bands, but H D and H A are the exchange energies matrix elements due to the donor and acceptor states, and for them it is clear be presented the following expressions: EÀε l Àσ l δ ij ; where D ki and A li are coupling between k-th donor state and l-th acceptor state, but the value σ k, l ¼ 1 are the selfenergy corrections which depend on locality of sites, and θ k, l ¼ E À ε k, l , and γ k, l are 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 SiO x , 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 crystalamorphous 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.

Experimental researches of field-assisted destruction of silicon nanocrystals
Nonlinear polarization associated with the phonons can be written as P ¼ 1 Using these equations it is possible to obtain the wave equation for field: Raman effect is result from the interaction of an electromagnetic field and optical phonon mode. The vibration wave < Q >¼ Q 2 exp i ω v t À k v z ð Þþc:c: ð The polarization at the Stokes The microcrystal wave function is a superposition of Eigen functions with k vectors centered at k 0 . We suppose that k 0 ¼ 0. If we are using the weighting function as Gaussian we will have the normalized first-order correlation function 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]: k ¼ 2π a k , ais lattice constant, k is dimensionless value,Γ 0˜3 :6cm À1 line width of the Si LO phonons in c-Si. The dispersion of the LO phonon in c-Si ω 2 k ð Þ ¼ A þ Bcos πk =2 ð Þ, where A = 1.714 Â 10 5 cm À2 and B = 10 5 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 N SiO = 210 21 cm À3 and for crystalline silicon the density of silicon bonds N SiÀSi = 510 22 cm À3 the lattice constants are following: a 2 = 0.98 a 0 аnd a 1 = 0.996 a 0 .
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 The deformation potential causes the appearance of effective acoustical charge ( Figure 2).
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 SiO 2 [18] the Raman spectrum of SiO 2 has the variation modes D1 (at 490 cm À1 ) with defects and activation energies 0.14 eV and pure mode w 1 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: where N SiM , δ SiM are the density of the bonds and partial charge for SiM bond (O or H), respectively ( Table 1).
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 SiH 2 [20]. The two atoms of hydrogen in the SiH 2 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 10 19 to 10 21 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 510 22 cm À3 , but densities of SiO and SiH bonds are 10 21 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 P SiO = 0.04 D and P SiH = 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 SiO x creation by using the Raman scattering data which correspond to the results reported in work [21].

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: 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 Si 1 and Si 3 atoms without oxygen) can be written asĤ ¼ . This matrix can be transformed into next triangu- Here, the matrix elementsÂ 22 Because, for the creation operator . Therefore, the dispersion curve is written as ω phonon = ω 0 AE ξ, where ω 0 dispersion curve without applied field (for energy 0.14 eV the estimated additional frequency is less than 3.4 Â 10 13 Hz, and additional frequency ξ= α 2 λ11 by 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).

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 ¼ PE ℏ is 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).
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 10 20 -10 21 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.
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(ω, r) value is a local field factor of film with crystalline volume fraction equals to r = 70%.
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 E 1 (for unperturbed state) and E 2 (perturbed state): The changes in Eigen values from the E 1 and E 2 by the field E = 0 to the new values of energies are following Figure 6. SHG spectra for P in ➔P out laser scheme and χ (2) xxx component of susceptibility.
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.
. Therefore, by a small perturbation of system ε ≤ E 1 À E 2 ð Þ , 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 ¼ 1 4πε0 0:2e 2 r 2 ; 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 ¼ HM kT k∇T, where H M is an enthalpy or energy of transport by heat, ∇T is 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∇σ À σ∇u ð Þ ds. 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 ¼ 1 2 Mω 2 x 2 , the energy reduction due to the linear electron-lattice interaction as E ¼ Ax.

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.
It is assumed that the electric fields of external field and local polarized field can be written as ; 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: , where N j is a number of neighbor atoms, E j ! and x j ! are 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: E2ÀδQ 1 E12x12ÀδQ 2 E22x22 ; 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: Accordingly, the ratio between the polarization charges for two dipoles inside the electric field is given by E11x11þE12x12 ¼ ÀAδQ 1 : Figure 8. Scheme of arrangement of dipole vectors along the axis.
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.
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).

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.  S ¼ 1 2 HL:, For such transformation the equations for triangular quantities X N = 2X N + 1 and for areas S N + 1 = 0.38*S N . The estimated value of Hausdorf dimension for such mathematical set is d ¼ ln2 ln2:63 ¼ 0:72.

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 r in classical theory we usually use the formula For r value of non-homogeneous anisotropic media it is possible to use the expression: 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 Ψ ≈ exp i kx À wt ð Þ ð Þthe 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 ¼ h , which by reemission result in radiation with various harmonics such as 2ω AE Ω 1 , 2ω AE Ω 2 , 2ω AE Ω 1 AE Ω 2 . The estimates for the values а and b as a levels' widths are following a ¼ ΔEa ΔEaþΔEa g a ; b ¼ ΔE b ΔEaþΔEa g b , and by using the Erginsoy formula for semiconductors impurities it is possible to determine the sizes of local area for such kind of effect~R 3 , 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 ¼ F 0 À κVE 2 2 : 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]: By equating of two right parts from expression ((11)) and (12), it is easy to obtain the result: The estimated value of free energy to destroy the silicon crystal phase is following F < 1.23 *10 À4 J/cm 2 . 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 ¼ α 4 3 πr 0 3 ∇Tr σ ð Þ, where Tr(σ) = σ 11 + σ 22 + σ 33 , σ is a stress, r 0 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: 8. MGT model and mechanism of crystal phase destruction by local fields where δq Si and δq O are partial charges for atoms, and B Si ,andB Si are 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 v drift ¼ D kT F: 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 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 N H is the density of silicon-hydrogen bonds, N DB is the density of dangling bonds, and N WB is a density of weak bonds, W 1 is a rate of dangling bonds elimination and weak bonds creation due to the annealing of silicon, mainly, but the W 2 is a rate of generating of dangling bonds due to the hydrogen diffusion, W 3 is a rate of floating bonds generating due to the decreasing of density of weak bonds. It is assumed, that the values W 2 and W 3 are the same order of magnitude. For atomic diffusion coefficient the following formula is D ¼ D 0 exp ; where E A is activation energy value and theoretical equation for D 0 can be written as [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, ΔS is an entropy ΔS ¼ À ∂ΔF ∂T , ΔF difference 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 10 14 cm À3 by the melting temperature. Therefore, the diffusion as interstitial as dopant diffusion plays mainly if we observe only the hydrogen diffusion with 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. 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 where the values are given.
By the relation between values W 2 and W 3 as W 2 = 0.001* W 3 , or W 2 ➔0, the coefficients A, B, C have view.
A W2¼0 ¼ 2W 1 þ W 3 ; B W2¼0 ¼ ÀW 1 W 3 ; C W2¼0 ¼ W 1 þ W 3 ð Þ W 1 N H þ W 1 W 3 N FB and the solution in this case can be written as. 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.
N H þ N FB ; where W 1 > W 3 =3 should 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 N DB t ð Þ, N WB t ð Þ as normalized values of density of dangling and weak bonds.
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 W 2 =0.001* W 3 the evolution changes its sharp increasing to the slow behavior of density of dangling bonds function as it is seen in Figure 14 . By the ratio between rates W 2 and W 3 equals to 100 the calculating evolution of density of dangling bonds can be presented as that was shown in Figure 15.
By the ratio between rates W 2 and W 3 equals to 10 the calculating evolution of density of dangling bonds can be presented as that was shown in Figure 16.