Latest ArticlesCement production accounts for approximately 12% of China's total CO2 emissions, having a significant challenge to achieving the national "dual carbon goals". Carbon capture, utilization, and storage (CCUS) represent a pivotal innovative technology for mitigating these emissions. However, conventional amine-based CO2 capture requires an energy-intensive high-temperature desorption, hindering its industrial implementation in cement plants. Also, the limited utilization pathways for captured CO2 pose another challenge for cement CCUS. Diethanolamine (DEA) offers a promising solution as it functions both as a CO2 absorber and a cement additive. This dual capability enables a potential carbonation utilization of CO2 absorbed DEA solutions without requiring the desorption step within cementitious systems. This study was thus to investigate the effect of CO2-absorbed diethanolamine (DEAC) on the early hydration behavior and strength development of cementitious systems. The findings could propose a novel approach for low-energy CO2 capture coupled with efficient in-situ utilization within cement industry.
Cement mortars with a water-to-cement ratio (W/C) of 0.50 were prepared with P·I 42.5 Portland cement (GB 8076) and ISO standard sand. The specimens were designated as REF, D0.1%C0%, D0.1%C0.02%, D1.0%C0%, and D1.0%C0.22%, respectively. DEA and its equivalent CO2 admixture were added as percentages of cement mass. All associated cement paste mixtures were prepared at a W/C ratio of 0.3. DEAC was prepared by continuously bubbling CO2 gas (≥99% purity) at a flow rate of 200 mL/min through a 5 mol/L DEA solution maintained at 40 ℃ until saturation was achieved. An eight-channel microcalorimeter recorded the hydration heat of cement paste specimens at 25 ℃ for 72 h. The phase composition of hardened cement pastes was determined by X-ray diffraction (XRD). The contents of bound water, CH, and CaCO3 were analyzed by thermogravimetric analysis (TGA). The cumulative porosity and pore size distribution of hardened paste samples at 7 d were characterized by mercury intrusion porosimetry (MIP). The compressive strength was measured on mortar specimens at 1, 3 d and 7 d of curing in accordance with the standard GB/T 17671.
The hydration calorimetry results demonstrate that D0.1%C0.02% and D0.1%C0% both accelerate the hydration rate of silicate phases, as evidenced by an increased second exothermic peak rate, while leaving the induction period duration unaffected. Conversely, D1.0%C0% and D1.0%C0.22% significantly reduce the second exothermic peak rate. D1.0%C0.22% extends the hydration induction period to 240 min, while D1.0%C0% has a negligible effect on its duration. The XRD patterns and TG analyses reveal that the impact of DEAC on the cement hydration depends critically on its specific DEA and CO2 dosage. At a low dosage (i.e., D0.1%C0.02%), a mild carbonation promotes a concurrent hydration of silicate and aluminate phases. However, a high dosage (i.e., D1.0%C0.22%) substantially inhibits early hydration of silicates. The MIP results indicate that DEAC and DEA both refine the pore structure of hardened cement paste. The pores below 20 nm are significantly reduced in D0.1%C0% and D0.1%C0.02% systems, aligning with their enhanced early hydration kinetics. This refinement also occurres in D1.0%C0% and D1.0%C0.22% systems despite inhibited silicate hydration. The results of compressive strength tests show that D0.1%C0.02% and D0.1%C0% can enhance mortar strengths at 1, 3 d, and 7 d, respectley. The strengths of D0.1%C0.02% systems can be increased by 8.4%, 10.2%, and 16.8% at these ages, respectively, primarily due to the DEA component with the weak carbonation contributing minimal additional enhancement. The 3-day and 7-day strengths of D1.0%C0.22% and D1.0%C0% systems both are increased (more significantly in the carbonated system), indicating a synergistic hydration-carbonation effect. However, the 1-day strength of D1.0%C0.22% system drastically is reduced by 52.2%, with silicates hydration inhibition by D1.0%C0% identified as a primary factor underpinning early strength reduction. According to the analysis of bound water content, calcium hydroxide (CH) content, porosity, and compressive strength relationships, a linear correlation between CH content and mortar strength is proposed. This demonstrates that silicate phase hydration kinetics can be modulated differently by DEAC and DEA formulations-fundamentally governed compressive strength development.
The addition of 0.1%DEA with 0.02% CO2 (D0.1%C0.02%) as DEAC enhanced the flexural and compressive strengths of cement mortar at 1, 3 d, and 7 d. In contrast, the addition of 1.0% DEA with 0.22% CO2 (D1.0%C0.22%) significantly reduced the 1-day strength. In D0.1%C0.02% system, the CO2 component reacted with dissolved Ca2+ released from cement minerals to precipitate CaCO3. This reaction promoted cement hydration, refined the pore structure of the hardened paste by reducing the volume of harmful pores, and facilitated a synergistic enhancement of hydration and carbonation. D1.0%C0.22% addition significantly retarded cement hydration within the first 24 h, primarily by inhibiting the dissolution of silicate phases and extending the induction period, leading to the reduced early strength. Although carbonation exacerbated the retardation of silicate phase hydration via DEA interaction, the hydration process recovered normal kinetics after 7 d.
The Faraday effect is one of the magneto-optical phenomena and refers to the conversion of linearly polarized light passing through a magnetic material into elliptically polarized light with the main axis-containing polarization plane rotated around the propagation vector. The angle by which the polarization plane is rotated, i.e., the Faraday rotation angle, is an important parameter determining the applicability of magnetic materials in devices such as electric-current and magnetic-field sensors, optical isolators, and optical circulators. Since the Faraday effect deals with a transmitted light, the transmittance of the magnetic materials is another important factor for applications. Thus, the materials are required to show a great magneto-optical figure of merit, that is defined as Faraday rotation angle or Verdet constant divided by absorbance or optical absorption coefficient. Here, the Verdet constant is defined as the Faraday rotation angle divided by external magnetic field and light path length inside the magnetic materials. It is well known that single crystals of garnet-type ferrites such as Y3Fe5O12 and (Gd,Bi)3Fe5O12 exhibit a large Faraday effect and a low optical absorption in the infrared region, especially in a wavelength range from 1.3 μm to 1.5 μm, and that they are effectively utilized as an optical isolator for optical telecommunications. However, compared to the garnet-type ferrites in the infrared region, magneto-optical materials with the superior performance, are lacking in the visible to ultraviolet region. Hence, the development of such materials is still in progress.
Oxide glasses rich in rare-earth ions exhibit a great Faraday effect, especially in the visible to ultraviolet range. Although these glasses feature magnetizations smaller than those of ferro- or ferri-magnetic oxide crystals such as abovementioned Y3Fe5O12 because the rare-earth-containing glasses are usually paramagnetic at room temperature, the transmittance of these glasses notably exceeds that of ferrite crystals in the visible to ultraviolet range. In addition, oxide glass has an advantage that it is feasible to tune continuously the composition so that optimized properties are attained and to fabricate large-sized and specific-shaped materials. In addition to the paramagnetic glasses, the Faraday effect of diamagnetic glasses is intensively investigated as well. The magnetization of diamagnetic glasses is further smaller than that of paramagnetic glasses, but the Faraday rotation angle or the Verdet constant of diamagnetic glasses is almost independent of temperature. This is an advantageous point of diamagnetic glasses, which cannot be realized in ferro-magnetic, ferri-magnetic, and para-magnetic materials. Furthermore, for wide-band gap oxide glass like SiO2 glass, which is diamagnetic, the Faraday effect can occur even in a very short wavelength range such as the deep and vacuum ultraviolet.
This review represents recent development on oxide glasses exhibiting large Faraday rotation. The macroscopic and microscopic mechanism of the Faraday effect are explained. The microscopic mechanism is very important to select magneto-optically active elements and to design glass compositions.Also, the Faraday effect of diamagnetic glasses is described. Heavy-metal oxide glasses and sulfide glasses are intensely exploited because the magnetic susceptibility of diamagnetic materials depends on the constituent atoms (ions) and the susceptibility is proportional to the squared atomic (ionic) radius and the number of electrons contained in the atom (ion). The Verdet constants of these glasses are summarized. The applications of diamagnetic glasses are briefly mentioned.
Subsequently, the Faraday effect of paramagnetic oxide glasses containing large amounts of rare-earth ions is reviewed. The pioneering work in this field has been carried out in the mid-1960s, showing that some ions like Ce3+, Pr3+, Tb3+, Dy3+, and Eu2+ give rise to larger Verdet constants in the visible range. A description is given to explain why these rare-earth ions exhibit larger Faraday effects than other ones. Recent researches seem to mainly pay attention to Tb3+-rich oxide glasses, for which higher concentrations of Tb3+ ions simply enhance the Verdet constant. In particular, Tb3+-rich oxide glasses fabricated via containerless processing, which is an emerging method and effective to expand the glass-forming region, showing the larger Verdet constant than single-crystalline Tb3Ga5O12 used as a commercially available optical isolator in the visible range. Furthermore, EuO-based amorphous oxides that have an unexpected ferromagnetism exhibit rather large Faraday effect.
In addition to the abovementioned diamagnetic oxide glasses and rare-earth-rich oxide glasses, a brief review concerns the Faraday effect of oxide glasses containing large amounts of 3d transition metal ions as well as glass-ceramics comprising ferro- or ferri-magnetic nano-sized crystalline particles embedded in transparent glass matrices.
The Faraday effect was discovered 180 years ago, but this phenomenon has been still utilized for practical applications as mentioned above. In particular, Tb3+-rich and Eu2+-rich oxide glasses are important for both fundamentals and applications. The Tb3+-rich glasses show a high transparency even in blue to ultraviolet region, so that the magneto-optical figure of merit is large enough to apply for an optical isolator. The Eu2+-rich glasses are ferromagnetic, so that they notably show a large Faraday effect. A new technique of glass formation such as containerless processing is effective to produce new glass compositions with further higher concentrations of rare-earth ions that are expected to exhibit a larger Verdet constant. Besides, the possible enhancement of Faraday effect based on plasmonics and Mie-tronics, i.e., the usage of localized surface plasmon resonance of metal nanoparticles and the Mie resonance of dielectric nanoparticles to increase the Verdet constant, becomes an important subject in the near future. With the development of high-power lasers, the demand for optical isolators that can operate in a wide wavelength range must increase. The oxide glasses have a promising application in such fields.
Solid oxide fuel cells (SOFCs) are high-efficient solid-state energy conversion devices. However, all-ceramic self-supported SOFCs face several challenges such as high brittleness, difficulty in mechanical processing, poor thermal shock resistance, and limited weldability, which result in high manufacturing costs and restrict the applications in mobile power systems. In contrast, metal-supported SOFCs (MS-SOFCs) with metal materials as the external structural support, exhibit remarkable mechanical strength, low cost, and rapid start-up capability, making it highly promising for mobile applications. The anode is a critical component of MS-SOFCs, serving as the site where fuel oxidation occurs to generate electrons. Its microstructure significantly influences the density and effectiveness of the triple-phase boundaries (TPB), where the gas phase, the ionic phase, and the electronic phase intersect. The TPB density largely determines the polarization resistance, with its low-frequency component being inversely related to anode gas diffusion. A common strategy to enhance gas transport is the incorporation of pore-formers, such as graphite, into the anode raw materials. Most studies focus on the type, particle size, and content of pore-formers, which directly affect the pores number, size, and distribution. In this work, atmospheric plasma spraying (APS) was employed to fabricate three types of anodes and corresponding cells. APS reduces thermal input to the metal substrate, effectively preventing oxidation, deformation, and elemental interdiffusion between the metal support and the anode at high temperatures. This study systematically investigates the influence of graphite incorporation methods on the microstructural evolution of the anode and the resulting cell performance, providing important theoretical insights into the operational mechanisms of SOFC anodes.
Porous 430L stainless steel substrates were used as supports. Three different NiO-GDC (Gd0.2Ce0.8O1.9) anode powders were prepared. C1 is the baseline without a pore-former. C2 contains 40% (in volume fraction) graphite, which is mixed by spray granulation process to produce composite particles. C3 is made by mechanically mixing 40% of the same graphite with C1 powder. All powders are spherical with good fluidity. Anode layers were deposited via APS. Subsequently, the C2 and C3 anodes were heat-treated in air to remove the graphite pore-former at 750 ℃ for 2 h. A ScSZ (Sc2O3-ZrO2) electrolyte and an LSCF (La0.6Sr0.4Co0.2Fe0.8O3-δ) cathode were subsequently deposited by APS to build single cells. The microstructure and porosity of the anodes were characterized using scanning electron microscopy (SEM) and image analysis software. The surface roughness was measured by profilometer. The electrochemical performance, including the open-circuit voltage (OCV), current-voltage-power (I-V-P) and electrochemical impedance spectroscopy (EIS), were evaluated in the range of 600-750 ℃ using humidified H2 as fuel and air as oxidant. The equivalent circuit fitting of EIS data is carried out to quantitatively analyze the contribution of charge transfer, surface adsorption/dissociation and gas diffusion in polarization impedance.
The incorporation of graphite pore-former significantly modified the pore size distribution and total porosity within the anodes. The measured porosity of C1, C2, and C3 was 26%±2%, 37%±3.1%, and 42%±2.3%, respectively. The C1 anode featured a relatively dense structure with uniformly distributed pores, which primarily consisted of submicron cracks and fine pores originating from the thermal stress inherent to the APS process. In contrast, the C2 anode showed a notable increase in both the number and size of pores, which were homogeneously dispersed without significant agglomeration. The C3 anode, however, contained a substantial amount of large pores, mostly ~5 μm in diameter, attributed to the agglomeration of graphite particles during mechanical mixing, resulting in coarse and irregular pore structures after heat treatment. Furthermore, the addition of graphite modified the thermal response characteristics of the agglomerated powder during the spraying process, promoting the formation of a more uniformly melted microstructure in the anode layer. The average surface roughness (Ra) values for C1, C2, and C3 were 6.64, 7.06 nm, and 7.66 μm, respectively, indicating that graphite addition increased anode surface roughness. This phenomenon is due to the thermal decomposition of graphite during the plasma spraying process, where high temperatures cause partial oxidation of graphite in the open atmosphere, thereby generating CO2 gas. The release of this gas from the incompletely solidified anode surface etches irregular pits and protrusions, ultimately leading to increased surface roughness. The cell without graphite pore-former (C1) consistently demonstrated the highest OCV and maximum power density, reaching 1.0 V and 957 mW·cm-2 at 750 ℃, respectively. EIS analysis revealed that the anodes with graphite pore-former (C2 and C3) exhibited improved charge transfer capability, thereby reducing the high-frequency polarization resistance. Despite this, the overall output performance of C2 and C3 did not show effective enhancement, which is attributed to their increased ohmic resistance (Ro) and lower OCV. The elevated surface roughness and inherent porosity in the C2 and C3 anodes adversely affected the quality of the subsequently sprayed electrolyte layer, introducing microcracks and gas permeation pathways. This resulted in increased Ro and reduced OCV, ultimately weakening the benefits gained from the reduced polarization resistance.
The method of graphite pore-former addition significantly affects the anode microstructure and overall cell performance. Compared with mechanical mixing, spray granulation produces a superior and uniform pore structure. However, contrary to conventional expectation, the introduction of graphite pore-forming agent into the APS anode reduces overall cell performance due to induced electrolyte defects, which elevated Ro and lowered OCV. Future optimization should focus on strategies to reduce anode surface roughness and refine pore structure without affecting the quality of electrolyte deposition. It is anticipated that this will further enhance the output performance of MS-SOFCs.
Scintillators are functional materials that immediately emit low-energy photons after absorbing high-energy ionizing radiations. Glasses have advantages as a material form for scintillator use, such as low manufacturing costs, ease of formability, high optical transparency, and high compositional tunability. A 6Li-glass scintillator doped with Ce is a commercial glass scintillator, and it has been used for thermal neutron detections owing to 6Li(n,α)3H neutron capture reactions. In general, the 3He counter is used in practice for neutron detections; however, the depletion of the 3He supply has driven vigorous research and development of alternative materials. 10B-contained glasses have an advantage of their large thermal neutron capture cross section (3840 barn). In this study, we focused on MgO-P2O5-B2O3 glass systems. They are composed of light elements; hence, neutron detection signals can be easily distinguished from noise due to X- and γ-rays. Furthermore, alkali-earth-embedded P2O5-B2O3 glasses realize high optical transparency with good chemical durability. As luminescence centers, Eu was selected. Eu exists in two different states: Eu2+ in reduction states and Eu3+ in oxidation states. The former shows broad emission bands with fast decay times of ~μs, whereas the latter shows sharp emission bands with slow decay times of ~ms. In general, Eu3+ is governed in the glasses prepared in an air atmosphere; however, alkali-earth-embedded P2O5-B2O3 glasses can exhibit Eu2+ owing to the localized reducing atmosphere generated by NH3 derived from the raw material (NH4H2PO4). Here, Eu-doped MgO-P2O5-B2O3 glasses with various concentrations of Eu were synthesized by the conventional melt quenching technique in air, and their optical, photoluminescence (PL), and scintillation properties were examined.
Undoped and 0.1%, 0.3%, 1.0%, and 2.0% (in mole fraction) Eu-doped 25MgO-30P2O5-45B2O3 (MPB) glasses were synthesized by the melt quenching method. Raw materials were homogeneously mixed and transferred into an alumina crucible. The powders were melted at 1100 ℃ for 1 h with an electrical furnace in air. The melt was flowed onto a preheated stainless-steel plate to quench, and pressed into batches. After that, the glasses were annealed at 400 ℃ for 1 h to remove thermal and mechanical strains. The surfaces of the glasses were polished for the following optical and scintillation measurements. The glass transition temperature (Tg) of the undoped sample was measured with a TG-DTA system. The powder X-ray diffraction (XRD) patterns were measured using a diffractometer. Diffuse transmission spectra were measured using a spectrophotometer (Shimadzu, SolidSpec-3700). PL excitation and emission spectra, PL quantum yields (QYs), and PL decay curves were measured using a Quantaurus-QY and Quantaurus-τ. Scintillation spectra, scintillation decay curves, and afterglow curves under X-ray irradiations, and pulse height spectra of 241Amα-rays and 137Cs γ-rays were measured with our original setups.
The appearances of undoped and Eu-doped MPB glasses were transparent and included some bubbles. Under ultraviolet light at 360 nm, Eu-doped samples showed bluish-red light. Some parts of the glasses were crushed into powder and used for the XRD measurements. Precipitations of crystalline phases were not observed in the XRD patterns; hence, the prepared samples formed glass phases with no periodical structures. Tg of the undoped glass was estimated to be 535 ℃. The transmittances of all the glasses were 70%-90% at 400-700 nm. Both the undoped and Eu-doped samples showed absorption peaks at 200-250 nm; hence, this can be due to the hosts. In addition, absorption peaks emerged at 250-400 nm in Eu-doped samples. The absorption bands at 250-280 nm and 280-400 nm originated from the 4f7-4f65d1 (T2g and Eg) transitions of Eu2+, respectively. An absorption peak due to the 7F1-5D3 transitions of Eu3+ was confirmed at 415 nm in the spectra of 2% Eu. A broad emission band due to the electronic transitions of Eu2+ was observed at 400-600 nm under excitation bands at 250-410 nm in the Eu-doped samples. Generally, emissions due to Eu2+ are difficult to observe in the glasses prepared in air because of the oxidizing atmosphere. The emissions appeared owing to the localized reduction atmosphere, which derived from NH4H2PO4 and (NH4)2O·5B2O3·8H2O. When excitation and monitored wavelengths were respectively set to 340 nm and 420 nm, the PL decay curves were fitted by a sum of two exponential functions. The PL decay time constants of both the fast (0.2-0.5 μs) and slow (0.7-1.2 μs) components were similar to those of other Eu-doped phosphors; hence, they are reasonable as 4f65d1(T2g)-4f7 and 4f65d1(Eg)-4f7 transitions of Eu2+. Eu-doped samples showed emission peaks at 350-500 nm under X-ray irradiations. These emission wavelengths were consistent with those of the PL spectra. Afterglow levels (AL) of undoped and 0.1%, 0.3%, 1.0%, and 2.0% Eu-doped glasses at 20 ms passed after X-ray irradiation of 2 ms were respectively estimated to be 1700×10-6, 4500×10-6, 3400×10-6, 800×10-6, and 160×10-6. Pulse height spectra of 241Am α-rays (5.5 MeV) were measured using the prepared MPB glasses. Eu-doped glasses showed clear full energy absorption peaks. Light yields (LY) of 0.3%, 1.0%, and 2.0% Eu-doped MPB glasses were respectively calculated to be 70, 150 photons/5.5 MeV, and 40 photons/5.5 MeV.
Eu-doped MPB glasses were synthesized by the conventional melt quenching method. All the samples showed halo peaks with no periodic patterns in XRD measurements. The transmittances were 70%-90%, and absorptions due to electronic transitions of Eu2+ and Eu3+ were observed. All the samples showed luminescence, which originated from the 4f65d1-4f7 transitions of Eu2+. PL QYs of 0.1%, 0.3%, 1.0%, and 2.0% Eu-doped samples were respectively calculated to be 46.7%, 37.2%, 20.4%, and 9.3% when monitored at 350-560 nm under excitation at 320 nm. ALs at 20 ms passed after X-ray irradiations of 2 ms were obtained to be 1700×10-6, 4500×10-6, 3400×10-6, 800×10-6, and 160×10-6 in undoped, 0.1%, 0. 3%, 1.0%, and 2.0% Eu-doped samples, respectively. Pulse height spectra of 241Am α-rays were measured using the prepared samples. LYs of 0.3%, 1.0%, and 2.0% Eu-doped samples were respectively determined to be 70, 150 photons/5.5 MeV, and 40 photons/5.5 MeV.
Crystallization control is crucial during glass production. In glass crystallization theory, it has been hypothesized that there is a direct relationship between glass-crystal interfacial energy and nucleation rate. Since it is difficult to measure interfacial energy directly, classical nucleation theory is used to obtain it. However, the estimation process is complex; various measurement data such as the nucleation rate, viscosity, and Gibbs free energy barriers between the glass and crystal are required. Overall, estimating the glass-crystal interfacial energy is time-consuming. In this regard, the interfacial energy can be directly obtained using molecular dynamics (MD) simulations. In this study, we implemented an interface model to estimate the glass-crystal interfacial energy based on the theory of Tielemann et al., which was recently introduced to generate crystal orientation planes using minimum-energy cuts. To the best of our knowledge, this is the first work in which the crystal orientation plane determined based on a minimum-energy cutting process has been used to build a glass-crystal interface in a more realistic environment (i.e., hypothesizes that a minimum-energy structure may occur during crystal nucleation) and calculate the interfacial energy. To achieve this, we considered stoichiometric alkali (Li, Na, and K) disilicate (2SiO2) glasses and crystals, as some early-stage experimental data have been reported, which are highly beneficial for validating the MD results. The effects of different potential models were also investigated and found that they had a significant impact on the reproduction of the experimental trend.
Simulations were performed using the LAMMPS package. The temperature and pressure were controlled with the Nose-Hoover thermostat and barostat, respectively. After setting up the glass-crystal interface using NPT, we ran MD simulations with NVT for another 200 ps to calculate the interfacial energy. The MD simulations were performed with a time step of 1 fs. Periodic boundary conditions were applied in all directions. Cutoff distances were applied according to references. We used three interatomic potential (SHIK, Du, and Pedone) models to estimate the interfacial energies of Li2O-2SiO2 glass-Li2O-2SiO2 (001), Na2O-2SiO2 glass-Na2O-2SiO2 (010), and K2O-2SiO2 glass-K2O-2SiO2 (001). All three of these potential functions are widely used. The glass structure was fabricated from the crystal structures by simulating a melt-quenching process. First, the Li2O-2SiO2 (001), Na2O-2SiO2 (010), and K2O-2SiO2 (001) crystal structures were prepared. Half the crystal was kept fixed, while the other half was melted at 3500 K with a canonical ensemble (NVT) for 300 ps and then quenched to 300 K at a cooling rate of 5 K/ps. After the melt-quenching process, the crystal part was unfixed and relaxed for 200 ps with a glassy structure using an isobaric-isothermal ensemble (NPT) at a temperature of 300 K. An interfacial model was developed by the theory of Tielemann et al, where the crystal plane was determined through the minimum energy cut in the crystal structure.
We first evaluated the glass density by calculating the local density profile along the z-direction of Li2O-2SiO2 glass-Li2O-2SiO2 (001), Na2O-2SiO2 glass-Na2O-2SiO2 (010), and K2O-2SiO2 glass-K2O-2SiO2 (001). The glass structures show density results comparable to the experimental data.
The calculated interfacial energy values using the SHIK potential show a similar experimental trend, while the Du and Pedone potentials are unable to reproduce the experimental trend. However, the potentials of Du and Pedone show better values of interfacial energy for the glass-crystal interface of Li2O-2SiO2 than the SHIK. Our MD results demonstrate that among the potential models, SHIK is a good candidate for calculating interfacial energy in terms of experimental reproduction.
We analyzed the interfacial coordination number (i.e., cation-anion) in the contact area between the glass and crystal surfaces. These coordination numbers are difficult to estimate experimentally. Typically, the coordination number is determined in the bulk region of a glass or crystal structure. We found that the interfacial coordination number varies significantly at the glass-crystal interface of Li2O-2SiO2, Na2O-2SiO2, and K2O-2SiO2. Among the interfacial systems, coordination number values of two, three, four, and five were observed, the most frequently observed value was one. The values were significantly lower at four and five for K2O-2SiO2 and Na2O-2SiO2, respectively. Our MD results demonstrated that the cation (Li, Na, and K)-anion environments were not the same in the interfacial domain. For further evaluation, we calculated the bond strength-coordination number.
The bond strength-coordination number has been reported for chalcogenide glasses. In the present work, the bond strength-coordination number is introduced for glass-crystal interface systems. The calculated results of bond strength-coordination number show a decreasing trend in the order Li > Na > K, which is similar to interfacial energy. Overall, analysis revealed that the alkali (Li, Na, and K)—O bond plays a crucial role in the interfacial strength.
We developed an interface model using molecular dynamics simulations based on the minimum energy cut of crystals with glassy structures. The minimum energy cut of the crystal was determined by applying Tielemann theory. The modeled interfaces were between glass and crystals in stoichiometric alkali (Li, Na, and K) disilicates. The interface model reproduced the experimental interfacial energy trend of Li > Na > K; the interatomic potential model was found to have a significant effect on reproduction. The MD results indicate that the SHIK-based MD model could reproduce the experimental results better than the other potential models. In addition, the interfacial coordination number was also calculated, which varies depending on the alkalinity of Li, Na, and K. Using the interfacial coordination number, the bond strength-coordination number was estimated, and found that the interfacial energy trend Li > Na > K is related to the bond strength-coordination number.