# Investigation of the Atomic Structure of Ge-Sb-Se Chalcogenide Glasses

Authors: M. Fabian; N. Dulgheru; K. Antonova; A. Szekeres; M. Gartner
Journal: Advances in Condensed Matter Physics
 (2018)

---

## Abstract

Glasses with composition ofGexSb40-xSe60 (x= 40, 35, 32, 27, 20, 15 at. %) have been synthesized. Neutron and X-ray diffraction techniques were used to study the atomic glassy structure, and Reverse Monte Carlo (RMC) simulations were applied to model the 3-dimensional atomic configurations and thorough mapping of the atomic parameters, such as first and second neighbour distances, coordination numbers, and bond-angle distributions. The results are explained with formation of GeSe4 and SbSe3 structural units, which correlate with the Ge/Sb ratio. For all the studied compositions, the Ge-Se, Sb-Se, Ge-Ge, and Se-Se bonds are significant. RMC simulations reveal the presence of Ge-Sb and Sb-Sb bonds, being dependent on Ge/Sb ratio. All atomic compositions satisfy formal valence requirements, i.e., Ge is fourfold coordinated, Sb is threefold coordinated, and Se is twofold coordinated. By increasing the Sb content, both the Se-Ge-Se bonds angle of 107±3° and Se-Sb-Se bonds angle of 118±3° decrease, respectively, indicating distortion of the structural units. Far infrared Fourier Transform spectroscopic measurements conducted in the range of 50-450 cm-1 at oblique (75°) incidence radiation have revealed clear dependences of the IR band’s shift and intensity on the glassy composition, showing features around x=27 at.% supporting the topological phase transition to a stable rigid network consisting mainly of SbSe3 pyramidal and GeSe4 tetrahedral clusters. These results are in agreement with the Reverse Monte Carlo models, which define the Ge and Sb environment.


---

## Body

## 1. Introduction

In recent years a great deal of interest has been devoted to the optical studies of chalcogenide glasses because of their wide area of applications as electrical and optical components as well as optical fibers in the infrared optical region, which are the basis of many applications [1–4]. Ternary Ge-Sb-Se glasses are promising materials for infrared optical fibers application because of their high transparency in the key region and good thermal, mechanical, and chemical properties. The physicochemical, elastic-plastic, and optical properties of the ternary Ge-Sb-Se chalcogenide system have been intensively studied [5–13], and the influence of Ge/Sb ratio on the optical, electronic, and microstructural properties of Ge-Sb-Se system has been considered. A common and useful parameter for the explanation of the compositional dependence on the physical properties is the coordination number Z of covalent bonds per atom, characterizing the structural atomic units [5, 13, 14]. For the compositions under investigation, a value of Z=2.65-2.67 corresponding to x=25-27 at.% has been obtained as critical for a transition from 2D cross-linked chains arrangements to formation of 3D network of tetrahedral GeSe4 and pyramidal SbSe3 when the x value increases [5, 13].Still one of the main questions remains how the structure of Ge environment changes, when the Ge atoms are replaced by Sb atoms. A comprehensive answer can be obtained by combining different methods and techniques, such as neutron diffraction (ND), high-energy X-ray diffraction (XRD), and far infrared (IR) Fourier Transform spectroscopy, which offers comprehensive information about the atomic structure. In addition, theoretical modeling of such experimental results using Reverse Monte Carlo (RMC) method allows building three-dimensional structure models and, thus, obtaining a more detailed description of the atomic-scale glassy structure. In connection with this, we have successfully performed experiments combining the ND and XRD techniques and RMC simulation method for the study of ternary As-Se-Te [15] and quaternary Ge-Sb-S(Se)-Te glasses [16–18]. Using the experience gained, we have started to investigate new Ge-Sb-Se system. The compositional dependence of the structure and physical and optical properties of ternary Ge-Sb-Se system is of interest, especially considering the effects of average coordination number,Z, of covalent bonds per atom.In this paper, we focused on the atomic-scale structure characterization of nonstoichiometricGexSb40-xSe60 compositions, wherex varied between 15 and 40 at. %. For examining the alteration of the glassy structure with changes in composition, we applied an approach in which ND and high-energy XRD techniques with Reverse Monte Carlo (RMC) modeling were combined. The procedure of RMC modeling is considered in details in Section 3.1. By this way, characteristic parameters, such as structure factors, partial atomic pair correlation functions, coordination numbers, and bond-angle distribution functions, are established. Recently, the chemical short order in Ge20SbxSe80-x (x=5, 15, 20) compositions has been studied by the nuclear techniques and extended X-ray absorption fine structure measurements [19] as one of the compositions, namely, Ge20Sb20Se40, is among the compositions we study herein. Our results obtained by the combination of different methods would widen the knowledge about these materials and would contribute to a better understanding of the structural changes by composition.In addition, we performed Fourier Transform infrared (FTIR) reflectance measurements in the terahertz spectral range in order to get more information about the microstructure and basic chemical bonding of the studied compositions. So far, few far IR spectral measurements have been carried out for studying the long-wavelength multiphonon edge of these materials. Raman and FTIR spectra have been reported in [5] and [8, 10, 20, 21], respectively, as in the latter PE pellets are used resulting in a not so good spectrum even when they are converted in absorption units. To avoid this, we carried out the FTIR measurements on thin (~ 1 μm) films deposited on quartz substrates. Using incidence radiation under a large angle of 75° the obtained spectra are more informative about the composition microstructure with variation ofxvalue. The results of this paper correspond to the conclusions presented by other authors and are in agreement with our earlier publications [5] on bulk samples served here as parent materials for thin films preparation.


## 2. Experimental Section

### 2.1. Sample Preparation

The binary Ge40Se60 and ternary GexSb40-xSe60 glasses with x=15, 20, 27, 32, and 35 at. % were prepared from elements of 99.999 % purity. Appropriate quantities from Ge, Se, and Sb were sealed into quartz ampoules after evacuation down to a pressure of 10-3 Pa. The bulk glasses were synthesized by the conventional melt-quenching method in a rotary furnace at 950°C. After homogenization for 24 h, the melts were quenched in air.Part of the bulk samples was powdered. At the ND and XRD studies, powdered material was used as specimen with ~ 3 g/each composition. The powdered glasses served also as parent material at evaporation of thin films for Fourier Transform infrared (FTIR) spectroscopic studies.The thin films were deposited onto quartz substrates by vacuum thermal evaporation of the powders from the corresponding parentGexSb40-xSe60 materials at a residual pressure of 10-4 Pa in the chamber. The thickness of the films (~1 μm) was controlled in situ by MIKI FFV quartz sensor device.


### 2.2. Neutron and X-Ray Diffraction Experiments

