Next Article in Journal
Recent Advances in Tactile Sensory Systems: Mechanisms, Fabrication, and Applications
Previous Article in Journal
Fish Scale for Wearable, Self-Powered TENG
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods

by
Chunhe Jiang
1,
Kejiang Li
2,*,
Zhisheng Bi
2,
Shufang Ma
3,
Jianliang Zhang
2,4,
Bo Liu
5 and
Jiaqi Li
5
1
Technical Support Center for Prevention and Control of Disastrous Accidents in Metal Smelting, University of Science and Technology Beijing, Beijing 100083, China
2
School of Metallurgical and Ecological Engineering, University of Science and Technology Beijing, Beijing 100083, China
3
Aetna Tianlong Tungsten Molybdenum Technology Co., Ltd., No. 11 Fenghui Middle Road, Haidian District, Beijing 100094, China
4
School of Chemical Engineering, The University of Queensland, Saint Lucia, QLD 4072, Australia
5
School of Advanced Engineering, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2024, 14(5), 464; https://doi.org/10.3390/nano14050464
Submission received: 9 February 2024 / Revised: 27 February 2024 / Accepted: 29 February 2024 / Published: 3 March 2024

Abstract

:
Molten slag has different properties depending on its composition. The relationship between its composition, structure, and properties has been the focus of attention in industrial manufacturing processes. This review describes the atomistic scale mechanisms by which oxides of different compositions affect the properties and structure of slag, and depicts the current state of research in the atomic simulation of molten slag. At present, the research on the macroscopic properties of molten slag mainly focuses on viscosity, free-running temperature, melting point, and desulphurization capacity. Regulating the composition has become the most direct and effective way to control slag properties. Analysis of the microevolution mechanism is the fundamental way to grasp the macroscopic properties. The microstructural evolution mechanism, especially at the atomic and nanoscale of molten slag, is reviewed from three aspects: basic oxides, acidic oxides, and amphoteric oxides. The evolution of macroscopic properties is analyzed in depth through the evolution of the atomic structure. Resolution of the macroscopic properties of molten slag by the atomic structure plays a crucial role in the development of fundamental theories of physicochemistry.

1. Introduction

Molten slag is a byproduct formed in the metallurgical process when oxide melts are involved. In the metallurgical process, slag serves a crucial role by facilitating the separation of gangue from liquid iron. Additionally, harmful impurities like phosphorus and sulfur can be effectively eliminated through the careful adjustment of slag compositions. Moreover, slag acts as a protective barrier, shielding the liquid steel from direct contact with elements such as hydrogen, nitrogen, and oxygen. Its multifaceted functions are pivotal in ensuring product quality, optimizing iron recovery rates, maintaining stable smelting operations, and achieving various technical and economic benchmarks [1,2,3,4,5].
Common oxides in molten slag can be divided into three types, including acid oxide, basic oxide, and amphoteric oxide [6,7]. As a representative of acidic oxides, SiO2 is a strong network former, and its atomic structure basically exists as a [SiO4]4− tetrahedron. SiO2 can increase the stability and polymerization degree of slag due to its network formation ability. It requires higher energy to destroy the network structure with [SiO4]4− as the basic unit. In addition, P2O5 is also a common acid oxide in aluminosilicate melt. CaO is a typical basic oxide and the main component for adjusting basicity. Its content has a vital impact on desulfurization and dephosphorization. In addition, increasing the content of CaO can improve the viscosity of slag. The properties of MgO are similar to CaO, and it can also improve fluidity by reducing the viscosity of the slag [8,9]. Al2O3 is a representative amphoteric oxide, and there are obvious differences in its properties in different contents. Excessive Al2O3 generally leads to an increase in slag viscosity, playing the role of acid oxide [10]. When the basic oxide is insufficient, Al2O3 will play the role of a basic oxide to ensure the balanced charge of the system [11]. B and Al elements belong to the same main group; they often have similar physical and chemical properties [12,13]. B2O5 can reduce the viscosity and the melting point of slag. In the micro scale, B3+ ions can increase the polymerization degree of the slag system. However, B2O3 has the ability to reduce slag viscosity [14,15]. Therefore, the adjustment mechanism of B2O3 on the properties of slag is worthy of further fundamental research. TiO2 usually has amphoteric properties [16]. In addition to the above main ingredients, due to the impact of ore resources, it may also contain a small number of other compounds, such as K2O, Na2O, MnO, FeO, etc.
The research on the macroscopic properties of slag is mainly focused on aspects such as melting characteristics, temperature, and fluidity [17,18,19,20]. Over the decades, a considerable amount of research has provided a comprehensive understanding of the high-temperature properties of slag. However, due to the limitations of high-temperature experiments, observing the microscopic structure has become a challenging aspect in the field of experimental research. With the development of modern physical chemistry and computational science, molecular dynamics simulation has become an effective method for researchers to study nanostructures under extreme conditions. Since the 21st century, molecular dynamics simulation methods have gradually been introduced into the field of metallurgy. Research on molten slag is the most widely applied application of molecular dynamics in metallurgy [21]. Through molecular dynamics simulation, researchers have vividly defined the types of oxygen atoms, characterized the network structure, and depicted the evolution of atomic properties. The transition from low-temperature qualitative characterization to high-temperature quantitative research has been achieved. Analyzing the microstructural evolution mechanisms of slag is crucial for controlling its macroscopic properties. Currently, traditional molecular dynamics can simulate studies on oxides such as SiO2, Al2O3, CaO, MgO, TiO2, P2O5, Na2O, K2O, B2O3, FeO, Fe2O3, MnO, etc. Basic microscopic properties, including radial distribution functions; coordination numbers; microscopic structural units; diffusion coefficients; etc., can be accurately computed. However, there is still a significant lack of information on more in-depth electronic structures. Meanwhile, classical molecular dynamics does not have the ability to simulate trace components in slag due to the limitations of potential parameters. Although first-principles molecular dynamics can study all elements, it currently cannot achieve the efficiency of molecular dynamics simulation due to its strong dependence on computational resources. Machine learning molecular dynamics is addressing both accuracy and computational efficiency issues, and is currently being applied to slag research.
This paper aims to review the research progress of molten slag, especially the relationship between the nanostructure and the macroproperties. Clarifying the action mechanism of oxides on the slag is critical to guiding the metallurgical process and the production process of other silicate-related industries.

2. The Atomic Simulation Method

Molecular dynamics simulation is a branch of computational chemistry based on Newtonian classical mechanics. It calculates the trajectories of atoms by computing the interatomic forces, and employs statistical methods to analyze the thermodynamic/dynamic properties of the system. The fundamental steps of molecular dynamics are illustrated in Figure 1. Through the analysis of atomic trajectories, information such as the structure, energy, thermodynamics, and mechanics of the system over a period of time can be obtained. This allows researchers to obtain desired computational results tailored to their specific research directions.
Molecular dynamics methods, based on different approaches to handling interatomic forces, can be categorized into classical molecular dynamics (CMD); ab initio molecular dynamics (AIMD); and machine learning molecular dynamics (MLMD). CMD typically employs a classical potential energy function expressed in a fixed mathematical form to describe the interatomic forces, while AIMD calculates interatomic interactions through quantum chemistry methods. MLMD, on the other hand, utilizes machine learning techniques to train ab initio simulation results to obtain a machine-learned potential. Although these three methods aim to describe interactions differently, they ultimately analyze the physical and chemical properties of a system based on the motion patterns of atoms.

2.1. Classic Molecular Dynamic

In CMD simulations, the choice of the force field is a crucial factor determining the accuracy of the simulation. The form of the force field is often fixed, generally describing the interactions of atoms based on the atomic distance and location. For these different categories of descriptions, classical force fields widely applied in CMD simulations currently include Lennard–Jones (LJ); Born–Mayer–Huggins (BMH); the embedded atom method (EAM) [22]; CHARMM [23]; AMBER [24]; ReaxFF [25]; and others. The emphases of these force fields varies based on the characteristics of different systems.
CMD is the primary method used for studying the atomistic structure of molten slag, as outlined in Table 1. BMH is the main force field employed in slag research, accurately describing the interactions between atoms, and thereby resolving atomic structural relationships. Tersoff [26] appears specifically relevant in the description of SiO2 single crystals. COMPASS [27] is predominantly utilized in molecular simulations using Materials Studio software. MORSE is mainly used in conjunction with the BMH to correct short-range interactions. Current CMD researches on slag have broadly covered common components. However, research on trace elements such as BaO, Cr2O3, etc., is still lacking, primarily due to the absence of force fields or parameters that can accurately describe their interatomic interactions.

2.2. Ab Initio Molecular Dynamics

Compared to CMD, AIMD can compute potential energy and forces using quantum chemistry methods [92,93]. It solves the Schrödinger equation with different kinds of approximations and then determines the trajectories of atomic motion based on Newton’s second law. Since each step in AIMD involves quantum chemical calculations, it heavily relies on computational resources, and the obtained simulation results are typically in the picosecond range. However, AIMD’s advantage lies in its independence from the choice of force field. It can accurately describe interatomic interactions even for less common elements, enabling precise simulations.
With the continuous development of computer technology, AIMD has gradually been applied to slag research, as shown in Table 2. AIMD can not only produce results similar to CMD, but also provide additional information at the electronic scale, such as charge distribution and band structure. Currently, there are relatively few research achievements using AIMD, but it is expected to become an important method in future researches.

2.3. Machine Learning Molecular Dynamics

Machine learning is not a new concept in the field of natural sciences [104], and its primary applications include material structure prediction [105], new material development, and material discovery [106]. However, in the field of molecular dynamics, machine learning is primarily used to obtain accurate interatomic interaction potentials. Machine learning potentials can, on one hand, compensate for the dependence of CMD on parameterized empirical force fields, and on the other hand, address the unacceptable computing resource consumption of large-scale AIMD simulations. MLMD involves training on energy and force information obtained from DFT calculations [107]. The generation of a machine learning potential typically requires three essential steps: (1) a dataset consisting of reference structures and corresponding quantum mechanics (QM) information; (2) a descriptor of the atomic neighbor structure information, which can be input into the machine learning algorithm; (3) a regressor to fit a multilayer perceptron. Based on quantum chemical descriptions and the extension of machine learning, MLMD has become another hot research direction in molecular dynamics simulations.
AIMD, due to its excessive dependence on computational resources, results in smaller system simulations and shorter simulation times compared to CMD. MLMD, on the other hand, allows for easy enlargement of the system, better capturing of accidental events during simulations, and thus the discovery of special phenomena in nature. In slag systems, MLMD has already found practical applications [108]. By thoroughly analyzing large-scale molecular dynamics data, machine learning algorithms can uncover microscopic features hidden in slag, providing a more comprehensive understanding of slag properties. The performance of MLMD is limited by the quality and quantity of training data. Obtaining a high-quality model requires a significant amount of real and diverse molecular dynamics data, which can be a substantial challenge.

