Abstract

The microstructure evolution of liquid silicon carbide (SiC) during rapid solidification under different pressure values is simulated with the Tersoff potential using molecular dynamics. The structure evolution characteristics of SiC are analyzed by considering the pair distribution function, bond angle distribution, coordination number, and the diagrams of the microstructure during rapid solidification. The results show that the average energy of atoms gradually increases with pressure. When the pressure reaches 100 GPa, the average energy of the atom is greater than the average energy of the atom in the initial liquid state. Under different pressures, the diffusion of atoms tends to remain stable at a temperature of about 3700 K. The application of pressure has a major impact on the arrangement of atoms, except on the third-nearest neighbor, while the impact on the nearest neighbor and the second-nearest neighbor is relatively small. The pressure increases the medium-range order of the system. The coordination numbers of Si and C atoms gradually decrease with the decrease in temperature and increase in pressure. Pressure changes the microstructure of the SiC amorphous system after solidification, and the density can be increased by adjusting the coordination number of atoms. As the pressure increases, the SiC amorphous system exhibits a dense structure with coordination numbers of 4, 5, 6, and 7.

1. Introduction

Silicon carbide (SiC) has excellent physical and mechanical properties, such as low density, high specific strength, and high-temperature resistance. It has great application prospects in fabricating weapon armor, coatings, engine turbine blades, and friction and electronic devices [1]. Extreme electronic devices, i.e., those based on wide-band-gap semiconductors, such as SiC, GaN-based nitride, and diamond, that are operated under complex and harsh conditions—such as at high temperatures and pressures, high strain rates, and even under strong radiation—are attracting substantial attention [2].

Extensive experimental research has been conducted to study the elastic deformation and phase transition of SiC under high pressures. Yoshida et al. [3] found that the 3C-SiC structure transforms from sphalerite to rock salt at about 100 GPa at room temperature (300K). Zhuravlev et al. [4] performed X-ray diffraction (XRD) and Raman experiments at room temperature. They reported that the 3C-SiC lattice remains stable up to 75 GPa, and they also noted a tendency for destabilization above 40 GPa. Daviau and Lee [5] explored the stability of 3C-SiC at high pressures and temperatures in a laser-heated diamond anvil cell. The threshold pressure for the phase transition from zincblende to rock salt structure was found to significantly decrease at high temperatures (1500 K and 1700 K). Daviau and Lee [6] reported that the decomposition that can occur at high temperatures and pressures halted once the RS phase was formed. However, previous studies lacked effective dynamic in situ detection technology and could not identify the process of dynamic deformation of materials in real time. Recently, advanced in situ detection technology has been used worldwide to observe the coexistence of low- and high-pressure phases in a mixed-phase system under laser shock, and the completion of the phase change has been found to require a high pressure of >200 GPa [7]. Nevertheless, the complete phase transition kinetics in real-time in situ conditions have still not been determined.

Molecular dynamics (MD) simulations have unique merits in revealing the atomic structure changes and micro-/nano-scale deformations, and they have been successfully used in studying the mechanical properties [8], surface [9, 10], defects [11], crystal growth [12], and nanotubes and nanowires [1316] of SiC. MD is especially suitable for the study of dynamic physical and mechanical property changes of materials under extreme loading [17, 18]. Recently, Li et al. [19] used large-scale MD to simulate the elastic-plastic and structural phase transitions of SiC under compression loading. Results show that the plastic deformation and phase transition are strongly dependent on the strain rate.