The neutron diffraction measurements were performed at room temperature in the 2-axis ‘PSD’ monochromatic neutron diffractometer (λ0=1.068 Å;Q=0.45-9.8 Å-1) [22] at the 10 MW Budapest Research Reactor, using thermal neutrons. The powder specimens (~3 g/each) were filled in cylindrical vanadium sample holder of 8 mm diameter. The structure factors,S(Q),were evaluated from the raw experimental data. The overall run of the ND experimental curves is similar for the investigated samples, but characteristic differences especially in the lowQ-range can be observed. Figure 1 presents theS(Q) from neutron diffraction experiments, where the results of the Reverse Monte Carlo model calculation are also drawn, but this will be discussed later, in Section 3.Figure 1
Neutron diffraction structure factor,S(Q), of the chalcogenide samples with composition of GexSb40-xSe60 (x=40, 35, 32, 27, 20, 15 at.%) is displayed: experimental data (colour marks) and RMC simulation (solid lines). For better clarity, the curves are shifted vertically.The X-ray diffraction experiments were carried out at the beam line BW5 at HASYLAB, DESY [23]. The powdered samples were filled into quartz capillary tubes of 2 mm in diameter (wall thickness of ~0.02 mm) and measured at room temperature. The energy of the radiation was 109.5 keV (λ0=0.113 Å). The high-energy synchrotron X-ray radiation makes it possible to reach diffraction data up to high-Q values. In this study the XRD structure factors were obtained up to 18 Å-1, as for higherQ-values the experimental data became rather noisy. Because the atomic parameters of the components are similar, the XRD structure factors were similar to ND ones. We consider here the XRD data for four samples which shows well the concentration dependence within the series. The XRD experimentalS(Q) data together with the RMC simulation results (details of the RMC modeling will be discussed in Section 3) are presented in Figure 2.Figure 2
X-ray diffraction structure factor,S(Q), of the ternary GexSb40-xSe60 with composition of x=35, 27, 15 at.% and the binary Ge40Se60 glasses is displayed: experimental data (colour marks) and RMC simulation (solid lines). For better clarity, the curves are shifted vertically.The shape and character of the ND (Figure1) and XRD (Figure 2) structure factors testify to a fully amorphous glassy structure in nature. However, there are some fine differences between the ND and XRD spectra. The differences in the overall run of these spectra lie in the different values of the weighting factors, wij, of the partial structure factors, Sij(Q), defined as(1)SQ=∑i,jkwijSijQ,(2)wijND=cicjbibj∑i,jkcibj2,(3)wijXRDQ=cicjfiQfjQ∑i,jkcifiQ2,where ci,cj are the molar fractions of the components, bi,bj are the coherent neutron and fi(Q),fj(Q) are the X-ray scattering amplitudes, andk is the number of elements in the sample. For the binary Ge40Se60 samplek is 2, thusk(k+1)/2=3, and for the ternary Ge-Sb-Se samples k is 3, thusk(k+1)/2=6, and, therefore, 6 different atom pairs are present in the studied series. The neutron scattering amplitude of an element is constant in the entireQ-range [24], while the X-ray scattering amplitude isQ-dependent [25] and for each atom it differs in a somewhat different way. To illustrate the differences, the comparison of the corresponding weighting factors, wij, for the two radiations is shown in Table 1, where wijND and wijXRD(Q) atQ=1.05 Å-1 are given for the Ge15Sb25Se60 sample.Table 1
Neutron and X-ray diffraction weighting factors,wij(%) at Q=1.05 Å-1, of the partial atomic pairs in Ge15Sb25Se60 sample.
G e 1 5 S b 2 5 S e 60 Ge-Ge Ge-Sb Ge-Se Sb-Sb Sb-Se Se-Se  ND weighting factor (%)  2.75 6.24 21.43 3.54 24.30 41.74  XRD weighting factor (%) atQ=1.05 Å -1  1.59 8.49 13.56 11.32 36.16 28.88It can be seen that the Ge-Se and Se-Se atom pairs have a significant contribution in the neutron experiment, while the Sb-Sb and Sb-Se atom pairs have a dominant weight in the X-ray experiment. On the other hand, the differences appeared for the Ge-Ge and Ge-Sb atom pairs are not more significant than those for the other ones. Taking into consideration all these characteristics, we have concluded that the two radiations give complementary information, and both types of measurements are needed to obtain a real structure for the investigated samples.


### 2.3. Fourier Transform Infrared Spectrophotometry

In order to detect the basic chemical bonding in theGexSb40-xSe60 compositions Fourier Transform infrared measurements were performed in the terahertz spectral range of 450-50 cm-1 of the polarized light. FTIR spectra were measured in close to reflectance-absorbance geometry with an incident radiation angle of 75°, permitting maximal optical path through the sample. A Bruker Vertex 70 instrument and a reflectance accessory Bruker A513/Q were used. The resolution was 2 cm-1 and the number of scans was 100. A gold standard mirror was used for the spectral normalization.


## 2.1. Sample Preparation

The binary Ge40Se60 and ternary GexSb40-xSe60 glasses with x=15, 20, 27, 32, and 35 at. % were prepared from elements of 99.999 % purity. Appropriate quantities from Ge, Se, and Sb were sealed into quartz ampoules after evacuation down to a pressure of 10-3 Pa. The bulk glasses were synthesized by the conventional melt-quenching method in a rotary furnace at 950°C. After homogenization for 24 h, the melts were quenched in air.Part of the bulk samples was powdered. At the ND and XRD studies, powdered material was used as specimen with ~ 3 g/each composition. The powdered glasses served also as parent material at evaporation of thin films for Fourier Transform infrared (FTIR) spectroscopic studies.The thin films were deposited onto quartz substrates by vacuum thermal evaporation of the powders from the corresponding parentGexSb40-xSe60 materials at a residual pressure of 10-4 Pa in the chamber. The thickness of the films (~1 μm) was controlled in situ by MIKI FFV quartz sensor device.


## 2.2. Neutron and X-Ray Diffraction Experiments

The neutron diffraction measurements were performed at room temperature in the 2-axis ‘PSD’ monochromatic neutron diffractometer (λ0=1.068 Å;Q=0.45-9.8 Å-1) [22] at the 10 MW Budapest Research Reactor, using thermal neutrons. The powder specimens (~3 g/each) were filled in cylindrical vanadium sample holder of 8 mm diameter. The structure factors,S(Q),were evaluated from the raw experimental data. The overall run of the ND experimental curves is similar for the investigated samples, but characteristic differences especially in the lowQ-range can be observed. Figure 1 presents theS(Q) from neutron diffraction experiments, where the results of the Reverse Monte Carlo model calculation are also drawn, but this will be discussed later, in Section 3.Figure 1
Neutron diffraction structure factor,S(Q), of the chalcogenide samples with composition of GexSb40-xSe60 (x=40, 35, 32, 27, 20, 15 at.%) is displayed: experimental data (colour marks) and RMC simulation (solid lines). For better clarity, the curves are shifted vertically.The X-ray diffraction experiments were carried out at the beam line BW5 at HASYLAB, DESY [23]. The powdered samples were filled into quartz capillary tubes of 2 mm in diameter (wall thickness of ~0.02 mm) and measured at room temperature. The energy of the radiation was 109.5 keV (λ0=0.113 Å). The high-energy synchrotron X-ray radiation makes it possible to reach diffraction data up to high-Q values. In this study the XRD structure factors were obtained up to 18 Å-1, as for higherQ-values the experimental data became rather noisy. Because the atomic parameters of the components are similar, the XRD structure factors were similar to ND ones. We consider here the XRD data for four samples which shows well the concentration dependence within the series. The XRD experimentalS(Q) data together with the RMC simulation results (details of the RMC modeling will be discussed in Section 3) are presented in Figure 2.Figure 2
X-ray diffraction structure factor,S(Q), of the ternary GexSb40-xSe60 with composition of x=35, 27, 15 at.% and the binary Ge40Se60 glasses is displayed: experimental data (colour marks) and RMC simulation (solid lines). For better clarity, the curves are shifted vertically.The shape and character of the ND (Figure1) and XRD (Figure 2) structure factors testify to a fully amorphous glassy structure in nature. However, there are some fine differences between the ND and XRD spectra. The differences in the overall run of these spectra lie in the different values of the weighting factors, wij, of the partial structure factors, Sij(Q), defined as(1)SQ=∑i,jkwijSijQ,(2)wijND=cicjbibj∑i,jkcibj2,(3)wijXRDQ=cicjfiQfjQ∑i,jkcifiQ2,where ci,cj are the molar fractions of the components, bi,bj are the coherent neutron and fi(Q),fj(Q) are the X-ray scattering amplitudes, andk is the number of elements in the sample. For the binary Ge40Se60 samplek is 2, thusk(k+1)/2=3, and for the ternary Ge-Sb-Se samples k is 3, thusk(k+1)/2=6, and, therefore, 6 different atom pairs are present in the studied series. The neutron scattering amplitude of an element is constant in the entireQ-range [24], while the X-ray scattering amplitude isQ-dependent [25] and for each atom it differs in a somewhat different way. To illustrate the differences, the comparison of the corresponding weighting factors, wij, for the two radiations is shown in Table 1, where wijND and wijXRD(Q) atQ=1.05 Å-1 are given for the Ge15Sb25Se60 sample.Table 1
Neutron and X-ray diffraction weighting factors,wij(%) at Q=1.05 Å-1, of the partial atomic pairs in Ge15Sb25Se60 sample.
G e 1 5 S b 2 5 S e 60 Ge-Ge Ge-Sb Ge-Se Sb-Sb Sb-Se Se-Se  ND weighting factor (%)  2.75 6.24 21.43 3.54 24.30 41.74  XRD weighting factor (%) atQ=1.05 Å -1  1.59 8.49 13.56 11.32 36.16 28.88It can be seen that the Ge-Se and Se-Se atom pairs have a significant contribution in the neutron experiment, while the Sb-Sb and Sb-Se atom pairs have a dominant weight in the X-ray experiment. On the other hand, the differences appeared for the Ge-Ge and Ge-Sb atom pairs are not more significant than those for the other ones. Taking into consideration all these characteristics, we have concluded that the two radiations give complementary information, and both types of measurements are needed to obtain a real structure for the investigated samples.