3. The Properties of Molten Slag

The properties of slag are closely related to its chemical composition, which generally directly affects its performance—such as the melting point, free running temperature, viscosity thermal conductivity, desulfuration ability, etc. The melting point and free running temperature are two important indicators used to evaluate the slag melting process [109,110]. The free running temperature means that although the slag has not reached the melting point, the viscosity has reached a certain threshold, and the slag can flow freely. It is defined according to the viscosity (Pa·s)-temperature (°C) curves, as shown in Figure 2a. The procedure for calculating entails drawing a tangent line with a slope of -0.02, with the temperature at the tangent point representing the slag’s free-running temperature [109,110], while the melting point refers to the temperature that the object melts completely from solid state to liquid. As shown in Figure 2b, the composition of slag directly influences the variation in slag’s melting point, which is directly related to the alteration of slag’s microstructure. The acid oxides represented by SiO2 can increase the melting point of the slag due to their strong network formation capabilities [111]. However, the effect of basic oxides is exactly the opposite, and within a certain range of compositions, they can often reduce the melting point of slag [111]. In addition, the amphoteric oxides, such as Al2O3, also have the ability to increase the melting point of slag [10,112]. Therefore, to fully understand the changes in macroscopic properties, it is necessary to conduct an in-depth analysis from the microstructure of the slag.
In general, the operating temperatures of metallurgical processes exceed the melting point of slag, making the fluidity of slag a key parameter of concern. During metallurgical processes, viscosity changes are typically controlled by adjusting the composition of slag. However, as high-quality mineral resources become increasingly scarce, the variety of ores entering the furnace has become richer, resulting in complex compositions of slag. In recent years, there has been extensive research on the impact of changes in different oxide components on slag viscosity, as shown in Figure 2c. Generally, basic oxides can improve slag fluidity by reducing the structural complexity. Previous research has already demonstrated that basic oxides, such as CaO [113], MgO [114], FeO [115], Na2O [116], etc., possess this characteristic. However, if the content is too high, there is a possibility of increasing the melting point, consequently leading to a decrease in the fluidity of slag. Xing et al. has found that the addition of BaO to slag leads to the formation of more complex structure, thereby increasing the viscosity. Acidic oxides generally only play the role of increasing structural complexity in slag. Therefore, SiO2 and P2O5 [117,118,119,120] will directly cause the viscosity of the slag to increase, while the role of amphoteric oxides will have different effects on the fluidity of slag according to changes in relative composition and temperature. Figure 2c shows that the influence of Al2O3 on the viscosity of slag is obviously different under different composition conditions [121]. The roles played by amphoteric oxides such as Cr2O3 [122] are similar to those of Al2O3. Therefore, to clarify the mechanism of action of amphoteric oxides, detailed analysis is required based on different conditions.
During the metallurgical process, changes in slag properties are crucial to the adjustment of impurity elements in molten metal. The desulfurization process is the most intuitive process of regulating slag properties. Sulfur is a harmful impurity in steel products. Its presence in steel materials will seriously affect their macroscopic properties. Therefore, regulating the sulfur content in materials is crucial to the production of high-quality steel materials [123]. In thermodynamics, sulfur capacity is usually used to characterize the desulfurization ability of slag, as shown in Figure 2d. Public researches showed that basic oxide components, such as CaO [124], MgO [125], MnO [126], and BaO [127] can improve the desulfurization ability of slag, while components such as SiO2 [128], Al2O3 [129], and TiO2 [130] have the opposite effect. The sulfur capacity of slag is also directly related to its microstructure evolution [125,126,127,129,130].
Figure 2. (a) The viscosity–temperature curve of slag ([110,111]); (b) the evolution of melting point with different compositions ([10,112,113]); (c) the variation of viscosity with different content ([114,115,116,117,118,119,120,121,122,123,131]); (d) the evolution of sulfur capacity with different contents ([126,127,128,130,131]).
Figure 2. (a) The viscosity–temperature curve of slag ([110,111]); (b) the evolution of melting point with different compositions ([10,112,113]); (c) the variation of viscosity with different content ([114,115,116,117,118,119,120,121,122,123,131]); (d) the evolution of sulfur capacity with different contents ([126,127,128,130,131]).
Nanomaterials 14 00464 g002
During the metallurgical process, regulating the composition of slag is the most intuitive method to adjust the evolution of the macroscopic properties of slag. In current research, most of the composition changes have been investigated in detail. However, high-temperature experiments are dangerous and complicated, and current research methods cannot detect the microstructural variation of slag under high-temperature conditions in real time. Therefore, the microscopic analysis of slag is still lacking. Developments in science have proven that microstructure is the fundamental factor affecting the evolution of macroscopic properties. Moreover, the ingredient adjustment method summarized through experience lacks scientificity. Therefore, analyzing the microstructure evolution mechanism of slag with different components is crucial to understanding the evolution of properties.

4. The Evolution of Molten Slag’s Microstructure

Molten slag is mainly composed of aluminosilicate, consisting of [SiO4]4− and [AlO4]5− tetrahedron. In the nanoscale network structure, the existence of oxygen atoms is critical to the properties of the aluminosilicate system. According to the role of oxygen, as shown in Figure 3, oxygen atoms can be divided into four types: free oxygen (FO), non-bridge oxygen (NBO), bridge oxygen (BO), and tricluster oxygen (TO). The network structure of molten slag is mainly connected by BO [53,64]. BO refers to the oxygen atom that connects the two tetrahedron structures, typically Si-O-Si, Si-O-Al, or Al-O-Al. Generally, three tetrahedron structures are connected by TOs, which is usually less in the system, but also plays an important role in improving the degree of polymerization (DOP) of the system. The oxygen atom occupied by only one tetrahedron structure is NBO, which exists only at the edge of the tetrahedron structure. The higher the NBO content, the lower the DOP of the system. In addition, there will be some FOs in the system. As usual, the ratio of BO to NBO can be used to represent the DOP of the system. Moreover, the number of BO on the tetrahedron can also represent the DOP of the system, such as QSin and QAln [42] (the proportion represented by SiO4 or AlO4 tetrahedra with non-bridging oxygen atoms). DOP represents the complexity of the microstructure of the slag. The more complicated the structure, the more difficult it is to change its original properties, and the more stable. DOP is directly related to the properties, such as viscosity, melting temperature, thermal conductivity, and desulfurization capacity, of molten slag.

4.1. The Effect of Acid Oxides

SiO2 and P2O5 are common acidic oxides that make up the tetrahedron structure in the slag. Usually, SiO2 is the main component in metallurgical molten slag, and P2O5 only exists in small amounts. In the silicate system, the bond length of the Si-O is about 1.67Å [101]. Si mainly forms a tetrahedron structure with 4 O atoms. The bond angle of O-Si-O is generally maintained at 108°, as shown in Figure 4a [28]. Si is the atom with the smallest diffusion coefficient in the slag, and the worst movement. Li et al. [30] investigated the influence of SiO2 content on the local structure and viscosity of slag. It was found that the main structural unit of the formation of the Si-O network was a [SiO4]4− tetrahedron. As shown in Figure 4b,c, with the decrease of SiO2 the BO content decreases, while the TO increases, in which the amount of O(Al, Al, Al) increases and the amount of O(Al, Al, Si) and O(Al, Al, Si) decreases. It indicates that [AlO4]5− tetrahedra are more likely to be connected by TO, while [SiO4]4− tetrahedra prefer to be connected by BO. Due to the stable structure of a [SiO4]4− tetrahedron, the chemical composition has little effect on the bond length of Si-O. With the decrease of SiO2 content, the diffusion coefficients of Si, Al, Ca, and O atoms increase significantly, which is mainly because the reduction of the Si-O bond weakens the stability of the network structure. Therefore, the viscosity of the system is reduced, which proves that the acid oxide has an obvious effect on the fluidity of the system. Sajid et al. investigated the electronic structure of SiO2-based slag, as shown in Figure 4d. Oxygen atoms typically gain more electrons, while Si and Al atoms primarily lose electrons, resulting in the formation of Si-O and Al-O bonds. Due to the greater electron loss by silicon atoms, the Si-O bond is the most stable chemical bond in the furnace slag system [96].
P2O5 is an acidic oxide with a lesser content in the molten slag [76], and it can also act as the network former in the slag, as shown in Figure 5a [77]. The bond length of P-O is about 1.53Å [48], which is slightly shorter than the Si-O. The stability of the P-O bond is relatively weaker than the Si-O bond. Fan et al. [72] found that the bond length of P-O will decrease as the P2O5 content increases. P atoms are more inclined to form an almost perfect [PO4]3− tetrahedron structure. The bond angle of O-P-O is about 109.3°, slightly larger than that of O-Si-O, as shown in Figure 5b. Du et al. [75] explored the microstructure of different P2O5-based binary melt. The basic unit of the P2O5-based binary melt is a [PO4]3− tetrahedron, which is very close to the perfect tetrahedron structure. With the increase in P2O5 content, Q0, Q1, and Q2 will gradually convert to Q3 structure, resulting in the more complex network structure of each phosphate system. Qian et al. [98] found that the role of P2O5 in aluminosilicate varies with the content. When the P2O5 content is minimal, most of the O atoms in [PO4]3− are NBOs, which can play a role in repolymerizing a [SiO4]4− based network structure and promote the increasing of the Si-O-Si proportion. With the further increase in the P2O5 content, the number of BO in the [PO4]3− tetrahedron structure continue to increase, and P atoms gradually participate in the construction of the network structure.
B2O3 is an excellent surface-active substance which can reduce the viscosity of slag. As a flux, B2O3 can be substituted for CaF2 to reduce the free-running temperature [15]. Structurally, B2O3-containing slag is more likely to form [BO3]3− trigonal planar structures, which affects the structure and properties of the slag. Wang et al. [132]. demonstrated that B2O3 has a significant inhibitory effect on CaO-SiO2 activity. Within a certain content, increasing the B2O3 content can reduce the strength of the network structure, and decrease the slag viscosity and apparent activation energy. Lai et al. [133] found that B2O3 leads to the formation of [SiO4]4− structures that increase the polymerization of a CaO-Al2O3-SiO2-TiO2-B2O3 system. The B and O atoms mainly form two structures, [BO3]3− trigonal planar and [BO4]5− tetrahedron. Zhang et al.’s research [83] found that the addition of B2O3 results in the formation of a B-O network structure, which becomes randomly mixed with the initial Si-O network structure. However, the properties exhibited by B2O3 vary significantly under different slag conditions. In low basicity conditions, B2O3 reduces the slag’s polymerization, whereas in high basicity conditions, the opposite occurs. Bi et al. [35] further revealed that B2O3 can exhibit amphoteric properties, like Al2O3. In the SiO2-CaO-B2O3 system, it transforms from acidic to basic when the B2O3 content is greater than 21%. When the content of acidic oxides is too high, B2O3 will show basic properties, but the role played is completely different from Ca2+. The Ca ions reduce the system’s polymerization by destroying the BO, while B3+ tends to capture O2+ to form chemical bonds. In turn, it changes from decreasing the molten slag microstructure polymerization, to reducing the slag complexity [37]. However, generally molten slag has a low content of B2O3, which can only show the nature of acidic oxides.