Studying the influence of solidification process control on the preparation and performance of materials can form a more unified understanding of the relationship between the preparation process, structure, and performance. The study of the solidification process of materials has always been an important research topic in the field of materials science. Pressure and temperature are two important parameters that affect the solidification process of materials. The effect of temperature on the solidification process has been extensively studied [20, 21]. In contrast, a few studies have explored the influence of pressure on the solidification process [22], especially the influence of high pressure on the solidification process, because there is no equipment that can directly observe the solidification process under high-pressure conditions. SiC is one of the materials necessary for fabricating extreme electronic devices. The study of the rapid solidification process of SiC under high pressure is necessary. Therefore, the MD simulation method is used to simulate the rapid solidification process of SiC under different pressures and to study the influence of pressure on the dynamic behavior and microstructure of SiC. The evolution of the microstructure provides theoretical guidance for the study of the physical and mechanical properties of SiC under extreme service conditions and can be the basis of the design of related materials and structures.

2. Simulation Conditions and Methods

SiC composed of 8000 atoms is placed in a cubic box with periodic boundary conditions, and the Tersoff potential function [23] and Nosé–Hoover algorithm are used [24, 25] to control the temperature and pressure. The simulation time step is set to 1 fs, and the macro constraint condition is the constant-pressure, constant-temperature (NPT) ensemble. Firstly, the system is heated to a temperature of 4500 K, which is higher than the melting point of SiC, and 40,000 steps are isothermally run to obtain the liquid SiC at equilibrium. Then, at different pressures of 0, 1, 5, 10, 20, 30, 40, 50, and 100 GPa, the system is rapidly cooled from 4500 K to 100 K at a cooling rate of 1 × 1012 K/s, and the atomic configuration is recorded for every 100 K decrement to analyze the microstructure.

3. Results and Analysis

3.1. Average Energy of Atoms and Volume Ratio

Figure 1 shows the average energy of atoms and the volume ratio as functions of temperature during the rapid solidification of the liquid SiC under different pressures.

Figure 1 shows that the volume ratio V/V0 (V and V0 are the liquid and solid volumes of SiC, respectively) of liquid SiC during rapid solidification decreases with the increase in pressure. The average energy of atoms gradually increases with pressure, and the extent of change gradually becomes smaller. Furthermore, when the pressure reaches 100 GPa, the average energy of the atoms becomes greater than the average energy of the atoms in the initial liquid state. The reason for this change is that the interatomic distance in the liquid SiC system gradually reduces with the increase in pressure. Hence, the volume correspondingly decreases, and the value of V/V0 becomes smaller. The distance reduces, resulting in a decrease in the mutual attraction potential energy between the atoms. As the temperature continues to decrease, the kinetic energy of the atoms also decreases. Therefore, during the solidification process, the average energy of the atoms shows a downward trend. When the applied pressure is 100 GPa, since the distance between atoms decreases and becomes smaller than the equilibrium distance r0 between the atoms, the repulsive potential energy increases. Therefore, the average energy of atoms under high pressure of 100 GPa is greater than that of liquids.

3.2. Mean Square Displacement (MSD)

To analyze the influence of different pressures on the dynamics of the rapid solidification process, the relationship between the mean square displacement (MSD) and the temperature during the rapid solidification of liquid SiC was analyzed (Figure 2).

Figure 2 shows that during the rapid solidification of SiC, as the pressure increases, MSD gradually decreases, and the slope of MSD also shows a decreasing trend, indicating that the diffusion coefficient of atoms decreases with the increase in pressure. The application of pressure hinders the movement of atoms. Under different pressures, the diffusion of atoms tends to remain stable at a temperature of about 3700 K.

3.3. Radial Distribution Function (RDF)

RDF is a function that describes the density of particles as a function of the distance from the reference atom. It is an important method to describe the ordering degree of liquid, crystalline, and amorphous structures. The formula of RDF is as follows:

To analyze the influence of pressure on the microstructure during the rapid solidification of SiC, the curve of the radial distribution function (RDF) as a function of temperature under different pressures (Figure 3) is considered.

