Pyrolysis Characteristics and Non-Isothermal Kinetics of Integrated Circuits

Due to the complexity of components and high hazard of emissions, thermochemical conversions of plastics among waste-integrated circuits (ICs) are more favorable compared with the common treatment options of electronic waste (E-waste), such as chemical treatment and burning. In this study, the waste random-access memory, as the representative IC, was used to investigate the thermal degradation behaviors of this type of E-waste, including a quantitative analysis of pyrolysis characteristics and non-isothermal kinetics. The results show that the pyrolysis of the ICs can be divided into three different decomposition stages. The pyrolysis temperature and gas atmosphere play an important role in the pyrolysis reaction, and the heating rate greatly affects the rate of the pyrolysis reaction. The non-isothermal kinetic parameters and reaction mechanisms of ICs are determined using the Friedman method, Coats and Redfern (CR) method, and Kissinger method. The results show that the actual average activation energy of the pyrolysis reaction of ICs should be between 170 and 200 kJ·mol−1. The optimally fitting model for the ICs pyrolysis is the three-step parallel model consisting of the random nucleation model (Am) and reaction order model (Cn).


Introduction
Over the past few decades, the electronic industries have been pushing factors of the modern industry's advancement. With rapid technological advancement and increasing customer demand, the use of consumer electronic items has increased dramatically, and the products are becoming obsolete much faster than ever before [1,2]. Approximately 45 million tons of electronic waste (E-waste) were produced in 2012 around the world [3,4] and according to UN statistics, in 2018, there were 48.5 million tons of E-waste worldwide. E-waste is growing rapidly [5].
Even though a great deal of effort has been directed toward the promotion of harmless recycling of E-waste [6][7][8], recycling E-waste without harming the environment and human health is still a problem. E-waste generated by developed countries is usually exported to developing countries, where recycling techniques include burning and dissolution in strong acids with few measures to protect human health and the environment [9,10]. Over the past few decades, developing countries, especially China and India, have recycled the most Ewaste worldwide at the expense of the environment and the health of local residents [11,12]. The amount of research on recycling E-waste without harming the environment and human health is not enough, and more research is required on the recycling technology of E-waste.
Integrated circuits (ICs), which play a vital role in electronic products as processors or memorizers, are widely used in electronic products such as computers, mobile phones, and The recovery of ICs must consider not only the recovery of resources but also the environmental impact of the process. ICs are encapsulated by epoxy resin with added PBDEs. Serious environmental problems will be caused by improper recycling processes [14][15][16].
The following E-waste treatment options can be used for reference. Common E-waste treatment options are burning [17,18], chemical treatment [7,19,20], and mechanical separation [6,8,21,22]. Burning is the most common method to recycle E-waste. However, the process of treatment releases a large number of harmful gasses, such as polycyclic aromatic hydrocarbons (PAHs), hydrogen bromide (HBr), and even dioxins bromide. Burning residue can contaminate the soil and groundwater near the incineration site. Emissions of harmful gas and burning residue pose serious health risks to local residents [11,12,15]. Chemical treatment is efficient and relatively clean. However, a large amount of waste residue and liquid that are difficult to handle are caused by chemical treatment [7]. Lastly, mechanical separation is a good method in theory, but in fact, mechanical separation can only be used for the primary separation of the metal fraction and non-metal fraction from E-waste. The high local temperature of E-waste caused by mechanical separation can also produce harmful gases [6].
Supercritical fluid extraction [23][24][25] and pyrolysis [13,[26][27][28] are promising methods for the harmless treatment of E-waste. Pyrolysis refers to the thermochemical decomposition of carbonaceous organic material in the absence or lack of oxygen through indirect heating to produce valuable fuels, such as pyrolysis gas, pyrolysis oil, and carbon black [28]. Wang and Zhang investigated the degradation process of brominated flame retardant (BFR) and BFR-containing waste computer housing plastic in various supercritical fluids (water, methanol, isopropanol, and acetone). Li and Xu used supercritical water (SCW) to decompose BER and recover metals from Waste Memory Modules (WMMs). Niu et al. used SCW to decompose the organics and recover Ta from Waste tantalum capacitors (WTCs). However, the technology of supercritical fluid extraction has strict requirements on the material of the equipment. More importantly, it is very expensive. Pyrolysis is a promising thermal technique, which can recover and remove organic material for resource recovery effectively and environmentally friendly. Diao et al. performed a series of molecular dynamics simulations (MDSs) to study the thermal decomposition characteristics of epoxy resin. Chen et al. researched vacuum pyrolysis characteristics and kinetic analysis of liquid crystal from scrap liquid crystal display panels. Liu et al. applied vacuum pyrolysis to degrade the organic material of waste ICs. Pyrolysis is considered a better treatment option than other options. Pyrolysis has the advantages of low environmental costs and high efficiency [13].
Little research has been conducted on the pyrolysis of waste ICs [13]. For this reason, the present study aims to investigate the thermal degradation characteristics and nonisothermal pyrolysis kinetics of ICs. Thermogravimetric analysis (TGA) is used to analyze the pyrolysis characteristics of organic material in Ics in a non-isothermal environment. The influence of pyrolysis-related factors on pyrolysis is studied. The non-isothermal kinetic parameters and kinetic reaction mechanisms of Ics are obtained by the Friedman method, Coats and Redfern (CR) method, and Kissinger method. The results of this study can provide process parameters and guidance for further industrial application of Ics pyrolysis technology to recycle ICs efficiently and environmentally friendly.

Materials
In this study, random-access memory (RAM), which was a representative IC, was used for the experiment. The RAM was selected as the sample for the following reasons. Firstly, the compositional characteristics of RAM are representative, and the pyrolysis characteristics of RAM can provide references for other types of ICs. Secondly, the application of RAM is widespread, and RAM generally has a higher content of precious metals and therefore has a higher recycling value than other types of ICs, such as a central processing unit, external memory, and so on. The samples were provided by the Zhongming Environmental Protection Group, which is an environmental company specializing in hazardous waste disposal and recycling. Prior to experiments, the materials were cut into small pieces by a shearing machine, and then ground into powder by a ball mill. The homogeneity of experimental samples can be guaranteed by turning the RAM into powder.
The composition of samples was determined by an Elemental Analyzer (Vario EL, Elementar Analysensysteme GmbH, Germany), X-Ray Fluorescence Spectrometer (Zetium, PANalytical B.V., Holland), Inductively Coupled Plasma-Atomic Emission Spectrometer (Prodigy 7, Leeman, America), and High-Performance Liquid Chromatography (HP1100, Agilent, America). The main components of ICs are listed in Table 1. Cu, C, Br, Si, and O are the main elements of ICs, accounting for more than 80 % of the components of ICs.

Experimental Methods
The pyrolysis experiments were performed in a tubular furnace on a laboratory scale. In order to control the atmosphere of the pyrolysis experiment, TGA was employed to deal with harmful waste gas, and a gas control system, an analytical balance system, and a tail gas treatment system were added, as shown in Figure 1. The analytical balance system consisted of an analytical balance, data line, and computer. During the experiment, the change in sample quality with time was transmitted to the computer through the data line. The change in temperature with time in the experiment was realized by measuring the temperature in the furnace with a thermocouple. deal with harmful waste gas, and a gas control system, an analytical balance system, and a tail gas treatment system were added, as shown in Figure 1. The analytical balance system consisted of an analytical balance, data line, and computer. During the experiment, the change in sample quality with time was transmitted to the computer through the data line. The change in temperature with time in the experiment was realized by measuring the temperature in the furnace with a thermocouple. During the experiment, the sample was put into the crucible, and then the crucible was lifted by stainless-steel wire with oxidation resistance and corrosion resistance. The other end of the stainless-steel wire was connected to the analytical balance. The crucible was located in the middle of the quartz tube of a tubular furnace. The temperature of the experiment was controlled according to the setting procedure of the control panel of the tubular furnace. At the same time, the gas pressure relief valve and the gas flow control valve were opened to control the atmosphere of the pyrolysis experiment. Then the vacuum pump was opened to control the gas to pass through the tail gas treatment system. The airtightness of the whole experimental device was guaranteed.

Thermogravimetric Analysis
The sample quantity measured by the thermogravimetric analyzer was 2~5 mg, which is suitable for the thermogravimetric analysis of the pure substance. The relative error among the measurements of the thermogravimetric analyzer was approximately 3~5 % [29,30]. While the powder sample made from the RAM was a mixture with an uneven composition, the relative error would be greater if the thermogravimetric analyzer was directly used for thermogravimetric analysis. So, the above experimental equipment was designed for thermogravimetric analysis of samples, as is shown in Figure 1. During the experiment, the sample was put into the crucible, and then the crucible was lifted by stainless-steel wire with oxidation resistance and corrosion resistance. The other end of the stainless-steel wire was connected to the analytical balance. The crucible was located in the middle of the quartz tube of a tubular furnace. The temperature of the experiment was controlled according to the setting procedure of the control panel of the tubular furnace. At the same time, the gas pressure relief valve and the gas flow control valve were opened to control the atmosphere of the pyrolysis experiment. Then the vacuum pump was opened to control the gas to pass through the tail gas treatment system. The airtightness of the whole experimental device was guaranteed.

Thermogravimetric Analysis
The sample quantity measured by the thermogravimetric analyzer was 2~5 mg, which is suitable for the thermogravimetric analysis of the pure substance. The relative error among the measurements of the thermogravimetric analyzer was approximately 3~5% [29,30]. While the powder sample made from the RAM was a mixture with an uneven composition, the relative error would be greater if the thermogravimetric analyzer was directly used for thermogravimetric analysis. So, the above experimental equipment was designed for thermogravimetric analysis of samples, as is shown in Figure 1.
For each experimental run, approximately 5 g of the sample was used. The samples were heated from 50 to 1050 • C at a constant heating rate of 10 • C/min. The system pressure is atmospheric pressure, and the gas atmosphere during pyrolysis is nitrogen and air. The flow rate of the carrier gas was fixed at 0.5 L/min. The relative error was controlled to less than 1%. In order to investigate the factors affecting the pyrolysis of ICs, the independent variables and control variables are shown in Table 2. Approximately 10 g of waste ICs were used in each experiment. The system pressure was one atmospheric pressure, and the flow rate of the carrier gas was fixed at 0.5 L/min. Samples before and after pyrolysis were analyzed by elemental analysis (C, H, N) to better understand the influence of different pyrolysis factors on the pyrolysis of ICs. Under any given condition, the experiment was usually carried out more than twice.

Pyrolysis Kinetics
The nomenclature used is shown in Table 3. The reaction equation of the pyrolysis of ICs can be written as follows: The non-isothermal kinetics for the pyrolysis of ICs is written as follows [30]: where α is the conversation rate and is given by In the above equation, m, m i , and m f represent the instantaneous, initial (at 50 • C), and final (at 1050 • C) weights of the sample. f(α) is the mechanical function of the reaction kinetics model and can be written as Common reaction kinetics models and mechanical functions are summarized in Table S1 in the Supplementary Materials. The reaction rate constant k is expressed in terms of the Arrhenius equation as: where A is the pre-exponential factor, E a is the activation energy, R is the universal gas constant, T is the absolute temperature, n is the reaction order. For a constant heating rate β: Substituting Equations (4)-(6) into Equation (2) gives: Different methods, namely the Friedman method of Equation (8), CR method of Equations (9) and (10), and Kissinger method of Equation (11), are used in this work, as shown in Equations (8)- (11) [29,30].

Pyrolysis Properties of ICs
There are few direct studies on the pyrolysis properties of ICs [13], but there are many studies on the pyrolysis properties of waste PCBs, which are similar to those of ICs [5,14,18,31,32]. These studies and the component analysis of ICs can provide some basic information about their behavior as well as potential products and applications in their thermochemical conversion processes.
In the pyrolysis processes, carbon and hydrogen are the main reaction elements, as the organic components of ICs are mainly composed of carbon and hydrogen. Organic components in ICs will be decomposed into small organic molecules and released by heating. Carbon monoxide and carbon dioxide will also be generated in an air atmosphere. The pyrolysis of elemental nitrogen under an air atmosphere can produce nitrogen oxides (NO x ). Unlike waste PCBs, ICs do not contain chlorine, and the bromine content is relatively high (9.17 wt%). Pyrolysis will produce bromine-based toxic gases, including hydrogen bromide (HBr), polybrominated biphenyls (PBBs), and even bromide dioxin. This requires methods to reduce the emissions of these toxic gases, and more research is needed on the migration and transformation of bromine in ICs. Copper, tin, and other metal elements are oxidized under the air atmosphere, but not under the nitrogen atmosphere. Noble metal elements will not react under any gas atmosphere.
Pyrolysis residue under the air atmosphere is mainly composed of metal oxides and glass cloth. Pyrolysis residue under a nitrogen atmosphere consists of metal, inorganic carbon, and glass cloth. The method of separating the glass cloth from the pyrolysis residue is the key to the recovery of ICs.

TG Analysis
The TGA and DTG curves of ICs under nitrogen and air atmospheres are shown in Figure 2. In general, the pyrolysis processes of ICs are characterized by a three-stage thermal degradation. In the first stage, the quality of the ICs changes negligibly. Some research shows there is PBDEs emission in the first stage [14,31]. The decomposition of ICs mainly takes place in the second stage. In this stage, the organic macromolecular polymer is thermally decomposed into a large number of volatile organic small molecules. In the third stage, as the temperature increases, the benzene ring begins to decompose.
due is the key to the recovery of ICs.

TG Analysis
The TGA and DTG curves of ICs under nitrogen and air atmospheres are shown in Figure 2. In general, the pyrolysis processes of ICs are characterized by a three-stage thermal degradation. In the first stage, the quality of the ICs changes negligibly. Some research shows there is PBDEs emission in the first stage [14,31]. The decomposition of ICs mainly takes place in the second stage. In this stage, the organic macromolecular polymer is thermally decomposed into a large number of volatile organic small molecules. In the third stage, as the temperature increases, the benzene ring begins to decompose. The third-stage reactions develop at approximately 410~874 °C. After the pyrolysis reaction, there is a weight loss of 19.92 wt% for ICs under a nitrogen atmosphere and 19.54 wt% for ICs under an air atmosphere. The weight loss in pyrolysis of ICs is lower than the weight loss in pyrolysis of PCBs [13,32]. Because the organic content of ICs is usually lower than that of PCBs, the weight loss between ICs under nitrogen and ICs under air is similar. The results of elemental analysis before and after pyrolysis are shown in Figure 3. The pyrolysis residue of the reaction in nitrogen contains pyrolysis residual carbon, while the pyrolysis residue of the reaction in the air does not contain residual carbon, and the metal part of ICs will be oxidized by oxygen. From the similar pyrolysis mass loss, it can be concluded that the mass increase in the oxidized metal part is approximately equal to the mass of the pyrolysis residual carbon. The peak temperature of ICs under a nitrogen atmosphere (338 °C) is lower than that of ICs under an air atmosphere (358 °C), but the decomposition intensity of the former (−0.4529 wt.%•°C −1 ) is substantially higher than that of the latter (−0.2024 wt.%•°C −1 ). The pyrolysis reaction speed under the nitrogen condition is faster than that under the air condition. The pyrolysis reaction time under the nitrogen condition is shorter than that under the air condition. The initial temperature of the pyrolysis reaction under the nitrogen condition is lower than that under the air condition. The micro morphology of the samples before and after pyrolysis were analyzed by SEM, The third-stage reactions develop at approximately 410~874 • C. After the pyrolysis reaction, there is a weight loss of 19.92 wt% for ICs under a nitrogen atmosphere and 19.54 wt% for ICs under an air atmosphere. The weight loss in pyrolysis of ICs is lower than the weight loss in pyrolysis of PCBs [13,32]. Because the organic content of ICs is usually lower than that of PCBs, the weight loss between ICs under nitrogen and ICs under air is similar. The results of elemental analysis before and after pyrolysis are shown in Figure 3. The pyrolysis residue of the reaction in nitrogen contains pyrolysis residual carbon, while the pyrolysis residue of the reaction in the air does not contain residual carbon, and the metal part of ICs will be oxidized by oxygen. From the similar pyrolysis mass loss, it can be concluded that the mass increase in the oxidized metal part is approximately equal to the mass of the pyrolysis residual carbon. The peak temperature of ICs under a nitrogen atmosphere (338 • C) is lower than that of ICs under an air atmosphere (358 • C), but the decomposition intensity of the former (−0.4529 wt.%· • C −1 ) is substantially higher than that of the latter (−0.2024 wt.%· • C −1 ). The pyrolysis reaction speed under the nitrogen condition is faster than that under the air condition. The pyrolysis reaction time under the nitrogen condition is shorter than that under the air condition. The initial temperature of the pyrolysis reaction under the nitrogen condition is lower than that under the air condition. The micro morphology of the samples before and after pyrolysis were analyzed by SEM, and the results are shown in Figure 3b-d. Before pyrolysis, the samples are mainly composed of metal particles and glass fibers, filled with organic materials. After pyrolysis under an N 2 atmosphere, the organic materials are completely pyrolyzed when the metal particles change slightly. On the surface of the sample, there is a small amount of residual carbon. After pyrolysis under a nitrogen atmosphere, the surface of the sample becomes rougher, and the metal particles are coated with an oxide film. Therefore, compared with pyrolysis under an air atmosphere, pyrolysis under a nitrogen atmosphere can preserve the original state of the metal parts in the waste ICs, which is conducive to further refining for metal recovery.
der an N2 atmosphere, the organic materials are completely pyrolyzed when the metal particles change slightly. On the surface of the sample, there is a small amount of residual carbon. After pyrolysis under a nitrogen atmosphere, the surface of the sample becomes rougher, and the metal particles are coated with an oxide film. Therefore, compared with pyrolysis under an air atmosphere, pyrolysis under a nitrogen atmosphere can preserve the original state of the metal parts in the waste ICs, which is conducive to further refining for metal recovery.

Effect of Related Factors on Pyrolysis
The temperature, heating rate, and gas atmosphere are important factors affecting pyrolysis. In this study, the effect of temperature, heating rate, and gas atmosphere on the pyrolysis progress and products were studied. The results show that pyrolysis temperature and gas atmosphere have an important influence on the pyrolysis reaction, and the heating rate only affects the rate of the pyrolysis reaction and has no effect on the reaction results.

Effect of Temperature
In order to investigate the effect of temperature on the pyrolysis of ICs, constanttemperature pyrolysis experiments were carried out at 250, 400, 600, and 800 °C. In the experiment, the tube heater is heated to a specified temperature, and then the sample is suspended in the tube heater. The selection of temperature-independent variables is based on the TGA and DTG curves of ICs under a nitrogen atmosphere.
The results are shown in Figure 4. At 250 °C, there is a considerable decrease in quality, although the rate of change is relatively slow. The PBDEs emission can explain this change [14,31]. At 400 °C, the reaction in the second stage of pyrolysis occurs. The organic macromolecular polymer is thermally decomposed into a large number of volatile organic small molecules. At the same time, most aromatic compounds containing benzene rings formed by pyrolysis are also released. Most organic components in ICs will be emitted at 400 °C. However, there are still a few organic compounds that do not react from elemental analysis of pyrolysis residue. The reaction rate of the second stage at 600 °C is quicker than that at 400 °C. Additionally, the reaction of the third stage can occur at 600 °C, but the reaction rate is slow. All organic components become residual carbon and volatile organic compounds from elemental analysis of pyrolysis residue. At 800 °C, the result of pyrolysis is similar. The reaction rate of the second stage is also similar compared to that at 600 °C. Nevertheless, the reaction rate in the third stage of pyrolysis at 800 °C is much faster than that at 600 °C. In conclusion, pyrolysis temperature has an important influence

Effect of Related Factors on Pyrolysis
The temperature, heating rate, and gas atmosphere are important factors affecting pyrolysis. In this study, the effect of temperature, heating rate, and gas atmosphere on the pyrolysis progress and products were studied. The results show that pyrolysis temperature and gas atmosphere have an important influence on the pyrolysis reaction, and the heating rate only affects the rate of the pyrolysis reaction and has no effect on the reaction results.

Effect of Temperature
In order to investigate the effect of temperature on the pyrolysis of ICs, constanttemperature pyrolysis experiments were carried out at 250, 400, 600, and 800 • C. In the experiment, the tube heater is heated to a specified temperature, and then the sample is suspended in the tube heater. The selection of temperature-independent variables is based on the TGA and DTG curves of ICs under a nitrogen atmosphere.
The results are shown in Figure 4. At 250 • C, there is a considerable decrease in quality, although the rate of change is relatively slow. The PBDEs emission can explain this change [14,31]. At 400 • C, the reaction in the second stage of pyrolysis occurs. The organic macromolecular polymer is thermally decomposed into a large number of volatile organic small molecules. At the same time, most aromatic compounds containing benzene rings formed by pyrolysis are also released. Most organic components in ICs will be emitted at 400 • C. However, there are still a few organic compounds that do not react from elemental analysis of pyrolysis residue. The reaction rate of the second stage at 600 • C is quicker than that at 400 • C. Additionally, the reaction of the third stage can occur at 600 • C, but the reaction rate is slow. All organic components become residual carbon and volatile organic compounds from elemental analysis of pyrolysis residue. At 800 • C, the result of pyrolysis is similar. The reaction rate of the second stage is also similar compared to that at 600 • C. Nevertheless, the reaction rate in the third stage of pyrolysis at 800 • C is much faster than that at 600 • C. In conclusion, pyrolysis temperature has an important influence on the pyrolysis reaction. It affects the final composition of reactants and reaction rate and can be controlled to obtain different pyrolysis results. Figure 4c-f shows the micro morphology of the samples before and after pyrolysis at different temperatures. For pyrolysis at 250 • C, a preliminary melting phenomenon appeared on the surface of the sample. Compared with pyrolysis at 250 • C, the surface reaction of samples pyrolyzed at 400 • C was more complete, and most of the organic components also reacted. After pyrolysis at 800 • C, the sample almost completely melted, and the organic components were fully reacted. However, the surface of the residue was seriously damaged. The higher the temperature, the more serious the damage to the part of the RAM that cannot be pyrolyzed. Hence, an appropriate temperature of 400 • C is suggested. Under this condition, the rate of reaction can be guaranteed while protecting the metal part from destruction.
Materials 2022, 15, x FOR PEER REVIEW 9 on the pyrolysis reaction. It affects the final composition of reactants and reaction rate can be controlled to obtain different pyrolysis results. Figure 4c-f shows the micro morphology of the samples before and after pyroly different temperatures. For pyrolysis at 250 °C, a preliminary melting phenomenon peared on the surface of the sample. Compared with pyrolysis at 250 °C, the surface tion of samples pyrolyzed at 400 °C was more complete, and most of the organic com nents also reacted. After pyrolysis at 800 °C, the sample almost completely melted, the organic components were fully reacted. However, the surface of the residue was ously damaged. The higher the temperature, the more serious the damage to the pa the RAM that cannot be pyrolyzed. Hence, an appropriate temperature of 400°C is gested. Under this condition, the rate of reaction can be guaranteed while protecting metal part from destruction.

Effect of Heating Rate
Heating rate is an important factor affecting the pyrolysis reaction, and 5 k/mi k/min, and 15 k/min are set to study the effect of heating rate on pyrolysis of ICs. results are shown in Figure 5.

Effect of Heating Rate
Heating rate is an important factor affecting the pyrolysis reaction, and 5 k/min, 10 k/min, and 15 k/min are set to study the effect of heating rate on pyrolysis of ICs. The results are shown in Figure 5. Figure 4c-f shows the micro morphology of the samples before and after pyrolysis a different temperatures. For pyrolysis at 250 °C, a preliminary melting phenomenon ap peared on the surface of the sample. Compared with pyrolysis at 250 °C, the surface reac tion of samples pyrolyzed at 400 °C was more complete, and most of the organic compo nents also reacted. After pyrolysis at 800 °C, the sample almost completely melted, an the organic components were fully reacted. However, the surface of the residue was ser ously damaged. The higher the temperature, the more serious the damage to the part o the RAM that cannot be pyrolyzed. Hence, an appropriate temperature of 400°C is sug gested. Under this condition, the rate of reaction can be guaranteed while protecting th metal part from destruction.

Effect of Heating Rate
Heating rate is an important factor affecting the pyrolysis reaction, and 5 k/min, 1 k/min, and 15 k/min are set to study the effect of heating rate on pyrolysis of ICs. Th results are shown in Figure 5.  At the heating rate of 5 k/min, the reaction time is the longest, while at the heating rate of 20 k/min, the reaction time is the shortest. The heating rate has a great influence on the reaction time of pyrolysis of ICs. The temperature at which the pyrolysis begins does not vary with the heating rate. Obviously, the heating rate has little influence on the reaction rate of the second stage of pyrolysis. This means that the reaction rate of the second stage is affected very little by the heating rate. The reaction rate of the third stage of pyrolysis increases as the heating rate increases. This result reflects that the reaction rate of the third stage is greatly influenced by the heating rate. The reaction time of the third stage can be shortened by increasing the heating rate of the third stage in future industrial applications. The results of elemental analysis of pyrolysis residue at different heating rates indicate the heating rate does not influence the final composition of reactants. In conclusion, the heating rate affects the rate of pyrolysis reaction. The desired reaction time can be obtained by controlling the heating rate, and the heating rate will not affect the results of the pyrolysis reaction.

Effect of Gas Atmosphere
Pyrolysis experiments in the different gas atmospheres often mean different reaction processes. It is important to study the influence of the gas atmosphere on pyrolysis experiments. The independent variables and control variables considering the influence of the gas atmosphere on the pyrolysis of ICs are shown in Table 2. In order to study the effect of the different gas atmospheres on the pyrolysis reaction at different temperatures, the temperature is added as another independent variable.
The results are shown in Figure 6. The pyrolysis of ICs at different temperatures under nitrogen is different from the pyrolysis of ICs at different temperatures under air. At 250 • C, the reaction rate in air is quicker than that in nitrogen. There is a weight loss of 7.04 wt.% for ICs under nitrogen and 9.8 wt.% for ICs under air. The difference in weight loss can be explained by the difference in the carbon content of pyrolysis residue, and the pyrolysis residue under air has a lower carbon content (8.42 wt.%) compared with the carbon content (12.84 wt.%) of pyrolysis residue under nitrogen. Whether ICs have a pyrolysis reaction at 250 • C under air needs further study. At 400 • C, the reaction rate in air is slower than that in nitrogen. From the results of elemental analysis, the residue under air contains little carbon (0.54 wt.%). However, residue under nitrogen has an 8.41 wt.% carbon content. At 800 • C, the reaction rate is almost the same between air and nitrogen. Pyrolysis residue under air does not have carbon and pyrolysis residue under nitrogen also contains a considerable amount of carbon (6.43 wt.%).
the reaction time of pyrolysis of ICs. The temperature at which the pyrolysis begins d not vary with the heating rate. Obviously, the heating rate has little influence on the re tion rate of the second stage of pyrolysis. This means that the reaction rate of the seco stage is affected very little by the heating rate. The reaction rate of the third stage of rolysis increases as the heating rate increases. This result reflects that the reaction rat the third stage is greatly influenced by the heating rate. The reaction time of the third st can be shortened by increasing the heating rate of the third stage in future industrial plications. The results of elemental analysis of pyrolysis residue at different heating ra indicate the heating rate does not influence the final composition of reactants. In con sion, the heating rate affects the rate of pyrolysis reaction. The desired reaction time be obtained by controlling the heating rate, and the heating rate will not affect the res of the pyrolysis reaction.

Effect of Gas Atmosphere
Pyrolysis experiments in the different gas atmospheres often mean different react processes. It is important to study the influence of the gas atmosphere on pyrolysis exp iments. The independent variables and control variables considering the influence of gas atmosphere on the pyrolysis of ICs are shown in Table 2. In order to study the ef of the different gas atmospheres on the pyrolysis reaction at different temperatures, temperature is added as another independent variable.
The results are shown in Figure 6. The pyrolysis of ICs at different temperatures der nitrogen is different from the pyrolysis of ICs at different temperatures under air 250 °C, the reaction rate in air is quicker than that in nitrogen. There is a weight los 7.04 wt.% for ICs under nitrogen and 9.8 wt.% for ICs under air. The difference in wei loss can be explained by the difference in the carbon content of pyrolysis residue, and pyrolysis residue under air has a lower carbon content (8.42 wt.%) compared with carbon content (12.84 wt.%) of pyrolysis residue under nitrogen. Whether ICs have a rolysis reaction at 250 °C under air needs further study. At 400 °C, the reaction rate in is slower than that in nitrogen. From the results of elemental analysis, the residue un air contains little carbon (0.54 wt.%). However, residue under nitrogen has an 8.41 w carbon content. At 800 °C, the reaction rate is almost the same between air and nitrog Pyrolysis residue under air does not have carbon and pyrolysis residue under nitro also contains a considerable amount of carbon (6.43 wt.%).  In summary, the gas atmosphere has little effect on the pyrolysis reaction rate. The pyrolysis of ICs in different atmospheres has different reaction processes and results. The differences in pyrolysis reaction at different atmospheres need further investigation.

Reaction Activation Energy
The experimental data of the pyrolysis of ICs under a nitrogen atmosphere at different heating rates were used to analyze the non-isothermal kinetics. Figure 7a shows the relationship between ln dα dt and 1 T obtained by the Friedman method. When α is less than 0.65, the linear correlation between ln dα dt and 1 T is good, and the slope changes very little, indicating that the activation energy does not change much. When α is greater than 0.65, the linear correlation between A and B is low, and the slope of the dotted line becomes larger, which means that the activation energy becomes larger, and the accuracy of this method decreases. This shows the complexity of the pyrolysis reaction. When the reaction conversion rate is low, the pyrolysis reaction is dominated by one type of reaction, and the activation energy of the reaction varies slightly. When the reaction conversion rate is high, the pyrolysis reaction is jointly dominated by various types of reactions and thus the activation energy becomes larger. In addition, the drastic changes in the activation energy of the reaction can be seen at both ends of the reaction conversion rate, indicating that this method amplifies the random error here. Figure 7b shows the relationship between ln β T 2 m and 1 T m obtained by the Kissinger method. The activation energy can be obtained from the slope and is 170.90 kJ·mol −1 , which is close to the activation energy obtained by the Friedman equation when the reaction conversion rate is lower than 0.7. The pyrolysis kinetics model of the CR equation method is a first-order reaction. The relationship between ln( − ln(1−α) T 2 ) and 1 T is shown in Figure 7c. It can be seen that the relationship is not completely linear. The curve can be divided into two sections for linear fitting, and the results are shown in Figure S1 of Supplementary Materials. The linear correlation of the piecewise fitting is good, and R 2 is above 0.9. It can be seen that the activation energy obtained by fitting in the first section is smaller than that obtained by the first two methods, while the activation energy obtained in the second section is very small, which deviated from the normal range of common organic pyrolysis activation energy. Therefore, the CR equation may not be suitable for the complex kinetics of the solid thermal decomposition reaction.

Reaction Activation Energy
The experimental data of the pyrolysis of ICs under a nitrogen atmosphere at different heating rates were used to analyze the non-isothermal kinetics. Figure 7a shows the relationship between ln ( ) and 1 T obtained by the Friedman method. When α is less than 0.65, the linear correlation between ln ( ) and 1 T is good, and the slope changes very little, indicating that the activation energy does not change much. When α is greater than 0.65, the linear correlation between A and B is low, and the slope of the dotted line becomes larger, which means that the activation energy becomes larger, and the accuracy of this method decreases. This shows the complexity of the pyrolysis reaction. When the reaction conversion rate is low, the pyrolysis reaction is dominated by one type of reaction, and the activation energy of the reaction varies slightly. When the reaction conversion rate is high, the pyrolysis reaction is jointly dominated by various types of reactions and thus the activation energy becomes larger. In addition, the drastic changes in the activation energy of the reaction can be seen at both ends of the reaction conversion rate, indicating that this method amplifies the random error here. Figure 7b shows the relationship between ln ( ) and is shown in Figure 7c. It can be seen that the relationship is not completely linear. The curve can be divided into two sections for linear fitting, and the results are shown in Figure S1 of Supplementary Materials. The linear correlation of the piecewise fitting is good, and R 2 is above 0.9. It can be seen that the activation energy obtained by fitting in the first section is smaller than that obtained by the first two methods, while the activation energy obtained in the second section is very small, which deviated from the normal range of common organic pyrolysis activation energy. Therefore, the CR equation may not be suitable for the complex kinetics of the solid thermal decomposition reaction. ) and obtained from CR method. Table 4 shows the values of activation energy obtained by the three methods. It is worth noting that the adopted activation energy of the CR equation is the average of the ) and 1 T obtained from CR method. Table 4 shows the values of activation energy obtained by the three methods. It is worth noting that the adopted activation energy of the CR equation is the average of the activation energy of the first section. The reaction activation energy obtained by the CR method is relatively small, while that obtained by the Friedman method is relatively large. When considered comprehensively, the actual average activation energy of the pyrolysis reaction of the waste memory module should be between 170 and 200 kJ·mol −1 . The obtained results on the activation energy correspond to those reported in the study of Jadhao et al. [33]. The comparison of activation energy in this study with the literature data [33][34][35][36] is displayed in detail in Table S2 in Supplementary Materials.