4.2. The Effect of Basic Oxides

CaO and MgO are the most common basic oxides in metallurgical molten slags, and are generally used as network modifiers to adjust the macroscopic properties. Microscopically, CaO and MgO dissociate into Ca2+, Mg2+, and O2−, as shown in Figure 6a. The cations will interact with the network structure formed by the acidic and amphoteric oxides in the slag, thus affecting the macroscopic properties of the slag. Jiang et al. [40,41] comparatively analyzed the mechanism of microscopic effects of CaO and MgO on the slag properties. The results show that the changes in the relative contents of CaO and MgO do not affect the interatomic chemical bond length in the slag. However, Ca2+ and Mg2+ ions are playing the role of network modifiers, which will damage the intrinsic network structure, as shown in Figure 6b. The increase in CaO content will lead to more conversion of the TO into BO and NBO, which in turn reduces the structural complexity of the system. MgO has the same function of reducing the structural complexity of the slag. Nevertheless, when CaO and MgO are only present alone and in the same amount, MgO is significantly more capable of destroying the microscopic network structure. Ca2+ ions have a stronger charge compensation ability than Mg2+ ions, so Ca2+ is more effective in maintaining the stability of the network structure. Meanwhile, both ions do not significantly affect the bond angles of the network structure. Zhang et al. [100] utilized AIMD simulations to discover that only a small amount of electron accumulation occurs around Ca atoms. Therefore, Ca-O is the ionic bond in the molten slag. Gao et al. [114] further demonstrated by Fourier transform infrared (FTIR) spectroscopy that the complex network structure of slag can be depolymerized into simpler units by CaO or MgO.
FeO is a special component of molten slag, which both dissolves into the slag to influence the evolution of slag properties, and changes the slag-iron interaction through reduction process. Macroscopically, FeO has the ability to reduce slag viscosity. Fan et al. [57] investigated the effect of FeO on the structure and transport properties of the TiO2-SiO2 slag system and found that FeO has the ability to destroy the Ti-O-Ti structure. The Fe has the effect of reducing the complexity of the slag by replacing one of the Ti ions in the structure, which in turn improves the fluidity, as shown in Figure 7a. Ma et al. [58] further investigated the effect of FeO on the SiO-Al2O3-CaO slag system, and further resolved that FeO can improve the slag by promoting the conversion of TO and BO to NBO, as shown in Figure 7b. However, CaO destroys the network structure faster than FeO under the same condition [45]. Fe atoms can exist in the slag in both Fe2+ (FeO) and Fe3+ (Fe2O3) ions. It was found that both Fe2+ and Fe3+ have the ability to reduce the BO, and the effect of Fe3+ is relatively stronger. Fe3+ can form the FeO4 structure by replacing Si or Al atoms, which in turn has the effect of reducing the degree of polymerization, as shown in Figure 7c. Fe2+, on the other hand, tends to exert only a similar effect as Ca2+ and Mg2+ [59].
With the shortage of advantageous mineral resources, the utilization of multicomponent minerals gradually leads to composition diversification in metallurgical molten slag. MnO is increasingly becoming a more common component. As expressed previously, MnO has the ability to improve slag viscosity and enhance desulphurization. Ma et al. [71] further investigated the mechanism of the effect of MnO on the slag under high temperature conditions. It was found that the change in MnO content does not directly affect the bond length of the original chemical bond, but it leads to a decrease in the coordination number of Al and O atoms. MnO can promote the transformation of BO and TO to NBO and FO, and the higher the basicity, the more obvious the effect of MnO. Also, the self-diffusion coefficient of each atom increases with the rise of MnO content, which in turn improves the slag viscosity. Yuan et al. [51] revealed that Mn ions tend to interact with BO structural units, thereby modifying the network structure. This behavior is notably different from Mg ions. Mn can serve as charge compensated ions, but its effectiveness is weaker compared to Mg ions. He et al. [99] comparatively analyzed the effect of CaO and MnO on the desulphurization properties of slag using the ab initio MD method. The results showed that MnO is also an ionic oxide, with the same ionic bond between Mn and O atoms, but weaker than Ca-O. A small number of Mn-Mn clusters will be present in the MnO-SiO2 melt, while MnO leads to a significant alteration of the desulphurization mechanism. As shown in Figure 8, S will form a stable bonding structure with Mn atoms when incorporated into the liquid MnO-SiO2, and the Si-S bond is hardly found in the S-doped MnO-SiO2 slag. MnO will result in the decrease of BO and the increase in NBO during desulphurization. MnO will increase the sulphur capacity by increasing the oxygen ion activity and decreasing the activity coefficient of S ions [99].
Na is a typical alkali element, which has been treated as a hazardous element in metallurgical processes. Although its content in molten slag has always been low, its action has been of wide interest. Jiang et al. [64] found that Na2O can increase the coordination number of Si-O as well as Al-O in aluminosilicate, with a more pronounced effect on Al-O. Generally, both Si-O and Al-O structures are most stable when they are tetra-coordinated, so the stability of the system decreases successively. Na2O has the ability to promote the conversion of BO to NBO, TO, and FO, thereby increasing the diffusion ability and leading to a reduction in viscosity. It indicates that Na2O also plays the role of network modifier in the molten slag system, but research has shown that its effect is significantly weaker than that of CaO [44]. In the phosphate glasses system, Na2O also plays a role in depolymerizing the network structure of the system and can significantly increase the ratio of NBO oxygen [63]. The number of chemically durable P-O-Fe bridges decreases and they are replaced by P-O-Na. The iron cations will be separated from PO4 units and form iron oxides clusters. K2O and Na2O are typical alkali oxides, and since Na and K belong to the same main group of elements, they are often generally treated as having the same properties. However, they play different roles in terms of microstructural mechanisms [34]. Li et al. [61] comparatively investigated the effects of K2O and Na2O on the fluidity of coke ash and found that they play opposite roles. K2O will lead to an increase in viscosity, while Na2O does just the opposite, and yet the effect of both on the microstructure is lacking minimally. Significant differences were however observed in the location of these ions in the oxygen bonding networks. While Na+ ions were preferentially located in the bridging/non-bridging oxygen networks, K+ ions tended to be present in various oxygen tri-clusters. However, when both basic alkali are present in the slag at the same time, it leads to a sharp BO and NBO increase and a sharp TO decrease [62].

4.3. The Effect of Amphoteric Oxide

Amphoteric oxides in different environments can show mutual transformation of both acidic and basic properties. Al2O3 is a common component of metallurgical molten slag, which in general plays the same role as SiO2 and is a constituent unit of the microscopic network structure. However, given the specificity of amphoteric oxides, the active role of Al2O3 is equally a focus of research. Chen et al. [11] investigated the effect of Al2O3 on the structure and properties of a CaO-SiO2-Al2O3 system using molecular dynamics simulations. Al2O3 is considered to be acidic and a network former in the environment with sufficient basic cations. In contrast, Al2O3 becomes basic with an insufficient basic cations environment, thus providing oxygen atoms and balancing charge. The increase in Al2O3 concentration leads to a relative shortage of Ca atoms and therefore to an increase in the content of TO and five-coordinated Al, resulting in structural instability. Bi et al. [35] similarly found that Al2O3 exhibits a significant amphoteric transformation based on the variation in content. When Al2O3 is relatively low, it mainly exhibits as acidic oxide. In this case, the BO content and polymerization degree of the system increase with the increase in Al2O3 content. When the content exceeds the critical point, the effect of basic oxides is exerted. It will make the BO content of the system decrease and the growth rate of polymerization degree decrease. Xuan et al. [36] analyzed the effects of the amphoteric transformation process of Al2O3 on the viscosity and microstructure, based on the results of previous research. As Al2O3 increases, the coordination number of the basic cations increases, changing from a network modifier to a charge compensation ion. The transition is determined by the charge balancing ability of the slag. The addition of Al2O3 leads to the replacement of [SiO4]4− by [AlO4]5−, as shown in Figure 9, with a decrease in the contents of Si-O-Si and Si-O-M, and an increase in Si-O-Al, Al-O-Al, and TO. Therefore, the network strength increases first and then decreases, which is consistent with the tendency of viscosity.
TiO2 can exhibit amphoteric properties, but there is currently a debate over whether it is truly an amphoteric oxide. Ti ions primarily exist in three forms in the molten slag: [TiO4], [TiO5], and [TiO6]. Most Ti-O-Ti bond angles are 100°, while a small proportion of Ti-O-Ti bond angles vary from 126.67° to 131.64° [57]. Yao et al. [68] found that Si–O–Ti linkages are more favorable than either Si–O–Si or Ti–O–Ti linkages in a CaO-SiO2-TiO2 system, as Si-O-Ti linkages account for the overwhelming majority of BO groups. Fan and Chen et al. [38,69,70] found that Ti4+ ions act as an amphoteric component in slag, i.e., it behaves as a network-forming component with a tetrahedral structure (acidic), or a network-disrupting component with an octahedral structure (basic). Basicity is the key parameter affecting the amphoteric transition of Ti4+. As the basicity decreases, the [TiO4]4− tetrahedral structure evolves to [TiO6]8− octahedral structure, implying that the addition of SiO2 favors the formation of [TiO6]8− octahedra. Ti-O stretching and Ti-O-Al vibrations can form short-range ordered simple network structures that reduce the effect of SiO2 on viscosity. TiO2 can reduce system viscosity by replacing SiO2, but also consumes BO to form short-range ordered [TiO4]4− tetrahedral networks.
In summary, it is shown that amphoteric oxides are capable of both atomic structure network formation and destruction, and the corresponding slag properties will be affected by the amphoteric transition.

5. Conclusions and Prospective