Figure 3 shows that as the temperature decreases, the first and second peaks of RDF continue to rise and become sharp, indicating that the short-range order in the solidification process is enhanced when different pressures are applied. At the same temperature, the first peak near 1.50 Å gradually decreases with increasing pressure. Two new peaks appear near 2.50 Å and near 3.00 Å, and the peaks increased with increasing pressure, and the position of the peaks gradually decreased with decreasing temperature. Two small peaks near 2.85 Å and 3.20 Å decrease with increasing pressure. To clarify the composition of RDF peaks after solidification, the RDF after solidification of each bond of SiC under different pressures was calculated and analyzed, as shown in Figure 4.

Figures 4(b)–4(d) show that the first peak near 1.5 Å corresponds to the C–C bond, the second peak near 1.9 Å corresponds to the C–Si bond, and the third peak near 2.35 Å corresponds to the Si–Si bond. When the pressure is less than 20 GPa, the peak near 2.5 Å is mainly ascribed to the C–C bond, and the peak decreased as the pressure increased. When the pressure is greater than 20 GPa, a peak rises around 2.5 Å, which reflects that C-Si atoms are replaced at their new positions to construct new short-range order. When the pressure is less than 10 GPa, the peak near 3.0 Å corresponds mainly to the C–C and C–Si bonds, and the peak value increases with pressure. When the pressure is greater than 10 GPa, the second peak corresponding to the Si–Si bond is compressed from 3.2 Å to 3.0 Å. Therefore, this peak corresponded mainly to the Si–Si bond, and its height increases significantly with pressure, while the peak near 3.2 Å gradually decreased. From the above analysis, it can be seen that the application of pressure has a greater impact on the arrangement of atoms, except for the third-nearest neighbor, while the impact on the nearest neighbor and the second-nearest neighbor is relatively small.

The analysis of RDF in our previous study [26] showed that the C atoms distributed around 1.5 Å mainly exhibited a coordination number of 3. After applying pressure, the peak of RDF at 1.5 Å decreases with the increase of pressure, which corresponds to the decrease in the number of three-coordinated C atoms, indicating that C atoms form other coordination structures. The peaks near 2.5 Å and 3.0 Å forward shift, and the increase corresponds to the increase in the number of atoms in this range, indicating that the increase in pressure makes the atomic arrangement of the system more compact and improves the medium-range order of the system.

3.4. Coordination Number (CN)

The coordination number (CN) reflects the arrangement of the nearest-neighbor atoms, and it is a statistical parameter used for studying the arrangement of the nearest-neighbor atoms. Figure 5 shows the variation in the CN of Si and C atoms with temperature under different pressures.

It can be seen from Figure 5(a) that during the rapid solidification of liquid SiC, the CN of all Si atoms is greater than four. When the pressure is higher than 1 GPa, the CN of Si atoms gradually decreases with decreasing temperature and increases with increasing pressure. When the pressure reaches 100 GPa, the CN of Si atoms exceeds six. It can be seen from Figure 5(b) that when the pressure is less than 5 GPa, the CN of C atoms increases with the decrease in temperature. When the pressure is greater than 5 GPa, the CN of C atoms increases with increasing pressure and gradually decreases with decreasing temperature when the pressure is greater than 40 GPa. The CN of C atoms in the solidified SiC is greater than four, indicating that as the pressure increases, the C atoms no longer form a three-coordinated structure.

3.5. Bond Angle Distribution Function (BADF)

The bond angle distribution function (BADF) reflects the positional relationship between the neighboring atoms and the central atom and provides information about the structural evolution of the system during solidification under different pressures. Figure 6 shows the BADF of SiC crystal at different temperatures under different pressures.

Figure 6 shows that the main peaks become higher and narrower as the temperature decreases. When the pressure is 1 GPa, the main peak is at around 120°, and the bond angles are mainly distributed between 70° and 170°. When the pressure is 5 GPa, the main peak at 120° decreases, and another peak starts to appear at 110°. As the pressure increases to 10 GPa, a new peak appears near 60°, the 120° peak further decreases, and the 110° peak’s height increases. Furthermore, when the pressure is 20 GPa, the 110° peak is higher than that of 120°. A new peak appears at 80°, and a shoulder peak appears in the range of 140°to 180°. With a further increase in pressure, three main peaks are formed near 60°, 80°, and 110°. The range of the shoulder peak and the bond angle distribution also increase with pressure, and the bond angle distribution range becomes 40°to 180° at 100 GPa.