Kinetics Model for ICs Pyrolysis
The kinetics model for pyrolysis was analyzed using the Master plot method. As shown in Figure 8a-c, the corresponding g(α), (b) y(α), and z(α) curves were determined. It can be seen that the rising speed of the g(α) curve of the ICs pyrolysis firstly increases, then slows down, and finally accelerates, of which the characteristics are similar to that of random nucleation models A2 and A3. The front part of the y(α) curve of the ICs pyrolysis coincides with that of random nucleation models A2 and A3 on the whole. At the end, the y(α) curve decreases faster than that of the random nucleation model and falls between the reaction order curves F1 and C2. The z(α) curve of the ICs pyrolysis is also most similar to random nucleation models A2 and A3 in the front, and the end of the curve is close to the reaction order model C2. Therefore, it can be concluded that the process of ICs pyrolysis can be described by the hybrid dynamic model consisting of the random nucleation model and reaction order model. Based on the information, the model parameters can be calculated to verify the correctness of the model. et al. [33]. The comparison of activation energy in this study with the literature data [33][34][35][36] is displayed in detail in Table S2 in Supplementary Materials.

Kinetics Model for ICs Pyrolysis
The kinetics model for pyrolysis was analyzed using the Master plot method. As shown in Figure 8a-c, the corresponding g(α), (b) y(α), and z(α) curves were determined. It can be seen that the rising speed of the g(α) curve of the ICs pyrolysis firstly increases, then slows down, and finally accelerates, of which the characteristics are similar to that of random nucleation models A2 and A3. The front part of the y(α) curve of the ICs pyrolysis coincides with that of random nucleation models A2 and A3 on the whole. At the end, the y(α) curve decreases faster than that of the random nucleation model and falls between the reaction order curves F1 and C2. The z(α) curve of the ICs pyrolysis is also most similar to random nucleation models A2 and A3 in the front, and the end of the curve is close to the reaction order model C2. Therefore, it can be concluded that the process of ICs pyrolysis can be described by the hybrid dynamic model consisting of the random nucleation model and reaction order model. Based on the information, the model parameters can be calculated to verify the correctness of the model.