## 2.3. Fourier Transform Infrared Spectrophotometry

In order to detect the basic chemical bonding in theGexSb40-xSe60 compositions Fourier Transform infrared measurements were performed in the terahertz spectral range of 450-50 cm-1 of the polarized light. FTIR spectra were measured in close to reflectance-absorbance geometry with an incident radiation angle of 75°, permitting maximal optical path through the sample. A Bruker Vertex 70 instrument and a reflectance accessory Bruker A513/Q were used. The resolution was 2 cm-1 and the number of scans was 100. A gold standard mirror was used for the spectral normalization.


## 3. Results and Discussion

### 3.1. Reverse Monte Carlo Simulation

The experimental diffractionS(Q) data has been simulated by the RMC method [26], which is a widely used effective tool to model disordered structures. The RMC minimizes the squared difference between the experimentalS(Q) and the calculated one from a 3D atomic configuration. From the RMC configuration the partial pair correlation functions, (gij(r)), the coordination numbers, and the bond-angle distributions were calculated using the software package RMC++ [27].In order to investigate the possible effect of the starting configuration on the results, the following model was built up. The initial configuration was generated by random distribution of 10 000 atoms in a cubic simulation box. The experimentally measured density was 0.032, 0.030, 0.031, 0.038, 0.039, and 0.038 atoms/Å3 [28] corresponding to box edges of 31.90 Å, 34.66 Å, 34.29 Å, 32.05 Å, 31.76 Å, and 32.05 Å for the Ge40Se60, Ge35Sb5Se60, Ge32Sb8Se60, Ge27Sb13Se60, Ge20Sb20Se60, and Ge15Sb25Se60 samples, respectively.The Ge, Sb, and Se elements possess comparable neutron and X-ray scattering amplitudes, which makes it rather difficult to identify the atomic positions. In order to avoid unreasonable short-range orders for all atom pairs, in the RMC simulation procedure constraints were applied to the minimum interatomic distances between atom pairs (cut-off distances). The starting cut-off distances were taken from the previous studies for the similar contents [15, 17] and from the literature based on the results of similar or very close glassy compositions [29, 30]. Several RMC runs have been performed by modifying slightly the cut-off distances for each atom pairs. The final cut-off distances, used further in the RMC modeling, are summarized in Table 2. From our previous study and from literature data, we expect that the main glassy network is built up by the GeSe4/2 tetrahedra and the SbSe3/2 pyramid units. Therefore, it is suggested that the atomic structure of the ternary glasses is built up by chain connection of GeSe4/2 tetrahedra and SbSe3/2 pyramid through a 2-coordinated Se atoms. Thus, the 4-coordinated Ge unit connects through a Se atom to a 3-coordinated Sb unit and, therefore, the first neighbour distance of Ge-Sb is less realistic. In the literature, there are studies [29, 31–33], which support this suggestion. Considering the values of the weighting factors in Table 1, one can see that Ge-Sb bond’s weighting factors (wND,Ge-Sb=6.24% and wXRD,Ge-Sb=8.49%) are lower than those for Ge-Se (wND,Ge-Se=21.43% and wXRD,Ge-Se=13.56%) and Sb-Se (wND,Sb-Se=24.30% and wXRD,Sb-Se=36.16%), which bonds have well known first neighbour distances. From these follows that Ge-Sb bonds most probably have the second neighbour distance. Nevertheless, in the RMC modeling, we did not exclude any first neighbour distances and Figure 3 presents all partial pair correlation functions obtained from the RMC calculations.Table 2
Cut-off distances (Å) for atom pairs used in the final RMC run.
Atom pairs Cut-off distances (Å) G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60  Ge-Se 2.15 2.05 2.15 2.15 2.10 2.125  Sb-Se - 2.10 2.125 2.15 2.175 2.17  Ge-Sb - 2.25 2.25 2.20 2.25 2.25  Ge-Ge 2.30 2.20 2.20 2.20 2.20 2.25  Se-Se 2.00 2.20 2.20 2.15 2.15 2.175  Sb-Sb - 2.35 2.35 2.345 2.35 2.30Figure 3
Partial atom pair correlation functions,gij(r) for the Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle) compositions, obtained by RMC modeling of the Ge-Se (a), Sb-Se (b), Ge-Sb (c), Ge-Ge (d), Sb-Sb (e), and Se-Se (f) atom pairs. For better clarity, the curves are shifted vertically.
(a) (b) (c) (d) (e) (f)The experimental structure factors (S(Q)) are compared with those accumulated from the data of the RMC models of the corresponding glass compositions. For each sample, close to 20 RMC configurations were obtained in the intervals of more than 1 500 000 accepted moves of atoms inside of the simulation box. Good agreement between the experimental and calculated diffraction data was achieved. The comparison between the experimental and calculated structure factors are illustrated in Figures 1 and 2. The measured neutron diffraction data in Figure 1 shows significant dependence on the composition, especially at low-Q values. For all compositions, the first and the most intensive peaks were situated at 1.1 and 2.1 Å-1, while the next less intensive ones are at 3.55 Å-1 and at 5.5 Å-1. Figure 2 displays the XRDS(Q) spectra, where the first intensive peak is observed at 0.6 Å-1, a second intensive peak at 1.8 Å-1, and two other broaden and much less intensive peaks at 3.5 Å-1 and at 5.7 Å-1.In order to get comprehensive structural information about the studied series, we have investigated all homopolar and heteropolar atomic distances. The partial atomic pair correlation functions,gij(r), were revealed from the RMC simulation with a good reproducibility and acceptable statistics. The obtained partial pair distribution functions are presented in Figure 3. The first peaks are well defined and, thus, we could perform an unambiguous evaluation of distances.It was found that Ge-Se, Sb-Se, Ge-Ge, and Se-Se bonds are significant in the studied samples. The peaks related to the first interatomic distances appear at 2.35±0.01 Å, 2.54±0.01 Å, 2.47±0.03 Å, and 2.32±0.02 Å, respectively, as these peaks appear at the same position within the error for all the studied compositions. Because of the allowed cut-off distances (Table2), first interatomic distances between Ge-Sb and Sb-Sb atoms were found at 2.55±0.05 Å and 2.57±0.03 Å, but the corresponding first peak’s intensity was very low (Figures 3(c) and 3(e)). This suggests that second neighbour distance between Ge and Sb is more realistic. Second characteristic peaks were obtained for Ge-Sb, Ge-Ge, Sb-Sb, and Se-Se atom pairs, which are displayed at 3.80±0.05 Å, 3.85±0.05 Å, 3.80±0.05 Å, and 3.80±0.05 Å, respectively.The partial atomic pair correlation functions are summarized in Table3. These values, with exceptions of Ge-Sb and Ge-Ge, are close to the data given in [19], but the differences could be due to the annealing of samples in that work, which altered their glassy structure.Table 3
Interatomic distances,rij (Å), in chalcogenide samples obtained from the RMC simulation.
Atom pairs Interatomic distances, r i j (Å) G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60 Ge-Se 2.34/3.85±0.02 2.35/3.82±0.02 2.36/3.90±0.02 2.35/3.80±0.02 2.35/3.85±0.02 2.35/3.82±0.02  Sb-Se - 2.54/3.88±0.02 2.55/3.90±0.05 2.55/3.80±0.02 2.53/3.90±0.02 2.55/3.80±0.03  Ge-Sb - 2.60/3.80±0.05 2.60/3.55/3.80±0.05 2.50/3.80±0.05 2.50/3.85±0.05 2.52/3.80±0.05  Ge-Ge 2.50/3.80±0.02 2.50/3.80±0.02 2.47/3.90±0.03 2.45/3.80±0.02 2.45/3.90±0.03 2.50/3.80±0.02  Se-Se 2.30/3.85±0.02 2.30/3.80±0.05 2.30/3.60/3.90±0.05 2.30/3.75±0.05 2.30/3.83±0.05 2.33/3.80±0.05  Sb-Sb - 3.80/4.35±0.05 2.57/3.85±0.05 2.60/3.80±0.05 2.60/3.75±0.05 2.55/3.80±0.05The great advance of the RMC method is that the coordination number,CNij, can be obtained from the final particle configurations. From the partial pair distribution functions we calculated the number of nearest neighbours for Ge, Sb, and Se atoms using the corresponding bond cut-off distances. It is necessary to specify the range ofr within which atoms are counted as neighbours and which determine the coordination shells. We introduced the quantitiesr1 andr2, wherer1 andr2 are the positions of minimum values on the lower and upper side of the corresponding peaks, respectively. From the literature, we know some details related to the coordination numbers of the present atoms: concerning the Ge atoms, they are 4-coordinated, the Sb atoms are 3-coordinated, and Se atoms expected to be close to 2-coordinated sites. To obtain more information about the atomic structure, we carried out constrained simulations, in which Ge and Sb atoms were forced to have four and three neighbours, respectively. The corresponding average coordination numbers, CNij, obtained from the RMC simulations, are tabulated in Table 4. The coordination numbers show a clear dependence with the Ge/Sb ratio. The results suggest the presence of two basic structural units, namely, the GeSe4 tetrahedral units with fourfold coordinated Ge atom and SbSe3 pyramids with threefold coordinated Sb atom.Table 4
Average coordination numbers,CNij, calculated from the RMC simulation. In brackets are indicated the intervals of r within which the atoms are counted as neighbours and, thus, the actual coordination numbers are calculated. The error is ~5% for Ge-Se, Sb-Se, Se-Ge and ~10% for Ge-Sb, Sb-Ge, Se-Sb, Ge-Ge, Sb-Sb, Se-Se. The total coordination distributions are signed by Ge-X for Ge, Sb-X for Sb, and Se-X for Se atom.
Atom pairs Coordination Number, C N i j G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60 Ge-Se 2.97(r1:2.15-r2:2.90) 3.10(r1:2.2-r2:2.85) 3.22(r1:2.15-r2:2.8) 3.45(r1:2.15-r2:2.95) 3.63(r1:2.15-r2:2.95) 3.84(r1:2.1-r2:2.9)  Ge-Sb - 0.05(r1:2.4-r2:2.85) 0.02(r1:2.4-r2:2.8) 0.03(r1:2.2-r2:2.75) 0.05(r1:2.3-2.75) 0.07(r1:2.3-r2:2.7)  Sb-Se - 2.97(r1:2.2-r2:3.0) 2.80(r1:2.15-r2:3.1) 2.78(r1:2.15-r2:3.0) 2.64(r1:2.2-r2:2.95) 2.59(r1:2.15-r2:3.0)  Sb-Ge - 0.03(r1:2.4-r2:2.85) 0.05(r1:2.4-r2:2.8) 0.04(r1:2.2-r2:2.75) 0.09(r1:2.3-r2:2.75) 0.11(r1:2.3-r2:2.7)  Se-Ge 1.98(r1:2.15-r2:2.90) 1.81(r1:2.2-r2:2.85) 1.72(r1:2.15-r2:2.8) 1.55(r1:2.15-r2:2.95) 1.21(r1:2.2-r2:2.85) 0.96(r1:2.1-r2:2.9)  Se-Sb - 0.21(r1:2.2-3.0) 0.31(r1:2.15-3.1) 0.58(r1:2.15-3.0) 0.88(r1:2.15-2.95) 1.08(r1:2.15-3.0)  Ge-Ge 0.74(r1:2.25-r2:2.90) 0.55(r1:2.2-r2:2.95) 0.52(r1:2.15-r2:2.8) 0.40(r1:2.2-r2:2.8) 0.23(r1:2.2-r2:2.8) 0.05(r1:2.25-r2:2.75)  Sb-Sb - 0.00 0.03(r1:2.35-r2:2.75) 0.02(r1:2.45-r2:2.75) 0.3(r1:2.35-r2:2.85) 0.15(r1:2.2-r2:2.75)  Se-Se 0.04(r1:2.1-r2:2.6) 0.02(r1:2.1-r2:2.65) 0.02(r1:2.15-r2:2.6) 0.03(r1:2.15-r2:2.6) 0.05(r1:2.1-r2:2.65) 0.08(r1:2.1-r2:2.6)  Ge-X 3.71 3.69 3.76 3.88 3.91 3.96  Sb-X - 3.00 2.88 2.84 2.85 2.85  Se-X 2.02 2.02 2.05 2.16 2.14 2.12In Figure4 the average coordination number distributions of the possible chemical bonds, calculated from the RMC modeling, are presented. It can be seen that the average coordination number around Ge atoms is close to four atoms, as is expected by the formation of tetrahedral units in the network.Figure 4
Coordination number distributions,CNij, for the Ge-Se and Se-Ge bonds (a) and Sb-Se and Se-Sb bonds (b) in the corresponding compositions of Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle).
(a) (b)However, with increasing Sb content the Ge-Se coordination number,CNGe-Se, slightly increases from 3.10 to 3.84 atoms, which indicates a small degree of angle distortion affecting the tetrahedral surrounding. This may be caused by the formation of Se-Ge bonds, as the coordination of Se atoms increases from 2.02 to 2.16. The average CNGe-X coordination number increases from 3.69 to 3.96 with increasing Sb content. This suggests that the chalcogenide network consists of stable tetrahedral Ge units. The CNSb-Se coordination number continuously decreases from 2.97 to 2.50 with increasing Sb concentration. The Sb-X average coordination number slightly varied from 3.00 to 2.85. The Se-Ge coordination number continuously decreases from 1.98 to 0.96 with increasing Sb concentration, while the Se-Sb coordination number increases from 0.21 to 1.08; nevertheless, the concentration of Se atoms is kept constant.Figure5 shows the bond-angle distributions for the Ge and Sb atoms with Se-Ge-Se and Se-Sb-Se three-particle bond angles. We calculated the three-particle bond-angle distributions using the final atomic configuration of the RMC algorithm, plotted both as the function of cos⁡(Θ) (scale below) and Θ (upper scale), where Θ represents the actual bond angle.Figure 5
Bond-angle distributions of Se-Ge-Se (a) and Se-Sb-Se (b) configurations for the corresponding compositions of Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle).
(a) (b)For the Se-Ge-Se bonds, the peak positions are at 107±3° (for the ternary samples) [in agreement with [34]] and at 109.6±2° (for the binary sample) [35] which are very close to the tetrahedral angle of 109.47°. The peak distribution is very similar for all the studied samples. The broad distribution is quite asymmetric, but with the increase of Sb concentration, a shift down can be observed suggesting that the tetrahedral environment has become distorted. The distribution of Se-Sb-Se bonding angles is broad and asymmetric, the average angles being 118±3° and 111±5°. The broad distribution suggests that 3-fold Sb atoms are present. The Se-Sb-Se bond angles distribution also shows similar characteristics for all samples, as with the increase of Sb concentration a shift down to the 111±5° can be observed implying considerable distortion in Sb3 planar geometry. In most of the cases, a weak peak around 60±3° appears, strongly dependent on glassy composition, which is due to the occasional presence of threefold rings [35].The present RMC simulation of the ND and XRD data points out that the Sb atoms incorporated in Ge40Se60 glass bind covalently to Se and form SbSe3 - trigonal units. On the other hand, Se atoms, which are also connected to Ge, form GeSe4 tetrahedral structural units. These findings are in agreement with several spectroscopic studies [5, 8, 36].The atomic parameters support the formation of Ge-tetrahedral and Sb-trigonal units in Ge40Se60 glass structure. When the Sb content gradually increases, strong Sb units are formed with light influence of the tetrahedral Ge network, which causes only negligible changing of the partial pair distribution functions, coordination numbers, and bond-angle distributions. The Ge-Sb bonds show a broad distribution at 3.8 Å, which suggests the existence of connection between Ge-tetrahedral and Sb-pyramidal units through Se atoms. Nevertheless, the Ge-Ge bonds indicate connections between the Ge-tetrahedral, even if in the literature only theoretical bonding is reported. Based on these findings we suppose that a few numbers of Ge-clusters exist, which are connected to the Sb-pyramids. From the values of bond lengths and coordination numbers, we can conclude that the network structure is constructed through connection of tetrahedral GeSe4 and trigonal SbSe3 structural units. The average coordination numbers vary with glassy composition, as it follows from Figure 4. In all the studied samples, the concentration of the Se atoms is constant. With increasing of Sb concentration, coordination of Ge atoms becomes close to 4, while the coordination of Sb atoms becomes close to 3, both coordination sites being close to theoretical ones. These results are in agreement with the literature and support formation of tetrahedral Ge and pyramidal Sb units [19, 29, 37]. For the Se atoms, the coordination distribution changes in opposite direction, and with increasing of Sb concentration the Se-Ge coordination numbers decrease and the Se-Sb coordination numbers increase. This indicates that Se atom prefers to bind to Sb atom rather than Ge atom.From the above considerations we can conclude that the changes in the distribution of the chemical bonds from theGexSb40-xSe60 compositions with average coordination numbers Z < 2.67 to the compositions with Z > 2.67 are responsible for the structural rearrangement of the glassy network [5, 6, 13]