To analyze the influence of pressure on the microstructure of liquid SiC after solidification, the bond angle distribution function at 100 K and under different pressures was analyzed (Figure 7).

Figure 7 shows that in the pressure range of 0–10 GPa, the bond angles after the solidification of SiC are mainly distributed around 120° and 110°, indicating that an amorphous random network structure dominated by three-and four-coordinated structures is formed after solidification. With the increase in pressure, the BADF around 120° gradually decreased, indicating that the number of three-coordinated structures decreases with the increase of pressure. The results of the analysis of the change in CN with pressure show that the CN of the C atom increases with pressure, i.e., the three-coordinated structure of C gradually changes to the four-coordinated structure, and the number of the four-coordinated structure increases, which is corresponding to the increasing of BADF around 110°. In the pressure range of 20 GPa to 100 GPa, the BADF around 60°, 80°, 110°, and 140°to 180° increases with pressure, indicating that a new coordination structure is formed in the system.

Figure 8 shows the coordination of C and Si atoms and their bond angle distribution. For example, cn3_C and cn4_C represent the three-coordinated bond angle and the four-coordinated bond angle, respectively, with the C atom at the center, and so on.

It can be seen from Figure 8(a) that the bond angle of the three-coordinated C atom is mainly distributed around 120°, and the peak value gradually decreases with increasing pressure. Figure 8(b) shows that the number of four-coordinated C atoms increases with increasing pressure. It can be seen from Figure 8(c) that the bond angles corresponding to the four-coordinated Si atoms are mainly between 70° and 140°, and the number decreases with increasing pressure. At 100 GPa, the number of four-coordinated Si atoms is almost zero. It can be seen from Figure 8(d) that the bond angles corresponding to the five-coordinated Si atoms are mainly distributed between 70° and 180°, and the number of five-coordinated Si atoms increases with pressure (0–30 GPa). In the range of 30 GPa to 100 GPa, the number of five-coordinated Si atoms decreases with increasing pressure. Figure 8(e) shows that when the pressure is higher than 20 GPa, the bond angles corresponding to the six-coordinated Si atoms are mainly between 40° and 180°, and there is a strong peak near 80° and another peak between 130° and 170°. The number of plateaus increases with pressure. It can be seen from Figure 8(f) that the bond angles corresponding to the seven-coordinated Si atoms are mainly between 40° and 180°. There are two obvious peaks near 60° and 80°, and there is a plateau between 110° and 150°. The peak value increases with pressure.

In summary, the pressure changes the microstructure of the SiC amorphous system after solidification, and the density is increased as a result of the changes in CN atoms. As the pressure increases, the SiC amorphous system exhibits a dense structure with a CN of 4, 5, 6, and 7.

3.6. Microstructure Analysis
3.6.1. Voronoi Polyhedron

The above analysis results show that during the solidification of liquid SiC, the pressure has an impact on the nearest neighbor, second-nearest neighbor, third-nearest neighbor, and on the CN and bond angle distribution. To identify the effect of pressure on the arrangement of atoms, the Voronoi polyhedron method [27] was used to analyze the number of Voronoi polyhedrons during the solidification process. The Voronoi index method is represented by (n3, n4, n5, n6), where n3, n4, n5, and n6 represent the number of triangles, quadrilaterals, pentagons, and hexagons in the polyhedron, respectively. For example, (0 4 0 0) is a tetrahedron composed of four quadrilaterals, (0 6 0 0) is a hexahedron composed of six quadrilaterals, and (2 3 0 0) is a pentahedron composed of two triangles and three quadrilaterals. The minimum CN in the Voronoi index method is four, which corresponds to tetrahedral clusters. This method can accurately describe the crystalline and amorphous structures of semiconductor materials.