Kinetic Parameters of ICs Pyrolysis
Based on the obtained information above, different parallel reaction models were developed, including the two-step parallel model Am-Cn, three-step parallel model Am-Cn-Cn, three-step parallel model Am-Cn-D1, and four-step parallel model Am-Cn-Cn-Cn. For two, three, and four-step parallel models, the relationship between conversion rates can be described as follows α = w α , w = 1 (12)

Kinetic Parameters of ICs Pyrolysis
Based on the obtained information above, different parallel reaction models were developed, including the two-step parallel model A m -C n , three-step parallel model A m -C n -C n , three-step parallel model A m -C n -D 1 , and four-step parallel model A m -C n -C n -C n . For two, three, and four-step parallel models, the relationship between conversion rates can be described as follows where n is the number of parallel reaction steps. Nonlinear fitting was conducted to obtain the model parameters. The fitting process and variances are shown in Figure 9a-e, and the obtained model parameters are listed in Table 5. On the whole, in the four models, the A m -C n -C n model shows the optimal fitting performance with the root-sum square (RSS) of 0.0021, 0.00831, and 0.01682 at the heating rates of 5 K/min, 10 K/min, and 20 K/min, respectively. The comparison of fitting variances (S 2 ) also confirms this result. Hence, the three-step parallel model A m -C n -C n is the optimal fitting model for the ICs pyrolysis.
where n is the number of parallel reaction steps. Nonlinear fitting was conducted to obtain the model parameters. The fitting process and variances are shown in Figure 9a-e, and the obtained model parameters are listed in Table 5. On the whole, in the four models, the Am-Cn-Cn model shows the optimal fitting performance with the root-sum square (RSS) of 0.0021, 0.00831, and 0.01682 at the heating rates of 5 K/min, 10 K/min, and 20 K/min, respectively. The comparison of fitting variances (S 2 ) also confirms this result. Hence, the three-step parallel model Am-Cn-Cn is the optimal fitting model for the ICs pyrolysis.