### 3.2. FTIR Spectra Analysis

The FTIR reflectance spectra in dependence on the indexx (Ge content) are presented in Figure 6. They show considerably rich set of bands because unlike the crystal structure, the lack of symmetry in amorphous material changes the selection rules and, thus, some forbidden IR modes become active. The spectra present well pronounced reflectance-absorbance bands of collective low-frequency vibrations. Although the spectral change with the Ge content is very small, there are bands which can be assigned to the vibrations affected by thex variation.Figure 6
Far FTIR reflectance spectra ofGexSb40-xSe60 films at different x values.The band at 250cm-1 is characteristic for GexSe100-x binary materials [12, 38]. With increasing Ge and with introducing Sb in our GexSb40-xSe60 samples the band becomes larger and more intensive. As we have already established that in GexSb40-xSe60 glasses structural transformation occurs at x= 25-27 at.% [6], we assume that GeSe2 and GeSe4 units coexist [5, 10, 13] together with Se-Se bond vibration (shoulder at 288 cm-1) [8]. With increasing the Ge content and, respectively, decreasing the content of heavier Sb atoms, these bands move to shorter wavelengths, as seen in Figure 7, where a peculiarity appears at x=27 at.%, corresponding to an average coordination number of Z=2.67.Figure 7
Shift of the bands at 250cm-1 (a) and 200 cm-1 (b) with variation of x value.
(a) (b)The broad and structured band at 160-210cm-1 is also sensitive to thex variation. Its components at 168, 187, and 203 cm-1 are reduced in intensity by increasing the Ge content. In [5, 39] the authors have discussed Raman spectra in this spectral region assigning the measured frequencies to the dynamics of the tetrahedral GeSe4 vibrations with addition of Sb atoms and to the appearance of pyramidal SbSe3 vibrations. Moreover, the heteropolar Sb-Se and Ge-Se bonds vibrational modes have close position at 190 cm-1 and 200 cm-1, respectively [5, 38]. The distinct peak around 168 cm-1 can be connected with formation of Ge-Ge bonds by cross-linked Ge atoms in the Ge2Se6/2 chains (Se3Ge-GeSe3 units) [5, 40, 41]. We are convinced that these modes are responsible for forming the complex band as IR active due to the lack of symmetry in our amorphous samples, as it was obtained by FTIR transmission measurements in [10, 21].The band picked at 150cm-1 with a shoulder at 157 cm-1 does not shift withx but its intensity varies (see Figure 8). It could be related to the Sb-Se (at smallx) and to the Ge-Se (at higherx) nonsymmetric vibrations, respectively, [8, 21]. At the highest value ofx=35 at.%, the band is slightly enlarged thus showing the coexistence of these vibrations.Figure 8
Dependence of bands intensity on the x values.The main bands in the terahertz range, 106cm-1 and 74 cm-1, belong to the vibrations of GeSe4 tetrahedral and respectively to SbSe3 pyramidal units [8, 21]. Because of these complex structures, they are surrounded with many weak peaks and shoulders. In Figure 8, it can be seen that the band intensities as function of x evolve in opposite directions, following the variation of the Ge and Sb content, respectively.Bands in the FTIR spectra, located in the 250-320cm-1 region (Figure 6) have very low intensity and a tendency to vanish by increasing the Ge at.% and decreasing the Sb at.%, respectively. They may be attributed to several vibrations of Se-Se, Ge-Ge, and Ge-Sb bonds in modified GeSe4 - structural units [21, 38, 41].In the 450-350cm-1 spectral range, a very broad band, centred at 420 cm-1 and with a weak shoulder around 340 cm-1, is observed. It is independent on the glassy composition and is related to vibration of oxygen impurity atoms in Ge-O bonding [42, 43]. The presence of oxygen-related impurity bonds in these Ge-Sb-Se films has been detected by infrared ellipsometric measurements performed in the 400- 4000 cm-1 spectral range [44].From the study of the basic chemical bonds vibrational modes in chalcogenideGexSb40-xSe60 films considering the FTIR reflection spectra at oblique (75°) incidence radiation, we obtained specific dependences of the band’s shift and intensity on the film composition. Figures 7 and 8 show features around x=27 at.%, thus confirming the earlier published model of the topological phase transition to a stable rigid network consisting mainly of SbSe3 pyramidal and GeSe4 tetrahedral clusters. These results support the RMC modeling, which reveals that through connection of these GeSe4 and SbSe3 structural units, it is constructed a Ge-Sb-Se network structure.


