Local electric fields appeared in dielectric and semiconductors due to the destruction of symmetry, creating the vacancies, point defects, and chemical impurities in material. By increasing external electric field value, numerous structural changes will be generated. Point defects in silicon films were characterized by using electron-paramagnetic resonance spectroscopy and laser picosecond spectroscopy. EPR spectroscopy provides the detailed microscopic information about the point defects in silicon films, which have unpaired electrons. The switching effect for A or (VO)− defects was observed by applied bias voltage for nanocrystalline silicon films. The intensities of EPR signal change according to the switching of position of oxygen incorporated in silicon (111) from one pair of silicon atom to another pair of silicon atoms. In this case, the interaction between px orbital is changed to pz orbital. The proposed mechanism of states’ switching for the vacancy (VO)− is based on a model of coupled dangling bonds (DBs) of silicon atoms and interaction between atomic orbitals. The crystal phase destruction in nanocrystalline silicon film by applying external electric field was investigated by Raman spectroscopy. The possible mechanism of phase creation was proposed.
- local field
- point defect
- dangling bonds
- Raman spectroscopy
- second-harmonic generation
- silicon films
Local electric fields appeared in dielectric and semiconductors due to the destruction of symmetry, creating the vacancies, point defects, and chemical impurities in material. By increasing external electric field value, 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. Studying the nature of local electric fields will open new tendencies in electronic device producing, from one side, and help to change materials’ properties according to our needs, from another side. The quantum size effect for nanocrystalline silicon films with broad size distribution was observed by PL and SHG spectral measurements. The optical properties of films are likely to depend on the size distribution of nanocrystals, along with an effect of hydrogen termination on the surface of nanocrystals. The studying of physical properties of such nanoscale objects as nanocrystals, point defects and chemical bonds is very important for understanding the appearance of local fields in silicon films, as it is necessary to investigate them and to detect the physical phenomena which appear due to the local fields. For small particles such as electron, the wave features play a significant role because the de Broglie wave length is given by where p is impulse and h is Planck constant. For electron rotating around nucleus of atom with circulation rate 3.54 × 106 m/s, Å. For oxygen atom at room temperature and the rate of motion 250 m/s, the de Broglie wavelength is Å. The probability of occupancy of each energy state is defined by wave function as a solution of the Schrodinger equation where H is full energy of system named as Hamiltonian. For example, for the quantum well with the size L = 100 Å, the solution of this equation is satisfied to a discrete energy according to the following expression where numerical data are 3.37 meV for n = 1, 13.48 meV for n = 2, 30.33 meV for n = 3, and 53.92 meV for n = 4.
2. Quantum objects in silicon films
The probability to find the particle in each state can be calculated as Hamiltonian of system can be written in following form: where
2.1. Nanocrystals in silicon films
The combination of laser picosecond spectroscopy with Raman scattering measurements is perspective for nondestructive analysis of composite thin films and silicon surface with complex relief. The method of laser picosecond spectroscopy is a powerful tool for surface analysis with high sensitivity. The pump-probe scheme is similar to the interferometric schemes with great possibility to measure the refractive index changes. For composite thin film measurements, such as poly-Si or porous silicon, the local factor field plays a great role because the crystalline volume fraction (or porosity) and dielectric functions’ values of film’s components are important. Also, due to the great amount of scattering centers (e.g., nanocrystallites or microcrystallites spatial distribution in film according to Gaussian distribution), the role of photo-acoustic distortion is negligible. However, the non-equilibrium perturbation of career density and its decay due to the diffusion or recombination (in case of great porosity or small value of ρ) are detected surely by using picosecond technique. In addition, it is necessary to note that the electron relaxation from surface and defect states can be reflected by optical response evolution, too. Also, by carrying these measurements, it is important to know the surface morphology. For elimination of the surface scattering background the correlation function value was measured.
2.1.1. Nanocrystal detection by Raman spectroscopy
The Raman spectra for poly-Si films consist of a narrow line near 520 cm−1 arising from a crystalline phase and a broad line around 480 cm−1 from an amorphous phase. The ρ value was estimated from the ratio of the Raman integrated intensity for the crystalline component to the total intensity, using the ratio of the integrated Raman cross-section for the crystalline phase to that of an amorphous phase. The penetration depth of the incident Ar-ion laser radiation (λ = 488 nm) into silicon is within about 0.3 μm.
The crystalline volume fraction was estimated using Raman scattering measurements from the ratio : , where
It was assumed that this change is because the bond angle variation value decreases by the starting crystallization process in a-Si film. We assume that the Raman data of hydrogenized amorphous silicon film resulting in the spectral peak around 445–447 cm−1 corresponds the LO mode but it is 480 cm−1 for TO mode. For the higher structural relaxed silicon thin film, by high level of hydrogen dilution, the TO mode reflects in 490 cm−1 value of peak position. The spectral peak width changes from the value of 40 cm−1 for a-Si to the 70 cm−1 for a-Si:H.
Figure 1 shows the various types of decompositions of Raman spectrum for microcrystalline silicon film with thickness around 100 nm. This film was deposited at 380°C. By using this procedure we use 3, 4, and even 9 peak approximations. The figure illustrates the three-peak decomposition with c-Si-related spectral line around 520 cm−1, nanocrystalline or intermediate spectral line in the range 500–510 cm−1 and for amorphous silicon it is 480 cm−1. The next figure, Figure 1b, shows the Raman data with noise signal for comparison, and fourth spectral component (~490 cm−1) is related to hydrogenated amorphous silicon network. The estimation of crystalline volume fraction was 27%. The hydrogen terminates all the dangling bonds (DBs). The spectral characteristics of amorphous-related peaks are changed by the hydrogen saturation of dangling bonds.
2.1.2. Time-resolved picoseconds laser spectroscopy of silicon nanocrystalline films
A mode-locked YAG:Nd3+ laser radiation with wavelength 532 nm was used as an optical pump of media, but the second-harmonic radiation (λ = 1064 nm) was used for probing the sample’s Figure 2. The correlation between the reflected signal intensity and probe laser beam for a-Si and p-Si films surface is given. The pulse duration was 120 ps. The pulse repetition rate was 100 MHz and frequency of Q-switched modulation of second-harmonic radiation was 6.2 MHz.
The correlation function of reflected signal intensity of probe laser beam was detected by means of pump-probe laser scheme where
where E is transition energy and < i| and |j > are the initial and final states of electron transition. The time evolution reflects the superposition of ensemble of different kinds of transitions such as band-to-band transition, transition through the surface or defect states. The width of the energy level depends on the chemical bond structure of silicon surface. The most important fact is appearance of quantum beats, which reflect the optical response from two-neighbor level. We suppose that the different fractions of oxygen incorporation in silicon film cause the production of neighbor levels inside the band gap with the spectral width: a-Si reflects the presence of complex SiO configuration. The figure shows the correlation function as a function of delay time for oxidized Si (111) surface and poly-Si silicon film with <δ> = 9.7 nm. It is expected that the oxide defect levels and surface state level (for Si (111) surface) are responsible for the observed oscillation during the exponential decay of electron density. We suggest that such oscillations are the same as quantum beats of levels of occupation by the laser time-resolved spectroscopy of molecular levels with the period of ΔΤ. Also, it is seen that exponential decrement for poly-Si film is weaker than the other. It is assumed that the correspondence between the values of decrement and width of oxygen-related levels in the band gap of silicon is proved and the magnitude of correlation function oscillations illustrates oxide amount on the surface. Also, our assumption has strong evidence from the FTIR data for aged poly-Si films.
Intervalley scattering rate is given by 
Dipole moment can be written in form 
In conclusion, the new method of laser spectroscopy for testing nanocrystal silicon structures and measurements of band gap energy and energy position of defect levels inside the band gap is proposed. Also, the competition between the electron diffusion and recombination processes in poly-Si films by the transmitted probe-pump laser scheme is analyzed. It is clear that tunneling process between the silicon nanocrystals through the thin oxidized layer plays a great role and it is possible to design a new electronic device.
2.2. Point defects in silicon
For several decades, a-Si:H films are intensively studied by different spectroscopic and microscopic methods. Most of the a-Si films contain some amount of nanocrystals. This is because the PECVD method to prepare nc-Si films and a-Si films is the same but differs only in deposition conditions such as substrate temperature, gas mixture, working pressure, and RF power. Figure 5 illustrates the CVD-prepared a-Si film with gas silane by 1 sccm flow rate. The RF power was 20 W and temperature of substrate was 25°C. It is clear that the sharp peaks relate to the various fractions of silane dissociation: (SiH)x; (SiH2)x; (SiH3)x; (SiH4)x. It is supposed that their responses are around spectral lines: 1.56, 1.65, 1.73, and 1.83 eV. Also, the broad spectral peak around 1.75 eV is related to the tiniest nanocrystals in a-Si film.
EPR spectroscopy provides the detailed microscopic information about the point defects in silicon films, which have unpaired electrons. It is seen in Figure 6, that there is a dangling bond signal on EPR spectrum. Therefore, the a-Si films contain dangling bond defects and are disordered due to the hydrogen termination of dangling bonds. The surface diffusion and desorption processes by growth are low at low temperatures of deposition.
2.2.1. Dangling bonds
For detail explanation of dangling bond (DB) defects several kinds of models were proposed such as defect-pool model of Powell and Dean . The main principle of this model is the following. There is equilibrium between the quantity of dangling bonds and weak bonds. The chemical reactions that cause the appearance of dangling bonds are the following: WB (2DB); SiH + WB (DB + SiH) + DB; 2SiH + WB (Si-H-H-Si) + 2DB. Thin film deposition at low temperatures of substrates results in the appearance of dangling bonds defects that can be terminated by hydrogen. The process of nucleation of crystal grains required the minimum quantity of silicon atoms.
These paramagnetic centers appear when there is an unpaired electron of the dangling bond of silicon atom that bonded with three silicon atoms (see Figure 7).
For coupled dangling bonds of silicon atoms on the surface (111) there is an interaction between atomic orbitals (as it is seen in Figure 8). The degeneracy of the level according to rules for wave functions of one silicon atom with free atom
2.2.2. Vacancies and oxygen incorporated in silicon (111)
At first time the term A defect (or center) in silicon was used by Watkins and Corbett in 1961 . They studied the irradiated by 1.5 MeV electron beam silicon with current density 2.5 μa/cm2. The role of oxygen in A center creation is significant due to its definite bonding with silicon atoms. Spin resonance is caused by unpaired electrons which is trapped by splitting of the atomic orbital of pair-coupled silicon atoms. The figure illustrates the appearance of A defect paramagnetic centers in inter-grain area of nanocrystalline silicon films experimentally detected by means of the EPR spectrometer technique.
Figures 9 and 10 show the switching effect for A centers by applied bias voltage for fluorinated nanocrystalline silicon films. The intensities of EPR signal change according to the switching of position of oxygen incorporated in silicon (111) from one pair of silicon atom to another pair of silicon atoms. In this case the interaction between
Figure 11 presents the various pictures of potential energy for oxygen atom bonded with couple Si atoms as a function of a distance by classical approach. It is clear that the distances that characterize the motion of oxygen atom are compared with its de Broglie wave length (2–4 Å) because there is a possibility to use quantum mechanical treatment of moving oxygen atom caused by the electric field from one silicon atom pair to another. Figure 12 shows the scheme of Si-O-Si bridge with coupled Si-Si atoms as the quantum oscillator with eigenvalues of energy of the oxygen atom.
The oxygen atoms and dimers are incorporated in the silicon grain boundary and have weak covalent bonds with silicon. But the activation energy for molecular diffusion is low, 0.3 eV, in contrast with the activation energy value for atomic diffusion (1.3 eV). The energy of Si-O-Si bridge interaction can be written in form of Morse function
where U0 is the energy of Si-O (4.5 eV) and α is coefficient. The oxygen atom in A center oscillates between two points of stable positions. Transition of probability for oxygen atom through the potential barrier by applied external electrical field U(r) = U0-Er is given by
where ν is a frequency of oscillations of oxygen atom. Wave function of the oxygen atom is the following
where φ1and φ2 are the wave functions of pure states 1 and 2 and Г1 and Г2 are the widths of the levels, respectively. We assume that all the Г values are approximately equal to each other. Raman scattering data help to determine the dipole Si-Si orientation along the laser E field ax of incident radiation on the silicon surface by applied electric field switching of O atom spatial position. Applied electric field causes the tunnel of oxygen atom from one pair of coupled silicon atoms to another with perpendicular axis of dipole orientation. The annealing procedure can assist to restore spectral characteristics related to the primary-ordered atomic position due to the minimums in potential energy diagrams.
2.2.3. Physical model of vacancy switching by applied electric field
The role of applied electric field in order–disorder transition can be surely understood by using several mechanisms of electric dipole creation as it is seen in Figure 10. The first mechanism is the changing length of silicon-oxygen Si-O bond or silicon fluorine (Si-F) bonding and atomic transfer of oxygen inside the vacancy defect which causes the reorientation of VO or VF defects. Such bonds are strongest in silicon film because of high-electron affinity values for O and F atoms. The dipoles Si-O and Si-F have the great value in their polarized charges because the interaction between them and electromagnetic field is sufficient, and it results in their efficient reorientations because the dipole in external field changes its orientation to compensate it.
The second mechanism of order–disorder transition in the silicon film was initially proposed by Watkins and Corbett , which is devoted to hydrogen atom reorientation around the vacancy by the temperature above 200°C as the possible EPR data explanation. For nanostructured silicon thin film that contains a great amount of tiniest nanocrystals, the hydrogen atoms are distributed in grain boundary of nanocrystal silicon grains. But by electric field application to the thin film along the surface, the smallest nanocrystals disappear and, partially, the crystal structure is destroyed. By annealing, the opposite tendency is observed. The third assumption is the creation of excited hydrogen molecule inside the silicon film by annealing to the temperatures 150–200°C: Si-H + Si-H→Si-Si + H*2 . It is clear that the silicon nanocrystal fraction decreases but the broad spectral line for amorphous phase significantly increases. The fourth mechanism is a deep valence hole and negatively charged OH- created by laser beam irradiation  (Figure 13).
The pump-probe laser picosecond spectroscopy data with phonon-stimulated stress pulses  confirm the existence of the quantum interference of closed neighbor levels resulting in quantum beats phenomena (see Figure 3). The estimation of the energy difference of level position is around 5 μeV. Picosecond acoustic longitudinal stress pulses cause the small displacement of interface media surrounding the silicon nanocrystals because the Si-O-Si bridges were slackened and the oxygen atom oscillates between two stable excess states of four coupled pairs of silicon atoms. The main principle of switching is a motion between two silicon coupled atoms with oxygen from one pair of silicon atoms to another. R. Biswas and coworkers proposed the flip model for hydrogen incorporation in amorphous silicon with two stable opposite positions for hydrogen atoms with axial symmetry relative to the bonding silicon atom .
The energy of Si-O-Si bridge interaction can be written in form of Morse function
where U0 is the energy of Si-O (4.5 eV) and α is the coefficient. The oxygen atom in A center oscillates between two points of stable positions. The change of the spatial position of oxygen atom can be described by using simple model of oscillator
where the eigenfrequency of oscillating oxygen atom can be expressed as ; , k is Boltzmann constant. The solution of equation by the resonance conditions (by the ) can be written as
The eigen frequency value for Si-O bond was estimated by means of the FTIR spectroscopy technique. It equals approximately to 3 * 1013 Hz. The resonance will be by the equality between the phonon frequency Ω and ω0 value.
2.2.4. Model of vacancy with oxygen as two coupled oscillators
For the vacancy-oxygen atom (VO) defect center the model of two coupled oscillators is much more convenient:
and the energy of atom interaction can be expressed as the following for two pairs of silicon atoms which are oriented along the x and z axis:
where x1 and x2 are the displacement values for oxygen atom positions relatively, Si-SiA and Si-SiB pairs.
The solution of the system of two two-order differential equations contains two eigen frequencies w1 and w2. The oxygen atom oscillates in A center with two different frequencies corresponding two energy levels for oxygen atom. The quantum beats are observed in picosecond laser spectroscopy experiments shown on Figure 19 caused by quantum interference of two energetic states of oxygen atom.
2.2.5. Dipole moments of Si-O bond and Si-O-Si bridge
By using the Slater atomic orbital for Si and O atoms which were written in my previous work we can easily estimate dipole moment for S-O configuration. Dipole moment by the zero displacement is ; where , 0 < r ≤ ξ, ζ is covalent radius of oxygen atom, x is bond length. For 3pSi-2pO orbital bonding we can ; ;
By r = 1.6 Å, r/a0 = 3, we can estimate the dipole moment value as μSi-O = 1.12 D.
The comparison of calculated overlapping integrals as a function of distance is shown in Figure 14. It is seen that the Si-Si orbital is more significant even at large distances than Si-O.
For Si-O-Si bridge we can obtain
We can evaluate the dipole moment by zero displacement . For Si-O-Si bridge the angle (see Figure 15) between two bonds is approximately 90–180°. It is clear that by α➔180o the vector value of dipole moment is zero, but by α➔90° the dipole moment can be estimated as .
By applied electrical field the dipole moment value for Si-O-Si bridge is defined by displacement evolution: . In the case of μ < μcritical there is no localized bonded state. It is clear that for μ0Si-O > μcritical there is localized bonded state for electron, but for μ0Si-O-Si < μcritical such a state cannot exist because the Si-O bonding by zero displacement generates localized electron states on Si (111) surface which can cause troubles for writing and storing procedures of memory cell. By the annealing of Si (111) substrate to a temperature of more than 250°C we can eliminate the Si-O bonds and produce the strong siloxane bonding. Therefore, it is possible to prepare the oxidized Si (111) surface with strict siloxane Si-O-Si bridges without localized electron states on the surface. Applying the electrical or even acoustical field we can stimulate A center switching inside the film that causes the appearance of localized electron states. The polarization wave generated by laser radiation is given by
We can estimate the yield of interference of two polarized waves with shift of phase’s π.
The spatial distribution of charges for polarized molecular bridges Si-O-Si on the clean Si (111) surface is used for electrical charge trapping. The potential energy barrier for electron trapping can be estimated as
where Δ d is the cell size and λ is the wavelength of polarization wave. By applied electrical field 107 V/m W ≈ 100 μeV and by the value of field 108 V/cm is 1 meV.
Figure 16 illustrates the point defect on silicon (111) surface caused by a silicon atom vacancy and oxygen atom incorporation into the silicon network. There are three silicon bonds which are connected with oxygen atom through the
184.108.40.206. Si3H-σ vacancy defect
The molecular-like spectrum of electron states for hydrogenated cluster Si3H−σ has a series of electron states in the range from 0.7 to 1.0 eV. . The eigen frequencies of oscillator Si3-H are conjugated with the electron states for the same oscillators and results in appearance of broadening in the electron spectra. The defect can be described as tied three Si3-H oscillators and one Si3- oscillator. It is clear that the quantum properties of defect are reflected in their molecular orbital spectral positions. It should be observed surely for nanostructured silicon films with structured cell less than 1 nm or for amorphous hydrogenated silicon film. In such kind of media we surely can change the spectral positions of electron states by using bias voltage.
220.127.116.11. Si3-F vacancy defect
Silicon functional groups can be determined as the following: Si-Si, Si-H, Si-H2, Si-H3, SiH2-F. The silicon bonding is arranged according to
18.104.22.168. Si-Ox clusters
We assume that the defects are localized in Si-SiO2 interface with complex composition of SiOx. The oxygen atoms and dimers are incorporated in the silicon grain boundary and have weak covalent bonds with silicon. But the activation energy for molecular diffusion is low, 0.3 eV, in contrast with the activation energy value for atomic diffusion (1.3 eV). The increase in oxygen concentration results in increasing SiO bond length. For the sample with smaller grain size the hydrogen termination of dangling bonds is greater, and density of oxygen is low. But the film with large grains contains a great amount of oxygen in Si-O-Si and SiO2 compositions. We assume that the exponential electron density decay is due to the exciton state decay. The short relaxation time can be described as possible relaxation process through the oxygen incorporation-related defect states, where density of states is the following:
The oxygen atom bonding with silicon atoms inside crystallites forms the dioxide composition along with defect production. Figure 17 shows the two kinds of oxidized surfaces: Si (100) surface annealed by 1100 K and oxidized poly-Si film with a great value of <δ> = 17.6 nm. The oxide-silicon configuration for both surfaces is different. The Si (100) surface contains mainly the SiO2 component with a high density of oxygen concentration but the poly-Si film is weaker oxidized with the presence of SiO0.5 configuration.
Figure 18 illustrates the oxygen incorporation into silicon grain boundary after transformation of the silanol groups into strong siloxane bonds.
3. Field-assisted crystal phase creation from amorphous
3.1. Fractal character of restored crystal phase
The first mechanism that can be used for explanation of crystal phase restoration phenomenon in silicon film is a difference in structural orders of silicon crystal structure and polysilane chains which are arranged into conglomerates after hydrogen migration, and further incorporation is transformed into the silicon net. The hydrogen surface diffusion by annealing also changes the phase fraction ratio in the silicon film. Figure 19 shows the microphotography of silicon microcrystal that contains fractal structures of silicon bonding. By applying the external electrical field, the dipole orientations change to compensate the external field. Because the dipoles’ reorientations cause numerous dislocations inside silicon nanocrystal, it is supposed that the order of crystalline structure inside the silicon film as grown usually is the same as order of fractals which is shown in Figure 19. The microphotography of the microcrystal of silicon on the glass substrate that is shown in Figure 19a demonstrates the hexagonal structure of the crystal along with fractal structures which are illustrated in Figure 19b. Because there is the necessity to analyze the mechanism of field-assisted phase transformation as destruction or phase restoration of native fractal nanoscale structures of silicon nanocrystals, for fractals, the dimension of Hausdorf can be surely evaluated as , where N(s) is the number of cubic volumes from which it is possible to substitute the rectangular nanocrystal with crystal orientation (111). Such a value of dimension equals 2.5. For determination of any structure the formula for Hausdorf metrics that is given by.
is suitable, where ε is radius of the open ball with centers in points X and Y. The space limit for the fractal structure can be evaluated using the silicon bonding size 2.5 Ǻ. For small 3D objects such as nanocrystals, the density of substitution can be written as
The second proposed mechanism of silicon crystal phase restoration can be illustrated by means of Figure 19. It is shown that the dislocations and dipoles in grain boundary of nanocrystal under the influence of external electric field cause the phase destruction into the arc series according to the modular group scenario. Such transformations cause the appearance of numerous arcs which will be destroyed after further triangular Knop scenario of transformation of the 2D silicon structure to 1 D silicon structure or polysilane chain. For the fractal structure the Hausdorf dimension can be varied and is less than 1. The down picture on Figure 8 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. For such transformation the equations for triangular quantities XN = 2XN + 1 and for areas SN = 2.63d*SN + 1. Comparative analysis of last two equalities results in 2.63d = 2. The estimated value of Hausdorf dimension for such mathematical set is 0.72, because, there is a decrease in dimension by the destruction mechanism for media. It is obviously that the field-assisted transformation of silicon nanocrystals into objects with low dimension can be followed by the way of a polysilane chain creation .
3.2. Raman spectra of silicon film with the created crystal phase
Figure 20 illustrates the Raman spectra of amorphous silicon film which is deposited by magnetron sputtering atoms of Pt and further applying the external electric field. The scheme of experiments is shown in the figure. There is a crystal phase that was created by fixed value of applied electric field due to the charges accumulated by Pt clusters. Change in free energy by local field results in the appearance of silicon stripes of the crystalline phase in amorphous silicon.
Free energy functional can be expressed as follows
, where is a phenomenological parameter and is a dimensionless average density. The solution of equation is
or finding the minimum of the functional results in the investigations of crystal phase creation. For 1D in the uniform phase the expression can be used as a possible solution. Figure 21 shows free energy as a function of reduced particle density at the melting point of Fe according to the phase-field crystal model estimations.
3.3. Renormalization group of translation for crystal phase creation
The resolving the equations for locality function by the free energy functional minimizing :
for locally conserved, and for locally non-conserved, from which we can easily find the equation for estimating the value of energy for order–disorder transition:
also, it is assumed that the evaluation value of free energy (of heat) equals the energy of local field force to arrange in order all the atoms which should be placed on their sites: .
The solution of nonlinear Swift-Hohenberg equation also can be presented as a linear combination of solutions if we derive the terms of as partial solutions and further combine them. The local function can be transformed into a series according to perturbation theory as follows , where the value is a steady-state solution, and next values are the next order of perturbed decomposing procedure. For such a purpose it can be used as a first order of the solution that can be expressed as: and amplitude , where is a renormalization constant .
Tsu R, Gonzalez-Hernandez J, Chao SS, Lee SC, Tanaka K. Critical volume fraction of crystallinity for conductivity percolation in phosphorus-doped Si: F: H alloys. Applied Physics Letters. 1982; 40:534. ISSN: 0003-6951
Lucovsky G, Nemanich RJ, Knight JC. Structural interpretation of the vibration spectra of a-Si:H alloys. Physical Review B. 1979; 19:2064-2073. ISSN: 0163-1829
Yu P, Cardona M. Fundamentals of Semiconductors. Berlin: Springer; 1996. p. 617. ISBN: 3-540-61461-3
Kahan VD. High frequency hopping transport in solids and dielectric transparency in compensated semiconductors. JETP. 2000; 117:452-456. ISSN: 1063-7761
Powell MJ, Deane SC. Improved defect-pool model for charged defects in amorphous silicon. Physical Review B. 1993; 48:10815-10826. ISSN: 0163–1829
Watkins G, Corbet JW. Defects in irradiated silicon I: Electron spin resonance of the Si-A center. Physical Review. 1961; 121:1001-1014. ISSN: 0163-1829
Watkins G, Corbett JW. Physics Review. 1965; 138:A543. ISSN: 0163-1829
Jones B, Coomer BJ, Goss JP, Hourahine B, Recende A. The interaction of hydrogen with deep level defects in silicon. Solid State Phenomena. 1999; 71:173-248
Bakos T. Defects in amorphous SiO2: reactions, dynamics and optical properties [PhD Thesis]. Nashville; 2003. p.70
Milovzorov D. Electronic structure of nanocrystalline silicon and oxidized silicon surfaces. Electrochemical and Solid State Letters. 2001; 4(7):G61-G63. ISSN: 0013-4651
Biswas R, Li Q, Pan BC, Yoon Y. Mechanism for hydrogen diffusion in amorphous silicon. Physical Review B. 1998; 57:2253-2256. ISSN: 0163-1829
Milovzorov D. Memory cell with photoacoustic switching. In: Lai WY, Pau S, Daniel Lopez O, editors. Proceedings of the SPIE 5592. Nanofabrication: Technologies, Devices, and Applications. 2005. pp. 427–437. DOI: 10.1117/12.568290
Xu C, Taylor T, Burton G, Neumark D. The Journal of Chemical Physics. 1998; 108:7645. ISSN: 0021-9606
Voronkov M. Topics in Current Chemistry. 1986; 131:99. ISSN: 0340-1022
Lanoo M, Allan G. A cluster plus effective tight-binding study of SiOx systems. Solid State Communications. 1978; 28:733-739. ISSN: 0038-1098
Hausdorff F. Grundzuge der Mengenlehre. Vol. 184. Berlin; 1914. ISBN: 3-540-42224-2
Athreya B, Goldenfeld N, Dantzig J. Renormalization-group theory for the phase-field crystal equation. Physical Review E. 2006; 74:011601. ISSN: 0163-1829
Emmerich H, Granasy L, Lowen H. Selected issued of phase-field crystal simulations. The European Physics Journal Plus. 2011; 126:102. ISSN: 2190-5444