# Ge-Sb-S-Se-Te amorphous chalcogenide thin films towards on-chip nonlinear photonic devices

J.-B. Dory, C. Castro-Chavarria, A. Verdy, J.-B. Jager, M. Bernard, C.

Sabbione, M. Tessaire, J.-M. Fedeli, A. Coillet, B. Cluzel, et al.

► To cite this version:

J.-B. Dory, C. Castro-Chavarria, A. Verdy, J.-B. Jager, M. Bernard, et al., Ge-Sb-S-Se-Te amorphous chalcogenide thin films towards on-chip nonlinear photonic devices. Scientific Reports, 2020, 10, pp.11894. 10.1038/s41598-020-67377-9. cea-02903381

HAL Id: cea-02903381

https://cea.hal.science/cea-02903381

Submitted on 3 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Ge-Sb-S-Se-Te amorphous chalcogenide thin films towards on-chip nonlinear photonic devices

J.-B. Dory<sup>D1</sup>, C. Castro-Chavarria<sup>1</sup>, A. Verdy<sup>D1</sup>, J.-B. Jager<sup>2</sup>, M. Bernard<sup>1</sup>, C. Sabbione<sup>1</sup>, M. Tessaire<sup>1</sup>, J.-M. Fédéli<sup>1</sup>, A. Coillet<sup>D3</sup>, B. Cluzel<sup>3</sup> & P. Noé<sup>D1✉</sup>

Thanks to their unique optical properties Ge-Sb-S-Se-Te amorphous chalcogenide materials and compounds offer tremendous opportunities of applications, in particular in near and mid-infrared range. This spectral range is for instance of high interest for photonics or optical sensors. Using co-sputtering technique of chalcogenide compound targets in a 200 mm industrial deposition tool, we show how by modifying the amorphous structure of GeSb<sub>w</sub>S<sub>x</sub>Se<sub>y</sub>Te<sub>z</sub> chalcogenide thin films one can significantly tailor their linear and nonlinear optical properties. Modelling of spectroscopic ellipsometry data collected on the as-deposited chalcogenide thin films is used to evaluate their linear and nonlinear properties. Moreover, Raman and Fourier-transform infrared spectroscopies permitted to get a description of their amorphous structure. For the purpose of applications, their thermal stability upon annealing is also evaluated. We demonstrate that depending on the GeSb<sub>w</sub>S<sub>x</sub>Se<sub>y</sub>Te<sub>z</sub> film composition a trade-off between a high transparency in near- or mid-infrared ranges, strong nonlinearity and good thermal stability can be found in order to use such materials for applications compatible with the standard CMOS integration processes of microelectronics and photonics.