Figure 9 shows the large number of clusters that appeared during the solidification process.

Figure 10 shows the variation in the number of Voronoi polyhedrons with temperature under different pressures.

It can be seen from Figure 10 that when the pressure is 1–100 GPa, as the temperature decreases, the (4 0 0 0) structure in the system gradually increases, however, the trend of change in the number of (4 0 0 0) structures with the increase in the pressure is different. When the pressure is 1–20 GPa, the number of (4 0 0 0) structures increases with increasing pressure. Furthermore, for pressure between 20–50 GPa, the number of (4 0 0 0) structures decreases with increasing pressure, and the number of (4 0 0 0) structures at 100 GPa is slightly increased than that of 50 GPa. When the pressure is 20 GPa, the number of (2 3 0 0) structures increases slowly with decreasing temperature, while that of (0 6 0 0) structure decreases with increasing temperature. When the pressure is 30–50 GPa, the number of (2 3 0 0) structures increases with the decrease in temperature and decreases with the increase in pressure. At 30 GPa, the number of (2 3 0 0) structures is the largest. The number of (0 6 0 0) structures increases with decreasing temperature and increases with increasing pressure. At 50 GPa, the number of (2 3 0 0) structures and that of (0 6 0 0) structures in the system are almost the same. When the pressure increases to 100 GPa, the number of (2 3 0 0) structures is greatly reduced, and the number of (0 6 0 0) structures and that of (0 5 2 0) structures exceed the number of (2 3 0 0) structures, which is corresponding to the reduction of the number of five-coordination atomic bond angles. During the high-pressure solidification process, the number of types of Voronoi polyhedrons increases with pressure, resulting in a complex structure with high coordination. The greater the pressure, the more coordinated the Voronoi polyhedrons, and the more complex the structure.

3.6.2. Cross-Section Atomic Distribution

To visually present the influence of pressure on the microstructure of amorphous SiC after solidification, partial cross-sectional views of the amorphous SiC system at 100 K after rapid solidification under different pressures are drawn, as shown in Figure 11: the cross-sectional thickness is 1.5 atomic layers. It can be seen from the figure that the rapid solidification of the system under different pressures results in the formation of an amorphous structure without periodical translation characteristics. With the increase in pressure, the number of bonds between the atoms is different. When the pressure is 0–5 GPa, the atoms in the cross-section are partially bonded to form large rings, and some atoms are far apart and scattered in the system. When the pressure is increased to 10 GPa, the system structure becomes denser with some medium rings, and the number of scattered atoms decreases compared with the situation of 5 GPa. With a further increase in pressure, the number of rings in the system gradually increases, and the cavities formed by the rings become smaller. A three-membered ring appears in the system at 20 GPa. When the pressure is increased to 100 GPa, the density of the system increases significantly, and the number of three-membered rings and that of four-membered rings increase. Except for a few atoms, most of the atoms are bonded, which improves the medium-range order of the system. These results correspond to the results of the analysis based on RDF, in which the influence of pressure on the degree of medium-range order was investigated.

3.6.3. Distribution of Central Atoms with Different CN

According to the results of the analysis of the bond angle distribution function, at high pressure, the SiC amorphous system exhibits a dense structure with CN of 4, 5, 6, and 7. Figure 12 shows the pressure-dependent changes in the number of central atoms with different CN in the amorphous SiC solidified under different pressures.

It can be seen from Figure 12 that in the pressure range of 5–100 GPa, the number of four-coordinated atoms is always the largest. Over the entire pressure range studied, the number of three-coordinated atoms decreases sharply with the increase in pressure, and the magnitude of the decrease decreases with the increase in pressure. Furthermore, the number of five-coordinated atoms varies with pressure in the range of 0–30 GPa, however, at pressures higher than 30 GPa, it decreases with increasing pressure. The number of atoms of CN of six and seven increases with pressure. These results agree with the pressure-dependent variation in CN. Figure 13 depicts the distribution of central atoms with different CN in the system.