## 3.1. Reverse Monte Carlo Simulation

The experimental diffractionS(Q) data has been simulated by the RMC method [26], which is a widely used effective tool to model disordered structures. The RMC minimizes the squared difference between the experimentalS(Q) and the calculated one from a 3D atomic configuration. From the RMC configuration the partial pair correlation functions, (gij(r)), the coordination numbers, and the bond-angle distributions were calculated using the software package RMC++ [27].In order to investigate the possible effect of the starting configuration on the results, the following model was built up. The initial configuration was generated by random distribution of 10 000 atoms in a cubic simulation box. The experimentally measured density was 0.032, 0.030, 0.031, 0.038, 0.039, and 0.038 atoms/Å3 [28] corresponding to box edges of 31.90 Å, 34.66 Å, 34.29 Å, 32.05 Å, 31.76 Å, and 32.05 Å for the Ge40Se60, Ge35Sb5Se60, Ge32Sb8Se60, Ge27Sb13Se60, Ge20Sb20Se60, and Ge15Sb25Se60 samples, respectively.The Ge, Sb, and Se elements possess comparable neutron and X-ray scattering amplitudes, which makes it rather difficult to identify the atomic positions. In order to avoid unreasonable short-range orders for all atom pairs, in the RMC simulation procedure constraints were applied to the minimum interatomic distances between atom pairs (cut-off distances). The starting cut-off distances were taken from the previous studies for the similar contents [15, 17] and from the literature based on the results of similar or very close glassy compositions [29, 30]. Several RMC runs have been performed by modifying slightly the cut-off distances for each atom pairs. The final cut-off distances, used further in the RMC modeling, are summarized in Table 2. From our previous study and from literature data, we expect that the main glassy network is built up by the GeSe4/2 tetrahedra and the SbSe3/2 pyramid units. Therefore, it is suggested that the atomic structure of the ternary glasses is built up by chain connection of GeSe4/2 tetrahedra and SbSe3/2 pyramid through a 2-coordinated Se atoms. Thus, the 4-coordinated Ge unit connects through a Se atom to a 3-coordinated Sb unit and, therefore, the first neighbour distance of Ge-Sb is less realistic. In the literature, there are studies [29, 31–33], which support this suggestion. Considering the values of the weighting factors in Table 1, one can see that Ge-Sb bond’s weighting factors (wND,Ge-Sb=6.24% and wXRD,Ge-Sb=8.49%) are lower than those for Ge-Se (wND,Ge-Se=21.43% and wXRD,Ge-Se=13.56%) and Sb-Se (wND,Sb-Se=24.30% and wXRD,Sb-Se=36.16%), which bonds have well known first neighbour distances. From these follows that Ge-Sb bonds most probably have the second neighbour distance. Nevertheless, in the RMC modeling, we did not exclude any first neighbour distances and Figure 3 presents all partial pair correlation functions obtained from the RMC calculations.Table 2
Cut-off distances (Å) for atom pairs used in the final RMC run.
Atom pairs Cut-off distances (Å) G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60  Ge-Se 2.15 2.05 2.15 2.15 2.10 2.125  Sb-Se - 2.10 2.125 2.15 2.175 2.17  Ge-Sb - 2.25 2.25 2.20 2.25 2.25  Ge-Ge 2.30 2.20 2.20 2.20 2.20 2.25  Se-Se 2.00 2.20 2.20 2.15 2.15 2.175  Sb-Sb - 2.35 2.35 2.345 2.35 2.30Figure 3
Partial atom pair correlation functions,gij(r) for the Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle) compositions, obtained by RMC modeling of the Ge-Se (a), Sb-Se (b), Ge-Sb (c), Ge-Ge (d), Sb-Sb (e), and Se-Se (f) atom pairs. For better clarity, the curves are shifted vertically.
(a) (b) (c) (d) (e) (f)The experimental structure factors (S(Q)) are compared with those accumulated from the data of the RMC models of the corresponding glass compositions. For each sample, close to 20 RMC configurations were obtained in the intervals of more than 1 500 000 accepted moves of atoms inside of the simulation box. Good agreement between the experimental and calculated diffraction data was achieved. The comparison between the experimental and calculated structure factors are illustrated in Figures 1 and 2. The measured neutron diffraction data in Figure 1 shows significant dependence on the composition, especially at low-Q values. For all compositions, the first and the most intensive peaks were situated at 1.1 and 2.1 Å-1, while the next less intensive ones are at 3.55 Å-1 and at 5.5 Å-1. Figure 2 displays the XRDS(Q) spectra, where the first intensive peak is observed at 0.6 Å-1, a second intensive peak at 1.8 Å-1, and two other broaden and much less intensive peaks at 3.5 Å-1 and at 5.7 Å-1.In order to get comprehensive structural information about the studied series, we have investigated all homopolar and heteropolar atomic distances. The partial atomic pair correlation functions,gij(r), were revealed from the RMC simulation with a good reproducibility and acceptable statistics. The obtained partial pair distribution functions are presented in Figure 3. The first peaks are well defined and, thus, we could perform an unambiguous evaluation of distances.It was found that Ge-Se, Sb-Se, Ge-Ge, and Se-Se bonds are significant in the studied samples. The peaks related to the first interatomic distances appear at 2.35±0.01 Å, 2.54±0.01 Å, 2.47±0.03 Å, and 2.32±0.02 Å, respectively, as these peaks appear at the same position within the error for all the studied compositions. Because of the allowed cut-off distances (Table2), first interatomic distances between Ge-Sb and Sb-Sb atoms were found at 2.55±0.05 Å and 2.57±0.03 Å, but the corresponding first peak’s intensity was very low (Figures 3(c) and 3(e)). This suggests that second neighbour distance between Ge and Sb is more realistic. Second characteristic peaks were obtained for Ge-Sb, Ge-Ge, Sb-Sb, and Se-Se atom pairs, which are displayed at 3.80±0.05 Å, 3.85±0.05 Å, 3.80±0.05 Å, and 3.80±0.05 Å, respectively.The partial atomic pair correlation functions are summarized in Table3. These values, with exceptions of Ge-Sb and Ge-Ge, are close to the data given in [19], but the differences could be due to the annealing of samples in that work, which altered their glassy structure.Table 3
Interatomic distances,rij (Å), in chalcogenide samples obtained from the RMC simulation.
Atom pairs Interatomic distances, r i j (Å) G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60 Ge-Se 2.34/3.85±0.02 2.35/3.82±0.02 2.36/3.90±0.02 2.35/3.80±0.02 2.35/3.85±0.02 2.35/3.82±0.02  Sb-Se - 2.54/3.88±0.02 2.55/3.90±0.05 2.55/3.80±0.02 2.53/3.90±0.02 2.55/3.80±0.03  Ge-Sb - 2.60/3.80±0.05 2.60/3.55/3.80±0.05 2.50/3.80±0.05 2.50/3.85±0.05 2.52/3.80±0.05  Ge-Ge 2.50/3.80±0.02 2.50/3.80±0.02 2.47/3.90±0.03 2.45/3.80±0.02 2.45/3.90±0.03 2.50/3.80±0.02  Se-Se 2.30/3.85±0.02 2.30/3.80±0.05 2.30/3.60/3.90±0.05 2.30/3.75±0.05 2.30/3.83±0.05 2.33/3.80±0.05  Sb-Sb - 3.80/4.35±0.05 2.57/3.85±0.05 2.60/3.80±0.05 2.60/3.75±0.05 2.55/3.80±0.05The great advance of the RMC method is that the coordination number,CNij, can be obtained from the final particle configurations. From the partial pair distribution functions we calculated the number of nearest neighbours for Ge, Sb, and Se atoms using the corresponding bond cut-off distances. It is necessary to specify the range ofr within which atoms are counted as neighbours and which determine the coordination shells. We introduced the quantitiesr1 andr2, wherer1 andr2 are the positions of minimum values on the lower and upper side of the corresponding peaks, respectively. From the literature, we know some details related to the coordination numbers of the present atoms: concerning the Ge atoms, they are 4-coordinated, the Sb atoms are 3-coordinated, and Se atoms expected to be close to 2-coordinated sites. To obtain more information about the atomic structure, we carried out constrained simulations, in which Ge and Sb atoms were forced to have four and three neighbours, respectively. The corresponding average coordination numbers, CNij, obtained from the RMC simulations, are tabulated in Table 4. The coordination numbers show a clear dependence with the Ge/Sb ratio. The results suggest the presence of two basic structural units, namely, the GeSe4 tetrahedral units with fourfold coordinated Ge atom and SbSe3 pyramids with threefold coordinated Sb atom.Table 4
Average coordination numbers,CNij, calculated from the RMC simulation. In brackets are indicated the intervals of r within which the atoms are counted as neighbours and, thus, the actual coordination numbers are calculated. The error is ~5% for Ge-Se, Sb-Se, Se-Ge and ~10% for Ge-Sb, Sb-Ge, Se-Sb, Ge-Ge, Sb-Sb, Se-Se. The total coordination distributions are signed by Ge-X for Ge, Sb-X for Sb, and Se-X for Se atom.
Atom pairs Coordination Number, C N i j G e 40 S e 60 G e 35 S b 5 S e 60 G e 32 S b 8 S e 60 G e 27 S b 13 S e 60 G e 20 S b 20 S e 60 G e 15 S b 25 S e 60 Ge-Se 2.97(r1:2.15-r2:2.90) 3.10(r1:2.2-r2:2.85) 3.22(r1:2.15-r2:2.8) 3.45(r1:2.15-r2:2.95) 3.63(r1:2.15-r2:2.95) 3.84(r1:2.1-r2:2.9)  Ge-Sb - 0.05(r1:2.4-r2:2.85) 0.02(r1:2.4-r2:2.8) 0.03(r1:2.2-r2:2.75) 0.05(r1:2.3-2.75) 0.07(r1:2.3-r2:2.7)  Sb-Se - 2.97(r1:2.2-r2:3.0) 2.80(r1:2.15-r2:3.1) 2.78(r1:2.15-r2:3.0) 2.64(r1:2.2-r2:2.95) 2.59(r1:2.15-r2:3.0)  Sb-Ge - 0.03(r1:2.4-r2:2.85) 0.05(r1:2.4-r2:2.8) 0.04(r1:2.2-r2:2.75) 0.09(r1:2.3-r2:2.75) 0.11(r1:2.3-r2:2.7)  Se-Ge 1.98(r1:2.15-r2:2.90) 1.81(r1:2.2-r2:2.85) 1.72(r1:2.15-r2:2.8) 1.55(r1:2.15-r2:2.95) 1.21(r1:2.2-r2:2.85) 0.96(r1:2.1-r2:2.9)  Se-Sb - 0.21(r1:2.2-3.0) 0.31(r1:2.15-3.1) 0.58(r1:2.15-3.0) 0.88(r1:2.15-2.95) 1.08(r1:2.15-3.0)  Ge-Ge 0.74(r1:2.25-r2:2.90) 0.55(r1:2.2-r2:2.95) 0.52(r1:2.15-r2:2.8) 0.40(r1:2.2-r2:2.8) 0.23(r1:2.2-r2:2.8) 0.05(r1:2.25-r2:2.75)  Sb-Sb - 0.00 0.03(r1:2.35-r2:2.75) 0.02(r1:2.45-r2:2.75) 0.3(r1:2.35-r2:2.85) 0.15(r1:2.2-r2:2.75)  Se-Se 0.04(r1:2.1-r2:2.6) 0.02(r1:2.1-r2:2.65) 0.02(r1:2.15-r2:2.6) 0.03(r1:2.15-r2:2.6) 0.05(r1:2.1-r2:2.65) 0.08(r1:2.1-r2:2.6)  Ge-X 3.71 3.69 3.76 3.88 3.91 3.96  Sb-X - 3.00 2.88 2.84 2.85 2.85  Se-X 2.02 2.02 2.05 2.16 2.14 2.12In Figure4 the average coordination number distributions of the possible chemical bonds, calculated from the RMC modeling, are presented. It can be seen that the average coordination number around Ge atoms is close to four atoms, as is expected by the formation of tetrahedral units in the network.Figure 4
Coordination number distributions,CNij, for the Ge-Se and Se-Ge bonds (a) and Sb-Se and Se-Sb bonds (b) in the corresponding compositions of Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle).
(a) (b)However, with increasing Sb content the Ge-Se coordination number,CNGe-Se, slightly increases from 3.10 to 3.84 atoms, which indicates a small degree of angle distortion affecting the tetrahedral surrounding. This may be caused by the formation of Se-Ge bonds, as the coordination of Se atoms increases from 2.02 to 2.16. The average CNGe-X coordination number increases from 3.69 to 3.96 with increasing Sb content. This suggests that the chalcogenide network consists of stable tetrahedral Ge units. The CNSb-Se coordination number continuously decreases from 2.97 to 2.50 with increasing Sb concentration. The Sb-X average coordination number slightly varied from 3.00 to 2.85. The Se-Ge coordination number continuously decreases from 1.98 to 0.96 with increasing Sb concentration, while the Se-Sb coordination number increases from 0.21 to 1.08; nevertheless, the concentration of Se atoms is kept constant.Figure5 shows the bond-angle distributions for the Ge and Sb atoms with Se-Ge-Se and Se-Sb-Se three-particle bond angles. We calculated the three-particle bond-angle distributions using the final atomic configuration of the RMC algorithm, plotted both as the function of cos⁡(Θ) (scale below) and Θ (upper scale), where Θ represents the actual bond angle.Figure 5
Bond-angle distributions of Se-Ge-Se (a) and Se-Sb-Se (b) configurations for the corresponding compositions of Ge40Se60 (black square), Ge35Sb5Se60 (red circle), Ge32Sb8Se60 (blue triangle), Ge27Sb13Se60 (green crosses), Ge20Sb20Se60 (orange square), and Ge15Sb25Se60 (purple reverse triangle).
(a) (b)For the Se-Ge-Se bonds, the peak positions are at 107±3° (for the ternary samples) [in agreement with [34]] and at 109.6±2° (for the binary sample) [35] which are very close to the tetrahedral angle of 109.47°. The peak distribution is very similar for all the studied samples. The broad distribution is quite asymmetric, but with the increase of Sb concentration, a shift down can be observed suggesting that the tetrahedral environment has become distorted. The distribution of Se-Sb-Se bonding angles is broad and asymmetric, the average angles being 118±3° and 111±5°. The broad distribution suggests that 3-fold Sb atoms are present. The Se-Sb-Se bond angles distribution also shows similar characteristics for all samples, as with the increase of Sb concentration a shift down to the 111±5° can be observed implying considerable distortion in Sb3 planar geometry. In most of the cases, a weak peak around 60±3° appears, strongly dependent on glassy composition, which is due to the occasional presence of threefold rings [35].The present RMC simulation of the ND and XRD data points out that the Sb atoms incorporated in Ge40Se60 glass bind covalently to Se and form SbSe3 - trigonal units. On the other hand, Se atoms, which are also connected to Ge, form GeSe4 tetrahedral structural units. These findings are in agreement with several spectroscopic studies [5, 8, 36].The atomic parameters support the formation of Ge-tetrahedral and Sb-trigonal units in Ge40Se60 glass structure. When the Sb content gradually increases, strong Sb units are formed with light influence of the tetrahedral Ge network, which causes only negligible changing of the partial pair distribution functions, coordination numbers, and bond-angle distributions. The Ge-Sb bonds show a broad distribution at 3.8 Å, which suggests the existence of connection between Ge-tetrahedral and Sb-pyramidal units through Se atoms. Nevertheless, the Ge-Ge bonds indicate connections between the Ge-tetrahedral, even if in the literature only theoretical bonding is reported. Based on these findings we suppose that a few numbers of Ge-clusters exist, which are connected to the Sb-pyramids. From the values of bond lengths and coordination numbers, we can conclude that the network structure is constructed through connection of tetrahedral GeSe4 and trigonal SbSe3 structural units. The average coordination numbers vary with glassy composition, as it follows from Figure 4. In all the studied samples, the concentration of the Se atoms is constant. With increasing of Sb concentration, coordination of Ge atoms becomes close to 4, while the coordination of Sb atoms becomes close to 3, both coordination sites being close to theoretical ones. These results are in agreement with the literature and support formation of tetrahedral Ge and pyramidal Sb units [19, 29, 37]. For the Se atoms, the coordination distribution changes in opposite direction, and with increasing of Sb concentration the Se-Ge coordination numbers decrease and the Se-Sb coordination numbers increase. This indicates that Se atom prefers to bind to Sb atom rather than Ge atom.From the above considerations we can conclude that the changes in the distribution of the chemical bonds from theGexSb40-xSe60 compositions with average coordination numbers Z < 2.67 to the compositions with Z > 2.67 are responsible for the structural rearrangement of the glassy network [5, 6, 13]