Chalcogenides are commonly defined as non-oxide compounds containing at least one chalcogen element such as S, Se and/or Te (belonging to group 16 of O) alloyed with electropositive elements (more often elements of group 15 (As, Sb, Bi) and/or group 14 (Si, Ge, Sn, Pb)). Chalcogenide exhibit a unique portfolio of properties which has led to their wide use for non-volatile memory applications such as optical data storage (CD-RW and DVD-RAM), Conductive-Bridging Random Access Memory or Phase-Change Random Access Memory<fnref n="1" />. More recently, thanks to the huge electronic nonlinearities and discontinuity [in particular the Ovonic Threshold Switching (OTS) mechanism] observed in some chalcogenide glasses (CGs) under electrical field application, the latter are considered as the most promising materials to be used as innovative selector element in 3D memory arrays²,³. Besides, thanks to a high transparency window in the infrared range and large optical nonlinearities⁴,⁵, chalcogenide glasses offer also opportunities for elaboration of innovative mid-infrared (MIR) components such as MIR super-continuum (SC) laser sources, optical sensors, IR micro-lens arrays and all-optical integrated circuits⁶,⁷. The potential of amorphous semiconductors for optical applications is already demonstrated but the progresses in glass science is still behind crystal science due to inherent complexity of highly disordered systems hindering thus their structural description and theoretical modelling⁸. Up to now, state-of-the-art of MIR SC generation have been achieved by using mainly chalcogenide compounds containing Arsenic such as As₂S₃ and As₂Se₃ fibres⁹,¹⁰ or GeAsSe rib waveguides¹¹. However, the R.E.A.C.H. European recommendation (https://echa.europa.eu/regulations/reach) as well as a recent publication from the World Health Organization (https://www.who.int/ipcs/assessment/public_health/chemicals_phc/en/) have both identified Arsenic as one of the ten most harmful chemicals for human health. As a result, developing new materials using less and less harmful or rare element is a huge challenge for the materials science.

In that context, the aim of the present study is to develop As-free amorphous chalcogenide thin films compatible with CMOS technology of microelectronics and photonics. The ultimate purpose would be to give clear clues to select the most suitable chalcogenide compositions exhibiting the best trade-off between stability and optical properties in order to enable future achievement of passive and active on-chip MIR components such

<footnote n="1">¹Université Grenoble Alpes, CEA, LETI, MINATEC Campus, 17 Avenue des Martyrs, 38000 Grenoble, France. ²Université Grenoble Alpes, CEA, IRIG, MINATEC Campus, 17 Avenue des Martyrs, 38000 Grenoble, France. ³ICB, UMR CNRS 6303, Université de Bourgogne Franche Comté, 9, Avenue Alain-Savary, BP 47870, 21078 Dijon cedex, France. ✉email: pierre.noe@cea.fr</footnote>
as MIR optical waveguides or portable SC laser sources. To this aim, a wide composition range of amorphous GeSb<sub>w</sub>S<sub>x</sub>Se<sub>y</sub>Te<sub>z</sub> chalcogenide thin films were deposited on 200 nm Si-based substrates by means of magnetron co-sputtering technique. Despite the ability of such technique to permit a fast and easy study of a wide range of chalcogenide thin films' compositions, only very few studies of the optical properties of co-sputtered chalcogenide films have been reported in literature yet<sup>12-14</sup>. Since the local structure of glasses determines their physical properties, a first description of the amorphous structure of the (co-)sputtered chalcogenide films was probed by means of Raman and Fourier-transform infrared (FTIR) spectroscopies. Then, modelling of the spectroscopic ellipsometry data acquired on as-deposited thin films permits to determine their linear optical constants in the visible to near infrared range. From these experimental data and by using the well-known Sheik-Bahae model, a first evaluation of their n<sub>2</sub> Kerr refractive index is given in order to get better insight on their nonlinear optical properties. Moreover, the thermal stability of all deposited thin films was probed by monitoring their optical reflectivity during annealing under N<sub>2</sub> atmosphere. Finally, some of the GeSb<sub>w</sub>S<sub>x</sub>Se<sub>y</sub>Te<sub>z</sub> compositions are shown to exhibit an optimized compromise between the good glass stability (but limited transparency window) of S-based chalcogenide and the high 3rd order nonlinear refractive index n<sub>2</sub> (but low thermal stability of the amorphous phase) of Te-based compositions<sup>15,16</sup>. Besides, we reveal the major impact of Sb introduction in GeS<sub>x</sub> and GeSe<sub>x</sub> binary as well as GeS<sub>x</sub>Se<sub>y</sub> and GeSe<sub>x</sub>Te<sub>y</sub> ternary compounds on the amorphous structure and therefore on the optical properties of chalcogenide thin films obtained by co-sputtering deposition technique.

## Methods

**Chalcogenide thin films deposition and basic characterization.** The amorphous chalcogenide thin films were obtained by means of (co-)sputtering on 200 nm Si(001) substrates from either a single or up to simultaneously three targets of pure chalcogenide compounds. Before (co-)sputtering deposition, we did not perform any particular surface preparation of the Si substrates. Indeed, the substrates are standard high-quality microelectronic Si wafers with high purity and very low surface contamination and covered, as usual, by a few nm thick thin layer of native silicon oxide. The films' thicknesses were either fixed to 100 nm or 200 nm. The compositions of sputtering targets were Ge<sub>40</sub>S<sub>60</sub>, Ge<sub>37</sub>S<sub>63</sub>, Ge<sub>26</sub>Se<sub>74</sub>, Ge<sub>30</sub>Se<sub>70</sub>, Ge<sub>52</sub>Te<sub>48</sub>, pure Sb and pure Ge. The targets were sputtered by magnetron sputtering with an Ar plasma either achieved by using radio frequency (RF) at 13.56 MHz or direct courant (DC) plasma discharges. It is well known that deposition parameters of magnetron sputtering technique, such as Argon gas flow, substrate temperature, deposition vessel geometry, RF or DC plasma, applied power on targets and deposition pressure especially can play a major role on the density, composition and structure of the thin films<sup>17-19</sup>. In this study, all deposition were made using the same Ar flow and the deposition pressure was kept at 5.10<sup>-3</sup> mBar. Note that due to geometry and deposition conditions in our multi-cathode chamber, the substrate temperature was kept very close to room temperature during deposition of chalcogenide films.

The different compositions of the sputtered films were obtained by varying the power applied to the targets between 10 and 220 W (see Sect. I of the Supplementary Information). Note that a slight composition deviation between targets and deposited films can occur. This is the result of the different sputtering yields between chemical elements due to a selective atoms-ions interactions. This effect also depends on the deposition conditions such as applied power or aging of the targets. For instance, Ge<sub>38</sub>S<sub>62</sub> and Ge<sub>36</sub>S<sub>64</sub> thin film samples were obtained by means of sputtering of the same Ge<sub>37</sub>S<sub>63</sub> target and in the same deposition conditions. The slight deviation of composition is probably due to aging of the target after a high number of deposition cycles. This variation is of the order of the accuracy level of characterization methods used to determine films' compositions (~ 1 at. %). However, this trend on compositions is also confirmed by changes of properties and amorphous structure of these films as shown in the following.

The different compositions of the films were adjusted by varying the deposition rates ratio in between the different sputtering targets during the co-sputtering deposition. Then, the obtained compositions were verified by means of Wavelength Dispersive X-Ray Fluorescence (WDXRF) and/or ion beam analysis (Rutherford Back Scattering (RBS) and Particle Induced X-Ray Emission (PIXE) measurements).

Immediately after deposition all films were transferred under high vacuum in a second deposition chamber in order to be capped by a 10 nm thick SiN<sub>x</sub> protective layer. This nitride layer deposited by reactive RF magnetron sputtering of a pure Si target under Ar/N<sub>2</sub> reactive atmosphere is shown to be highly efficient to protect chalcogenide thin films from surface oxidation since chalcogenides are particularly prone to oxidation<sup>18,20,21</sup>. We must note that during SiN<sub>x</sub> deposition the temperature of the chalcogenide thin film samples can reach up to 180 °C due to thermal heating resulting from RF sputtering of the Si target. We note that since this temperature is below the glass transition temperature of our chalcogenide samples, it is expected to have no detrimental impact on the films, as already studied for similar chalcogenide bulk glasses and thin films<sup>4,22-27</sup>.

In order to get a first evaluation of the limit of stability of the deposited amorphous chalcogenides upon annealing, the reflectivity of the films at 670 nm was monitored during an annealing under N<sub>2</sub> atmosphere. The heating ramp rate was fixed to 10 °C/min. The limit temperature was defined as the temperature for which an irreversible change of sample reflectivity was detected indicating an irreversible structural modification of the amorphous material (phase crystallization, phase segregation, layer delamination ...).

**FTIR and Raman scattering measurements.** Fourier-Transform Infrared spectroscopy (FTIR) analysis of local order of the amorphous films was performed in transmission mode in the 100–500 cm<sup>-1</sup> range. All absorbance spectra were acquired in the same experimental conditions (average over 32 scans and spectral resolution of 2 cm<sup>-1</sup>). A reference absorbance spectrum collected on a Si(100) substrate covered by a 10 nm SiN<sub>x</sub> capping layer was used as background subtraction for all FTIR spectra acquired on the chalcogenide thin
film samples. Then, all raw spectra were normalized to the chalcogenide film thicknesses in order to get a more accurate comparison between the different thin film samples.

Raman scattering spectra were acquired in a micro-Raman spectrometer in the range from 100 to 500 cm⁻¹ using a laser probe at 532 nm wavelength. The acquisition conditions (laser power, magnification and exposure time) were adjusted for each films in order to optimize the signal-to-noise ratio but with a particular emphasis to keep no or very limited impact on the material's structure.

**Spectroscopic ellipsometry measurements.** Spectroscopic ellipsometry (SE) measurements were performed in the 400–1,700 nm range. Data were collected at three incidence angles (55, 65 and 75°). Analysis of the raw data was performed using WVASE 32° software. A 10 nm SiNₓ layer deposited on a Si substrate was also measured separately in order to take into account any possible influence of capping layer when modelling of chalcogenide films' data. For chalcogenide thin film samples, the film thicknesses, dielectric functions, optical constants (refractive index n and extinction coefficient k) and absorption coefficient α as a function of the photon energy in the 0.73–3.1 eV range were obtained by means of modelling of the SE data with a Cody-Lorentz (CL) model (see also Sect. 3 of the Supplementary Information).

The optical bandgaps of the films were estimated by using the bandgap values obtained from the CL fitting model (E<sub>g</sub><sup>CL</sup>) as well as by considering the energy for which the absorption reaches 10<sup>4</sup> cm<sup>-1</sup> (E<sub>g</sub><sup>04</sup>).

Using the M-line technique<sup>28</sup> at two wavelengths (1,313 and 1,548 nm), the effective refractive indices of the films as well as their thicknesses (by prism coupling technic) were also accurately determined (not shown). The obtained values of n at 1,313 and 1,548 nm were compared with those obtained from SE modelling and were used to validate the accuracy of SE results.

In order to get an estimation of the optical nonlinearities of the studied glasses the well-known Sheik–Bahae model was used<sup>29</sup>. This method allows to estimate nonlinear Kerr refractive index n<sub>2</sub> by means of an analytical approach using linear refractive index and optical band gap energy values. This model takes into account contributions from several physical origins: two-photon and Raman transitions, linear Stark and quadratic Stark effects. A divergent term is also added in order to subtract the unphysical behaviour resulting from the formula used to adjust these contributions. For these calculations, we used a Sellmeier fit of the refractive index in the material's transparency range obtained by using the Cody-Lorentz modelling and the optical band gap energy E<sub>g</sub><sup>04</sup>. The results of Sellmeier fits were extrapolated to wavelengths beyond spectral range of the ellipsometry measurements.

## Results and discussion

**As-deposited chalcogenide thin films' composition map.** All compositions of the studied chalcogenide thin films are reported on the ternary diagrams of Fig. 1 (see also Sect. 1 of the Supplementary Information). As shown in Fig. 1 the wide composition range of amorphous chalcogenide thin films that can be obtained exhibits a minimal thermal stability of 250 °C and up to higher than 400 °C depending on film composition (see "Methods" and Sect. 2 of the Supplementary Information).

**Amorphous structure of chalcogenide thin films by FTIR and Raman spectroscopies.** A summary of the bonds nature and main structural motifs detected in the amorphous structure of chalcogenide thin films deposited by co-sputtering is listed hereafter. All the details are given in Sect. 3 of the Supplementary Information. Table 1 summarizes the main structural motifs for each type of chalcogenide thin film compounds of the present study deduced from analysis of the Raman and FTIR vibration modes as well as the references from literature supporting our conclusions. As follows, the main conclusions on the amorphous structure drawn from the analysis of the Raman and FTIR spectra of each film is summarized system by system based on the extended discussion of the Supplementary Information.

**Ge<sub>1-x</sub>Se<sub>x</sub>, Ge<sub>1-x</sub>S<sub>x</sub> and [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> thin films.** In Fig. 2a are shown the FTIR and Raman spectra acquired on the Ge<sub>1-x</sub>Se<sub>x</sub>, thin films, with x varying in 0.63–0.74 range. Se enrichment in Ge<sub>1-x</sub>Se<sub>x</sub> films leads to reduction of homopolar inter-tetrahedral Ge–Ge bonds (Ge–Ge<sub>ETH</sub>) in favour of Se–Se ones accompanied by an increase of the number of corner-sharing (CS) compared with edge-sharing (ES) GeSe<sub>4/2</sub> tetrahedra in excellent agreement with previous works<sup>30</sup>.

The FTIR and Raman spectra of Ge<sub>1-x</sub>S<sub>x</sub> films are presented in Fig. 2b. In amorphous Ge<sub>1-x</sub>S<sub>x</sub> films, vibrational modes corresponding to GeS<sub>4/2</sub> tetrahedra are observed along with contributions attributed to Ge–Ge and S–S homopolar bonds. The decrease of disorder in the amorphous Ge<sub>1-x</sub>S<sub>x</sub> films evidenced by a hardening of the Raman modes is observed when going from Ge<sub>40</sub>S<sub>60</sub> toward Ge<sub>36</sub>S<sub>64</sub> composition as the amorphous network evolves toward composition being more and more close to that of the stoichiometric GeS<sub>2</sub> glass. The intensity of the visible contribution of the crystalline Si substrate (c-Si) depends on the thickness and transparency of the chalcogenide films at 532 nm.

In Fig. 2c, the FTIR and Raman spectra of the [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> thin films obtained by co-sputtering of Ge<sub>40</sub>S<sub>60</sub> and Ge<sub>26</sub>Se<sub>74</sub> targets exhibit the main modes related to Ge–S and Ge–Se bonds with broadening and a slight frequency shift of their intensity maxima compared with those detected in pure amorphous Ge<sub>40</sub>S<sub>60</sub> and Ge<sub>26</sub>Se<sub>74</sub> compounds<sup>31</sup>. The relative intensities of these two main contributions depend on the Ge<sub>40</sub>S<sub>60</sub>/Ge<sub>26</sub>Se<sub>74</sub> concentration ratio introduced in the film during co-sputtering. Besides, random incorporation of sulphur and selenium in mixed GeS<sub>4-m</sub>Se<sub>m</sub> tetrahedra (with m = 1, 2 and 3) is observed as supported by previous experimental and simulation studies<sup>32–35</sup>.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Ternary and pseudo-ternary diagrams showing compositions of the chalcogenide thin films deposited by magnetron (co)-sputtering: (a) Ge<sub>1-x-y</sub>S<sub>x</sub>Se<sub>y</sub>, (b) [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x-y</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub>Sb<sub>y</sub> and (c) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub> films. The dots' size denotes the error bars estimated on the compositions of the films and the dots' colour corresponds to an estimation of their limit temperature of stability defined by means of the temperature-resolved reflectivity measurements (see "Methods" and the Supplementary Information). Compounds' compositions indicated in square bracket correspond to those of the sputtering targets which may differ from those obtained for the deposited films.</figcaption></figure>

[Ge₃₀Se₇₀]₁₋ₓSbₓ and [Ge₃₇S₆₃]₁₋ₓSbₓ thin films. Figure 3a shows the Raman and FTIR spectra acquired on the [Ge₃₀Se₇₀]₁₋ₓSbₓ films. Addition of Sb in Ge₃₀Se₇₀ films results in formation of SbSe₃/₂ pyramids, which are reminiscent of the main structural units of amorphous Sb₂Se₃, connected to GeSe₄/₂ tetrahedral units by the Se atoms, and a decreasing number of Se–Se homopolar bonds in favour of a more and more prominent amount of homopolar Sb–Sb and wrong Ge–Sb bonds. This observation is of a major importance since homopolar Sb–Sb and wrong Ge–Sb bonds are shown to play a major role on electronic and hence optical properties of chalcogenide glasses³,³⁶,³⁷. Apart the non-negligible amount of Sb–Sb and Ge–Sb bonds present in our films which are chalcogen-deficient compared with the GeSe₂ and Sb₂Se₃ stoichiometric compositions, a similar trend has been previously observed in [GeSe₂]₁₋ₓ[Sb₂Se₃]ₓ bulk glasses³⁸.

The FTIR and Raman spectra of the [Ge₃₇S₆₃]₁₋ₓSbₓ thin films are shown in Fig. 3b. S-deficient [Ge₃₇S₆₃]₁₋ₓSbₓ films are mainly composed of GeS₄/₂ tetrahedral units and, upon addition of pure Sb, Sb–Sb and Sb–Ge bonds are formed with appearance of a small amount of SbS₃/₂ pyramids. Therefore, in other words the amorphous network of our S-poor glasses can be depicted as a mix of GeSb₄₋ₘSₘ tetrahedral motifs (with m = {0, 1, 2, 3, 4}) and few SbGe₃₋ₘnSbₘSₙ pyramidal units (with n + m = {0, 1, 2, 3}) in favour of homopolar Sb–Sb and wrong Sb–Ge bonds³⁹⁻⁴².

In Fig. 3c, the FTIR and Raman spectra of the [Ge₃₇S₆₃]₁₋ₓ[Ge₃₀Se₇₀]ₓ thin films upon Sb incorporation show vanishing of the Raman modes of mixed GeS₄₋ₘSeₘ tetrahedral units at least in favour of an increase number of GeSe₄/₂ tetrahedra. This indicates differences between Ge and Sb atoms in chemical bonding affinity with S and Se chalcogen elements. This has a different impact on the amorphous structure depending on the S/Se ratio. Increasing Sb concentration in [Ge₃₇S₆₃]₁₋ₓ[Ge₃₀Se₇₀]ₓSbᵧ films leads to preferential formation of Ge–Se and Sb–S bonds by detriment to Ge–S and Sb–Se ones in films exhibiting a lack of chalcogen element compared with the stoichiometric compositions in agreement with a previous study⁴¹. Raman and in a less manner FTIR
<table><caption>**Table 1.** List of bonds and their main vibration frequencies in the corresponding structural motifs as detected by Raman and FTIR analysis of chalcogenide thin films. The frequencies of the Raman and FTIR modes correspond to position of experimental peaks reported for stoichiometric glasses in the literature. The analysis of these peaks is extensively discussed in Sect. 3 of the Supplementary Information. CS, ES and ETH denotes respectively corner-sharing, edge-sharing and ethane-like motifs as commonly labelled in the literature.</caption><thead><tr><th rowspan="2">System</th><th rowspan="2">Bond</th><th rowspan="2">Main structural motifs</th><th rowspan="2">Raman modes (cm<sup>-1</sup>)</th><th rowspan="2">References</th><th colspan="2">IR modes</th></tr><tr><th>(cm<sup>-1</sup>)</th><th>References</th></tr></thead><tbody><tr><td rowspan="6">Ge<sub>1-x</sub>Se<sub>x</sub></td><td rowspan="3">Ge-Se</td><td>GeSe<sub>4/2</sub> tetrahedra</td><td>115; 140</td><td>68,69</td><td>115; 260; 285; 310</td><td>69</td></tr><tr><td>CS GeSe<sub>4/2</sub> tetrahedra</td><td>195</td><td>70</td><td>–</td><td>–</td></tr><tr><td>ES GeSe<sub>4/2</sub> tetrahedra</td><td>218; 310</td><td>70</td><td>–</td><td>–</td></tr><tr><td rowspan="3">Ge-Ge</td><td>ETH Ge<sub>2</sub>Se<sub>6</sub> units</td><td>–</td><td>–</td><td>220</td><td>69</td></tr><tr><td>ETH Ge<sub>2</sub>Se<sub>6</sub> units</td><td>179; 270</td><td>24,70</td><td>–</td><td>–</td></tr><tr><td>Ge tetrahedra in amorphous Ge phase</td><td>275</td><td>71</td><td>–</td><td>–</td></tr><tr><td rowspan="6">Ge<sub>1-x</sub>S<sub>x</sub></td><td>Se-Se</td><td>Se-Se bridge between GeSe<sub>4/2</sub> tetrahedra</td><td>265</td><td>14</td><td>–</td><td>–</td></tr><tr><td rowspan="3">Ge-S</td><td>GeS<sub>4/2</sub> tetrahedra</td><td>115; 150</td><td>72,73,74</td><td>147; 367; 388</td><td>75,76</td></tr><tr><td>CS GeS<sub>4/2</sub> tetrahedra</td><td>343; 425</td><td>70,72,76</td><td>343</td><td>75</td></tr><tr><td>ES GeS<sub>4/2</sub> tetrahedra</td><td>370; 437</td><td>70,72,76</td><td>437</td><td>75</td></tr><tr><td rowspan="3">Ge-S-Ge</td><td>GeS<sub>4/2</sub> tetrahedra</td><td>–</td><td>–</td><td>265</td><td>32</td></tr><tr><td>ETH Ge<sub>2</sub>Se<sub>6</sub> units</td><td>250</td><td>30,70</td><td>–</td><td>–</td></tr><tr><td>S rings</td><td>220; 475</td><td>73</td><td>–</td><td>–</td></tr><tr><td rowspan="6">Ge<sub>1-x-y</sub>S<sub>x</sub>Se<sub>y</sub></td><td rowspan="3">Ge-Se</td><td>S chains</td><td>485</td><td>32</td><td>–</td><td>–</td></tr><tr><td>GeS<sub>1</sub>Se<sub>3</sub> mixed tetrahedra</td><td>218</td><td>34,35</td><td>–</td><td>–</td></tr><tr><td>GeS<sub>2</sub>Se<sub>2</sub> mixed tetrahedra</td><td>232</td><td>34,35</td><td>–</td><td>–</td></tr><tr><td rowspan="3">Ge-S</td><td>GeS<sub>3</sub>Se<sub>1</sub> mixed tetrahedra</td><td>265</td><td>35</td><td>–</td><td>–</td></tr><tr><td>GeS<sub>1</sub>Se<sub>3</sub> mixed tetrahedra</td><td>392</td><td>35</td><td>–</td><td>–</td></tr><tr><td>GeS<sub>2</sub>Se<sub>2</sub> mixed tetrahedra</td><td>383</td><td>35</td><td>–</td><td>–</td></tr><tr><td rowspan="6">Ge<sub>1-x-y</sub>S<sub>x</sub>S<sub>y</sub></td><td rowspan="3">Sb-Se</td><td>GeS<sub>3</sub>Se<sub>1</sub> mixed tetrahedra</td><td>367</td><td>35</td><td>–</td><td>–</td></tr><tr><td>SbSe<sub>3/2</sub> pyramids</td><td>190</td><td>14</td><td>180; 200; 250</td><td>32,77</td></tr><tr><td>ETH Sb<sub>2</sub>Se<sub>4</sub> units</td><td>159</td><td>14</td><td>156</td><td>78</td></tr><tr><td rowspan="2">Sb-Sb</td><td>amorphous Sb phase</td><td>140</td><td>79,80</td><td>–</td><td>–</td></tr><tr><td>amorphous Ge<sub>15</sub>Sb<sub>85</sub> phase</td><td>140</td><td>81</td><td>–</td><td>–</td></tr><tr><td rowspan="3">Sb-S</td><td>SbS<sub>3/2</sub> pyramids</td><td>280; 308</td><td>82</td><td>285; 330</td><td>83,84</td></tr><tr><td>ETH Sb<sub>2</sub>Se<sub>4</sub> units</td><td>170</td><td>40</td><td>–</td><td>–</td></tr><tr><td>S<sub>2</sub>Ge-SbS<sub>2</sub> units</td><td>205</td><td>85</td><td>–</td><td>–</td></tr><tr><td rowspan="6">Ge<sub>1-x-y</sub>Se<sub>x</sub>Te<sub>y</sub></td><td>Ge-Te</td><td>GeTe<sub>4/2</sub> tetrahedra</td><td>–</td><td>–</td><td>150; 220</td><td>86</td></tr><tr><td>Ge-Te</td><td>Ge-GeTe<sub>3</sub> tetrahedra or GeTe defective octahedral motifs</td><td>120</td><td>87</td><td>–</td><td>–</td></tr><tr><td>Se-Te</td><td>Se-Se bridge between GeSe<sub>4/2</sub> tetrahedra</td><td>200</td><td>22,88</td><td>–</td><td>–</td></tr><tr><td rowspan="2">Te-Te</td><td>Te chains</td><td>150</td><td>89</td><td>–</td><td>–</td></tr><tr><td>amorphous Te phase</td><td>157</td><td>90,91</td><td>–</td><td>–</td></tr><tr><td>Ge<sub>1-x-y</sub>Sb<sub>x</sub>Se<sub>y</sub>Te<sub>z</sub></td><td>Sb-Se</td><td>SbTe<sub>3/2</sub> pyramids</td><td>145</td><td>92</td><td>–</td><td>–</td></tr><tr><td>Si substrate</td><td>Si-Si</td><td>c-Si modes</td><td>300; 520</td><td>93</td><td>–</td><td>–</td></tr></tbody></table>

spectra also evidences the presence of a significant amount of Ge-Ge, Sb-Sb and Ge-Sb bonds in our sputtered [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x-y</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub>Sb<sub>y</sub> films.

**[Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub> and [Ge<sub>1-2x</sub>Se<sub>x</sub>Te<sub>x</sub>]<sub>1-y</sub>Sb<sub>y</sub> thin films.** Figure 4a shows the FTIR and Raman spectra acquired on the [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> thin films obtained by co-sputtering of Ge<sub>30</sub>Se<sub>70</sub> and Ge<sub>52</sub>Te<sub>48</sub> targets. Ge-GeTe<sub>3</sub> tetrahedra and GeTe defective octahedral motifs present in a-GeTe phase<sup>43,44</sup> are expected also in Ge<sub>52</sub>Te<sub>48</sub>-rich [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> films. The structure of our [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> thin films obtained by co-sputtering can be depicted as GeSe<sub>4-n</sub>Te<sub>n</sub> tetrahedral motifs with n = {0, 1, 2, 3, 4} forming a disordered network connected, for small x values, by means of the chalcogen elements and, for higher x values, coexisting with a non-negligible amount of Ge-Ge homopolar bonds and possibly some Te-Te bonds as well as Se-Te bonds as proposed in a previous work<sup>22</sup>.

In Fig. 4b,c are shown the FTIR and Raman spectra of [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub> and [Ge<sub>1-2x</sub>Se<sub>x</sub>Te<sub>x</sub>]<sub>1-y</sub>Sb<sub>y</sub> thin films obtained by co-sputtering of Ge<sub>30</sub>Se<sub>70</sub>, Ge<sub>52</sub>Te<sub>48</sub> and Sb targets. The comparison between both figure tends to indicate that in our amorphous [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> films, which are getting more and more Ge-rich as x value is increased, Ge-Ge homopolars are found and Ge-Se(Te) bonds represent the vast majority of
<figure><img src="image_4.png" /><figcaption>**Figure 2.** FTIR and Raman spectra of as-deposited (a) Ge₁-xSex (b) Ge₁-xSx and (c) [Ge₄₀S₆₀]₁-x[Ge₂₆Se₇₄]ₓ thin films. The different contribution of vibration modes related to GeSe₂ and GeS₂ amorphous phases are indicated by orange and green dash-lines, respectively. The contribution of the c-Si substrate to the Raman signal appearing on some spectra is also indicated by a red dash-line at 300 cm⁻¹ as well as a shoulder visible after 450 cm⁻¹ that corresponds to a c-Si phonon mode near 520 cm⁻¹.</figcaption></figure>

bonds involving the chalcogen elements. Upon introduction of Sb in the [Ge₃₀Se₇₀]₁-x[Ge₅₂Te₄₈]ₓ films, whereas the contribution of Sb-Se bonds to the Raman spectra is weak their presence in IR absorption is more visible in FTIR spectra probably thanks to the high IR cross section of such bonds compared with others. Sb doping of the Ge₄₃Se₂₉Te₂₈ compound lead to formation of Sb-Sb homopolar bonds and probably Ge-Sb wrong bonds as already discussed above for Sb-doped Ge₃₀Se₇₀ films. As a result, [Ge₃₀Se₇₀]₁-x[Ge₅₂Te₄₈]ₓSbₓ films obtained by co-sputtering appear to be a highly disordered system. It can be depicted as a mix of Ge-GeTe₃ and GeSb₄-m-nSeₘTeₙ tetrahedra (with m + n = {0, 1, 2, 3, 4}) as well as SbGe₃-m-nSeₘTeₙ motifs (with m + n = {0, 1, 2, 3}). The relative fraction between these motifs varies as a function of the Ge/Se/Te atomic ratio of the films.

To conclude, the effect of Sb incorporation in the amorphous structure of [Ge₃₇S₆₃]₁-x[Ge₃₀Se₇₀]ₓ and [Ge₃₀Se₇₀]₁-x[Ge₅₂Te₄₈]ₓ thin films is different depending on the nature of the chalcogen elements present in the film. A clear difference is observed between films containing S/Se compared with Se/Te as chalcogen elements. In particular, this is evidenced by formation of GeSb₄-m-nSeₘTeₙ tetrahedra (with m + n = {0, 1, 2, 3, 4}) as well
<figure><img src="image_5.png" /><figcaption>**Figure 3.** FTIR and Raman spectra of (a) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>Sb<sub>x</sub> (b) [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x</sub>Sb<sub>x</sub> and (c) [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x-y</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub>Sb<sub>y</sub> thin films. The position of the main vibration modes related to GeSe<sub>2</sub> and GeS<sub>2</sub> stoichiometric glasses are indicated by orange and green dash-lines, respectively. The contribution of the c-Si substrate to the Raman signal appearing on some spectra is also indicated by a red dash-line at 300 cm<sup>-1</sup> as well as a shoulder visible after 450 cm<sup>-1</sup> that corresponds to a c-Si phonon mode near 520 cm<sup>-1</sup>.</figcaption></figure>

as SbGe<sub>3-m-n</sub>Se<sub>m</sub>Te<sub>n</sub> motifs (with m + n = {0, 1, 2, 3}) in Sb-doped [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> thin films whereas no such mixed Ge tetrahedra nor mixed Sb pyramids can be observed in [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub> compounds doped with Sb. Besides, in both case Sb addition [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub> and [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> films does not lead to a random and homogeneous distribution of the chalcogen atoms in Ge-centered tetrahedra as observed in Ge-based chalcogenide films. From all the above FTIR and Raman study, one can conclude that the local order and the structure of the (co-)sputtered chalcogenide thin films are shown to largely vary with in particular a very different amount of homopolar bonds and the latter is shown to depend on thin films' atomic composition. Since the properties of materials being intimately linked to their structure, probing the link between structure and optical properties is an invaluable clue in order to propose design rules aiming at fabricating chalcogenide compounds thin films with optimized properties for the applications in photonics. In the following, the linear and nonlinear optical constants of the films, such as the real and imaginary part of refractive index as well as Kerr nonlinear refractive index, are studied.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** FTIR and Raman spectra of (a) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>, (b) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub> and (c) [Ge<sub>1-2x</sub>Se<sub>x</sub>Te<sub>x</sub>]<sub>1-y</sub>Sb<sub>y</sub> thin films. The position of the main modes related to amorphous GeTe or GeSe<sub>2</sub> compounds as well as those related to Sb atoms are indicated by orange, purple and black dash-lines, respectively.</figcaption></figure>

**Optical properties of the chalcogenide thin films.** *Linear optical constants.* The optical constants (refractive index n and extinction coefficient k) of the films were deduced from the spectroscopic ellipsometry measurements from visible to near-IR (NIR) range (see "Methods"). The refractive indices n and Tauc's plots of (a,E)<sup>1/2</sup> vs energy (eV) obtained from the extinction coefficient k (see "Methods") are plotted in Figs. 5, 6, 7 for each composition of the chalcogenide films. First, Figs. 5, 6, 7 call for a general comment. For wavelength range located above inter band absorption range, which corresponds to energies lower than the bandgap energy (see the Tauc's plots of Figs. 5, 6, 7), the refractive indices tend progressively to a kind of plateau. The latter gives therefore an estimation of the refractive indices in the MIR range up to multi-phonons absorption appearing at significantly higher wavelengths than NIR<sup>45,46</sup>. Therefore, the study of optical properties in the visible-NIR range is the best compromise in order to get an estimation of refractive indices from visible to MIR range as well as giving an estimation of the optical band gap energy from absorption measurements or ellipsometry data fitting models (see "Methods" and Sect. 4 of the Supplementary Information).
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Refractive index n (left panels) and Tauc's plots (right panels) of (a) Ge<sub>1-x</sub>Se<sub>x</sub> (b) Ge<sub>1-x</sub>S<sub>x</sub> and (c) [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> thin films deposited by (co-)sputtering.</figcaption></figure>

**Ge<sub>1-x</sub>Se<sub>x</sub>, Ge<sub>1-x</sub>S<sub>x</sub> and [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> thin films.** In Fig. 5a, the refractive index of the sputtered Ge<sub>1-x</sub>Se<sub>x</sub> (0.63 < x < 0.74) thin films in the NIR range (above 850 nm) does not evolve linearly with composition but reaches a maximum for x = 0.63 as well as a minimum for x = 0.7. This effect can be related to the amorphous structure of Ge<sub>37</sub>Se<sub>63</sub> film which exhibits the highest amount of Ge-Ge homopolar bonds in ethane-like structures among all Ge<sub>1-x</sub>Se<sub>x</sub> (0.63 < x < 0.74) films (see the Raman of Fig. 2a). Besides, ab initio molecular dynamics simulations reported that the presence of numerous distorted Se-Ge-Se angles in Ge<sub>1-x</sub>Se<sub>x</sub> (0.60 < x < 0.66) permits to reduce the stress induced by the increase of the Ge content and thus the mean atomic coordination number or network connectivity<sup>47</sup>. By opposite, the Ge<sub>30</sub>Se<sub>70</sub> film which has the lowest refractive index also presents the lowest amount of Ge-Ge and Se-Se homopolars and the narrowest and most well-defined modes corresponding to CS and ES GeSe<sub>4/2</sub> tetrahedra (195 and 212 cm<sup>-1</sup>) as shown in Fig. 2a. This trend on the refractive index is in good agreement with previous literature for bulk glasses in which compositions with the highest Ge-Se/(Ge-Ge+Se-Se) bonding ratio exhibited the lowest refractive indices<sup>26</sup>. The absorption of Ge<sub>30</sub>Se<sub>70</sub>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Refractive index n and Tauc's plots of (a) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>Sb<sub>x</sub> (b) [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x</sub>Sb<sub>x</sub> and (c) [Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x-y</sub>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub>Sb<sub>y</sub> thin films deposited by (co-)sputtering.</figcaption></figure>

film evidences its highest band gap value among all Ge<sub>1-x</sub>Se<sub>x</sub> (0.63 < x < 0.74) films as seen on the Tauc's plots of Fig. 5a. Besides, the Ge<sub>37</sub>Se<sub>63</sub> composition with the highest refractive index has the lowest band gap energy. The decreasing number of Ge-Ge bonds may be at origin of the strong increase of absorption between Ge<sub>34</sub>Se<sub>66</sub> and Ge<sub>37</sub>Se<sub>63</sub> films.

Ge<sub>1-x</sub>S<sub>x</sub> thin films experience a significant increase of their refractive index in the NIR as the Ge fraction is slightly increased from x = 0.64 to x = 0.60 (Fig. 5b). This trend is consistent with previous literature studies<sup>48,49</sup>. This could be attributed to an increasing number of distorted S-Ge-S bonds similarly to the previously mentioned presence of numerous distorted Se-Ge-Se angles in Ge<sub>1-x</sub>Se<sub>x</sub> (0.60 < x < 0.66)<sup>47</sup> related to the increase of the number of Ge-Ge homopolar bonds (see Fig. 2b) with the decrease of the S content. This leads to a reduction of the band gap energy in agreement with previous literature<sup>49</sup> due to creation of new electronic states under the conduction band explaining the progressive shift of absorption toward lower energy (Fig. 5b)<sup>50</sup>. The absorption and refractive index change monotonously as x increases in our Ge<sub>1-x</sub>S<sub>x</sub> films (Fig. 5b).
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Refractive index n and Tauc's plots of (a) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub> (b) [Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub> and (c) [Ge<sub>1-2x</sub>Se<sub>x</sub>Te<sub>x</sub>]<sub>1-y</sub>Sb<sub>y</sub> thin films deposited by (co-)sputtering.</figcaption></figure>

The refractive index of the [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> films also exhibits a minimum for the Ge<sub>33</sub>S<sub>37</sub>Se<sub>30</sub> composition (Fig. 5c). This is not surprising since, similarly to Ge<sub>1-x</sub>Se<sub>x</sub> and Ge<sub>1-x</sub>S<sub>x</sub> films, the amount of Ge-Ge homopolar bonds is expected to be minimum for this near-stoichiometric composition with a Ge/(S + Se) atomic ratio close to 1/2 (Fig. 2c). As a result, the general trends on refractive indices of studied Ge-based chalcogenide films are driven by the ratio between heteropolar and homopolar bonds.

Among the [Ge<sub>40</sub>S<sub>60</sub>]<sub>1-x</sub>[Ge<sub>26</sub>Se<sub>74</sub>]<sub>x</sub> thin films, the band gap energy reaches a maximum for the Ge<sub>33</sub>S<sub>37</sub>Se<sub>30</sub> composition as evidenced by the absorption curves in Fig. 5c. This results from the composition of Ge<sub>33</sub>S<sub>37</sub>Se<sub>30</sub> film which is close to the GeCh<sub>2</sub> stoichiometric compound (with Ch referring to S or Se chalcogen element) expected to exhibit almost no homopolar bonds. Note that for a same Ge concentration, replacing S by Se atoms results in a progressive decrease of band gap energy<sup>31</sup>.
[Ge₃₀Se₇₀]₁₋ₓSbₓ, [Ge₃₇S₆₃]₁₋ₓSbₓ and [Ge₃₇S₆₃]₁₋ₓ[Ge₃₀Se₇₀]ₓSbᵧ thin films. In Fig. 6, replacing tetravalent Ge atoms, as observed in structural motifs of stoichiometric GeSe₂ glass, by trivalent Sb atoms, as found in amorphous stoichiometric Sb₂Se₃ compound, within the GeSₓSe₁₋ₓ films increases the non-chalcogen/chalcogen element ratio required to keep the stoichiometric composition. This has been clearly evidenced for composition tie-lines crossing the stoichiometric GeSe₂/Sb₂Se₃ pseudo-binary tie-line as Ge₄₀₋ₓSbₓSe₆₀⁵¹, GeₓSb₁₀Se₉₀₋ₓ⁵²,⁵³ GeₓSb₁₅Se₈₅₋ₓ or GeₓSb₂₀Se₈₀₋ₓ⁵³ glasses. These compositions are again those with the smallest amount of homopolar bonds and corresponding to a limit of topological phase transition as reported in Ge₄₀₋ₓSbₓSe₆₀ glasses⁵⁴,⁵⁵. Besides, in bulk glasses as well as in thin films, the decrease of the GeSe₂/Sb₂Se₃ ratio was shown to result in an increase of the refractive index and a decrease of the band gap energy since electronic polarizability of Sb–Se bonds is much higher than that of Ge-Se bonds²⁴,⁵⁶,⁵⁷.

In Fig. 6 are shown the change of optical constants upon Sb addition by means of co-sputtering in the Ge₃₀Se₇₀, Ge₃₃S₆₇ and [Ge₃₃S₆₇]₁₋ₓ[Ge₃₀Se₇₀]ₓ thin films in order to increase films' refractive indices by drifting away from compositions located on the GeCh₂/Sb₂Ch₃ pseudo-binary tie-line. In Figs. 6a–c the incorporation of Sb in Ge₃₀Se₇₀, Ge₃₃S₆₇ and [Ge₃₃S₆₇]₁₋ₓ[Ge₃₀Se₇₀]ₓ films significantly increases films' refractive index and at the same time reducing the optical band gap energy of the material. This effect is attributed to the higher electronic polarizability of Sb atoms in particular when forming Sb–Ch and Sb–Sb bonds. These highly polarizable bonds are revealed to appear in an increasing level as the incorporated Sb amount is increased in films (see amorphous structure analysis detailed in the Amorphous Structure section and in the Supplementary Information). One can also note that among all these compositions, some films exhibit compositions very close to those of well-studied commercial glasses such as for instance AMTIR-3 glass (also commercially called IG5 or IRG-5 or BD-2 or OPTIR-3) of Ge₂₈Sb₁₂Se₆₀ composition. Thus, this well-known glass is close to the Ge₂₆Sb₁₃Se₆₁ film of the present study and can be used as a point of comparison. In Fig. 6a, the optical constant values of our [Ge₃₀Se₇₀]₁₋ₓSbₓ films are in good agreement with the literature¹⁴,⁵²,⁵³.

[Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓ, [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓSbᵧ and [Ge₁₋₂xSeₓTeₓ]₁₋ᵧSbᵧ thin films. Figure 7a shows that, similarly to Sb, introduction of Ge₅₂Te₄₈ in [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓ films results in a huge increase of refractive index of the Ge₃₀Se₇₀ glass. As in the Sb doping, this effect could be attributed to the higher electronic polarizability of atoms of Ge–Te motifs compared with those in Ge-Se ones, and that form in a growing concentration in [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓ films upon increasing the Ge₅₂Te₄₈ concentration (see Fig. 4a). Moreover, the refractive index of [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓ films can be further increased by Sb incorporation as shown in Fig. 7b,c. Such an increase is near directly correlated to the concentration of Te and Sb atoms in the [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓSbᵧ and [Ge₁₋₂xSeₓTeₓ]₁₋ᵧSbᵧ films with a more and more significant effect when highly polarizable environments such as Sb–Te, Sb–Sb bonds or in a less extent Ge–Te bonds form in the amorphous material⁵⁸.

In Fig. 7a, the increase of absorption of [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓ films upon increasing x can be related to the smallest bandgap of a-GeTe compared to that of Se-based compound²⁹. The absorption progressively decreases as the band gap of the film increases upon moving toward the Ge₃₀Se₇₀ composition. By opposite, in Fig. 7b,c for [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓSbᵧ films, a more complex trend is observed. First, in Fig. 7b an increase of the (Sb + Te)/Se ratio in [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓSbᵧ films lead to a shift of the absorption curve toward lower energy due to a decrease of band gap energy. More surprising, in Fig. 7c the absorption curves and optical band gap of [Ge₁₋₂xSeₓTeₓ]₁₋ᵧSbᵧ films remain almost constant for the three compounds containing Sb concentration ranging from 15 to 30 at. %. Increasing the concentration of homopolar bonds, mainly Sb-Sb ones, and Te-related bonds due to an increase of the (Sb + Te)/Se ratio in [Ge₃₀Se₇₀]₁₋ₓ[Ge₅₂Te₄₈]ₓSbᵧ films result in an increase of band tails as well as creation of electronic states in the bandgap, and hence to a reduction of the optical bandgap energy.

However, in Fig. 7c upon incorporation of Sb in the ternary Ge₁₋₂xSeₓTeₓ compound, the refractive index of the films significantly increase while the decrease in bandgap energy is negligible. The explanation is related to the changes in the amorphous structure. First, the increasing concentration of Sb leads to an increase of the material polarizability and thus that of the refractive index. However, at the same time the amount of Te–Ge heteropolar bonds decreases leading to an increase of band gap energy of the material hence counterbalancing the effect of Sb. Therefore in [Ge₁₋₂xSeₓTeₓ]₁₋ᵧSbᵧ films, the incorporation of highly polarizable Sb bonds leads to an increase of the refractive index but does not significantly affect the optical band gap energy value. Indeed, the replacement of Ge–Ge homopolars by Sb–Sb bonds upon Sb introduction is expected to impact mostly the density and nature of localized electronic defect states in the material band gap but with no or limited effect on the bandgap energy.

To summarize, the observed trend of refractive indices as a function of chalcogenide thin films' compositions can be reasonably correlated to changes in films' amorphous structure as extensively described in the FTIR/ Raman experiments of the Amorphous Structure section and in Sect. 3 of the Supplementary Information. The electronic polarizability of local chemical environments and bonding configurations in the amorphous is shown to play key role aiming at controlling the refractive indices of the films. In particular, introduction of highly polarizable bonds such as Sb–Sb homopolars, Sb–Ch (with Ch = Te or Se) or Te–X (with X = Sb or Ge) and in a less extent Ge–Ge bonds leads to a significant increase of the refractive index concomitant to a decrease of the optical bandgap energy. The optical band gap energy is also an utmost parameter that should be taken into account for MIR and nonlinear photonic applications regarding for instance transparency window, two-photon absorption (TPA) losses and so on. The band gap energy values of the chalcogenide thin films are listed in Table 2. The obtained band gap energies and refractive index values are found to follow the well-known Moss rule⁶⁰, which relates the refractive index (or optical dielectric constant) to the optical band gap energy of semiconductors⁶¹.

In order to evaluate the potential interest for applications of such chalcogenide thin films in the emerging field of on-chip nonlinear photonics, determining the Kerr refractive indices, which quantify the nonlinear frequency conversion efficiencies, can give a first interesting insight. For that purpose, in the following the Kerr
refractive indices of the films are estimated by means of their linear optical constants as well as using estimated values of the band gap energy.

Kerr nonlinear refractive indices $n_2$. The $n_2$ Kerr nonlinear refractive indices were evaluated mainly using the Sheik-Bahae model²⁹. We considered this model as enough accurate in order to get qualitative trends on $n_2$ values (see Sect. 5 of the Supplementary Information). It has to be emphasized that a systematic error of about 20% exists in the experimental determination of non-linear optical coefficients.

The Kerr refractive index is strongly correlated with the linear refractive index and the optical band gap energy value. Using the Sheik-Bahae model, the maximum of the third order nonlinear parameter is found for an energy close to $0,534 \times E_g^{\text{opt}}$. Therefore, the values of $n_2$ can vary considerably depending on the wavelength and must be taken into consideration depending on the value of the wavelength that will be used in the application. However, previous experimental work reported a maximum of Kerr index for photon energy values higher than that corresponding to $0,534 \times E_g^{\text{opt}}$⁵³. Table 2 shows the maximum values of the $n_2$ Kerr refractive indices for all the films studied in this work and their corresponding wavelengths, as well as the values obtained at 1,550 nm, the standard wavelength for telecommunications. In order to validate the Sheik-Bahae model, this method was also applied to samples of silica and silicon nitride thin film. The $n_2$ values at 1,550 nm for these two reference materials have been calculated to be $3 \times 10^{-20}$ and $2 \times 10^{-19} \text{ m}^2/\text{W}$, respectively. These values are in excellent agreement with those experimentally measured in the literature⁶⁰. Therefore, one can conclude that in our (co-)sputtered amorphous chalcogenide thin films, the Kerr indices are of two to three orders of magnitude higher than those obtained for silica or silicon nitride materials. These values are also in excellent agreement with previous experimental results reported for similar $Ge_{1-x}Sb_xSe_y$, $Ge_{1-x}Sb_xS_y$ compounds³⁸,⁶³⁻⁶⁵ as well as those deduced using the Sheik-Bahae model for $Ge_{1-x}S_x$ and $Ge_{1-x}Sb_xSe_y$ glasses⁴⁹,⁶⁶.

For some of the chalcogenide thin films, the Sheik-Bahae model gives negative Kerr indices at 1,550 nm. In previous works, only positive experimental $n_2$ values were observed in the transparency window of chalcogenide glasses⁶¹,⁶⁷. Nevertheless, negative Kerr indices values were reported in $Ge_{1-x}Sb_xSe_y$ glasses at 800 nm⁶⁶ where absorption starts to be non-negligible. Therefore, we can emphasize that the negative Kerr indices values at 1.55 µm (~0.8 eV) are found for co-sputtered films for which the optical band gap energy is below 1.6 eV and as a result when TPA becomes significant.

In literature, $n_2$ values have been related to the amorphous structure. In particular, a clear correlation was found between $n_2$ values and the concentration of highly polarizable heteropolar bonds³⁸. Herein, one observe a clear correlation between material polarizability, which is proportional to the linear refractive index, and the nonlinearities evidenced by the $n_2$ values calculated by using the Sheik-Bahae model. Thus, one can relate the enhancement of electronic polarizability and the resulting increase of $n_2$ values to presence of peculiar local atomic motifs and bonding configurations found in the amorphous structure of films. For instance, the amount of homopolar bonds play a key role in the increase of linear and nonlinear refractive indices as described above.

Finally, in Table 2 are also reported for each films the limit of temperature after which a degradation of the material could be observed by monitoring optical reflectivity at 670 nm upon annealing (see "Methods" and Sect. 2 of the Supplementary Information). This is also very instructive since temperature limit in between 250 and 400 °C were obtained. This emphasizes that a compromise between optical properties and thermal stability has to be found for selection of a particular composition depending on final application as well as taking into account the thermal budget required for integration process flow in devices. Indeed, an annealing after deposition is also expected to significantly affect the optical properties of these Ge-based amorphous chalcogenide thin films due to the structural relaxation. Thus, aging of these metastable amorphous materials will have to be studied in the future aiming at ensuring durability of MIR components that would integrate some of these compounds.

## Conclusion

To conclude, industrial co-sputtering deposition method is a powerful tool in order to fastly study a wide compositions range of amorphous chalcogenide thin films aiming at ultimately achieving highly nonlinear on-chip MIR components. By means of a systematic study of the amorphous structure correlated with the trend on optical properties (linear optical constants, optical band gap and Kerr nonlinear refractive index) of the as-deposited chalcogenide thin films one can get invaluable clues in order to optimize the materials optical properties towards future applications. The materials' polarizability and thus linear and nonlinear refractive indices increase significantly when moving from light to heavier chalcogen element such as S to Se and toward Te-based chalcogenide compounds but accompanied with a decrease of the thermal stability. The ratio of homopolar/wrong on heteropolar bonds in the amorphous chalcogenide is shown to play a main role on the electronic and thus optical properties of the films. For instance, the introduction of Sb atoms deeply modifies the amorphous structure as well as introducing Sb-Sb homopolar bonds. As a result, the electronic polarizability of the glass is significantly enhanced as evidenced by the significant increase of the material refractive indices. However, this also leads to a detrimental decrease of the material's thermal stability and bandgap. Moreover, we show that the outstanding and state-of-the-art Kerr refractive indices in the MIR range are found for chalcogenide thin films deposited by means of industrial sputtering technique enabling fast transfer to applications. Finally, we demonstrate that a good trade-off between high nonlinearity, good thermal stability and optimized working wavelength in the IR can be found opening wide range of opportunities for future on-chip photonic applications fully compatible with CMOS large-scale integration technologies.
<table><caption>**Table 2.** Summary of chalcogenide thin films obtained by (co-)sputtering deposition: composition, nature of the sputtering targets used for (co-)sputtering deposition, refractive index at 1.55 µm, band gap energy E<sub>g</sub><sup>04</sup> and E<sub>g</sub><sup>CL</sup>, n<sub>2</sub> Kerr nonlinear refractive index calculated either by means of the Sheik-Bahae model at 1.55 µm (n<sub>2</sub>) and maximal n<sub>2</sub> values (n<sub>2 max</sub>) at energy near 0.534 × E<sub>g</sub><sup>opt</sup> eV (the corresponding wavelength value is also indicated into brackets) and a first evaluation of the limit temperature for material's stability (see text for details as well as "Methods" and Sect. 2 of the Supplementary Information).</caption><thead><tr><th>Groups</th><th colspan="2">Co-sputtered Targets</th><th>Compositions</th><th>E<sub>g</sub><sup>04</sup><br>(eV)</th><th>E<sub>g</sub><sup>CL</sup><br>(eV)</th><th>n<br>λ = 1.55μm</th><th>n<sub>2</sub><br>(10<sup>-18</sup> m<sup>2</sup>.W<sup>-1</sup>)<br>λ = 1.55μm</th><th>n<sub>2 max</sub><br>(10<sup>-18</sup> m<sup>2</sup>.W<sup>-1</sup>)</th><th>Limit of<br>thermal<br>stability<br>(°C)</th></tr></thead><tbody><tr><td rowspan="3">[Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x</sub><br>Sb<sub>x</sub></td><td rowspan="3"></td><td>Ge<sub>37</sub>S<sub>63</sub></td><td>Ge<sub>31</sub>Sb<sub>15</sub>S<sub>54</sub></td><td>1.94</td><td>1.61</td><td>2.53</td><td>8.3</td><td>19.5 (λ = 1200 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>29</sub>Sb<sub>20</sub>S<sub>51</sub></td><td>1.76</td><td>1.46</td><td>2.67</td><td>13.6</td><td>25.9 (λ = 1320 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>27</sub>Sb<sub>25</sub>S<sub>48</sub></td><td>1.61</td><td>1.34</td><td>2.79</td><td>25.2</td><td>33.9 (λ = 1450 nm)</td><td>300</td></tr><tr><td rowspan="3">Ge<sub>1-x-y</sub>S<sub>x</sub>Se<sub>y</sub></td><td rowspan="3"></td><td>Ge<sub>40</sub>S<sub>60</sub></td><td>Ge<sub>36</sub>S<sub>64</sub></td><td>2.82</td><td>2.46</td><td>2.12</td><td>1.8</td><td>6.1 (λ = 820 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>38</sub>S<sub>62</sub></td><td>2.7</td><td>2.37</td><td>2.19</td><td>2.1</td><td>6.8 (λ = 860 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>40</sub>S<sub>60</sub></td><td>2.39</td><td>1.98</td><td>2.29</td><td>3.4</td><td>10.2 (λ = 970 nm)</td><td>400</td></tr><tr><td rowspan="5">Ge<sub>1-x-y</sub>Se<sub>x</sub></td><td rowspan="5"></td><td>Ge<sub>25</sub>Se<sub>75</sub></td><td>Ge<sub>38</sub>S<sub>50</sub>Se<sub>12</sub></td><td>2.43</td><td>2.03</td><td>2.27</td><td>3.2</td><td>9.7 (λ = 960 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>33</sub>S<sub>37</sub>Se<sub>30</sub></td><td>2.57</td><td>2.3</td><td>2.24</td><td>2.5</td><td>8.0 (λ = 900 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>32</sub>S<sub>28</sub>Se<sub>40</sub></td><td>2.5</td><td>2.26</td><td>2.26</td><td>2.8</td><td>8.7 (λ = 930 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>28</sub>S<sub>15</sub>Se<sub>57</sub></td><td>2.35</td><td>2.15</td><td>2.33</td><td>3.6</td><td>10.5 (λ = 990 nm)</td><td>400</td></tr><tr><td></td><td>Ge<sub>26</sub>Se<sub>74</sub></td><td>2.22</td><td>2.11</td><td>2.42</td><td>4.4</td><td>12.3 (λ = 1050 nm)</td><td>300</td></tr><tr><td rowspan="4">[Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub><br>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub></td><td rowspan="4"></td><td>Ge<sub>30</sub>Se<sub>70</sub></td><td>Ge<sub>30</sub>Se<sub>70</sub></td><td>2.29</td><td>2.19</td><td>2.39</td><td>3.8</td><td>11.1 (λ = 1010 nm)</td><td>350</td></tr><tr><td></td><td>Ge<sub>34</sub>Se<sub>66</sub></td><td>2.13</td><td>2.05</td><td>2.42</td><td>5.3</td><td>14.4 (λ = 1090 nm)</td><td>350</td></tr><tr><td></td><td>Ge<sub>37</sub>Se<sub>63</sub></td><td>1.88</td><td>1.64</td><td>2.55</td><td>9.9</td><td>21.8 (λ = 1240 nm)</td><td>350</td></tr><tr><td></td><td>Ge<sub>34</sub>Se<sub>60</sub>Te<sub>6</sub></td><td>2.01</td><td>1.71</td><td>2.48</td><td>7.1</td><td>17.6 (λ = 1160 nm)</td><td>400</td></tr><tr><td rowspan="4">[Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x-y</sub><br>[Ge<sub>52</sub>Te<sub>48</sub>]<sub>x</sub>Sb<sub>y</sub></td><td rowspan="4"></td><td>Ge<sub>52</sub>Te<sub>48</sub></td><td>Ge<sub>37</sub>Se<sub>50</sub>Te<sub>13</sub></td><td>1.73</td><td>1.46</td><td>2.67</td><td>15.2</td><td>27.7 (λ = 1350 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>40</sub>Se<sub>39</sub>Te<sub>21</sub></td><td>1.47</td><td>1.26</td><td>2.92</td><td>43.7</td><td>44.4 (λ = 1580 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>44</sub>Se<sub>28</sub>Te<sub>28</sub></td><td>1.23</td><td>1.12</td><td>3.21</td><td>19.4</td><td>75.7 (λ = 1890 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>28</sub>Sb<sub>15</sub>Se<sub>52</sub>Te<sub>5</sub></td><td>1.54</td><td>1.29</td><td>2.94</td><td>35.6</td><td>36.6 (λ = 1510 nm)</td><td>250</td></tr><tr><td rowspan="6">[Ge<sub>30</sub>Se<sub>70</sub>]<sub>1-x</sub><br>[Sb]<sub>x</sub></td><td rowspan="6"></td><td>Ge<sub>30</sub>Sb<sub>25</sub>Se<sub>30</sub>Te<sub>15</sub></td><td>Ge<sub>29</sub>Sb<sub>20</sub>Se<sub>40</sub>Te<sub>11</sub></td><td>1.3</td><td>1.14</td><td>3.28</td><td>30.1</td><td>58.6 (λ = 1790 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>37</sub>Sb<sub>15</sub>Se<sub>24</sub>Te<sub>24</sub></td><td>1.13</td><td>0.97</td><td>3.6</td><td>-10.2</td><td>85.7 (λ = 2060 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>34</sub>Sb<sub>20</sub>Se<sub>23</sub>Te<sub>23</sub></td><td>1.12</td><td>0.95</td><td>3.6</td><td>-13.7</td><td>88.7 (λ = 2080 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>32</sub>Sb<sub>25</sub>Se<sub>22</sub>Te<sub>21</sub></td><td>1.09</td><td>0.93</td><td>3.69</td><td>-26.0</td><td>95.0 (λ = 2140 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>26</sub>Sb<sub>13</sub>Se<sub>61</sub></td><td>1.77</td><td>1.54</td><td>2.77</td><td>12.2</td><td>23.5 (λ = 1320 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>23</sub>Sb<sub>22</sub>Se<sub>55</sub></td><td>1.44</td><td>1.27</td><td>3.03</td><td>42.5</td><td>45.4 (λ = 1620 nm)</td><td>250</td></tr><tr><td rowspan="4">[Ge<sub>37</sub>S<sub>63</sub>]<sub>1-x-y</sub><br>[Ge<sub>30</sub>Se<sub>70</sub>]<sub>x</sub>Sb<sub>y</sub></td><td rowspan="4"></td><td>Ge<sub>30</sub>Sb<sub>11</sub>S<sub>33</sub>Se<sub>26</sub></td><td>Ge<sub>21</sub>Sb<sub>30</sub>Se<sub>49</sub></td><td>1.2</td><td>1.03</td><td>3.3</td><td>11.7</td><td>78.8 (λ = 1930 nm)</td><td>250</td></tr><tr><td></td><td>Ge<sub>30</sub>Sb<sub>16</sub>S<sub>31</sub>Se<sub>25</sub></td><td>2.02</td><td>1.79</td><td>2.53</td><td>6.7</td><td>16.6 (λ = 1150 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>28</sub>Sb<sub>22</sub>S<sub>29</sub>Se<sub>23</sub></td><td>1.82</td><td>1.59</td><td>2.67</td><td>10.7</td><td>22.4 (λ = 1270 nm)</td><td>300</td></tr><tr><td></td><td>Ge<sub>24</sub>Sb<sub>27</sub>S<sub>28</sub>Se<sub>21</sub></td><td>1.61</td><td>1.46</td><td>2.83</td><td>24.7</td><td>33.1 (λ = 1450 nm)</td><td>250</td></tr><tr><td></td><td></td><td></td><td></td><td>1.42</td><td>1.37</td><td>2.96</td><td>45.3</td><td>50.6 (λ = 1640 nm)</td><td>250</td></tr><tr><td></td><td></td><td></td><td>SiN<sub>x</sub></td><td>4.85</td><td>X</td><td>2.07</td><td>0.2</td><td>0.8 (λ = 480 nm)</td><td>&gt;&gt;400</td></tr></tbody></table>

Received: 6 June 2019; Accepted: 13 May 2020

Published online: 17 July 2020

## References

1. Noé, P., Vallée, C., Hippert, F., Fillot, F. & Raty, J.-Y. Phase-change materials for non-volatile memory devices: From technological challenges to materials science issues. *Semicond. Sci. Technol.* **33**, 1 (2018).

2. Ovshinsky, S. R. Reversible electrical switching phenomena in disordered structures. *Phys. Rev. Lett.* **21**, 1450 (1968).

3. Verdy, A. *et al*. Improved electrical performance thanks to Sb and N doping in Se-rich GeSe-based OTS selector devices. *IEEE 9th Int. Mem. Work., IMW* (2017).
4. Popescu, M. A. *Non Crystalline Chalcogenides* (Springer, Berlin, 2000).

5. Eggleton, B. J. *et al.* Chalcogenide photonics. *Nature Phot.* **5**, 141 (2011).

6. Tsiulyanu, D. & Ciobanu, M. Room temperature A.C. operating gas sensors based on quaternary chalcogenides. *Sensor. Actuat. B-Chem.* **223**, 95–100 (2016).

7. Lyubin, V., Klebanov, M., Feigel, A. & Sfez, B. Films of chalcogenide glassy semiconductors: New phenomena and new applications. *Thin Solid Films* **459**, 183–186 (2004).

8. Tanaka, K. Optical nonlinearity in photonic glasses. *J. Mater. Sci. Mater. Electron.* **16**, 633–643 (2005).

9. Mouawad, O. *et al.* Multioptave midinfrared supercontinuum generation in suspended-core chalcogenide fibers. *Opt. Lett.* **39**, 2684 (2014).

10. Petersen, C. R. *et al.* Mid-infrared supercontinuum covering the 1.4–13.3 µm molecular fingerprint region using ultra-high NA chalcogenide step-index fibre. *Nat. Phot.* **8**, 830–834 (2014).

11. Yu, Y. *et al.* Mid-infrared supercontinuum generation in chalcogenides. *Opt. Mat. Express* **3**, 1075 (2013).

12. Wang, H. *et al.* Compositional dependence of crystallization and structural stability in Ge-Sb-Se chalcogenide films. *J. Non-Cryst. Solids* **453**, 108–112 (2016).

13. Lin, L. *et al.* Photo-induced structural changes in Ge-Sb-Se films. *Infrared Phys. Techn.* **81**, 59–63 (2017).

14. Halenkovič, T. *et al.* Amorphous Ge-Sb-Se thin films fabricated by co-sputtering: Properties and photosensitivity. *J. Am. Ceram. Soc.* **101**, 1–13 (2018).

15. Petit, L. *et al.* Correlation between the nonlinear refractive index and structure of germanium-based chalcogenide glasses. *Mater. Res. Bull.* **42**, 2107–2116 (2007).

16. Sharma, N., Sharda, S., Sharma, V. & Sharma, P. Far-infrared investigation of ternary Ge–Se–Sb and quaternary Ge–Se–Sb–Te chalcogenide glasses. *J. Non-Cryst. Solids* **375**, 114–118 (2013).

17. Baudet, E. *et al.* Experimental design approach for deposition optimization of RF sputtered chalcogenide thin films devoted to environmental optical sensors. *Sci. Rep.* **7**, 3500 (2017).

18. Noé, P. *et al.* Impact of interfaces on scenario of crystallization of phase change materials. *Acta Mater.* **110**, 142–148 (2016).

19. Fillot, F. *et al.* Crystallization of Ge₂Sb₂Te₅ and nitrogren-doped Ge₂Sb₂Te₅ phase-change-materials thin films studied by in situ combined X-ray scattering techniques. *J. Appl. Crystallogr.* **51**, 1691–1705 (2018).

20. Smith, C., Jackson, J., Petit, L., Rivero-Baleine, C. & Richardson, K. Processing and characterization of new oxy-sulfo-telluride glasses in the Ge–Sb–Te–S–O system. *J. Solid stat.* **Chem.** **183**, 1891–1899 (2010).

21. Knotek, P. *et al.* Ageing of Ge₂Sb₁₁.₆S₆₃.₅ thin films under various conditions. *Mat. Chem. Phys.* **195**, 236–246 (2017).

22. Černošek, Z. *et al.* The properties and structure of Ge–Se–Te glasses and thin film. *J. Non-Cryst. Solids* **460**, 169–177 (2017).

23. Suriñach, S., Baró, M. D., Clavaguera-Mora, M. T. & Clavaguera, N. Glass forming ability and cristallization kinetics of allows in the GeSe₂–GeTe–Sb₂Te₃ system. *J. Non-Cryst. Solids* **111**, 113–119 (1989).

24. Olivier, M. *et al.* Photosensitivity of pulsed laser deposited Ge–Sb–Se thin films. *Opt. Mater. Express* **5**, 781–793 (2015).

25. Svoboda, R., Málek, J. & Liška, M. Correlation between the structure and relaxation dynamics of (Ge₂S)₃(Sb₂S₃)₃⁻ glassy matrices. *J. Non-Cryst. Solids* **479**, 113–119 (2018).

26. Yang, G. *et al.* Physical properties of the GeₓSe₁₋ₓ glasses in the 0<x<0.42 range in correlation with their structure. *J. Non-Cryst. Solids* **377**, 54–59 (2013).

27. Borisova, Z. *Glassy Semiconductors* (Springer, Berlin, 1981).

28. Tien, P. K., Ulrich, R. & Martin, R. J. Modes of propagating light waves in thin deposited semiconductor films. *Appl. Phys. Lett.* **14**, 291 (1969).

29. Sheik-Bahae, M., Hutchings, D. C., Hagan, D. J. & Van Stryland, E. W. Dispersion of bound electronic nonlinear refraction in solids. *J. Quantum Electron.* **27**, 1296–1308 (1991).

30. Boolchand, P., Grothaus, J., Tenhover, M., Hazle, M. A. & Grasselli, R. K. Structure of GeS₂ glass: Spectroscopic evidence for broken chemical order. *Phys. Rev. B* **33**, 5421 (1986).

31. Gu, S., Zhang, Q. & Pan, R. Microstructure and mid-infrared refractive index dispersion of Ge(SₓSe₁₋ₓ)₄ glasses. *Chalcogenide Lett.* **12**, 257–262 (2015).

32. Petit, L. *et al.* Effect of the substitution of S for Se on the structure of the glasses in the system Ge₀.₂₃Sb₀.₀₇S₀.₇₀₋ₓSeₓ. *J. Phys. Chem. Solids* **66**, 1788–1794 (2005).

33. Dongol, M., Elhady, A. F., Ebied, M. S. & Abuelwafa, A. A. Impact of sulfur content on structural and optical properties of Ge₂Se₀.₈₋ₓSₓ chalcogenide glasses thin films. *Opt. Mat.* **78**, 266–272 (2018).

34. Griffiths, J. E., Espinosa, G. P., Phillips, J. C. & Remeika, J. P. Raman spectra and athermal laser annealing of Ge(SₓSe₁₋ₓ)₂ glasses. *Phys. Rev. B* **28**, 4444–4453 (1983).

35. Xuecai, H., Guangying, S., Yu, L., Hongbo, Y. & Yonghua, L. Structure and vibrational modes of Ge–S–Se glasses: Raman scattering and ab initio calculations. *Chalcogenide Lett.* **9**, 465–474 (2012).

36. Noé, P. *et al.* Toward ultimate nonvolatile resistive memories: The mechanism behind ovonic threshold switching revealed. *Sci. Adv.* **6**, eaay2830 (2020).

37. Lee, J. H. *et al.* Unravelling interrelations between chemical composition and refractive index dispersion of infrared-transmitting chalcogenide glasses. *Sci. Rep.* **8**, 15482 (2018).

38. Olivier, M. *et al.* Wide-range transmitting chalcogenide films and development of micro-components for infrared integrated optics applications. *Opt. Mater. Express.* **4**, 525–540 (2014).

39. Pethes, I. *et al.* Atomic level structure of Ge-Sb-S glasses: Chemical short range order and long Sb–S bonds. *J. Alloy. Compd.* **774**, 1009–1016 (2019).

40. Watanabe, I. *et al.* Study on local structure in amorphous Sb–S films by Raman scattering. *J. Non-Cryst. Solids* **58**, 35–40 (1983).

41. Guery, G. *et al.* Evolution of glass properties during a substitution of S by Se in Ge₂Sb₁₂S₆₀₋ₓSeₓ glass network. *J. Non-Cryst. Solids* **358**, 1740–1745 (2012).

42. d'Acapito, F. *et al.* Role of Sb dopant in Ag:GeSₓ-based conducting bridge random access memories. *Phys. Status Solidi A* **213**, 311–315 (2016).

43. Upadhyay, M., Murugavel, S., Anbarasu, M. & Ravindran, T. R. Structural study on amorphous and crystalline state of phase change material. *J. Appl. Phys.* **110**, 083711 (2011).

44. Raty, J.-Y. Aging in phase change materials: Getting insight from simulation. *Phys. status solidi Rapid Res. Lett.* **1800590**, 1800590 (2019).

45. Klocek, P. & Colombo, L. Index of refraction, dispersion, bandgap and light scattering in GeSe and GeSbSe glasses. *J. Non-Cryst. Solids* **93**, 1–16 (1987).

46. Fang, Y. *et al.* Methods for determining the refractive indices and thermos-optic coefficients of chalcogenide glasses at MIR wavelengths. *Opt. Mater. X* **2**, 1000030 (2019).

47. Micoulaut, M. *et al.* Structure, topology, rings, and vibrational and electronic properties of GeₓSe₁₋ₓ glasses across the rigidity transition: A numerical study. *Phys. Rev. B* **88**, 054203 (2013).

48. Munzar, M. *et al.* Some optical properties of Ge–S amorphous thin film. *Curr. Appl. Phys.* **2**, 181–185 (2002).

49. Tasseva, J., Lozanova, V. & Todorov, R. Linear and non-linear optical properties of GeS₂ doped with the elements from III and V group of the periodic table. *Bulg. Chem. Commun.* **45**, 33–37 (2013).

50. Raty, J.-Y. *et al.* Aging mechanisms in amorphous phase-change materials. *Nat. Commun.* **6**, 7467 (2015).
51. Ivanova, Z. G. et al. On the structural phase transformations in Ge<sub>2</sub>Sb<sub>40-x</sub>Se<sub>x</sub> glasses. *J. Non-Cryst. Solids* **293–295**, 580–585 (2001).

52. Sharda, S., Sharma, N., Sharma, P. & Sharma, V. Band gap and dispersive behavior of Ge alloyed a-SbSe thin films using single transmission spectrum. *Mater. Chem. Phys.* **134**, 158–162 (2012).

53. Wang, T. et al. Systematic z-scan measurements of the third order nonlinearity of chalcogenide glasses. *Opt. Mater. Express* **4**, 1011–1022 (2014).

54. Pamukchieva, V., Szekeres, A., Sharlandjiev, P., Alexieva, Z. & Gartner, M. Effet of composition on refractive index dispersion in GeSbS thin films. *Proc. Int. Semiconductor. Con.* **1**, 101–104 (1998).

55. Fabian, M., Dulgheru, N., Antonova, K., Szekeres, A. & Gartner, M. Investigation of the atomic structure of Ge–Sb–Se chalcogenide glasses. *Adv. Condens. Matter Phys.* **2018**, 1–11 (2018).

56. Abdel-Wahab, F., Ali Karar, N. N., El Shaikh, H. A. & Salem, R. M. Effect of Sb on the optical properties of the Ge–Se chalcogenide thin films. *Phys. B* **422**, 40–46 (2013).

57. Shaaban, E. R., Kaid, M. A., El Sayed, M. & Adel, A. Effect of compositional variations on the optical properties of Sb–Ge–Se thin films. *J. Phys. D Appl. Phys.* **41**, 125301 (2008).

58. Shportko, K. et al. Resonant bonding in crystalline phase-change materials. *Nat. Mater.* **7**, 653–658 (2008).

59. Noé, P. & Hippert, F. Structure and properties of chalcogenide materials for PCM. In *Phase Change Memory* 125–179 (Springer, Cham, 2018).

60. Moss, T. S. A relationship between the refractive index and the infrared threshold of sensitivity for photoconductors. *Proc. Phys. Soc. Sect. B* **63**, 167–176 (1950).

61. Zakery, A. & Elliott, S. R. *Optical Nonlinearities in Chalcogenide Glasses and their Applications* (Springer, Berlin, 2007).

62. Wang, L. et al. Nonlinear silicon nitride waveguides based on a PECVD deposition platform. *Opt. Express* **26**, 9645–9654 (2018).

63. Choi, J. W., Chen, G. F. R., Ng, D. K. T., Ooi, K. J. A. & Tan, D. T. H. Wideband nonlinear spectral broadening in ultra-short ultra-silicon rich nitride waveguides. *Sci. Rep.* **6**, 27120 (2016).

64. Kuriakose, T. et al. Measurement of ultrafast optical Kerr effect of Ge-Sb-Se chalcogenide slab waveguides by the beam self-trapping technique. *Optics Comm.* **403**, 352–357 (2017).

65. Serna, S. et al. Nonlinear optical properties of integrated GeSbS chalcogenide waveguides. *Photon. Res.* **6**, B37–B42 (2018).

66. Chen, Y. et al. Optical and structure properties of amorphous Ge–Sb–Se films for ultrafast all-optical signal processing. *J. Alloy. Compd.* **580**, 578–583 (2013).

67. Romanova, E. et al. Measurement of non-linear optical coefficients of chalcogenide glasses near the fundamental absorption band edge. *J. Non-Cryst. Solids* **480**, 13–17 (2018).

68. Sugai, S. Stochastic random network model in Ge and Si chalcogenide glasses. *Phys. Rev. B* **35**, 1345–1361 (1987).

69. Fukunaga, T., Tanaka, Y. & Murase, K. Glass formation and vibration properties in the (Ge, Sn)-Se system. *Solid State Commun.* **42**, 513–516 (1982).

70. Jackson, K. et al. Raman-active modes of a-GeSe<sub>2</sub> and a-GeS<sub>2</sub>: A first-principles study. *Phys. Rev. B* **60**, 14985–14989 (1999).

71. Lannin, J. S., Maley, N. & Kshirsagar, S. T. Raman scattering and short range order in amorphous germanium. *Solid State Commun.* **53**, 939–942 (1985).

72. Tao, H., Mao, S., Dong, G., Xiao, H. & Zhao, X. Raman scattering studies of the Ge-In sulphide glasses. *Solid State Commun.* **137**, 408–412 (2006).

73. Guo, H. T., Zhang, M.-J., Xu, Y.-T., Xia, X.-S. & Yang, Z.-Y. Structural evolution study of additions of Sb<sub>2</sub>S<sub>3</sub> and CdS into GeS<sub>2</sub> chalcogenide glass by Raman spectroscopy. *Chin. Phys. B* **26**, 104208 (2017).

74. Lucovsky, G., Galeener, F. L., Keezer, R. C., Geils, R. H. & Six, H. A. Structural interpretation of the infrared and Raman spectra of glasses in the alloy system Ge<sub>1-x</sub>S<sub>x</sub>. *Phys. Rev. B* **10**, 5134–5146 (1974).

75. Lucovsky, G., Nemanich, R. J. & Solin, S. A. Coordination dependent vibrational properties of amorphous semiconductor alloys. *Solid State Commun.* **17**, 1567–1572 (1975).

76. Julien, C. et al. Raman and infrared spectroscopic studies of Ge–Ga–Ag sulphide glasses. *Mat. Sci. Eng. B* **22**, 191–200 (1994).

77. Quiroga, I. et al. Infrared studies of a Ge<sub>0.26</sub>Sb<sub>0.05</sub>Se<sub>0.75</sub> glassy semiconductor. *J. Non-Cryst. Solids* **196**, 183–186 (1996).

78. Petkov, K., Vassilev, G., Todorov, R., Tasseva, J. & Vassilev, V. Optical properties and structure of thin films from the system GeSe<sub>2</sub>–Sb<sub>2</sub>Se<sub>3</sub>–AgI. *J. Non-Cryst. Solids* **357**, 2669–2674 (2011).

79. Lannin, J. S. Raman scattering properties of amorphous As and Sb. *Phys. Rev. B* **15**, 3863–3871 (1977).

80. Rossow, U. et al. Growth mode of ultrathin Sb layers on Si studied by spectroscopic ellipsometry and Raman scattering. *Appl. Surf. Sci.* **63**, 35–39 (1993).

81. Shakhvorostov, D. et al. Evidence for electronic gap-driven metal-semiconductor transition in phase-change materials. *PNAS* **106**, 10907–10911 (2009).

82. Lin, C. et al. Network structure in GeS<sub>2</sub>–Sb<sub>2</sub>S<sub>3</sub> chalcogenide glasses: Raman spectroscopy and phase transformation study. *J. Phys. Chem. C* **116**, 5862–5867 (2012).

83. Droichi, M. S., Vaillant, F., Bustarret, E. & Jousse, D. Study of localized states in amorphous chalcogenide Sb<sub>2</sub>S<sub>3</sub> films. *J. Non-Cryst. Solids* **101**, 151–155 (1988).

84. Kamitsos, E. I., Kapoutsis, J. A., Culeac, I. P. & Iovu, M. S. Structure and bonding in As-Sb-S chalcogenide glasses by infrared reflectance spectroscopy. *J. Phys. Chem. B* **101**, 11061–11067 (1997).

85. Nazabal, V. et al. Sputtering and pulsed laser deposition for near- and mid-infrared applications: A comparative study of Ge<sub>25</sub>Sb<sub>10</sub>S<sub>65</sub> and Ge<sub>25</sub>Sb<sub>10</sub>Se<sub>65</sub> amorphous thin films. *J. Appl. Ceram. Technol.* **8**, 990–1000 (2011).

86. Raty, J.-Y. et al. Vibrational properties and stabilization mechanism of the amorphous phase of doped GeTe. *Phys. Rev. B* **88**, 014203 (2013).

87. Andrikopoulos, K. S., Yannopoulos, S. N., Kolobov, A. V., Fons, P. & Tominaga, J. Raman scattering study of GeTe and Ge<sub>2</sub>Sb<sub>2</sub>Te<sub>5</sub> phase-change materials. *J. Phys. Chem. Solids* **68**, 1074–1078 (2007).

88. Mendoza-Galvan, A., García-García, E., Vorobiev, Y. V. & González-Hernández, J. Structural, optical and electrical characterization of amorphous Se<sub>x</sub>Te<sub>1-x</sub> thin film alloys. *Microelectron. Eng.* **51–52**, 677–687 (2000).

89. Varma, G. S. et al. Thermally reversing window in Ge<sub>15</sub>Te<sub>85-x</sub>In<sub>x</sub> glasses: Nanoindentation and micro-Raman studies. *J. Non-Cryst. Solids* **358**, 3103–3108 (2012).

90. Uemura, O., Hayasaka, N., Tokairin, S. & Usuki, T. Local atomic arrangement in Ge–Te and Ge–S–Te glasses. *J. Non-Cryst. Solids* **205–207**, 189–193 (1996).

91. Brodsky, M. H., Gambino, R. J., Smith, J. E. & Yacoby, Y. The Raman spectrum of amorphous tellurium. *Phys. Stat. sol. B* **52**, 609 (1972).

92. Němec, P. et al. Amorphous and crystallized Ge–Sb–Te thin films deposited by pulsed laser: Local structure using Raman scattering spectroscopy. *Mater. Chem. Phys.* **136**, 935–941 (2012).

93. Hu, J. et al. Exploration of waveguide fabrication from thermally evaporated Ge–Sb–S glass films. *Opt. Mater.* **30**, 1560–1566 (2008).

## Acknowledgements

E. Henaff is gratefully acknowledged for his availability and support during sputtering deposition of the chalcogenide thin films. We also sincerely acknowledge L. Duraffourg, head of the LCO laboratory of the DOPT department of CEA-LETI, for providing funding and related financial support for PhD thesis of J.-B. Dory. This
work has been partially supported by the EIPHI Graduate School (under contract"ANR-17-EURE-0002") and
the FEDER Program of the Burgundy Region.

**Author contributions**

P.N. conceptualized and managed the project with help of C.S., B.C. and J.-B. J. J.-B. D. performed all experi-
ments and treated all experimental data under supervision of P.N. and B.C. J.-B. D., C. C.-C., C.S., M.T. and
M.B. prepared all the thin film samples in LETI clean rooms. J.-M. F. supported J.-B. D. for the supervision of
the technological lots in LETI 200 mm technological platform. J.-B. D. performed spectroscopic ellipsometry
data modelling under supervision of P.N., B.C. and A.C. J.-B.D. performed all Raman and FTIR experiments
with invaluable help of A.V. The paper was written by J.-B. D. and P.N. with help of B.C. and A.V. All authors
have given their approval to the final version of the manuscript.

Competing interests

The authors declare no competing interests.

Additional information

**Supplementary information** is available for this paper at https://doi.org/10.1038/s41598-020-67377-9.

Correspondence and requests for materials should be addressed to P.N.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020