In recent years, the development of computers and molecular dynamics has greatly expanded the level and depth of metallurgical process research. The combination of metallurgical physicochemical and atomic simulation can significantly increase the safety, efficiency, and scientificity of the metallurgical experimental process. It is conducive to improving the theory of physicochemistry at the atomic scale and promoting the development of basic research on the properties of oxides in molten slag. Molten slag and its derivatives have been widely valued and studied in the fields of metallurgy, materials, coal chemical industry and earth science. At present, research on the mechanisms of the influence of different natural oxides on the structure of metallurgical slag mainly focuses on the following aspects:
(1)
Acidic oxides, represented by SiO2 and P2O5, mainly form [SiO4]4− tetrahedra and [PO4]5− tetrahedra in molten slag. The [SiO4]4− tetrahedra and [PO4]5− tetrahedra are interconnected with other polyhedral structures, leading to the formation of complex chain-like or network structures. Acidic oxides often increase the viscosity, liquidus temperature, and enthalpy of molten slag by increasing the degree of polymerization of the its atomic structure.
(2)
Basic oxides, represented by CaO and MgO, dissociate into free cations and free oxygen ions in the molten slag. The free cations of basic oxides undergo dynamic equilibrium vibrations around polyhedral structures to compensate for the negative charge overflow. Meanwhile, cations and oxygen ions destroy the network structures of the system during migration and interaction, thereby improving the properties of the melt.
(3)
Amphoteric oxides, represented by Al2O3, possess both acidic and basic oxide properties. When the content of acidic oxides is high, they tend to exhibit the property of basic oxides destroying the network structure. When the content of acidic oxides is low, they generally exhibit the property of forming acidic oxide network structures.
CMD simulations have demonstrated significant capabilities in the field of molten slag research. However, the current composition of metallurgical molten slag has gradually become more complex, and the existing potential functions of CMD can no longer meet the needs of all compositional research. Furthermore, due to the limitations of current force fields, it is difficult for CMD to reproduce the chemical bond properties between different types of atoms in the molten slag, including ionic bonds, covalent bonds, etc., which is also a limitation of CMD. Moreover, current CMD methods can only qualitatively reproduce the regularities of changes in macroscopic properties, and the calculation of macroscopic properties still lacks accuracy. There is an urgent need to develop more accurate potential functions to supplement existing research in depth, aiming to enhance the research of molten slag properties in various scenarios. AIMD is a powerful method for addressing simulations with complex compositions in molten slag, but its computational efficiency limits its applicability to large-scale simulations. However, AIMD can accurately depict the interactions between atoms, characterize the types of chemical bonds between atoms, and elucidate the mechanisms of electron interactions. It serves as an effective complement to CMD methods. Therefore, the synergistic use of CMD and AIMD can address most of the issues in slag simulation. MLMD holds great potential in balancing computational efficiency and accuracy. By ensuring that the dataset used for training the potential functions is diverse enough, the detailed exploration of the structures and properties of molten slag with the coexistence of multiple complex components can be expected. In the future, CMD, AIMD, and MLMD methods will each play crucial roles in different emphasized areas.

Funding

This research received the support of the Young Elite Scientist Sponsorship Program by CAST (YESS20210090).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