## 3.2. FTIR Spectra Analysis

The FTIR reflectance spectra in dependence on the indexx (Ge content) are presented in Figure 6. They show considerably rich set of bands because unlike the crystal structure, the lack of symmetry in amorphous material changes the selection rules and, thus, some forbidden IR modes become active. The spectra present well pronounced reflectance-absorbance bands of collective low-frequency vibrations. Although the spectral change with the Ge content is very small, there are bands which can be assigned to the vibrations affected by thex variation.Figure 6
Far FTIR reflectance spectra ofGexSb40-xSe60 films at different x values.The band at 250cm-1 is characteristic for GexSe100-x binary materials [12, 38]. With increasing Ge and with introducing Sb in our GexSb40-xSe60 samples the band becomes larger and more intensive. As we have already established that in GexSb40-xSe60 glasses structural transformation occurs at x= 25-27 at.% [6], we assume that GeSe2 and GeSe4 units coexist [5, 10, 13] together with Se-Se bond vibration (shoulder at 288 cm-1) [8]. With increasing the Ge content and, respectively, decreasing the content of heavier Sb atoms, these bands move to shorter wavelengths, as seen in Figure 7, where a peculiarity appears at x=27 at.%, corresponding to an average coordination number of Z=2.67.Figure 7
Shift of the bands at 250cm-1 (a) and 200 cm-1 (b) with variation of x value.
(a) (b)The broad and structured band at 160-210cm-1 is also sensitive to thex variation. Its components at 168, 187, and 203 cm-1 are reduced in intensity by increasing the Ge content. In [5, 39] the authors have discussed Raman spectra in this spectral region assigning the measured frequencies to the dynamics of the tetrahedral GeSe4 vibrations with addition of Sb atoms and to the appearance of pyramidal SbSe3 vibrations. Moreover, the heteropolar Sb-Se and Ge-Se bonds vibrational modes have close position at 190 cm-1 and 200 cm-1, respectively [5, 38]. The distinct peak around 168 cm-1 can be connected with formation of Ge-Ge bonds by cross-linked Ge atoms in the Ge2Se6/2 chains (Se3Ge-GeSe3 units) [5, 40, 41]. We are convinced that these modes are responsible for forming the complex band as IR active due to the lack of symmetry in our amorphous samples, as it was obtained by FTIR transmission measurements in [10, 21].The band picked at 150cm-1 with a shoulder at 157 cm-1 does not shift withx but its intensity varies (see Figure 8). It could be related to the Sb-Se (at smallx) and to the Ge-Se (at higherx) nonsymmetric vibrations, respectively, [8, 21]. At the highest value ofx=35 at.%, the band is slightly enlarged thus showing the coexistence of these vibrations.Figure 8
Dependence of bands intensity on the x values.The main bands in the terahertz range, 106cm-1 and 74 cm-1, belong to the vibrations of GeSe4 tetrahedral and respectively to SbSe3 pyramidal units [8, 21]. Because of these complex structures, they are surrounded with many weak peaks and shoulders. In Figure 8, it can be seen that the band intensities as function of x evolve in opposite directions, following the variation of the Ge and Sb content, respectively.Bands in the FTIR spectra, located in the 250-320cm-1 region (Figure 6) have very low intensity and a tendency to vanish by increasing the Ge at.% and decreasing the Sb at.%, respectively. They may be attributed to several vibrations of Se-Se, Ge-Ge, and Ge-Sb bonds in modified GeSe4 - structural units [21, 38, 41].In the 450-350cm-1 spectral range, a very broad band, centred at 420 cm-1 and with a weak shoulder around 340 cm-1, is observed. It is independent on the glassy composition and is related to vibration of oxygen impurity atoms in Ge-O bonding [42, 43]. The presence of oxygen-related impurity bonds in these Ge-Sb-Se films has been detected by infrared ellipsometric measurements performed in the 400- 4000 cm-1 spectral range [44].From the study of the basic chemical bonds vibrational modes in chalcogenideGexSb40-xSe60 films considering the FTIR reflection spectra at oblique (75°) incidence radiation, we obtained specific dependences of the band’s shift and intensity on the film composition. Figures 7 and 8 show features around x=27 at.%, thus confirming the earlier published model of the topological phase transition to a stable rigid network consisting mainly of SbSe3 pyramidal and GeSe4 tetrahedral clusters. These results support the RMC modeling, which reveals that through connection of these GeSe4 and SbSe3 structural units, it is constructed a Ge-Sb-Se network structure.