Conclusions
The pyrolysis characteristics of ICs were investigated in this study through thermogravimetric analyses. Three different decomposition stages during the pyrolysis of ICs were clearly identified. The pyrolysis reaction speed under the nitrogen condition is faster than that under the air condition. The pyrolysis reaction time under the nitrogen condition is shorter than that under the air condition. It is also worth mentioning that the gas cost of pyrolysis under nitrogen is higher than that under air. The influence of pyrolysis-related factors on the pyrolysis of ICs has been studied by the control variable method. Pyrolysis temperature affects the final composition of reactants and reaction rate and can be controlled to obtain different pyrolysis results. The heating rate has little influence on the reaction rate of the second stage of pyrolysis. The reaction rate of the third stage of pyrolysis increases as the heating rate increases. Different gas atmospheres, or the presence or absence of oxygen, lead to different reaction processes and results of pyrolysis of ICs. Non-isothermal kinetics of ICs has been examined through the Friedman method, CR method, and Kissinger method. The reaction activation energy obtained by the CR method is relatively small, while that obtained by the Friedman method is relatively large. Considered comprehensively, the actual average activation energy of the pyrolysis reaction of the waste memory module should be between 170 and 200 kJ·mol −1 . The process of ICs pyrolysis can be described by the three-step parallel model (A m -C n -C n ). Overall, an efficient and environmentally friendly process technology (pyrolysis) for organic material degradation in waste ICs was determined when pyrolysis gas is considered for use as a circulating fuel. This study provided thermodynamic data and kinetic parameters for further industrial applications.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ma15134460/s1, Figure S1: The linear regression for the extraction of the kinetic parameters from CR method; Table S1: Common reaction kinetics models and mechanism functions.; Table S2: Comparison of activation energy obtained in this study with studies reported in literature.