3.6.4. Pressure-Dependent Changes in the Main Structure with Central Atom

To analyze the microstructure images of SiC after solidification under different pressures, the temperature-dependent changes in the four-, five-, six-, and seven-coordination structures of Si and the three- and four-coordination structures of C, which are dominant in the rapid solidification process, are analyzed.

AmBn was used to represent the type and the number of atoms in the structure. For example, C1Si3 represents a structure composed of 1 C and 3 Si. cn_C_3 represents a three-coordinated structure with C atoms at the center, respectively. Figure 14 shows the pressure-dependent changes in the main structure with C or Si atoms as the central atom at 100 K.

It can be seen from Figure 14(a) that the number of three-coordinated structures with the C central atom decreases with increasing pressure, which corresponds to the decrease in the peak value of the C–C bond at ∼1.5 Å shown in Figure 4(b). It also corresponds to the decrease in the peak value of the bond angle distribution function of cn3_C at ∼120°, as shown in Figures 8(a). From Figure 14(b), it can be seen that the number of four-coordinated structures with the C atom at the center increases with increasing pressure, which corresponds to the increase in the bond angle distribution function of cn4_C at ∼110°, as shown in Figure 8(b). Furthermore, it can be seen from Figure 14(c) that the number of four-coordinated structures with the Si atom at the center decreases with increasing pressure, which corresponds to the decrease of the Si-Si bond peak at ∼ 2.35 Å in Figure 4(d) and corresponds to the decrease of the BADF of cn4_Si at ∼110° in Figure 8(c). It can be seen from Figure 14(d) that the number of five-coordinated structures with Si atoms at the center increases first and then decreases with increasing pressure. The number of Si4C2 structures, which is the largest number in cn_Si_5, reaches the maximum value at a pressure of 20 GPa, and then gradually decreases with increasing pressure. The number of Si3C3 and Si2C4 reached the maximum value at 30 GPa, and then gradually decrease. Thus, the change in the five-coordination structure corresponds to the decrease in the peak value of the bond angle distribution function of cn5_Si at 70°–120°, as shown in Figure 8(d), and it also corresponds to the change in the number of five-coordination (2 3 0 0) structures with pressure in Figure 10. From Figures 14(e) and 14(f), it can be seen that the number of six-coordinated and seven-coordinated structures with Si atoms at the center increases with increasing pressure, which corresponds to the peak increase of the BADF of cn6_Si and cn7_Si in Figures 8(e) and 8(f) at 60°–180°, and it also corresponds to the change in the number of six-coordinated (0 6 0 0) and seven-coordinated (0 5 2 0) with pressure in Figure 10.

4. Conclusion

From the results of the analysis of the rapid solidification process of liquid SiC under different pressures, the main conclusions obtained are as follows:

During the solidification process, the application of pressure has a major impact on the arrangement of atoms, except on the third-nearest neighbor. Furthermore, the impact on the nearest neighbor and the next nearest neighbor is small. The increase in pressure makes the atomic arrangement of the system more compact and improves the medium-range order of the system. At higher pressures, the system exhibits a complex bond angle distribution, with three main peaks around 60°, 80°, and 110° and a shoulder peak at 140°–180°. The pressure changes the microstructure of the SiC amorphous system after solidification, and the density is improved by adjusting the CN atoms. As the pressure increases, the SiC amorphous system exhibits a dense structure with a CN of 4, 5, 6, and 7.

Data Availability

The underlying data supporting the results of our study are submitted along with the manuscript.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 51671004) and the funds from the Key Laboratory of Materials Simulation and Computing of Anshun University (Grant no. Asxyxkpt201803). The authors thank LetPub (https://letpub.com) for its linguistic assistance during the preparation of this manuscript.