Shufang Ma was employed by the company Aetna Tianlong Tungsten Molybdenum Technology Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Wu, T.; Wang, Q.; Yu, C.F.; He, S.P. Structural and viscosity properties of CaO–SiO2–Al2O3–FeO slags based on molecular dynamic simulation. J. Non Cryst. Solids 2016, 450, 23–31. [Google Scholar] [CrossRef]
  2. Ozbay, E.; Erdemir, M.; Durmus, H.l. Utilization and efficiency of ground granulated blast furnace slag on concrete properties–A review. Constr. Build. Mater. 2016, 105, 423–434. [Google Scholar] [CrossRef]
  3. Taylor, R.; Richardson, I.; Brydson, R. Composition and microstructure of 20-year-old ordinary Portland cement ground granulated blast furnace slag blends containing 0 to 100% slag. Cem. Concr. Res. 2010, 40, 971–983. [Google Scholar] [CrossRef]
  4. Karamanova, E.; Avdeev, G.; Karamanov, A. Ceramics from blast furnace slag, kaolin and quartz. J. Eur. Ceram. Soc. 2011, 31, 989–998. [Google Scholar] [CrossRef]
  5. Greeley, R.; Schneid, B.D. Magma generation on Mars: Amounts, rates, and comparisons with Earth, Moon, and Venus. Science 1991, 254, 996–998. [Google Scholar] [CrossRef] [PubMed]
  6. Yan, Z.-M.; Lv, X.-W.; Li, Z.-S. Physicochemical properties and structure of titania-containing metallurgical slags: A review. J. Iron Steel Res. Int. 2022, 29, 187–206. [Google Scholar] [CrossRef]
  7. Waseda, Y.; Toguri, J.M. The Structure and Properties of Oxide Melts: Application of Basic Science to Metallurgical Processing; World Scientific: London, UK, 1998. [Google Scholar]
  8. Liu, W.; Pang, Z.; Wang, J.; Zuo, H.; Xue, Q. Investigation of viscosity and structure of CaO–SiO2–MgO–Al2O3–BaO–B2O3 slag melt. Ceram. Int. 2022, 48, 17123–17130. [Google Scholar] [CrossRef]
  9. Qiu, G.; Miao, D.; Wei, X.; Bai, C.; Li, X. Effect of MgO/Al2O3 and CaO/SiO2 on the Metallurgical Properties of CaO–SiO2–Al2O3–MgO–TiO2 Slag. J. Non Cryst. Solids 2022, 585, 121545. [Google Scholar] [CrossRef]
  10. Yao, L.; Ren, S.; Wang, X.; Liu, Q.; Dong, L.; Yang, J.; Liu, J. Effect of Al2O3, MgO, and CaO/SiO2 on viscosity of high alumina blast furnace slag. Steel Res. Int. 2016, 87, 241–249. [Google Scholar] [CrossRef]
  11. Chen, Y.; Pan, W.; Jia, B.; Wang, Q.; Zhang, X.; Wang, Q.; He, S. Effects of the amphoteric behavior of Al2O3 on the structure and properties of CaO–SiO2–Al2O3 melts by molecular dynamics. J. Non Cryst. Solids 2021, 552, 120435. [Google Scholar] [CrossRef]
  12. Pang, Z.; Lv, X.; Jiang, Y.; Ling, J.; Yan, Z. Blast furnace ironmaking process with super-high TiO2 in the slag: Viscosity and melting properties of the slag. Metall. Mater. Trans. B 2020, 51, 722–731. [Google Scholar] [CrossRef]
  13. Yan, H.-j.; Liu, L.; Zhuang, J.-c.; Zhou, P.; Zhou, C.Q. Molecular Dynamics Simulation of Carbon Effect on the Thermal Physical Properties of the Molten Iron. ISIJ Int. 2018, 59, 221–226. [Google Scholar] [CrossRef]
  14. Yan, W.; Chen, W.; Yang, Y.; Lippold, C.; McLean, A. Evaluation of B2O3 as replacement for CaF2 in CaO–Al2O3 based mould flux. Ironmak. Steelmak. 2016, 43, 316–323. [Google Scholar] [CrossRef]
  15. Kim, G.H.; Sohn, I. Effect of CaF2, B2O3 and the CaO/SiO2 mass ratio on the viscosity and structure of B2O3 containing calcium silicate based melts. J. Am. Ceram. Soc. 2019, 102, 6575–6590. [Google Scholar] [CrossRef]
  16. Liao, J.; Li, J.; Wang, X.; Zhang, Z. Influence of TiO2 and basicity on viscosity of Ti bearing slag. Ironmak. Steelmak. 2012, 39, 133–139. [Google Scholar] [CrossRef]
  17. Park, H.; Park, J.Y.; Kim, G.H.; Sohn, I. Effect of TiO2 on the viscosity and slag structure in blast furnace type slags. Steel Res. Int. 2012, 83, 150–156. [Google Scholar] [CrossRef]
  18. Shankar, A.; Görnerup, M.; Lahiri, A.; Seetharaman, S. Estimation of viscosity for blast furnace type slags. Ironmak. Steelmak. 2007, 34, 477–481. [Google Scholar] [CrossRef]
  19. Lv, X.; Lv, X.; Qiu, J.; Liu, M. Viscosity and Structure Evolution of Slag in Ferronickel Smelting Process from Laterite. J. Min. Metall. Sect. B Metall. 2017, 53, 147–154. [Google Scholar] [CrossRef]
  20. Min, D.J.; Tsukihashi, F. Recent advances in understanding physical properties of metallurgical slags. Met. Mater. Int. 2017, 23, 1–19. [Google Scholar] [CrossRef]
  21. Zhou, C.; Li, J.; Wang, S.; Zhao, J.; Ai, L.; Chen, Q.; Chen, Q.; Zhao, D. Development of Molecular Dynamics and Research Progress in the Study of Slag. Materials 2023, 16, 5373. [Google Scholar] [CrossRef]
  22. Finnis, M.; Sinclair, J. A simple empirical N-body potential for transition metals. Philos. Mag. A 1984, 50, 45–55. [Google Scholar] [CrossRef]
  23. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar] [CrossRef] [PubMed]
  24. Cornell, W.D.; Cieplak, P.; Bayly, C.I.; Gould, I.R.; Merz, K.M.; Ferguson, D.M.; Spellmeyer, D.C.; Fox, T.; Caldwell, J.W.; Kollman, P.A. A second generation force field for the simulation of proteins, nucleic acids, and organic molecules. J. Am. Chem. Soc. 1995, 117, 5179–5197. [Google Scholar] [CrossRef]
  25. Chenoweth, K.; Van Duin, A.C.; Goddard, W.A. ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. J. Phys. Chem. A 2008, 112, 1040–1053. [Google Scholar] [CrossRef] [PubMed]
  26. Tersoff, J. New empirical approach for the structure and energy of covalent systems. Phys. Rev. B 1988, 37, 6991. [Google Scholar] [CrossRef] [PubMed]
  27. Sun, H. COMPASS: An ab initio force-field optimized for condensed-phase applications overview with details on alkane and benzene compounds. J. Phys. Chem. B 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  28. Hoang, V.V. Molecular dynamics simulation of amorphous SiO2 nanoparticles. J. Phys. Chem. B 2007, 111, 12649–12656. [Google Scholar] [CrossRef] [PubMed]
  29. Nath, S.D. Study of the effect of sizes on the structural properties of SiO2 glass by molecular dynamics simulations. J. Non Cryst. Solids 2013, 376, 50–60. [Google Scholar] [CrossRef]
  30. Li, K.; Bouhadja, M.; Khanna, R.; Zhang, J.; Liu, Z.; Zhang, Y.; Yang, T.; Sahajwalla, V.; Yang, Y.; Barati, M. Influence of SiO2 reduction on the local structural order and fluidity of molten coke ash in the high temperature zone of a blast furnace. A molecular dynamics simulation investigation. Fuel 2016, 186, 561–570. [Google Scholar] [CrossRef]
  31. Liu, Y.; Bai, C.; Lv, X.; Wei, R. Molecular dynamics simulation on the influence of Al2O3 on the slag structure at 1873 K. Mater. Today Proc. 2015, 2, S453–S459. [Google Scholar] [CrossRef]
  32. Wu, T.; Wang, Q.; Yao, T.; He, S. Molecular dynamics simulations of the structural properties of Al2O3-based binary systems. J. Non Cryst. Solids 2016, 435, 17–26. [Google Scholar] [CrossRef]
  33. Yuelin, Q.; Hao, L.; Yanhua, Y. Structure evolution of blast furnace slag with high Al2O3 Content and 5 mass% TiO2 via molecular dynamics simulation and fourier transform infrared spectroscopy. Metall. Res. Technol. 2018, 115, 113. [Google Scholar] [CrossRef]
  34. Zhao, Y.; Du, J.; Qiao, X.; Cao, X.; Zhang, C.; Xu, G.; Liu, Y.; Peng, S.; Han, G. Ionic self-diffusion of Na2O–Al2O3–SiO2 glasses from molecular dynamics simulations. J. Non Cryst. Solids 2020, 527, 119734. [Google Scholar] [CrossRef]
  35. Bi, Z.; Li, K.; Jiang, C.; Zhang, J.; Ma, S. Effects of amphoteric oxide (Al2O3 and B2O3) on the structure and properties of SiO2–CaO melts by molecular dynamics simulation. J. Non Cryst. Solids 2021, 559, 120687. [Google Scholar] [CrossRef]
  36. Xuan, W.; Zhang, Y. Exploration of the amphoteric transition of Al2O3 on melt structure and viscosity of silicate slag. Ceram. Int. 2023, 49, 25815–25822. [Google Scholar] [CrossRef]
  37. Bi, Z.; Li, K.; Jiang, C.; Zhang, J.; Ma, S.; Sun, M.; Wang, Z.; Li, H. Performance and transition mechanism from acidity to basicity of amphoteric oxides (Al2O3 and B2O3) in SiO2–CaO–Al2O3–B2O3 system: A molecular dynamics study. Ceram. Int. 2021, 47, 12252–12260. [Google Scholar] [CrossRef]
  38. Chen, Y.; Yang, J.; Zhang, X.; Wang, Q.; Wang, Q.; He, S. Amphoteric behavior of component and microstructure feature on CaO-Al2O3-TiO2 ternary melt by molecular dynamics simulation. Comput. Mater. Sci. 2022, 205, 111223. [Google Scholar] [CrossRef]
  39. Liu, J.; Yu, Y.; Kong, W.; Wang, Q.; He, Z.; Yang, X. Effect of Al2O3 on the molten slag microstructure information by molecular dynamics simulation and experimental analysis. J. Non Cryst. Solids 2022, 578, 121354. [Google Scholar] [CrossRef]
  40. Jiang, C.; Li, K.; Zhang, J.; Qin, Q.; Liu, Z.; Liang, W.; Sun, M.; Wang, Z. The effect of CaO (MgO) on the structure and properties of aluminosilicate system by molecular dynamics simulation. J. Mol. Liq. 2018, 268, 762–769. [Google Scholar] [CrossRef]
  41. Jiang, C.; Li, K.; Zhang, J.; Liu, Z.; Niu, L.; Liang, W.; Sun, M.; Ma, H.; Wang, Z. The effect of CaO and MgO on the structure and properties of coal ash in the blast furnace: A molecular dynamics simulation and thermodynamic calculation. Chem. Eng. Sci. 2019, 210, 115226. [Google Scholar] [CrossRef]
  42. Zhao, H.; Li, J.; Yang, S.; Liu, J.; Liu, W. Molecular dynamics study of structural properties of refining slag with various CaO/Al2O3 ratios. Minerals 2021, 11, 398. [Google Scholar] [CrossRef]
  43. Fan, H.; Wang, R.; Duan, H.; Chen, D.; Xu, Z. Structural and transport properties of TiO2–SiO2–MgO–CaO system through molecular dynamics simulations. J. Mol. Liq. 2021, 325, 115226. [Google Scholar] [CrossRef]
  44. Gao, L.; Liu, X.; Bai, J.; Kong, L.; Bai, Z.; Li, W. Structure and flow properties of coal ash slag using ring statistics and molecular dynamics simulation: Role of CaO/Na2O in SiO2–Al2O3–CaO–Na2O. Chem. Eng. Sci. 2021, 231, 116285. [Google Scholar] [CrossRef]
  45. Ma, S.; Li, K.; Zhang, J.; Jiang, C.; Bi, Z.; Sun, M.; Wang, Z.; Li, H. The effects of CaO and FeO on the structure and properties of aluminosilicate system: A molecular dynamics study. J. Mol. Liq. 2021, 325, 115106. [Google Scholar] [CrossRef]
  46. Gu, C.; Lyu, Z.; Bao, Y. Size-dependent dissolution behavior of CaO in the CaO–SiO2–FeO slag system: A molecular dynamics study. J. Mol. Liq. 2022, 367, 120488. [Google Scholar] [CrossRef]
  47. Diao, J.; Ke, Z.; Jiang, L.; Zhang, Z.; Zhang, T.; Xie, B. Structural properties of molten CaO–SiO2–P2O5–FeO system. High. Temp. Mater. Process. 2017, 36, 871–876. [Google Scholar] [CrossRef]
  48. Liu, Y.; Lv, X.; Li, B.; Bai, C. Relationship between structure and viscosity of CaO–SiO2–MgO–30.00 wt% Al2O3 slag by molecular dynamics simulation with FT-IR and Raman spectroscopy. Ironmak. Steelmak. 2018, 45, 492–501. [Google Scholar] [CrossRef]
  49. Sadaf, S.; Wu, T.; Zhong, L.; Liao, Z.Y.; Wang, H.C. Effect of basicity on the structure, viscosity and crystallization of CaO–SiO2-B2O3 based mold fluxes. Metals 2020, 10, 1240. [Google Scholar] [CrossRef]
  50. Yuan, H.; Wang, Z.; Zhang, Y.; Wang, C. Roles of MnO and MgO on structural and thermophysical properties of SiO2–MnO–MgO-B2O3 welding Fluxes: A molecular dynamics study. J. Mol. Liq. 2023, 386, 122501. [Google Scholar] [CrossRef]
  51. Jiang, C.; Li, K.; Zhang, J.; Qin, Q.; Liu, Z.; Sun, M.; Wang, Z.; Liang, W. Effect of MgO/Al2O3 ratio on the structure and properties of blast furnace slags: A molecular dynamics simulation. J. Non Cryst. Solids 2018, 502, 76–82. [Google Scholar] [CrossRef]
  52. Jiang, C.; Li, K.; Zhang, J.; Qin, Q.; Liu, Z.; Liang, W.; Sun, M.; Wang, Z. Molecular dynamics simulation on the effect of MgO/Al2O3 ratio on structure and properties of blast furnace slag under different basicity conditions. Metall. Mater. Trans. B 2019, 50, 367–375. [Google Scholar] [CrossRef]
  53. Kong, W.-G.; Liu, J.-H.; Yu, Y.-W.; Hou, X.-M.; He, Z.-J. Effect of w (MgO)/w (Al2O3) ratio and basicity on microstructure and metallurgical properties of blast furnace slag. J. Iron Steel Res. Int. 2021, 28, 1223–1232. [Google Scholar] [CrossRef]
  54. Fan, H.; Chen, D.; Liu, P.; Duan, H.; Huang, Y.; Long, M.; Liu, T. Structural and transport properties of FeO–TiO2 system through molecular dynamics simulations. J. Non Cryst. Solids 2018, 493, 57–64. [Google Scholar] [CrossRef]
  55. Siakati, C.; Macchieraldo, R.; Kirchner, B.; Tielens, F.; Peys, A.; Seveno, D.; Pontikes, Y. Unraveling the nano-structure of a glassy CaO-FeO-SiO2 slag by molecular dynamics simulations. J. Non Cryst. Solids 2020, 528, 119771. [Google Scholar] [CrossRef]
  56. Fan, H.; Wang, R.; Xu, Z.; Duan, H.; Chen, D. Structural and transport properties of FeO–TiO2–SiO2 systems: Insights from molecular dynamics simulations. J. Non Cryst. Solids 2021, 571, 121049. [Google Scholar] [CrossRef]
  57. Ma, S.; Li, K.; Zhang, J.; Jiang, C.; Sun, M.; Li, H.; Wang, Z.; Bi, Z. Structural Characteristics of CaO–SiO2–Al2O3–FeO Slag with Various FeO Contents Based on Molecular Dynamics Simulations. JOM 2021, 73, 1637–1645. [Google Scholar] [CrossRef]
  58. Ma, S.; Ren, S.; Li, K.; Zhang, J.; Jiang, C.; Bi, Z.; Sun, M. The Effects of FeO and Fe2O3 on the Structure and Properties of Aluminosilicate System: A Molecular Dynamics Study. JOM 2022, 74, 4162–4173. [Google Scholar] [CrossRef]
  59. Yang, H.; Liu, Y.; Wan, X.; Zhang, T.-a.; Lin, S.; Wang, K. The Effect of FeO on Transport Properties of Dephosphorization Slag from Microstructure: A Molecular Dynamics Simulation Study. Trans. Indian. Inst. Met. 2023, 76, 3165–3173. [Google Scholar] [CrossRef]
  60. Li, K.; Khanna, R.; Bouhadja, M.; Zhang, J.; Liu, Z.; Su, B.; Yang, T.; Sahajwalla, V.; Singh, C.V.; Barati, M. A molecular dynamic simulation on the factors influencing the fluidity of molten coke ash during alkalization with K2O and Na2O. Chem. Eng. J. 2017, 313, 1184–1193. [Google Scholar] [CrossRef]
  61. Li, K.; Khanna, R.; Zhang, J.; Bouhadja, M.; Sun, M.; Barati, M.; Liu, Z.; Singh, C.V. Molecular dynamics investigation on coke ash behavior in the high-temperature zones of a blast furnace: Influence of alkalis. Energy Fuels 2017, 31, 13466–13474. [Google Scholar] [CrossRef]
  62. Goj, P.; Stoch, P. Molecular dynamics simulations of P2O5-Fe2O3-FeO-Na2O glasses. J. Non Cryst. Solids 2018, 500, 70–77. [Google Scholar] [CrossRef]
  63. Jiang, C.; Zhang, H.; Xiong, Z.; Chen, S.; Li, K.; Zhang, J.; Liang, W.; Sun, M.; Wang, Z.; Wang, L. Molecular dynamics investigations on the effect of Na2O on the structure and properties of blast furnace slag under different basicity conditions. J. Mol. Liq. 2020, 299, 112195. [Google Scholar] [CrossRef]
  64. Zhang, X.; Liu, C.; Jiang, M. Effect of Na Ions on Melt Structure and Viscosity of CaO–SiO2–Na2O by Molecular Dynamics Simulations. ISIJ Int. 2021, 61, 1389–1395. [Google Scholar] [CrossRef]
  65. Zhang, S.; Zhang, X.; Bai, C.; Wen, L.; Lv, X. Effect of TiO2 content on the structure of CaO–SiO2–TiO2 system by molecular dynamics simulation. ISIJ Int. 2013, 53, 1131–1137. [Google Scholar] [CrossRef]
  66. Zhang, S.; Zhang, X.; Peng, H.; Wen, L.; Qiu, G.; Hu, M.; Bai, C. Structure analysis of CaO–SiO2–Al2O3–TiO2 slag by molecular dynamics simulation and FT-IR spectroscopy. ISIJ Int. 2014, 54, 734–742. [Google Scholar] [CrossRef]
  67. Yao, T.-H.; He, S.-P.; Wu, T.; Wang, Q. Molecular dynamics simulations of microstructural properties of CaO–SiO2–TiO2 fluorine-free slag systems. Ironmak. Steelmak. 2017, 44, 551–558. [Google Scholar] [CrossRef]
  68. Fan, H.; Zhu, Y.; Xu, Z.; Wang, R. Structure of the TiO2−MgO–Al2O3 system: Insights from molecular dynamics simulations. J. Non Cryst. Solids 2022, 584, 121482. [Google Scholar] [CrossRef]
  69. Kim, Y.; Nam, C.; Kim, S.; Jeon, H. Rheological properties and structure of molten FeO-TiO2-B2O3 ilmenite smelting slag. J. Non Cryst. Solids 2021, 552, 120308. [Google Scholar] [CrossRef]
  70. Ma, S.; Li, K.; Zhang, J.; Jiang, C.; Bi, Z.; Sun, M.; Wang, Z. Effect of MnO content on slag structure and properties under different basicity conditions: A molecular dynamics study. J. Mol. Liq. 2021, 336, 116304. [Google Scholar] [CrossRef]
  71. Fan, G.; Diao, J.; Jiang, L.; Zhang, Z.; Xie, B. Molecular dynamics analysis of the microstructure of the CaO–P2O5–SiO2 slag system with varying P2O5/SiO2 ratios. Mater. Trans. 2015, 56, 655–660. [Google Scholar] [CrossRef]
  72. Saiko, I.; Saetova, N.; Raskovalov, A.; Il’ina, E.; Molchanova, N.; Kadyrova, N. Hopping conductivity in V2O5–P2O5 glasses: Experiment and non-constant force field molecular dynamics. Solid. State Ion. 2020, 345, 115180. [Google Scholar] [CrossRef]
  73. Nguyen, M.A.; Nguyen, V.H. Structural Properties of Liquid CaO–SiO2–P2O5 System. VNU J. Sci. Math. Phys. 2023, 39. [Google Scholar] [CrossRef]
  74. Du, Y.; Yuan, Y.; Li, L.; Long, M.; Duan, H.; Chen, D. Insights into structure and properties of P2O5-based binary systems through molecular dynamics simulations. J. Mol. Liq. 2021, 339, 116818. [Google Scholar] [CrossRef]
  75. Zhao, Y.; Du, J.; Cao, X.; Zhang, C.; Xu, G.; Qiao, X.; Liu, Y.; Peng, S.; Han, G. A modified random network model for P2O5–Na2O–Al2O3–SiO2 glass studied by molecular dynamics simulations. RSC Adv. 2021, 11, 7025–7036. [Google Scholar] [CrossRef]
  76. Sun, H.; Yang, J.; Zhang, R.; Xu, L. Insight into the structure and transport properties of CaO-SiO2-P2O5 system during the phosphorus enrichment process: A molecular dynamics simulation. J. Non Cryst. Solids 2024, 627, 122818. [Google Scholar] [CrossRef]
  77. Zhang, Z.; Xie, B.; Zhou, W.; Diao, J.; Li, H.-Y. Structural characterization of FeO–SiO2–V2O3 slags using molecular dynamics simulations and FT-IR spectroscopy. ISIJ Int. 2016, 56, 828–834. [Google Scholar] [CrossRef]
  78. Saetova, N.; Raskovalov, A.; Antonov, B.; Denisova, T.; Zhuravlev, N. Structural features of Li2O–V2O5–B2O3 glasses: Experiment and molecular dynamics simulation. J. Non Cryst. Solids 2020, 545, 120253. [Google Scholar] [CrossRef]
  79. Feng, X.; Yao, W.; Li, J. Effect of B2O3 on the structure of CaO–Al2O3–B2O3 ternary melts: A molecular dynamics simulation. J. Non Cryst. Solids 2021, 574, 121141. [Google Scholar] [CrossRef]
  80. Bi, Z.; Li, K.; Jiang, C.; Zhang, J.; Ma, S.; Sun, M.; Wang, Z.; Li, H. Effects of B2O3 on the structure and properties of blast furnace slag by molecular dynamics simulation. J. Non Cryst. Solids 2021, 551, 120412. [Google Scholar] [CrossRef]
  81. Fan, H.; Wang, R.; Xu, Z.; Duan, H.; Chen, D. The effect of B2O3 on the structure and properties of titanium slag melt by molecular dynamics simulations. J. Mater. Res. Technol. 2021, 15, 1046–1058. [Google Scholar] [CrossRef]
  82. Zhang, X.; Liu, C.; Jiang, M. Effect of B2O3 on the Melt Structure and Viscosity of CaO–SiO2 System. Steel Res. Int. 2022, 93, 2100520. [Google Scholar] [CrossRef]
  83. Yang, R.; Zhang, Y.; Zu, Q.; Huang, S.; Zhang, L.; Deng, L.; Zeng, H. Molecular dynamics simulations study on structure and properties of CaO–MgO–B2O3–Al2O3–SiO2 glasses with different B2O3/MgO. J. Non Cryst. Solids 2023, 616, 122458. [Google Scholar] [CrossRef]
  84. Asada, T.; Yamada, Y.; Ito, K. The estimation of structural properties for molten CaO–CaF2–SiO2 system by molecular dynamics simulations. ISIJ Int. 2008, 48, 120–122. [Google Scholar] [CrossRef]
  85. Zhang, X.; Liu, C.; Jiang, M. Effect of Fluorine on Melt Structure for CaO–SiO2–CaF2 and CaO–Al2O3–CaF2 by Molecular Dynamics Simulations. ISIJ Int. 2020, 60, 2176–2182. [Google Scholar] [CrossRef]
  86. Zhang, X.; Liu, C.; Jiang, M.; Sun, J.; Zheng, X. The Relationship Between Composition and Structure: The Influence of Different Methods of Adding CaF2 on the Melt Structure of CaO–SiO2 Slags by Molecular Dynamics Simulations. J. Sustain. Metall. 2022, 8, 1269–1275. [Google Scholar] [CrossRef]
  87. Jia, B.; Li, M.; Yan, X.; Wang, Q.; He, S. Structure investigation of CaO–SiO2–Al2O3–Li2O by molecular dynamics simulation and Raman spectroscopy. J. Non Cryst. Solids 2019, 526, 119695. [Google Scholar] [CrossRef]
  88. Jiao, S.; Guo, P.; Chen, F.; Min, Y.; Liu, C. Molecular Dynamics Analysis of the Structural Behavior of Aluminum Ion in the Slag of CaO–SiO2–Al2O3–Li2O System. Steel Res. Int. 2024, 95, 2300310. [Google Scholar] [CrossRef]
  89. Wang, Z.; Li, Z.; Zhong, M.; Li, Z.; Wang, C. Elucidating the effect of Al2O3/SiO2 mass ratio upon SiO2-MnO-CaF2-Al2O3-based welding fluxes: Structural analysis and thermodynamic evaluation. J. Non Cryst. Solids 2023, 601, 122071. [Google Scholar] [CrossRef]
  90. Zheng, K.; Zhang, Z.; Yang, F.; Sridhar, S. Molecular dynamics study of the structural properties of calcium aluminosilicate slags with varying Al2O3/SiO2 ratios. ISIJ Int. 2012, 52, 342–349. [Google Scholar] [CrossRef]
  91. Xu, B.; Cao, Y.; Wang, Z.; Du, P.; Long, Y. Molecular Dynamics Study on the Effect of SiO2/Al2O3 Mass Ratio on the Structural Properties and Viscosity of Molten Fused Red Mud. Minerals 2022, 12, 925. [Google Scholar] [CrossRef]
  92. Marx, D.; Hutter, J. Ab initio molecular dynamics: Theory and implementation. Mod. Methods Algorithms Quantum Chem. 2000, 1, 141. [Google Scholar]
  93. Iftimie, R.; Minary, P.; Tuckerman, M.E. Ab initio molecular dynamics: Concepts, recent developments, and future trends. Proc. Natl. Acad. Sci. USA 2005, 102, 6654–6659. [Google Scholar] [CrossRef]
  94. Scherer, C.; Schmid, F.; Letz, M.; Horbach, J. Structure and dynamics of B2O3 melts and glasses: From ab initio to classical molecular dynamics simulations. Comput. Mater. Sci. 2019, 159, 73–85. [Google Scholar] [CrossRef]
  95. Sajid, M.; Bai, C.; Yu, W.; You, Z.; Tan, M.; Aamir, M. First principle study of electronic structural and physical properties of CaOSiO2Al2O3 ternary slag system by using Ab Initio molecular and lattice dynamics. J. Non Cryst. Solids 2020, 550, 120384. [Google Scholar] [CrossRef]
  96. He, X.; Ma, S.; Wang, L.; Dong, H.; Chou, K. An Ab Initio Molecular Dynamics Simulation of Liquid FeO–SiO2 Silicate System with Sulfur Dissolving. Metall. Mater. Trans. B 2021, 52, 3346–3353. [Google Scholar] [CrossRef]
  97. Qian, Y.; Song, B.; Jin, J.; Prayogo, G.I.; Utimula, K.; Nakano, K.; Maezono, R.; Hongo, K.; Zhao, G. Ab initio molecular dynamics simulation of structural and elastic properties of SiO2–P2O5–Al2O3–Na2O glass. J. Am. Ceram. Soc. 2022, 105, 6604–6615. [Google Scholar] [CrossRef]
  98. He, X.; Ma, S.; Wang, L.; Dong, H.; Chou, K. Comparison of desulfurization mechanism in liquid CaO–SiO2 and MnO–SiO2: An ab initio molecular dynamics simulation. J. Alloys Compd. 2022, 896, 163008. [Google Scholar] [CrossRef]
  99. Zhang, C.; Kong, Y.-Q.; Wu, T.; Lei, J.; Wang, H.-C. First-principles study on microstructure of CaO–Al2O3–B2O3 slag. J. Mol. Liq. 2022, 368, 120738. [Google Scholar] [CrossRef]
  100. Sarnthein, J.; Pasquarello, A.; Car, R. Structural and electronic properties of liquid and amorphous si o 2: An ab initio molecular dynamics study. Phys. Rev. Lett. 1995, 74, 4682. [Google Scholar] [CrossRef]
  101. Jiang, C.; Li, K.; Barati, M.; Guo, P.; Danaei, A.; Liang, W.; Zhang, J. The interaction mechanism between molten SiO2–Al2O3–CaO slag and graphite with different crystal orientations: Experiment and Ab initio molecular dynamics simulation. Ceram. Int. 2023, 49, 8295–8301. [Google Scholar] [CrossRef]
  102. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. [Google Scholar] [CrossRef]
  103. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef]
  104. Gasteiger, J.; Zupan, J. Neural networks in chemistry. Angew. Chem. Int. Ed. Engl. 1993, 32, 503–527. [Google Scholar] [CrossRef]
  105. Senior, A.W.; Evans, R.; Jumper, J.; Kirkpatrick, J.; Sifre, L.; Green, T.; Qin, C.; Žídek, A.; Nelson, A.W.; Bridgland, A. Improved protein structure prediction using potentials from deep learning. Nature 2020, 577, 706–710. [Google Scholar] [CrossRef]
  106. Dral, P.O. Quantum chemistry in the age of machine learning. J. Phys. Chem. Lett. 2020, 11, 2336–2347. [Google Scholar] [CrossRef]
  107. Zhang, L.; Han, J.; Wang, H.; Car, R.; Weinan, E. Deep potential molecular dynamics: A scalable model with the accuracy of quantum mechanics. Phys. Rev. Lett. 2018, 120, 143001. [Google Scholar] [CrossRef] [PubMed]
  108. Sun, Y.; Tan, M.; Li, T.; Li, J.; Shang, B. Study on the structural properties of refining slags by molecular dynamics with deep learning potential. J. Mol. Liq. 2022, 353, 118787. [Google Scholar] [CrossRef]
  109. Jiang, T.; Wang, S.; Guo, Y.F.; Chen, F.; Zheng, F.Q. Effects of Basicity and MgO in Slag on the Behaviors of Smelting Vanadium Titanomagnetite in the Direct Reduction-Electric Furnace Process. Metals 2016, 6, 107. [Google Scholar] [CrossRef]
  110. Wei Zhang, K.Z. jiating Rao, Weiguo Fu and Mansheng Chu, Discussion on standardization of Melt-property Temperature of slag. J. IRon Steel Res. 2011, 23, 16–19. [Google Scholar]
  111. Liang, H.; Chu, M.; Feng, C.; Tang, J.; Liu, Z.; Wang, W. Optimisation study and affecting mechanism of CaO/SiO2 and MgO on viscous behaviours of titanium-bearing blast furnace slag. Ironmak. Steelmak. 2018, 47, 106–117. [Google Scholar] [CrossRef]
  112. Gao, Y.; Bian, L.; Liang, Z. Influence of B2O3 and TiO2 on viscosity of titanium-bearing blast furnace slag. Steel Res. Int. 2015, 86, 386–390. [Google Scholar] [CrossRef]
  113. Gao, Y.-M.; Wang, S.-B.; Hong, C.; Ma, X.-J.; Yang, F. Effects of basicity and MgO content on the viscosity of the SiO2–CaO–MgO–9 wt% Al2O3 slag system. Int. J. Miner. Metall. Mater. 2014, 21, 353–362. [Google Scholar] [CrossRef]
  114. Kim, H.; Kim, W.H.; Sohn, I.; Min, D.J. The effect of MgO on the viscosity of the CaO–SiO2–20 wt% Al2O3–MgO slag system. Steel Res. Int. 2010, 81, 261–264. [Google Scholar] [CrossRef]
  115. Jiao, K.X.; Zhang, J.L.; Wang, Z.Y.; Chen, C.L.; Liu, Y.X. Effect of TiO2 and FeO on the viscosity and structure of blast furnace primary slags. Steel Res. Int. 2017, 88, 1600296. [Google Scholar] [CrossRef]
  116. Xu, R.Z.; Zhang, J.L.; Han, W.X.; Chang, Z.Y.; Jiao, K.X. Effect of BaO and Na2O on the viscosity and structure of blast furnace slag. Ironmak. Steelmak. 2018, 47, 168–172. [Google Scholar] [CrossRef]
  117. Ye, G.; Yang, J.; Zhang, R.; Yang, W.; Sun, H. Behavior of phosphorus enrichment in dephosphorization slag at low temperature and low basicity. Int. J. Miner. Metall. Mater. 2021, 28, 66–75. [Google Scholar] [CrossRef]
  118. Schwitalla, D.H.; Guhl, S.; Körner, J.; Laabs, M.; Bai, J.; Meyer, B. Meta-study on the effect of P2O5 on single phase slag viscosity and the effect of P2O5 induced liquid phase immiscibility on dispersion viscosity. Fuel 2021, 305, 121501. [Google Scholar] [CrossRef]
  119. Xia, Y.; Li, J.; Fan, D.; Hou, G. Effects of interfacial oxygen potential and slag phase changing during slag formation process on dephosphorization behavior. ISIJ Int. 2019, 59, 1519–1526. [Google Scholar] [CrossRef]
  120. Sun, Y.; Zheng, K.; Liao, J.; Wang, X.; Zhang, Z. Effect of P2O5 addition on the viscosity and structure of titanium bearing blast furnace slags. ISIJ Int. 2014, 54, 1491–1497. [Google Scholar] [CrossRef]
  121. Zhang, S.; Li, H.; Ran, M.; Jiang, Z.; Zheng, L.; Feng, H.; Yu, J.; Dai, Y. Effect of Al2O3 on Viscosity and Refining Ability of High Basicity Slag for Heat-resistant Austenitic Stainless Steel. ISIJ Int. 2022, 62, 2207–2216. [Google Scholar] [CrossRef]
  122. Xu, R.Z.; Zhang, J.L.; Wang, Z.Y.; Jiao, K.X. Influence of Cr2O3 and B2O3 on viscosity and structure of high alumina slag. Steel Res. Int. 2017, 88, 1600241. [Google Scholar] [CrossRef]
  123. Cherry, J.; Davies, H.; Mehmood, S.; Lavery, N.; Brown, S.; Sienz, J. Investigation into the effect of process parameters on microstructural and physical properties of 316L stainless steel parts by selective laser melting. Int. J. Adv. Manuf. Technol. 2015, 76, 869–879. [Google Scholar] [CrossRef]
  124. Zhang, J.; Lv, X.; Yan, Z.; Qin, Y.; Bai, C. Desulphurisation ability of blast furnace slag containing high Al2O3 and 5 mass% TiO2 at 1773 K. Ironmak. Steelmak. 2016, 43, 378–384. [Google Scholar] [CrossRef]
  125. Shen, F.; Jiang, X.; Wu, G.; Wei, G.; Li, X.; Shen, Y. Proper MgO addition in blast furnace operation. ISIJ Int. 2006, 46, 65–69. [Google Scholar] [CrossRef]
  126. Liu, W.; Zuo, H. Effect of MnO and substituting CaO with BaO on the desulfurization ability of blast furnace slag. Metall. Mater. Trans. B 2021, 52, 2275–2282. [Google Scholar] [CrossRef]
  127. Ju, J.; Liu, W.; Xing, X.; Wang, J.; An, J. Sulphide capacity of CaO–SiO2–8% MgO–Al2O3–BaO slags ranging from 0% to 5% in BaO. Metall. Res. Technol. 2019, 116, 106. [Google Scholar] [CrossRef]
  128. Park, J.H.; Park, G.-H. Sulfide capacity of CaO–SiO2–MnO–Al2O3–MgO Slags at 1873 K. ISIJ Int. 2012, 52, 764–769. [Google Scholar] [CrossRef]
  129. Talapaneni, T.; Yedla, N.; Sarkar, S. Study on desulfurization capacity of high alumina blast furnace slag at 1773 K using slag-metal equilibrium technique. Metall. Res. Technol. 2018, 115, 502. [Google Scholar] [CrossRef]
  130. Liu, W.; Xing, X.; Zuo, H. Effect of TiO2 on viscosity and sulfide capacity of blast furnace slag containing barium. ISIJ Int. 2020, 60, 1886–1891. [Google Scholar] [CrossRef]
  131. Xu, C.-y.; Wang, C.; Xu, R.-Z.; Zhang, J.-l.; Jiao, K.-X. Effect of Al2O3 on the viscosity of CaO–SiO2–Al2O3–MgO–Cr2O3 slags. Int. J. Miner. Metall. Mater. 2021, 28, 797–803. [Google Scholar] [CrossRef]
  132. Wang, H.; Li, G.; Zhu, X.; Zhao, Y. Research on Activities of CaO–SiO2–B2O3(Al2O3, Fe2O3) Slag System. Appl. Mech. Mater. 2012, 217–219, 492–495. [Google Scholar]
  133. Lai, F.; Yao, W.; Li, J. Effect of B2O3 on structure of CaO–Al2O3–SiO2–TiO2–B2O3 glassy systems. ISIJ Int. 2020, 60, 1596–1601. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of molecular dynamics simulation.