## 4. Conclusions

Binary Ge40Se60 and ternary GexSb40-xSe60 glasses with x=35, 32, 27, 20, and 15 at.% compositions have been synthesized. Their structural properties have been investigated by neutron and X-ray diffraction measurements and RMC simulation procedures, in which all the possible chemical bonds in the glasses are taken into consideration. All specimens are fully amorphous in nature as confirmed by the ND and XRD patterns. It is revealed that addition of Sb atom to the binary Ge-Se glass does not change significantly the basic Ge-tetrahedral structural units but leads to their distortion. With increasing of Sb concentration, well-defined Sb-pyramidal units are formed in the ternary glass structure. The GexSb40-xSe60 glassy network builds up from GeSe4/2 tetrahedra and the SbSe3/2 pyramid units, connected through 2-coordinated Se atoms. Therefore, we may conclude that the “structural unit”-based description of the RMC results is one possible way of interpreting the measured diffraction data.The FTIR reflection spectra of chalcogenideGexSb40-xSe60 films studied in the terahertz spectral range at oblique (75°) incidence radiation showed characteristic bands related to Sb-Se and Ge-Se chemical bonds vibration in SbSe3 pyramidal and GeSe4 tetrahedral units, respectively. The peculiarities in the compositional dependences of the band positions and intensities confirm the existence of a topological phase transition in the GexSb40-xSe60 glasses around a critical Ge content of x=27 at.%., corresponding to an average coordination number Z=2.67.



---

*Source: 7158079-2018-12-02.xml*