Open Access Article
Kum-Hui Kima,
Kyong-Sik Ju
*a,
Song-Hui Kima,
Hyon-Hui Ria,
Su-Bom Ria and
Kwang-Phil Kimb
aHigh-Tech Research and Development Center, Kim Il Sung University, Pyongyang, Democratic People's Republic of Korea. E-mail: ks.ju1025@ryongnamsan.edu.kp
bFaculty of Life Science, Kim Il Sung University, Pyongyang, Democratic People's Republic of Korea
First published on 24th April 2026
In this study, a silk fiber/graphene oxide (GO) composite was directly obtained by feeding silkworm larvae with GO. Methods such as temperature-dependent Fourier transform infrared (FT-IR) spectroscopy, second-derivative analysis and two-dimensional (2D) correlation FT-IR spectroscopy were used to investigate the bonding interactions between GO and silk fibers in the silk fiber/GO composite. In addition, the silk fiber/GO composite was physically characterized by scanning electron microscopy (SEM) and thermogravimetric analysis (TGA). Temperature-dependent second-derivative FT-IR spectra of the reinforced silk fiber/GO composite showed a characteristic absorption band of GO at 1053 cm−1 and 40 °C, confirming the incorporation of GO into the silk fiber matrix. 2D correlation analysis of the temperature-dependent FT-IR spectra indicated that the asynchronous correlation intensity of the silk fiber/GO composite decreased 1.7-fold compared with that of the silk fiber, and the number of correlation peaks was significantly reduced, which suggests that GO present in the silk fiber/GO composite is thermally reduced with increasing temperature, resulting in enhanced thermally stability. The results of TGA showed that the silk fiber/GO composite exhibits higher decomposition temperatures than the silk fiber, which is in good agreement with the temperature-dependent 2D correlation analysis.
Several methods such as temperature-dependent FT-IR spectroscopy, first- and second-derivative analysis and 2D correlation spectroscopy have been used to study graphene-based polymers. The conformational changes of polyvinyl alcohol (PVA) have been studied by temperature-dependent FT-IR spectroscopy at different temperatures.18 Mahendia19 studied the increase in glass transition temperature of reduced GO–PVA composites using the first derivative curves of temperature-dependent FT-IR spectra of reduced GO–PVA composites at different concentrations of reduced GO incorporated into PVA. In general, first and second-derivative techniques are both commonly used as spectral processing methods to resolve overlapping peaks.20
2D correlation FT-IR spectroscopy can successfully separate overlapping peaks and has been widely used to investigate the bonding mechanism in graphene-based composites, polymers and carbon nanomaterials.24–29 Su et al.24 applied temperature-dependent FT-IR spectroscopy with generalized 2D correlation spectroscopy to investigate the phase transition mechanism of poly(N-isopropylacrylamide-co-2-hydroxyethyl methacrylate) (PNIPAM-co-HEMA) hydrogels. They suggested that the phase transition can be divided into two stages (21.8–31.4 °C and 31.4–36.5 °C), and the two phase-transition processes can be further divided into nine steps based on generalized 2D correlation analysis. Shao et al.25 used second-derivative analysis and 2D FT-IR correlation spectroscopy to investigate conformation changes in glucose oxidase induced by adsorption onto GO sheets. The thermal reduction process of GO to graphene was already shown to occur in two stages by temperature-dependent 2D correlation and moving-window correlation, and the reduction mechanism was discussed in detail.30
In the study, a second-derivative method was combined with temperature-dependent FT-IR spectroscopy to detect the characteristic absorption bands of the epoxide and ether groups of GO in the silk fiber/GO composite. In addition, 2D correlation analysis was performed to investigate the bonding interactions between GO and silk fiber in the silk fiber/GO composite.
Bombyx mori silkworm larvae in the fifth instars were divided into two groups and raised at 25 °C and 20% relative humidity, being fed either with GO-treated leaves or normal mulberry leaves.
The cocoons obtained from the two groups were cut into appropriate sizes and washed with distilled water. They were then boiled in a 1-g per L Na2CO3 solution at 85 °C for 30 min to remove sericin, followed by washing with distilled water to remove residual Na2CO3. The degummed silk fibers were dried in an oven at 60 °C for 8 h.
The silk fibers obtained from silkworms fed with normal mulberry leaves and those fed with GO-treated leaves were denoted as the silk fiber and silk fiber/GO composite, respectively.
The synchronous and asynchronous 2D correlation FTIR spectra were calculated and processed using the MATLAB functions.23
Raman spectra were recorded on a LabRam HR800 Raman analyzer using a 532.17 nm Nd-YAG laser. For Raman measurements, 100 mg of degummed silk fiber/GO composite was dissolved in 20 mL of concentrated nitric acid, and this solution was diluted 30 times using distilled water. Finally, the dissolved silk fibroin was removed by vacuum filtration, and the resulting precipitate was used for Raman spectroscopic analysis.
Mechanical testing was performed using an AG-X Plus instrument (Shimadzu, Japan) at a stretching speed of 1 mm min−1, and the stress–strain curves were collected automatically by a computer.
![]() | ||
| Fig. 1 SEM images of the silk fiber (a) and silk fiber/GO composite (b). (c) Stress–strain curves of silk fibers. TG and DTG curves of silk fiber (d), TG and DTG curves silk fiber/GO composite (e). | ||
The mechanical testing was performed on two groups. The typical stress–strain curves for the two groups are shown in Fig. 1c, which indicate higher values for the silk fiber/GO composite compared to the silk fiber. Table 1 shows a summary of the average diameters and mechanical properties of the two groups.
| Types | Diameter (µm) | Tensile strength (MPa) | Elongation at break (%) | Young's modulus (GPa) |
|---|---|---|---|---|
| Silk fiber | 12.2 ± 1.7 | 430 ± 25 | 14.7 ± 1.3 | 29.4 ± 1.2 |
| Silk fiber/GO composite | 11.3 ± 1.5 | 580 ± 32 | 16.4 ± 1.6 | 35.9 ± 1.4 |
The tensile strength and Young's modulus of the silk fiber were 430 MPa and 29.4 GPa, respectively, while those of the silk fiber/GO composite were 580 MPa and 35.9 GPa.
TG and DTG (derivative thermogravimetric) curves of both fibers were recorded from room temperature up to 800 °C under an N2 atmosphere to evaluate their thermal degradation properties, and the measurement results are shown in Fig. 1d and e. The initial weight loss observed around 100 °C for the two groups was due to water evaporation. The silk fiber exhibited significant weight loss at about 280.7 °C and showed a peak degradation temperature at about 305.6 °C. For the silk fiber/GO composite, these weight losses occurred at 291.2 °C and 313.4 °C, respectively, which can be attributed to the presence of GO in the silk fiber/GO composite.
The TG and DTG results demonstrate that the incorporation of GO improves the thermal stability of silk fibers, resulting in a slow thermal degradation process in the silk fiber/GO composite.
![]() | ||
| Fig. 2 Temperature-dependent FT-IR spectral profiles of the silk fiber (a) and silk fiber/GO composite (b). | ||
Fig. 3 shows the temperature-dependent spectral intensity changes of the amide I (1635 cm−1), amide II (1519 cm−1) and amide III (1230 cm−1) bands over the temperature range of 20 °C to 330 °C. The intensity changes of the amide I, II and III bands in both fibers were significantly different with increasing temperature. The spectra of silk fiber showed an absorbance change peaks at 252 °C, whereas the silk fiber/GO composite displayed three absorbance change peaks at 252 °C, 272 °C and 312 °C.
The maximum temperature changes of the spectral intensities corresponding to the amide bands in the silk fiber and silk fiber/GO composite are shown in Table 2. The temperature-dependent FT-IR spectra show that the thermal behavior of the silk fibers/GO composite differs significantly from that of the silk fiber. These differences between the two fibers were investigated using second-derivative treatment of the temperature-dependent FT-IR spectra. Overlapping peaks can be resolved into two or more peaks by applying the second-derivative treatment. The second-derivative spectra were obtained using a Savitzky–Golay filter through 7-point smoothing.
| Wavenumber (cm−1) | Silk fiber | Silk fiber/GO composite | ||
|---|---|---|---|---|
| T1 (°C) | T1 (°C) | T2 (°C) | T3 (°C) | |
| 1635 (amide I) | 252 | 252 | 272 | 312 |
| 1520 (amide II) | 252 | 252 | 272 | 312 |
| 1230 (amide III) | 252 | 252 | 272 | 312 |
The second-derivative spectra are more complex and significantly narrower than the first-derivative spectra, which is preferred for analyzing the thermal behavior of silk fibers/GO composite. The temperature-dependent second-derivative FT-IR spectra of the two fibers, recorded upon heating from 30 °C to 320 °C in the 3000–2800 and 1200–1000 cm−1 regions, are illustrated in Fig. 4 and 5.
![]() | ||
| Fig. 4 Temperature-dependent second-derivative FT-IR spectra of the (a) silk fiber and (b) silk fiber/GO composite in the 3000–2800 cm−1 region. | ||
![]() | ||
| Fig. 5 Temperature-dependent second-derivative FT-IR spectra of the (a) silk fiber and (b) silk fiber/GO composite in the 1200–1000 cm−1 region. | ||
The temperature-dependent second-derivative FT-IR spectra of both fibers in the 3000–2800 cm−1 region were evidently different. The bands at 2932, 2915 and 2850 cm−1 can be assigned to the stretching vibrations of the methylene group of glycine, and the bands at 2981, 2861 and 2905 cm−1 can be attributed to the stretching vibrations of the methyl group of alanine and the tertiary CH stretching vibrations of alanine. The second-derivative FT-IR spectra of the silk fiber (Fig. 4a) showed no changes between 30 °C and 70 °C. In contrast, in the silk fiber/GO composite (Fig. 4b), a new peak at 2905 cm−1 appeared at 50 °C, and the band intensity increased.
The most obvious differences between the temperature-dependent second-derivative FT-IR spectra in Fig. 5a and b were observed at 40 °C and 50 °C. A peak at 1053 cm−1 appeared at 40 °C in the silk fiber/GO composite, representing the characteristic absorption band of the C–O–C groups in GO.30 In addition, the peaks at 1053 cm−1 and 1123 cm−1 disappeared simultaneously at 50 °C. This can be attributed to the weakening or breaking of intermolecular bonds between GO and the silk fiber with increasing temperature.
Fig. 6 shows the temperature-dependent second-derivative infrared spectra of both fibers in the 1650–1500 cm−1 region recorded upon heating from 30 °C to 80 °C. The C
C skeletal vibrations in GO appear near 1570 cm−1. As shown in Fig. 6a, the absorbance intensity at 1570 cm−1 for the silk fiber hardly changes with increasing temperature. On the other hand, for the silk fiber/GO composite (Fig. 6b), the absorbance intensity at 1570 cm−1 decreases at 30 °C and 40 °C, accompanied by a wavenumber shift to 1566 cm−1. The absorbance intensity increases and shifts towards higher wavenumbers above 50 °C. Hence, it indicates an increase in the band intensity corresponding to the C
C skeletal vibrations at 1570 cm−1, relative to the decrease in the C–O bands at 1053 cm−1, suggesting that the trace amount of GO incorporated into the silk fibers is thermally reduced upon heating.
![]() | ||
| Fig. 6 Temperature-dependent second-derivative FT-IR spectra of the (a) silk fiber and (b) silk fiber/GO composite in the 1650–1500 cm−1 region. | ||
The above results can confirm the incorporation of GO into the silk fiber/GO composite, but they provide limited information about the bonding interactions between the silk fiber and GO.
Fig. 7 shows the Raman spectra of GO and the precipitate obtained from the silk fiber/GO composite. As shown in Fig. 7(1), the Raman spectrum of GO exhibits a G band at 1584 cm−1 and a D band at 1346 cm−1. The G band at 1584 cm−1 is related to the sp2 graphite structure, and the D band at 1346 cm−1 is assigned to the local defects in the carbon material,17 which is also evident in the Raman spectrum of the precipitate obtained from the silk fiber/GO composite (Fig. 7(2)). The results indicate that GO was incorporated into the silk fibers through the diet. However, using Raman spectroscopy to confirm the presence of GO in silk fiber composites has the disadvantage of very complex sample pretreatment.
![]() | ||
| Fig. 7 Raman spectra of the (1) GO and (2) precipitate obtained from the silk fiber/GO composite in the 2000–1000 cm−1 region. | ||
2D correlation spectroscopy has been used as a powerful tool to analyze changes in physical and chemical processes in response to external perturbations such as temperature.32 2D correlation analysis of temperature-dependent FT-IR spectra was carried out to compare the temperature behavior and bonding modes of the silk fiber and silk fiber/GO composite.
In Fig. 8a, the auto-correlation spectrum shows prominent auto-correlation peaks at 1643, 1540 and 1041 cm−1; on the other hand, the asynchronous correlation spectrum (Fig. 8b) shows the cross-correlation peaks between them. The asynchronous correlation spectrum (Fig. 8d) shows the cross-correlation peaks between 1643, 1540, 1508 and 1643 cm−1.
The correlation wavenumbers in the 2D correlation spectra, along with the signs of the synchronous and asynchronous correlation peaks and the direction of intensity changes at each corresponding wavenumber in the dynamic spectra, are summarized in Table 3.
The assignments of the correlation wavenumbers listed in Table 3 are summarized in Table 4.
| Frequency (cm−1) | Assignments |
|---|---|
| 1041 | C–N stretching vibration of phenylalanine |
| 1187 | C–N stretching vibration of tyrosine |
| 1288 | Amide II in β-sheet conformation |
| 1415 | CH2 bending in asparagine |
| 1508 | Amide II [N–H bending and C–N bending] |
| 1540 | Amide II [peptide bond 60% N–H bending and 40% C–N bending] |
| 1643 | Amide I in β-sheet conformation |
Based on the results in Tables 3 and 4, the directions of intensity changes and the order of spectral variations at individual wavenumbers can be derived as follows: the spectral intensities at 1643 and 1508 cm−1 decrease prior to the spectral intensity change (↘↗↘) observed at 1540, 1415, 1288, and 1187 cm−1. Subsequently, a further decrease in spectral intensity at 1187 cm−1 occurs with increasing temperature.
The 2D correlation spectra of the silk fiber (Fig. 8) and the silk fiber/GO composite (Fig. 9) show significant differences in the synchronous and asynchronous correlation spectra. The intensity ranges of the synchronous and asynchronous correlation spectra for the silk fiber are −5.6 × 10−4–0.0123 and −0.0011–0.0011, respectively, while for the silk fiber/GO composite, they are −5.14 × 10−4–0.013 and −6.402 × 10−4–6.402 × 10−4, respectively. The intensities of the synchronous correlation spectra are similar between the two fibers, but the intensities of the asynchronous correlation spectra differ markedly. The asynchronous correlation strength of the silk fiber/GO composite decreases by about 1.7-fold compared to that of the silk fiber. The silk fiber/GO composite shows the smaller correlation intensities compared to the silk fiber upon heating; i.e., it was thermally stable 1.7-fold. In addition, the synchronous correlation spectra of the silk fiber/GO composite show new peaks at 1261 cm−1 and 926 cm−1 compared to the silk fiber (Fig. 9b). The asynchronous correlation spectra of the silk fiber/GO composite also differed greatly, with a significant decrease in the number of asynchronous correlation peaks compared to that of the silk fiber (Fig. 9d). The decrease in the correlation intensity and the decrease in the number of these correlation peaks in the asynchronous correlation spectra of the silk fiber/GO composite indicate that the correlation intensity changes in FT-IR absorption bands upon external heating is small compared to that of the silk fiber.
Table 5 shows the correlation wavenumbers of the synchronous and asynchronous temperature-dependent correlation spectra for the silk fiber/GO composite and their intensity change directions.
The assignments of the interested FT-IR absorption bands of GO are summarized in Table 6.
Based on the signs and directions of the correlation signals in the 2D correlation spectra of the silk fiber/GO composite, the order of intensity changes of the absorption bands with increasing temperature can be derived as follows:
The spectral intensities at 1542, 1643 and 1666 cm−1 decrease, while those at 990 and 1203 cm−1 show a decrease → increase → decrease pattern. On the other hand, the spectral intensities at 990 and 1203 cm−1 decrease prior to the spectral intensity changes at 1542 and 1643 cm−1. Subsequently, the spectral intensity at 921 cm−1 decreases, followed by a spectral intensity decreases at 1666, 1573, 1566 and 1554 cm−1 with increasing temperature.
In summary, the order of intensity changes is 990 and 1203 cm−1 (↘↗↘) → 1542, 1643, and 1666 cm−1 (↓) → 921 cm−1 (↘↗↘) → 1554 cm−1 (↓) → 1566 cm−1 (↓) → 1666 cm−1 (↓).
Table 7 lists the assignments of the correlation wavenumbers shown in Table 5.
The comparison of the 2D correlation analysis results for silk fiber and silk fiber/GO composites shows significant differences. In silk fibers, changes in the amide I and II bands, corresponding to the silk fiber structure, occur first, and then structural changes occur in order. In silk fiber/GO composites, spectral intensity changes occur first at the 990 cm−1 band, corresponding to the C–C stretching vibration of the glycine–alanine units in silk fiber, and at the 1203 cm−1 band, corresponding to the C–O stretching vibration of the epoxide group in GO, followed by sequential structural changes in the silk fiber. Therefore, it can be seen that the thermal stability of the silk fiber/GO composite is improved compared with that of silk fiber. In addition, as shown in Tables 3 and 4, the number of asynchronous correlation peaks in the silk fiber/GO composite decreased compared to silk fiber. From these results, it can be expected to be due to the rearrangement of the silk fiber structure by GO, resulting in a more stable structure.
These results demonstrate that the relationship between the structure and thermal stability of the silk fiber and the silk fibers/GO composite can be effectively investigated by temperature-dependent 2D correlation FT-IR spectroscopy.
The characteristic peaks of GO were observed in the 2D correlation FT-IR spectra of the silk fiber/GO composite, which confirmed the incorporation of GO into the silk fiber/GO composite. It was also found that the asynchronous correlation intensity of the silk fiber/GO composite decreased by around 1.7-fold compared to that of the silk fiber with increasing temperature, and the number of correlation peaks was significantly reduced, indicating enhanced thermal stability due to the presence of GO in the silk fiber/GO composite.
The results of TGA are in good agreement with the temperature-dependent 2D correlation FT-IR analysis. This study shows the usefulness of temperature-dependent 2D correlation FT-IR spectroscopy for detecting trace amounts of GO in GO composites and for investigating their bonding interactions with the matrix. The optimal content of GO required to produce silk fibers with the desired mechanical strength should be determined in future work.
| This journal is © The Royal Society of Chemistry 2026 |