Figure 1. Schematic diagram of molecular dynamics simulation.
Nanomaterials 14 00464 g001
Figure 3. Atomic configuration of oxygen structural unit (Reproduced with permission from [64]. Elsevier, 2020).
Figure 3. Atomic configuration of oxygen structural unit (Reproduced with permission from [64]. Elsevier, 2020).
Nanomaterials 14 00464 g003
Figure 4. (a) Distribution of bond angle of O-Si-O with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (b) evolution of non-bridging oxygen, bridge oxygen and tricluster oxygen with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (c) fraction of different types of O triclusters amongst total triclusters with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (d) electron density difference calculated at 1773 K (Reproduced with permission from [96]. Elsevier, 2020).
Figure 4. (a) Distribution of bond angle of O-Si-O with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (b) evolution of non-bridging oxygen, bridge oxygen and tricluster oxygen with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (c) fraction of different types of O triclusters amongst total triclusters with change of SiO2 content (Reproduced with permission from [30]. Elsevier, 2016); (d) electron density difference calculated at 1773 K (Reproduced with permission from [96]. Elsevier, 2020).
Nanomaterials 14 00464 g004
Figure 5. (a) Silicate network and phosphate network structure in molten slag system (Reproduced with permission from [77]. Elsevier, 2024); (b) schematic diagram of average bond angles of O–Si–O, and O–P–O (Reproduced with permission from [72]. J-STAGE, 2015).
Figure 5. (a) Silicate network and phosphate network structure in molten slag system (Reproduced with permission from [77]. Elsevier, 2024); (b) schematic diagram of average bond angles of O–Si–O, and O–P–O (Reproduced with permission from [72]. J-STAGE, 2015).
Nanomaterials 14 00464 g005
Figure 6. (a) Structure of molten slag (purple tetrahedron represents [SiO4], orange tetrahedron represents [AlO4], green and yellow balls represent Ca and Mg cations, respectively) (Reproduced with permission from [40]. Elsevier, 2019); (b) depolymerization of network structure with MgO and CaO (Reproduced with permission from [52]. Elsevier, 2018).
Figure 6. (a) Structure of molten slag (purple tetrahedron represents [SiO4], orange tetrahedron represents [AlO4], green and yellow balls represent Ca and Mg cations, respectively) (Reproduced with permission from [40]. Elsevier, 2019); (b) depolymerization of network structure with MgO and CaO (Reproduced with permission from [52]. Elsevier, 2018).
Nanomaterials 14 00464 g006
Figure 7. (a) The oxygen connection state in FeO-based slag with varying FeO content (Reproduced with permission from [57]. Elsevier, 2021); (b) the evolution of bridging oxygen, non-bridging oxygen, tricluster oxygen and free oxygen with the FeO content (Reproduced with permission from [58]. Springer Nature, 2021); (c) a schematic diagram of local atomic structure, which is also used to explain the variation of BO, NBO, FO, and TO in the system (Reproduced with permission from [59]. Springer, 2022).
Figure 7. (a) The oxygen connection state in FeO-based slag with varying FeO content (Reproduced with permission from [57]. Elsevier, 2021); (b) the evolution of bridging oxygen, non-bridging oxygen, tricluster oxygen and free oxygen with the FeO content (Reproduced with permission from [58]. Springer Nature, 2021); (c) a schematic diagram of local atomic structure, which is also used to explain the variation of BO, NBO, FO, and TO in the system (Reproduced with permission from [59]. Springer, 2022).
Nanomaterials 14 00464 g007
Figure 8. Schematic diagram of desulfurization mechanism of liquid MnO-SiO2 system (Reproduced with permission from [99]. Elsevier, 2022).
Figure 8. Schematic diagram of desulfurization mechanism of liquid MnO-SiO2 system (Reproduced with permission from [99]. Elsevier, 2022).
Nanomaterials 14 00464 g008
Figure 9. Amphoteric transformation of Al2O3 in the melt (Reproduced with permission from [36]. Elsevier, 2023).
Figure 9. Amphoteric transformation of Al2O3 in the melt (Reproduced with permission from [36]. Elsevier, 2023).
Nanomaterials 14 00464 g009
Table 1. Representative research of CMD in molten slag.
Table 1. Representative research of CMD in molten slag.
VariableSystemPotentialTemperatureReference
SiO2SiO2BMH and MORSE 1700–1800 K[28]
SiO2Tersoff 300 K and 5000 K[29]
SiO2-CaO-Al2O3BMH 2223 K[30]
Al2O3SiO2-CaO-Al2O3-MgOBMH 1873 K[31]
Na2O-Al2O3, K2O-Al2O3, MgO-Al2O3, and CaO-Al2O3,BMH 2473 K[32]
SiO2-CaO-Al2O3-MgO-TiO2BMH 1773 K[33]
Na2O-Al2O3-SiO2BMH 4000 K, 3000 K, 2500 K, 2000 K, 1500 K, 1000 K, 750 K, 700 K, 650 K, 600 K, 550 K, 500 K, 30 K[34]
SiO2-CaO-Al2O3BMH 1773 K[35]
SiO2-CaO-Al2O3-MgO-FeOBMH 1873 K[36]
SiO2-CaO-Al2O3BMH 1873 K[11]
SiO2-CaO-Al2O3-B2O3BMH 1773 K[37]
CaO-Al2O3-TiO2BMH 1873 K[38]
SiO2-CaO-Al2O3-MgOCOMPASS 1773 K[39]
CaOSiO2-CaO-Al2O3BMH 1773 K[40]
SiO2-CaO-Al2O3-MgOBMH 1773 K[41]
SiO2-CaO-Al2O3-MgOBMH 1800 K[42]
TiO2-SiO2-MgO-CaOBMH 1973 K[43]
SiO2-CaO-Al2O3-Na2OBMH 1873 K[44]
SiO2-CaO-Al2O3BMH 1773 K[45]
SiO2-CaO-FeOBMH 1873 K[46]
BasicitySiO2-CaO-Al2O3-FeOBMH and LJ 1873 K[1]
SiO2-CaO-FeO-P2O5BMH 1673 K[47]
SiO2-CaO-Al2O3-MgOBMH 1873 K[48]
SiO2-CaO-Al2O3-MgO-TiO2BMH 1773 K[33]
SiO2-CaO-B2O3BMH 1873 K[49]
MgOSiO2-CaO-Al2O3-MgO-TiO2BMH 1773 K[33]
SiO2-MgO-Al2O3BMH 1773 K[40]
SiO2-CaO-Al2O3-MgOBMH 1773 K[41]
TiO2-SiO2-MgO-CaOBMH 1973 K[43]
SiO2-MnO-MgO-B2O3BMH and LJ 1873 K[50]
MgO/Al2O3CaO-MgO-Al2O3-SiO2BMH 1773 K[51]
CaO-MgO-Al2O3-SiO2BMH 1773 K[52]
CaO-MgO-Al2O3-SiO2COMPASS 1773 K[53]
FeOFeO-TiO2BMH and LJ 1973 K[54]
SiO2-CaO-FeOBMH 300 K[55]
FeO-TiO2-SiO2BMH and LJ 1973 K[56]
CaO-SiO2-Al2O3-FeOBMH 1773 K[57]
SiO2-CaO-Al2O3BMH 1773 K[45]
CaO-SiO2-Al2O3-FeOBMH 1873 K[58]
CaO-SiO2-FeO-P2O5BMH 1673 K[59]
Na2O or K2OCaO-Al2O3-SiO2-Na2O(K2O)BMH and MORSE 2273 K[60]
CaO-Al2O3-SiO2-Na2O-K2OBMH and MORSE 2273 K[61]
P2O5-Fe2O3-FeO-Na2OBMH 300 K[62]
SiO2-CaO-Al2O3-MgO-Na2OBMH and MORSE 1773 K[63]
CaO-SiO2-Na2OBMH 1773 K[64]
SiO2-Al2O3-CaO-Na2OBMH 1873 K[44]
TiO2CaO-SiO2-TiO2BMH 1723 K[65]
CaO-SiO2-Al2O3TiO2BMH 1773 K[66]
CaO-SiO2-TiO2BMH 1873 K[67]
FeO-TiO2-SiO2BMH and LJ 1973 K[57]
CaO-SiO2-TiO2BMH 1873 K[38]
TiO2− MgO-Al2O3BMH 1973 K[68]
FeO-TiO2-B2O3BMH 1823 K[69]
MnOCaO-SiO2-Al2O3-MnOBMH and LJ 1773 K[70]
SiO2-MnO-MgO-B2O3BMH and LJ 1873 K[51]
P2O5CaO-P2O5-SiO2BMH 1673 K[71]
V2O5-P2O5BMH 298 K[72]
CaO-SiO2-P2O5BMH 3000 K[73]
P2O5-based binary systemBMH 1873 K[74]
SiO2-Al2O3-Na2O-P2O5Modified Teter 300 K[75]
CaO-SiO2-P2O5BMH 1673 K[76]
V2O5FeO-SiO2-V2O3BMH 1823 K[77]
B2O3Li2O-V2O5-B2O3BMH and MORSE 298 K[78]
CaO-Al2O3-B2O3BMH 1973 K[79]
SiO2-CaO-B2O3BMH 1773 K[35]
SiO2-CaO-Al2O3-B2O3BMH 1773 K[80]
SiO2-CaO-Al2O3-B2O3BMH 1773 K[37]
B2O3-TiO2-CaO-MgO-SiO2BMH 1973 K[81]
SiO2-CaO-B2O3BMH 2573 K[82]
CaO-MgO-B2O3-Al2O3-SiO2BMH and short order correction300 K[83]
CaFCaO-CaF2-SiO2BMH 1773 K[84]
CaO-SiO2(Al2O3)-CaF2BMH 1873 K[85]
CaO-CaF2-SiO2BMH 1823 K[86]
Li2OCaO-SiO2-Al2O3-Li2OBMH 1673 K[87]
CaO-SiO2-Al2O3-Li2OBMH 1673 K[88]
SiO2/Al2O3 ratioSiO2-MnO-CaF2-Al2O3BMH 1823 K[89]
CaO-SiO2-Al2O3BMH 2000 K[90]
CaO-SiO2-Al2O3-MgO-FeO-TiO2-Na2OBMH 1800 K[91]
Table 2. Representative research of AIMD in molten slag.
Table 2. Representative research of AIMD in molten slag.
SystemExchange–Correlation Functionals
and Pseudopotential
TemperatureReference
B2O3GGA-PBE, PAW Method800 K, 2300 K, 3600 K[94]
CaO-SiO2-Al2O3GGA-PW91, PAW Method1773 K[95]
FeO-SiO2-SGGA-PBE, PAW Method2000 K[96]
SiO2-P2O5-Al2O3-Na2OGGA-PBE, PAW Method300 K, 1000 K, 1500, 2000 K[97]
CaO(MnO)-SiO2GGA-PBE, PAW Method2000 K[98]
CaO-Al2O3-B2O3GGA-PBE, GHT-DZVP1873 K[99]
Liquid and Amorphous SiO2LDA, ultrasoft PP300K, 3000 K[100]
SiO2-Al2O3-CaOGGA-PBE, PAW Method2000 K[101]
Projector-augmented wave (PAW) [102]; Perdew–Burke–Ernzerhof (PBE) [103]; generalized gradient approximation (GGA) [103].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jiang, C.; Li, K.; Bi, Z.; Ma, S.; Zhang, J.; Liu, B.; Li, J. Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods. Nanomaterials 2024, 14, 464. https://doi.org/10.3390/nano14050464

AMA Style

Jiang C, Li K, Bi Z, Ma S, Zhang J, Liu B, Li J. Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods. Nanomaterials. 2024; 14(5):464. https://doi.org/10.3390/nano14050464

Chicago/Turabian Style

Jiang, Chunhe, Kejiang Li, Zhisheng Bi, Shufang Ma, Jianliang Zhang, Bo Liu, and Jiaqi Li. 2024. "Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods" Nanomaterials 14, no. 5: 464. https://doi.org/10.3390/nano14050464

APA Style

Jiang, C., Li, K., Bi, Z., Ma, S., Zhang, J., Liu, B., & Li, J. (2024). Developments in Atomistic and Nano Structure Evolution Mechanisms of Molten Slag Using Atomistic Simulation Methods. Nanomaterials, 14(5), 464. https://doi.org/10.3390/nano14050464

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop